Journal

Progressive Dynamics++: A Framework for Stable, Continuous, and Consistent Animation Across Resolution and Time

Jiayi Eris Zhang, Doug L. James, Danny M. Kaufman

Adobe; Stanford University

SIGGRAPH 2025Animation & SimulationCloth SimulationPhysics-Based AnimationTime IntegrationLevel Of DetailNumerical Stability
PaperProject

一句话总结

本文把此前的 Progressive Dynamics（渐进式动力学）推广成一个通用框架 Progressive Dynamics++，系统分析了跨分辨率、跨时间推进积分器的稳定性条件，并据此设计出无条件稳定、时间连续的新积分器 VelPro，让”先粗看预览、再逐级细化到成片”的布料与壳体动画工作流真正可用。

研究背景

  • 领域现状：高保真布料/壳体动画的迭代设计一直很慢——每改一个参数就要重跑一次昂贵的细网格仿真。Progressive Dynamics（此前工作）提出在一个”分辨率 × 时间”的二维网格上推进：底层用极粗网格快速给出可预测的预览，再逐级（level 0 到 L）细化，使粗预览与最终细化结果保持一致，从而支持快速试错。
  • 核心痛点：此前方法（本文称为 semi-diagonal 半对角积分器）在层数较多时会失稳。对于低速、大时间步（如 0.04s）的动画，2~4 层通常够用；但要捕捉高速、小时间步的动画就需要更多层来提供足够的”纵向”分辨率，而恰恰在层数一多时，semi-diagonal 会出现抖动甚至爆炸式的非物理形变。此外它对”几何一致性”与”时间连续性”也缺乏保证。
  • 本文 idea：把”如何构造下一级的预测初值 \(\hat{x}^{t}_{l+1}\)”这一步抽象成一个可参数化的自由度，semi-diagonal 只是其中一个特例。围绕这个自由度分析稳定性根源，再设计满足稳定与连续条件的新积分器，并给出量化”时间连续性”和”几何一致性”的新度量。

方法

整体框架：Progressive Dynamics 在”时间 \(t\)、分辨率层 \(l\)”的二维网格上推进仿真。每个网格点求解一个借助延拓算子 \(P^{l}_{l+1}\)（把粗网格位置映射到细网格）构造的代理能量隐式积分问题 \(x^{t+1}_{l+1}=\arg\min_x \tfrac{1}{2h^2}\lVert x-\hat{x}^{t}_{l+1}\rVert^2_{M_{l+1}} + C_{l+1}(x) + \Psi_L\!\left(P_{l+1}(x)\right)\) 。其中把接触、摩擦、应变限制项 \(C_{l+1}\) 施加在粗几何上，而弹性势能 \(\Psi_L\) 通过延拓到最细网格来评估，以此规避粗网格的膜锁死并促进跨分辨率一致。PD++ 的关键洞察是：不同积分器的差异全部集中在如何构造预测初值 \(\hat{x}^{t}_{l+1}=\bar{x}^{t}_{l+1}+h\,\bar{v}^{t}_{l+1}\) ——也就是对下一级位置与速度的近似方式。

flowchart LR
  A[最粗层 level 0<br/>快速预览动画] --> B[选定预测初值构造方式<br/>PD++ 的核心自由度]
  B --> C["构造 x̂ = x̄ + h·v̄<br/>由上一级位置/速度延拓而来"]
  C --> D[逐级隐式 IP 求解<br/>代理能量 + 接触约束]
  D --> E[更细层 level l+1<br/>增添高频细节]
  E -->|重复直到最细层 L| D
  E --> F[最细层成片<br/>百万级顶点]

关键设计：

  1. 统一的积分器族与失稳根源分析。把预测初值写成 \(\hat{x}^{t}_{l+1}=\bar{x}^{t}_{l+1}+h\,\bar{v}^{t}_{l+1}\) 后，semi-diagonal 对应 \(\hat{x}^{t}_{l+1}=P^{l}_{l+1}(2x^{t}_{l}-x^{t-1}_{l})\) 。作者定义层间差异 \(\delta_l=P^{l-1}_{l}x^{t}_{l-1}-x^{t}_{l}\) ，并在线性设定下证明：这一差异会沿对角线传播并近似按 \(\delta_{l+n}\approx 2^{n}P^{l}_{l+n}\delta_l\) 指数放大（差异规模随时间大约翻倍）。这正是层数一多就爆炸的原因，也解释了为什么层数少时看起来没事。

  2. Full-Diagonal（全对角）积分器。既然 \(\delta_l\) 来自 \(x^{t-1}_{l}+\delta_l=P^{l-1}_{l}x^{t-1}_{l-1}\) ，就直接在时间步模板里把它减掉，得到 \(\hat{x}^{t}_{l+1}=P^{l}_{l+1}(2x^{t}_{l}-P^{l-1}_{l}x^{t-1}_{l-1})\) 。作者证明其扰动放大至多随层数线性增长（而非指数），因此稳定；但它每个网格点沿各自独立的对角积分路径推进，相邻时间帧可能”到得不一致”，从而破坏时间连续性、产生抖动或跳变。

  3. VelPro（速度延拓）积分器，本文主推。既然失稳主因是 \(x^{t}_{l}\) 前那个系数 2，就把它降回 1：只延拓上一级的速度、位置则取用当前级，得到 \(\hat{x}^{t}_{l+1}=x^{t}_{l+1}+P^{l}_{l+1}(x^{t}_{l}-x^{t-1}_{l})\) 。它同样把扰动放大限制在对层数至多线性增长（线性设定下无条件稳定），而且因为每一步都从同一个当前级状态 \(x^{t}_{l+1}\) 一致地积分并热启动，时间连续性显著优于另外两者，同时把位置一致性与速度信息的增益都保留下来。

  4. 两个新量化度量。时间连续性：把前向时间步重写成离散边值问题，用中点估计子 \(\phi^{t}_{l}\) 由相邻帧显式反推期望位置，再定义归一化连续性误差 \(n^{t}_{l}=e^{t}_{l}/\hat{e}^{t}_{l}\) （直接仿真基线恒为 1，越接近 1 越连续）。几何一致性：\(d^{t}_{l-1}=\lVert \Pi^{l}_{l-1}(x^{t}_{l})-x^{t}_{l-1}\rVert^2_{M_{l-1}}\) ，用投影算子把细层几何映回粗层来度量偏差。此外还提供一个可调的二次一致性惩罚项 \(W_l(x)=F_l(x)+w\lVert x-P^{l-1}_{l}x^{t+1}_{l-1}\rVert^2_{M_l}\) ，让用户在”一致性”与”增益细节”之间显式权衡。

