Shell Elements Overview
The historical shell element in Radioss is a simple bilinear Mindlin plate element coupled with a reduced integration scheme using one integration point. It is applicable in a reliable manner to both thin and moderately thick shells.
This element is very efficient if the spurious singular modes, called “hourglass modes”, which result from the reduced integration are stabilized.
The stabilization approach consists of providing additional stiffness so that the spurious singular modes are suppressed. Also, it offers the possibility of avoiding some locking problems. One of the first solutions was to generalize the formulation of Kosloff and Frazier 1 for brick element to shell element. It can be shown that the element produces accurate flexural response (thus, free from the membrane shear locking) and is equivalent to the incompatible model element of Wilson et al. 2 without the static condensation procedure. Taylor 3 extended this work to shell elements. Hughes and Liu 4 employed a similar approach and extended it to nonlinear problems.
Belytschko and Tsay 5 developed a stabilized flat element based on the projections developed by Flanagan and Belytschko 6. Its essential feature is that hourglass control is orthogonal to any linear field, thus preserving consistency. The stabilized stiffness is approached by a diagonal matrix and scaled by the perturbation parameters which are introduced as a regulator of the stiffness for nonlinear problems. The parameters are generally chosen to be as small as possible, so this approach is often called, perturbation stabilization.
- The parameters are user-inputs and are generally problem-dependent.
- Poor behavior with irregular geometries (in-plane, out-of-plane). The stabilized stiffness (or stabilized forces) is often evaluated based on a regular flat geometry, so they generally do not pass either the Patch-test or the Twisted beam test.
Belytschko 7 extended this perturbation stabilization to the 4-node shell element which has become widely used in explicit programs.
Belytschko 8 improved the poor behavior exhibited in the warped configuration by adding a coupling curvature-translation term, and a particular nodal projection for the transverse shear calculation analogous to that developed by Hughes and Tezduyar 6, and MacNeal 10. This element passes the Kirchhoff patch test and the Twisted Beam test, but it cannot be extended to a general 6 DOF element due to the particular projection.
Belytschko and Bachrach 11 used a new method called physical stabilization to overcome the first drawback of the quadrilateral plane element. This method consists of explicitly evaluating the stabilized stiffness with the help of 'assumed strains', so that no arbitrary parameters need to be prescribed. Engelmann and Whirley 12 have applied it to the 4-node shell element. An alternative way to evaluate the stabilized stiffness explicitly is given by Liu et al. 13 based on Hughes and Liu's 4-node selected reduced integration scheme element 4, in which the strain field is expressed explicitly in terms of natural coordinates by a Taylor-series expansion. A remarkable improvement in the one-point quadrature shell element with physical stabilization has been performed by Belytschko and Leviathan 14. The element performs superbly for both flat and warped elements especially in linear cases, even in comparison with a similar element under a full integration scheme, and is only 20% slower than the Belytschko and Tsay element. More recently, based on Belytschko and Leviathan's element, Zhu and Zacharia 15 implemented the drilling rotation DOF in their one-point quadrature shell element; the drilling rotation is independently interpolated by the Allman function 16 based on Hughes and Brezzi's 17 mixed variational formulation.
The physical stabilization with assumed strain method seems to offer a rational way of developing a cost effective shell element with a reduced integration scheme. The use of the assumed strains based on the mixed variational principles, is powerful, not only in avoiding the locking problems (volumetric locking, membrane shear locking, as in Belytschko and Bindeman 18; transverse shear locking, as in Dvorkin and Bathe 19), but also in providing a new way to compute stiffness. However, as highlighted by Stolarski et al. 20, assumed strain elements generally do not have rigorous foundations; there is almost no constraint for the independent assumed strains interpolation. Therefore, a sound theoretical understanding and numerous tests are needed in order to prove the legitimacy of the assumed strain elements.
The greatest uncertainty of the one-point quadrature shell elements with physical stabilization is with respect to the nonlinear problems. All of these elements with physical stabilization mentioned above rely on the assumptions that the spin and the material properties are constant within the element. The first assumption is necessary to ensure the objectivity principle in geometrical nonlinear problems. The second was adapted in order to extend the explicit evaluation of stabilized stiffness for elastic problems to the physical nonlinear problems. It is found that the second assumption leads to a theoretical contradiction in the case of an elastoplastic problem (a classic physical nonlinear problem), and results in poor behavior in case of certain crash computations.
Zeng and Combescure 21 have proposed an improved 4-node shell element named QPPS with one-point quadrature based on the physical stabilization which is valid for the whole range of its applications (see Shell Formulations). The formulation is based largely on that of Belytschko and Leviathan.
Based on the QPPS element, Zeng and Winkelmuller have developed a new improved element named QEPH which is integrated in Radioss 44 version (see 3-Node Shell Elements).