• 소성∙가공
• /
• 제3권2호
• /
• pp.156-166
• /
• 1994

An arbitrary Lagrangian-Eulerian (ALE) finite element method has been developed. The finite element formation is derived and implemented for rigid-viscoplastic materials. The developed computer program is applied to the analysis of axisymmetric square die extrusion, which has many difficulties with updated Lagrangian approach. The results are compared with those from updated Largrangian approach. The results are compared with those from updated Lagrangian finite element program. Updating scheme of time dependent variables and mesh control are also examined.

• Structural Engineering and Mechanics
• /
• 제51권4호
• /
• pp.601-625
• /
• 2014

Based on continuum mechanics and the principle of virtual displacements, incremental total Lagrangian formulation (T.L.) and incremental updated Lagrangian formulation (U.L.) were presented. Both T.L. and U.L. considered the large displacement stiffness matrix, which was modified to be symmetrical matrix. According to the incremental updated Lagrangian formulation, small strain, large displacement, finite rotation of three dimensional Timoshenko fiber beam element tangent stiffness matrix was developed. Considering large displacement and finite rotation, a new type of tangent stiffness matrix of the beam element was developed. According to the basic assumption of plane section, the displacement field of an arbitrary fiber was presented in terms of nodal displacement of centroid of cross-area. In addition, shear deformation effect was taken account. Furthermore, a nonlinear finite element method program has been developed and several examples were tested to demonstrate the accuracy and generality of the three dimensional beam element.

• 전산구조공학
• /
• 제10권2호
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

• 대한기계학회논문집
• /
• 제13권5호
• /
• pp.820-830
• /
• 1989

본 논문에서는 타이어의 각 부분의 물성치 계산을 위한 식을 유한요소법에 적용할 수 있도록 제안하였다. 이 식은 강철 코드의 굽힙효과를 고려 하였으며, 특히 각 요소에서 전단변형이 일어나는 동안의 굽힘효과를 고려하였다. 유한요소 공식화는 가상일의 원리에 의하여 평형 방정식으로부터 유도하였고, Updated refer- ence coordinate에 대해 증분해석을 적용하여 Updated Lagrangian공식화를 하였다. 그리고 차량하중에 의하여 타이어가 노면에 접지될때의 응력상태를 게산할 수 있도록 접촉문제 공식화를 유한요소 공식화에 첨가 하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2002년도 추계 학술대회논문집
• /
• pp.47-52
• /
• 2002

We present an improved Lagrangian vortex method in 2-D incompressible unsteady viscous flows, which is based on a mesh-free integral approach of the velocity-vorticity formulation. Vorticity fields are represented by discrete vortex blobs that are updated by the Lagrangian vorticity transport with the particle strength exchange scheme. Velocity fields are expressed in a form of the Helmholtz decomposition, which are calculated by a fast algorithm of the Biot-Savart integration with a smoothed kernel and by a well-established panel method. No-slip condition is enforced through viscous diffusion of vorticity from a solid body into field. The vorticity flux is determined in such a way that spurious slip velocity vanishes. Through the comparison with the existing finite volume scheme for the transient vortical flows around an impulsively started cylinder at Reynolds number Re=550, we would obtain a more accurate scheme for vortex methods in complicated flows.

• Coupled systems mechanics
• /
• 제7권5호
• /
• pp.583-610
• /
• 2018

This paper discusses the issues associated with modeling frictional contact between solid bodies undergoing large deformations. The most common model for friction on contact interfaces in solid mechanics is the Coulomb friction model, in which two distinct responses are possible: stick and slip. Handling the transition between these two phases computationally has been a source of algorithmic instability, lack of convergence and non-unique solutions, particularly in the presence of large deformations. Most computational models for frictional contact have used penalty or updated Lagrangian approaches to enforce frictional contact conditions. These two approaches, however, present some computational challenges due to conditioning issues in penalty-type implementations and the iterative nature of the updated Lagrangian formulation, which, particularly in large simulations, may lead to relatively slow convergence. Alternatively, a plasticity-inspired implementation of frictional contact has been shown to handle the stick-slip conditions in a local, algorithmically efficient manner that substantially reduces computational cost and successfully avoids the issues of instability and lack of convergence often reported with other methods (Laursen and Simo 1993). The formulation of this approach, however, has been limited to the small deformations realm, a fact that severely limited its application to contact problems where large deformations are expected. In this paper, we present an algorithmically consistent formulation of this method that preserves its key advantages, while extending its application to the realm of large-deformation contact problems. We show that the method produces results similar to the augmented Lagrangian formulation at a reduced computational cost.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 추계학술대회논문집
• /
• pp.1253-1258
• /
• 2007