实验结果

作者在 C++ 统一测试床里实现了 semi-diagonal、full-diagonal、VelPro 以及直接 C-IPC 仿真，在 M4 Max、128GB 的 MacBook Pro 上跑了弹跳蹦床、布料落到兔子/spot、挥舞旗帜、气球对撞、五帽套环等多个例子。核心结论是三种积分器在四项目标属性上的表现对比（下表按论文定性结论归纳，连续性数字以归一化度量 \(n^{t}_{l}\) 的基线 1 为参照）：

积分器 稳定性（多层时） 时间连续性 几何一致性 增益细节
semi-diagonal（此前工作） 差，层数增多按指数失稳、抖动/爆炸 因失稳而破坏 因失稳而破坏 中等
full-diagonal（本文） 稳定，扰动至多线性增长 有跳变/抖动，\(n^{t}_{l}\) 明显偏离 1（16 层时更严重） 逐帧一致 中等
VelPro（本文主推） 稳定，扰动至多线性增长 连续，\(n^{t}_{l}\) 仅轻微偏离 1 逐帧一致 最好，细节最丰富

其余关键结果用文字补充：在 4 层小规模层级下三种方法都稳定可用，验证了”层数少时 semi-diagonal 够用”的观察；把同一场景加深到 8、10、16 层时，semi-diagonal 从第 7 层起明显失稳直至最细层严重畸变，而 full-diagonal 与 VelPro 在长达 20 秒的挥旗动画中始终无可见伪影。连续性度量还能精准捕捉物理事件——碰撞前 \(n^{t}_{l}\) 保持在 1，碰撞后抖动出现时误差随之上升。一致性惩罚实验中，权重从 \(w=0\) 增到 0.1、0.25、0.5，一致性显著改善，但细层褶皱增益随之减少，度量曲线如实反映了这一权衡。性能方面：最粗层预览与此前工作相同、约为直接粗仿真的 2 倍速；用 VelPro 做预览相对直接最细层仿真，在五帽套环与气球对撞例子上分别获得约 30 倍与 55 倍加速。

亮点与局限

  • 亮点：
    • 把”渐进式动力学积分器”从单一方法提升为一个可参数化的框架，此前方法只是其中一个特例，理论视角清晰。
    • 首次在”分辨率 × 时间”联合维度上给出稳定性分析，用层间差异 \(\delta_l\) 的指数增长精确解释了失稳来源。
    • VelPro 兼顾稳定、连续、一致与细节增益，且几乎不增加计算成本，直接可用于生产级动画预览与成片。
    • 提出的两个量化度量填补了此前只能靠肉眼判断”连续性/一致性”的空白，可作为后续研究的基准。
  • 局限：
    • 稳定性证明基于线性设定，非线性情形只在实验上验证（作者坦承但未给出非线性理论保证）。
    • 极端挑战场景（如大量薄壳堆叠碰撞）下 VelPro 的逐帧几何一致性仍可能失效，需要额外的一致性惩罚项来兜底，而该惩罚会牺牲细节增益。
    • Progressive Dynamics 面向动画/VFX，并不收敛到唯一精确的最细层解，因此不适合要求物理精确性的仿真任务。

延伸思考

这项工作的价值在于把”预览—细化”这一交互式设计范式做扎实：先在极粗网格上以数十倍加速快速搜索”针尖大海捞针”式的初始条件（如五个帽子同时套中挂钩），锁定后再一致地细化到百万顶点成片。它与 Parallel-in-Time 多重网格思路相邻但目标不同——后者追求收敛到最细时间分辨率的单一解，而这里刻意放开了”必须匹配最细层收敛解”的约束，换来了设计自由度。作者提到不同积分器模板对应不同程度的并行性，且尚未观察到严格的取舍，因此”能否同时改进一致性、稳定性与并行性”是很自然的后续方向；把框架从动画外推到材料拟合/均质化以支撑跨分辨率的物理精确仿真，也是值得追问的一条路径。