Gelöste Aufgaben/FEB1

Aufgabenstellung

Auch wenn es nicht so aussieht: für das rotierende Rotorblatt suchen wir eine statische Lösung - das Problem heißt "quasistatisch".

Lageplan

Ein Hubschrauber-Rotor dreht mir der konstanten Winknelgeschwindigkeit Ω. Das Rotor-Blatt ist aus Aluminium. Gesucht ist die FEM-Lösung für der Verschiebung der Querschnitte und die Dehnung der Querschnitte.

Lösung mit Maxima

Lorem Ipsum ....

tmp

Zu diesem Problem finden wir einfach eine analytische Lösung - die nutzen wir später, um unsere Ergebnisse dimensionslos zu machen.

Ausgangspunkt ist die Gleichgewichtsbeziehung des Dehnstabes. Die Streckenlast n (hier n(r)) ist die Zentrifugalkraft - oder die D'Alembert'sche Trägheitskraft der Zentripetalbeschleunigung, also

$n(x)=-\varrho Ar{\mathrm{\Omega }}^2$ .

Die Bewegungsgleichung lautet dann

$EA{u}^{″}(r)=\varrho Ar{\mathrm{\Omega }}^2$ .

Die allgemeine Lösung ist

$AEu(r)=-\frac{\varrho A{\mathrm{\Omega }}^2{r}^3}{6}+c_{1}r+c_0$ ,

angepasst an die Randbedingungen

$\begin{array}{l}u(0)=0\\ EA\frac{d}{dr}u{|}_{r={\ell }_0}=0\end{array}$

erhalten wir

$u(r)=\frac{3{\ell }_0^{2}r-r^3}{6AE}\cdot \varrho A{\mathrm{\Omega }}^2$ .

Als Bezugslänge gewinnen wir hier

$\begin{array}{ll}{u}_s& =u({\ell }_0)\\ & =\frac{{\ell }_0^3{\mathrm{\Omega }}^2\varrho }{3E}\end{array}$===Header=== Text

1+1

tmp

Für die Komposition der Bewegungsgleichungen brauchen wir die Element-Steifigkeitsmatrix und die Element-Lastmatrix.

Wir gehen von gleichlangen Elementen aus, hier von

$I=4$

Elementen, also

${\ell }_i=\frac{{\ell }_0}{I}$.

Wie im Abschnitt "Finite Elemente Methode" verwenden wir die linearen Ansatzfunktionen

${\varphi }_1=1-\xi \text{ und }{\varphi }_1=\xi$,

also

$u_i({\xi }_i)=U_{i-1}\cdot {\varphi }_1({\xi }_1)+U_i\cdot {\varphi }_2({\xi }_i)$.

Die Element-Steifigkeitsmatrix ergibt sich für den Dehnstab zu

${\underset{\_}{\underset{\_}{K}}}_i=\frac{EA}{{\ell }_i}\cdot \left(\begin{array}{rr}1& -1\\ -1& 1\end{array}\right)$.

Die Element-Lastmatrix ist etwas schwieriger. Für ein einzelnes Element i ist

$\delta W^a=\underset{r_i-1}{\overset{r_i}{\int }}\varrho Ar{\mathrm{\Omega }}^2\cdot \delta u(r)dr$.

Mit

$r=r_{i-1}+{\ell }_i{\xi }_i\text{ und }{r}_{i-1}=(i-1)\cdot {\ell }_i$

schreiben die Beziehung um zu

$\delta W^a=\varrho A{\mathrm{\Omega }}^2{\ell }_i^2\cdot \underset{0}{\overset{1}{\int }}(i+\xi )\cdot \delta u(\xi )d\xi$,

wir erhalten

${\underset{\_}{P}}_i=\varrho A{\mathrm{\Omega }}^2{\ell }_i^2\left(\begin{array}{c}\frac{3i+1}{6}\\ \frac{3i+2}{6}\end{array}\right)$. ===Declarations=== Text

1+1

tmp

Die Gleichgewichtsbeziehungen des Gesamtsystems erhalten wir durch aus der Addition aller virtuellen Arbeiten des Systems, praktisch durch das Hinzuaddieren der Anteile je Element in die System-Matrizen für K und P:

$\underset{\_}{\underset{\_}{K}}=\left(\begin{array}{ccccc}\frac{AE}{{\ell }_i}& -\frac{AE}{{\ell }_i}& 0& 0& 0\\ -\frac{AE}{{\ell }_i}& \frac{2AE}{{\ell }_i}& -\frac{AE}{{\ell }_i}& 0& 0\\ 0& -\frac{AE}{{\ell }_i}& \frac{2AE}{{\ell }_i}& -\frac{AE}{{\ell }_i}& 0\\ 0& 0& -\frac{AE}{{\ell }_i}& \frac{2AE}{{\ell }_i}& -\frac{AE}{{\ell }_i}\\ 0& 0& 0& -\frac{AE}{{\ell }_i}& \frac{AE}{{\ell }_i}\end{array}\right)$

und

$\underset{\_}{P}=\varrho A{\ell }_i^2{\mathrm{\Omega }}^2\cdot \left(\begin{array}{c}\frac{1}{6}\\ 1\\ 2\\ 3\\ \frac{11}{6}\end{array}\right)$===Equlibrium Conditions=== Text

1+1

tmp

Die geometrische Randbedingung U₀ = 0 arbeiten wir ein, indem wir die erste Zeile des Gleichungssystem streichen sowie die erste Spalte der Steifigketsmatrix.===Boundary Conditions=== Text

1+1

tmp

Die Lösung des linearen Gleichungssystems

$\underset{\_}{\underset{\_}{K}}\cdot \underset{\_}{U}=\underset{\_}{P}$

ist

$\underset{\_}{U}=\left(\begin{array}{c}\frac{47l_0^3{\mathrm{\Omega }}^2\rho }{384E}\\ \frac{11l_0^3{\mathrm{\Omega }}^2\rho }{48E}\\ \frac{39l_0^3{\mathrm{\Omega }}^2\rho }{128E}\\ \frac{l_0^3{\mathrm{\Omega }}^2\rho }{3E}\end{array}\right)$.

Oder - in dimensionsloser Form

$\frac{\underset{\_}{U}}{u_s}=\left(\begin{array}{c}\frac{47}{128}\\ \frac{11}{16}\\ \frac{117}{128}\\ 1\end{array}\right)$.===Solving=== Text

1+1

tmp

Auslenkung u(x).

Die Ergebnisse der FE-Rechnung und der analytischen Lösung können wir jetzt übereinander auftragen:

Die Dehnungen (und Spannungen) im Bauteil sind

${\epsilon }_{rr}=\frac{du}{dr}$.

Auch dieses Ergebnis können wir auftragen:

Dehnung ε(x).

===Post-Processing===

✔ Konstante Dehnung je Element:

Was ausschaut wie ein Fehler - nämlich die "Treppenfunktion" für die Dehnung im FE-Modell - ist in Wirklichkeit die Folge unserer linearen Ansatzfunktionen. Diese abgeleitet liefern eine konstante Dehnung je Element.

1+1

Links

• ...

Literature

• ...