Soal Tabrakan Bola 2d: tidak ada kekekalan energi

Saya mencoba menulis simulasi fisika sederhana di mana bola dengan jari-jari dan massa yang berbeda-beda memantul dalam lingkungan yang elastis sempurna dan tanpa gesekan. Saya menulis kode saya sendiri dengan mengikuti sumber ini: http://www.vobarian.com/collisions/2dcollisions2.pdf dan saya juga menguji kode dari sini: Ball to Ball Tabrakan - Deteksi dan Penanganan

PERTANYAAN DIEDIT

Dengan bantuan Rick Goldstein dan Ralph, kode saya berfungsi (ada kesalahan ketik..). Terima kasih banyak atas bantuan Anda. Namun saya masih bingung mengapa algoritma lain tidak bekerja untuk saya. Bola memantul ke arah yang benar, namun energi total sistem tidak pernah kekal. Kecepatannya semakin cepat hingga bola mulai berkedip dalam posisi statis di layar. Saya sebenarnya ingin menggunakan kode ini dalam program saya, karena lebih ringkas daripada yang saya tulis.

Berikut adalah algoritma fungsional yang saya tulis (walaupun saya mengambil bit pertama dari sumber lain itu). Ini ada di kelas Bubble:

public void resolveCollision(Bubble b)
{
    // get the minimum translation distance
    Vector2 delta = (position.subtract(b.position));
    float d = delta.getMagnitude();
    // minimum translation distance to push balls apart after intersecting
    Vector2 mtd = delta.multiply(((getRadius() + b.getRadius())-d)/d); 

    // resolve intersection --
    // inverse mass quantities
    float im1 = 1 / getMass(); 
    float im2 = 1 / b.getMass();

    // push-pull them apart based off their mass
    position = position.add(mtd.multiply(im1 / (im1 + im2)));
    b.position = b.position.subtract(mtd.multiply(im2 / (im1 + im2)));

    //get the unit normal and unit tanget vectors
    Vector2 uN = b.position.subtract(this.position).normalize();
    Vector2 uT = new Vector2(-uN.Y, uN.X);

    //project ball 1 & 2 's velocities onto the collision axis
    float v1n = uN.dot(this.velocity);
    float v1t = uT.dot(this.velocity);
    float v2n = uN.dot(b.velocity);
    float v2t = uT.dot(b.velocity);

    //calculate the post collision normal velocities (tangent velocities don't change)
    float v1nPost = (v1n*(this.mass-b.mass) + 2*b.mass*v2n)/(this.mass+b.mass);
    float v2nPost = (v2n*(b.mass-this.mass) + 2*this.mass*v1n)/(this.mass+b.mass);

    //convert scalar velocities to vectors
    Vector2 postV1N = uN.multiply(v1nPost);
    Vector2 postV1T = uT.multiply(v1t);
    Vector2 postV2N = uN.multiply(v2nPost);
    Vector2 postV2T = uT.multiply(v2t);

    //change the balls velocities
    this.velocity = postV1N.add(postV1T);
    b.velocity = postV2N.add(postV2T);
}

Dan inilah yang tidak berhasil

public void resolveCollision(Bubble b)
{
    // get the minimum translation distance
    Vector2 delta = (position.subtract(b.position));
    float d = delta.getMagnitude();
    // minimum translation distance to push balls apart after intersecting
    Vector2 mtd = delta.multiply(((getRadius() + b.getRadius())-d)/d); 

    // resolve intersection --
    // inverse mass quantities
    float im1 = 1 / getMass(); 
    float im2 = 1 / b.getMass();

    // push-pull them apart based off their mass
    position = position.add(mtd.multiply(im1 / (im1 + im2)));
    b.position = b.position.subtract(mtd.multiply(im2 / (im1 + im2)));

    // impact speed
    Vector2 v = (this.velocity.subtract(b.velocity));
    float vn = v.dot(mtd.normalize());

    // sphere intersecting but moving away from each other already
    if (vn > 0.0f) return;

    // collision impulse (1f is the coefficient of restitution)
    float i = (-(1.0f + 1f) * vn) / (im1 + im2);
    Vector2 impulse = mtd.multiply(i);

    // change in momentum
    this.velocity = this.velocity.add(impulse.multiply(im1));
    b.velocity = b.velocity.subtract(impulse.multiply(im2));
}

Beri tahu saya jika Anda menemukan sesuatu. Terima kasih