A Locking-Free Weak Galerkin Finite Element Method for Linear Elasticity Problems

09/11/2023

∙

by   Fuchang Huo, et al.

∙

0

∙

share

In this paper, we introduce and analyze a lowest-order locking-free weak Galerkin (WG) finite element scheme for the grad-div formulation of linear elasticity problems. The scheme uses linear functions in the interior of mesh elements and constants on edges (2D) or faces (3D), respectively, to approximate the displacement. An H(div)-conforming displacement reconstruction operator is employed to modify test functions in the right-hand side of the discrete form, in order to eliminate the dependence of the Lamé parameter λ in error estimates, i.e., making the scheme locking-free. The method works without requiring λ∇·𝐮_1 to be bounded. We prove optimal error estimates, independent of λ, in both the H^1-norm and the L^2-norm. Numerical experiments validate that the method is effective and locking-free.