# Bobylew-Steklow-Lösung

Die Bobylew-Steklow-Lösung (englisch case of Bobylev–Steklov)$^{:1173–1185}$$^{:92}$ ist in der Kreiseltheorie eine von Wladimir Steklow und Dmitri Bobyljow, der in der Schreibung Bobylew publizierte, 1896 unabhängig voneinander gefundene Lösung der Euler-Poisson-Gleichungen, die die Drehung von schweren Kreiseln um einen festen Punkt bestimmen.
Bobyljow und Steklow zeigten, dass bei einem unsymmetrischen Kreisel mit Hauptträgheitsmomenten A, B und C, bei dem C = 2A ist und der Massenmittelpunkt auf der zu C gehörenden Hauptträgheitsachse (kurz C-Achse) liegt, die Euler-Poisson-Gleichungen integriert werden können, wenn der Kreisel anfänglich weder eine Winkelgeschwindigkeit noch eine Winkelbeschleunigung um die B-Achse aufweist.
Ein solcher Kreisel kommt nie in Drehung um die B-Achse, hat eine konstante Winkelgeschwindigkeit um die C-Achse und alle seine anderen Zustandsgrößen sind Jacobische elliptische Funktionen der Zeit. Der schwere unsymmetrische Kreisel führt analytisch darstellbare periodische Präzessionsbewegungen aus. Obwohl das Hauptträgheitsmoment B nicht in die analytische Beschreibung eingeht, beeinflusst es die Stabilität der entstehenden Bewegung.$^{:1180ff.}$
Steklow fand drei Jahre später (1899) eine weitere Lösung, die 1952 von P. A. Kuz'min$^{:96}$ oder 2005 von A. P. Markeev ergänzt wurde.

## Bobylew-Steklow-Lösung
### Lösung der Bewegungsgleichungen

Der Kreisel unterliegt in der Bobylew-Steklow-Lösung ausschließlich der Schwerkraft, womit seine Bewegung von den Euler-Poisson-Gleichungen bestimmt wird. Diese Gleichungen besitzen zwei Integrale der Bewegung, die Gesamtenergie und den Drehimpuls in Lotrichtung. Im Bobylew-Steklow-Fall sind weiterhin die Winkelgeschwindigkeiten um die 2- und die 3-Achse konstant. Diese Achsen besitzen in vertikaler z-Richtung Komponenten n$_{2,3}$, die sich als Funktionen der z-Komponente n$_1$ der 1-Achse erweisen. Für n$_1$ existiert infolgedessen eine mit elliptischen Funktionen erfüllbare autonome Differentialgleichung, mit deren Lösung alle weiteren Zustandsgrößen des Kreisels dargestellt werden können.

#### Euler-Poisson-Gleichungen

Die Euler-Poisson-Gleichungen lauten im Bobylew-Steklow-Fall:
    ${\begin{aligned}A{\dot {p}}=&(B-2A)qr+c_{0}n_{2}\\B{\dot {q}}=&Arp-c_{0}n_{1}\\2A{\dot {r}}=&(A-B)pq\\{\dot {n}}_{1}=&rn_{2}-qn_{3}\\{\dot {n}}_{2}=&pn_{3}-rn_{1}\\{\dot {n}}_{3}=&qn_{1}-pn_{2}\end{aligned}}$
Darin sind n$_{1,2,3}$ die Koordinaten des antiparallel zur Gewichtskraft nach oben weisenden Einheitsvektors und p, q, r = ω$_{1,2,3}$ die Winkelgeschwindigkeiten im Hauptachsensystem. Das Stützpunktmoment c$_0$ = mgs ist das Produkt aus der Gewichtskraft mg mit dem Abstand s des Massenmittelpunkts vom Stützpunkt entlang der 3-Achse. Der Überpunkt bildet die Zeitableitung.
Wie bei jedem schweren Kreisel sind in der Bobylew-Steklow-Lösung die Gesamtenergie und der Drehimpuls L$_z$ in Lotrichtung konstant. Der Präzessionswinkel um die Vertikale kann aus diesen Gleichungen nicht ermittelt werden, was im Hauptartikel nachzuschlagen ist.

#### Anfangsbedingungen

