Berechnung der spezifischen Umlaufenergie, der großen Halbachse und der Umlaufzeit eines umlaufenden Körpers

Ist es möglich, die spezifische Orbitalenergie zu berechnen? ϵ , die große Halbachse A , und die Umlaufzeit T (oder P ) ohne dass Ihnen einer davon zur Verfügung steht? Die Werte, die mir zur Verfügung stehen, sind die Geschwindigkeit des umkreisenden Körpers relativ zum Schwerpunkt, seine aktuelle Position (ebenfalls relativ zum Schwerpunkt) und die zentrale Masse, die die Quelle der Schwerkraft darstellt M . Ich habe auch die Masse des umkreisenden Körpers, aber sie ist vernachlässigbar.

Ist es angesichts all dieser Dinge und ohne äußere Faktoren möglich, einen der oben aufgeführten Werte zu berechnen? Nach dem dritten Keplerschen Gesetz ist die Umlaufzeit im Verhältnis gegeben 4 π 2 / T 2 = G M / R 3 Wo T ist die Umlaufzeit, G ist Newtons universelle Gravitationskonstante, M ist die Masse des größeren Körpers (da die Masse des umkreisenden Körpers vernachlässigbar ist) und R ist der Abstand zwischen dem Schwerpunkt und dem umlaufenden Körper. Dies hilft nicht sehr, da es sich um eine Proportion handelt und nicht mit Algebra umgangen werden kann, um einen echten Wert zu erhalten T (Ich finde?).

Wie auch immer, ich habe das Internet und Wikipedia durchsucht, um einen Weg zu finden, diese Werte zu berechnen, aber ich bin ratlos. Ich versuche zu sehen, ob es eine Möglichkeit gibt, diese Dinge für ein kleines Programmierprojekt / einen kleinen Simulator zu berechnen. Andernfalls müsste das Programm eine Zeit lang simuliert werden, um einen dieser Werte zu ermitteln und die anderen zu berechnen.

Ich gebe Ihnen nur einen ersten Hinweis: Die Gesamtenergie der Umlaufbahn ist gerecht 1 2 M v 2 G M M R . Für den Rest gibt es eine genau bekannte Lösung für das Kepler-Problem, die im Internet verfügbar sein sollte. Passen Sie dies an Ihre Anfangsbedingungen an, und Sie sind goldrichtig.
Nun, ich kann nicht glauben, dass ich diese Formel für die Energie einer Umlaufbahn verpasst habe. Ich bin immer wieder auf Formulierungen gestoßen, die sich auf die große Halbachse beziehen. Nun, um das Internet nach der Lösung zu durchsuchen, von der Sie sprechen ...
Interessante Situation - wenn Ihre "Geschwindigkeit", wie ich vermute, ein Skalar und keine Vektorgröße ist, haben Sie nicht genügend Informationen, um die vollständige Umlaufbahn zu bestimmen (keine Exzentrizitätsinformationen), aber Sie haben immer noch genug, um die gewünschten Mengen zu erhalten . Kleiner Punkt: Das ist die neue und verbesserte Version von Keplers Gesetz, die Ihnen von Newton präsentiert wurde – es ist eine Gleichheit, keine Proportion.
Die Geschwindigkeit ist durch einen Vektor gegeben, ja. Danke aber für die Kleinigkeit! Das wusste ich (offensichtlich) nicht.

Antworten (1)

Ja, Sie können alle diese Größen ableiten. Die spezifische Orbitalenergie E Ist

E = 1 2 v 2 μ R = μ 2 A ,
Wo μ = G M 3 / ( M + M ) 2 , Und A ist die große Halbachse. Die Umlaufzeit folgt aus Keplers drittem Gesetz:
T 2 = ( 2 π ) 2 A 3 μ .
Wenn Sie auch die Radialgeschwindigkeit kennen v R und die Tangentialgeschwindigkeit v T getrennt bei R , dann kannst du auch den spezifischen relativen Drehimpuls berechnen H und die orbitale Exzentrizität e :
H 2 = R 2 v T 2 = μ A ( 1 e 2 ) .


Bearbeiten

Mehrere Leute haben versucht, sich zu ändern μ hinein μ = G ( M + M ) . Das ist falsch, denn das ist die Formel für Relativbewegung statt Bewegung in Bezug auf den Massenmittelpunkt . Die Bewegungsgleichungen des Zweikörperproblems lauten

M R ¨ M = G M M | R M R M | 3 ( R M R M ) , M R ¨ M = G M M | R M R M | 3 ( R M R M ) ,
Wo R M Und R M sind die Positionen des kleinen und großen Körpers in Bezug auf den Massenmittelpunkt. Was wir wollen, ist die Bewegung des kleinen Körpers in Begriffen von auszudrücken R M . Per Definition bleibt die Lage des Massenmittelpunkts konstant,
M R M + M R M = 0 ,
so dass
R M R M = M + M M R M .
Deshalb,
M R ¨ M = G M M M 3 ( M + M ) 3 R M 3 ( M + M M R M ) ,
oder
R ¨ M = μ R M 3 R M ,
mit μ = G M 3 / ( M + M ) 2 . In meiner Antwort R = R M . Ich hoffe, das klärt die Dinge auf.

Nochmals vielen Dank für die Hilfe, Pulsar. Ich schätze die Formeln!
Sie können sich auch ausdrücken A Und e in Abhängigkeit von Ihren vorgegebenen Werten v Und R :
A = μ R 2 μ R v 2
e = μ 2 + cos 2 ϕ R v 2 ( R v 2 2 μ ) μ
Wo ϕ ist der Winkel der Geschwindigkeit relativ zum Vektor senkrecht zum Positionsvektor.