Gelöste Aufgaben/Kw24

Aufgabenstellung

Ein Mechanismus besteht aus einem starren Stab AC (Länge ℓ₁, Massenmoment J) und einer Feder BC (Länge ℓ₂ im entspannten Zustand, Steifigkeit k).

Die Körper sind drehbar gelagert, der Abstand AB ist ℓ₃.

Lageplan

Hier geht es um die verschiedenen Lösungstypen des Systems - in Abhängigkeit von den Anfangsbedingungen. Gesucht sind Lösungen als Anfangswertproblem, nichtlineare Schwingungen und Aussagen über die Stabilität der Bewegungen.

Lösung mit Maxima

Lorem Ipsum ....

tmp

Gesucht ist die Lösung zu dem nichtlinearen Anfangswertproblem.

Header

Text

📑 toggle code

1+1

tmp

Wir wählen die Längenverhältnisse

${\ell }_2={\ell }_1,{\ell }_3=\sqrt{2}\cdot {\ell }_1$

so dass die Feder für φ(t) = 45° entspannt ist.

Declarations

Text

📑 toggle code

1+1

tmp

Der Mechanismus hat genau einen Freiheitsgrad, den wir entweder durch die Federlängung w oder den Drehwinkel φ erfassen. Den Zusammen zwischen den beiden stellen wir über

${\ell }_1^2\cdot {\mathrm{s}\mathrm{i}\mathrm{n}\left(\mathrm{\phi }\left(t\right)\right)}^2+{({\ell }_3-{\ell }_1\cdot \mathrm{c}\mathrm{o}\mathrm{s}\left(\mathrm{\phi }\left(t\right)\right))}^2={({\ell }_2+w)}^2$

her. Nun steht uns frei, ob wir diese Beziehung nach w oder φ auflösen. Wir wählen φ als Minimal-Koordinate und erhalten

$w=\frac{{\ell }_1\cdot \sqrt{12-2^{\frac{7}{2}}\cdot \cos (\phi )}-2\cdot {\ell }_1}{2}$

Nach den Regeln der Variationsmethoden erhalten wir daraus:

$\delta w=\frac{2^{\frac{3}{2}}\cdot {\ell }_1\cdot \sin (\phi )\cdot \delta \phi }{\sqrt{12-2^{\frac{7}{2}}\cdot \cos (\phi )}}$.

Kinematics

Text

📑 toggle code

1+1

tmp

Mit dem Prinzip der virtuellen Verrückungen schreiben wir

$\delta W=-(\frac{d^2}{d{t}^2}\cdot \mathrm{\phi }\left(t\right))\cdot J\cdot \delta \phi -k\cdot w\cdot \delta w$

und erhalten die Bewegungsgleichung

$J\cdot \ddot{\phi }+\sqrt{2}\cdot l_1^2\cdot k\cdot \sin \phi -\frac{2^{\frac{3}{2}}\cdot l_1^2\cdot k\cdot \sin \phi }{\sqrt{12-2^{\frac{7}{2}}\cdot \cos \phi }}=0$

Equilibrium Conditions

Text

📑 toggle code

1+1

tmp

Die Bewegungsgleichung machen wir dimensionslos, indem wir charakteristische Größen aus der linearisierten Bewegungsgleichung verwenden.

Die drei Gleichgewichts-Lagen bzgl. derer wir die Bewegungsgleichung linearisieren können, sind

$\begin{array}{ll}{\phi }_{0,1}& =-\frac{\pi }{4}\\ {\phi }_{0,2}& =0\\ {\phi }_{0,3}& =+\frac{\pi }{4}\end{array}$.

Eine Taylor-Reihenentwicklung des Rückstellmoments bzgl. φ_0,3 ,also

$\phi ={\phi }_{0,3}+\varphi$

liefert die linearisierte Bewegungsgleichung

$\ddot{\varphi }+{\omega }_0^2\cdot \varphi =0\text{ mit }{\omega }_0^2=\frac{{\ell }_1^2\cdot k}{J}$.

Dabei ist ω₀ die Eigenkreisfreuquenz der linearisierten Bewegungsgleichung.

