A simple modification to mitigate locking in conforming FEM for nearly incompressible elasticity

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 λ∞, or equivalently as the Poisson ratio 1/2. This effect is known as locking or non-robustness. For the piecewise linear case, the error in the L2-norm of the standard Galerkin conforming FEM is bounded by~Cλ h2, resulting in poor accuracy for practical values of~h if λ is sufficiently large. In this short paper, we show that the locking phenomenon can be reduced by replacing λ with~λh=λμ/(μ+λ h/L)<λ in the stiffness matrix, where μ is the second Lam\'e parameter and L is the diameter of the body . We prove that with this modification, the error in the L2-norm is bounded by Ch for a constant C that does not depend on λ. Numerical experiments confirm this convergence behaviour and show that, for practical meshes, our method is more accurate than the standard method if λ is larger than about μ L/h. Our analysis also shows that the error in the H1-norm is bounded by Cλh1/2\,h, which improves the Cλ1/2\,h estimate for the case of conforming FEM.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…