统一薄壳-固体-杆件接触求解器：PPF 的屏障势能方法与工程实践

在物理仿真领域，接触求解一直是计算复杂度最高的子问题之一。传统方案往往针对不同几何类型（薄壳、固体、杆件）分别设计求解器，导致代码冗余且难以处理异构物体间的耦合。ZOZO 技术团队开源的 ppf-contact-solver 提供了一种统一视角：基于增量势接触（Incremental Potential Contact, IPC）方法，在位置级别直接求解无穿透约束，同时覆盖 codimension-0（固体）、codimension-1（薄壳）、codimension-2（杆件）三类几何体。

核心设计：屏障势能而非惩罚力

PPF 求解器的关键技术决策是放弃传统的惩罚力（penalty force）或基于速度互补的 LCP 框架，转而采用 IPC 的屏障势能（barrier potential）模型。该方法将非穿透约束转化为一个光滑可微的能量项，嵌入到隐式时间积分的优化框架中。具体而言，对于任意两个接近的图元对（primitive pair），当它们之间的距离低于阈值时，屏障函数会产生趋于无穷大的排斥势能，从而在数学上保证求解结果严格无穿透。

这种位置级（position-level）约束处理相比速度级（velocity-level）方法具有显著优势：后者虽然能保证接触瞬间速度一致，但在时间推进后仍可能因 advection 的非线性而产生几何穿透。PPF 通过在每个时间步求解一个约束优化问题，直接得到满足无穿透条件的构型，从根本上避免了这一问题。

有限元基础与严格布料约束

求解器采用有限元方法（FEM）作为可变形体的物理基础。对于薄壳，使用三角形离散化并引入膜能量和基于二面角的弯曲能量；对于杆件，则采用边拉伸能量和中心线弯曲能量。这种离散方式与工业界标准的 FEM 流程兼容，便于与现有资产 pipeline 集成。

特别值得注意的是 PPF 对布料仿真的处理。传统布料模拟常面临 "橡胶化" 问题 —— 即布料在张力下产生非物理的弹性拉伸。PPF 通过引入严格的应变限制（strain-limited cloth）约束，确保布料在受力时保持近似不可拉伸的特性，这对服装仿真等应用场景至关重要。

GPU 执行与性能权衡

PPF 支持单精度 GPU 执行，能够处理大规模接触场景。求解器采用牛顿型优化器配合线性搜索（line search）和连续碰撞检测（CCD）过滤器，确保在复杂接触情况下仍能收敛。然而，开发者坦诚地指出当前实现存在性能局限：求解器相对较慢，UI 仍需大量手动配置，在极端情况下可能出现停滞。

从工程实践角度，建议将 PPF 定位为离线仿真工具而非实时引擎组件。对于需要高精度接触处理的特效镜头或服装设计验证场景，其无穿透保证和统一几何处理能力值得额外的计算开销。

Blender 集成与部署模式

ZOZO 团队已发布 Blender 5.0 插件，将求解器能力引入 DCC 工具链。插件支持三种部署模式：

1. 本地执行：Windows/Linux 环境下的单机 GPU 计算
2. 远程执行：通过局域网连接计算节点
3. 云端 GPU：支持 vast.ai、AWS 等云平台的弹性计算资源

这种灵活的部署架构使小型工作室也能获得此前仅大型 VFX 公司具备的高精度仿真能力。

可落地参数配置清单

基于项目文档和工程实践，整理以下关键参数配置建议：

局限与未来方向

PPF 当前实现存在若干限制：首先，求解器仅处理粒子 - 流体穿透，对于边 - 流体或三角形 - 流体的距离度量尚未定义，这要求固体粒子分布相对密集以避免穿透；其次，方法目前仅支持疏水接触，无法处理其他接触角场景；最后，半拉格朗日 advection 引入了数值耗散，对于需要长期能量守恒的场景可能需要额外处理。

尽管如此，PPF 在统一多几何类型接触求解方面的探索具有重要意义。它将 IPC 方法从固体 - 固体接触扩展到固体 - 流体耦合领域，为复杂物理场景的仿真提供了新的技术路径。对于需要严格无穿透保证的离线仿真任务，PPF 是值得评估的开源选项。

资料来源

• GitHub 仓库: st-tech/ppf-contact-solver — ZOZO 开源的接触求解器实现
• 80 Level: "New Open-Source Physics Engine For Blender Released" — Blender 插件发布详情与部署模式说明
• arXiv:2505.12539 "Penetration-free Solid-Fluid Interaction on Shells and Rods" — 屏障势能方法在固 - 流耦合中的理论基础

systems