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.

Inverse deformation analysis: an experimental and numerical assessment using the FEniCS Project

Fast Corotational Simulation for Example-Driven Deformation.

Geometrically nonlinear deformation reconstruction of based on Euler–Bernoulli beam theory using a nonlinear iFEM algorithm

Full-field deformation reconstruction of beams using the inverse Finite Element Method: Application to thin-walled structures

An Efficient Large Deformation Method Using Domain Decomposition

Two-Dimensional Deformation Estimation of Beam-Like Structures Using Inverse Finite-Element Method: Theoretical Study and Experimental Validation

Nonlinear deformation monitoring of elastic beams based on isogeometric iFEM approach

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

FEniCS implementation of the Virtual Fields Method (VFM) for nonhomogeneous hyperelastic identification

Geometrically Nonlinear Deformation Reconstruction Based on iQS4 Elements Using a Linearized Iterative iFEM Algorithm

Inverse finite element method with energy-based regularization for deformation reconstruction and structural health monitoring

Geometrically nonlinear shape sensing of anisotropic composite beam structure using iFEM algorithm and third-order shear deformation theory

A Variational Bayesian Inference Theory of Elasticity and Its Mixed Probabilistic Finite Element Method for Inverse Deformation Solutions in Any Dimension

An Explicit Total Lagrangian Fragile Points Method for Finite Deformation of Hyperelastic Materials

Multi-compartment poroelastic models of perfused biological soft tissues: implementation in FEniCSx

Sparse Surface Constraints for Combining Physics-based Elasticity Simulation and Correspondence-Free Object Reconstruction

Accelerated complex-step finite difference for expedient deformable simulation

Absolute Nodal Coordinate Formulation-Based Shape Sensing Approach for Large Deformation: Plane Beam

Application of Inverse Finite Element Method to Shape Sensing of Curved Beams

Comparative study of face-based smoothed point interpolation method and linear finite element method for elastoplastic and large deformation problems in geomaterials

Variable Thickness Strain Pre-Extrapolation for the Inverse Finite Element Method