Subspace-Preconditioned GPU Projective Dynamics with Contact for Cloth Simulation

一句话总结

用”子空间求解低频运动 + 全空间 Jacobi 松弛高频细节”的预条件框架加速投影动力学（PD），并以时间分裂方式把 IPC 无穿透接触整合进来，实现高分辨率布料在 GPU 上的高效、无穿透仿真。

研究背景

• 领域现状：布料仿真为了呈现细腻的褶皱、折痕，需要高分辨率网格离散；投影动力学（PD）通过”局部投影 + 全局二次求解”的交替迭代成为常用的快速求解框架，但高分辨率带来的自由度爆炸让全阶预分解和求解变得昂贵。
• 核心痛点：一是非线性布料动力学在高分辨率下计算量巨大；二是碰撞与自碰撞无处不在，而当前最强的接触建模——增量势接触（IPC）——会显著增大非线性系统的刚度，使求解更难收敛。PD 的全局步依赖固定系统矩阵的预分解，一旦引入随步变化的接触项就难以复用预分解。
• 本文 idea：观察到布料整体运动由低频模态主导，可以把大部分运动放进一个小规模子空间里高效求解（可预分解），再用全阶并行 Jacobi 松弛补上子空间无法表达的高频细节；接触则用时间分裂拆成”弹性步 + 穿透校正步”，并用子空间 BFGS 把接触二次代理逐步并入子空间矩阵，避免反复重分解。

方法

整体框架：在优化式后向欧拉时间积分下用 PD 求解。每个全局步先做一次子空间求解拿到低频主运动，再在全阶上跑固定次数的 Jacobi 松弛压制高频残差；接触通过时间分裂被拆成一个用 Newton-PCG 求解的穿透校正步和一个 PD 弹性步，并借 CCD 保证无穿透。

flowchart LR
  A["预测位置(惯性)"] --> B["弹性步: 子空间预条件 PD"]
  B --> C["局部步: 三角形并行投影"]
  C --> D["全局步 I: 子空间求解(低频)"]
  D --> E["全局步 II: 全阶 Jacobi 松弛(高频)"]
  E --> F{"收敛?"}
  F -->|否| C
  F -->|是| G["穿透校正步: Newton-PCG + CCD"]
  G --> H["无穿透新状态"]

关键设计：

1. 子空间预条件的 PD 全局步。全局步本是一个系统矩阵固定的二次问题。作者用一个基矩阵把位移投影到低维子空间，子空间系统规模小、可预分解，从而高效捕获低频主运动。但子空间无法表达褶皱这类高频特征，于是紧接着在全阶上做若干次 Jacobi 松弛来削减高频残差。这种”子空间求解 + 全阶松弛”的组合正是它相较纯 Chebyshev 加速的全阶 PD 能大幅提速的核心（示例场景从 26.2s/帧 降到 3.0s/帧）。

2. 面向布料的 B 样条子空间。布料由缝纫裁片组成，每片有二维 UV 平面参数化。作者只需在二维平面上构造基：用二次 B 样条构造一维基再张成三维基，天然满足单位分解（partition of unity）。单位分解对仿真质量很重要，缺失它会带来明显伪影。

3. 接触的时间分裂与子空间 BFGS。接触沿用带二次代理的分裂方案：把动量方程拆成弹性步与穿透校正步，穿透校正用 Newton-PCG 求解并配合 CCD 保证无穿透。引入接触二次代理后全局步不再是纯线性系统，为避免每步重分解，作者用子空间 BFGS 做低秩更新去近似接触带来的 Hessian 变化，并让这些低秩更新在同一时间步的所有全局步之间共享，从而把接触信息渐进并入子空间矩阵而不破坏可预分解结构。整套算法在 GPU 上高效实现。

实验结果

在多个高分辨率布料场景上与 GPU CIPC、PD-IPC 对比单帧耗时（越低越好），本文方法在保持无穿透的同时显著更快：

场景 (顶点数 / 基数)	本文↓	GPU CIPC↓	PD-IPC↓
252K V / 8427 基	27s	241s	49s
104K V / 7665 基	10s	90s	18s
285K V / 7119 基	46s	1033s	134s
232K V / 9432 基	23s	362s	253s

可以看到本文相对 GPU CIPC 普遍快约一个数量级（最大场景 46s vs 1033s，约 22 倍），相对 PD-IPC 也快约 2–10 倍，且优势在顶点数与接触复杂度更高的场景更明显。

亮点与局限

• 亮点：
  • 把子空间积分与全阶并行松弛结合，同时压制低频与高频残差，兼顾效率与细节表现。
  • 借二维 B 样条子空间贴合布料裁片结构，并满足单位分解，保证仿真质量。
  • 通过时间分裂 + 子空间 BFGS 低秩更新，将 IPC 无穿透接触无缝并入 PD，且避免昂贵的系统重分解。
  • 高分辨率场景下相较现有 GPU 求解器有约一个数量级的提速。
• 局限：
  • 高速运动时时间分裂误差会引入阻尼（damping）效应，需要自适应子步（adaptive substepping）来缓解。
  • 穿透校正步中的 Newton 法可能出现过冲（overshooting），需要为该步设计专门的求解器。

延伸思考

• 子空间 + 全阶松弛的思路与经典多重网格/预条件求解相通，B 样条子空间相当于人为构造的低频光滑算子；能否用数据驱动或几何自适应的方式学习更优基，进一步压低基数与残差是值得探索的方向。
• 时间分裂带来的阻尼与穿透校正的过冲，本质都源于把强耦合的弹性—接触问题近似解耦；如何在保持可预分解优势的同时更精确地耦合两者，是该类方法的共同难点。
• 该框架天然依赖缝纫裁片的 UV 参数化来构造子空间，对无明显 UV 结构或拓扑复杂的可变形体如何推广，是把方法从布料扩展到通用弹性体时需要回答的问题。