Wie modelliert man hyperbolische Bahnen?

Ich programmiere einen Orbitalsimulator und benötige Hilfe bei der Modellierung hyperbolischer Umlaufbahnen (oder Trajektorien, wie auch immer Sie sie nennen möchten). Bisher kann ich elliptische Bahnen mit den Standardbahnelementen (Perigäum, Exzentrizität, große Halbachse, mittlere Anomalie usw.) modellieren; Ich kann die mittlere Anomalie in wenigen Sekunden berechnen und daraus die Position und Geschwindigkeit des Fahrzeugs ableiten. Gibt es eine Möglichkeit, dies für hyperbolische Bahnen zu tun? Wenn Sie nur mit dem Geschwindigkeits- und Positionsvektor eines Fahrzeugs beginnen, können Sie vorhersagen, wo es sich von jetzt an jede Sekunde befinden wird, während es sich in einer hyperbolischen Umlaufbahn / Flugbahn befindet?

Antworten (1)

Ja. (Es klingt, als hätten Sie es mit einem Zwei-Körper-Problem zu tun, das analytische Lösungen hat.)

Zunächst einmal, angesichts der Position, r , die Geschwindigkeit, v , sowohl in einem Koordinatensystem mit dem Zentralkörper im Ruhezustand als auch dem GM des Zentralkörpers, μ , können Sie leicht die spezifische Energie (Energie pro Masseneinheit) berechnen, die Ihnen sagt, ob die Flugbahn elliptisch (negative Energie) oder hyperbolisch (positive Energie) ist. Oder vielleicht parabolisch (genau Nullenergie), was aber nur von mathematischem Interesse ist. Aus der Energie können Sie die große Halbachse erhalten, a . Hier v = | v | und r = | r | :

E = v 2 2 μ r

a = μ 2 E

Sie propagieren wahrscheinlich elliptische Bahnen mit Sinus und Cosinus. Eine hyperbolische Trajektorie wird auf ähnliche Weise mit hyperbolischen Sinus und hyperbolischen Cosinus durchgeführt. Sie können den spezifischen Drehimpuls nur aus Position und Geschwindigkeit berechnen und daraus die Exzentrizität erhalten, e .

M = r × v

e = 1 + M M μ a

Dann eine einfache Form der Flugbahn in der x j Ebene, mit engster Annäherung an die + x Achse ist:

x = a ( e cosch τ )
j = a e 2 1 Sünde τ

wo τ ist die exzentrische Anomalie, bezogen auf die Zeit, wo τ = 0 und t = 0 ist am nächsten:

t = a 3 μ ( e Sünde τ τ )

Würde dies für eine mittlere Anomalie von mehr als 90 Grad funktionieren?
Die mittlere Anomalie ist in diesem Fall kein Winkel. Es ist ein anderes Tier, definiert als M = e Sünde τ τ (zu sehen in der Gleichung für die Zeit).