The MAC scheme for linear elasticity in displacement-stress formulation on non-uniform staggered grids

TL;DR

This paper introduces a marker-and-cell finite difference scheme for 2D and 3D linear elasticity problems that ensures stability, avoids spurious oscillations, and achieves super-convergence, suitable for various elastic materials.

Contribution

The paper develops a novel staggered grid finite difference method for elasticity that guarantees stability, super-convergence, and is locking-free for all elastic material types.

Findings

The scheme is stable and super-convergent for displacement and stress.

It is locking-free for all Lame constants, suitable for compressible and nearly incompressible materials.

Numerical results confirm the theoretical accuracy and robustness.

Abstract

A marker-and-cell finite difference method is developed for solving the two dimensional and three dimensional linear elasticity in the displacement-stress formulation on staggered grids. The method employs a staggered grid arrangement, where the displacement components are approximated on the midpoints of cell edges, the normal stresses are defined at the cell centers, and the shear stresses are defined at the grid points. This structure ensures local conservation properties and avoids spurious oscillations in stress approximation. A rigorous mathematical analysis is presented, establishing the stability of the scheme and proving the second-order L2-error super-convergence for both displacement and stress. The proposed method is locking-free with respect to the Lame constant, making it suitable for both compressible and nearly incompressible elastic materials. Numerical experiments demonstrate the efficiency and robustness of the finite difference scheme, and the computed results show excellent agreement with the theoretical predictions.