• Geomechanics and Engineering
• /
• 제33권4호
• /
• pp.369-380
• /
• 2023

The aim of this paper is to investigate stress analysis of semi-infinite soils consisting of two layers with twin rectangular tunnels under static loads. The region close to the ground surface and tunnel modelled within finite elements. In order to use a more realistic model, the far region is modelled within infinite elements. The material model of the layered soil is considered as elastic and isotropic. In the finite element solution of the problem, two dimensional (2D) plane solid elements are used with sixteen-nodes rectangular finite and eight-nodes infinite shapes. Finite and infinite elements are ordered to be suitable for the tunnel and the soils. The governing equations of the problem are obtained by using the virtual work principle. In the numerical process, the five-point Gauss rule is used for the calculation of the integrations. In order to validate using methods, comparison studies are performed. In the numerical results, the stress distributions of the two layered soils containing twin rectangular tunnels presented. In the presented results, effects of the location of the tunnels on the stress distributions along soil depth are obtained and discussed in detail. The obtained results show that the locations of the tunnels are very effective on the stress distribution on the soils.

• KCI Concrete Journal
• /
• 제13권2호
• /
• pp.67-74
• /
• 2001

In this study, a numerical model for the simulation of reinforced concrete columns subject to cyclic loading is presented. The model consists of three separate models representing concrete, reinforcing steel bars and bond-slip between a reinforcing bar and ambient concrete. The concrete model is represented by the plane stress plastic-damage model and quadrilateral finite elements. The nonlinear steel bar model embedded in truss elements is used for longitudinal and transverse reinforcing bars. Bond-slip mechanism between a reinforcing bar and ambient concrete is discretized using connection elements in which the hysteretic bond-slip link model defines the bond stress and slip displacement relation. The three models are connected in finite element mesh to represent a reinforced concrete structure. From the numerical simulation, it is shown that the proposed model effectively and realistically represents the overall cyclic behavior of a reinforced concrete column. The present plastic-damage concrete model is observed to work appropriately with the steel bar and bond-slip link models in representing the complicated localization behavior.

• Journal of Mechanical Science and Technology
• /
• 제15권11호
• /
• pp.1499-1506
• /
• 2001

Various transition elements are used in general for the effective finite element analysis of complicated mechanical structures. In this paper, a solid-to-beam transition finite element, which can b e used for connecting a C1-continuity beam element to a continuum solid element, is proposed. The shape functions of the transition finite element are derived to meet the compatibility condition, and a transition element equation is formulated by the conventional finite element procedure. In order to show the effectiveness and convergence characteristics of the proposed transition element, numerical tests are performed for various examples. As a result of this study, following conclusions are obtained. (1) The proposed transition element, which meets the compatibility of the primary variables, exhibits excellent accuracy. (2) In case of using the proposed transition element, the number of nodes in the finite element model may be considerably reduced and the model construction becomes more convenient. (3) This formulation method can be applied to the usage of higher order elements.

• Structural Engineering and Mechanics
• /
• 제23권5호
• /
• pp.487-504
• /
• 2006

Several new versatile two-dimensional p-version finite elements are developed. The element matrices are integrated analytically to guarantee the accuracy and monotonic convergence of the predicted solutions of the proposed p-version elements. The analysis results show that the convergence rate of the present elements is very fast with respect to the number of additional Fourier or polynomial terms in shape functions, and their solutions are much more accurate than those of the linear finite elements for the same number of degrees of freedom. Additionally, the new p-version plate elements without the reduced integration can overcome the shear locking problem over the conventional h-version elements. Using the proposed p-version elements with fast convergent characteristic, the elasto-plastic impact of the structure attached with the absorber is simulated. Good agreement between the simulated and experimental results verifies the present p-version finite elements for the analyses of structural dynamic responses and the structural elasto-plastic impact. Further, using the elasto-plastic impact model and the p-version finite element method, the absorber of the T structure on the Qiantang River is designed for its collision protection.

• 한국CDE학회논문집
• /
• 제2권3호
• /
• pp.175-185
• /
• 1997

Given the widespread use of the Finite Element Method in strength analysis, automatic mesh generation is an important component in the computer-aided design of parts and assemblies. For a given resolution of geometric accuracy, the purpose of mesh generators is to discretize the continuous model of a part within this error limit. Sticking to this condition often produces many small elements around small features in spite that these regions are usually of little interest and computer resources are thus wasted. Therefore, it is desirable to selectively suppress small features from the model before discretization. This can be achieved by low-pass filtering a CAD model. A spatial function of one dimension higher than the model of interest is represented using the Fourier basis functions and the region where the function yields a value greater than a prescribed value is considered as the extent of a shape. Subsequently, the spatial function is low-pass filtered, yielding a shape without the small features. As an undesirable effect to this operation, all sharp corners are rounded. Preservation of sharp corners is important since stress concentrations might occur there. This is why the LPF (low-pass filtered) model can not be directly used. Instead, the distances of the boundary elements of the original shape from the LPF model are calculated and those that are far from the LPF model are identified and removed. It is shown that the number of mesh elements generated on the simplified model is much less than that of the original model.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2003년도 가을 학술발표회 논문집
• /
• pp.444-451
• /
• 2003

