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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.3
• /
• pp.195-202
• /
• 2012

Two nonlinear frame elements taking into account geometric nonlinearity is presented and compared based on the Lagrangian co-rotational formulation. The first frame element is believed to be geometrically-exact because not only tangent stiffness matrices is exactly evaluated including stiffness matrices due to initial deformation but also total member forces are directly determined from total deformations in the deformed state. Particularly two exact tangent stiffness matrices based on total Lagrangian and updated Lagrangian formulation, respectively, are verified to be identical. In the second frame element, the deformed curved shape is regarded as the polygon and current flexural deformations in iteration process are neglected in evaluating tangent stiffness matrices and total member forces. Two numerical examples are given to demonstrate the accuracy and the good performance of the first frame element compared with the second element. Furthermore it is shown that the first frame element can be used in tracing nonlinear behaviors of cable members.

• 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.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.485-490
• /
• 2015

This paper presents development of a geometrically nonlinear triangular planar element including rotational degrees of freedom, based on the co-rotational(CR) formulation. The CR formulation is one of the efficient geometrically nonlinear formulations and it is based on the assumptions on small strain and large rotation. In this paper, modified CR formulation is suggested for the developemnt of a triangular planar element. The present development is validated regarding the static and time transient problems. The present results are compared with the results predicted by the previous researchers and those obtained by the existing commercial software.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.123-131
• /
• 2005

According to ow previous study, we confirmed That the Petrov-Galerkin natural element method(PG-NEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method(BG-NEM). This paper is an extension of PG-NEM to two-dimensional geometrically nonlinear problem. For the analysis, a linearized total Lagrangian formulation is approximated with the PS-NEM. At every load step, the grid points ate updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates The large deformation problem.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.788-791
• /
• 2005

Nodal displacements are referred to the Initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian formulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid fer structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacements and traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One closed loop structure undergoing large deformations is analyzed to demonstrate the efficiency and validity of the proposed method.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2001.05a
• /
• pp.163-167
• /
• 2001

A two-dimensional progressive failure analysis method is presented for the strength characterization of the composite joints under pin loading. The eight-nodes laminated she]1 element is utilized based on the updated Lagrangian formulation. The criteria by Yamada-Sun, Tsai-Wu, and the maximum stress are used for the failure estimation. The stiffness of failed layer is degraded by the complete unloading method. No factor depending on test is included in the finite element analysis except for the material strength and stiffness. Total 20 plate specimens with and without hole are tested to validate the finite element prediction. The Tsai-Wu failure criterion most conservatively estimates the strength of laminate, and the maximum stress criterion yields the highest strength because it does not consider the coupling of the failure modes. The strength by Yamada-Sun method neglecting the matrix failure effect are located between other two methods and shows best agreement with test result for laminate with hole.

• Selected Papers of The Society of Naval Architects of Korea
• /
• v.1 no.1
• /
• pp.91-100
• /
• 1993

A displacement-based finite element method Is presented for the geometrically nonlinear analysis of eccentrically stiffened plates. A nonlinear degenerated shell element and a nonlinear degenerated eccentric isoparametric beam (isobeam) element are formulated on the basis of Total Agrangian and Updated Lagrangian descriptions. In the formulation of the isobeam element, some additional local decrees of freedom are implementd to describe the stiffener's local plate buckling modes. Therefore this element can be effectively employed to model the eccentric stiffener with fewer D.O.F's than the case of a degenerated shell element. Some detailed buckling and nonlinear analyses of an eccentrically stiffened plate are performed to estimate the critical buckling loads and the post buckling behaviors including the local plate buckling of the stiffeners discretized with the degenerated shell elements and the isobeam elements. The critical buckling loads are found to be higher than the analytical plate buckling load but lower than Euler buckling load of the corresponding column, i.e, buckling strength requirements of the Classification Societies for the stiffened plates.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.4 s.235
• /
• pp.534-539
• /
• 2005

Nodal displacements are referred to the initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian furmulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid for structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacement sand traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One open loop and one closed loop structure undergoing large deformations are analyzed to demonstrate the efficiency and validity of the proposed method.