Eulersche Gleichungen (Kreiseltheorie)
# Eulersche Gleichungen (Kreiseltheorie)

Die Euler’schen Kreiselgleichungen oder uneindeutig Euler’schen Gleichungen sind Bewegungsgleichungen für die Rotation eines starren Körpers. Es sind die Komponenten des für den Starrkörper in seinem Hauptachsensystem aufgeschriebenen Drallsatzes und stellen die wichtigste Grundgleichung der Kreiseltheorie dar.
Wird der Körper einem Drehmoment ausgesetzt, entwickeln sich Kreiselwirkungen, die versuchen, die Eigendrehung mit der erzwungenen Drehung in Deckung zu bringen. Die Kreiselwirkungen sind die summierten Drehmomente der Eulerkräfte und Zentrifugalkräfte an allen Massenpunkten des Körpers. Das Moment und die Kreiselwirkungen befinden sich im dynamischen Gleichgewicht, was die Kreiselgleichungen ausdrücken:
    ${\begin{aligned}M_{1}=&\Theta _{1}{\dot {\omega }}_{1}+(\Theta _{3}-\Theta _{2})\omega _{2}\omega _{3}={\dot {L}}_{1}+\left({\frac {1}{\Theta _{2}}}-{\frac {1}{\Theta _{3}}}\right)L_{2}L_{3}\\M_{2}=&\Theta _{2}{\dot {\omega }}_{2}+(\Theta _{1}-\Theta _{3})\omega _{3}\omega _{1}={\dot {L}}_{2}+\left({\frac {1}{\Theta _{3}}}-{\frac {1}{\Theta _{1}}}\right)L_{3}L_{1}\\M_{3}=&\Theta _{3}{\dot {\omega }}_{3}+(\Theta _{2}-\Theta _{1})\omega _{1}\omega _{2}={\dot {L}}_{3}+\left({\frac {1}{\Theta _{1}}}-{\frac {1}{\Theta _{2}}}\right)L_{1}L_{2}\end{aligned}}$
Darin sind jeweils für $k=1,2,3$
    $M_{k}$ die von außen angreifenden Drehmomente,
    $\Theta _{k}$ die Hauptträgheitsmomente,
    $L_{k}=\Theta _{k}\omega _{k}$ die Drehimpulse,
    $\omega _{k}$ die Winkelgeschwindigkeiten und
    ${\dot {\omega }}_{k}$ die Winkelbeschleunigungen
im Hauptachsensystem. Gelegentlich wird auch die dazu gehörige Vektorgleichung
    ${\vec {M}}=\mathbf {\Theta } \cdot {\dot {\vec {\omega }}}+{\vec {\omega }}\times \mathbf {\Theta } \cdot {\vec {\omega }}={\frac {\mathrm {d} _{r}}{\mathrm {d} t}}{\vec {L}}+{\vec {L}}\cdot \mathbf {\Theta } ^{-1}\times {\vec {L}}$
mit dem Trägheitstensor $\mathbf {\Theta }$ als Euler’sche Kreiselgleichung angegeben. Hier bildet „·“ die Vektortransformation, „×“ das Kreuzprodukt und ${\tfrac {\mathrm {d} _{r}}{\mathrm {d} t}}$ die relative Zeitableitung im Hauptachsensystem.
Die Drehmomente, Hauptträgheitsmomente und Drehimpulse werden mit einem Bezugspunkt berechnet, für den sich der Massenmittelpunkt oder ein unbeschleunigter, in einem Inertialsystem ruhender Stützpunkt eignen, siehe Drallsatz am starren Körper.
Die ersten Summanden auf den rechten Seiten, bestehend aus den Winkelbeschleunigungen und Drehimpulsänderungen, resultieren aus den Kreiselwirkungen der Euler-Kräfte und die anderen, in den Winkelgeschwindigkeiten und Drehimpulsen quadratischen Terme berücksichtigen die Kreiselwirkungen der Zentrifugalkräfte. Wenn die Bewegung bekannt ist, dann können aus diesen Gleichungen die Momente berechnet werden, die im Bezugspunkt eingeleitet werden müssen, damit der Körper die vorgegebene Bewegung ausführt.
Die Kreiselgleichungen wurden von Leonhard Euler 1750 aufgestellt und später zum Drallsatz weiterentwickelt.

