Numerical analysis of a locking-free primal hybrid method for linear elasticity with $H(\mathrm{div})$-conforming stress recovery

TL;DR

This paper introduces a locking-free primal hybrid finite element method for linear elasticity that ensures optimal convergence and robust stress recovery, particularly effective for nearly incompressible materials.

Contribution

The work develops a general analysis for stable approximation spaces and proposes a stress recovery strategy, enhancing robustness and accuracy in linear elasticity simulations.

Findings

Achieves optimal convergence orders on triangular and quadrilateral meshes.

Provides a stress recovery method that yields $H(\mathrm{div})$-conforming, locally equilibrated stress fields.

Demonstrates robustness in nearly incompressible elasticity problems.

Abstract

In this work, we study a primal hybrid finite element method for the approximation of linear elasticity problems, posed in terms of displacement, an auxiliary pressure field, and a Lagrange multiplier related to the traction. We develop a general analysis for the existence and uniqueness of the solution for the discrete problem, which is applied to the construction of stable approximation spaces on triangular and quadrilateral meshes. The use of these spaces lead to optimal convergence orders, resulting in a locking-free method capable of providing robust approximations for nearly incompressible problems. Finally, we propose a strategy for recovering the stress field from the hybrid solution by solving element-wise sub-problems. The resulting stress approximation is $H(\mathrm{div})$-conforming, locally equilibrated, weakly symmetric, and robust to locking.