Achieving Material Robustness via Symmetric Stress Finite Element Discretizations
For computational mechanics researchers using finite element methods, this work identifies a critical flaw in weak symmetry enforcement that can cause catastrophic errors in stress predictions for anisotropic materials.
The paper shows that weak enforcement of stress tensor symmetry in finite element discretizations can yield arbitrarily poor stress approximations for anisotropic materials, while strong enforcement ensures material-robust accuracy.
When discretizing symmetric stress tensors in variational problems arising in continuum mechanics, one has to choose how to enforce the symmetry of the stress tensor: (i) strongly by requiring the discrete tensors to be pointwise symmetric or (ii) weakly by introducing a Lagrange multiplier. For $H(\mathrm{div})$-conforming finite element discretizations of Hellinger--Reissner elasticity and velocity--stress formulations of incompressible flow, where symmetry of the Cauchy stress tensor is tied to the conservation of angular momentum, we show that this choice may substantially impact the accuracy of the numerical scheme. Through a series of benchmark problems featuring anisotropic constitutive laws inspired by fiber reinforced material, liquid crystal polymer networks, and polar fluids, we show that schemes enforcing symmetry weakly can yield arbitrarily poor stress approximations -- even for zero-stress configurations. However, schemes enforcing symmetry strongly deliver accurate stress approximations independently of the constitutive law, a property we term material robustness. We present a unifying theory that rigorously explains this behavior.