• Structural Engineering and Mechanics
• /
• 제35권6호
• /
• pp.677-697
• /
• 2010

This paper focuses on geometrically non-linear static analysis of a simply supported beam made of hyperelastic material subjected to a non-follower transversal uniformly distributed load. As it is known, the line of action of follower forces is affected by the deformation of the elastic system on which they act and therefore such forces are non-conservative. The material of the beam is assumed as isotropic and hyperelastic. Two types of simply supported beams are considered which have the following boundary conditions: 1) There is a pin at left end and a roller at right end of the beam (pinned-rolled beam). 2) Both ends of the beam are supported by pins (pinned-pinned beam). In this study, finite element model of the beam is constructed by using total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In order to use the solution procedures of Newton-Raphson type, there is need to linearized equilibrium equations, which can be achieved through the linearization of the principle of virtual work in its continuum form. In the study, the effect of the large deflections and rotations on the displacements and the normal stress and the shear stress distributions through the thickness of the beam is investigated in detail. It is known that in the failure analysis, the most important quantities are the principal normal stresses and the maximum shear stress. Therefore these stresses are investigated in detail. The convergence studies are performed for various numbers of finite elements. The effects of the geometric non-linearity and pinned-pinned and pinned-rolled support conditions on the displacements and on the stresses are investigated. By using a twelve-node quadratic element, the free boundary conditions are satisfied and very good stress diagrams are obtained. Also, some of the results of the total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element are compared with the results of SAP2000 packet program. Numerical results show that geometrical nonlinearity plays very important role in the static responses of the beam.

• 한국압력기기공학회 논문집
• /
• 제19권2호
• /
• pp.75-83
• /
• 2023

For three-dimensional finite element welding residual stress simulation, several methods are available. Two widely used methods are the moving heat source model using heat flux and the temperature boundary condition model using the temperature profile of the welded beads. However, each model has pros and cons in terms of calculation times and difficulties in determining welding parameters. In this paper, a new method using the moving temperature profile model is proposed to perform efficiently 3-D FE welding residual stress analysis for large structures. Comparison with existing experimental residual stress measurement data of two-pass welding pipe and SNL(Sandia National Laboratories) mock-up canister shows the accuracy and efficiency of the proposed method.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2000년도 봄 학술발표회 논문집
• /
• pp.669-676
• /
• 2000

The estimation of consolidation rate is one of the important factors in the construction on soft clayey deposits. A number of researches are carried out to predict the consolidation behavior in field, however, most of the results show the discrepancies between the prediction and observation. This paper analyzes consolidation behavior of normally consolidated clay in K/sub o/ condition with 2-dimensional drainage by use of the numerical methods. Elastic and elastic-plastic finite element analyses are compared in terms of the dissipation of excess pore pressure. These results are also compared with Terzaghi-Rendulic's equation that is implemented by finite difference method. The consolidation time calculated by using elastic model is found to be similar to the result of Terzaghi-Rendulic's equation. The consolidation predicted by MCC model takes more time than other cases. Initial increase of excess pore pressure in radial drainage can be shown, however, this phenomenon does not have a significant effect on tile final consolidation time.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• 제3B권2호
• /
• pp.67-73
• /
• 2003

In this paper, the impedance of a piezoelectric transducer is calculated using the three-dimensional finite element method. The validity of numerical routine is confirmed experimentally. Using this numerical routine, the effects of material coefficients on piezoelectric actuators characteristics are analyzed. The material constants, which make significant effects, are selected and the relations between material constants are studied. Using these processes, three variables of material constants for a piezoelectric transducer are selected and the design sensitivity method is adopted as an inversion scheme. The validity of the inversion scheme for a piezoelectric transducer is confirmed by applying the proposed method to the sample piezoelectric transducer.

• 한국해양공학회지
• /
• 제8권2호
• /
• pp.80-92
• /
• 1994

The dynamic response characteristics of a floating ocean city are examined for presenting the basic data for the design of huge offshore structures supported by a large number of floating bodies in waves. The numerical approach which is accurate in linear system is based on combination of a three dimensional source distribution method, wave interaction theory and the finite element method of using the space frame element. The hydrodynamic interactions among the floating bodies are taken into account in their exact form within the context of linear potential theory in the motion and structural analysis. The method is applicable to an arbitrary number of three dimensional bodies having any individual body geometries and geometrical arrangement with the restriction that the circumscribed, bottom-mounted. Imaginary vertical cylinder for each body does not contain any part of the other body. The validity of this procedure was verified by comparing with experimental results obtained in the literature.

• 한국컴퓨터정보학회논문지
• /
• 제17권9호
• /
• pp.9-16
• /
• 2012

