• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.5
• /
• pp.141-153
• /
• 1999

A finite element procedure for the analysis of rubber-like hyperelastic material is developed. The volumetric incompressiblity conditions of the rubber deformation is included in the formulation by using penalty method. In this paper, the behavior of the rubber deformation is represented by hyperelastic constitutive relations based on a generalized Mooney-Rivlin model. The principle of virtual work is used to derive nonlinear finite element equation for the large displacement problem and presented in total-Lagrangian description. The finite element procedure using analytic differentiation resulted in very close solution to the result of the well known commercial packages NISAII AND ABAQUS. Numerical tests show that the results from the numerical differentiation method coincide very well with those from the analytic method and the well known commercial packages in static analysis. The convergency of rubber usingν iteration method is also discussed.

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

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.10a
• /
• pp.249-253
• /
• 1997

A numerical method for evaluating the equilibrium path of cylindrical shell subject to axial load and eccentrically axial load is presented. The effects of both material and geometric nonlinearities were also considered in the analysis. The nonlinear formulation was based on the total Lagrangian description and nonlinear equtions were solved by the Newton-Raphson method with load increment procedures. Degenerate shell elements with layered approach were employed for the analysis. The elasto-plastic deformation can be found in several examples and a large eccentricity of the axial load reduces the stress level at the time of the local buckling of the pipe considerably.

• Computational Structural Engineering
• /
• v.2 no.2
• /
• pp.53-62
• /
• 1989

In recent times, the subterranean caverns for storing crude oils and oil products are increasingly needed. The elasto-VIScoplastic DYNamic finite element analysis program(VISDYN) has been developed in order to investigate dynamic responses of the storage cavity. And validity of the program is studied through a numerical example. Mohr-Coulomb yield criterion is adopted and associated flow rule is assumed. Geometrically nonlinear behaviour is taken into account using a total Lagrangian formulation. In dynamic deformation reponses, the difference between the steady state displacements and the unsteady state ones by the static analysis can be neglected.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.335-342
• /
• 2002

A continuum-based design sensitivity analysis (DSA) method fur non-shape problems is developed for geometrically nonlinear elastic structures. The non-shape problem is characterized by the design variables that are not associated with the domain of system like sizing, material property, loading, and so on. Total Lagrangian formulation with the Green-Lagrange strain and the second Piola-Kirchhoff stress is employed to describe the geometrically nonlinear structures. The spatial domain is discretized using the 4-node isoparametric plane stress/strain elements. The resulting nonlinear system is solved using the Newton-Raphson iterative method. To take advantage of the derived analytical sensitivity In topology optimization, a fast and efficient design sensitivity analysis method, adjoint variable method, is employed and the material property of each element is selected as non-shape design variable. Combining the design sensitivity analysis method and a gradient-based design optimization algorithm, an automated design optimization method is developed. The comparison of the analytical sensitivity with the finite difference results shows excellent agreement. Also application to the topology design optimization problem suggests a very good insight for the layout design.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1995.05a
• /
• pp.141-153
• /
• 1995

An incremental Total Lagrangian Formulation is implemented for the finite element analysis of laminated composite pressure vessel with consideration of the material and geometric nonlinearities. For large displacements/large rotations due to geometric nonlinearities, the incremental equations are derived using a quadratic approximation for the increment of the reference vectors in terms of the nodal rotation increments. This approach leads to a complete tangent stiffness matrix. For material nonlinearity, the analysis is performed by using the piecewise linear method, taking account of the nonlinear shear stress-strain relation. The results of numerical tests include the large deflection behavior of the selected composite shell problem. When compared with the previous analysis, tile results are in good agreement with them. As a practical example, filament wound pressure vessel is analyzed with consideration of the geometrically and materially nonlinearity. The numerical results agree fairly well with the existing experimental results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.4
• /
• pp.375-387
• /
• 2004

A new finite element model will be presented to analyze the nonlinear behavior of RC beams and slabs strengthened by a patch repair. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on hardening rule, crushing condition, plate-end debonding strength model and so on. The Gauss-Lobatto numerical quadrature is applied to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version nonlinear finite element model is demonstrated through the load-deflection curves, the ultimate loads, and the failure modes of RC beams or slabs bonded with steel plates or FRP plates compared with available result of experiment and other numerical methods.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

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

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.7 s.100
• /
• pp.836-842
• /
• 2005

Nonlinear dynamic characteristics of active piezolaminated plates are investigated with respect to the thermopiezoelastic behaviors. For largely deformed structures with small strain, the incremental total Lagrangian formulation is presented based on the virtual work principles. A multi-field layer-wise finite shell element is proposed for assuring high accuracy and non-linearity of displacement, electric and thermal fields. For dynamic consideration of thermopiezoelastic snap-through phenomena, the implicit Newmark's scheme with the Newton-Raphson iteration is implemented for the transient response of various piezolaminated models with symmetric or eccentric active layers. The bifurcate thermal buckling of symmetric structural models is first investigated and the characteristics of piezoelectric active responses are studied for finding snap-through piezoelectric potentials and the load-path tracking map. The thermoelastic stable and unstable postbuckling, thermopiezoelastic snap-through phenomena with several attractors are proved using the nonlinear time responses for various initial conditions and damping loss factors. Present results show that thermopiezoelastic snap-through phenomena can result in the difficulty of buckling and postbuckling control of intelligent structures.