2D质量守恒

时间:2019-02-05 03:04:48

标签: java math physics

我有这个程序,球在屏幕周围弹跳。我希望这些球能够彼此弹起。

所以我有这个     boolean intersects(ball b2) {. . .} 方法,它检查b。与b2相交。如果b.intersects(x)返回true,则调用calcDX()calcDY()

void calcDX(Ball b2) {
    double b1Momentum = mass * dx;
    double b2Momentum = b2.mass * b2.dx;

    double b2FinalVelocity = (b1Momentum + b2Momentum + mass * (dx + b2.dx)) / (mass + b2.mass);

    dx = b2FinalVelocity - dx - b2.dx;
    b2.dx = b2FinalVelocity; //hi
}

void calcDY(Ball b2) {
    double b1Momentum = mass * dy;
    double b2Momentum = b2.mass * b2.dy;

    double b2FinalVelocity = (b1Momentum + b2Momentum + mass * (dy + b2.dy)) / (mass + b2.mass);

    dy = b2FinalVelocity - dy - b2.dy;
    b2.dy = b2FinalVelocity;
}

这是我正在使用的弹性碰撞方程,其中dxdyx的变化和y的变化。

所以我有这个循环(我知道,while(true),忽略它。)每秒更新60次,每次调用updateBallPositions()

@Override
public void run() {
    int ticksPerSecond = 60;
    long lastTime = System.nanoTime();
    double nanoSecondsPerTick = 1000000000.0 / ticksPerSecond; // How many nano-seconds in a tick
    double delta = 0.0;

    while (true) {  // main game loop
        long now = System.nanoTime();
        delta += (now - lastTime) / nanoSecondsPerTick;
        lastTime = now;

        while (delta >= 1) { // tick
            delta -= 1;

            updateBallPositions();
            f.panel.repaint();
        }
    }
}

private void updateBallPositions() {
    for (Ball b : balls) {
        b.x += b.dx;
        b.y += b.dy;

        for (Ball b2 : balls) {
            if(b != b2 & b.intersects(b2)) {
                b.calcDX(b2);
                b.calcDY(b2);
            }
        }

        if(b.isTouchingHorizontalWall(f))
            b.dy *= -1;
        if(b.isTouchingVerticalWall(f))
            b.dx *= -1;
    }
}

这一切都很棒,但是没有用。球经常粘在一起并永远碰撞,同时应允许球平稳弹跳。 (如果您需要更多代码,请这样说)

1 个答案:

答案 0 :(得分:1)

您的碰撞处理看起来不对。

让球A的质量为ma,速度矢量为VA0
令球B的质量为mb,速度矢量为VB0
它们在C点发生碰撞
要计算碰撞后的速度,必须使用能量守恒定律和动量守恒定律

考虑临时轴:
T-通过点C与两个球相切,并且通过径向R连接球中心

enter image description here

T上的速度投影是切向矢量(例如,VAt),轴R上的速度投影是径向矢量(例如,VAr

重要-碰撞期间切向动量不变:VAt1=VAt, VBt1=VBt

所以我们在T-R坐标系中有方程式

 ma * VAr + mb * VBr = ma * VAr1 + mb * VAr2
 ma * VA0^2 + mb*VB0^2 = ma * VA1^2 + mb*VB1^2   // I removed division by two
 but
 VA1^2 = VAr1^2 + VAt^2
 VB1^2 = VBr1^2 + VBt^2

最后是由两个未知数为VAr1, VBr1的两个方程组成的系统

 ma * VAr + mb * VBr = ma * VAr1 + mb * VAr2
 ma * VA0^2 + mb*VB0^2 = ma * (VAr1^2 + VAt^2) + mb*(VBr1^2 + VBt^2)

解决它,将速度分量转换为世界坐标系OXY