Benutzer:FridtjofNansen/DoppelPendel
Herleitung der Bewegungsgleichung des Doppelpendels
massenlos, reibungsfrei, starr! Winkel sind wie definiert? Skizze!
Um die Bewegungsgleichung des Doppelpendels herzuleiten muss zuerst das Koordinatensystem definiert werden. Der Ursprung des Systems liegt bei (0,0) und ist auch der Ursprung des ersten Pendelarmes 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 l_1 } . Die y-Achse ist als positiv nach oben definiert.
Die Bewegungsgleichungen werden mittels der Positionen der Massen 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_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 m_2 } hergeleitet. Jedoch werden dann die Winkel betrachtet, da diese die Bewegung genau so gut beschreiben jedoch die Anzahl der Gleichungen reduzieren. Die Winkel werden genutzt um die Geschwindigkeit und damit die kinetische und potentielle Energie der Massen zu bestimmen. Die beiden Energien ergeben zusammen die Lagrange-Funktion des Systems welche dann in die Euler-Lagrange Gleichung gesetzt wird und damit zwei Gleichungen geben, die Entwicklung der Winkel mit der Zeit angeben.
Gleichungen
Position m1
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_1 = l_1 sin(\theta_1) }
Fehler beim Parsen (Konvertierungsfehler. Der Server („https://wikimedia.org/api/rest_“) hat berichtet: „Cannot get mml. Server problem.“): {\displaystyle y_{1}=-l_{1}cos(\theta _{1})}
Position m2
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_1 = l_1 sin(\theta_1) + l_2 sin(\theta_2)}
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_1 = -l_1 cos(\theta_1) -l_2 cos(\theta_2)}
Geschw m1
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 u_1 = \frac{\partial x_1}{\partial \theta_1} = \dot{\theta_1} l_1 cos(\theta_1)}
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 v_1 = \frac{\partial y_1}{\partial \theta_1} = \dot{\theta_1} l_1 sin(\theta_1)}
Geschw m2
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 u_2 = \frac{\partial x_2}{\partial \theta_2} + u_1 = \dot{\theta_2} l_2 cos(\theta_2) + u_1}
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 L = T - V} (der Lagrange-Gleichung) 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 T_1 = \frac{1}{2} m_1 (u_1^2 + v_1^2)} 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 T_2 = \frac{1}{2} m_2 (u_2^2 + v_2^2)} 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 V_1 = m_1 g y_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 V_2 = m_2 g y_2} folgt:
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 T = T_1 + T_2 = \frac{1}{2} \dot{\theta_1^2} l_1^2+ \frac{1}{2} (m_2 \dot{\theta_2^2}l_2^2 + m_2 \dot{\theta_1^2} l_1^2 + 2 m_2 \dot{\theta_1} l_1 \dot{\theta_2} l_2 cos(\theta_1 - \theta_2))} 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 V = V_1 + V_2 = -m_1 g l_1 cos(\theta_1) -m_2 g l_1 cos(\theta_1) - m_2 g l_2 cos(\theta_2)} .
Zusammen ergibt das für die Lagrange-Funktion:
Mittels der Euler-Lagrange-Gleichung (ist das der richtige Link???) :
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 \frac{\partial}{\partial t} \frac{\partial L}{\partial \dot{\theta_i}} - \frac{\partial L}{\partial \theta_i} = 0}
lässt sich dann erhalten:
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_1} = -\frac{m_2}{m_1 + m_2} \frac{l_2}{l_1} (\ddot{\theta_2} cos(\theta_1 - \theta_2) + \dot{\theta_2^2} sin(\theta_1 - \theta_2)) - \frac{g}{l_1} sin(\theta_1)} 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_2} = - \frac{l_1}{l_2} (\ddot{\theta_1} cos(\theta_1 - \theta_2) - \dot{\theta_1^2} sin(\theta_1 - \theta_2)) - \frac{g}{l_2} sin(\theta_2)} .
Dies sind die Winkelbeschleunigungen für die 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_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 \theta_2} und beschreibung die Evolution des Systems mit der Zeit in Abhängigkeit von den Längen, Massen und der Erdbeschleunigung.
Evaluierung des chaotischen Verhaltens
Zur Betrachtung des chaotischen Verhaltens des Doppelspendels gibt es eine Reihe von Möglichkeiten. Oft kann mittels einfachster Berechnungen eine Aussage über chaotisches Verhalten getroffen werden.
MLE
Der MLE ist der sog. maximaler/größter Ljapunow-Exponent (maximum Lyapunov exponent) und beschreibt die "Stärke" des chaotischen Verhaltens. Er ist Bestandteil des Ljapunow-Spektrums welches alle Ljapunow-Exponenten (je einer pro Freiheitsgrad). Man geht davon aus, dass das System eine Störung in der Richtung des MLE hat und da er das größte Wachstum zeigt, ist zu erwarten, dass der MLE nach einer gewissen Zeit die Evolution des Systems dominiert. Ein positiver MLE zeigt normalerweise ein chaotisches System an. Er wird approximiert 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 \lambda_{max}(t) = \frac{1}{t} \ln \left(\frac{||\delta (t)||}{|| \delta (t=0)||} \right)} .
Bei zwei Experimenten mit einer anfänglichen Separation 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/“:): {\textstyle \delta(t=0) \backsim \mathcal{O}(10^{-8} )} in den Anfangsbedingungen oder sogar weniger, verstärkt sich diese Differez exponetiell und lässt die Trajektorien divergieren.[1] Die Separation (der natürliche Logarithmus der obigen Gleichung) kann dann in einem halblogarithmischen Diagramm gegen die Zeit aufgetragen werden. Dann wird mittels linearer Regression die Steigung bestimmt und diese gibt dann den approximierten MLE.
Bifurkationsdiagram
Bifurkationsdiagramme sind eine Möglichkeit komplexe Informationen über den Phasenraum eines dynamischen Systems in einen zweidimensionalen, visualisierbaren Plot zu komprimieren. Üblicherweise wird die qualitative Änderungen des Verhaltens eines Systems mittels der Variation eines geeigneten Parameters untersucht. So können für das Doppelpendel bspw. das Verhältnis der Massen, das Verhältnis der Längen, die Erdbeschleunigung oder die Anfangsbedingungen herangezogen werden. Durch die kontinuierliche Veränderung des gewählten Bifurkationsparameters wird das System auf Stabilität (periodische, quasi-periodische Lösungen) bzw. auf Chaos geprüft.
Wenn man die anfänglichen 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/“:): {\textstyle \theta_1 = \theta_2 } als den Bifurkationsparameter wählt, lässt sich das qualitativ veränderliche Verhalten des Doppelpendels sehr gut veranschaulichen. Dazu werden die beiden Winkel simultan Stück für Stück erhöht und für jedes Inkrement wird das Doppelpendel erneut integriert (berechnet). Mit diesen Daten lässt sich dann veranschaulichen, wie das System schwingt. Man hat also einen vier-dimensionalen Phasenraum der sich aus 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/“:): {\textstyle \theta_1, \theta_2, \dot{\theta_1}, \dot{\theta_2}} zusammensetzt. Praktischerweise oszillieren die Winkelgeschwindigkeiten 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/“:): {\textstyle \dot{\theta_1}, \dot{\theta_2}} , obwohl mit unbestimmter Amplitude, um Null. Daher ist zu erwarten, dass beide immer wieder die Null überqueren. Für ein harmonisch schwingendes System (periodische Lösung) sind die Nullüberquerungen 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 \dot{\theta_1}, \dot{\theta_2}} an festen Punkten, da das System immer an bestimmten Punkten () seine Auf- und Abwärtsbewegung beendet und zurückschwingt. Das ist vergleichbar mit einem normalen starren Pendel. Daher ist im Umkehrschluss zu erwarten, dass das chaotisch schwingende System an allen möglichen Punkten (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/“:): {\textstyle \theta_1, \theta_2} ) die Winkelgeschwindigkeit Null zeigt. Wenn man dann eine "Scheibe" aus dem Phasenraum gesondert betrachtet, bspw. 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/“:): {\textstyle \dot{\theta_2} = 0} , kann man die Bifurkation des Verhaltens zwei-dimensional darstellen indem man die 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/“:): {\textstyle \dot{\theta_1}} gegen die veränderliche Anfangsbedingung aufträgt (siehe rechts).
- ↑ Wolf, A., Swift, J. B., Swinney, H. L., & Vastano, J. A. (1985). Determining lyapunov exponents from a time series. Physica D: Nonlinear Phenomena, Vol. 16 No. 3: 285 - 317
