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.

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 nicht vor und ist daher eine zyklische Variable. Dies ist ein Beispiel für das Noether-Theorem.)

Damit lässt sich aus der Differentialgleichung für eliminieren:

Diese Bewegungsgleichung für 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 berücksichtigt wird, die sich daraus ergibt, dass keine explizite Abhängigkeit von der Zeit vorliegt. Daraus folgt:

• Die Bewegung ist auf einen Bereich eingeschränkt, findet also zwischen zwei Breitenkreisen statt.
• Die Auf- und Abbewegung zwischen den Breitenkreisen ist periodisch (aber nicht harmonisch).
• Die azimutale Winkelgeschwindigkeit ist entweder konstant Null (wenn ) oder hat das feste Vorzeichen von . Der Drehsinn der Pendelbewegung um die z-Achse kann sich daher nicht umkehren.
• Bei schwingt das sphärische Pendel exakt periodisch durch die Ruhelage wie ein ebenes mathematisches Pendel.
• Bei ist und . 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 (ebenso auch die Punkte mit ) um ein versetzt sind.
• Dann ist die Bewegung als ganze nur periodisch, wenn der Versatz ein rationaler Bruchteil der vollen Drehung um 360° ist.

Mithilfe der Berücksichtigung der Konstanz der Energie kann die Bewegungsgleichung für in eine Differentialgleichung erster Ordnung überführt werden, die allerdings auch nicht elementar lösbar ist:

Konisches Pendel

Das konische Pendel wird durch die Lösung mit

beschrieben. Dann ist und folglich nach der obigen Bewegungsgleichung

Demnach beschreibt das Pendel mit der konstanten Winkelgeschwindigkeit

einen Kegelmantel, wobei sein muss, der konstante Auslenkwinkel also auf den Bereich eingeschränkt ist.

Behandlung in der Newtonschen Mechanik

Allgemeiner Fall

Bahnkurve eines sphärischen Pendels

Zum Animieren klicken

---

Das Pendel (mit einer Stange statt eines Fadens zwischen Masse und Aufhängepunkt) wird aus horizontaler Lage mit einer nach oben und hinten gerichteten Geschwindigkeit angestoßen.

Die Bahnkurve der Pendelmasse ergibt sich nach der Newtonschen Mechanik als Lösung der vektoriellen Differentialgleichung für die Beschleunigung

Der erste Summand auf der rechten Seite ist die Beschleunigung aufgrund der Schwerkraft , 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 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 nötige Zentripetalkraft ausüben. Die Zwangskraft wirkt daher in radialer Richtung und ist gegeben durch:

bezeichnet den vom Aufhängepunkt weg gerichteten radialen Einheitsvektor.

Zusammen mit der eingeprägten Kraft kann man schreiben:

Hier zeigt sich, dass die gesamte Beschleunigung 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 und , die – wie oben angemerkt – nicht geschlossen gelöst werden können. Für eine numerische Lösung sind kartesische Koordinaten günstiger, weil in sphärischen Koordinaten der Winkel am Ort der Ruhelage nicht definiert ist. Die nebenstehende Animation, die einen komplizierten Bewegungsablauf zeigt, wurde auf diese Weise mit einem SciLab Skript erstellt.

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

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 unabhängig ist, mit zunehmendem Ausschlag aber sinkt. In niedrigster Näherung gilt (siehe in mathematisches Pendel)

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 sowie und . Es ergeben sich zwei gekoppelte Differentialgleichungen für und , die in Potenzreihen entwickelt werden können.

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

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 abspielt und die zur Ruhelage hin rücktreibende Kraft durch die tangentiale Komponente der Schwerkraft gegeben ist, wobei diese (für Auslenkung in x-Richtung) durch

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.

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 und mit der Kreisfrequenz

Die Ellipse rotiert dabei im Sinn des Umlaufs so, dass der Scheitelpunkt bei jedem Umlauf um den Winkel

versetzt wird. Das entspricht einer Drehung der Bahn mit einer Winkelgeschwindigkeit

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.

Einzelnachweise

1. Bergmann-Schaefer Lehrbuch der Experimentalphysik, Band 1: Mechanik, Akustik, Wärme, IV. Kapitel, Abschnitt 35
2. A. Budó: Theoretische Mechanik. 4. Auflage. VEB Deutscher Verlag der Wissenschaften, Berlin 1967, § 23, S. 117–118. 
3. M.G. Olsson: Spherical pendulum revisited, American Journal of Physics 49, 531 (1981); doi: 10.1119/1.12666
4. 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
5. Szostak, Roland: Ein permanent schwingendes Foucault-Pendel für Schulen, Der Mathematische und Naturwissenschaftliche Unterricht, PLUS LUCIS 2/2002-1/2003, S. 11–15, PDF

Weblinks

• Physik und Mathematik mit Maple, Kugelpendel, abgerufen am 22. Dezember 2014