Sphärisches Pendel
Ein sphärisches Pendel, auch Kugelpendel oder räumliches Pendel, ist ein Pendel, dessen Aufhängung Ausschläge in unterschiedliche Richtungen zulässt. Im Unterschied zum (ebenen) Kreispendel, bei dem die Bewegung der Pendelmasse auf einen vertikalen Kreis beschränkt ist, bewegt sich beim (räumlichen) Kugelpendel die Pendelmasse auf einer Kugelfläche.
Ein Spezialfall des Kugelpendels ist das konische Pendel, auch Kegelpendel, Kreispendel, Rundlaufpendel oder Zentrifugalpendel, bei dem sich die Pendelmasse auf einer horizontalen Kreisbahn bewegt und der Faden deshalb eine Kegelfläche beschreibt.[1]
In der theoretischen Behandlung des sphärischen Pendels wird häufig vereinfachend die Aufhängung als masselos und der Pendelkörper als punktförmig angenommen sowie der Einfluss der Reibung vernachlässigt. Neben der Energieerhaltung ist beim sphärischen Pendel auch die Drehimpulserhaltung von Bedeutung. In der Projektion auf eine horizontale Ebene überstreicht der Pendelfaden daher in gleichen Zeiten gleiche Flächen (siehe Flächensatz).
Eine Anwendung des sphärischen Pendels ist das Foucaultsche Pendel, mit dessen Hilfe ohne Bezug auf Beobachtungen am Himmel die Erdrotation anschaulich nachgewiesen werden kann.
Behandlung nach Lagrange
Allgemeiner Fall
Da sich die Pendelmasse des Kugelpendels auf einer Kugelfläche bewegt, lässt sich seine Bewegung am besten in Kugelkoordinaten beschreiben:
Der Aufhängepunkt ist der Ursprung und die z-Achse weist zur stabilen unteren Ruhelage. Dann ist
- die Länge des Pendels, die sich wegen der starren Verbindung zwischen Aufhängungspunkt und Pendelkörper nicht ändern kann
- der Polarwinkel die Auslenkung aus der unteren Gleichgewichtslage
- der Azimutwinkel die Rotation um die senkrechte -Achse.
Da die Länge konstant gehalten wird, sind die beiden Winkel die einzigen freien Variablen, also die generalisierten Koordinaten für dieses System. Es ist nun die Lagrange-Funktion
zu bilden, wobei die kinetische Energie und die potentielle Energie in Abhängigkeit von den beiden generalisierten Koordinaten und ihren Zeitableitungen bezeichnen.
Die potentielle Energie des Pendels bezüglich des Aufhängepunktes beträgt
und hat ihr Minimum bei . Die kinetische Energie beträgt
- .
Die Bewegungsgleichungen ergeben sich dann aus den Lagrangegleichungen 2. Art:
Die Lagrange-Gleichungen ergeben (nach Kürzen von ):
- .
- .
Die zweite Lagrange-Gleichung führt sofort auf
- .
Diese Gleichungen bilden ein System von zwei gekoppelten Differentialgleichungen 2. Ordnung, von denen die zweite allerdings sofort einmal integriert werden kann, wie man an der darüberstehenden Lagrange-Gleichung sieht, aus der sie hervorgegangen ist.
Nach dieser zweiten Lagrange-Gleichung ist der zu gehörige konjugierte Impuls nämlich konstant. Es ist die -Komponente des Drehimpulses
( kommt in Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L} nicht vor und ist daher eine zyklische Variable. Dies ist ein Beispiel für das Noether-Theorem.)
Damit lässt sich Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot\phi} aus der Differentialgleichung für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta} eliminieren:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \ddot \theta = \frac{L_z^2}{m^2 R^4} \frac{\cos \theta}{\sin^3 \theta} - \frac gR \sin \theta} .
Diese Bewegungsgleichung für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta} ist im Allgemeinen nicht elementar lösbar, und es können sich komplexe Bewegungen ergeben. Ein einfach lösbarer Fall ist das konische Pendel (s. u.).
Aussagen zu allgemeinen Eigenschaften der Bewegung lassen sich gewinnen, wenn zusätzlich die Konstanz der Gesamtenergie Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle E = T + V} berücksichtigt wird, die sich daraus ergibt, dass keine explizite Abhängigkeit von der Zeit vorliegt. Daraus folgt[2]:
- Die Bewegung ist auf einen Bereich Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta_\text{min} \le \theta \le \theta_\text{max}} eingeschränkt, findet also zwischen zwei Breitenkreisen statt.
- Die Auf- und Abbewegung zwischen den Breitenkreisen ist periodisch (aber nicht harmonisch).
- Die azimutale Winkelgeschwindigkeit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot \phi} ist entweder konstant Null (wenn Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L_z=0} ) oder hat das feste Vorzeichen von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L_z \ne 0} . Der Drehsinn der Pendelbewegung um die z-Achse kann sich daher nicht umkehren.
- Bei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L_z = 0} schwingt das sphärische Pendel exakt periodisch durch die Ruhelage wie ein ebenes mathematisches Pendel.
- Bei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle L_z \ne 0} ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 0 < \theta_\text{min} < 90^\circ } und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta_\text{max}<180^\circ} . Das Pendel hält sowohl vom tiefsten als auch vom höchsten Punkt der Kugel einen Mindestabstand ein. Der periodischen Auf- und Abbewegung überlagert sich eine azimutale Drehung, so dass die aufeinanderfolgende Punkte mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta=\theta_\text{min} } (ebenso auch die Punkte mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta=\theta_\text{max} } ) um ein Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Delta \phi } versetzt sind.
- Dann ist die Bewegung als ganze nur periodisch, wenn der Versatz Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Delta \phi } ein rationaler Bruchteil der vollen Drehung um 360° ist.
Mithilfe der Berücksichtigung der Konstanz der Energie kann die Bewegungsgleichung für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta} in eine Differentialgleichung erster Ordnung überführt werden, die allerdings auch nicht elementar lösbar ist:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot \theta = \pm \sqrt{\frac{2E}{mR^2} - \frac{L_z^2}{m^2 R^4 \sin^2 \theta} + 2 \frac gR \cos \theta}} ,
Konisches Pendel
Das konische Pendel wird durch die Lösung mit
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta = \text{const.} }
beschrieben. Dann ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \ddot \theta = 0 } und folglich nach der obigen Bewegungsgleichung
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot \phi^2 \cos \theta= \frac gR } .
Demnach beschreibt das Pendel mit der konstanten Winkelgeschwindigkeit
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \omega = \dot \phi= \pm \sqrt{\frac{g}{R \cos \theta}}}
einen Kegelmantel, wobei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \cos \theta > 0 } sein muss, der konstante Auslenkwinkel also auf den Bereich Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 0 \le \theta < 90^\circ} eingeschränkt ist.
Behandlung in der Newtonschen Mechanik
Allgemeiner Fall
Die Bahnkurve Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec r(t)} der Pendelmasse ergibt sich nach der Newtonschen Mechanik als Lösung der vektoriellen Differentialgleichung für die Beschleunigung
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \ddot\vec r=\vec g + \vec a_Z } .
Der erste Summand auf der rechten Seite ist die Beschleunigung aufgrund der Schwerkraft Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle m\vec g} , die hier die eingeprägte Kraft ist. Der zweite Summand rührt von der durch den Stab ausgeübten Zwangskraft her. Sie muss den Körper auf der Kugelschale mit dem Radius Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle R} halten, also – bei jeder Position und Geschwindigkeit des Körpers – die radiale Komponente der Schwerkraft aufheben und die für die Bahnkrümmung mit dem Krümmungsradius Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle R} nötige Zentripetalkraft ausüben. Die Zwangskraft wirkt daher in radialer Richtung und ist gegeben durch:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec a_Z = \left[ -(\vec g\cdot\vec e_r)-\frac{\dot\vec r^2}{R}\right] \vec e_r} .
Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec e_r=\vec r/R} bezeichnet den vom Aufhängepunkt weg gerichteten radialen Einheitsvektor.
Zusammen mit der eingeprägten Kraft kann man schreiben:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \ddot\vec r=\left[\vec g -(\vec g\cdot\vec e_r)\vec e_r\right] -\frac{\dot\vec r^2}{R} \vec e_r = - \frac{\vec r \times \left( \vec r \times \vec g \right)}{R^2} -\frac{\dot\vec r^2}{R^2} \vec r} ,
Hier zeigt sich, dass die gesamte Beschleunigung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \ddot\vec r} durch die tangentiale Komponente der Schwerkraft (Term in eckigen Klammern) und die radiale Zentripetalkraft verursacht wird.
Drückt man diese Gleichung in Kugelkoordinaten aus, ergeben sich wieder die Differentialgleichungen für die Winkel Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \phi} , die – wie oben angemerkt – nicht geschlossen gelöst werden können.[3] Für eine numerische Lösung sind kartesische Koordinaten günstiger, weil in sphärischen Koordinaten der Winkel Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \phi} am Ort der Ruhelage nicht definiert ist. Die nebenstehende Animation, die einen komplizierten Bewegungsablauf zeigt, wurde auf diese Weise mit einem SciLab Skript erstellt.[4]
Harmonische und anharmonische Näherung bei kleinen Ausschlägen
Qualitative Beschreibung
Bei kleinen Ausschlägen sind die Bewegungen des sphärischen Pendels einfach: Sind die Ausschläge infinitesimal klein, schwingt es wie ein isotroper zweidimensionaler harmonischer Oszillator mit derselben Frequenz
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \omega_0 = \sqrt{\frac{g}{R}}}
wie das ebene mathematische Pendel in harmonischer Näherung. Das heißt, die Bahnkurven sind raumfeste Ellipsen, einschließlich der Grenzfälle der linearen Schwingung und des Kreises. Hierzu siehe den speziellen Abschnitt im Artikel Harmonischer Oszillator. Bei kleinen, aber endlichen Ausschlägen treten anharmonische Effekte auf, die eine Verringerung der Umlauffrequenz und eine Präzession der Bahnellipse (im Drehsinn des Umlaufs) nach sich ziehen. Beides rührt daher, dass die Frequenz des ebenen mathematischen Pendels nur im infinitesimalen Bereich von der Größe des Ausschlags Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \theta_\text{max}} unabhängig ist, mit zunehmendem Ausschlag aber sinkt. In niedrigster Näherung gilt (siehe in mathematisches Pendel)
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \omega(\theta_\text{max}) = \omega_0 \left( 1-\frac{1}{16} \sin^2(\theta_\text{max}) \right)}
Eine elliptische Schwingung kann als Überlagerung von zwei linearen Schwingungen gleicher Frequenz mit verschieden großen Ausschlägen angesehen werden, die um eine Viertelperiode versetzt und rechtwinklig zueinander entlang der großen und kleinen Halbachse der Ellipse erfolgen. Diese Möglichkeit ist bei infinitesimal kleinen Ausschlägen gegeben, so dass die linearen Schwingungen synchron bleiben und eine raumfeste Ellipse bilden. Bei realen Auslenkungen ist aber beim sphärischen Pendel die Schwingung längs der kleinen Halbachse etwas schneller als die Schwingung längs der großen Halbachse, so dass sie schon über ihren Nullpunkt hinaus ist, wenn die andere erst bei ihrer maximalen Auslenkung, d. h. am Scheitelpunkt der Bahnkurve, ankommt. Zusammengesetzt ergibt sich, dass der Scheitelpunkt auf einem Kreis herumwandert.
Berechnung
Die Bewegungen bei kleinen Ausschlägen werden am einfachsten in kartesischen Koordinaten durch eine Entwicklung nach Potenzen behandelt. Der Ursprung liegt im Aufhängepunkt und die z-Achse ist nach unten gerichtet. Kleine Abweichungen von der Ruhelage sind definiert durch Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle |x|/R \ll 1} sowie Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle |y|/R \ll 1} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle |z-R|/ R \ll 1} . Es ergeben sich zwei gekoppelte Differentialgleichungen für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \ddot x} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \ddot y} , die in Potenzreihen entwickelt werden können.[3]
- 1. Näherung – Linearisierung
Berücksichtigt man für infinitesimale Ausschläge nur die Glieder niedrigster Potenz, erhält man zwei entkoppelte Differentialgleichungen für ein Paar harmonischer Oszillatoren gleicher Frequenz
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \begin{array}{ll} \ddot x &= - \omega_0^2 \,x \\ \ddot y &= - \omega_0^2 \,y \end{array}}
Für Lösungsweg und Lösung siehe harmonischer Oszillator. Dieselben Differentialgleichungen erhält man für kleine Auslenkungen aus der physikalisch begründeten Näherung, dass die Bewegung sich nur in der Ebene Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z=R} abspielt und die zur Ruhelage Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle x=y=0} hin rücktreibende Kraft durch die tangentiale Komponente der Schwerkraft gegeben ist, wobei diese (für Auslenkung in x-Richtung) durch
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle -m\, g\, \sin(\theta) \approx -m\, g\, \theta \approx - \frac{m\, g}{R} x}
genähert wird (für y-Richtung entsprechend). Die Bahnkurven sind raumfeste Ellipsen mit beliebiger Orientierung der Achsen in der Schwingungsebene, einschließlich der Grenzfälle Strecke und Kreis.
- 2. Näherung – kubische Glieder, rotierende Ellipse
In nächster Näherung treten kubische Glieder auf, über die die beiden Differentialgleichungen auch gekoppelt sind. Eine geschlossene Lösung ist nicht möglich. Eine Näherungslösung geht, der obigen qualitativen Diskussion entsprechend, vom Ansatz einer langsam rotierenden Ellipsenbahn aus. Demnach durchläuft der Pendelkörper eine Ellipse mit den Halbachsen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle a} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle b} mit der Kreisfrequenz
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \omega(a,b) = \omega_0 \left( 1-\frac18 \frac{a^2 + b^2}{R^2} \right)} .
Die Ellipse rotiert dabei im Sinn des Umlaufs so, dass der Scheitelpunkt bei jedem Umlauf um den Winkel
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Delta \phi(a,b) = \frac34 \frac{\pi a b}{R^2}}
versetzt wird. Das entspricht einer Drehung der Bahn mit einer Winkelgeschwindigkeit
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Omega(a,b) = \frac38 \omega_0 \frac{a b}{R^2}} .
Diese Präzessionsbewegung ist zum Beispiel eine häufige Störung beim Foucaultschen Pendel, weil sie leicht die Größe der Präzession aufgrund der Erddrehung erreicht.[5]
Einzelnachweise
- ↑ Bergmann-Schaefer Lehrbuch der Experimentalphysik, Band 1: Mechanik, Akustik, Wärme, IV. Kapitel, Abschnitt 35
- ↑ A. Budó: Theoretische Mechanik. 4. Auflage. VEB Deutscher Verlag der Wissenschaften, Berlin 1967, § 23, S. 117–118.
- ↑ a b M.G. Olsson: Spherical pendulum revisited, American Journal of Physics 49, 531 (1981); doi: 10.1119/1.12666
- ↑ Hier ist das SciLab-Skript angegeben.
//SKRIPTANFANG. Dieses Skript ist mit der freien Mathematik-Software "Scilab" als sce-Datei ausführbar
// Autor Modalanalytiker 20.08.2018 //Numerische Lösung der vektoriellen Bewegungsgleichung des sphärischen Pendels //Bahnkurve für Wikimedia-Bild "SphaerischesPendelxyz.svg" //Durchgängig SI-Einheiten! d2r=%pi/180;//degree to radian
//########## [ EINGABE ############################################## g=[0;0;-9.81]//Erdbeschleunigung, kartesische Koordinaten L=1//Pendellänge TSim=3// Dauer der Simulation deltaT=0.01//Zeitschritt der Ausgabe //Anfangszustand in sph. Koord., um unrealistische Zustände auszuschließen Th0=90*d2r; //Polarwinkel dTh0=-200*d2r; //Polarwinkelgeschwindigkeit Ph0=0*d2r; //Azimut dPh0=90.*d2r; //Azimutwinkelgeschwindigkeit //########## ] EINGABE ENDE #########################################
function dz=xyzPendel(t,z,g)// //Ableitung des Zustandsvektors z=[x;y;z;dx;dy;dz], 6 x 1 //g=[gx; gx; gz]: Schwerebeschleunigung dz(1:3)=z(4:6)//Geschwindigkeit 3 x 1 r=norm(z(1:3))//Abstand der Masse vom Ursprung, Länge der Pendelstange, 1 x 1 er=z(1:3)/r //Radialer Einheitsvektor der Massenposition, 3 x 1 dz(4:6)=g-er*(g'*er + z(4:6)'*z(4:6)/r) //Beschleunigung, 3 x 1 endfunction
t=0:deltaT:TSim;//Zeitspanne der Simulation //Umrechnung der AW in kartesische Koordinaten //Anfangsposition x0=L*sin(Th0)*cos(Ph0); y0=L*sin(Th0)*sin(Ph0); z0=L*cos(Th0); Z0(1:3)=[x0;y0;z0];//Anfangsposition der Masse, Großbuchstabe Z0! //Anfangsgeschwindigkeit eTh0=[cos(Th0)*cos(Ph0);cos(Th0)*sin(Ph0);-sin(Th0)]//polarer Einheitsvektor ePh0=[-sin(Ph0);cos(Ph0);0]//azimutaler Einheitsvektor Z0(4:6)=L*dTh0 *eTh0 + L*sin(Th0)*dPh0 *ePh0//Anfangsgeschwindigkeit der Masse
//Integration der Dgl. //Solver automatically selects between nonstiff predictor-corrector Adams method and //stiff Backward Differentiation Formula (BDF) method zk=ode(Z0, t(1), t, list(xyzPendel,g));//###### Dgl.-Lösung ###########
//Graphik xdel(); param3d(zk(1,:),zk(2,:),zk(3,:))//#################################### title('Sphärisches Pendel mit Aufhängung im Ursprung ""o""') ce=gce(); ce.foreground=5; ce.thickness=3; //rote Kurve param3d(zk(1,1),zk(2,1),zk(3,1))//ausgefüllter roter Kreis für Anfang ce=gce();ce.mark_mode="on"; ce.mark_style=0; ce.mark_size_unit = "point"; ce.mark_size = 10; ce.mark_foreground = 5; ce.mark_background = 5; param3d(zk(1,:),zk(2,:),-L*ones(zk(1,:)))//Projektion z=konst ce=gce(); ce.foreground=13; ce.line_style=9;ce.thickness=2; param3d(zk(1,1),zk(2,1),-L*ones(zk(1,1)))//Projektion z=konst, Anfang ce=gce();ce.mark_mode="on"; ce.mark_style=0; ce.mark_size_unit = "point"; ce.mark_size = 5; ce.mark_foreground = 5; ce.mark_background = 5; param3d(-L*ones(zk(1,:)),zk(2,:),zk(3,:))//Projektion x=konst ce=gce(); ce.foreground=5; ce.line_style=9;ce.thickness=2; param3d(-L*ones(zk(1,1)),zk(2,1),zk(3,1))//Projektion, Anfang ce=gce();ce.mark_mode="on"; ce.mark_style=0; ce.mark_size_unit = "point"; ce.mark_size = 5; ce.mark_foreground = 5; ce.mark_background = 5; param3d(zk(1,:),+L*ones(zk(1,:)),zk(3,:))//Projektion y=konst ce=gce(); ce.foreground=5; ce.line_style=9;ce.thickness=2; param3d(zk(1,1),+L*ones(zk(1,1)),zk(3,1))//Projektion, Anfang ce=gce();ce.mark_mode="on"; ce.mark_style=0; ce.mark_size_unit = "point"; ce.mark_size = 5; ce.mark_foreground = 5; ce.mark_background = 5; param3d(0,0,0) //schwarzes O-Zeichen für Ursprung ce=gce();ce.mark_mode="on"; ce.mark_style=9; ce.mark_size_unit = "point"; ce.mark_size = 5; ce.mark_foreground = 1; ca=gca(); ca.rotation_angles=[83.25 280]; ca.data_bounds = 1.03*[-1,-1,-1.;1,1,1]; ca.tight_limits = ["on","on","on"];
//SKRIPTENDE - ↑ Szostak, Roland: Ein permanent schwingendes Foucault-Pendel für Schulen, Der Mathematische und Naturwissenschaftliche Unterricht, PLUS LUCIS 2/2002-1/2003, S. 11–15, html
Weblinks
- Physik und Mathematik mit Maple, Kugelpendel, abgerufen am 22. Dezember 2014