Reduced Integration and Hourglass Stabilization
COMSOL Multiphysics uses by default a full and consistent quadrature order when numerically evaluating integrals such as the virtual work. However, a full quadrature order does not necessarily lead to an optimal formulation of the resulting finite element matrices in structural mechanics, where higher-order terms can introduce locking and other detrimental effects. A remedy can be to reduce the quadrature order for selected parts of the virtual work, typically the strain energy density. When full integration is used, the quadrature order is twice the shape order of the displacement field. And when a reduced integration scheme is used, the quadrature order is reduced by half.
Using a reduced integration scheme not only affects the numerical integration of the weak contribution related to the displacement field, but for consistency it also changes the shape order of auxiliary dependent variables such as the out-of-plane strain in plane stress, the auxiliary pressure in a mixed formulation, and other auxiliary dependent variables added by nonlinear material models. This also applies to state variables used in material models such as plasticity, damage, and viscoelasticity. Lowering the order of auxiliary dependent variables and state variables significantly reduces the computation time due to the reduced assembly time.
Using reduced integration affects the stationary part of the virtual work; a full quadrature order is always applied to inertial term contributions in implicit solvers, and mass lumping techniques can be used for explicit dynamics analysis.
Furthermore, reduced integration is applied to weak contributions added by material models. Other contributions to the virtual work such as body loads and other boundary conditions use a full quadrature scheme.
Hourglass Stabilization
The hourglass stabilization method adds a stabilization term Wstb to the strain energy density and to the principle of virtual work.
The exact formulation of Wstb differs from case to case, but it is implemented so that it is inexpensive to evaluate. For example, the expression does not include evaluation of nonlinear materials such as plasticity.
Energy Sampling Method
When using reduced integration along with the energy sampling method, the stationary part of the virtual work is written as
The additional terms that include Wstb are integrated using different quadrature orders. The first and second integrals are evaluated using reduced integration, while the third integral is evaluated using full integration. Wstb is chosen such that the combined contribution of the second and third integrals evaluates to zero every time the reduced integration in the first integral is accurate, while it adds a non zero contribution in presence of hourglass modes.
Flanagan–Belytschko Method
The Flanagan–Belytschko method (Ref. 53) introduces an artificial stiffness to resist hourglass modes in quadrilateral and hexahedral elements. Hourglass modes are represented as a linear combination of nodal displacement basis.
The method is efficient for controlling hourglass modes without significantly over-stiffening the elements.
Hessian Method
The strain energy contribution Wstb is derived from the norm of the second gradient of the displacement field, also called the Hessian of the field.
See also Using Reduced Integration in the Structural Mechanics Modeling chapter.