Mit den Anfangsbedingungen $q={\dot {q}}=0$ spezialisieren sich die Euler-Poisson-Gleichungen zu:
    ${\begin{aligned}A{\dot {p}}=&c_{0}n_{2}&(1)\\B{\dot {q}}=&Arp-c_{0}n_{1}=0&(2)\\2A{\dot {r}}=&0&(3)\\{\dot {n}}_{1}=&rn_{2}&(4)\\{\dot {n}}_{2}=&pn_{3}-rn_{1}&(5)\\{\dot {n}}_{3}=&-pn_{2}&(6)\end{aligned}}$
Die Zeitableitung der Gleichung (2) verschwindet wegen der Bedingungen (1), (3) und (4), weswegen (2), (3) und q = 0 permanent gilt. Die sechs Gleichungen sind dadurch nicht nur anfänglich, sondern dauerhaft in Kraft. Aus (3) ist ersichtlich, dass r konstant ist. Bei r = 0 ergibt sich eine pendel­artige Drehung um die horizontale erste Hauptachse,$^{:1177}$ was hier nicht interessiert. Hier wird
    $r=r_{0}={\sqrt {\frac {c_{0}}{A}}}\rho \neq 0$
angenommen. Der an dieser Stelle eingeführte dimensionslose Parameter $\rho :=r_{0}{\sqrt {\tfrac {A}{c_{0}}}}$ ist proportional zum Drehimpuls um die 3-Achse.

#### Herleitung des Elliptischen Integrals

Mit den Gleichungen (6), (2) und (4) wird n$_3$ berechnet:
    ${\dot {n}}_{3}=-pn_{2}=-{\frac {c_{0}n_{1}}{Ar_{0}}}{\frac {{\dot {n}}_{1}}{r_{0}}}=-{\frac {(n_{1}^{2}{\dot\;}}}{2\rho ^{2}}}\quad \rightarrow \quad n_{3}={\frac {1}{2\rho ^{2}}}(j\rho -n_{1}^{2})$
Die Integrationskonstante j ist dimensionslos und proportional zum lotrechten Drehimpuls L$_z$. Die fünfte Gleichung liefert mit dem Gefundenen
    ${\dot {n}}_{2}=r_{0}\left({\frac {j}{2\rho ^{3}}}-1\right)n_{1}-{\frac {r_{0}}{2\rho ^{4}}}n_{1}^{3}$
und die zeitabgeleitete Gleichung (4) eine autonome Differentialgleichung
    ${\ddot {n}}_{1}=r_{0}{\dot {n}}_{2}=r_{0}^{2}\left({\frac {j}{2\rho ^{3}}}-1\right)n_{1}-{\frac {r_{0}^{2}}{2\rho ^{4}}}n_{1}^{3}$
in n$_1$. Multiplikation mit $2{\dot {n}}_{1}$ erlaubt eine Zeitintegration:
    ${\dot {n}}_{1}^{2}=D+r_{0}^{2}\left({\frac {j}{2\rho ^{3}}}-1\right)n_{1}^{2}-{\frac {r_{0}^{2}}{4\rho ^{4}}}n_{1}^{4}$
Wegen ${\dot {n}}_{1}=r_{0}n_{2},\,n_{3}={\tfrac {1}{2\rho ^{2}}}(j\rho -n_{1}^{2}),\,n_{1}^{2}+n_{2}^{2}+n_{3}^{2}=1$ ist die Integrationskonstante D nicht beliebig, sondern es ergibt sich nach algebraischen Umformungen
    ${\dot {n}}_{1}=\pm {\frac {r_{0}}{2\rho ^{2}}}{\sqrt {4\rho ^{4}-j^{2}\rho ^{2}+(2j\rho -4\rho ^{4})n_{1}^{2}-n_{1}^{4}}}$
Trennung der Veränderlichen führt auf ein elliptisches Integral
    $\pm {\frac {r_{0}}{2\rho ^{2}}}\int \mathrm {d} t=\int {\frac {\mathrm {d} n_{1}}{\sqrt {4\rho ^{4}-j^{2}\rho ^{2}+(2j\rho -4\rho ^{4})n_{1}^{2}-n_{1}^{4}}}}$
womit das Problem im Prinzip gelöst ist. Die darin vorkommenden Werte
    $\rho ={\sqrt {\frac {A}{c_{0}}}}r_{0},\quad j={\frac {L_{z}}{\sqrt {Ac_{0}}}}=2\rho \left(n_{3}+{\frac {n_{1}^{2}}{2\rho ^{2}}}\right)$
