Gelöste Aufgaben/GYRO

Aufgabenstellung

Kreisel und Ihre Bewegungsgleichungen sind immer wieder eine Herausforderung für uns Ingenieure. Hier nähern wir uns hier dem Thema mit dem Aufstellen und Lösen der Bewegungsgleichungen für große Winkel - und damit nichtlinearen Differentialgleichungen. Die sonst üblichen Linearisierungen führen wir nicht durch - und schauen, ob wir die Gleichungen noch gelöst bekommen.

Kreisel mit gelenkiger Fußpunktlagerung.

Gesucht ist die Präzessions-Bahn der eines Kreisels. Unser Kreisel ist ein Kegel der Höhe H und Radius R. Die Bewegungsgleichungen sollen mit dem Prinzip der virtuellen Verrückungen aufgestellt werden. Wir suchen nach der Trajektorie des Mittelpunktes des Kreisels für große Kippwinkel.

Lösung mit Maxima und Matlab^®

Wir arbeiten mit Maxima und Matlab.

Maxima brauchen wir zum Aufstellen der Bewegungsgleichungen, die zunächst mal recht komplex aussehen. Dabei nehmen wir die Bewegungsgleichungen wie sie aus den Gleichgewichtsbedingungen kommen und lösen sie als Anfangswertproblem.

Geomertrie des Kreisels

Header

Wir lösen hier das Anfangswertproblem mit den nichtlinearen Bewegungsgleichungen zu den Koordinaten

${\underline {\varphi }}=\left({\begin{array}{c}\varphi _{1}(t)\\\varphi _{2}(t)\\\varphi _{3}(t)\end{array}}\right)$,

wobei φ₁, φ₂ Kippwinkel bezüglich der x₁- und x₂-Achsen sind, φ₃ ist der Rotationswinkel um die Kreisel-Symmetrieachse.

Die Gleichgewichtsbedingungen kommen aus dem Prinzip der virtuellen Verrückungen. Für diese Gleichungen brauchen wir die virtuelle Arbeit der D'Alembert'schen Trägheitskräfte und die virtuelle Arbeit der Gewichtskraft des Kreisels.

/*******************************************************/
/* MAXIMA script                                       */
/* version: wxMaxima 21.05.2                           */
/* author: Andreas Baumgart                            */
/* last updated: 2022-03-27                            */
/* ref: gyroscopic precession                          */
/* description: solve using principle of virt. work    */
/*******************************************************/

Declarations

Wir brauchen zunächst die Transformationsmatrizen der Euler-Rotation, also die Drehmatrizen, die unser inertiales Koordinatensystem in das Kreiselfeste Koordinatensystem überführt. Wir kippen zunächst um die x₂-Achse, dann um die resultierende x'₁-Achse und schließlich um die Rotationsachse x₂. Die Koordinatentransformation ist dann

${\underline {\underline {D}}}_{0\to K}={\underline {\underline {D}}}(\varphi _{2}(t))\cdot {\underline {\underline {D}}}(\varphi _{1}(t))\cdot {\underline {\underline {D}}}(\varphi _{3}(t))$.

Und wir führen noch Abkürzungen für die Systemparameter

${I_{y}}={\frac {\displaystyle 3\left({{R}^{2}}+4{{H}^{2}}\right)m}{20}},{I_{z}}={\frac {\displaystyle 3{{R}^{2}}m}{10}},m={\frac {\displaystyle {\pi }H{{R}^{2}}\rho }{3}}$

und die Zustandsgrößen

${\frac {d}{dt}}{{\varphi }_{i}}(t)={{\omega }_{i}}(t)$

ein.

/* Euler-Drehmatrizen */
D[1](φ) := matrix([1,0,0],[0,cos(φ),-sin(φ)],[0,sin(φ),cos(φ)]);
D[2](φ) := matrix([cos(φ),0,sin(φ)],[0,1,0],[-sin(φ),0,cos(φ)]);
D[3](φ) := matrix([cos(φ),-sin(φ),0],[sin(φ),cos(φ),0],[0,0,1]);

coord: [[ φ[1](t), φ[2](t), φ[3](t)],
        [δφ[1]   ,δφ[2]   ,δφ[3]   ]];
q: funmake('matrix,makelist([coord[1][i]],i,1,3));

abbr: makelist(diff(φ[i](t),t)=ω[i](t),i,1,3);

params: [
  I[y] =      m       * 3/20*( R^2+4*H^2),
  I[z] =      m       * 3/10*  R^2,
   m   = ρ*%pi*R^2*H/3           ];
params: append(subst(params[3], [params[1], params[2]]),[params[3]]);
params: [params, makelist(solve(params[i], [H^3,R^4,R^2][i])[1],i,1,3)];

Equations of Motion

Die Gleichgewichtsbedingungen für den Kreisel kommen aus

Fehler beim Parsen (Unbekannte Funktion „\begin{array}“): {\displaystyle \begin{array}{lll}\delta W &=& \sum_i (\vec{F}_i\cdot\delta r_i ) - \int_V \rho \ddot{\vec{r}}_P\cdot \delta\vec{r}}_P\\ &=& 0 \end{array}} ,

wobei wir nur eine Kraft F - die Gewichtskraft $mg(-{\vec {e}}_{3})$ - haben und V das Kreisel-Volumen ist.

Der Ortsvektor zu einem Punkt des Kreisels ist

${\vec {r}}_{P}=\left({\vec {e}}_{1}{\vec {e}}_{2}{\vec {e}}_{3}\right)\cdot {\underline {\underline {D}}}_{0\to K}(\varphi _{1},\varphi _{2},\varphi _{3})\cdot {\underline {\underline {D}}}_{3}(\alpha )\cdot \left({\begin{array}{c}r\\0\\z\end{array}}\right)$

1+1

Solving

Wir brauchen

1+1

Postprocessing

Wir brauchen

1+1

---

Links

• ...

Literature

• ...