The media transport systems, such as printers, copy machines, facsimiles, ATMs, cameras, etc. have been widely used and being developed rapidly. In the development of those sheet-handling machineries, it is important to predict the static and dynamic behavior of the sheet with a high degree of reliability because the sheets are fed and stacked at such a high speed. Flexible media are very thin, light and flexible, so they behave in geometric nonlinearity with large displacement and large rotation but small strain. In the flexible media analysis, aerodynamic effect from the surrounding air must be included because any small force can make large deformation. In this paper, only the flexible media analysis is performed as early stage of analysis including aerodynamic effect. Through formulations and simulations for total Lagrangian(TL), updated Lagrangian (UL) and co-rotational(CR) method which are widely used for geometric nonlinear analysis, usefulness and reliability of each methods are investigated.

• 한국공간구조학회논문집
• /
• 제9권2호
• /
• pp.91-97
• /
• 2009

본 논문은 Total Lagrangian(TL)과 Updated Lagrangian(UL)을 사용하여 구조물 하부에 케이블을 연결하고, 그 케이블에 인장력을 가하여 구조물의 처짐을 제어하는 언더텐션 시스템에 있어서 케이블의 인장력에 따른 구조물의 처짐 거동의 변화와 각 부재력을 비교함으로서 언더텐션 시스템의 효율성을 검증하는데 있다. 일반적인 빔과 거더로 이루어진 구조와는 달리 언더텐션 시스템에서는 상부에서는 하중을 하부 케이블의 인장력을 이용하여 그 하중을 양 단부로 전달하게 된다. 언더텐션 시스템은 스트럿의 개수와 길이, 케이블의 초기 인장력의 크기에 따라서 그 효과가 다르게 나타날 수 있다. 또한 케이블이 설치된 장스팬 구조의 경우, 그 거동에 있어서 비선형성적 거동이 크게 나타나게 된다. 따라서 본 논문에서는 선형과 비선형해석 결과를 비교함으로써 비선형해석의 필요성에 대하여 논의하고자 한다.

• Structural Engineering and Mechanics
• /
• 제12권1호
• /
• pp.35-50
• /
• 2001

A numerical methodology is presented in this paper for the geometrically non-linear analysis of slender uni-dimensional structural elements under unilateral contact constraints. The finite element method together with an updated Lagrangian formulation is used to study the structural system. The unilateral constraints are imposed by tensionless supports or foundations. At each load step, in order to obtain the contact regions, the equilibrium equations are linearized and the contact problem is treated directly as a minimisation problem with inequality constraints, resulting in a linear complementarity problem (LCP). After the resulting LCP is solved by Lemke's pivoting algorithm, the contact regions are identified and the Newton-Raphson method is used together with path following methods to obtain the new contact forces and equilibrium configurations. The proposed methodology is illustrated by two examples and the results are compared with numerical and experimental results found in literature.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2003년도 봄 학술발표회 논문집
• /
• pp.235-242
• /
• 2003

The solids and the fluids in porous media have a relative velocity to each other. Due to physically and chemically different material properties and their relative velocity, the behavior of saturated porous media is extremely complicated. Thus, in order to describe and clarify the deformation behavior of saturated porous media, constitutive models for deformation of porous media coupling several effects such as flow of the fluids or thermodynanical change need to be developed in frame of Arbitrary Lagrangian Eulerian (ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian elements, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media that are considered for the behavior of the solids and the fluids. In this work, governing equations of porous media based on ALE description are obtained from governing equations in frame of updated Lagrangian description. Then, weak forms of these equations are derived using arbitrary weighting functions.