sind Konstanten und liegen mit den Anfangsbedingungen vor, in denen
    $q=0,\quad {\dot {q}}={\frac {1}{B}}(Arp-c_{0}n_{1})=0$
eingestellt sein muss.
Aus den zusammen getragenen Beziehungen ergibt sich die Gesamtenergie des Kreisels zu
    $E:={\frac {1}{2}}(Ap^{2}+Cr^{2})+c_{0}n_{3}={\frac {2\rho ^{3}+j}{2\rho }}c_{0}$.

### Darstellung mit Jacobis elliptischen Funktionen

Das Polynom unter der Wurzel im elliptischen Integral für n$_1$ hat die Nullstellen
    ${n_{1}}_{1,2,3,4}=\pm {\sqrt {\rho j-2\rho ^{4}\pm 2\rho ^{2}{\sqrt {\rho \left(\rho ^{3}+{\frac {1}{\rho }}-j\right)}}}}$
Je nachdem zwei oder vier dieser Nullstellen reell sind, entstehen unterschiedliche Lösungsfunktionen:
* Region I mit -2ρ < j < 2ρ erlaubt nur zwei reelle Nullstellen.
* Region II mit $\rho <1,\;2\rho <j<\rho ^{3}+{\tfrac {1}{\rho }}$ ermöglicht vier reelle Nullstellen.
Die elliptischen Moduli der elliptischen Lösungsfunktionen sind im Bild dargestellt. Bei k nahe 0 weichen die elliptischen Funktionen nur gering vom Sinus und Cosinus ab. Die weiß belassenen Flächen sind dem Kreisel nicht zugänglich, wenn $n_{3}={\tfrac {1}{2\rho ^{2}}}(j\rho -n_{1}^{2})$, eine Auflage, die aus den speziellen Anfangsbedingungen resultiert. Vergleichbare Lösungen entstehen, wenn ρ und j beide gleichzeitig ihr Vorzeichen wechseln. Der Einfachheit halber werden im Folgenden positive Werte für ρ angenommen.

#### Region I

In der Region I ist |j| < 2ρ und die analytische Lösung lautet:
    ${\begin{aligned}p=&{\hat {p}}\operatorname {cn} (u,k),&n_{1}=&{\hat {n}}_{1}\operatorname {cn} (u,k)\\q=&0,&n_{2}=&-2dk\operatorname {sn} (u,k)\operatorname {dn} (u,k)\\r=&r_{0},&n_{3}=&{\frac {j}{2\rho }}-2dk^{2}\operatorname {cn} ^{2}(u,k)\end{aligned}}$
mit
    ${\begin{aligned}{\hat {p}}=&2k\omega ,&{\hat {n}}_{1}=&{\sqrt {\frac {A}{c_{0}}}}\rho {\hat {p}}\\u=&\pm \omega (t-t_{0})&\omega =&{\sqrt {{\frac {c_{0}}{A}}d}},\\k=&{\sqrt {{\frac {1}{2}}\left(1+{\frac {j-2\rho ^{3}}{2\rho d}}\right)}},&d=&{\sqrt {\rho \left(\rho ^{3}+{\frac {1}{\rho }}-j\right)}}\end{aligned}}$
Die Konstanten j und ρ sind bei der #Herleitung des Elliptischen Integrals angegeben und sn, cn sowie dn sind drei grundlegende Jacobische elliptische Funktionen. Die Winkelgeschwindigkeit p hat die Periode ${\tfrac {4}{\omega }}K(k)$, worin K(k) das vollständige elliptische Integral der I. Art ist. Bei dieser Lösung schwingen p sowie n$_{1,2}$ periodisch um 0.
Im Grenzfall j = -2ρ oder j = 2ρ > 2 ist k = 0 und es stellen sich Staude-Drehungen mit der Winkelgeschwindigkeit $r={\sqrt {\tfrac {c_{0}}{A}}}\rho$ um die vertikale 3-Achse mit n$_3$ = -1 bzw. n$_3$ = +1 ein.$^{:1176}$
Falls j = 2ρ < 2 ist k = 1 und es ergibt sich ein mit der Region II gemeinsamer #Aperiodischer Grenzfall.

#### Region II