In this study, a new three-dimensional finite element analysis model of high-speed railway bridges considering train-bridge interaction, in which various improved finite elements are used for modeling structural members, is proposed. The box-type bridge deck of a railway bridge is modeled by the NFS(Nonconforming Flat Shell) elements with 6 degrees of freedom. Track structures are idealized using the beam finite elements with the offset of beam nodes and those on Winkler foundation with two parameters. And, the vehicle model devised for a high-speed train is employed, which has an articulated bogie system. By Lagrange's equations of motion, the equations of motion of a bridge-train system can be formulated. Finally, by deriving the equations of the forces acting on a bridge considering bridge-train interaction the complete system matrices of total bridge-train system can be constructed. As numerical examples of this study, 2-span PC box-girder bridge is analyzed and results are compared with experimental results.

• Coupled systems mechanics
• /
• 제9권2호
• /
• pp.125-145
• /
• 2020

In this paper, we present a 3D thermo-hydro-mechanical coupled discrete beam lattice model of structure built of the nonisothermal saturated poro-plastic medium subjected to mechanical loads and nonstationary heat transfer conditions. The proposed model is based on Voronoi cell representation of the domain with cohesive links represented as inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities in axial and both transverse directions. The enhanced Timoshenko beam finite element is capable of modeling crack formation in mode I, mode II and mode III. Mode I relates to crack opening, mode II relates to in-plane crack sliding, and mode III relates to the out-of-plane shear sliding. The pore fluid flow and heat flow in the proposed model are governed by Darcy's law and Fourier's law for heat conduction, respectively. The pore pressure field and temperature field are approximated with linear tetrahedral finite elements. By exploiting nodal point quadrature rule for numerical integration on tetrahedral finite elements and duality property between Voronoi diagram and Delaunay tetrahedralization, the numerical implementation of the coupling results with additional pore pressure and temperature degrees of freedom placed at each node of a Timoshenko beam finite element. The results of several numerical simulations are presented and discussed.

• Journal of Welding and Joining
• /
• 제21권6호
• /
• pp.71-76
• /
• 2003

For the purpose of clarifying the influence of welding residual stress to the fatigue crack propagations behaviour, an analytical investigation based on finite element method is performed to examine the opening behaviour of tip-closed crack in the compressive residual stress. A finite element model comprised of contact elements for the crack plane and plane stress elements for the base material is used to evaluate crack opening stress of the crack existing in the residual stress field. Also an analytical method based on the superposition principle to estimate the length of opened part of tip closed crack and the stress distribution adjacent to the crack during uploading is applied to the finite element model. The software for the analysis is ABAQUS, which is a general purpose finite element package. The results show that stresses distributed on the crack surfaces are reduced and approached to zero as the applied stresses are increased up to crack tip opening stress and no mechanical discontinuity is found at the boundary of contact elements and plane stress elements. It is verified that the opening behavior of the fatigue crack in the residual stress can be predicted by finite element method with the proposed analytical method.

• 한국수자원학회논문집
• /
• 제42권12호
• /
• pp.1053-1067
• /
• 2009

본 연구의 목적은 하도의 형상이 불규칙한 자연하천에서 2차원 흐름 특성을 해석하고 예측하기 위해 2차 요소를 이용한 정확하고 효과적인 상향가중 유한요소모형의 개발에 있다. 모형의 개발을 위해 선형 삼각형 요소, 선형 사각형 요소와 혼합요소를 적용하였고 2차 삼각형, 사각형 요소와 혼합요소를 적용하여 모형을 개발하였으며, 지배방정식의 수치적분식으로 Gauss Quadrature 방법을 사용하였다. 개발된 모형의 적용성 검증을 위해 하상융기가 있는 수로, U자형 수로 등에 모의를 실시하여 해석해 및 실측치와 비교 검토하였다. 모의 결과 2차 요소가 선형 요소에 비해 보다 정확한 해를 제공하는 것으로 판단되었으며 2차요소를 적용한 상용모형인 RMA-2 모형과 비교한 결과 본 연구 개발 모형이 보다 정확한 해를 나타내는 것을 확인할 수 있었다. 개발된 모형을 향후 자연하천에 적용할 경우 기존의 모형에 비해 향상된 결과를 얻을 수 있을 것으로 판단된다.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2008년도 추계학술대회 논문집
• /
• pp.1379-1388
• /
• 2008

In this paper, we present the equations of motion by which the natural vibration of a rotating annular disk can be accurately analyzed. These equations are derived from the theory of finite deformation and the principle of virtual work. The radial displacements of annular disk which is rotating at constant angular velocity are determined by non-linear equations formulated using 1-dimensional finite elements in radial direction. The equations of the in-plane vibrations at disturbed state are also formulated using 1-dimensional finite elements in radial direction along the number of nodal diameters. They are expressed as in functions of the radial displacements at the steady state and the disturbed displacements about the steady state. In-plane static deformation modes of the annular disk are used as the interpolation functions of 1-dimensional finite elements in radial direction. The natural vibrations of an annular disk with different boundary conditions are analyzed by using the presented model and the 3-dimensional finite element model to verify accuracy of the presented equations of motion. Its results are compared and discussed.