Gelöste Aufgaben/FEC1/FEC1 - Teil II: Das Modell erstellen

Mathematisches Modell formulieren

Die Bewegungsgleichungen leiten wir mit dem Prinzip der virtuellen Verrückungen her, die Gleichgewichtsbedingung lautet also

$\begin{array}{lcl}\delta W& =& \delta W^a-\delta \mathrm{\Pi }\\ & \stackrel{!}{=}& 0\end{array}$.

Die Anteile ergeben sich jeweils aus der Summe der Teilsysteme

$\delta W^a=\sum _i{W}_i^a,\delta \mathrm{\Pi }=\sum _i{\mathrm{\Pi }}_i$

mit

$i\in \{S,I,T,B_B,B_C,U_A,U_D\}$ .

Auf keines der Teilsysteme wirkt eine äußere, eingeprägte Last (keine Kraft, kein Moment) - es treten also nur die Anteile der d'Alembert'schen Kräfte und der Formänderungsenergie auf.

Impeller- und Turbinen-Rad sowie die Unwuchten erfassen wir als Starrkörper - ihre virtuelle Formänderungsenergie ist Null:

$\begin{array}{lc}\delta {\mathrm{\Pi }}_I& =0\\ \delta {\mathrm{\Pi }}_T& =0\\ \delta {\mathrm{\Pi }}_{U,A}& =0\\ \delta {\mathrm{\Pi }}_{U,D}& =0\end{array}$ .

Für die beiden Lager berücksichtigen wir keine Massen, hier ist also

$\begin{array}{lc}\delta W_{B,B}& =0\\ \delta W_{B,C}& =0\end{array}$ .

Die verbleibenden virtuellen Arbeiten des Systems schreiben wir an, sortieren nach den Koordinaten und deren Variation und formen daraus das System von Bewegungsgleichungen.

tmp

Die Welle ist ein elastisches Kontinuum. Wir können sie als Stab (als eindimensionales Kontinuum) modellieren. Dabei passen

• ein Euler-Bernoulli-Balken oder
• ein Timoshenko-Balken.

Der Einfachheit halber wählen wir den klassischen Euler-Bernoulli-Balken. Die Ansatzfunktionen und Element-Steifigkeits-Matrizen nehmen wir aus der FEM-Formulierung für den Euler-Bernoulli-Balken.

Für die Komplexität des Modells ist die Anzahl der Finiten Elemente entscheidend. Jeder Knoten eines Finite Elements für diesen 3D-Euler-Bernoulli-Ballken hat die vier Koordinaten

$\begin{array}{c}{V}_i(t)\\ {\mathrm{\Psi }}_i(t)\\ W_i(t)\\ {\mathrm{\Phi }}_i(t)\end{array}$.

Die Anzahl der Koordinaten K ist also mit N als Anzahl der Elemente

$K=4\cdot (N+1)$

Uns geht es hier primär um die Anschaulichkeit - wir approximieren die Bewegung der Welle mit einem Finiten Element.

Damit sind die Koordinaten der Welle - und damit des Gesamt-Systems

Hier steht der Index 0 für Punkt A und Index 1 für Punkt B.===Diskretisierung der Welle=== Text

1+1

tmp

Allgemein ist die Virtuelle Arbeit von D'Alembert'schen Trägheitkräften definiert als

$\delta W^a=\underset{V}{\int }(-\rho \ddot{\overrightarrow{r}})\cdot \delta \overrightarrow{r}dV$.

Wir brauchen also nur einen Ortsvektor zu den materiellen Punkten der Körper als Funktion der Koordinaten der Bewegung - den wir dann in die Beziehung oben einsetzen und über das Volumen des Körpers integrieren.

Virtuelle Arbeit von D'Alembert'schen Trägheitkräften

Text

1+1

tmp

Virtuelle Formänderungsenergie

Text

1+1

tmp

Komponieren der System-Matrizen

{{{text}}}

1+1

Links

• ...

Literature

• ...