• Structural Engineering and Mechanics
• /
• v.24 no.2
• /
• pp.227-245
• /
• 2006

The performance of a three dimensional non-linear finite element model for hyperelastic material considering the effect of compressibility is studied by analyzing rubber blocks under different modes of deformation. It includes simple tension, pure shear, simple shear, pure bending and a mixed mode combining compression, shear and bending. The compressibility of the hyperelastic material is represented in the strain energy function. The nonlinear formulation is based on updated Lagrangian (UL) technique. The displacement model is implemented with a twenty node brick element having u, ${\nu}$ and w as the degrees of freedom at each node. The results obtained by the present numerical model are compared with the analytical solutions available for the basic modes of deformation where the agreement between the results is found to be satisfactory. In this context some new results are generated for future references since the number of available results on the present problem is not sufficient enough.

• Journal of the korean Society of Automotive Engineers
• /
• v.12 no.2
• /
• pp.29-42
• /
• 1990

Three dimensional frame structures with such nonlinearities as large displacements, medium rotations, plastic hinges and local defects are efficiently analyzed by introducing the nonlinear rotational spring. Formulations are based on the incremental updated Lagrangian descriptions and the virtual work principle, Axial displacement and twisted angle in beam elements are interpolated linearly, while bending displacements are approximated by the Hermite polynomials. The modified are length method is used as a solution method. The moment-angle of rotation relationship obtained analytically or experimentally can be easily incorporated into the solution procedure. Several examples tested show that the present method can be used efficiently in analyzing nonlinear frame structures with plastic hinges or local defect.

• Journal of Hydrospace Technology
• /
• v.1 no.1
• /
• pp.111-124
• /
• 1995

An improved analysis model for material nonlinearity induced by elasto-plastic deformation and damage including a large strain response was proposed. The elasto-plastic-damage constitutive model based on the continuum damage mechanics approach was adopted to overcome limitations of the conventional plastic analysis theory. It can manage the anisotropic tonsorial damage evolved during the time-independent plastic deformation process of materials. Updated Lagrangian finite element formulation for elasto-plastic damage coupling problems including large deformation, large rotation and large strain problems was completed to develop a numerical model which can predict all kinds of structural nonlinearities and damage rationally. Finally a finite element analysis code for two-dimensional plane problems was developed and the applicability and validity of the numerical model was investigated through some numerical examples. Calculations showed reasonable results in both geometrical nonlinear problems due to large deformation and material nonlinearity including the damage effect.

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

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.8
• /
• pp.2113-2122
• /
• 1994

The 3-D FEM analysis for simulating the stamping operation of planar anisotropic sheet metals with arbitrarily-shaped tools is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The consistent full set of governing relations, comprising equilibrium equation and mesh-normal geometric constraints, is appropriately linearized. The linear triangular elements are used for depicting the formed sheet, based on membrane approximation. Barlat's non-quadratic anisotropic yield criterion(strain-rate potential) is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and non-quadratic function parameter. The planar anisotropic finite element formulation is tested with the numerical simulations of the stamping of an automotive hood inner panel and the drawing of a hemispherical punch. The in-plane anisotropic effects on the formability of both mild steel and aluminum alloy sheet metals are examined.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.4
• /
• pp.513-528
• /
• 2013

A floating Oscillating Water Column (OWC) wave energy converter, a Backward Bent Duct Buoy (BBDB), was simulated using a state-of-the-art, two-dimensional, fully-nonlinear Numerical Wave Tank (NWT) technique. The hydrodynamic performance of the floating OWC device was evaluated in the time domain. The acceleration potential method, with a full-updated kernel matrix calculation associated with a mode decomposition scheme, was implemented to obtain accurate estimates of the hydrodynamic force and displacement of a freely floating BBDB. The developed NWT was based on the potential theory and the boundary element method with constant panels on the boundaries. The mixed Eulerian-Lagrangian (MEL) approach was employed to capture the nonlinear free surfaces inside the chamber that interacted with a pneumatic pressure, induced by the time-varying airflow velocity at the air duct. A special viscous damping was applied to the chamber free surface to represent the viscous energy loss due to the BBDB's shape and motions. The viscous damping coefficient was properly selected using a comparison of the experimental data. The calculated surface elevation, inside and outside the chamber, with a tuned viscous damping correlated reasonably well with the experimental data for various incident wave conditions. The conservation of the total wave energy in the computational domain was confirmed over the entire range of wave frequencies.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.5
• /
• pp.85-98
• /
• 1993

A method of analysis for the material and geometric nonlinear analysis of planar prestressed concrete cable-stayed bridges including the time-dependent effects due to load history, creep, shrinkage, aging of concrete and relaxation of prestress is described. The analysis procedure, based on the finite element method, is capable of predicting the response of these structures through elastic, cracking, inelastic and ultimate ranges. The nonlinear formulation for the description of motion is based on the updated Lagrangian approach. To account for the material nonlinearity, nonlinear stress-strain relationship and cracking of concrete, nonlinear stress-strain relationships of reinforcing steel, prestressing steel, and cable, including load reversal are given. Results from a numerical examples on ultimate analyses of cable-stayed bridges are presented to illustrate the analysis method.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1995.03a
• /
• pp.196-202
• /
• 1995

Drawbead formulation is the first process together with a binder wrap process in a sheet metal forming process. The purpose of a drawbead is to control the flow of the metal into the die in panel press forming. To simulate the drawbead formation process, an elasto-plastic finite element formulation is derived from the equilibrium equation an drelated boundary conditions considering the proper contact conditons. The developed finite element program is applied to drawbead formation in the plane strain condition. The simulation of drawbead formation produces the distribution fo stress and strain along the bead and the resultant elongation of the sheet in the cavity region with respect to various cavity dimensions of the sheet as well as the punch force of a drawbead and the amount of draw-in with respect to the stroke fo a drawbead. The numerical resutls provides the fundamental information as a boundary condition to analyze the complex binder wrap phenomena and panel press forming in simple way.

• Structural Engineering and Mechanics
• /
• v.56 no.6
• /
• pp.917-938
• /
• 2015

The nonlinear numerical analysis of the impact response of reinforced concrete/mortar beam incorporated with the updated Lagrangian method, namely the Smoothed Particle Hydrodynamics (SPH) is carried out in this study. The analysis includes the simulation of the effects of high mass low velocity impact load falling on beam structures. Three material models to describe the localized failure of structural elements are: (1) linear pressure-sensitive yield criteria (Drucker-Prager type) in the pre-peak regime for the concrete/mortar meanwhile, the shear strain energy criterion (Von Mises) is applied for the steel reinforcement (2) nonlinear hardening law by means of modified linear Drucker-Prager envelope by employing the plane cap surface to simulate the irreversible plastic behavior of concrete/mortar (3) implementation of linear and nonlinear softening in tension and compression regions, respectively, to express the complex behavior of concrete material during short time loading condition. Validation upon existing experimental test results is conducted, from which the impact behavior of concrete beams are best described using the SPH model adopting an average velocity and erosion algorithm, where instability in terms of numerical fragmentation is reduced considerably.