Triangular, tetrahedral, and wedge elements

The library of solid elements in Abaqus includes first- and second-order triangles, tetrahedra, and wedge elements for planar, axisymmetric, and three-dimensional analysis.

The following topics are discussed:

Interpolation
Integration

Interpolation

The interpolation is defined in terms of the element coordinates g, h, and r shown in Figure 1. Since Abaqus is a Lagrangian code for most applications, these are also material coordinates. They each span a range from 0 to 1 in an element but satisfy the constraint that $g + h \leq 1$ for triangles and wedges and $g + h + r \leq 1$ for tetrahedra. The node numbering convention used in Abaqus for these elements is also shown in Figure 1. Corner nodes are numbered first, and then the midside nodes for second-order elements. The interpolation functions are as follows.

First-order triangle (3 nodes):

u = (1 - g - h) u_{1} + g u_{2} + h u_{3}

Second-order triangle (6 nodes):

\begin{matrix} u = & 2 (\frac{1}{2} - g - h) (1 - g - h) u_{1} + 2 g (g - \frac{1}{2}) u_{2} + 2 h (h - \frac{1}{2}) u_{3} \\ + 4 g (1 - g - h) u_{4} + 4 g h u_{5} + 4 h (1 - g - h) u_{6} \end{matrix}

First-order tetrahedron (4 nodes):

u = (1 - g - h - r) u_{1} + g u_{2} + h u_{3} + r u_{4}

Second-order tetrahedron (10 nodes):

\begin{matrix} u = & (2 (1 - g - h - r) - 1) (1 - g - h - r) u_{1} + (2 g - 1) g u_{2} + (2 h - 1) h u_{3} \\ + (2 r - 1) r u_{4} + 4 (1 - g - h - r) g u_{5} + 4 g h u_{6} + 4 (1 - g - h - r) h u_{7} \\ + 4 (1 - g - h - r) r u_{8} + 4 g r u_{9} + 4 h r u_{10} \end{matrix}

Figure 1. Isoparametric master elements.

First-order wedge (6 nodes):

\begin{matrix} u = & \frac{1}{2} (1 - g - h) (1 - r) u_{1} + \frac{1}{2} g (1 - r) u_{2} + \frac{1}{2} h (1 - r) u_{3} \\ + \frac{1}{2} (1 - g - h) (1 + r) u_{4} + \frac{1}{2} g (1 + r) u_{5} + \frac{1}{2} h (1 + r) u_{6} \end{matrix}

Second-order wedge (15 nodes):

\begin{matrix} u = & \frac{1}{2} ((1 - g - h) (2 (1 - g - h) - 1) (1 - r) - (1 - g - h) (1 - r^{2})) u_{1} \\ + \frac{1}{2} (g (2 g - 1) (1 - r) - g (1 - r^{2})) u_{2} + \frac{1}{2} (h (2 h - 1) (1 - r) - h (1 - r^{2})) u_{3} \\ + \frac{1}{2} ((1 - g - h) (2 (1 - g - h) - 1) (1 + r) - (1 - g - h) (1 - r^{2})) u_{4} \\ + \frac{1}{2} (g (2 g - 1) (1 + r) - g (1 - r^{2})) u_{5} + \frac{1}{2} (h (2 h - 1) (1 + r) - h (1 - r^{2})) u_{6} \\ + 2 (1 - g - h) g (1 - r) u_{7} + 2 g h (1 - r) u_{8} + 2 h (1 - g - h) (1 - r) u_{9} \\ + 2 (1 - g - h) g (1 + r) u_{10} + 2 g h (1 + r) u_{11} + 2 h (1 - g - h) (1 + r) u_{12} \\ + (1 - g - h) (1 - r^{2}) u_{13} + g (1 - r^{2}) u_{14} + h (1 - r^{2}) u_{15} \end{matrix}

Second-order variable 15–18 node wedge (assuming all 18 nodes are defined):

