Energy-Controllable Time Integration for Elastodynamic Contact

TL;DR

This paper introduces A-search, a novel energy-controllable integrator for elastodynamic simulations that maintains stability and energy fidelity, outperforming traditional methods like BDF2 in complex collision and deformation scenarios.

Contribution

The paper presents a new integrator, A-search, which allows user-defined energy control while ensuring stability and physical accuracy in elastodynamic simulations.

Findings

A-search effectively controls energy to match user targets.

A-search outperforms BDF2 in stability and energy preservation.

The method handles large deformations and complex collisions reliably.

Abstract

Dynamic simulation of elastic bodies is a longstanding task in engineering and computer graphics. In graphics, numerical integrators like implicit Euler and BDF2 are preferred due to their stability at large time steps, but they tend to dissipate energy uncontrollably. In contrast, symplectic methods like implicit midpoint can conserve energy but are not unconditionally stable and fail on moderately stiff problems. To address these limitations, we propose a general class of numerical integrators for Hamiltonian problems which are symplectic on linear problems, yet have superior stability on nonlinear problems. With this, we derive a novel energy-controllable time integrator, A-search, a simple modification of implicit Euler that can follow user-specified energy targets, enabling flexible control over energy dissipation or conservation while maintaining stability and physical fidelity. Our method integrates seamlessly with barrier-type energies and allows for inversion-free and penetration-free guarantees, making it well-suited for handling large deformations and complex collisions. Extensive evaluations over a wide range of material parameters and scenes demonstrate that A-search has biases to keep energy in low frequency motion rather than dissipation, and A-search outperforms traditional methods such as BDF2 at similar total running times by maintaining energy and leading to more visually desirable simulations.