본 논문에서는 유한요소법에서 희소행렬의 효율적인 저장을 위한 2차원 가변길이 벡터 저장구조를 제안한다. 제안한 저장구조는 유한요소 전체 방정식의 거대희소행렬 $N{\times}N$ 대신, 전체 행의 개수 N의 상삼각행렬에서 0이 아닌 실제 필요한 값들만 2차원 가변길이 벡터를 이용하여 저장하는 방법이다. 이 방법을 이용하면, 해석대상의 2차원 격자구조에서는 각 절점당 최소 1개에서 최대 5개까지의 저장 공간이 필요하게 되고, 3차원 격자구조에서는 각 절점당 최소 1개에서 최대 14개까지의 저장 공간이 필요하게 된다. 인덱스를 포함해도 2배 이상을 넘지 않는다. 본 논문의 실험 결과에 의해, 제안한 저장구조는 총 절점 개수가 많아질수록 기존의 최대칼럼 높이를 저장하는 스카이 라인 저장구조보다 메모리 공간을 효과적으로 줄일 수 있는 구조임을 알 수 있었다.

• 대한조선학회지
• /
• 제12권2호
• /
• pp.35-42
• /
• 1975

This paper investigates the distribution of stress and its behavior at the Root Toe in fillet welding joint. Furthermore, the stress components and principal stresses in the fillet welds are calculated by the finite element method. The distribution of stresses obtained numerically by means of the finite element method is also compared with the experimental results of two dimensional photoelasticity. A Cover plate type and Center block type of fillet welds are used as models for the numerical calculations covering the variations of 2 W/M(thickness of main plate/thickness of cover plate)=1 through 2W/M=4. The results obtained in these studies are summarized as follows; 1) When W2/M values become small, the stress concentration factors of the Root are larger than of the Toe in a C-type. Its critical value is 2W/M=3.00. However, no critical value exists in a T-type. 2) For 2W/M Values being avove 3.5 in a C-type and above 4.0 in a T-type, $K_R$ and $K_{\tau}$ become 1. 3) According to the differences of 2W/M values, the differences in stress become increasing in the Root but become decreasing in the Toe. These differences, however, disappear as the free boundary surface is approached. 4) The stress concentration factors of both the Root and Toe obtained by means of the finite element method have somewhat lower values than obtained by the photoelasiticity. But their principal stress directions coincide in either method. 5) It proves beneficial to employ the finite element method for two-dimensional plane stress analysis in front fillet welding joint.

• 대한기계학회논문집A
• /
• 제29권2호
• /
• pp.293-302
• /
• 2005

In this paper, the thick-walled laminated, orthotropic as well as bimaterial, composite hollow spheres under impact pressure are analyzed in detail by using the semi-discrete finite element method with the Houbolt time-integration scheme which results in unconditionally stable transient numerical results. Numerical results are obtained by using the self-constructed spherically symmetric (one-dimensional) and axially symmetric (two-dimensional) finite element programs, and compared with the previous solutions by other researchers, being shown some of which are incorrect. The finite element package Nastran is also adopted for numerical comparison.

• Earthquakes and Structures
• /
• 제26권2호
• /
• pp.163-174
• /
• 2024

Dynamic equilibrium equations for finite element analysis were derived for the free field one-dimensional shear wave propagation through the horizontally layered soil deposits with the elastic half-space. We expressed Rayleigh's viscous damping consisting of mass and stiffness proportional terms. We considered two cases where damping matrices are defined in the total and relative displacement fields. Two forms of equilibrium equations are presented; one in terms of total motions and the other in terms of relative motions. To evaluate the performance of new equilibrium equations, we conducted two sets of site response analyses and directly compared them with the exact closed-form frequency domain solution. Results show that the base shear force as earthquake load represents the simpler form of equilibrium equation to be used for the finite element method. Conventional finite element procedure using base acceleration as earthquake load predicts exact solution reasonably well even in soil deposits with unrealistically high damping.

• 대한토목학회논문집
• /
• 제30권2D호
• /
• pp.135-141
• /
• 2010

국내 줄눈콘크리트 포장 다웰바의 경우 합리적 근거가 없이 경험에 의하거나 국외의 기준에 의해 설계시공되고 있다. 다웰바를 합리적으로 설계하기 위해서는 다웰바 거동에 대한 고찰이 요구되며 이에 3차원 유한요소해석을 이용하여 분석하였다. 다웰바의 거동을 평가하기 위해서 다웰바의 3차원 유한요소해석과 다웰바 거동에 대한 대표적인 모델인 Timoshenko 이론을 비교분석하였다. 본 연구에서는 다웰바의 사이즈 및 배치에 대한 합리적인 설계를 하기 위해서 3차원 유한요소해석을 사용하여 다웰바의 그룹작용(Dowel Group Action) 범위와 차량바퀴의 주행경로(Wheel Path)분포를 고려하여 슬래브 처짐을 감소시킬 수 있는 다웰바 최적배치를 제시하였다.