## Spezialfälle
### Euler-Poisson-Gleichungen

Die Euler-Poisson-Gleichungen sind die spezifischen Kreiselgleichungen für den schweren Kreisel, bei dem das äußere Moment von der Schwerkraft herrührt. Die klassische Kreiseltheorie ist fast ausschließlich dem schweren Kreisel mit Stützpunkt gewidmet und es wurde viel Aufwand in das Auffinden exakter Lösungen gesteckt. Eine Auflistung einiger dieser Lösungen findet sich im Hauptartikel.

### Kugelkreisel

Ein Kugelkreisel ist ein Kreisel mit drei identischen Hauptträgheitsmomenten Θ, sodass sich die Kreiselgleichungen dann auf
    ${\vec {M}}=\Theta {\dot {\vec {\omega }}}$
reduzieren. Die Winkelbeschleunigung ist beim Kugelkreisel also parallel zum angreifenden Moment. Beim Kugelkreisel sind die Fliehkräfte im Körper immer im mechanischen Gleichgewicht. Ein Vergleich mit den Bewegungsgleichungen bei einer Translationsbewegung zeigt, dass der Kugelkreisel das genaue Analogon des Massenpunkts bei Rotationsbewegungen ist.

### Ebene Bewegungen

Bei einer ebenen Bewegung um eine Hauptträgheitsachse, beispielsweise die 3-Achse, entfallen Drehungen und Momente um die 1- und 2-Achsen, und die Gleichungen reduzieren sich auf
    $M_{3}=\Theta _{3}{\ddot {\varphi }}\,,$
wobei φ der Drehwinkel um die 3-Achse ist.

## Lösungen der Kreiselgleichungen bei ebenen Bewegungen

Im ebenen Fall sind die Kreiselgleichungen oftmals analytisch lösbar, wofür die beiden folgenden Fälle Beispiele sind.

### Anstoß einer Billardkugel

Parallel zur Tischplatte soll eine Billardkugel mit Radius r, Masse m und Massenträgheitsmoment Θ so angestoßen werden, dass sie nicht über den Tisch rutscht, siehe Abb. 1. Es stellt sich die Frage, in welcher Höhe h über der Platte die Kraft F eingeleitet werden muss, damit für das schlupflose Rollen keine Reibkraft am Tisch notwendig ist.
Die exzentrisch an der Kugel angreifende horizontale Kraft entwickelt ein Moment M = -(h-r)F um den Massenmittelpunkt, das die Kugel gemäß der Kreiselgleichung
    $M=-(h-r)F=\Theta {\ddot {\varphi }}$
in Drehung versetzt. Das Moment ist negativ, weil es entgegen der Zählrichtung des Drehwinkels φ wirkt. Außerdem beschleunigt die Kraft die Kugel gemäß dem Gesetz „Kraft gleich Masse mal Beschleunigung“:
    $F=m{\ddot {x}}\,.$
Die Beschleunigung ${\ddot {x}}$ ist parallel zum Tisch in Richtung der Kraft. Die Bedingung für schlupfloses Rollen
    ${\ddot {x}}=-r{\ddot {\varphi }}$
schließt das Gleichungssystem für die drei Unbekannten h, φ und x ab. Damit berechnet sich
    $-(h-r)F=\Theta {\ddot {\varphi }}=-{\frac {\Theta }{r}}{\ddot {x}}=-{\frac {\Theta }{mr}}F\quad \rightarrow \quad h=r+{\frac {\Theta }{mr}}\,.$
Die Kraft muss demnach im Stoßmittelpunkt angreifen, damit die Kugel schlupflos rollt. Bei einer massiven homogenen Kugel ist das Massenträgheitsmoment Θ = $^2$⁄$_5$mr² und somit
    $h={\frac {7}{5}}r={\frac {7}{10}}d\,,$
wobei d = 2r der Durchmesser der Kugel ist.

### Ein eine schiefe Ebene hinabrollendes Rad

