Python/Pygame中的弹性碰撞

Question

我在引力场中处理弹性碰撞时遇到了严重的问题。我试图按照能量守恒定律来实现它，但是结果出错，就像这个视频所示。首先，两个物体粘在一起，碰撞后几百帧之后，它们的速度变得非常快。 完整代码在此处在线上，但是负责在碰撞后提供输出速度的方法是world.py中的这个函数:

def collision(self, obj1, obj2):
    R = obj1.radius + obj2.radius   # code is used to "jump back in time" to avoid penetration when there's a collision
    dx = obj1.x - obj2.x            #
    dy = obj1.y - obj2.y            #
    K = math.hypot(dx, dy)          #   
    dvx = obj1.vx - obj2.vx         #   
    dvy = obj1.vy - obj2.vy         #
    dv = math.hypot(dvx, dvy)       #
    deltat = (R - K)/dv             #
    print dv
    print deltat
    obj1.x = obj1.rect.centerx = obj1.x - obj1.vx # *deltat
    obj2.x = obj2.rect.centerx = obj2.x - obj2.vx # *deltat
    obj1.y = obj1.rect.centery = obj1.y - obj1.vy # *deltat
    obj2.y = obj2.rect.centery = obj2.y - obj2.vy # *deltat
    dx = obj2.x - obj1.x
    dy = obj2.y - obj1.x
    delta = math.hypot(dx, dy)
    nx = dx/delta
    ny = dy/delta
    vx1bc = obj1.vx * nx
    vx2bc = obj2.vx * nx
    vy1bc = obj1.vy * ny
    vy2bc = obj2.vy * ny
    vx2ac = (obj2["energy_loss"]*(vx1bc - vx2bc) + vx1bc + (obj2["mass"]/obj1["mass"]*vx2bc))/((obj2["mass"]/obj1["mass"])+1)
    vy2ac = (obj2["energy_loss"]*(vy1bc - vy2bc) + vy1bc + (obj2["mass"]/obj1["mass"]*vy2bc))/((obj2["mass"]/obj1["mass"])+1)
    vx1ac = (vx1bc + obj2["mass"]/obj1["mass"]*vx2bc - obj2["mass"]/obj1["mass"]*vx2ac)*obj1["energy_loss"]
    vy1ac = (vy1bc + obj2["mass"]/obj1["mass"]*vy2bc - obj2["mass"]/obj1["mass"]*vy2ac)*obj1["energy_loss"]
    V1cx = obj1.vx * ny
    V1cy = obj1.vy * ny
    V2cx = obj2.vx * ny
    V2cy = obj2.vy * ny
    alfa = math.atan2(ny, nx)
    alfa_deg = math.degrees(alfa)
    v1a = math.hypot(vx1ac, vy1ac)
    v2a = math.hypot(vx2ac, vy2ac)
    obj1.vx = v1a*math.cos(alfa)+V1cx * math.sin(alfa)
    obj2.vx = v2a*math.cos(alfa)+V2cx * math.sin(alfa)
    obj1.vy = v1a*math.sin(alfa)+V1cx * math.cos(alfa)
    obj2.vy = v2a*math.sin(alfa)+V2cx * math.cos(alfa)

- Siekacz

---

首先，两个物体粘在一起，在碰撞几百帧后达到巨大的速度。向量？我也想使用它们，但Python没有任何标准的向量实现。 - Siekacz

你不必遵循Python标准库：你可以编写自己的代码！ - Gareth Rees

仍然无法解决问题。 - Siekacz

我沿着您的代码一直跟到 delta 声明之后，但是我不理解接下来发生了什么。 "energy_loss" 代表什么？此外，能否解释一下变量名称 nx、ny、vx1bc、vx2bc、vy1bc、vy2bc、vx2ac、vy2ac、vx1ac、vy1ac、V1cx、V1cy、V2cx、V2cy、v1a、v2a 的含义？ - Kevin

边界框是什么形状？是由“半径”属性建议的圆形吗？还是由“rect”建议的矩形？ - nneonneo

你也应该发布你的数学推导。在你的计算中考虑了哪些因素和输入变量？它与直接弹性碰撞有何不同？ - nneonneo

Answer 1

你的代码使用了很多不同的变量，它们的命名也不是很相关。这使得跟上所有内容变得困难，也让其他人难以看出你的程序在做什么。也许你应该将变量附加到列表中并访问它们？此外，我建议你查看一个很棒的物理教程：http://www.petercollingridge.co.uk/pygame-physics-simulation。

- Syklis