Energy-Controllable Time Integration for Elastodynamic Contact
作者: Kevin You, Juntian Zheng, Minchen Li
分类: cs.GR, math.NA
发布日期: 2026-02-08
💡 一句话要点
提出能量可控时间积分方法A-search，用于弹性体动态接触仿真，提升稳定性和视觉效果。
🎯 匹配领域: 支柱四：生成式动作 (Generative Motion)
关键词: 弹性体仿真 动态仿真 时间积分 能量控制 隐式方法
📋 核心要点
- 现有隐式积分方法（如隐式欧拉和BDF2）在弹性体动态仿真中虽稳定，但能量耗散不可控，辛方法虽能保持能量，但稳定性不足。
- 论文提出A-search，一种能量可控的时间积分器，通过修改隐式欧拉，使用户能够指定能量目标，灵活控制能量耗散或守恒。
- 实验表明，A-search在保持能量方面优于传统方法BDF2，能产生更具视觉吸引力的仿真效果，且适用于大型变形和复杂碰撞。
📝 摘要（中文）
弹性体的动态仿真一直是工程和计算机图形学中的一项长期任务。在图形学中，隐式欧拉和 BDF2 等数值积分器因其在大时间步长下的稳定性而备受青睐，但它们往往会不受控制地耗散能量。相比之下，隐式中点等辛方法可以保持能量，但并非无条件稳定，并且在适度刚性问题上会失效。为了解决这些局限性，我们提出了一类用于哈密顿问题的通用数值积分器，该积分器在线性问题上是辛的，但在非线性问题上具有卓越的稳定性。在此基础上，我们推导出一种新颖的能量可控时间积分器 A-search，它是隐式欧拉的一个简单修改，可以遵循用户指定的能量目标，从而灵活控制能量耗散或守恒，同时保持稳定性和物理保真度。我们的方法与 barrier-type 能量无缝集成，并允许无反转和无穿透保证，使其非常适合处理大型变形和复杂碰撞。对各种材料参数和场景的广泛评估表明，A-search 倾向于将能量保持在低频运动中而不是耗散，并且 A-search 在相似的总运行时间下优于 BDF2 等传统方法，因为它保持了能量并产生了更具视觉吸引力的模拟。
🔬 方法详解
问题定义：弹性体动态仿真中，如何在保证数值稳定性的前提下，精确控制能量的耗散与守恒是一个关键问题。传统的隐式积分方法（如隐式欧拉和 BDF2）虽然稳定，但会引入不必要的能量耗散，导致仿真结果不真实。而辛方法虽然能保持能量，但稳定性较差，难以处理复杂场景和高刚度问题。
核心思路：论文的核心思路是设计一种既具有隐式方法的稳定性，又能精确控制能量耗散的时间积分方案。通过对隐式欧拉方法进行改进，引入能量控制机制，使得用户可以根据需要指定能量目标，从而灵活地控制能量的耗散或守恒。
技术框架：A-search 方法基于隐式欧拉方法，并在其基础上进行修改。其整体流程如下：1. 使用隐式欧拉方法进行初步的时间步进；2. 计算当前状态下的能量；3. 根据用户指定的能量目标，调整下一步的迭代方向，以实现能量的精确控制；4. 重复迭代，直到满足能量目标或达到最大迭代次数。该方法可以与 barrier-type 能量函数无缝集成，保证无反转和无穿透。
关键创新：A-search 的关键创新在于其能量控制机制。它通过在隐式欧拉的迭代过程中引入能量目标，使得积分器能够根据用户的需求，主动地调整能量的耗散或守恒。这种能量控制机制使得 A-search 既具有隐式方法的稳定性，又能够精确地控制能量，从而获得更真实、更稳定的仿真结果。与现有方法相比，A-search 能够更好地平衡稳定性和能量守恒，适用于更广泛的场景。
关键设计：A-search 的关键设计在于能量目标的设定和迭代过程中的调整策略。能量目标可以根据具体的应用场景进行设定，例如，可以设定为保持能量不变，也可以设定为按照一定的速率耗散能量。在迭代过程中，A-search 通过调整迭代方向和步长，使得能量逐渐逼近目标值。具体的调整策略需要根据具体的能量函数和仿真场景进行优化。
🖼️ 关键图片
📊 实验亮点
实验结果表明，A-search 在保持能量方面优于传统的 BDF2 方法，能够产生更具视觉吸引力的仿真效果。在相似的总运行时间下，A-search 能够更好地保持能量，避免了 BDF2 方法中出现的过度耗散现象。此外，A-search 还能够处理大型变形和复杂碰撞，适用于更广泛的仿真场景。
🎯 应用场景
该研究成果可广泛应用于计算机图形学、游戏开发、工程仿真等领域。例如，在游戏开发中，可以使用 A-search 方法来模拟角色的运动和碰撞，从而获得更真实、更稳定的游戏体验。在工程仿真中，可以使用 A-search 方法来模拟结构的动态响应，从而评估结构的安全性。
📄 摘要（原文）
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.