Equilibrium Points and Stability

Text

📑 toggle code

1+1

tmp

Eine Bezugszeit für unsere nichtlineare Bewegungsgleichung kommt nun aus

${\omega }_0=\frac{2\pi }{t_{Bez}}$

Damit lautet die dimensionslose Bewegungsgleichung

$\frac{d^2}{d{\tau }^2}\mathrm{\Phi }+\frac{2^{\frac{5}{2}}\cdot {\pi }^2\cdot \sqrt{12-2^{\frac{7}{2}}\cdot \mathrm{c}\mathrm{o}\mathrm{s}\left(\mathrm{\Phi }\right)}\cdot \mathrm{s}\mathrm{i}\mathrm{n}\left(\mathrm{\Phi }\right)-2^{\frac{7}{2}}\cdot {\pi }^2\cdot \mathrm{s}\mathrm{i}\mathrm{n}\left(\mathrm{\Phi }\right)}{\sqrt{12-2^{\frac{7}{2}}\cdot \mathrm{c}\mathrm{o}\mathrm{s}\left(\mathrm{\Phi }\right)}}=0\text{ mit }\tau =t/t_{Bez}$

Rückstellmoment der Feder

Den nichtlinearen (zweiten) Summand können wir als dimensionsloses "Rückstellmoment auffassen und über dem Drehwinkel auftragen:

Nullstellen des Rückstellmoments zeigen die Gleichgewichtslagen des Systems an. Die Steigung dieses Rückstellmoments in den Gleichgewichtslagen deutet dabei auf die Stabilität der linearisierten Bewegung (oder Bewegung im Kleinen) hin. So ist die Lösung Φ = 0 eine instabile Gleichgewichtslage, die beiden anderen sind stabil.===Make Equations of Motion Dimensionless=== Text

📑 toggle code

1+1

tmp

Für die numerische Lösung müssen wir die Bewegungsgleichung

$\ddot{\mathrm{\Phi }}=f(\mathrm{\Phi })$

umschreiben auf eine Differentialgleichung erster Ordnung:

$\left(\begin{array}{c}\mathrm{\Phi }\\ \mathrm{\Omega }\end{array}\right)=\left(\begin{array}{c}\mathrm{\Omega }\\ f(\mathrm{\Phi })\end{array}\right)$

Diese lösen wir mit einer fertigen Routine - hier einem Runge-Kutta-Verfahren - als Anfangswertproblem. Anders als bei linearisierten Bewegungsgleichungen hängt die Charakteristik der Lösung von den Anfangsbedingungen ab.

Numerical Solution

Text

📑 toggle code

1+1

tmp

Das schauen wir uns am Verlauf von Lösungen im Zeitbereich an:

Lösung f) fällt aus dem Raster - hier gibt es keine periodische Lösung mehr: Bei einer anderen Skalierung der Φ-Achse erkennen wir, dass der Stab AC um den Punkt A rotiert:

Lösung f) im Zeitbereich.

Oft erkennt man das Charakteristische einer Lösung im Phasendiagramm. Hier wird die Geschwindigkeit Ω über dem Drehwinkel Φ aufgetragen:

Phasendiagramme der Lösungen

Periodische Lösungen - die oft gesucht sind - erkennt man an einem geschlossenen Verlauf. Die hellblaue Lösung gehört zu den Anfangswerten von f), bei dem der Drehwinkel immer größer wird - sich also keine periodische Lösung ergibt.

Wir sehen, dass die Lösung e) (schwarz) eine "gerade noch" periodische Lösung ist - etwas mehr Anfangsgeschwindigkeit im System und die Lösung sieht wie in f) aus. Diese Lösung nennt man auch Separatrix - sie separiert zwei Charakteristiken von Lösungen.

Analog sehen Sie in der Lösung c) (grün) eine Separatrix zwischen der periodischen Lösung der kleinen Schwingung um die Gleichgewichtslage und solchen periodischen Lösungen, bei denen der Stab die Horizontale durchbricht:

Phasendiagramme

Post-Processing

Text

📑 toggle code

1+1

a)	b)
c)	d)
e)	f)

---

Links

• ...

Literature

• ...