Abstract:In deformable simulation, an important computing task is to calculate the gradient and derivative of the strain energy function in order to infer the corresponding internal force and tangent stiffness matrix. The standard numerical routine is the finite difference method, which evaluates the target function multiple times under a small real-valued perturbation. Unfortunately, the subtractive cancellation prevents us from setting this perturbation sufficiently small, and the regular finite difference is doomed for computing problems requiring a high-accuracy derivative evaluation. In this paper, we graft a new finite difference scheme, namely the complex-step finite difference (CSFD), with physics-based animation. CSFD is based on the complex Taylor series expansion, which avoids subtractions in first-order derivative approximation. As a result, one can use a very small perturbation to calculate the numerical derivative that is as accurate as its analytic counterpart. We accelerate the original CSFD method so that it is also as efficient as the analytic derivative. This is achieved by discarding high-order error terms, decoupling real and imaginary calculations, replacing costly functions based on the theory of equivalent infinitesimal, and isolating the propagation of the perturbation in composite/nesting functions. CSFD can be further augmented with multicomplex Taylor expansion and Cauchy-Riemann formula to handle higher-order derivatives and tensor-valued functions. We demonstrate the accuracy, convenience, and efficiency of this new numerical routine in the context of deformable simulation - one can easily deploy a robust simulator for general hyperelastic materials, including user-crafted ones to cater specific needs in different applications. Higher-order derivatives of the energy can be readily computed to construct modal derivative bases for reduced real-time simulation. Inverse simulation problems can also be conveniently solved using gradient/Hessian-based optimization procedures.

Hyperelastic Material Modeling Using First- and Second-Order Differential Operations with Complex Step Derivative Approximation

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

Accelerated complex-step finite difference for expedient deformable simulation

Robust and efficient implementation of finite strain generalized continuum models for material failure: Analytical, numerical, and automatic differentiation with hyper-dual numbers

An additive Mori-Tanaka scheme for elastic-viscoplastic composites based on a modified tangent linearization

A simple interpolation-based approach towards the development of an accurate phenomenological constitutive relation for isotropic hyperelastic materials

Accelerated Nonlinear Finite Element Method for Analysis of Isotropic Hyperelastic Materials Nonlinear Deformations

A mixed FFT-Galerkin approach for incompressible or slightly compressible hyperelastic solids under finite deformation

High-performance unsymmetric 8-node hexahedral element in modeling nearly-incompressible soft tissues

An efficient and accurate method for modeling nonlinear fractional viscoelastic biomaterials

Plane stress finite element modelling of arbitrary compressible hyperelastic materials

A Review on Material Models for Isotropic Hyperelasticity

Model-free data-driven inelasticity in Haigh-Westergaard space -- a study how to obtain data points from measurements

Extracting non-equilibrium material properties in visco-hyperelastic polymers based on analytical stress solution

Improved Carroll's hyperelastic model considering compressibility and its finite element implementation

A generic formulation of anisotropic thermo-elastoviscoplasticity at finite deformations for finite element codes

A Method for Determining Elastic Constants and Boundary Conditions of Three-Dimensional Hyperelastic Materials

An Inverse Model and Mathematical Solution for Inferring Viscoelastic Properties and Dynamic Deformations of Heterogeneous Structures

An Exact Consistent Tangent Stiffness Matrix for a Second Gradient Model for Porous Plastic Solids: Derivation and Assessment

STATE OF THE ART OF CONSTITUTIVE RELATIONS OF HYPERELASTIC MATERIALS

Simplifying FFT-based methods for solid mechanics with automatic differentiation