现代游戏物理引擎刚体动力学与反弹几何解算：从连续碰撞检测(CCD)到四元数角动量衰减实战

在现代竞技射击与高拟真动作游戏中，投掷物（高爆手雷、战术烟雾弹、震爆弹）及高速刚体载具的碰撞解算，构成了游戏物理交互的核心护城河。无数初中级开发者在实现投掷物物理模拟时，常常陷入两大典型困境：其一，由于离散时间步长积分引发的“量子穿墙”隧道效应（Tunneling Effect）；其二，由于忽略反弹切向冲量与自旋角动量耦合，导致刚体反弹运动失真、生硬，缺乏真实动能传递感。

本文将从物理数学建模底层出发，深入剖析现代游戏引擎（如 PhysX、Havok 及自研确定性物理系统）中的刚体动力学、连续碰撞检测（Continuous Collision Detection, CCD）与动量衰减算法实现。

一、弹道轨迹数值积分：从显式欧拉到辛积分（Verlet / RK4）

在游戏主循环固定时间步长（Fixed Timestep $\Delta t$）下，离散求解刚体运动状态方程是物理管线的第一步。

1. 显式欧拉法（Explicit Euler）的能量发散缺陷

最直观的状态推进公式为：

$$x_{t+\Delta t} = x_t + v_t \cdot \Delta t$$

$$v_{t+\Delta t} = v_t + a_t \cdot \Delta t$$

在具有重力与阻力场的作用下，显式欧拉法属于一阶非辛积分器。随着模拟时间的推移，离散截断误差会导致系统能量无界单调递增，表现为弹跳高度越弹越高甚至飞向天际。

2. 辛欧拉法与半隐式积分（Semi-Implicit Euler）

在工业级物理引擎中，计算开销最小且能保证相空间辛几何守恒的方案是半隐式欧拉（速度优先更新）：

$$v_{t+\Delta t} = v_t + a_t \cdot \Delta t$$

$$x_{t+\Delta t} = x_t + v_{t+\Delta t} \cdot \Delta t$$

仅通过交换更新顺序，半隐式欧拉即具备良好的能量守恒特性，长期弹跳振荡衰减平稳，是绝大多数中低速投掷物模拟的首选基石。

二、连续碰撞检测（CCD）与扫掠体积（Swept Volume）

当高速刚体（速度达 $30 \sim 80 \, \text{m/s}$）在固定 Tickrate（如 64Hz，$\Delta t \approx 15.6\text{ms}$）下运动时，单帧位移量可达 $0.5 \sim 1.2$ 米，极易彻底跨越厚度较薄的静态墙体或铁丝网碰撞体，造成灾难性的穿透穿模。

1. 扫掠体积构造（Swept Sphere / Capsule）

为消除隧道效应，必须采用连续碰撞检测（Continuous Collision Detection, CCD）：

• 将刚体在 $[t, t+\Delta t]$ 时间段内的位移矢量 $\vec{D} = \vec{v} \cdot \Delta t$ 构造为时空几何连续体（扫掠胶囊体 Swept Capsule）；
• 通过射线检测或凸包扫掠（Convex Sweep），求解时空相交方程的首个碰撞时间因子（Time of Impact, $TOI \in [0, 1]$）。

2. GJK 与 EPA 算法实时相交解算

对于任意凸多面体刚体碰撞对，物理引擎通过 GJK（Gilbert-Johnson-Keerthi）算法判断 Minkowski 差集是否包含原点。一旦处于碰撞浸润状态，立即调用 EPA（Expanding Polytope Algorithm）求解最小穿透深度向量 $\vec{d}$ 与接触法线 $\vec{n}$，保证刚体在产生物理冲量前被精确弹出表面。

三、碰撞冲量解算与切向摩擦锥模型

刚体接触碰撞瞬间，遵循经典牛顿冲量恢复定律（Restitution Law）与库仑摩擦模型（Coulomb Friction）。

1. 法向冲量解算（Normal Impulse Resolution）

设碰撞前两物体接触点相对速度为 $\vec{v}_{\text{rel}} = \vec{v}_1 - \vec{v}_2$，接触面法线为 $\vec{n}$（指向刚体 1）。
碰撞恢复系数 $e \in [0, 1]$（0 为完全非弹性碰撞，1 为完全弹性碰撞）。
沿法线方向的分离冲量标量 $J_n$ 解算公式为：

$$J_n = -\frac{(1 + e) \, (\vec{v}_{\text{rel}} \cdot \vec{n})}{\frac{1}{m_1} + \frac{1}{m_2} + \left( I_1^{-1}(\vec{r}_1 \times \vec{n}) \times \vec{r}_1 + I_2^{-1}(\vec{r}_2 \times \vec{n}) \times \vec{r}_2 \right) \cdot \vec{n}}$$

其中：

• $m_1, m_2$ 为物体质量；
• $I_1, I_2$ 为惯性张量（Inertia Tensor）；
• $\vec{r}_1, \vec{r}_2$ 为质心到碰撞接触点的矢量。

2. 切向摩擦冲量（Tangent Friction Impulse）

沿切向方向，接触面相对滑动速度为：

$$\vec{v}_t = \vec{v}_{\text{rel}} - (\vec{v}_{\text{rel}} \cdot \vec{n})\vec{n}$$

归一化切向矢量 $\vec{t} = \frac{\vec{v}_t}{\|\vec{v}_t\|}$。

根据库仑摩擦定律，摩擦冲量标量 $J_t$ 受到摩擦锥（Friction Cone）的严格钳制：

$$|J_t| \le \mu \cdot J_n$$

其中 $\mu$ 为动摩擦因数。若切向滑动冲量超出上限，则进入动摩擦打滑状态；反之则进入静摩擦锁定状态。

四、四元数姿态演进与自旋角动量耦合

在战术手雷落地反弹时，仅有线速度的改变是不够的。切向摩擦冲量会产生相对于刚体质心的力矩 $\vec{\tau} = \vec{r} \times \vec{J}$，从而瞬时赋予手雷强烈的自旋角速度 $\vec{\omega}$。

1. 角动量更新

$$\vec{L}_{t+\Delta t} = \vec{L}_t + (\vec{r} \times \vec{J})$$

$$\vec{\omega}_{t+\Delta t} = I^{-1} \vec{L}_{t+\Delta t}$$

2. 基于单位四元数（Quaternion）的姿态积分

为杜绝欧拉角万向节死锁（Gimbal Lock），物理引擎采用单位四元数 $q = [s, \vec{v}]$ 记录空间朝向。
姿态一阶导数方程为：

$$\dot{q} = \frac{1}{2} \boldsymbol{\omega}_q \otimes q$$

其中 $\boldsymbol{\omega}_q = [0, \vec{\omega}]$。通过规范化投影保证四元数模长恒等于 1：

$$q_{t+\Delta t} = \text{Normalize}\left(q_t + \frac{\Delta t}{2} \boldsymbol{\omega}_q \otimes q_t\right)$$

五、工业级确定性物理核心实现代码（C++ / GLM）

以下是抽离自工业级射击游戏物理模块的核心反弹几何与冲量解算函数：

#include <glm/glm.hpp>
#include <glm/gtc/quaternion.hpp>
#include <algorithm>

struct RigidBodyState {
    glm::vec3 position;
    glm::vec3 linearVelocity;
    glm::quat orientation;
    glm::vec3 angularVelocity;
    float mass;
    float restitution; // 恢复系数 e
    float friction;    // 摩擦因数 mu
    glm::mat3 inertiaTensorInverse;
};

struct ContactManifold {
    glm::vec3 contactPoint;
    glm::vec3 normal;
    float penetrationDepth;
};

