• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

• Steel and Composite Structures
• /
• v.21 no.2
• /
• pp.373-394
• /
• 2016

The postbuckling behavior of laminated composite plates and shells, subjected to various shear loadings, is presented, using a modified 8-ANS method. The finite element, based on a modified first-order shear deformation theory, is further improved by the combined use of assumed natural strain method. We analyze the influence of the shell element with the various location and number of enhanced membrane and shear interpolation. Using the assumed natural strain method with proper interpolation functions, the present shell element generates neither membrane nor shear locking behavior even when full integration is used in the formulation. The effects of various types of lay-ups, materials and number of layers on initial buckling and postbuckling response of the laminated composite plates and shells for various shear loading have been discussed. In addition, the effect of direction of shear load on the postbuckling behavior is studied. Numerical results and comparisons of the present results with those found in the literature for typical benchmark problems involving symmetric cross-ply laminated composites are found to be excellent and show the validity of the developed finite element model. The study is relevant to the simulation of barrels, pipes, wing surfaces, aircrafts, rockets and missile structures subjected to intense complex loading.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.540-547
• /
• 1994

Contact and frictional locking conditions and the effect of shielding due to contact at the facet, which could be represented by the difference in energy release rate, as a function of phase angle of loading are analyzed in this study for the case of interfacial cracks by assuming single crack-kink model. The analysis of contact effects on interfacial fracture resistance shows that relative shielding increases as the shear component was increased, which indicates a qualitative agreement with the previous experimental results.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1993.10a
• /
• pp.195-199
• /
• 1993

In the analysis of metal forming processes by the finite element method, there are many numerical instabilities such as element locking, hourglass mode, shear locking. These instabilities may have a bad effect upon accuracy and convergence. The present work is concerned with improvement of stability and efficiency in two dimensional rigid-plastic finite element method using various type of elements and numerical integration schemes. AS metal forming examples, upsetting and backward extrusion are taken for comparison among the methods : various element types and numerical integration schemes. comparison is made in terms of stability and efficiency. As a result, it has been shown that the finite element computation is stabilized from the viewpoint of computational time, convergency, and numerical instability.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.367-371
• /
• 2001

In the shape optimization based on the finite element method, the accuracy of finite element analysis of a given structure is important to determine the final shape. In case of a bending dominant problem, finite element solutions by the full integration scheme are not reliable because of the locking phenomenon. Furthermore, in the process of shape optimization, the mesh distortion is large due to the change of the structure outline: therefore, we cannot guarantee the accurate result unless the finite element itself is accurate. We approach to more accurate shape optimization to diminish these inaccuracies by improving the existing finite element. The shape optimization using the modified finite element is applied to a two-dimensional simple beam. Results show that the modified finite element have improved the optimization results.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.11a
• /
• pp.173-176
• /
• 2000

Finite element model based on the Naghdi's shell theory in the general tensor-based form is formulated in the present study. Partial mixed variational functional for assumed strain is formulated in order to avoid the severe locking troubles known as transverse shear and membrane locking. The proposed assumed strain element in general tensor Naghdi's shell model provides very accurate solutions for thin shells in benchmark problems. In additions, linear elastic constitutive equations are given in the general curvilinear coordinate system including anisotropic layered structures. Thus laminated composited shell structures are easily analyzed in the present formulation.

• Structural Engineering and Mechanics
• /
• v.73 no.1
• /
• pp.1-16
• /
• 2020

Finite elements based on the partition of unity (PU) approximation have powerful capabilities for p-adaptivity and solutions with high smoothness without remeshing of the domain. Recently, the PU approximation was successfully applied to the three-node shell finite element, properly eliminating transverse shear locking and showing excellent convergence properties and solution accuracy. However, the enrichment with the PU approximation results in a significant increase in the number of degrees of freedom; therefore, it requires greater computational cost, thus making it less suitable for practical engineering. To circumvent this disadvantage, we propose a new strategy to decrease the total number of degrees of freedom in the existing PU-based shell element, without loss of optimal convergence and accuracy. To alleviate the locking phenomenon, we use the method of mixed interpolation of tensorial components and perform convergence studies to show the accuracy and capability of the proposed shell element. The excellent performances of the new shell elements are illustrated in three benchmark problems.

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

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.55-61
• /
• 2002

In the present study, a first order shear deformable Loop-subdivision triangular element which can handle transverse shear deformation of moderately thick shell and composite laminated or sandwich shells are developed. The developed element is more general than the previous one based on classical shell theory, since it includes the effect of transverse shell deformation and has standard five degrees of freedom per node. The quartic box spline function is employed as the interpolation basis function. Numerical examples for the benchmark static shell problems are analyzed to assess the performance of the developed subdivision shell element and locking trouble.

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