A Locking-free modified conforming FEM for planar elasticity (2407.06831v2)

Abstract: Due to the divergence-instability, the accuracy of low-order conforming finite element methods (FEMs) for nearly incompressible elasticity equations deteriorates as the Lam\'e parameter $\lambda\to\infty$, or equivalently as the Poisson ratio $\nu\to1/2$. This effect is known as \itshape{locking} or \itshape{non-robustness}. For the piecewise linear case, the error in the ${\bf L}^2$-norm of the standard Galerkin conforming FEM is bounded by $C\lambda h^2$, resulting in poor accuracy for practical values of $h$ if $\lambda$ is sufficiently large. In this paper, we show that for 2D problems the locking phenomenon can be controlled by replacing $\lambda$ with $\lambda_h=\lambda/(1+\lambda h/L)$ in the stiffness matrix, where $L$ is the diameter of the body $\Omega$. We prove that with this modification, the error in the ${\bf L}^2$-norm is bounded by $Ch$ for a constant $C$ that does not depend on $\lambda$. Numerical experiments confirm this convergence behaviour and show that, for practical meshes, our method is more accurate than the standard method if the material is nearly incompressible. Our analysis also shows that the error in the ${\bf H}^1$-norm is bounded by $Ch^{1/2}$, but our numerical experiments suggest that this bound is not sharp.