Abstract:We present a structure-preserving scheme based on a recently-proposed mixed formulation for incompressible hyperelasticity formulated in principal stretches. Although there exist Hamiltonians introduced for quasi-incompressible elastodynamics based on different variational formulations, the one in the fully incompressible regime has yet been identified in the literature. The adopted mixed formulation naturally provides a new Hamiltonian for fully incompressible elastodynamics. Invoking the discrete gradient formula, we are able to design fully-discrete schemes that preserve the Hamiltonian and momenta. The scaled mid-point formula, another popular option for constructing algorithmic stresses, is analyzed and demonstrated to be non-robust numerically. The generalized Taylor-Hood element based on the spline technology conveniently provides a higher-order, robust, and inf-sup stable spatial discretization option for finite strain analysis. To enhance the element performance in volume conservation, the grad-div stabilization, a technique initially developed in computational fluid dynamics, is introduced here for elastodynamics. It is shown that the stabilization term does not impose additional restrictions for the algorithmic stress to respect the invariants, leading to an energy-decaying and momentum-conserving fully discrete scheme. A set of numerical examples is provided to justify the claimed properties. The grad-div stabilization is found to enhance the discrete mass conservation effectively. Furthermore, in contrast to conventional algorithms based on Cardano's formula and perturbation techniques, the spectral decomposition algorithm developed by Scherzinger and Dohrmann is robust and accurate to ensure the discrete conservation laws and is thus recommended for stretch-based material modeling.

Automatic Consistency Integrations for the Stress Update of J(2) Elastoplastic Rate Constitutive Equations with Combined Hardening

An explicit integration algorithm for the unified rate-dependent elastoplastic damage model applied to high-rate dynamics

IMPLICIT INTEGRATION ALGORITHM FOR NUMERICAL SOLUTION OF PARTIALLY SATURATED SOILS WITH AN ENHANCED GENERALIZED PLASTICITY CONSTITUTIVE MODEL

A trust region algorithm for finite-strain plasticity with strongly coupled hardening

A robust stress update algorithm for elastoplastic models without analytical derivation of the consistent tangent operator and loading/unloading estimation

A multiplicative finite strain crystal plasticity formulation based on additive elastic corrector rates: Theory and numerical implementation

A return mapping algorithm based on the hyper dual step derivative approximation for elastoplastic models

Semi-Implicit Algorithm for Elastoplastic Damage Models Involving Energy Integration

Characteristics of elasto-brittle-plastic softening constitutive theory and its numerical calculation

A structure-preserving integrator for incompressible finite elastodynamics based on a grad-div stabilized mixed formulation with particular emphasis on stretch-based material models

Consistent Tangent Moduli for Multi-Yield-Surface J(2) Plasticity Model

Development of Elasto-Plastic Material Models in A Natural Coordinate System

An implicit stress update algorithm for the plastic nonlocal damage model of concrete

Implicit Numerical Integration of Highly Nonlinear Plasticity Models

Explicit, Comprehensive Modeling of Multi-Axial Finite Strain Pseudo-Elastic Smas Up to Failure

Implicit Damping Iterative Algorithm to Solve Elastoplastic Static and Dynamic Equations.

m‐AGC tangent visco‐plastic operator with hardening/softening, and application to the visco‐plastic relaxation analysis of stable and unstable problems using fracture‐based geomechanical interfaces

Stress rate and the lagrangian formulation of the finite-strain plasticity for a von mises kinematic hardening model

Choice of the stress integration scheme for accurate large-deformation finite element analysis

A large deformation theory for rate-dependent elastic–plastic materials with combined isotropic and kinematic hardening

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening