C++物理引擎碰撞检测实战指南

该函数检查 x、y 轴区间是否同时重叠，仅当所有轴重叠时才判定为碰撞。

2.2 常见碰撞体类型设计与 C++ 类封装

在物理引擎实现中，碰撞体的抽象建模是核心环节。为支持多种几何形状，通常采用基类定义通用接口，派生具体类型以实现多态判据。

基础类设计

定义抽象基类 Collider，封装公共属性与虚函数，便于运行时动态调用。

class Collider {
public:
    virtual bool intersects(const Collider* other) const = 0;
    virtual ~Collider() = default;
protected:
    Vec3 position; // 碰撞体中心位置
};

该类声明了纯虚函数 intersects，强制子类实现相交检测逻辑。position 统一管理空间坐标，降低耦合。

典型子类实现

常用碰撞体包括球体与轴对齐包围盒（AABB），其 C++ 封装如下：

2.3 轴对齐包围盒（AABB）检测算法实现

基本原理

轴对齐包围盒（AABB）通过为每个物体定义一个最小和最大的坐标边界框，判断两个物体的包围盒在各轴上是否重叠，从而快速检测碰撞。

代码实现

struct AABB {
    float minX, minY, minZ;
    float maxX, maxY, maxZ;
};

bool checkCollision(const AABB& a, const AABB& b) {
    return (a.minX <= b.maxX && a.maxX >= b.minX) &&
           (a.minY <= b.maxY && a.maxY >= b.minY) &&
           (a.minZ <= b.maxZ && a.maxZ >= b.minZ);
}

该函数通过比较两个 AABB 在 X、Y、Z 三个轴上的区间是否重叠来判断碰撞。只要任一轴无重叠，则判定无碰撞，逻辑简洁高效。

性能优势

• 计算复杂度低，仅需 6 次比较
• 适用于大规模场景的初步筛选
• 易于并行化处理

2.4 分离轴定理（SAT）原理与多边形碰撞检测

分离轴定理（Separating Axis Theorem, SAT）是判断两个凸多边形是否发生碰撞的核心算法之一。其核心思想是：若存在一条轴，使得两个多边形在该轴上的投影不重叠，则这两个多边形不相交。

投影与分离轴的选取

对于两个凸多边形，需检查每条边的法线方向作为潜在分离轴。若所有轴上的投影均重叠，则判定为碰撞。

• 仅适用于凸多边形
• 需归一化法向量以保证投影准确性
• 投影计算使用点积操作

代码实现示例

std::pair<float, float> projectPolygon(const std::vector<Vec3>& vertices, const Vec3& axis) {
    float min = dot(vertices[0], axis);
    float max = min;
    for (size_t i = 1; i < vertices.length(); i++) {
        float p = dot(vertices[i], axis);
        if (p < min) min = p;
        if (p > max) max = p;
    }
    return {min, max};
}

上述函数将多边形顶点集沿指定轴投影，返回最小和最大投影值。dot 表示向量点积，用于计算顶点在单位法向量上的投影长度。后续需比较两多边形在相同轴上的投影区间是否重叠。

2.5 碰撞信息反馈：法向量、穿透深度计算

在物理引擎中，碰撞检测不仅需要判断物体是否相交，还需提供精确的碰撞信息用于响应处理。其中，法向量与穿透深度是关键参数。

法向量的作用

法向量指明了两个物体接触面的垂直方向，决定了分离力的方向。它通常由体积较小的物体指向较大的物体，确保响应的一致性。

穿透深度的计算逻辑

当两凸体相交时，可通过 分离轴定理（SAT） 找到最小穿透方向与距离。以下为简化实现示例：

// 假设已通过 SAT 确定最小穿透轴 axis 与深度 depth
Vector3 contactNormal = axis; // 法向量
float penetrationDepth = depth; // 穿透深度

// 调整方向：确保从 A 指向 B
if (dot(centerB - centerA, contactNormal) < 0) {
    contactNormal = -contactNormal;
}

