2 Elements
2 Elements
2 Elements
1.Introduction
Abaqus provides a wide range of elements for different geometries and analysis types. An element is
characterized by the following:
Family: Continuum, Shell, Membrane, Rigid, Beam, Truss, etc.
Number of Nodes: Depends on element shape and its order of interpolation
Degrees of Freedom per Node: Depends on the solution field of the analysis (Displacement,
Rotation, Temperature, etc.) and the order of modeling space (1D, 2D and 3D).
Formulation: Applicable to some of the element families (eg small- and finite-strain shells)
otherwise most elements provide both automatically
Integration: Reduced and full integration
Fig 2- Example of Abaqus elements and naming convention for elements in Abaqus.
3. Element Selection
Element selection generally requires some knowledge on the theory and application of finite element
analysis (FEA). For different problems and analysis types there are some general (rule-of-thumb)
guidelines on the element selection. User should generally know that in Abaqus:
All elements include a variety of element-based loads, for example body forces (eg gravity) and
surface pressure loads on solid and shell elements, and forces distributed on a line for beam and
shell elements.
All elements are suitable for geometrically non-linear analysis (except small-strain shell).
There is generally no restriction on the type of material used for an element.
Most elements can be used for both Standard and Explicit solvers.
3.1.Solid Elements
The solid (or continuum) elements in Abaqus can be used for linear analysis and for complex nonlinear
analyses involving contact, plasticity, and large deformations. They are available for stress, heat
transfer, acoustic, coupled thermal-stress, coupled pore fluid-stress, piezoelectric, magnetostatic,
electromagnetic, and coupled thermal-electrical analyses. An Abaqus user should know that
Triangular and quadrilateral elements are commonly used for 2D modeling of solids.
For 3D modeling tetrahedral and hexahedral elements can be used.
Generally the use of quadrilateral and hexahedral elements is recommended for modeling of
continuum solids.
Solid element in Abaqus can provide first-order (linear) and second-order (quadratic)
interpolation.
3.1.1 Quadrilateral and Hexahedra Elements
Two-dimensional quadrilateral elements include CPE4R (plane strain), CPS4R (plane stress) and
CAX4R (axisymmetric) and three-dimensional hexahedral element include C3D8R (8-node trilinear
brick) and C3D20R (20-node quadratic brick). Application of quads and hex elements are generally
recommended whenever possible however these elements are geometrically less versatile and difficult
to be used for complicated geometries. It should be noted that by effective use of different domain
decompositions techniques (eg partitioning) it is possible use these elements even in complex
geometries.
3.1.2 Triangular and Tetrahedral elements
Tetrahedral elements provide geometrically versatility and most of automatic meshing algorithms
utilize these types of elements. User of Abaqus should generally know that
First-order (linear) tetrahedral elements (CPE3R, CPS3R, CAX3R and C3D4) are stiff and they
should not be used for structural analysis (except if an extremely fine mesh is used).
Modified second-order tetrahedral and triangles (CPE6M, CPS6M, CAX6M, and C3D10M) are
effective alternatives to linear types when the application requires use of these elements (like in
complex geometries). These elements can be used with automatic mesh generation and are
robust for large deformations and contact problems.
A good mesh of hexahedral elements (C3D8R) usually provides a solution of equivalent
accuracy at less cost.
3.2. Shell and Membrane Elements
Shell elements are developed based on the shell theory that approximates a thin 3D continuum (small
thickness compared to lateral dimensions) using a 2D formulation. Depending on the ration between
thickness and lateral directions dimensions different theories have been developed. The use of shell and
membrane elements requires a good level of understanding about the theories of shells and membranes
and proper selection of these elements should be consulted with experiences users. Abaqus offers
conventional shell elements (2D geometry and thickness is defined by section property) and continuum
shell elements (3D geometry).
They use of these elements is also recommended for composite laminate structures because they
capture through-thickness response. These elements have both wedge-shaped (SC6R) and hexahedral
(SC8R) geometries. The choice of thickness direction for hexahedral shell elements can be ambiguous
since any of the six faces can be the bottom face. By default Abaqus uses standard nodal connectivity
convention to define thickness direction, however non-default direction can be specifined using the
stack direction parameter on the shell section definition.
Trusses are rods that only carry tensile and compressive loads. Abaqus also provide truss elements
(T3D2).
3.4 Rigid Bodies/Elements
A rigid body is a collection of points that their motion is governed by the motion of a single reference
point. This means a part of body that can not deform can be considered as a rigid body. Rigid bodies
are very commonly used in finite element analysis because they are computational very efficient (only
require 6 degrees of freedom). Abaqus allows for definition of analytical rigid bodies and rigid
elements. Analytical rigid bodies are more commonly used because of their several advantages
including having potential to reduce noise in the solution, computational efficiency, and ease of
definition. Analytical rigid bodies have some disadvantages for example they can not be used to define
a general 3D rigid geometry. Three types of analytical rigid surfaces are available in Abaqus
(Segments, Cylinder, and Revolution).
uses a lower order integration (only one Gauss-point or more dependent on the geometry of element
and order of interpolation) to calculate the element matrices for example C3D20R only uses 8
integration points compared to C3D20 which uses 27 integration points therefore its computation is
almost 3.5 less costly.
Fig 7- Integration points used in reduced- (left) and full (right) integration in quadrilateral elements.
4.2 Hourglassing
The numerical algorithms used for the calculations and integration of the element matrices during
assembly step may influence how element behaves. Reduced-order integration (primarily used by
Abaqus Explicit) allow for fast and cheap calculation of the element matrices but the used of this
integration algorithm comes with some drawbacks. Hourglassing can be a problem with first-order,
reduced-integration elements (CPS4R, CAX4R, C3D8R, etc.) in stress/displacement analyses. Since
the elements have only one integration point, it is possible for them to distort in such a way that the
strains calculated at the integration point are all zero, which, in turn, leads to uncontrolled distortion of
the mesh. First-order, reduced-integration elements in Abaqus include hourglass control, but they
should be used with reasonably fine meshes. Hourglassing can also be minimized by distributing point
loads and boundary conditions over a number of adjacent nodes.
In Abaqus/Standard the second-order reduced-integration elements, with the exception of the 27-node
C3D27R and C3D27RH elements, do not have the same difficulty and are recommended in all cases
when the solution is expected to be smooth. The C3D27R and C3D27RH elements have three
unconstrained, propagating hourglass modes when all 27 nodes are present. These elements should not
be used with all 27 nodes, unless they are sufficiently constrained through boundary conditions. Firstorder elements are recommended when large strains or very high strain gradients are expected.
4.3 Shear and Volumetric Locking
Fully integrated elements in Abaqus/Standard and Abaqus/Explicit do not hourglass but may suffer
from locking behavior: both shear and volumetric locking. Shear locking occurs in first-order, fully
integrated elements (CPS4, CPE4, C3D8, etc.) that are subjected to bending. The numerical
formulation of the elements gives rise to shear strains that do not really existthe so-called parasitic
shear. Therefore, these elements are too stiff in bending, in particular if the element length is of the
same order of magnitude as or greater than the wall thickness.
Volumetric locking occurs in fully integrated elements when the material behavior is (almost)
incompressible. Spurious pressure stresses develop at the integration points, causing an element to
behave too stiffly for deformations that should cause no volume changes. If materials are almost
incompressible (elastic-plastic materials for which the plastic strains are incompressible), second-order,
fully integrated elements start to develop volumetric locking when the plastic strains are on the order of
the elastic strains. However, the first-order, fully integrated quadrilaterals and hexahedra use selectively
reduced integration (reduced integration on the volumetric terms). Therefore, these elements do not
Use second-order, fully integrated elements close to stress concentrations to capture the severe
gradients in these regions. However, avoid these elements in regions of finite strain if the
material response is nearly incompressible.
Use first-order quadrilateral or hexahedral elements or the modified triangular and tetrahedral
elements for problems involving large distortions. If the mesh distortion is severe, use reducedintegration, first-order elements.
If the problem involves bending and large distortions, use a fine mesh of first-order, reducedintegration elements.
Hybrid elements must be used if the material is fully incompressible (except when using plane
stress elements). Hybrid elements should also be used in some cases with nearly incompressible
materials.
6. References
[1] Abaqus 6.13 Documentation