• Structural Engineering and Mechanics
• /
• v.15 no.2
• /
• pp.199-223
• /
• 2003

The formulation of a non-linear shear deformable shell element is presented for the solution of stability problems of stiffened plates and shells. The formulation of the geometrical stiffness presented here is exactly defined on the midsurface and is efficient for analyzing stability problems of thick plates and shells by incorporating bending moment and transverse shear resultant force. As a result of the explicit integration of the tangent stiffness matrix, this formulation is computationally very efficient in incremental nonlinear analysis. The element is free of both membrane and shear locking behaviour by using the assumed strain method such that the element performs very well in the thin shells. By using six degrees of freedom per node, the present element can model stiffened plate and shell structures. The formulation includes large displacement effects and elasto-plastic material behaviour. The material is assumed to be isotropic and elasto-plastic obeying Von Mises's yield condition and its associated flow rules. The results showed good agreement with references and computational efficiency.

• Structural Engineering and Mechanics
• /
• v.18 no.1
• /
• pp.41-58
• /
• 2004

The postbuckling behaviour of thin shells has fascinated researchers because the theoretical prediction and their experimental verification are often different. In reality, shell panels possess small imperfections and these can cause large reduction in static buckling strength. This is more relevant in thin laminated composite shells. To study the postbuckling behaviour of thin, imperfect laminated composite shells using finite elements, explicit incremental or secant matrices have been presented in this paper. These incremental matrices which are derived using Marguerre's shallow shell theory can be used in combination with any thin plate/shell finite element (Classical Laminated Plate Theory - CLPT) and can be easily extended to the First Order Shear deformation Theory (FOST). The advantage of the present formulation is that it involves no numerical approximation in forming total potential energy of the shell during large deformations as opposed to earlier approximate formulations published in the literature. The initial imperfection in shells could be modeled by simply adjusting the ordinate of the shell forms. The present formulation is very easy to implement in any existing finite element codes. The secant matrices presented in this paper are shown to be very accurate in tracing the postbuckling behaviour of thin isotropic and laminated composite shells with general initial imperfections.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1992.03a
• /
• pp.201-212
• /
• 1992

The plane strain analysis for simulating the stretch/draw forming operation of arbitrarily-shaped tool profiles and arbitrarily draw-in conditions is introduced. An implicit, incremental, updated Lagrangian formulation is employed, introducing a rigid-viscoplastic constitutive equation. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshe without depending on the explicit spatial derivatives of tool surfaces. The FEM formulation is tested in the sections automotive inner panel and two-side draw-in. Not only the excellent agreement between measured and computed strains in the stretched section is obtained, but also the numerical stability of current formulation is verified in the two-side draw-in section.

• Structural Engineering and Mechanics
• /
• v.30 no.1
• /
• pp.119-129
• /
• 2008

A small strain and elastoplastic formulation of Polygonal Element Method (PEM) is developed for efficient analysis of elastoplastic solids. In this work, the polygonal elements are constructed based on traditional triangular finite meshes. The construction method of polygonal mesh can directly utilize the sophisticated triangularization algorithm and reduce the difficulty in generating polygonal elements. The Wachspress rational finite element basis function is used to construct the approximations of polygonal elements. The incremental variational form and a von Mises type model are used for non-linear elastoplastic analysis. Several small strain elastoplastic numerical examples are presented to verify the advantages and the accuracy of the numerical formulation.

• Transactions of Materials Processing
• /
• v.1 no.1
• /
• pp.95-103
• /
• 1992

The plane strain analysis for simulating the stretch/draw forming operation with an arbitrarily-shaped tool profile 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 linear line elements are used for depicting the formed sheet, based on membrane approximation. The FEM formulation is tested in the sections of automotive inner panel and two-side draw-in. Not only the excellent agreement between measured and computed strains is obtained in the stretched section, but also the numerical stability of formulation is verified in the draw-in section.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.3
• /
• pp.107-115
• /
• 1985