上述代码中，axis 是使投影重叠最小的分离轴，depth 表示沿该轴需分离的距离。通过中心点相对位置校正法向量方向，保证反馈一致性。

• 法向量决定碰撞响应方向
• 穿透深度影响物体分离程度
• 二者共同构成解决穿透状态的基础

第三章：空间划分与性能优化策略

3.1 网格哈希与空间分块加速检测

在大规模点云或三维场景中，直接遍历所有数据进行碰撞或邻近检测效率极低。网格哈希技术通过将空间划分为固定大小的体素块，仅对相邻体素内的点进行局部检测，大幅减少计算量。

空间分块策略

采用均匀网格划分三维空间，每个网格单元存储其内部点的索引。查询时只需访问目标点所在网格及其邻接网格，避免全局搜索。

哈希映射优化存储

使用三维坐标哈希函数将网格坐标映射为一维键值，实现稀疏网格的高效存储与快速查找：

int64_t hash_key = x * 73856093 ^ y * 19349663 ^ z * 83492791;

该哈希函数选用大质数进行异或运算，有效减少碰撞概率。参数 x、y、z 为网格整数坐标，输出唯一键用于哈希表索引。

1. 将点云坐标转换为体素索引
2. 通过哈希表聚合同格内点集
3. 仅检测目标点所在及邻近体素中的候选点

3.2 动态对象管理与粗检测阶段实现

在实时渲染系统中，动态对象的高效管理是性能优化的关键。通过引入空间分区结构，如四叉树或 BVH，可快速定位移动物体并减少碰撞检测的计算量。

对象注册与更新机制

每个动态对象在初始化时注册至管理器，并分配唯一 ID。每帧根据其运动状态更新包围盒位置。

class ObjectManager {
public:
    void Update(uint32_t id, const Vector3& pos) {
        if (auto obj = objects.find(id); obj != objects.end()) {
            obj->second.Center = pos; // 更新空间位置
        }
    }
private:
    std::unordered_map<uint32_t, BoundingVolume*> objects;
};

上述代码展示了对象位置的动态更新逻辑，通过哈希表实现 O(1) 级查找效率，确保高频调用下的响应速度。

粗检测流程

采用层次化剔除策略，先进行轴对齐包围盒（AABB）相交测试，过滤明显不相交的对象对。

• 遍历所有活动对象
• 构建当前帧的检测任务列表
• 执行批量粗检测并输出潜在碰撞对

3.3 层次包围体树（BVH）在 C++ 中的轻量实现

BVH 结构设计原则

层次包围体树通过递归划分空间加速碰撞检测与光线追踪。其核心在于构建紧凑的包围盒层级，降低场景遍历复杂度。

轻量节点定义

struct AABB {
    Vec3 min, max;
    bool intersects(const Ray& r) const;
};

struct BVHNode {
    AABB bounds;
    int left, right; // 子节点索引：负值表示叶子指向图元
    int splitAxis;   // 分割轴（0=x,1=y,2=z）
};

该结构采用数组存储节点，提升缓存友好性；叶子节点用负索引指向图元 ID，避免指针开销。

构建策略对比

• 中点分割：实现简单，适合均匀分布图元
• 表面面积启发式（SAH）：计算成本高但结构更优
• 均分法：平衡左右子树大小，适用于快速预览场景

第四章：复杂场景下的高精度检测实战

4.1 连续碰撞检测（CCD）防止高速物体穿透

在物理引擎中，高速移动的物体会因离散时间步长导致帧间穿透，传统离散碰撞检测（DCD）无法捕捉运动轨迹中的碰撞瞬间。连续碰撞检测（CCD）通过追踪物体在时间步内的运动路径，有效防止穿透问题。

CCD 核心原理

CCD 将物体运动视为连续过程，计算其在当前帧内的扫掠体积（Swept Volume），并检测该体积是否与其它物体相交。常见方法包括扫掠测试（Sweep Test）和时间步细分（Temporal Subdivision）。

实现示例：球体扫掠检测

