Gelöste Aufgaben/Kv52

Aufgabenstellung

Lageplan

Ein starrer Stab AB (Massenmoment J_A, Länge ℓ₁) wird aus dem Winkel φ₀ im Erdschwerefeld losgelassen und stößt in C auf einen Anschlag. Der Stoß zwischen Stab und Oberfläche sei ideal-elastisch.

Gesucht ist die nichtlineare Bewegungsgleichung und die numerische Lösung als Anfangswertproblem. Dabei denken wir uns den Anschlag als elastische Feder, die nur für φ<-π/2 Kontakt zum Stab hat. Die Federkraft ist also Null, solange der Stab die Oberfläche in C nicht berührt und sie ist proportional zur Federkompression w, wenn sich Kugel und Oberfläche berühren.

Lösung mit Maxima

Lorem Ipsum ....

tmp

Wir stellen die Bewegungsgleichungen des Systems als System con Differentialgleichungen erster Ordnung auf. Die Nichtlinearität kommt aus den großen Winkeln φ(t) und dem Kontakt mit der Wand.

Mit unterschiedlichen Steifigkeiten für den Kontakt testen wir die Möglichkeiten der numerischen Integration aus.

Header

Text

1+1

tmp

Die System-Parameter sind

${\begin{array}{l}\ell _{2}={\frac {3}{4}}\ell _{1}\\J_{A}={\frac {1}{3}}m\;\ell _{1}^{2}\end{array}}$.

Zum Dimensionslos-Machen der Bewegungsgleichungen brauchen wir später eine Bezugszeit t_B, die wir mit Hilfe der Eigenfrequenz der zugeordneten linearen Systems so wählen:

$t_{B}=1/\omega _{0}$

(Achtung: das macht 2π-periodische Lösungen)

Die Eigenkreisfrequenz bei Schwingungen um φ =- π/2 herum ist

$\displaystyle \omega _{0}^{2}={\frac {m\,g\,\ell }{2\,J_{A}}}$

Außerdem wählen wir eine dimensionslose Federsteifigkeit κ, so dass

$\displaystyle k=\kappa \cdot {\frac {m\,g}{\ell _{2}}}$

Kennlinie

Für den nichtlinearen Kontakt wählen wir eine Kennlinie wie in Kw23 zu

$C({\tilde {u}},\varepsilon )=\left\{{\begin{array}{ll}0&{\text{ wenn }}{\tilde {u}}<-\varepsilon \\{\frac {1}{2}}({\tilde {u}}+\varepsilon )^{2}&{\text{ wenn }}-\varepsilon <{\tilde {u}}<\varepsilon {\text{ und }}\\{\tilde {u}}&{\text{ sonst }}\end{array}}\right.$

Und so wie rechts im Bild sieht sie dann aus:

Declarations

Text

1+1

tmp

Die Gleichgewichtsbeziehungen erhalten wir aus einem Momentengleichgewicht um den Punkt A.

Das Freikörperbild

liefert die Bewegungsgleichung

$\displaystyle J_{A}{\ddot {\varphi }}-K(\varphi )\,\ell _{2}+mg{\frac {\ell _{1}}{2}}\cos(\varphi )=0$

Die Kontaktkraft konstruieren wir mit Hilfe der nichtlinearen Kennlinie C und

$\displaystyle {\tilde {u}}=-{\frac {\pi }{2}}-\varphi$

zu

$K(\varphi )=\ell _{2}\cdot C(-\pi /2-\varphi ,\varepsilon )$.

Mit der dimensionslosen Zeit

$\tau =t/t_{B}$

finden wir als Bewegungsgleichung des Systems

$\displaystyle {\frac {d^{2}}{d\tau ^{2}}}\varphi +\cos(\varphi )-{\frac {3}{2}}\,\kappa \cdot C(-\pi /2-\varphi ,\varepsilon )=0$.

Equilibrium Conditions

Text

1+1

tmp

Umschreiben in ein Differentialgleichungssystem erster Ordnung - wie wir es für die numerische Lösung brauchen - liefert

$\left({\begin{array}{c}{\dot {\varphi }}\\{\dot {\omega }}\end{array}}\right)=\left({\begin{array}{c}\omega \\\displaystyle {\frac {3}{2}}\,\kappa \cdot C(-\pi /2-\varphi ,\varepsilon )-\cos(\varphi )\end{array}}\right)$.

Mit den Anfangsbedingungen

${\begin{array}{l}\varphi (0)=0\\\omega (0)=0\end{array}}$

und einer Lösungsroutine nach dem Runge-Kutta-Verfahren 4.ter Ordnung liefert das numerische Ergebnis.

Solving

Text

1+1

tmp

Wir tragen hier zwei Ergebnisse auf:

... in dieser Spalte den Referenzfall für $\kappa =100.$	... und hier für $\kappa =1000.$
zunächst im Zeitbereich
und dann im Phasenraum:

Im Referenzfall erhalten wir - wie erhofft - eine periodische Lösung: die Phasenkurve ist geschlossen. Schon für eine um den Faktur 10 größere Steifigkeit des Anschlags jedoch weicht die Lösung sichtbar von einer periodischen ab. Das klassische Problem von numerischen Integrationsroutinen wird hier sichtbar.

Post-Processing

Text

1+1

---

Links

• ...

Literature

• ...