• Journal of the Earthquake Engineering Society of Korea
• /
• v.16 no.4
• /
• pp.37-52
• /
• 2012

In this study, a new 4-node general shell element with 6 DOFs per node is presented. Drilling rotational degrees of freedom are introduced by the variational principle with an independent rotation field. In formulation of the element, substitute transverse shear strain fields are used to avoid shear locking, while four nonconforming modes are applied in the in-plane displacement fields as a remedy for membrane locking. In addition, a direct modification method for nonconforming modes is employed in the numerical implementation of nonconforming modes to represent constant strain states. A 9-points integration rule is adopted for volume integration in the computation of the element stiffness matrix. With the combined use of these techniques, the developed shell element has no spurious zero energy modes, and can represent a constant strain state. Several numerical tests are carried out to evaluate the performance of the new element developed. The test results show that the behavior of the elements is satisfactory.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.3
• /
• pp.279-293
• /
• 2004

The analysis of linear static and free vibration problems of isotropic and laminated composite plates and shells is performed by the improved 9-node shell element with the new strain displacement relationship. In that relationship, the effect of new additional terms between the bending strain and displacement has been investigated in the warping problem. Natural co ordinate based strains, stresses and constitutive equations are used. The assumed natural strain method is used to alleviate both membrane and shear locking behavior from the element. The Lanczos method is employed in the calculation of the eigenvalues of laminated composite structures and the Gauss integration rule is adopted to evaluate the mass matrix. The numerical examples are compared with the analytical solutions to validate the current formulation and the results presented could be useful for the understanding of the behaviour of laminates under free vibration conditions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.9 s.180
• /
• pp.2292-2301
• /
• 2000

We propose an accurate and efficient estimation method of transverse shear stresses for analysis and design of laminated composite structures by 4-node quadrilateral degenerated shell elements. To get proper distributions of transverse shear stresses in each layer, we use 3-dimensional equilibrium equations instead of constitutive equations with shear correction factors which vary diversely according to the shapes of shell sections. Three dimensional equilibrium equations are integrated through the thickness direction with complete polynomial membrane stress fields, which are recovered by REP (Recovery by Equilibrium in Patches) recovery method. The 4-node quadrilateral degenerated shell element used in this paper has drilling degrees of freedom and shear stresses derived from assumed strain fields that are set up at natural coordinate systems. The numerical results demonstrate that the proposed estimation method attains reasonable accuracy and efficiency compared with other methods and FE analysis using 4-node degenerated shell elements.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.4
• /
• pp.107-114
• /
• 2004

In this paper, a simple improved 8-node finite element for the finite element analysis of fibrous composite plates is presented by using the direct modification. We drive explicit expressions of shape functions for the 8-node element with bilinear element geometry, which is modified so that it can represent any quadratic fields in Cartesian coordinates. The refined first-order shear deformation theory is proposed, which results in parabolic through-thickness distribution of the transverse shear strains and stresses from the formulation based on the third-order shear deformation theory. It eliminates the need for shear correction factors in the first-order theory. This finite element is further improved by combined use of assumed strain, modified shape function, and refined first-order theory. To show the effectiveness of our simple modification on the 8-node finite elements, numerical studies are carried out the static, buckling and free vibration analysis of fibrous composite plates.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.2
• /
• pp.25-34
• /
• 2003

In the present paper, formulation of a nine-node assumed strain shell element is modified and extended for analysis of actuator embedded/attached structures. The shell element can alleviate locking and has sic DOFs per node by discarding assumption of no thickness change. In modeling of the physicalquantities, we have assumed linear strain field through the whole thickness direction. The electric and mechanical quantities have been coupled through the constitutive equations. Unlike typical shell element, the present shell element allows thickness change. Thus, three-dimensional piezoelasticity can be accurately simulated. Base on the formulation, a finite element program is generated and the code is validated by solving numerical examples. The results from the present work are well agreed with those from other references.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.2982-2994
• /
• 1993

Low-velocity impact responses of composite laminates are investigated using the finite element method based on various theories. In two-dimensional nonlinear analysis, a displacement field considering higher order shear deformation and large deflection of the laminate is assumed and a finite element formulation is developed using a C$^{o}$-continuous 9-node plate element. Also, three-dimensional linear analysis based on the infinitesimal strain-displacement assumptions is performed using 8-node brick elements with incompatible modes. A modified Hertzian contact law is incorporated into the finite element program to evaluate the impact force. In the time integration, the Newmark constant acceleration algorithm is used in conjuction with successive iterations within each time step. Numerical results from static analysis as well as the impact response analysis are presented including impact force histories, deflections, strains in the laminate. Impact responses according to two typical low-velocity impact conditions are compared each other.

• Structural Engineering and Mechanics
• /
• v.15 no.1
• /
• pp.83-114
• /
• 2003

A new simple 3-node triangular flat-shell element with standard nodal DOF (6 DOF per node) is proposed for the linear and geometrically nonlinear analysis of very thin to thick plate and shell structures. The formulation of element GT9 (Long and Xu 1994), a generalized conforming membrane element with rigid rotational freedoms, is employed as the membrane component of the new shell element. Both one-point reduced integration scheme and a corresponding stabilization matrix are adopted for avoiding membrane locking and hourglass phenomenon. The bending component of the new element comes from a new generalized conforming Kirchhoff-Mindlin plate element TSL-T9, which is derived in this paper based on semiLoof constrains and rational shear interpolation. Thus the convergence can be guaranteed and no shear locking will happen. Furthermore, a simple hybrid procedure is suggested to improve the stress solutions, and the Updated Lagrangian formulae are also established for the geometrically nonlinear problems. Numerical results with solutions, which are solved by some other recent element models and the models in the commercial finite element software ABAQUS, are presented. They show that the proposed element, denoted as GMST18, exhibits excellent and better performance for the analysis of thin-think plates and shells in both linear and geometrically nonlinear problems.

• Journal of Korean Society of Steel Construction
• /
• v.23 no.4
• /
• pp.455-464
• /
• 2011

In this study, the buckling analysis of the vertically curved stiffened web plate was conducted through finite-element analysis, using an eight- and nine-node flat shell element with a substitute shear strain field. To investigate the buckling behavior of the vertically curved web plate with a longitudinal or vertical stiffener under in-plane moment loading, parametric studies were conducted for the variation of the width (b) and ratio of the bending stiffness of the stiffener to that of the plate (${\gamma}=EI/bD$). The static behavior of the vertically curved web plate without a stiffener was also investigated, and then the buckling abilities of the longitudinal and vertical stiffeners were compared under moment loading.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.3
• /
• pp.11-21
• /
• 2008

A linear and non-linear analysis for plates and shells with arbitrary edge supports subjected to various loading was presented. The 9-node ANS(Assumed Natural Strain) hell element was employed and the spring element, which could express an arbitrary edge support using the six degrees of freedom, was introduced. For the application of his analysis, the plates and shells with various edge supports were analyzed, and the ending behavior with these edge supports were obtained accurately. For these edge supports, particularly elastic edge support was simulated by six springs and reasonable results were obtained. The results show that the present method can be widely used to analyze the bending behavior of plates and shells with arbitrary edge conditions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.