Soll ich die Gravitationskonstante mit der Skala ändern und warum führen Änderungen der fps- und Zeitskala dazu, dass meine Umlaufbahn bricht?

Ich benutze Unity3D und schreibe C#

• Erdmasse: 5.972e+22
• Erdgeschwindigkeit: [0,0,0]
• Erdradius: 63,71
• ISS-Masse: 4194,55
• Startgeschwindigkeit der ISS: [76.50, 0, 0]
• ISS-Höhe: 4
• Gravitationskonstante: 6,67408e-15

Alle Werte sind die realen Werte / 100, da ich die Dinge verkleinere. (z. B. Gravitationskonstante war 6,67408e-11 / 11 = 6,67408e-15.

Hier ist der Code...

void Gravity(Gravity body) {
    float dSquared = (body.transform.position - transform.position).sqrMagnitude;
    float force = (body.G * shipMass * body.planetMass) / dSquared;
    Vector3 forceDirection = (body.transform.position - transform.position).normalized;
    Vector3 forceVector = (forceDirection * force);
    shipVelocity += forceVector * Time.deltaTime;

    transform.position += shipVelocity * Time.deltaTime;
}

Damit dies funktioniert, muss ich die Gravitationskonstante auf 1,59508e-21 einstellen

Hier ist ein Bild der Umlaufbahn, die rosa Linie ist der vergangene, nicht der zukünftige Pfad (die angezeigten berechneten Umlaufbahnelemente sind höchstwahrscheinlich falsch):

Hier ist ein Bild der Umlaufbahn, nachdem ich die Zeit beschleunigt habe:

Bearbeiten 1:

Mit dieser Beschleunigungsformel und dem ersten Teil von Russell Borogoves Antwort konnte ich mit diesen Werten eine kreisförmige Umlaufbahn erreichen:

• Erdmasse: 5.972e+22
• Erdgeschwindigkeit: [0,0,0]
• Erdradius: 63,71
• ISS-Masse: 4194,55
• Startgeschwindigkeit der ISS: [0,0765, 0, 0]
• ISS-Höhe: 4
• Gravitationskonstante: 6,67408e-24 (verschoben -15 dann +2 Größenordnungen)

und dieser Code:

void Gravity(Gravity body) {
    float dSquared = (body.transform.position - transform.position).sqrMagnitude;
    float M = body.planetMass;
    Vector3 forceDirection = (body.transform.position - transform.position).normalized;
    float g = G * M / dSquared;
    Vector3 forceVector = (forceDirection * g);
    shipVelocity += forceVector * Time.deltaTime;
}

Bearbeiten 2: Gravitationsberechnungen mit Russell Borogoves Akkumulator-Skript unten hinzugefügt ...

void Update() {
    Gravity(orbitee);
    ApplyForce();
}

void ApplyForce() {
    transform.position += shipVelocity * Time.deltaTime;
    transform.Rotate(transform.forward * shipTorque * Time.deltaTime);
}

void Gravity(Gravity body) {
    accumulator += Time.deltaTime;
    while (accumulator > 0.0f) {
        float dSquared = (body.transform.position - transform.position).sqrMagnitude;
        float M = body.planetMass;
        Vector3 forceDirection = (body.transform.position - transform.position).normalized;
        float g = G * M / dSquared;
        Vector3 forceVector = (forceDirection * g);
        shipVelocity += forceVector * simtime;
        accumulator -= simtime;
    }
}

Bearbeiten 3: Aktualisiertes Skript für Orbitalelemente (wahrscheinlich nicht 100% korrekt)

    //Mass of satellite
    float m1 = shipMass;

    //Mass of planet
    float m2 = orbiting.planetMass;
    float M = m2;

    //Relative Position Vector
    Vector3 r = FindRelativePosition(orbitee.transform.position, transform);

    //Relative Velocity Vector
    Vector3 v = shipCurVel;
    //#!# add relative to planet velocity

    //Specific Angular Momentum
    Vector3 h = Vector3.Cross(r, v);

    //Standard Gravitational Parameter
    float µ = G * (m1 + m2);
    //float µ = G * M;

    //Eccentricity Vector
    Vector3 evec = (Vector3.Cross(v, h) / µ) - (r / Vector3.Magnitude(r));

    //Eccentricity
    float e = Vector3.Magnitude(evec);

    //Vector to Ascending Node #!#
    Vector3 n = Vector3.Cross(new Vector3(0, 0, 1), h);

    //True Anomaly
    float t = Mathf.Acos((Vector3.Dot(evec, r)) / (e * Vector3.Magnitude(r)));
    if (Vector3.Dot(r, v) < 0)
        t = (2 * Mathf.PI) - t;

    //Eccentric Anomaly
    float E = 2 * Mathf.Atan(Mathf.Tan(t / 2) / Mathf.Sqrt((1 + e) / (1 - e)));

    //Longitude of Ascending Node (2D)
    float Ω = 0;

    //Inclination (2D)
    float i = 0;

    //Argument of Periapsis
    float ω = Mathf.Atan2(evec.y, evec.x);
    float ωdegrees = ω * (180 / Mathf.PI);
    //float ω = Mathf.Acos((Vector3.Dot(n, evec)) / (e * Vector3.Magnitude(n)));

    //Mean Anomaly 
    float MA = E - (e * Mathf.Sin(E));

    //Semi-Major Axis
    float a = 1 / ((2 / Vector3.Magnitude(r)) - (Mathf.Pow(Vector3.Magnitude(v), 2) / µ));

    //Apoapsis
    float ap = a * (1 + e);

    //Periapsis
    float pe = a * (1 - e);

    //Orbital Period
    float T = (2 * Mathf.PI) * Mathf.Sqrt(Mathf.Pow(a, 3) / µ)/60;

Letzte Änderung: Problem behoben, das ich mit falscher Exzentrizität hatte - die relative Position wurde falsch berechnet.