In einer ebenen Bewegung rolle ein Rad mit Radius r, Masse m und Massenträgheitsmoment Θ eine mit dem Winkel α geneigte Ebene unter Einfluss einer Schwerebeschleunigung g hinab, siehe Abb. 2. Weil sich das Rad dabei auch translatorisch bewegt, geht auch seine Masse in die Beschleunigung ein. Die Beschleunigung wächst jedoch, wenn das Massenträgheitsmoment abnimmt.
Aufgrund des schlupflosen Abrollens entsteht am Aufstandspunkt des Rades eine Reibkraft R, die das Rad in Drehung versetzt, denn es entspricht einem Moment M=-r R um den Massenmittelpunkt. Das Moment ist negativ, weil es entgegen der Zählrichtung des Drehwinkels φ arbeitet. Damit lautet die Kreiselgleichung im ebenen Fall:
    $-rR=\Theta {\ddot {\varphi }}\,.$
Die auf das Rad hangabwärts wirkende Komponente F der Gewichtskraft mg hat die Größe F=mgsin(α). Ihr entgegen steht die Reibkraft, sodass nach dem zweiten newtonschen Gesetz gilt:
    $F-R=m{\ddot {x}}\quad \rightarrow \quad R=mg\sin(\alpha )-m{\ddot {x}}\,,$
worin ${\ddot {x}}$ die hangabwärts zählende Beschleunigung des Rades und sin die Sinusfunktion ist. Die Bedingung für schlupfloses Rollen ${\ddot {x}}=-r{\ddot {\varphi }}$ schließt das Gleichungssystem für die drei Unbekannten R, φ und x ab und es ergibt sich
    ${\ddot {x}}={\frac {mr^{2}}{\Theta +mr^{2}}}g\sin(\alpha )$
Ein hangabwärts, reibungsfrei rutschender Klotz erfährt die Beschleunigung ${\ddot {x}}=g\sin(\alpha )$, die größer ist als die des Rades, denn beim Rad wird ein Teil der potentiellen Energie in Rotationsenergie umgesetzt, die dann für die Translation fehlt.

## Beispiel: Fliehkraftpendel

Ein bemerkenswerter Anwendungsfall der Euler’schen Kreiselgleichungen ist die "Perle auf rotierendem Kreisring im Schwerefeld" das konische oder Fliehkraftpendel, denn es hat überraschende Eigenschaften. In Worten der Kreiseltheorie handelt es sich um einen schweren Kreisel, der in einem Punkt festgehalten wird und bei dem die Eigendrehung um die Pendelachse durch äußere Momente verhindert und die Präzessionsgeschwindigkeit konstant gehalten wird. Komplexere Eigenschaften treten bei beliebiger Lage des Massenmittelpunkts und anfänglicher Ausrichtung der Hauptachsen auf. Vereinfachungen ergeben sich mit Massenmittelpunkt auf der 1-Achse und permanent horizontaler 3-Achse des Hauptachsen­systems.

### Kinematik

Die in der Kreiseltheorie üblichen Definitionen der Winkel und Winkelgeschwindigkeiten werden übernommen. Insbesondere wird die z-Komponente der 2-Achse
    $n_{2}=\sin(\vartheta )\cos(\varphi )=\cos(\varphi )$
bei horizontaler 3-Achse (ϑ ≡ 90°) benötigt, die identisch mit der 2-Komponente des Einheitsvektors ist, der vertikal nach oben weist. Die Funktionen sin und cos bilden den Sinus und Kosinus.
Die Winkelgeschwindigkeiten lauten unter den hiesigen Bedingungen:
    ${\begin{aligned}\mu :=&{\dot {\psi }}=\omega _{1}\sin(\varphi )+\omega _{2}\cos(\varphi )={\text{const.}}\\{\dot {\vartheta }}=&\omega _{1}\cos(\varphi )-\omega _{2}\sin(\varphi )\equiv 0\\\omega _{1}=&\mu \sin(\varphi ),\;\omega _{2}=\mu \cos(\varphi ),\;\omega _{3}={\dot {\varphi }}\end{aligned}}$
