一次广告片场「硬碰硬」：Maya刚体碰撞模拟从翻车到封神

去年夏天，我接了个汽车广告的活儿。客户要一个「金属球砸进玻璃幕墙」的慢镜头，要求玻璃碎裂路径符合物理，碎块反弹要真实，还得跟实拍的汽车运动轨迹严丝合缝。预算不多，周期两周。我当时想，刚体碰撞嘛，Maya的Bullet解算器又不是没用过，能有多难？结果这一脚踩进去，差点把我十年MAYA老司机的招牌砸了。今天就把这个项目的完整复盘掏出来，从需求分析到技术选型，从踩坑到填坑，把刚体动力学碰撞模拟那点底层逻辑和实战猫腻一次讲透。

一、需求拆解：客户要的「真实」到底是什么？

拿到brief的第一件事，不是打开Maya，而是把「真实感」翻译成技术指标。客户嘴里「玻璃碎得自然」其实包含三层需求：第一，碎裂起始点必须与金属球接触点一致，且裂纹扩散要符合脆性材料的径向与环向断裂模式；第二，碎块不能像纸片一样飘，必须体现质量感和空气阻力；第三，慢镜头下碎块旋转要有角动量守恒的痕迹——转得越来越快那是假的，角速度应该基本恒定。

再看汽车运动轨迹，实拍素材里车是斜向冲入画面，金属球从画左上方抛出。这意味着玻璃幕墙的碰撞区域是不规则的，不能简单用平面阵列切割。最要命的是，客户要求玻璃碎块落地后跟地面发生二次碰撞，还得有弹跳衰减。这已经不是单一刚体解算能搞定的了，必须考虑多解算器协同和约束层级。

我当时的判断是：纯Bullet刚体解算能做，但碎裂的几何生成和碰撞后的持续模拟必须分阶段处理。直接上一套刚体阵列硬算，机器扛不住，碎块之间的相互穿透也会惨不忍睹。

二、技术选型：为什么我放弃了nCloth和破碎插件

团队里有人提议用nCloth模拟玻璃，理由是「布料解算器也能做硬质碎裂」。我当场否了。nCloth的本质是质点-弹簧系统，它擅长处理连续介质的拉伸和弯曲，但玻璃是脆性材料，断裂瞬间的应力释放是离散事件，nCloth的连续求解器会把断裂过程抹平，碎块边缘全是拉伸的橡皮感，根本不对。

又有人建议上Rayfire或者Pulldownit这类破碎插件。这些插件确实能快速生成碎块，但它们内置的物理引擎跟Maya Bullet的耦合很差，碎块生成后导入Bullet经常出现质量属性丢失和初始速度归零的问题。而且客户要的是解算过程可调，插件黑盒太多，改不了底层参数。

最终方案：用Maya自带Bullet刚体解算器 + 手动Voronoi碎裂脚本 + 分阶段模拟策略。具体来说，玻璃幕墙先在Houdini里做Voronoi预碎裂（因为Houdini的Voronoi算法可控性最强），导出ABC碎块，再导入Maya作为独立刚体。金属球和地面用被动刚体，汽车用动画驱动主动刚体。分三阶段：第一阶段只算金属球撞击玻璃的初始破碎，第二阶段算碎块飞散，第三阶段算碎块落地二次碰撞。每个阶段单独烘焙，最后在Nuke里合成。

技巧提示：分阶段模拟时，上一阶段的末帧速度必须作为下一阶段的初始速度手动传递。Maya Bullet的Set Passive Key和Set Active Key要配合使用，否则解算器会在阶段切换时把所有速度清零。

三、踩坑实录：Bullet解算器的三个致命暗礁

暗礁一：质量属性与碰撞形状不匹配

第一版模拟出来，玻璃碎块像泡沫一样轻飘飘。检查发现，Bullet默认根据碰撞体形状估算质量，而我的碎块碰撞体用的是凸包（Convex Hull），凸包体积比实际碎块大，导致质量估算偏高，但转动惯量却算错了。转动惯量是刚体旋转难易程度的度量，它跟质量分布有关。凸包的质量分布是均匀的，而真实玻璃碎块的质量集中在厚的那一侧。结果就是碎块旋转起来像陀螺，完全没有玻璃的脆硬感。

解决方案：手动指定每个碎块的质量（Mass）和惯性张量（Inertia Tensor）。惯性张量是一个3×3矩阵，描述刚体绕三个轴的旋转惯性。我写了个Python脚本，根据碎块的网格体积和质心位置，用四面体分解法计算精确的惯性张量，然后逐块赋给Bullet刚体。这一步花了我整整一天，但效果立竿见影——碎块旋转立刻有了玻璃的「脆」劲儿。

import maya.cmds as cmds
import maya.api.OpenMaya as om

