Gelöste Aufgaben/Kw28

Aufgabenstellung

Das System besteht aus einer Kugel (Radius r, Masse m₂) und einer Plattform (Masse m₁). Wie skizziert ist die Plattform mit einer Parallelführung aus zwei Euler-Bernoulli-Balken (Biegesteifigkeit EI) elastisch gelagert. Aus der statischen Referenzkonfiguration wird die Kugel aus der Höhe H über der Plattform losgelassen. Kugel und Plattform stoßen also aufeinander und führen dann Ihre eigene vertikale Bewegung durch - bis zur nächsten Kollision. Der Stoß zwischen Kugel und Oberfläche sei ideal-elastisch.

Lageplan

Für das skizzierte System modellieren Sie die Kugel als elastisch, die elastisch gelagerte Plattform als starr. Gesucht ist eine numerische Lösung als Anfangswertproblem und die nichtlinearen Schwingungen der beiden Systemteile.

Lösung mit Maxima

Lorem Ipsum ....

tmp

Die Kugel können Sie sich im unteren Teil durch eine Feder der Steifigkeit k₂ ersetzt denken - das geht analog zu Beispiel Kw23.

Die Federkraft K ist also Null, solange die Kugel die Oberfläche nicht berührt und sie ist proportional zur Federkompression w, wenn sich Kugel und Oberfläche berühren.

Kugel-Modell: elast. Kontakt mit Einfederung w.

Die beiden Körper haben jeweils einen Freiheitsgrad in vertikale Richtung. Die Bewegungsgleichungen für die beiden Körper sind stückweise linear (Kontakt: K = k w / kein Kontakt K = 0). Den Kontakt erfassen wir durch eine Kennlinie, die wir zwischen den beiden linearen Bereichen ausrunden. Das Ausrunden macht die numerische Integration schneller.

Wir lösen das Anfangswertproblem zu der zugeordneten nichtlinearen Bewegungsgleichung. Die Nichtlinearität kommt hier aus der Kontaktbedingung zwischen der Kugel und der Plattform..

Header

Text

1+1

tmp

Die Koordinaten der Verschiebung der beiden Massen nennen wir

${\underline {Q}}(t)=\left({\begin{array}{l}u_{1}(t)\\u_{2}(t)\end{array}}\right){\text{ und damit }}\delta {\underline {Q}}=\left({\begin{array}{l}\delta u_{1}\\\delta u_{2}\end{array}}\right)$.

Für die beiden elastischen Balken verwenden die Ersatzfeder-Steifigkeit

$\displaystyle k_{1}=\displaystyle 2\cdot {\frac {12EI}{\ell ^{3}}}$.

Und als Abkürzungen verwenden wir

${\begin{array}{ll}\alpha &=\displaystyle {\frac {m_{2}}{m_{1}}}\\\kappa &=\displaystyle {\frac {k_{2}}{k_{1}}}\end{array}}$.

Declarations

Text

1+1

tmp

Die Kennlinie für den Kontakt definieren wir stückweise zu

