5.1 Element geometry

Two types of shell elements are available in ABAQUS: conventional shell elements and continuum shell elements. Conventional shell elements discretize a reference surface by defining the element's planar dimensions, its surface normal, and its initial curvature. The nodes of a conventional shell element, however, do not define the shell thickness; the thickness is defined through section properties. Continuum shell elements, on the other hand, resemble three-dimensional solid elements in that they discretize an entire three-dimensional body yet are formulated so that their kinematic and constitutive behavior is similar to conventional shell elements. Continuum shell elements are more accurate in contact modeling than conventional shell elements, since they employ two-sided contact taking into account changes in thickness. For thin shell applications, however, conventional shell elements provide superior performance.

In this manual only conventional shell elements are discussed. Henceforth, we will refer to them simply as “shell elements.” For more information on continuum shell elements, see Shell elements: overview, Section 15.6.1 of the ABAQUS Analysis User's Manual.


5.1.1 Shell thickness and section points

You define the shell thickness using either the *SHELL SECTION or *SHELL GENERAL SECTION option.

If you use the *SHELL SECTION option, ABAQUS uses numerical integration to calculate the stresses and strains independently at each section point (integration point) through the thickness of the shell, thus allowing nonlinear material behavior. For example, an elastic-plastic shell may yield at the outer section points while remaining elastic at the inner section points. The location of the single integration point in an S4R (4-node, reduced integration) element and the configuration of the section points through the shell thickness are shown in Figure 5–1.

Figure 5–1 Configuration of section points in a numerically integrated shell.

You can specify any odd number of section points through the shell thickness with the *SHELL SECTION option. By default, ABAQUS uses five section points through the thickness of a homogeneous shell, which is sufficient for most nonlinear design problems. However, you should use more section points in some complicated simulations, especially when you anticipate reversed plastic bending (nine is normally sufficient in this case). For linear problems three section points provide exact integration through the thickness. However, the *SHELL GENERAL SECTION option is more efficient for linear elastic shells.

If you use the *SHELL GENERAL SECTION option, the material behavior must be linear elastic, as the stiffness of the cross-section is calculated only once at the beginning of the simulation. In this case, all calculations are done in terms of the resultant forces and moments across the entire cross-section. If you request stress or strain output, ABAQUS provides default output for the bottom surface, the midplane, and the top surface.


5.1.2 Shell normals and shell surfaces

The connectivity of the shell element defines the direction of the positive normal, as shown in Figure 5–2.

Figure 5–2 Positive normals for shells.

For axisymmetric shell elements the positive normal direction is defined by a 90° counterclockwise rotation from the direction going from node 1 to node 2. For three-dimensional shell elements the positive normal is given by the right-hand rule going around the nodes in the order in which they appear in the element definition.

The “top” surface of a shell is the surface in the positive normal direction and is called the SPOS face for contact definition. The “bottom” surface is in the negative direction along the normal and is called the SNEG face for contact definition. Normals should be consistent among adjoining shell elements.

The positive normal direction defines the convention for pressure load application and output of quantities that vary through the shell thickness. A positive pressure load applied to a shell element produces a load that acts in the direction of the positive normal. (The pressure load convention for shell elements is opposite to that for continuum elements.)


5.1.3 Initial shell curvature

Most shell elements in ABAQUS are formulated as true curved shell elements and, therefore, require accurate calculation of the initial surface curvature. ABAQUS automatically calculates the surface normals at the nodes of every shell element to estimate the initial curvature of the shell. The surface normal at each node is determined using a fairly elaborate algorithm, which is discussed in detail in Defining the initial geometry of conventional shell elements, Section 15.6.3 of the ABAQUS Analysis User's Manual.

With a coarse mesh as shown in Figure 5–3, ABAQUS may determine several independent surface normals at the same node for adjoining elements. Physically, multiple normals at a single node mean that there is a fold line between the elements sharing the node. While it is possible that you intend to model such a structure, it is more likely that you intend to model a smoothly curved shell; ABAQUS will try to smooth the shell by creating an averaged normal at a node.

Figure 5–3 Effect of mesh refinement on the nodal surface normals.

The basic smoothing algorithm used is as follows: if the normals at a node for each shell element attached to the node are within 20° of each other, the normals will be averaged. The averaged normal will be used at that node for all elements attached to the node. If ABAQUS cannot smooth the shell, a warning message is issued in the data (.dat) file.

You may override the default algorithm. To introduce fold lines into a curved shell or to model a curved shell with a coarse mesh, use the *NODE and *NORMAL options to define the normals manually. With the *NODE option you specify the surface normal at a node as the 4th, 5th, and 6th values on the data line, following the nodal coordinates. A normal you define with *NODE is the normal used for all elements sharing that node, unless *NORMAL is also used. Use the *NORMAL option to specify a normal at a node for selected elements only. Normals defined with *NORMAL override normals defined with *NODE. See Defining the initial geometry of conventional shell elements, Section 15.6.3 of the ABAQUS Analysis User's Manual, for further details.


5.1.4 Reference surface offsets

The reference surface of the shell is defined by the shell element's nodes and normal definitions. When modeling with shell elements, the reference surface is typically coincident with the shell's midsurface. However, many situations arise in which it is more convenient to define the reference surface as offset from the shell's midsurface. For example, surfaces created in CAD packages usually represent either the top or bottom surface of the shell body. In this case it may be easier to define the reference surface to be coincident with the CAD surface and, therefore, offset from the shell's midsurface.

Shell offsets can also be used to define a more precise surface geometry for contact problems where shell thickness is important. Shell offsets can be useful when modeling a shell with continuously varying thickness. In this case defining the nodes at the shell midplane can be difficult. If one surface is smooth while the other is rough, as in some aircraft structures, it is easiest to use shell offsets to define the nodes at the smooth surface.

Offsets can be introduced by using the OFFSET parameter on the *SHELL SECTION and *SHELL GENERAL SECTION options. The offset value is defined as a fraction of the shell thickness measured from the shell's midsurface to the shell's reference surface, as shown in Figure 5–4.

Figure 5–4 Schematic of shell offset for an offset value of 0.5.