\begin{matrix} u = & (\frac{1}{2} ((1 - g - h) (2 (1 - g - h) - 1) (1 - r) - (1 - g - h) (1 - r^{2})) + \frac{1}{4} (N_{16} + N_{18})) u_{1} \\ + (\frac{1}{2} (g (2 g - 1) (1 - r) - g (1 - r^{2})) + \frac{1}{4} (N_{16} + N_{17})) u_{2} \\ + (\frac{1}{2} (h (2 h - 1) (1 - r) - h (1 - r^{2})) + \frac{1}{4} (N_{17} + N_{18})) u_{3} \\ + (\frac{1}{2} ((1 - g - h) (2 (1 - g - h) - 1) (1 + r) - (1 - g - h) (1 - r^{2})) + \frac{1}{4} (N_{16} + N_{18})) u_{4} \\ + (\frac{1}{2} (g (2 g - 1) (1 + r) - g (1 - r^{2})) + \frac{1}{4} (N_{16} + N_{17})) u_{5} \\ + (\frac{1}{2} (h (2 h - 1) (1 + r) - h (1 - r^{2})) + \frac{1}{4} (N_{17} + N_{18})) u_{6} \\ + (2 (1 - g - h) g (1 - r) - \frac{1}{2} N_{16}) u_{7} + (2 g h (1 - r) - \frac{1}{2} N_{17}) u_{8} \\ + (2 h (1 - g - h) (1 - r) - \frac{1}{2} N_{18}) u_{9} + (2 (1 - g - h) g (1 + r) + \frac{1}{2} N_{16}) u_{10} \\ + (2 g h (1 + r) + \frac{1}{2} N_{17}) u_{11} + (2 h (1 - g - h) (1 + r) + \frac{1}{2} N_{18}) u_{12} \\ + ((1 - g - h) (1 - r^{2}) - \frac{1}{2} (N_{16} + N_{18})) u_{13} + (g (1 - r^{2}) - \frac{1}{2} (N_{16} + N_{17})) u_{14} \\ + (h (1 - r^{2}) - \frac{1}{2} (N_{17} + N_{18})) u_{15} + + N_{16} u_{16} + N_{17} u_{17} + N_{18} u_{18}, \end{matrix}

where

\begin{matrix} N_{16} = 4 g (1 - g - h) (1 - r^{2}) \end{matrix}

\begin{matrix} N_{17} = 4 g h (1 - r^{2}) \end{matrix}

\begin{matrix} N_{18} = 4 h (1 - g - h) (1 - r^{2}) . \end{matrix}

Integration

The first-order triangle and tetrahedron are constant stress elements and use a single integration point for the stiffness calculation when used in stress/displacement applications. A lumped mass matrix is used for both elements, with the total mass divided equally over the nodes. For heat transfer applications a three-point integration scheme is used for the conductivity and heat capacity matrices of the first-order triangle, with the integration points midway between the vertices and the centroid of the element; and a four-point integration scheme is used for the first-order tetrahedron. Distributed loads are integrated with two and three points for first-order triangles and tetrahedrons, respectively.

The three-point scheme is also used for the stiffness of the second-order triangle when it is used in stress/displacement applications. The mass matrix is integrated with a six-point scheme that integrates fourth-order polynomials exactly (Cowper, 1973). Distributed loads are integrated using three points. The heat transfer versions of the element use the six-point scheme for the conductivity and heat capacity matrices.

For stress/displacement applications the second-order tetrahedron uses 4 integration points for its stiffness matrix and 15 integration points for its consistent mass matrix. For heat transfer applications the conductivity and heat capacity matrices are integrated using 15 integration points. The first-order wedge uses 2 integration points for its stiffness matrix but 6 integration points for its lumped mass matrix. The second-order wedge uses 9 integration points for its stiffness matrix but 18 integration points for its consistent mass matrix. The integration schemes used for the second-order tetrahedra and wedge elements can be found in Stroud (1971).