Four-field mixed finite elements for incompressible nonlinear elasticity

TL;DR

This paper introduces a stable four-field mixed finite element method for incompressible nonlinear elasticity that avoids stabilization, achieves high accuracy, and demonstrates superior robustness in 2D and 3D simulations.

Contribution

It proposes a novel four-field mixed formulation with discontinuous displacement, eliminating the need for stabilization and providing proven stability and error estimates.

Findings

Achieves optimal or super convergence rates in 2D and 3D.

Demonstrates enhanced robustness over existing methods like CSFEM.

Provides an efficient postprocessing technique for accurate displacement recovery.

Abstract

We present a stable finite element method for incompressible nonlinear elasticity based on a four-field mixed formulation involving the displacement, displacement gradient, first Piola--Kirchhoff stress and pressure. Unlike existing four-field mixed formulations, such as the compatible strain mixed finite element method (CSFEM), the proposed approach employs a discontinuous displacement field and requires no stabilisation in either 2D or 3D. A Newton--Raphson linearisation is derived and finite element pairs satisfying the relevant inf-sup conditions are identified. To recover accurate continuous displacement fields, an efficient postprocessing technique is further introduced. We establish the well-posedness of the linearised continuous problem together with a priori error estimates for the discrete formulation. Extensive numerical experiments in both 2D and 3D demonstrate optimal or even super convergence rates and enhanced robustness, particularly in 3D where CSFEM typically requires stabilisation.