SFML中Verlet积分碰撞求解器失效：球体碰撞粘连求助

1. 碰撞响应迭代次数不足 #

Verlet积分依赖位置驱动，单次碰撞修正经常无法彻底解决帧内穿透，导致下一帧检测时球体仍处于重叠状态，表现为粘连。

• 解决方法：每帧执行3-5次碰撞修正迭代，逐步把球体推到分离状态：

// 每帧多次迭代碰撞修正
for (int iter = 0; iter < 5; ++iter) {
    for (size_t i = 0; i < balls.size(); ++i) {
        for (size_t j = i + 1; j < balls.size(); ++j) {
            Ball& a = balls[i];
            Ball& b = balls[j];
            sf::Vector2f delta = a.position - b.position;
            float distSq = delta.x*delta.x + delta.y*delta.y;
            float minDist = a.radius + b.radius;
            float minDistSq = minDist * minDist;
            
            if (distSq < minDistSq && distSq > 0) {
                float dist = sqrt(distSq);
                float overlap = 0.5f * (dist - minDist);
                sf::Vector2f correction = delta / dist * overlap;
                // 修正当前位置
                a.position -= correction;
                b.position += correction;
            }
        }
    }
}

2. 未同步修正旧位置 #

Verlet中物体的速度由当前位置和旧位置推导而来，如果只改当前位置，旧位置没同步调整，下一帧物体会“回弹”回重叠区，造成粘连。

• 解决方法：修正当前位置时同步更新旧位置，保持速度逻辑一致：

// 在位置修正后添加这两行
a.oldPosition -= correction;
b.oldPosition += correction;

3. 碰撞精度与阈值问题 #

距离计算的精度误差或未处理极近距离（比如distSq为0），会导致修正逻辑失效，球体卡住。

• 解决方法：给最小距离加个微小的epsilon值，避免除以零或精度问题：

const float EPSILON = 0.001f;
float minDistSq = (minDist + EPSILON) * (minDist + EPSILON);

4. 缺少碰撞弹跳逻辑 #

如果没给碰撞添加弹跳系数（restitution），物体碰撞后动能快速损耗，看起来像粘在一起。

• 解决方法：位置修正后，基于碰撞法线调整旧位置来模拟弹跳：

// 计算碰撞法线
sf::Vector2f normal = delta / dist;
// 推导Verlet下的速度
sf::Vector2f velA = a.position - a.oldPosition;
sf::Vector2f velB = b.position - b.oldPosition;
float relVel = (velA.x - velB.x) * normal.x + (velA.y - velB.y) * normal.y;

// 只有物体靠近时才处理弹跳
if (relVel > 0) continue;
float restitution = 0.8f; // 0是完全非弹性，1是完全弹性
float impulse = -(1 + restitution) * relVel / (1/a.mass + 1/b.mass);

// 通过旧位置调整来改变速度
a.oldPosition += normal * impulse / a.mass;
b.oldPosition -= normal * impulse / b.mass;

5. 时间步长不稳定 #

帧时间步长波动大，会导致Verlet积分的位置更新误差积累，引发穿透和粘连。

• 解决方法：用固定时间步长更新物理，分离渲染与物理逻辑：

const float FIXED_TIMESTEP = 1.0f / 60.0f;
float accumulator = 0.0f;

void update(float deltaTime) {
    accumulator += deltaTime;
    while (accumulator >= FIXED_TIMESTEP) {
        // 执行物理更新（包括碰撞迭代）
        updatePhysics(FIXED_TIMESTEP);
        accumulator -= FIXED_TIMESTEP;
    }
    // 可选：渲染时用accumulator/FIXED_TIMESTEP做位置插值
}

内容的提问来源于stack exchange，提问作者Dreadward07