A locking-free face-centred finite volume (FCFV) method for linear elasticity

TL;DR

This paper introduces a novel face-centred finite volume (FCFV) method for linear elasticity that is robust, locking-free, and accurate on distorted meshes, demonstrated through various benchmark tests in 2D and 3D.

Contribution

The paper presents a new FCFV method that enforces stress symmetry explicitly and computes displacement and stress fields element-wise without reconstruction.

Findings

Optimal convergence of displacement and stress fields.

Method is robust against mesh distortion and stretching.

Effective in complex geometries and classical benchmark problems.

Abstract

A face-centred finite volume (FCFV) method is proposed for the linear elasticity equation. The FCFV is a mixed hybrid formulation, featuring a system of first-order equations, that defines the unknowns on the faces (edges in two dimensions) of the mesh elements. The symmetry of the stress tensor is strongly enforced using the well-known Voigt notation and the displacement and stress fields inside each cell are obtained element-wise by means of explicit formulas. The resulting FCFV method is robust and locking-free in the nearly incompressible limit. Numerical experiments in two and three dimensions show optimal convergence of the displacement and the stress fields without any reconstruction. Moreover, the accuracy of the FCFV method is not sensitive to mesh distortion and stretching. Classical benchmark tests including Kirch's plate and Cook's membrane problems in two dimensions as well as three dimensional problems involving shear phenomenons, pressurised thin shells and complex geometries are presented to show the capability and potential of the proposed methodology.