An efficient and accurate semi-implicit time integration scheme for dynamics in nearly- and fully-incompressible hyperelastic solids

TL;DR

This paper introduces two semi-implicit time integration schemes for large-scale dynamics in nearly- and fully-incompressible hyperelastic solids, demonstrating their stability, accuracy, and specific stability properties, while highlighting challenges in volume preservation.

Contribution

The paper develops and empirically verifies two second-order accurate semi-implicit schemes, analyzing their stability and effectiveness for nonlinear elasticity problems.

Findings

FEBDF2's time step size scales inversely with shear wave speed

Both schemes are second-order accurate and stable

Semi-implicit schemes have difficulty preserving volume globally

Abstract

The choice of numerical integrator in approximating solutions to dynamic partial differential equations depends on the smallest time-scale of the problem at hand. Large-scale deformations in elastic solids contain both shear waves and bulk waves, the latter of which can travel infinitely fast in incompressible materials. Explicit schemes, which are favored for their efficiency in resolving low-speed dynamics, are bound by time step size restrictions that inversely scale with the fastest wave speed. Implicit schemes can enable larger time step sizes regardless of the wave speeds present, though they are much more computationally expensive. Semi-implicit methods, which are more stable than explicit methods and more efficient than implicit methods, are emerging in the literature, though their applicability to nonlinear elasticity is not extensively studied. In this research, we develop and investigate the functionality of two time integration schemes for the resolution of large-scale dynamics in nearly- and fully-incompressible hyperelastic solids: a Modified Semi-implicit Backward Differentiation Formula integrator (MSBDF2) and a forward Euler / Semi-implicit Backward Differentiation Formula Runge-Kutta integrator (FEBDF2). We prove and empirically verify second order accuracy for both schemes. The stability properties of both methods are derived and numerically verified. We find FEBDF2 has a maximum time step size that inversely scales with the shear wave speed and is unaffected by the bulk wave speed -- the desired stability property of a semi-implicit scheme. Finally, we empirically determine that semi-implicit schemes struggle to preserve volume globally when using nonlinear incompressibility conditions, even under temporal and spatial refinement.