Unconditionally stable, second-order schemes for gradient-regularized, non-convex, finite-strain elasticity modeling martensitic phase transformations

TL;DR

This paper introduces two unconditionally stable, second-order time-integration schemes for simulating gradient-regularized, non-convex finite-strain elasticity in martensitic phase transformations, addressing numerical challenges like stiffness and microstructure resolution.

Contribution

The paper presents two novel unconditionally stable, second-order schemes specifically designed for gradient elasticity problems involving martensitic transformations, improving numerical stability and accuracy.

Findings

Schemes are unconditionally stable and second-order accurate.

Numerical examples demonstrate improved stability and microstructure resolution.

Each scheme offers distinct advantages in computational performance.

Abstract

In the setting of continuum elasticity martensitic phase transformations are characterized by a non-convex free energy density function that possesses multiple wells in strain space and includes higher-order gradient terms for regularization. Metastable martensitic microstructures, defined as solutions that are local minimizers of the total free energy, are of interest and are obtained as steady state solutions to the resulting transient formulation of Toupin's gradient elasticity at finite strain. This type of problem poses several numerical challenges including stiffness, the need for fine discretization to resolve microsstructures, and following solution branches. Stable and accurate time-integration schemes are essential to obtain meaningful solutions at reasonable computational cost. In this work we introduce two classes of unconditionally stable second-order time-integration schemes for gradient elasticity, each having relative advantages over the other. Numerical examples are shown highlighting these features.