20-Node Hexahedral Element (HEX20)
20-Node Hexahedral Element (HEX20): Theoretical Foundations
HEX20 — The Highest Accuracy 3D Element
Is HEX20 the most accurate 3D element in FEM?
It can be said to be the most accurate among practical elements. It is a 20-node quadratic hexahedral element, with 3 degrees of freedom per node, totaling 60 DOF. It has 2.5 times the DOF of HEX8 (24 DOF), but its accuracy is orders of magnitude higher.
Shape Functions
HEX20 has 8 vertex nodes and 12 edge mid-point nodes. The shape functions are of the serendipity type:
Vertex nodes:
Edge mid-point nodes (e.g., $\xi_i = 0$):
What is the Serendipity type? Is it different from the Lagrange type?
The Lagrange type (27-node HEX) has nodes on faces and inside, while the Serendipity type (20-node HEX) has mid-side nodes only on edges, with no nodes at face centers or inside. The Serendipity type has nearly equivalent accuracy to the Lagrange type with fewer DOF, making it overwhelmingly more common in practice.
Accuracy of HEX20
How accurate is HEX20?
Comparison with TET10:
| Accuracy for the same DOF count | HEX20 | TET10 |
|---|---|---|
| Displacement | Very High | High |
| Stress | Very High | High |
| Efficiency (Accuracy per DOF) | Highest | 50–70% of HEX20 |
- Displacement: $O(h^4)$ (HEX20) vs. $O(h^3)$ (TET10)
- Stress: $O(h^3)$ (HEX20) vs. $O(h^2)$ (TET10)
HEX20 converges one order faster. It achieves comparable accuracy even with a mesh twice as coarse...
Therefore, for geometries where HEX meshing is possible, HEX20 is the most efficient. The problem is the effort required for mesh generation. This trade-off between "accuracy efficiency vs. mesh generation cost" is the essence of choosing between HEX20 and TET10.
Integration Scheme
Integration for HEX20:
| Integration | Gauss Points | Characteristics |
|---|---|---|
| Full Integration (3×3×3) | 27 | Highest accuracy. Low risk of locking |
| Reduced Integration (2×2×2) | 8 | Avoids shear locking. Hourglassing is 1 mode |
Reduced integration for HEX20 is 2×2×2 = 8 points? That's the same number as full integration for HEX8.
Yes. Reduced integration for HEX20 (C3D20R) uses the same 8 integration points as full integration for HEX8, but delivers far higher accuracy. C3D20R is considered one of the most efficient elements in 3D FEM.
Does C3D20R have hourglassing issues?
Reduced integration for HEX20 has only one hourglass mode. This is extremely few compared to the 12 modes of HEX8R, posing almost no problem in practice. However, caution is needed for single-element patch tests.
When to Use HEX20
In what situations should HEX20 be used?
Conversely, in what situations is HEX20 disadvantageous?
Summary
Let me organize the theory of HEX20.
Key points:
- 20 nodes, Serendipity type, quadratic hexahedron — Highest accuracy among practical 3D elements
- C3D20R (reduced integration) is highly recommended — 8 integration points for maximum efficiency
- Fast convergence — One order faster than TET10
- Mesh generation is a challenge — Automatic HEX is difficult. Suitable for sweep/mapped meshing
- Not suitable for large deformation or explicit analysis — HEX8R or TET10 are more stable
So HEX20 is "the best element if you can create it," right?
Exactly. There is a wall of mesh generation between "the highest accuracy" and "if you can create it." For engineers who can overcome this wall, HEX20 is the most powerful weapon.
Theoretical Superiority of Quadratic Hexahedral Elements
The 20-node hexahedral element (Serendipity element) is a masterpiece of isoparametric formulation proposed by Ergatoudis et al. in 1966. Its shape functions, which include complete quadratic polynomials, provide nearly 8 times the stress accuracy within an element compared to linear elements. Adopted in structural analysis for space agencies since NASA's Apollo program in the 1970s, it was once called the "gold standard element."
Computational Methods for 20-Node Hexahedral Element (HEX20)
C3D20R — The Most Efficient 3D Element
Please explain why C3D20R (Abaqus's HEX20 reduced integration) is considered "the most efficient".
There are three reasons:
1. Highest accuracy/DOF ratio — 2–3 times more accurate than TET10 for the same DOF count
2. No shear locking — Avoided by reduced integration
3. Minimal volumetric locking — No problem up to about $\nu = 0.49$
So C3D20R (reduced integration) is recommended over C3D20 (full integration).
C3D20 uses 27 integration points, resulting in high computational cost, and can show signs of volumetric locking for materials with high $\nu$. C3D20R uses 8 integration points, is faster to compute, and has less locking. In practice, C3D20R is overwhelmingly more common.
Element Names by Solver
| Variation | Abaqus | Nastran | Ansys |
|---|---|---|---|
| Full Integration | C3D20 | CHEXA(20) | SOLID186(full) |
| Reduced Integration | C3D20R | — | SOLID186(red.) |
| Hybrid | C3D20H, C3D20RH | — | u-P formulation |
Does Nastran's CHEXA(20) not have a reduced integration option?
Nastran's CHEXA(20) uses 2×2×2 (reduced) integration by default. There is an option to switch to full integration, but the default is reduced. So, Nastran's CHEXA(20) ≈ Abaqus's C3D20R.
Handling of Mid-Side Nodes
HEX20's mid-side nodes are at the midpoints of edges, but for curved surfaces, they should be snapped to the CAD, right?