• Structural Engineering and Mechanics
• /
• v.12 no.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.

• Structural Engineering and Mechanics
• /
• v.14 no.2
• /
• pp.135-150
• /
• 2002

Numerical analysis of ultimate behaviour of thin bionics shells is treated in present paper. Interactive conditions in resonance and stability ultimate response are considered. Numerical treatment of nonlinear problems appearing is made using the updated Lagrangian formulation of motion. Each step of the iteration approaches the solution of linear problem and the feasibility of parallel processing FETM-technique with adaptive mesh refinement and substructuring for the analysis of ultimate action of thin bionics shells is established. Some numerical results are submitted in order to demonstrate the efficiency of the procedures suggested.

• Computational Structural Engineering
• /
• v.10 no.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.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• 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.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1995.10a
• /
• pp.150-156
• /
• 1995

In the finite element analysis of metal forming processes, the updated Lagrangian approach has been widely and effectively used to simulate the non-steady state problems. However some difficulties have arisen from abrupt flow change as in extrusion through square dies. In the present work, a ALE(arbitrary Lagrangian-Eulerian) finite element formulation for deformation analysis are presented for rigid viscoplastic materials. The developed finite element program is applied to the analysis of square die extrusion of a square section. The computational results are compared with those from the updated Lagrangian finite element analysis.

• Journal of Korean Society of Steel Construction
• /
• v.11 no.2 s.39
• /
• pp.153-165
• /
• 1999

In order to trace the lateral-torsional post-bucking behaviors of thin-walled space frames with non-symmetric cross sections, a geometrically non-linear finite element formulation is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for non-symmetric thin-walled cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, tangent stiffness matrices of thin-walled space frame element are derived by using the Hermition polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines and incremental member forces.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.205-216
• /
• 1997

In order to trace the lateral post-buckling behaviors of thin-wafled space frames, a geometrically nonlinear finite element formulation is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric thin-walled cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, tangent stiffness matrices of the thinwalled space frame element with 7 degrees of freedom including the restrained warping for each node are derived by using the Hermition polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformations. Finite element solutions for the spatial buckling and post-buckling analysis of thin-walled space frames are presented and compared with available solutions and other researcher's results.

• Computational Structural Engineering
• /
• v.10 no.1
• /
• pp.201-211
• /
• 1997

A clearly consistent finite element formulation for geometrically non-linear analysis of space frames is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, elastic and geometric stiffness matrices of the space frame element are derived by using the Hermitian polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformaions. Finite element solutions for the spatial buckling and post-buckling analysis of space frames are compared with available solutions and other researcher's results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.3
• /
• pp.29-38
• /
• 1990

A totally, new approach of Lagrangian formulation named 'Pseudo Lagrangian Formulation(PLF)' for large deformation analysis of continue and structures by the finite of element method has been presented, and the efficiency and accuracy of nonlinear analysis beam element formulated by PLF has been discussed by solving several numerical examples. In PLF, the deformation of a body is maeasured by assigning a nonphysical 'Pseudo' configuration as reference. The Lagrangian deformation and the finite element mapping of the traditonal Lagrangian approaches are then carried out directly at the same time, The result of numerical tests shows superior performance of PLF to the traditional Lagrangian methods, Applications of PLF to small and finite deformation problems indicate that PLF not only serves as an alternative but has certain implementational advantages over total or updated Lagrangian formulations.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.3
• /
• pp.25-37
• /
• 1992

A nonlinear finite element formulation which has the capability of handling very large geometrical changes is presented. The formulation is based on an updated material reference frame and hence true stress-strain test can be directly applied to properly characterize properties of materials which are subjected to very large deformation. For the large deformation, a consistent formulation based on the continuum mechanics approach is derived. The kinematics is referred to an updated material frame. Body equilibrium is also established in an updated geometry and the second Piola-Kirchhoff stress and the updated Lagrangian strain tensor are used in the formulation. Numerical examples for very large deformation of framed structures and plane solids are analyzed for verification purposes. The numerical solutions are obtained by an incremental numerical procedure. The importance of handing material properties properly is also demonstrated.