Gelöste Aufgaben/Kerb

Aufgabenstellung

Ein Satellit soll eine stabile Umlaufbahn um die Erde beschrieben.

Berechnen Sie mögliche Lösungen.

Lösung mit Maxima

Bei unserer Suche geht es um kartesische- und polare-Koordinatensysteme, Gleichgewichtsbedingungen formuliert mit dem Prinzip der virtuellen Verrückungen und er Stabilität von Lösungen.

tmp

Header

Ein Satellit im Orbit ist einer Erdbeschleunigung g ausgesetzt, die nichtlinear von seinem Abstand zum Erdmittelpunkt abhängt. Dieses nichtlineare Anfangswertproblem lösen wir hier als Anfangswertproblem numerisch.

tmp

Declarations

Wir brauchen die Systemparameter

Symbol	Bedeutung
R = 6378,137 km	Erdradius am Äquator
g₀ = 9.81 m/s2	Erdbeschleunigung an der Erdoberfläche
T = 1 d	Periodendauer einer Erdumdrehung (d .. day)
O_r = 35786 km	Höhe eines Satelliten über der Erdoberfläche in einem geostationären Orbit.

Kinematics

Wir führen ein kartesisches x, y, z - Koordinaten-System im Erdmittelpunkt (earth-centered inertial) ein, das nicht mit der Erde mitrotiert - wir setzen dieses als Intertialsystem an (obwohl wir wissen, dass die Erde um die Sonne rotiert und die Sonne selbst auch kein Intertialsystem in unserem Sinne ist).

Wir brauchen die Einheitsvektoren

{\vec{e}}_{x}, {\vec{e}}_{y}

entlang der x,y-Achsen, um den Ortsvektor zum Satelliten

\vec{r} = r (t) \cdot {\vec{e}}_{r} (t)

zu erfassen.

Dafür verwenden wir

{\vec{e}}_{r} (t) = \cos (φ (t)) \cdot {\vec{e}}_{x} + \sin (φ (t)) \cdot {\vec{e}}_{y}

.

Die Bewegung des Satelliten können wir also auch in kartesischen Koordinarten

\begin{array}{l} u_{x} (t) = r (t) \cdot \cos (φ (t)) \\ u_{y} (t) = r (t) \cdot \sin (φ (t)) \end{array}

,

anschreiben, was hier allerdings wenig anschaulich ist! Für die Variation von u_x, u_y erhalten wir

\begin{array}{l} δ u_{x} = \cos (φ (t)) \cdot δ r - r (t) \cdot \sin (φ (t)) \cdot δ φ \\ δ u_{y} = \sin (φ (t)) \cdot δ r + r (t) \cdot \cos (φ (t)) \cdot δ φ \end{array}

tmp

Equilibrium Conditions

Die Bewegungsgleichungen schreiben wir mit dem Prinzip der virtuellen Verrückungen an, hier mit den Anteilen der D'Alembert'sche Trägheitskräfte und der Gewichtskraft

δ W = \underset{d'Alembert'sche Trägheitskräfte}{\underset{⏟}{- m \cdot {\ddot{u}}_{x} (t) \cdot δ u_{x} - m \cdot {\ddot{u}}_{y} (t) \cdot δ u_{y}}} \underset{Gewichtskraft}{\underset{⏟}{- m g (r) \cdot δ r}}

,

dabei ist m die Masse des Satelliten - die sich aber sofort herauskürzt.

In die Gleichgewichtsbedingung

δ W \overset{!}{=} 0

setzen wir alle Größen ein und finden

δ W = - \frac{(g_{0} m R^{2} + m \cdot r^{2} \ddot{r} (t) - m \cdot r^{3} {\dot{φ}}^{2}) δ r + (2 m \cdot r^{3} \dot{φ} \dot{r} + m \cdot r^{4} \ddot{φ}) δ φ}{r^{2}}

.

Das Trennen nach den virtuellen Verrückungen liefert dann die beiden Bewegungsgleichungen

\begin{array}{cl} - r \cdot {\dot{φ}}^{2} + \ddot{r} + g_{0} \cdot {(\frac{R}{r})}^{2} & = 0, \\ - 2 r \cdot \dot{φ} \cdot \dot{r} - r^{2} \cdot \ddot{φ} & = 0 \end{array}

.

Für die dimensionslose Schreibweise setzen wir

\begin{array}{ll} t & = τ \cdot T \\ r (t) & = ρ (t) \cdot R \end{array}

an und verwenden die Abkürzungen

\begin{array}{ll} \frac{d}{d τ} φ & = Φ \\ \frac{d}{d ρ} φ & = P \end{array}

,

wobei das P für das große griechische ρ steht. Damit - und mit dem dimensionslosen Parameter μ statt

g = \frac{μ R}{T^{2}}

lesen sich die beiden Bewegungsgleichungen so:

\begin{array}{ll} \frac{d}{d τ} P & = \frac{Φ^{2} ρ^{3} - μ}{ρ^{2}} \\ \frac{d}{d ρ} Φ & = - 2 \frac{Φ P}{ρ} \end{array}

.

tmp

Stationary Solution

..

tmp

Solving

..

tmp

Post-Processing

..

Links

...

Literature

...

Trajectory animation