Euler-Poisson-Gleichungen
Die Euler-Poisson-Gleichungen nach Leonhard Euler und Siméon Denis Poisson sind die in der Kreiseltheorie benutzten Bewegungsgleichungen für den schweren Kreisel mit Stützpunkt. Sie stellen für diesen Kreisel die Komponenten des Drallsatzes und der Zeitableitung der Lotrichtung oder Gewichtskraft im Hauptachsensystem dar.
Die klassische Kreiseltheorie ist fast ausschließlich dem schweren Kreisel mit Stützpunkt gewidmet und es wurde und wird[1] viel Aufwand in das Auffinden exakter Lösungen gesteckt. Im Zeitalter leistungsfähiger Rechenmaschinen haben diese Lösungen nicht mehr die früher berechtigte, zentrale Bedeutung. Heute bereitet es keine Schwierigkeiten, die Euler-Poisson-Gleichungen mit beliebigen Anfangsbedingungen durch Numerische Simulation mit jeder gewünschten Genauigkeit zu berechnen.[2]
Wilhelm Hess konnte 1890 mit Hilfe der Integrale der Bewegung die Richtungskosinus eliminieren.[3] Mit Hilfe dieser Formulierung untersucht die Forschung Anfang des 21. Jahrhunderts die Topologie der Energieflächen des schweren Kreisels mit Stützpunkt.[4]
Geschichte
[Bearbeiten | Quelltext bearbeiten]Leonhard Euler stellte 1750 die Kreiselgleichungen auf und konnte 1758 bereits eine Lösung – den Euler-Kreisel – angeben. Joseph-Louis Lagrange leistete 1788 einen wichtigen Beitrag durch Lösung der Gleichungen für den symmetrischen schweren Kreisel mit Fixpunkt. Carl Gustav Jacob Jacobi veröffentlichte 1829 die Theorie der Jacobi’schen elliptischen Funktionen und der Theta-Funktionen, mit denen sich die Euler-Poisson-Gleichungen lösen lassen. Sofia Kowalewskaja entdeckte 1888 den letzten durch Theta-Funktionen lösbaren Fall, den schweren, symmetrischen, inhomogenen Kowalewskaja-Kreisel, und reformulierte das Problem mit analytischen Funktionen einer komplexen Zahl.
A. M. Ljapunow bewies 1894, dass die drei Fälle von Euler, Lagrange und Kowalewskaja die einzigen sind, in denen die allgemeine Lösung bei beliebigen Anfangsbedingungen eine eindeutige Funktion der Zeit ist. In anderen Fällen sind die Lösungen für bestimmte Anfangsbedingungen mehrdeutige Funktionen der Zeit. Édouard Husson zeigte 1905[5], dass der Euler-, Lagrange- und Kowalewskaja-Kreisel die einzigen mit algebraischen ersten Integralen lösbaren Fälle der Euler-Poisson-Gleichungen sind. Darüber hinaus können nur spezielle Bewegungen integrabel sein. Ein viertes algebraisches Integral existiert nur in den Fällen, wo die Lösung eine eindeutige Funktion der Zeit ist.[6]
Formulierung
[Bearbeiten | Quelltext bearbeiten]Vektorgleichungen
[Bearbeiten | Quelltext bearbeiten]Das Schweremoment ergibt sich aus dem Kreuzprodukt × des Hebelarms vom Stützpunkt zum Massenmittelpunkt mit der Gewichtskraft des Kreisels:
Darin ist m die Masse, g die Schwerebeschleunigung und êz der lotrecht nach oben weisende Einheitsvektor. Das Schweremoment wird in den Drallsatz eingesetzt[7]:
Darin ist Θ der Trägheitstensor, Θ –1 seine inverse, die Winkelgeschwindigkeit und der Drehimpuls des Kreisels bezüglich des Stützpunkts und bildet die relative Zeitableitung im Hauptachsensystem. Die Eulerʹschen Kreiselgleichungen sind die Komponenten der Vektorgleichung in ebendiesem System im Fall eines beliebigen äußeren Moments[8].
Die von Poisson entdeckten kinematischen Gleichungen[9] formulieren mathematisch die Konstanz der Lotrichtung im Hauptachsensystem und schreiben sich in Vektorform
wo wieder die relative Zeitableitung im Hauptachsensystem bildet. Die Komponenten dieser Gleichung werden zu Ehren ihres Entdeckers Poisson-Gleichungen genannt[10].
Euler-Poisson-Gleichungen
[Bearbeiten | Quelltext bearbeiten]Das autonome gewöhnliche Differentialgleichungssystem aus Euler- und Poisson-Gleichungen wird Euler-Poisson-Gleichungen genannt[10][11] und schreibt sich im Hauptachsensystem
Der Überpunkt bildet die Zeitableitung, mg ist die Gewichtskraft und für k = 1,2,3 ist jeweils
- ωk die Komponente der Winkelgeschwindigkeit,
- Lk die Komponente des Drehimpulses,
- sk die Komponente des Massenmittelpunkts,
- nk die Komponente des Einheitsvektors êz und
- Θk ein Hauptträgheitsmoment des Kreisels
im Hauptachsensystem. Häufig werden
- die Hauptträgheitsmomente Θ1,2,3 mit A, B bzw. C,
- die Winkelgeschwindigkeiten ω1,2,3 mit p, q bzw. r und
- die Richtungskosinus n1,2,3 mit γ1,2,3 oder γ, γ', γ", gelegentlich auch mit umgekehrtem Vorzeichen
bezeichnet. Manchmal werden die Kreiselgleichungen durch die Gewichtskraft dividiert oder der Faktor mg den Richtungskosinus zugeschlagen, sodass in den obigen Formeln dieser Faktor nicht mehr auftritt.
Berechnung des Präzessionswinkels
[Bearbeiten | Quelltext bearbeiten]Weil die Lotrichtung vom Präzessionswinkel ψ der Drehung um die Lotrichtung unabhängig ist, kann aus den Euler-Poisson-Gleichungen der Präzessionswinkel nicht berechnet werden. Er kann jedoch mit der Differentialgleichung
zeitgleich oder nachträglich ermittelt werden.[11]
Integrale der Bewegung
[Bearbeiten | Quelltext bearbeiten]Bei jedem schweren Kreisel mit Stützpunkt ist die Norm des Richtungsvektors der Lotrichtung, der Drehimpuls in Lotrichtung und die Gesamtenergie E konstant:
Die Konstanten sind Integrale der Bewegung und werden in der Kreiseltheorie kurz Integrale genannt, die ersten beiden auch Casimir-Invarianten. Das zweite Integral Lz heißt nach dem Drall- oder Flächensatz auch Drall- bzw. Flächenintegral. Es ist Konstant weil das Schweremoment keine Komponente in Lotrichtung hat. Die Gesamtenergie E wird in der analytischen Mechanik auch als Hamilton-Funktion bezeichnet. Sie ist konstant weil das Schwerefeld der Erde konservativ ist und somit die Kreiselbewegung den Energieerhaltungssatz befolgt. Die Unveränderlichkeit der Integrale lässt sich auch durch Zeitableitung und Einsetzen der Euler-Poisson-Gleichungen nachweisen.
Hess’sche Gleichungen
[Bearbeiten | Quelltext bearbeiten]Wilhelm Hess hat 1890 in seinem Aufsatz[3], in dem er das loxodromische Pendel einführte, alternative Formulierungen veröffentlicht. Es gelang ihm die Richtungskosinus n1,2,3 mit den Integralen und dem Drehimpuls auszudrücken:
Darin ist T die Rotationsenergie, die gleich dem ersten Summanden im obigen Integral E ist, und . Damit entsteht Hess’ Bewegungsgleichung (6):
die nur noch vom Drehimpuls oder der Winkelgeschwindigkeit abhängt. Es ist erwiesen, dass √f ein integrierender Faktor für diese Gleichung ist, weswegen schon ein weiteres Integral genügt, um die Bewegungsgleichungen zu lösen. Die Funktion f ist für die Topologie der Energieflächen bedeutsam.[4]
Hess konnte auch schon
angeben. In seinen Gleichungen (8) ersetzte er die Winkelgeschwindigkeiten durch das Drehimpulsbetragsquadrat , die Projektion des Drehimpulses auf die Schwereachse und die Rotationsenergie T:
Darin ist . Ähnliche Gleichungen fand P. A. Schiff 1904[12]. Die Gleichungen sind äquivalent zu den Euler-Poisson-Gleichungen, vorausgesetzt es werden keine Nebenbedingungen an die Variablen ν, ρ und T gestellt[13].
Zusammenhang mit der Riccatischen Differentialgleichung
[Bearbeiten | Quelltext bearbeiten]Bei gefundenen oder gegebenen Winkelgeschwindigkeiten ω1,2,3 stellt sich die Frage nach den entsprechenden Richtungskosinus n1,2,3. Nach Gaston Darboux empfiehlt sich die Lösung mittels der Riccatischen Differentialgleichung. Dazu wird mit der imaginären Einheit i und unbekannten Funktionen x und y
gesetzt. Die Funktionen x und y sind infolge dessen Integrale der Riccatischen Differentialgleichung
Diese lässt sich durch Quadraturen integrieren, sobald ein partikuläres Integral von ihr bekannt ist.[14]
Lösungen der Euler-Poisson-Gleichungen
[Bearbeiten | Quelltext bearbeiten]Für die technische Anwendung gibt es bedeutsame Spezialfälle, bei denen sich die Euler-Poisson-Gleichungen soweit vereinfachen, dass sie integrabel sind. In diesen Fällen weisen die Trajektorien des Kreisels einen zumindest quasi-periodischen Verlauf auf, können die verschiedenen Bewegungsmodi klassifiziert und die Zeitfunktionen der Variablen sowie ihre geometrische Bedeutung angegeben werden. Insbesondere beim Kowalewskaja-Kreisel und im Gorjatschew-Tschaplygin-Kreisel sind die analytischen Lösungen so kompliziert, dass die Herausarbeitung der vorgenannten typischen Eigenschaften der Bewegung äußerst aufwändig ist. Hier helfen topologische Analyse (Bifurkationsdiagramm), Stabilitätsanalyse, Phasenraum-Diagramme und Computeranimationen dabei, Einblicke in die Vorgänge im Kreisel zu erhalten und deren typischen Eigenschaften heraus zu arbeiten. Die so erzielten Ergebnisse können praktische Anwendungen motivieren[15].
Die folgende Tabelle enthält eine Auswahl räumlicher Bewegungen von Kreiseln, in denen bis Anfang des 21. Jahrhunderts exakte Lösungen der Euler-Poisson-Gleichungen gelungen sind[16].
Entdecker | Hauptträgheits- momente |
Lage des Schwerpunkts |
Anfangsbedingungen (t = 0) |
---|---|---|---|
Leonhard Euler siehe Euler-Kreisel |
beliebig | s1=s2=s3=0 | beliebig |
Joseph-Louis Lagrange siehe Lagrange-Kreisel |
A=B | s1=s2=0, s3≠0 | beliebig |
Sofia Kowalewskaja siehe Kowalewskaja-Kreisel |
A=B=2C | beliebig | |
Wilhelm Hess siehe Hesssches Pendel |
beliebig | s2=0 |
|
Gorjatschew und Chaplygin |
A=B=4C | Lz=0 | |
Merzalow | A=B=4C | r = 0 | |
Bobylev und Steklow siehe Bobylew-Steklow-Lösung |
2A=C | s1=s2=0, s3≠0 | |
Otto Staude siehe Staude-Drehung |
beliebig | beliebig | |
Giuseppe Grioli siehe Griolische Präzession |
beliebig | s2=0 |
Bis auf einen Freiheitsgrad eindeutig festgelegt |
In der Tabelle sind A, B, C = θ1,2,3 die Hauptträgheitsmomente, p, q, r = ω1,2,3 die Winkelgeschwindigkeiten und s1,2,3 die (konstanten) Koordinaten des Massenmittelpunkts im Hauptachsensystem.
Einzelnachweise
[Bearbeiten | Quelltext bearbeiten]- ↑ S. V. Ershakov: New exact solution of Euler’s equations (rigid body dynamics) in the case of rotation over the fixed point. In: Archive of Applied Mechanics. Band 84, Nr. 3. Springer-Verlag, 2014, ISSN 0939-1533, S. 385–389, doi:10.1007/s00419-013-0806-x.
- ↑ Magnus (1971), S. 109.
- ↑ a b Wilhelm Hess: Ueber die Euler’schen Bewegungsgleichungen und über eine neue particuläre Lösung des Problems der Bewegung eines starren Körpers um einen festen Punkt. In: Mathematische Annalen. Vol. 37, 1890, S. 153–181 (digizeitschriften.de [abgerufen am 2. Mai 2018]).
- ↑ a b I. G. Gashenenko, P. H. Richter: Enveloping Surfaces And Admissible Velocities Of Heavy Rigid Bodies. In: World Scientific Publishing Company (Hrsg.): International Journal of Bifurcation and Chaos. Band 14, Nr. 8, 2004, ISSN 0218-1274, S. 2525–2553, doi:10.1142/S021812740401103X (iamm.su [PDF; abgerufen am 2. Juni 2019] siehe S. 2537).
- ↑ Édouard Husson: Recherche des intégrales algébriques dans le mouvement d’un solide pesant autour d’un point fixe. In: Annales de la faculté des sciences de Toulouse 2e série. 1906, S. 73–152, doi:10.5802/afst.232 (französisch, numdam.org [PDF; abgerufen am 7. März 2018] Auf Seite 74 wird ein erster Beweisversuch von Roger Liouville 1897 als fehlerhaft aufgedeckt.).
- ↑ Leimanis (1965), S. 53 ff.
- ↑ In der Tensoralgebra kann auf Klammerungen verzichtet werden:
- ↑ Magnus (1970), S. 106.
- ↑ Siméon Denis Poisson: Traité de Méchanique. 3. Auflage. 1 bis 6. J. G. Garnier, Brüssel 1838 (französisch, archive.org [abgerufen am 3. November 2019]).
- ↑ a b Leimanis (1965), S. 7.
- ↑ a b Peter H. Richter, Holger R. Dullin, Andreas Wittek: Kovalevskaya Top. Hrsg.: Institut für den wissenschaftlichen Film (IWF). 1997, ISSN 0073-8433, S. 41 (englisch, researchgate.net [abgerufen am 28. März 2018] Dort hat φ die Bedeutung von ψ hier, siehe auch Kreiseltheorie#Bezugssysteme und Euler-Winkel).
- ↑ P. A. Schiff (П. А. Шиффъ): Über die Bewegungsgleichungen eines schweren starren Körpers mit einem festen Punkt. In: Матем. сб. Band 24, Nr. 2, 1904, S. 169–177 (russisch, mathnet.ru [abgerufen am 10. Juni 2019] Originaltitel: Объ уравненiяхъ движенiя тяжелаго твердаго тѣла, имѣющаго неподвижную точку.).
- ↑ Leimanis (1965), S. 104.
- ↑ Felix Klein, Conr. Müller: Encyklopädie der Mathematischen Wissenschaften mit Einschluss ihrer Anwendungen. Mechanik. Hrsg.: Akademien der Wissenschaften zu Göttingen, Leipzig, München und Wien. 4. Band, 1. Teilband. B. G. Teubner, Leipzig 1908, ISBN 3-663-16021-1, S. 565, doi:10.1007/978-3-663-16021-2 (eingeschränkte Vorschau in der Google-Buchsuche [abgerufen am 7. März 2020] siehe auch wikisource}).
- ↑ A. V. Borisov, I. S. Mamaev: Euler-Poisson Equations and Integrable Cases. 2001, doi:10.1070/RD2001v006n03ABEH000176, arxiv:nlin/0502030 (englisch, Enthält Lösungen der Euler-Poisson-Gleichungen, deren ausführliche Beschreibung und weiter führende Literaturangaben.).
- ↑ Magnus (1971), S. 108.
Literatur
[Bearbeiten | Quelltext bearbeiten]- K. Magnus: Kreisel: Theorie und Anwendungen. Springer, 1971, ISBN 3-642-52163-0, S. 109. (eingeschränkte Vorschau in der Google-Buchsuche [abgerufen am 7. Februar 2019]).
- Eugene Leimanis: The General Problem of the Motion of Coupled Rigid Bodies about a Fixed Point. Springer Verlag, Berlin, Heidelberg 1965, ISBN 3-642-88414-8, S. 53 f., doi:10.1007/978-3-642-88412-2 (englisch, eingeschränkte Vorschau in der Google-Buchsuche [abgerufen am 21. März 2018]).