Gelöste Aufgaben/ODE1

Aufgabenstellung

Hier untersuchen wir systematisch ein System aus einem starren Pendel im Erdschwerefeld.

Lageplan

Gesucht sind Gleichgewichtslagen, Schwingungen um diese und die numerische Lösung der nichtlinearen Bewegungsgleichung.

Lösung mit Maxima

bzgl A ist J_A, der Schwerpunkt des Stabes liegt in seiner Mitte.

Header

Wir nehemn an: das Massenmoment des homogenen Stabes bzgl A ist JA, der Schwerpunkt des Stabes liegt in seiner Mitte.

/*******************************************************/ /* MAXIMA script */ /* version: wxMaxima 15.08.2 */ /* author: Andreas Baumgart */ /* last updated: */ /* ref: TM-C */ /* description: solve the IVP, check for stability */ /*******************************************************/

tmp

Aus dem Prinzip der virtuellen Verrückungen kommt die Gleichgewichtsbeziehung für den Stab:

$\delta W=\displaystyle \underbrace {-J_{A}\cdot {\ddot {\varphi }}\cdot \delta \varphi -m\cdot g\cdot {\frac {\ell }{2}}\cdot \sin(\varphi )\cdot \delta \varphi } _{\displaystyle =\delta W^{a}}-\underbrace {0} _{\displaystyle =\delta \Pi }$,

wobei

$\varphi =\varphi (t)$

Mit der neuen, dimensionslosen Zeit

$\tau =\displaystyle {\frac {t}{t_{Bez}}}{\text{ und }}t_{Bez}^{2}={\frac {J_{A}}{m\cdot g\cdot {\frac {\ell }{2}}}}$

wird aus der Bewegungsgleichung

$\displaystyle {\frac {d^{2}}{d\tau ^{2}}}(\varphi )+\sin(\varphi )=0$

Declarations

Text

1+1

tmp

Die Gleichung

$\displaystyle {\frac {d^{2}}{d\tau ^{2}}}\varphi =0$

ist erfüllt für

$\sin(\varphi ){\stackrel {!}{=}}0{\text{, also }}\varphi =0,\pi ,2\,\pi ,\,\ldots$

Gleichgewichtslagen

Die Funktion der Rückstell'kraft' sin(φ) schneidet die y-Achse dabei mit der Steigung +1 (z.B. bei φ₁ = 0, φ₃ = 2π) und mit der Steigung  -1 (z.B. bei φ₂ = π).

Das Vorzeichen der Steigung wird über die Stabilität der Lösung entscheiden!

Equilibrium Conditions

Text

1+1

tmp

Equilibrium Conditions and Stability

Text

1+1

tmp

Solving

Text

1+1

tmp

Post-Processing

Text

1+1

Phasendiagramm

Lösung im Zeitbereich

Lösung im Zeitbereich.

Lösung im Zeitbereich

---

Links

• ...

Literature

• ...