For the finite deflection analysis of plates with initial deflections subjected to uniaxial compression, the formulation of incremental finite strip method is made and has been incorporated into a computer program. A new in plane displacement function varying along the load: direction has been derived from the out-of-plane displacement function by considering the curvature of a plate. Either incremental load type analysis or incremental displacement type analysis may be selected to solve incremental equibrium equations in the program. The following results have been obtained: 1. Incremental displacement type analysis is superior to incremental load type analysis in that the former converges more rapidly than the latter. 2. The finite strip method using the new displacement function gives as accurate results as analytical method and other finite element methods.

• Transactions of the KSME C: Technology and Education
• /
• v.3 no.2
• /
• pp.107-115
• /
• 2015

The analysis of multi-flexible-body dynamics (MFBD) has been an important issue in the area of the computational dynamics. This technique has been developed and improved in RecurDyn solver. This paper reviews the formulation which is applied in the RecurDyn solver. Basically, in order to solve the multi-flexible-body dynamics problem, an incremental finite element formulation using a corotational procedure is used. In particular, in order to solve the rigid and flexible bodies together, a constraint equation between a rigid body and a flexible body is applied, in which a virtual body and a flexible body joint are introduced.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.25 no.9
• /
• pp.938-943
• /
• 2014

The scattered field by a gap/crack on the PEC surface of a large object having low-observable RCS cannot be negligible, but may not be analyzed by the known high-frequency technique. If the electrical width of the crack is very small, the crack can be modeled by an impedance strip, whose scattering formulation can be analytically obtained based on a low-frequency approximation. The scattering solution is formulated for the 2D strip and TE(Transverse Electric) or TM(Transverse Magnetic) wave incidence, from which a 3D ILDC(Incremental Length Diffraction Coefficients) can be extracted. Using the ILDC formulation, the scattering by any arbitrary shaped crack can be estimated. In this paper, an improved ILDC equations are proposed, which combine the known TE and TM solutions. The improved accuracy of the proposed solution is numerically verified.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.197-208
• /
• 2000

Choi et al./sup 1)/ presented the total Lagrangian formulation based upon the degenerated shell element. Geometrically correct formulation is developed by updating the direction of normal vectors and taking into account the second order rotation terms in the incremental displacement field. Assumed strain concept is adopted in order to overcome the shear locking phenomena and to eliminate the spurious zero energy mode. In this paper, for the ultimate strength analysis of stiffened shell structures considering effects of residual stresses, the return mapping algorithm based on the consistent elasto-plastic tangent modulus is applied to anisotropic shell structures. In addition, the load/displacement incremental scheme is adopted for non-linear F.E. analysis. Based on such methodology, the computer program is developed and numerical examples to demonstrate the accuracy and the effectiveness of the proposed shell element are presented and compared with the results in literatures.

• Structural Engineering and Mechanics
• /
• v.26 no.6
• /
• pp.725-739
• /
• 2007

In recent years there are many plate bending elements that emerged for solving both thin and thick plates. The main features of these elements are that they are based on mix formulation interpolation with discrete collocation constraints. These elements passed the patch test for mix formulation and performed well for linear analysis of thin and thick plates. In this paper a member of this family of elements, namely, the Discrete Reissner-Mindlin (DRM) is further extended and developed to analyze both thin and thick plates with geometric nonlinearity. The Von K$\acute{a}$rm$\acute{a}$n's large displacement plate theory based on Lagrangian coordinate system is used. The Hu-Washizu variational principle is employed to formulate the stiffness matrix of the geometrically Nonlinear Discrete Reissner-Mindlin (NDRM). An iterative-incremental procedure is implemented to solve the nonlinear equations. The element is then tested for plates with simply supported and clamped edges under uniformly distributed transverse loads. The results obtained using the geometrically NDRM element is then compared with the results of available analytical solutions. It has been observed that the NDRM results agreed well with the analytical solutions results. Therefore, it is concluded that the NDRM element is both reliable and efficient in analyzing thin and thick plates with geometric non-linearity.