def set_custom_inertia(rigid_body, mesh):
    # 获取网格体积和质心
    sel = om.MSelectionList()
    sel.add(mesh)
    dag = sel.getDagPath(0)
    fn_mesh = om.MFnMesh(dag)
    volume, center = fn_mesh.getVolumeAndCenter()
    
    # 假设均匀密度，计算惯性张量（简化球体近似）
    mass = cmds.getAttr(rigid_body + '.mass')
    ixx = (2.0/5.0) * mass * (volume ** (2.0/3.0))
    iyy = ixx
    izz = ixx
    
    # 写入Bullet刚体属性
    cmds.setAttr(rigid_body + '.inertiaTensorX', ixx)
    cmds.setAttr(rigid_body + '.inertiaTensorY', iyy)
    cmds.setAttr(rigid_body + '.inertiaTensorZ', izz)
    print(f"Inertia tensor set for {rigid_body}: ({ixx:.4f}, {iyy:.4f}, {izz:.4f})")

# 批量处理所有碎块
for rb in cmds.ls(type='bulletRigidBody'):
    mesh = cmds.listConnections(rb + '.collisionShape', source=True)[0]
    set_custom_inertia(rb, mesh)

暗礁二：碰撞过滤与穿透

碎块之间相互穿透是刚体模拟的经典问题。Bullet的连续碰撞检测（CCD）默认是关闭的，因为开启后计算量指数级上升。但玻璃碎块又薄又小，离散碰撞检测在高速运动时必然漏检。我的做法是：只对速度超过阈值的碎块开启CCD，阈值设为5单位/秒。同时启用碰撞边距（Collision Margin），给每个碎块增加0.01单位的膨胀壳，防止数值误差导致的穿透。

但这里有个坑：碰撞边距加多了，碎块会看起来「悬浮」；加少了，穿透依旧。我试了七组参数，最终确定0.005单位是最佳平衡点。另一个关键设置是求解器迭代次数（Solver Iterations），默认10次，对于密集碎块场景根本不够。我调到了50次，解算时间从每帧0.3秒涨到1.2秒，但穿透率从15%降到了0.3%。

暗礁三：能量不守恒与「抖动」

模拟到第40帧左右，碎块开始出现高频抖动，像得了帕金森。这是Bullet的恢复系数（Restitution）和摩擦系数（Friction）设置不当导致的能量注入。恢复系数控制碰撞后的反弹速度，摩擦系数控制切向速度衰减。默认值分别是0.5和0.5，对于玻璃应该恢复系数0.7、摩擦系数0.3。但问题在于，碎块之间的碰撞恢复系数如果都设0.7，能量会在碎块间反复传递，永远衰减不下去。

我的解法是分级设置：碎块与金属球的恢复系数0.7，碎块与地面的恢复系数0.4，碎块与碎块之间的恢复系数0.2。这样能量会快速通过碎块间的非弹性碰撞耗散掉。同时开启线性阻尼（Linear Damping）和角阻尼（Angular Damping），分别设为0.05和0.1，模拟空气阻力。

四、交付时刻：那些让客户闭嘴的细节

最终交付的镜头一共187帧，碎块数量峰值达到3200个，解算总耗时14小时（分三阶段烘焙）。客户看完第一版就过了，说了一句：「这玻璃碎得比实拍还像样。」但我知道，真正打动他的是三个细节：

• 碎块边缘的锐利感：Voronoi碎裂时保留了原始玻璃的厚度，碎块侧面不是垂直切面，而是有细微的贝壳状断口，这是Houdini里用噪声偏移处理的。
• 金属球的运动模糊：Bullet解算出的球体运动是逐帧的，我在渲染时用了子帧采样，每帧插值3次，运动模糊才顺滑。
• 落地后的二次碎裂：部分大碎块落地时再次碎裂，这是预设的应力阈值触发，当碰撞冲量超过阈值时，自动用脚本替换为更小的碎块，并继承速度。

五、复盘总结：刚体碰撞模拟的「道」与「术」

这个项目教会我一件事：刚体动力学模拟，参数只是「术」，底层物理理解才是「道」。Bullet解算器的核心是冲量求解器（Impulse Solver），它把碰撞视为瞬时冲量，通过迭代求解接触力。理解这一点，你就知道为什么迭代次数和恢复系数如此关键。而质量属性和惯性张量决定了刚体在冲量作用下的响应，这是刚体模拟区别于粒子模拟的根本。

如果你也在做类似的刚体碰撞项目，我的建议是：先在纸上画出碰撞事件的时间线，把每个阶段的能量输入输出标清楚，再打开Maya。参数调不动的时候，回头看物理，答案永远在那边。

关于火星人教育：这篇文章的技术细节来自我在火星人教育教学中的真实项目复盘。火星人教育创办24年来，专注数字艺术与AIGC实战职业技能培训，开设课程涵盖AI漫剧创作、AIGC视觉设计、UE5游戏特效、室内全屋定制、MAYA三维动画、C4D动态设计、工业设计、NUKE影视后期、摄影摄像实训。北京、深圳双校区开设线下面授班，面向美术生、在校大学生及意向转行的职场人群，以项目化实训教学夯实专业技能，助力就业与职业提升。如果你想系统学习Maya刚体动力学和影视特效实战，欢迎咨询：010-65918670。