In der Region II ist ρ < 1, 2ρ < j < ρ$^3$ + 1/ρ und die analytische Lösung lautet:
    ${\begin{aligned}p=&{\hat {p}}\operatorname {dn} (u,k),&n_{1}=&{\hat {n}}_{1}\operatorname {dn} (u,k)\\q=&0,&n_{2}=&-2d\operatorname {sn} (u,k)\operatorname {cn} (u,k)\\r=&r_{0},&n_{3}=&{\frac {j}{2\rho }}-{\frac {2d}{k^{2}}}\operatorname {dn} ^{2}(u,k)\end{aligned}}$
mit
    ${\begin{aligned}{\hat {p}}=&\pm 2\omega ,&{\hat {n}}_{1}=&\pm \rho ^{2}{\frac {\hat {p}}{r_{0}}}=\pm {\sqrt {\frac {2\rho (j-2\rho ^{3})}{2-k^{2}}}}\\u=&\omega (t-t_{0}),&\omega =&{\sqrt {\frac {j-2\rho ^{3}}{2\rho ^{3}(2-k^{2})}}}r_{0}\\k=&{\sqrt {\frac {4\rho d}{j-2\rho ^{3}+2\rho d}}},&d=&{\sqrt {\rho \left(\rho ^{3}+{\frac {1}{\rho }}-j\right)}}\end{aligned}}$
Die Konstanten j sowie ρ sind bei der #Herleitung des Elliptischen Integrals angegeben, sn, cn sowie dn sind drei grundlegende Jacobische elliptische Funktionen. Die Winkelgeschwindigkeit p hat die Periode ${\tfrac {2}{\omega }}K(k)$, worin K(k) das vollständige elliptische Integral der I. Art ist. Bei dieser Lösung wechseln p sowie n$_1$ nie das Vorzeichen und schwingen periodisch um einen von 0 verschiedenen Wert.
Bei j = ρ$^3$ + 1/ρ ist k = 0 und der Kreisel führt eine Staude-Drehung um die Vertikale mit n$_1$ = ±√(1 - ρ$^4$), n$_2$ = 0 und n$_3$ = ρ$^2$ aus. Falls j = 2ρ kommt es zum aperiodischen Grenzfall.

#### Aperiodischer Grenzfall

Den beiden Regionen I und II ist der Grenzfall j = 2ρ < 2 gemeinsam. Dort ist k = 1 und die Jacobischen elliptischen Funktionen werden zu den aperiodischen Hyperbelfunktionen:
    ${\begin{aligned}\\p=&{\frac {2\omega }{\cosh(u)}},&n_{1}=&{\frac {2\rho {\sqrt {1-\rho ^{2}}}}{\cosh(u)}}\\q=&0,&n_{2}=&-2(1-\rho ^{2}){\frac {\sinh(u)}{\cosh ^{2}(u)}}\\r=&r_{0},&n_{3}=&1-{\frac {2(1-\rho ^{2})}{\cosh ^{2}(u)}}\end{aligned}}$
mit
    $u=\omega (t-t_{0})\quad {\text{und}}\quad \omega ={\sqrt {\frac {c_{0}}{A}}}{\sqrt {1-\rho ^{2}}}$
Die Bewegung geht asymptotisch in eine Staude-Drehung mit der Winkelgeschwindigkeit $r={\sqrt {\tfrac {c_{0}}{A}}}\rho$ um die vertikale 3-Achse über.$^{:1176}$

### Stabilität

Die Stabilität der Bewegung hängt nicht nur von den Parametern ρ und j ab, sondern auch vom mittleren Hauptträgheitsmoment B. Dieses liegt bei einem realen Kreisel im Intervall A < B < 3A. Die Zonen stabilen und instabilen Verhaltens sind im Parameterraum ρ, j, β = B/A durch dreidimensional gekrümmte Flächen getrennt. In Region I sind Kreisel mit zu großem B > A oft instabil. In Region II sind Kreisel mit B nahe 2A oder ρ > 0,8 oft instabil.$^{:1182}$

## Zweite Lösung von Steklow und Markeev