$K(w,\epsilon )=k\cdot \left\{{\begin{array}{ccl}0&{\text{für}}&\;\;\;\;\;\;\;\;w<-\epsilon \\\displaystyle {\frac {1}{4}}\left(1+{\frac {w}{\epsilon }}\right)^{2}&{\text{für}}&-\epsilon <w<+\epsilon \\w&{\text{für}}&+\epsilon <w\\\end{array}}\right.$

Kontakt-Kennlinie

Und so sieht die Kennlinie dann aus:

Die Parabel zwischen den beiden linearen Kennlinien-Stücken macht die Kraft K stetig differentierbar in w. Das macht die numerische Integration schneller - und genauer.

Contact Characteristic

Text

1+1

tmp

Koordinaten

Die Gleichgewichtsbedingung konstruieren wir mit dem Prinzip der virtuellen Verrückungen - wir brauchen also kein Freikörperbild.

Dazu verwenden wir die Verschiebungs-Koordinaten u₁, u₂ der beiden Körper wie skizziert.

Für die allgemeine Gleichgewichtsbedingung

${\begin{array}{ll}\delta W&=0\\&=\delta W^{a}-\delta \Pi \end{array}}$

setzen wir für die virtuelle Arbeit von äußeren, eingeprägte Kräften (inklusive der D'Alembert'sche Trägheitskraft) an:

$\displaystyle \delta W^{a}=-\sum _{i=1}^{2}m_{i}\cdot \left({\ddot {u}}_{i}+g\right)\cdot \delta u_{i}$

sowie für die virtuelle Formänderungsenergie der beiden Ersatz-Federn:

$\delta \Pi =k_{1}u_{1}(t)\cdot \delta u_{1}+\underbrace {K(u_{1}(t)-u_{2}(t))} _{\displaystyle {\text{Kontakt-Kraft}}}\cdot (\delta u_{1}-\delta u_{2})$.

Sortieren nach den virtuellen Verrückungen und anschrieben der Gleichgewichtsbedingungen in Matrix-Form liefert

${\begin{pmatrix}{{m}_{1}}&0\\0&{{m}_{2}}\end{pmatrix}}\cdot {\begin{pmatrix}{{\ddot {u}}_{1}}\\{{\ddot {u}}_{2}}\end{pmatrix}}+{\begin{pmatrix}{{k}_{1}}&0\\0&0\end{pmatrix}}\cdot {\begin{pmatrix}{{u}_{1}}\\{{u}_{2}}\end{pmatrix}}=-g\cdot {\begin{pmatrix}{{m}_{1}}\\{{m}_{2}}\end{pmatrix}}+K\left(u_{1}-u_{2}\right)\cdot {\begin{pmatrix}-1\\+1\end{pmatrix}}$.

Equilibrium Conditions

Text

1+1

tmp

Mit einer passenden Bezugszeit und -länge machen wir die Bewegungsgleichungen dimensionslos.

Hier liefert die Periodendauer der Plattform - allein, ohne Kugel - die Bezugszeit T. Die statische Absenkung der Plattform - ohne Einfluss der Kugel - liefert die Bezugslänge L.

Aus der Bewegungsgleichung oben lesen wir für K(u₁-u₂) = 0 ab:

$T=\displaystyle {\frac {2\pi }{\omega _{0}}}{\text{ mit }}\omega _{0}^{2}=\displaystyle {\frac {k_{1}}{m_{1}}}$

und

$L=\displaystyle {\frac {m_{1}}{k_{1}}}g$.

Nun ist also

$\displaystyle \tau ={\frac {t}{T}}{\text{ mit }}(.)':={\frac {d}{d\tau }}(.)$

und

$\displaystyle U_{i}={\frac {u_{i}}{L}}{\text{ für }}i=1,2$

Die Kennlinie überführen wir entsprechend in

$K(w)=k_{2}\cdot L\cdot C\left(\displaystyle {\frac {W}{\epsilon }}\right){\text{ mit }}W(\tau )=U_{1}-U_{2}$.

Einsetzen liefert die zwei Differentialgleichungen

${\begin{array}{llll}\displaystyle {\frac {1}{4\cdot \pi ^{2}}}U''_{1}\left(\tau \right)&+{{U}_{1}}\left(\tau \right)&+1&=-\kappa \cdot C\left({\frac {W}{\epsilon }}\right)\\\displaystyle {\frac {\alpha }{4\cdot \pi ^{2}}}U''_{2}\left(\tau \right)&&+\alpha &=+\kappa \cdot C\left({\frac {W}{\epsilon }}\right)\end{array}}$.

Dimensionless Equations of Motion

Text

1+1

tmp

Für die numerische Lösung schreiben wir die zwei Bewegungsgleichung zweiter Ordnung um als vier Bewegungsgleichungen erster Ordnung zu

$\displaystyle {\frac {d}{d\,\tau }}\left({\begin{array}{c}U_{1}\\U_{2}\\V_{1}\\V_{2}\end{array}}\right)={\underline {f}}(U_{1},U_{2}){\text{ mit }}V_{i}={\frac {dU_{i}}{d\,\tau }}$.

Die rechte Seite der gewöhnlichen Differentialgleichung ("Right-hand-side of ODE") lautet damit

${\underline {f}}=\left({\begin{array}{c}1\\1\\\displaystyle 4\pi ^{2}(-U_{1}-1-\kappa \cdot C\cdot \left({\frac {W}{\epsilon }}\right))\\\displaystyle 4\pi ^{2}(-1+{\frac {\kappa \cdot C}{\alpha }}\cdot \left({\frac {W}{\epsilon }}\right))\end{array}}\right)$

Die dimensionslosen Parameter wählen wir zu

$\displaystyle \kappa =10,\alpha ={\frac {1}{2}},\epsilon ={\frac {1}{100}}$

und lösen die Bewegungsgleichungen numerisch - hier mit dem Runge-Kutta-Verfahren 4.ter Ordnung.

Solving

Text

1+1

tmp

Als Lösung tragen wir die Zeitverläufe der Höhen U₁, U₂ über der Zeit τ

Zeitverlauf der Lösung

und U₁, U₂  jeweils über V₁, V₂  im Phasendiagramm auf:

Phasendiagramm

Post-Processing

Text

1+1

✔ Wie genau ist die Lösung?:

Überlegen Sie, wie Sie bei diesem Ergebnis die Genauigkeit prüfen können. Das ist hier besonders einfach, weil es um ein konservatives System geht - also eins, bei dem keine Energie hinzugefügt oder dissipiert wird.

---

Links

• ...

Literature

• ...