Wie berechnet man den Geschwindigkeitsvektor bei einer hyperbolischen Umlaufbahn?

Das Problem

Ich versuche, eine Formel zu bekommen, um die Zustandsvektoren zu berechnen R Und v = R ˙ auf einer Umlaufbahn, bei einer echten Anomalie v . Ich folge dem hier beschriebenen Prozess: https://downloads.rene-schwarz.com/download/M001-Keplerian_Orbit_Elements_to_Cartesian_State_Vectors.pdf . Der erste Berechnungsschritt beinhaltet die Berechnung der einfachen Zwischenzustandsvektoren Ö Und Ö ˙ Verlegung in der xy-Ebene (die dann im Raum gedreht werden, um zu bekommen R Und v bzw) :

Ö = R ( cos v Sünde v 0 ) ,     Ö ˙ = μ A R ( Sünde E 1 e 2 cos E 0 )

mit R = P 1 + e cos v .

Dies funktioniert gut bei einer elliptischen Umlaufbahn, ist aber bei einer hyperbolischen deswegen ungültig e > 1 Und A < 0 ergebend 1 e 2 Und μ A undefiniert sein.

Meine jetzige Lösung

Im Fall einer hyperbolischen Umlaufbahn habe ich die zweite Antwort aus diesem Beitrag angepasst: Berechnung des Geschwindigkeitszustandsvektors mit Umlaufbahnelementen in 2D zu berechnen Ö ˙ mit dem Flugbahnwinkel ϕ , den Drehimpuls kennen H = | | H | | . Wir berechnen zuerst den radialen Einheitsvektor u Ö ^ des Zwischenpositionsvektors und u S ^ der Einheitsvektor senkrecht zu u Ö ^ in der xy-Ebene:

u Ö ^ = Ö R = ( u Ö ,   X u Ö ,   j 0 ) ,     u S ^ = ( u Ö ,   j u Ö ,   X 0 )

Wir berechnen dann sin und cos des Flugbahnwinkels:

cos ϕ = H R v ,     Sünde ϕ = e Sünde v 1 + e cos v cos ϕ

mit v ist die Größe der Geschwindigkeit, berechnet aus der vis-viva-Gleichung . Und schließlich erhalten wir den Zwischengeschwindigkeitsvektor:

Ö ˙ = v ( Sünde ( ϕ ) u Ö ^ + cos ( ϕ ) u S ^ )

Eine bessere Lösung?

Gibt es einen besseren, einfacheren Weg, um diesen Zwischengeschwindigkeitsvektor im Fall einer hyperbolischen Umlaufbahn zu berechnen? Eine, die kein Wissen erfordert H . Gibt es beispielsweise eine ähnliche Formel wie im PDF, die die hyperbolische exzentrische Anomalie verwendet? H ?

Vielen Dank im Voraus.

@uhoh Danke für deinen Kommentar! Ich habe die Mathematik etwas tiefer ausgearbeitet und tatsächlich eine Lösung gefunden. Ich werde es in einer Antwort sehr bald beschreiben.
Das ist ausgezeichnet!

Antworten (1)

Lösung für hyperbolische Geschwindigkeit

Nach einigen mathematischen Manipulationen fand ich schließlich eine tatsächliche Lösung, die die hyperbolische Anomalie nutzt H .

Im Folgenden, e ist die Exzentrizität der Umlaufbahn, v ist die wahre Anomalie und A ist die große Halbachse.

Voraussetzung: Beweisen ich H = E in einer hyperbolischen Umlaufbahn

Dieser kleine Beweis soll nur zeigen, wie man die Gleichheit aus den bekannten Formeln der exzentrischen Anomalien herausholen kann.

  • Für eine Ellipsenbahn ( e < 1 ), die exzentrische Anomalie E ist definiert durch:

    (1) bräunen E 2 = 1 e 1 + e bräunen v 2

  • Für eine hyperbolische Umlaufbahn ( e > 1 ), die hyperbolische Anomalie H (auch geschrieben F ) ist definiert durch:

    (2) Tanh H 2 = e 1 e + 1 bräunen v 2

Bei einer hyperbolischen Umlaufbahn 1 e < 0 führt zu einer undefinierten Definition von E in (1) wegen des Quadratwurzelterms. Daher die Notwendigkeit, sein hyperbolisches Äquivalent (2) zu verwenden. Betrachtet man jedoch die Beziehung (1) in C durch Einbeziehung ich = 1 ermöglicht eine komplexe Definition von E :

bräunen E 2 = ich e 1 e + 1 bräunen v 2

wobei wir den rechten Term von (2) bemerken. Dies ist tatsächlich direkt verknüpft E C Und H R von:

(3) bräunen E 2 = ich Tanh H 2

Die Beziehungen zwischen hyperbolischen und trigonometrischen Funktionen ergeben:

X R       ich Tanh ( X ) = bräunen ( ich X )

was, wenn es auf (3) angewendet wird, zu Folgendem führt:

bräunen E 2 = bräunen ich H 2

Und da X bräunen ( ich X ) ist bijektiv X R weil es proportional ist X Tanh ( X ) , wir folgern daraus:

(4) ich H = E

im Fall einer hyperbolischen Umlaufbahn mit H R (und deshalb E ich R ).

Anpassung der Zwischengeschwindigkeitsvektorformel an hyperbolische Bahnen

Die von René Schwarz beschriebene Gleichung zur Berechnung des Zwischengeschwindigkeitsvektors (ohne Berücksichtigung der z-Komponente mit dem Wert 0) lautet:

(5) Ö ˙ = μ A R ( Sünde E 1 e 2 cos E )

Wir nehmen daher jetzt eine hyperbolische Umlaufbahn an e > 1 Und A < 0 . Daher kann (5) nicht direkt verwendet werden, weil 1 e 2 Und μ A sind undefiniert in R . Indem man die Tatsache nutzt, dass A = | A | Und 1 e 2 = ( e 2 1 ) , unter Berücksichtigung der Gleichung in C gibt:

( 5 )         Ö ˙ = μ | A | R ( Sünde E ( e 2 1 ) cos E ) = ich μ | A | R ( Sünde E ich e 2 1 cos E ) = μ | A | R ( ich Sünde E e 2 1 cos E ) = μ | A | R ( Sünde ich E e 2 1 cosch ich E )

Weil X R     ich Sünde X = Sünde ich X Und cos X = cosch ich X .
Endlich einbeziehen ich H = E ich E = H , und die Tatsache, dass cosch ist sogar und Sünde ungerade ist , erhalten wir:

(6) Ö ˙ = μ A R ( Sünde H e 2 1 cosch H )

mit | A | geschrieben als A zur Klarheit.

Diese Formel scheint in praktischen Fällen zu funktionieren (Bahnsimulation und -bestimmung). Bitte zögern Sie nicht, einen Kommentar abzugeben, um eventuelle Fehler zu korrigieren.