Gelöste Aufgaben/FEC1/FEC1 - Teil III: Simulieren

⬅ Zurück zu Teil II

Computerprogramm schreiben

Um unsere Bewegungsgleichungen

$\underset{\_}{\underset{\_}{M}}\underset{\_}{\ddot{Q}}+\underset{\_}{\underset{\_}{G}}\underset{\_}{\dot{Q}}+({\underset{\_}{\underset{\_}{K}}}_C\cos (2\mathrm{\Omega }t)+{\underset{\_}{\underset{\_}{K}}}_S\sin (2\mathrm{\Omega }t)+{\underset{\_}{\underset{\_}{K}}}_0)\underset{\_}{Q}=\underset{\_}{P}$

zu lösen, war der Computer und Maxima bis hierher eine gute Hilfe - wir hätten aber auch alle Gleichungen per Hand aufstellen können. Ab hier geht das nicht mehr. Für die Lösung der Bewegungsgleichungen brauchen wir Computer. Für einfache Probleme wie dieses reichen Interpreter-Programme wie Matlab©, wenn's komplizierter wird, müssten wir auf Kompiler wie C++ umsteigen. Für alle etablierten Programme gibt es Numerik-Bibliotheken, um die Probleme zu lösen, die wir jetzt angehen.

Bei allen Lösungsansätzen steht in der Bewegungsgleichung noch ein Parameter: die Dimensionslose Drehzahl η mit

$\mathrm{\Omega }=\eta \cdot {\omega }_{0,0}\text{ und der Eigenfrequenz }{\omega }_{0,0}^2=3\frac{EI}{{\ell }^3{m}_I}$

Partikulare Lösung für konstante Koeffizienten

Für die Bewegungsgleichung oben gibt es keine geschlossene partikulare Lösung - die periodischen Koeffizienten machen das unmöglich. Anschaulich berücksichtigen diese Koeffizienten die Drehwinkel-veränderliche Steifigkeit der Lager beim Drehen. 

Wir könnten von Δk = 0 ausgehen und die Lösung von

${\underset{\_}{\underset{\_}{K}}}_0{\underset{\_}{Q}}_P=\underset{\_}{P}$

berechnen. Aber hier steht nichts, was wir nicht mit einer Handrechnung auch hinbekommen würden: die Auslenkung des Systems wächst quadratisch mit der Drehzahl η(Ω) - also der Zentrifugalkraft.