// 球体沿方向 dir 移动，检测是否与静态平面碰撞
bool sweepSpherePlane(const Sphere& sphere, const Vec3& dir, const Plane& plane, float& outTime) {
    float denom = dot(plane.normal, dir);
    if (fabs(denom) < 1e-6) return false; // 平行无碰撞
    float t = (plane.distance - dot(plane.normal, sphere.center)) / denom;
    if (t >= 0 && t <= 1.0f) {
        outTime = t;
        return true;
    }
    return false;
}

该函数计算球体在单位时间内是否穿过平面，返回碰撞发生的时间点 t，用于调整位置以避免穿透。

性能对比

4.2 多点接触与碰撞响应的物理模拟

在复杂交互场景中，多点接触的物理模拟需精确处理多个接触点之间的力分布与响应。传统单点模型无法准确还原真实物理行为，因此引入基于冲量迭代的接触解算器成为关键。

接触点检测与法向力计算

系统首先通过 GJK 算法检测凸体间的穿透深度，并生成接触点集。每个接触点包含位置、法向与穿透深度信息。

struct Contact {
    vec3 position;
    vec3 normal;      // 碰撞法向
    float depth;      // 穿透深度
    float restitution;// 恢复系数
};

该结构体用于存储每个接触点的物理属性。法向用于分解冲量方向，深度参与约束求解，恢复系数决定反弹强度。

迭代式冲量解算

采用顺序脉冲法（Sequential Impulses）对多个接触点进行迭代求解，确保总动量守恒：

• 遍历所有接触点，计算相对速度在法向的投影
• 根据恢复系数生成目标分离速度
• 调整冲量值使当前接触点满足约束条件
• 重复迭代直至系统能量收敛

4.3 自定义材质交互与摩擦力模型集成

在物理仿真系统中，真实感的物体交互依赖于精确的材质响应。通过自定义材质属性，可为不同表面赋予独特的摩擦行为。

材质属性定义

每个材质类型需声明静态与动态摩擦系数，例如：

struct Material {
    float staticFriction;   // 静摩擦系数
    float dynamicFriction;  // 动摩擦系数
    std::string name;
};

上述结构体用于描述基础表面特性。静摩擦系数决定物体启动滑动的阈值，而动摩擦系数影响滑动过程中的减速行为。

摩擦力计算逻辑

两物体接触时，系统根据材质对查表或插值获取综合摩擦系数：

最终摩擦力由库仑摩擦模型计算：$ F = \mu \cdot N $，其中 $ \mu $ 为组合系数，$ N $ 为法向力。

4.4 实时调试可视化：绘制碰撞轮廓与法线方向

在物理仿真调试中，实时可视化碰撞轮廓与法线方向是定位问题的关键手段。通过图形化反馈，开发者能够直观判断物体间的交互是否符合预期。

绘制碰撞轮廓

使用调试渲染器遍历碰撞体的几何顶点，生成闭合线条表示轮廓。例如在 Unity 中可通过 Gizmos.DrawLine 实现：

void OnDrawGizmos() {
    Collider2D col = GetComponent();
    if (col is PolygonCollider2D poly) {
        Vector2[] points = poly.points;
        for (int i = 0; i < points.Length; i++) {
            Vector2 current = transform.TransformPoint(points[i]);
            Vector2 next = transform.TransformPoint(points[(i + 1) % points.Length]);
            Gizmos.DrawLine(current, next);
        }
    }
}

该代码将局部坐标转换为世界坐标并连接顶点，形成可视轮廓。

法线方向可视化

法线用于指示碰撞响应方向。通常以箭头形式从接触点向外绘制，长度标准化为单位值，便于观察朝向一致性。

总结

本文系统阐述了 C++ 物理引擎中碰撞检测的关键技术与实现细节。从基础的 AABB 与 SAT 算法，到 BVH 空间划分优化，再到 CCD 连续检测与材质交互模型，涵盖了从理论到实践的全流程。通过合理的架构设计与性能优化策略，可有效构建高精度的物理模拟系统，满足游戏及仿真领域对实时性与准确性的双重需求。