Der Winkel ϑ ist unveränderlich, wenn er anfangs keine Zeitableitung hat und dauerhaft keine Beschleunigung erfährt.

### Kinetik

Die 1- oder Pendelachse trägt im Abstand s > 0 vom Fixpunkt den Massenmittelpunkt, dessen Gewichtskraft mg somit ein Schweremoment auf das Pendel ausübt. Zusätzlich zu diesem Schweremoment wirken äußere Momente M$_{1,2}$ in 1- bzw. 2-Richtung des Hauptachsensystems, um Winkelbeschleunigungen von ϑ und μ zu unterdrücken. Die Kreiselgleichungen in 1- und 2-Richtung dienen nur der Bestimmung dieser Momente, die hier nicht interessieren. Bemerkenswert ist jedoch, dass diese Drehmomente eine zeitliche Änderung des vertikalen Drehimpulses verursachen, der bei schweren Kreiseln immer konstant ist.
Um die horizontale 3-Richtung kann das Pendel frei drehen und das Schweremoment in dieser Richtung ist -mgsn$_2$ = -mgs cos(φ), siehe Kinematik oben und Euler-Poisson-Gleichungen. Damit liefert die dritte der Kreiselgleichungen
    $M_{3}=-mgs\cos(\varphi )=C{\dot {\omega }}_{3}+(B-A)\omega _{1}\omega _{2}$
die Bewegungsgleichung
Darin sind A, B und C die Hauptträgheitsmomente um die 1-, 2- bzw. 3-Achse.

### Gleichgewichtslagen

Im Gleichgewicht verschwindet in (*) die Winkelbeschleunigung:
    $C{\ddot {\varphi }}=[(A-B)\mu ^{2}\sin(\varphi )-mgs]\cos(\varphi )=0$.
Anders als beim sich selbst überlassenen Lagrange-Kreisel ist eine Präzession mit waagerechter Achse (φ = 0) nicht möglich. Mangels Eigendrehung gibt es hier keine Kreiselwirkung, die das Schweremoment ausgleichen könnte. Das Gleichgewicht mit lotrechter 1-Achse (φ = ±90°) ist allen Kreiseln zugänglich – auch dem Kugelkreisel, der nur dort Gleichgewicht findet. Bei A ≠ B ergeben sich die Gleichgewichtslagen aus
    $\left[\sin(\varphi )-{\frac {mgs}{(A-B)\mu ^{2}}}\right]\cos(\varphi )=0$
Wenn A - B dasselbe Vorzeichen wie sin(φ) hat, kann der Term in der eckigen Klammer null werden. Dazu muss
    $|\sin(\varphi )|\leq 1\quad \Rightarrow \quad |\mu |\geq {\sqrt {\frac {mgs}{|A-B|}}}=:\mu _{c}$
erfüllt sein. Das Pendel muss eine kritische Drehgeschwindigkeit μ$_c$ überschreiten, damit diese Gleichgewichtslage abseits der Senkrechten existiert. Beim abgeplatteten Kreisel weist die 1-Achse im Gleichgewicht nach oben (ist A > B und sin(φ) > 0) und beim gestreckten nach unten (A < B und sin(φ) < 0).
Analog zum symmetrischen Kreisel wird das Pendel hier und im folgenden abgeplattet genannt, wenn A > B, und gestreckt, wenn B > A.

### Energiebetrachtung

Multiplikation der Bewegungsgleichung (*) mit ${\dot {\varphi }}$ ermöglicht eine Zeitintegration:
    $E={\frac {C}{2}}{\dot {\varphi }}^{2}+{\frac {B-A}{2}}\mu ^{2}\sin ^{2}(\varphi )+mgs\sin(\varphi )$