Steklow fand 1899 eine weitere Lösung, bei der er den Ansatz n$_2$ ~ p·q und n$_3$ ~ p·r machte. Der Massenmittelpunkt liegt diesmal auf der zu A gehörenden ersten Hauptachse im Abstand s vom Stützpunkt, womit sich die Euler-Poisson-Gleichungen mit den Bezeichnungen vom Abschnitt #Euler-Poisson-Gleichungen
    ${\begin{aligned}A{\dot {p}}=&(B-C)qr\\B{\dot {q}}=&c_{0}n_{3}+(C-A)rp\\C{\dot {r}}=&-c_{0}n_{2}+(A-B)pq\\{\dot {n}}_{1}=&rn_{2}-qn_{3}\\{\dot {n}}_{2}=&pn_{3}-rn_{1}\\{\dot {n}}_{3}=&qn_{1}-pn_{2}\end{aligned}}$
schreiben. Im Gebiet I, siehe Bild, lautet die Lösung
    ${\begin{aligned}p=&-\mu {\sqrt {\frac {A(2B-A)(A-2C)}{(A-C)^{2}(B-A)}}}\operatorname {cn} (u,k)\\q=&\mu {\sqrt {\frac {A^{2}(A-2C)}{(B-A)(A-C)(B-C)}}}\operatorname {sn} (u,k)\\r=&\mu {\sqrt {\frac {A^{2}(2B-A)}{(A-C)^{2}(B-A)}}}\operatorname {dn} (u,k)\\n_{1}=&1-{\frac {A}{A-C}}\operatorname {cn} ^{2}(u,k)\\n_{2}=&{\sqrt {\frac {A(2B-A)}{(A-C)(B-C)}}}\operatorname {sn} (u,k)\operatorname {cn} (u,k)\\n_{3}=&-{\sqrt {\frac {A(A-2C)}{(A-C)^{2}}}}\operatorname {cn} (u,k)\operatorname {dn} (u,k)\end{aligned}}$
mit
    ${\begin{aligned}u=&\mu {\sqrt {\frac {A(B-C)}{(A-C)(B-A)}}}(t-t_{0})\\\mu =&{\sqrt {\frac {c_{0}}{A}}}\\k=&{\sqrt {\frac {B-A}{B-C}}}\end{aligned}}$
Damit alle Koeffizienten reell sind, muss bei einem realen Starrkörper 2C < A < B < A + C sein. Es zeigt sich, dass der vertikale Drehimpuls null ist und die Gesamtenergie des Kreisels
    $E=c_{0}\left({\frac {A^{2}}{2(A-C)(B-A)}}+1\right)$
wird. Diese Bewegung ist immer instabil.
P. A. Kuz'min (1952)$^{:96}$ oder A. P. Markeev(2005) ergänzten die Lösung im Gebiet II:
    ${\begin{aligned}p=&\mu {\sqrt {\frac {A(2B-A)(A-2C)}{(A-C)^{2}(A-B)}}}\operatorname {sn} (u,k)\\q=&\mu {\sqrt {\frac {A^{2}(A-2C)}{(A-B)(A-C)(B-C)}}}\operatorname {cn} (u,k)\\r=&\mu {\sqrt {\frac {A^{2}(2B-A)}{(A-B)(A-C)(B-C)}}}\operatorname {dn} (u,k)\\n_{1}=&{\frac {A}{A-C}}\operatorname {sn} ^{2}(u,k)-1\\n_{2}=&{\sqrt {\frac {A(2B-A)}{(A-C)(B-C)}}}\operatorname {sn} (u,k)\operatorname {cn} (u,k)\\n_{3}=&-{\sqrt {\frac {A(A-2C)}{(A-C)(B-C)}}}\operatorname {sn} (u,k)\operatorname {dn} (u,k)\end{aligned}}$
mit
    ${\begin{aligned}u=&\mu {\sqrt {\frac {A}{A-B}}}(t-t_{0})\\\mu =&{\sqrt {\frac {c_{0}}{A}}}\\k=&{\sqrt {\frac {A-B}{A-C}}}\end{aligned}}$
Diesmal müssen die Hauptträgheitsmomente die Bedingungen 2C < A, B < A < B + C erfüllen. Der vertikale Drehimpuls ist wieder null und die Gesamtenergie des Kreisels lautet:
    $E=c_{0}\left({\frac {A^{2}}{2(A-B)(A-C)}}-1\right)$
Die Bewegung im Gebiet II ist fast überall stabil.
Da die Trajektorien ausschließlich von den Kreiseleigenschaften festgelegt werden, müssen die Anfangsbedingungen wie bei der Grioli’schen Präzession genau auf den Kreisel abgestimmt sein, denn es ist nur der Zeitparameter t$_0$ frei.