流体中可操控载具Three.js模拟崩溃原因及修复方案

基于《Physics for Game Developers》中的3D刚体飞机模拟器，我将其移植到Three.js的Vector3、Matrix3和Quaternion类，实现流体（空气、水等）中可操控载具（飞机、气垫船、船只）的模拟：

• 施加推力后，载具加速并在介质中运动，速度会根据舵面角度产生阻力和升力；
• 这些力偏离载具部件（主体与舵）的共同质心，产生力矩使载具产生角运动（转向）；
• 阻力和升力与速度平方成正比，阻力会限制给定推力下的最大速度；稳态时合力为零，阻力与升力的合力和推力大小相等、方向相反。

但模拟会随参数出现振荡（稳定或崩溃），振荡似乎发生在舵面法向量与阻力向量正交时——此时二者点积为零，用于计算攻角的数值趋近于零。调整参数（更小的流体密度RHO、更小的升力系数lc、更大的阻力系数dc）能延迟问题出现，但更新演示版本（将升力系数设为与攻角成正比、减小阻力系数并增大RHO）后，模拟在第227步崩溃。请问该问题的原因是什么？如何修复？

（注：相关演示代码可复现该问题，升力和阻力的计算逻辑从第215行左右开始）

原因分析 #

1. 攻角趋近零时的数值奇异点
   当舵面法向量与阻力向量正交，攻角接近零，若升力系数与攻角线性相关，此时升力会出现近乎突变的波动；同时点积为零会导致方向计算（如升力方向的推导）出现数值不稳定，甚至引发除以极小值的计算错误，直接导致力和力矩的输出异常。

2. 数值积分的稳定性不足
   若模拟使用基础欧拉积分，对于速度平方相关的非线性力系统，稳态附近的微小计算误差会被放大，引发力的方向反复切换，导致振荡不断加剧直至崩溃。这类积分方法对非线性系统的稳定性较差，无法有效约束稳态附近的扰动。

3. 参数失衡放大扰动
   增大RHO同时减小阻力系数，会让升力的相对影响大幅提升——当攻角接近零时，升力的微小波动会被放大为剧烈的力矩变化，破坏系统的动态平衡，加速模拟崩溃。

修复方案 #

1. 消除攻角的数值奇异点 #

• 在攻角计算中添加极小的偏移量（如ε = 1e-6），避免攻角精确为零，例如将攻角计算式改为angle = Math.acos(dotProduct) + ε，或在点积为零时直接赋予一个极小的默认值；
• 对舵面法向量、阻力向量先做归一化处理，确保点积计算的稳定性，避免因向量长度异常导致的点积偏差。

2. 替换为稳定的数值积分方法 #

将基础欧拉积分替换为半隐式欧拉积分或Verlet积分：

• 半隐式欧拉积分先更新速度再更新位置，对非线性力的稳定性更好；
• Verlet积分通过位置的二阶差分计算加速度，能自动抑制高频振荡，适合刚体动力学模拟。

3. 平滑升力系数的过渡逻辑 #

当攻角绝对值小于阈值（如0.01弧度）时，用平滑函数替代线性关系，避免升力突变：

// 示例：用二次函数平滑升力系数
const alphaThreshold = 0.01;
let liftCoeff;
if (Math.abs(alpha) < alphaThreshold) {
  liftCoeff = lc * (alpha * alpha / alphaThreshold);
} else {
  liftCoeff = lc * alpha;
}

也可使用Sigmoid函数或三角函数实现更自然的过渡。

4. 添加阻尼项抑制振荡 #

• 线运动阻尼：添加与速度成正比的阻尼力，公式为dampingForce = -linearDamping * velocity，降低速度的突变；
• 角运动阻尼：添加与角速度成正比的阻尼力矩，公式为dampingTorque = -angularDamping * angularVelocity，抑制转向时的振荡。

5. 平衡参数配置 #

避免过度增大RHO同时减小阻力系数：保证阻力的阻尼作用足够，让系统对微小扰动的敏感度降低；可通过测试找到RHO、阻力系数、升力系数的平衡区间，比如先固定RHO，逐步调整阻力系数至系统能稳定收敛到稳态。

内容的提问来源于Stack Exchange，提问作者匿名