Die Integrationskonstante E ist die Energie des Pendels im rotierenden System:
* der erste Summand trägt die Rotationsenergie um die 3-Achse bei,
* der zweite Summand ist das Zentrifugalpotential der Präzession und
* der dritte Summand steht für die Lageenergie im Schwerefeld der Erde.
Die Energie und ihre Ableitungen (·)' nach φ ergeben sich mit der Abkürzung $z:={\tfrac {mgs}{(A-B)\mu ^{2}}}$ zu
    ${\begin{aligned}E=&{\frac {C}{2}}{\dot {\varphi }}^{2}+(B-A)\mu ^{2}{\bigg{\frac {\sin ^{2}(\varphi )}{2}}-z\sin(\varphi ){\bigg )}\\E'=&(B-A)\mu ^{2}[\sin(\varphi )-z]\cos(\varphi )\\E''=&(B-A)\mu ^{2}[1+z\sin(\varphi )-2\sin ^{2}(\varphi )]\end{aligned}}$
Beim gestreckten Kreisel ist B > A und z < 0 und beim abgeplattenten Kreisel ist B < A sowie z > 0. Stabilität liegt in einem Energieminimum vor, wo E' = 0 und E" > 0 ist. Die Energie ist in Gleichgewichtslagen stationär: bei sin(φ) = z oder cos(φ) = 0, wo sin(φ) = -1 oder sin(φ) = +1 ist:

#### Gleichgewichtslage mit sin(φ) = z

    Hier ist E" = (A-B)μ²(z²-1). Diese Bewegung ist beim gestreckten Kreisel stabil, wenn z$^2$ < 1, also immer, und beim abgeplatteten, wenn z$^2$ > 1, also nie.

#### Hängendes Pendel sin(φ) = -1

    Hier ist E" = (A-B)μ²(1+z). Der untere Totpunkt ist beim abgeplatteten Kreisel stabil, wenn z > -1, also immer, und beim gestreckten, wenn z < -1, also |μ| < μ$_c$. Bei zunehmendem μ durchläuft das gestreckte Pendel bei μ = μ$_c$ eine subkritische Pitchfork-Bifurkation, wo die stabile Lage am unteren Totpunkt instabil wird und zwei stabile Gleichgewichtslagen mit sin(φ) = z entstehen. Im Gegensatz dazu ist der lotrecht hängende Lagrange-Kreisel immer stabil, siehe dort.

#### Aufrechtes Pendel sin(φ) = +1

    Hier ist E" = (A-B)μ²(1-z). Der obere Totpunkt ist beim gestreckten Kreisel stabil, wenn z > 1, also nie, und beim abgeplatteten, wenn z < 1 oder |μ| > μ$_c$. Bei zunehmendem μ durchläuft das abgeplattete Pendel bei μ = μ$_c$ eine superkritische Pitchfork-Bifurkation, wo die instabile Lage am oberen Totpunkt stabil wird und zwei neue instabile Gleichgewichtslagen mit sin(φ) = z entstehen.

### Schwingungen

Die Gleichung (*) kann um Gleichgewichtslagenφ$_0$, wo z = sin(φ$_0$) ist, linearisiert werden. Dazu wird φ = φ$_0$+δ mit konstantem φ$_0$ und kleiner Abweichung δ angenommen. Die Drehbeschleunigung
    $C{\ddot {\varphi }}=C{\ddot {\delta }}=(A-B)\mu ^{2}\sin(\varphi _{0}+\delta )\cos(\varphi _{0}+\delta )-mgs\cos(\varphi _{0}+\delta )$
entwickelt sich mit den Additionstheoremen und sin(δ) ≈ δ, cos(δ) ≈ 1 zur Schwingungsgleichung
    ${\ddot {\delta }}+{\frac {B-A}{C}}\mu ^{2}\cos ^{2}(\varphi _{0})\delta =0$ mit $\sin(\varphi _{0})={\frac {mgs}{(A-B)\mu ^{2}}}$
Beim gestreckten Kreisel ist B > A und der "Rückstellkoeffizient" vor δ immer positiv, weswegen der gestreckte Kreisel kleine Schwingungen um Gleichgewichtslagen ausführen kann. Beim abgeplattenten Kreisel ist wegen A > B der Rückstellkoeffizient negativ. Daher kann der abgeplattete Kreisel keine Schwingungen um Gleichgewichtslagen abseits der Senkrechten ausführen, was im Einklang mit den Ergebnissen aus der Energiebetrachtung ist.