The ABAQUS element library provides a complete geometric modeling capability. For this reason any combination of elements can be used to make up the model. Sometimes multi-point constraints are required for application of the necessary kinematic relations to form the model (for example, to model part of a shell surface with solid elements and part with shell elements or to model a pipe elbow with a mixture of beam and shell elements).
All elements use numerical integration to allow complete generality in material behavior. Shell and beam element properties can be defined as general section behaviors, or each cross-section of the element can be integrated numerically, so that nonlinear response can be tracked accurately when needed. A composite layered section can be specified, with different materials at different heights through the section. Some special elements (such as line springs) use an approximate analytical solution to model nonlinear behavior.
All of the elements in ABAQUS are formulated in a global Cartesian coordinate system except the axisymmetric elements, which are formulated in terms of r–z coordinates. In almost all elements, primary vector quantities (such as displacements and rotations ) are defined in terms of nodal values with scalar interpolation functions. For example, in elements with a two-dimensional topology the interpolation can be written as
All elements in ABAQUS are integrated numerically. Hence, the virtual work integral as described in Nonlinear solution methods in ABAQUS/Standard, Section 2.2.1, will be replaced by a summation:
The advantage of the reduced integration elements is that the strains and stresses are calculated at the locations that provide optimal accuracy, the so-called Barlow points (Barlow, 1976). A second advantage is that the reduced number of integration points decreases CPU time and storage requirements. The disadvantage is that the reduced integration procedure can admit deformation modes that cause no straining at the integration points. These zero-energy modes make the element rank-deficient and cause a phenomenon called “hourglassing,” where the zero energy mode starts propagating through the mesh, leading to inaccurate solutions. This problem is particularly severe in first-order quadrilaterals and hexahedra. To prevent these excessive deformations, an additional artificial stiffness is added to the element. In this so-called hourglass control procedure, a small artificial stiffness is associated with the zero-energy deformation modes. This procedure is used in many of the solid and shell elements in ABAQUS.
Most fully integrated solid elements are unsuitable for the analysis of (approximately) incompressible material behavior. The reason for this is that the material behavior forces the material to deform (approximately) without volume changes. Fully integrated solid element meshes, and in particular lower-order element meshes, do not allow such deformations (other than purely homogeneous deformation). For that reason ABAQUS uses “selectively reduced” integration in these elements: reduced integration is used for the volume strain and full integration for the deviatoric strains. As a consequence the lower-order elements give an acceptable performance for approximately incompressible behavior. For fully incompressible material behavior, another complication occurs: the bulk modulus and, hence, the stiffness matrix becomes infinitely large. For this case a mixed (hybrid) formulation is required, where the displacement field is augmented with a hydrostatic pressure field. In this formulation only the inverse of the bulk modulus appears, and, consequently, the contribution to the operator matrix vanishes. The hydrostatic pressure field plays the role of a Lagrange multiplier field enforcing the incompressibility constraints.
ABAQUS/Standard also provides elements for multifield problems. Examples are the pore pressure elements used for the analysis of porous solids with fluid diffusion, coupled temperature-displacement elements that couple heat transfer with stress analysis, and piezoelectric elements that couple electrical conduction with stress analysis. In these multifield elements the scalar variable (such as the temperature) is usually interpolated with different scalar functions as the displacement field; i.e.,