void ResolveRigidBodyCollision(RigidBodyState& rb, const ContactManifold& contact, float deltaTime) {
    // 1. 穿透位置校正 (Baumgarte Stabilization 投影法)
    const float slop = 0.005f; // 允许穿透误差容限 (5mm)
    const float percent = 0.4f; // 穿透校正系数
    float correctionMagnitude = std::max(contact.penetrationDepth - slop, 0.0f) * percent;
    rb.position += contact.normal * correctionMagnitude;

    // 2. 接触点相对速度矢量解算
    glm::vec3 r = contact.contactPoint - rb.position;
    glm::vec3 pointVelocity = rb.linearVelocity + glm::cross(rb.angularVelocity, r);
    float normalVelocity = glm::dot(pointVelocity, contact.normal);

    // 物体正在分离，跳过冲量解算
    if (normalVelocity > 0.0f) return;

    // 3. 计算旋转转动惯量标量项
    glm::vec3 rCrossN = glm::cross(r, contact.normal);
    glm::vec3 angularFactor = glm::cross(rb.inertiaTensorInverse * rCrossN, r);
    float effectiveMassNormal = (1.0f / rb.mass) + glm::dot(angularFactor, contact.normal);

    // 4. 法向冲量标量 Jn
    float targetNormalVel = -normalVelocity * rb.restitution;
    float deltaNormalVel = targetNormalVel - normalVelocity;
    float jn = deltaNormalVel / effectiveMassNormal;
    glm::vec3 normalImpulse = contact.normal * jn;

    // 5. 切向摩擦冲量 Jt
    glm::vec3 tangentVelocity = pointVelocity - contact.normal * normalVelocity;
    float tangentSpeed = glm::length(tangentVelocity);
    glm::vec3 frictionImpulse(0.0f);

    if (tangentSpeed > 1e-4f) {
        glm::vec3 tangentDir = tangentVelocity / tangentSpeed;
        glm::vec3 rCrossT = glm::cross(r, tangentDir);
        glm::vec3 angularFactorT = glm::cross(rb.inertiaTensorInverse * rCrossT, r);
        float effectiveMassTangent = (1.0f / rb.mass) + glm::dot(angularFactorT, tangentDir);

        float jt = -tangentSpeed / effectiveMassTangent;
        // 库仑摩擦锥约束
        float maxFriction = jn * rb.friction;
        jt = std::clamp(jt, -maxFriction, maxFriction);
        frictionImpulse = tangentDir * jt;
    }

    // 6. 应用合冲量到刚体线速度与角速度
    glm::vec3 totalImpulse = normalImpulse + frictionImpulse;
    rb.linearVelocity += totalImpulse / rb.mass;
    rb.angularVelocity += rb.inertiaTensorInverse * glm::cross(r, totalImpulse);
}

六、多端网络同步与确定性回放工程考量

在高精度多人物理同步中，若直接依赖浮点数物理引擎，跨平台（x86 与 ARM / NEON）浮点数舍入误差会在数十秒后导致两端弹道轨迹严重分歧。

1. 定点数数学库（Fixed-Point Math）：对弹道抛物线、碰撞点及四元数姿态统一采用 Q16.16 或 Q32.32 定点数解算，杜绝编译器优化指令集（FMA / AVX）造成的单比特差异；
2. 状态快照与回滚预测（Client-Side Prediction & Reconciliation）：本地客户端先行模拟投掷轨迹，服务器按固定 Tick 广播刚体物理快照（位置、速度、四元数）。一旦检测到偏差超出阈值，采用四阶 Hermite 样条插值平滑插值修正视觉网格，兼顾物理精准度与前台丝滑手感。

结语

高品质的游戏物理体验，绝不是简单调用一个 AddForce 或修改物理材质的 Bounciness 属性所能企及的。只有深入理解离散时间积分、CCD 时空扫掠、接触面摩擦锥约束与四元数姿态演进的数学物理本质，才能打造出经得起硬核竞技考验的工业级物理系统。