Gelöste Aufgaben/W8Zv

Aufgabenstellung

Zu den tabellierten Standardlösungen für den Euler-Bernoulli-Blaken berechnen wir eine Näherungslösung für einen beidseitig gelenkig gelagerten Euler-Bernoulli-Balken:

Caption

Hier soll mit dem Ansatz mit der Methode der Finiten Elemente gearbeitet werden. Gesucht ist das Verschiebungsfeld w(x) im Vergleich von FEM und analytischer Lösung.

Diese Aufgabe ist eng verwandt mit W8Zt und W8Zu.

Lösung mit Maxima

Mit dem Föppl-Symbol "<>", sowie

$\alpha =a/\ell$, $\beta =1-\alpha$ und $\xi =x/\ell$

ist die analytische Lösung:

$EIw(x)=\frac{F{\ell }^3}{6}[\beta \xi (1-{\beta }^2-{\xi }^2)+<\xi -\alpha {>}^3]$.

Bei dieser Lösung hat die unabhängige Koordinate x ihren Ursprung in A - wir verwenden unten einen anderen Ursprung!

Mit den passenden Ansatzfunktionen berechnen wir eine Näherungslösung des Problems nach dem Prinzip der virtuellen Verrückungen.

tmp

Bei der Methode der Finiten Elemente setzen wir die virtuelle Formänderungsenergie des Systems additiv aus den Anteilen je Element zusammen. Hier arbeiten wir mit zwei Elementen, die am Kraft-Angriffspunkt aneinander stoßen.

Also ist

$\delta \mathrm{\Pi }=\delta {\mathrm{\Pi }}_1+\delta {\mathrm{\Pi }}_2$ ,

die virtuelle Arbeit der äußeren Kraft F ist

$\delta W^a=F\cdot \delta w_F$ .

Header

Text

1+1

tmp

Declarations

Text

1+1

tmp

Aus der Finite Elemente Methode kennen wir die virtuelle Formänderungsenergie eines Balkenelements

$\delta {\mathrm{\Pi }}_i=(\delta W_{i-1},\delta {\mathrm{\Phi }}_{i-1},\delta W_i,\delta {\mathrm{\Phi }}_i)\cdot {\underset{\_}{\underset{\_}{K}}}_i\cdot \left(\begin{array}{l}{W}_{i-1}\\ {\mathrm{\Phi }}_{i-1}\\ W_i\\ {\mathrm{\Phi }}_i\end{array}\right)$

für das klassische x-z-Koordinatensystem des Euler-Bernoulli-Balkens und mit der Element-Steifigkeitsmatrix

${\underset{\_}{\underset{\_}{K}}}_i=\frac{{E}{I}}{{\ell }_E^3}\cdot \left(\begin{array}{cccc}12& 6\cdot {\ell }_E& -12& 6\cdot {\ell }_E\\ 6\cdot {\ell }_E& 4\cdot {\ell }_E^2& -6\cdot {\ell }_E& 2\cdot {\ell }_E^2\\ -12& -6\cdot {\ell }_E& 12& -6\cdot {\ell }_E\\ 6\cdot {\ell }_E& 2\cdot {\ell }_E^2& -6\cdot {\ell }_E& 4\cdot {\ell }_E^2\end{array}\right)$.

Virtual Strain-Energy per Element

Text

1+1

tmp

Das Gesamt-Gleichungssystem für die Koordinaten

$\underset{\_}{Q}=\left(\begin{array}{l}{W}_0\\ {\mathrm{\Phi }}_0\\ W_1\\ {\mathrm{\Phi }}_1\\ W_2\\ {\mathrm{\Phi }}_2\end{array}\right)$

ist

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

mit der 6x6 Matrix K₀. In diese Matrix müssen wir jetzt die Element-Steifigkeitsmatrizen für das Element 1 und das Element 2 hineinaddieren, also

${\underset{\_}{\underset{\_}{K}}}_0=EI\cdot \left(\begin{array}{cccccc}\frac{12}{{\ell }_1^3}& \frac{6}{{\ell }_1^2}& -\frac{12}{{\ell }_1^3}& \frac{6}{{\ell }_1^2}& 0& 0\\ \frac{6}{{\ell }_1^2}& \frac{4}{{\ell }_1}& -\frac{6}{{\ell }_1^2}& \frac{2}{{\ell }_1}& 0& 0\\ -\frac{12}{{\ell }_1^3}& -\frac{6}{{\ell }_1^2}& \frac{12}{{\ell }_2^3}+\frac{12}{{\ell }_1^3}& \frac{6}{{\ell }_2^2}-\frac{6}{{\ell }_1^2}& -\frac{12}{{\ell }_2^3}& \frac{6}{{\ell }_2^2}\\ \frac{6}{{\ell }_1^2}& \frac{2}{{\ell }_1}& \frac{6}{{\ell }_2^2}-\frac{6}{{\ell }_1^2}& \frac{4}{{\ell }_2}+\frac{4}{{\ell }_1}& -\frac{6}{{\ell }_2^2}& \frac{2}{{\ell }_2}\\ 0& 0& -\frac{12}{{\ell }_2^3}& -\frac{6}{{\ell }_2^2}& \frac{12}{{\ell }_2^3}& -\frac{6}{{\ell }_2^2}\\ 0& 0& \frac{6}{{\ell }_2^2}& \frac{2}{{\ell }_2}& -\frac{6}{{\ell }_2^2}& \frac{4}{{\ell }_2}\end{array}\right)$

Weil die Auslenkung des Kraft-Angriffspunktes

$\delta w_F=\delta W_1$

ist

$\underset{\_}{P}=\left(\begin{array}{c}0\\ 0\\ F\\ 0\\ 0\\ 0\end{array}\right)$.

Hier fehlen noch die ...

Equilibrium Conditions

Text

1+1

tmp

Einarbeiten der Randbedingungen in die Systemmatrix

Im Gesamt-Gleichungssystem müssen wir noch die Randbedingungen einarbeiten, nämlich

$\begin{array}{ll}{W}_0=0& (\delta W_0=0)\\ W_2=0& (\delta W_2=0)\end{array}$

Das machen wir durch Streichen der zugehörigen Zeilen und Spalten (1 und 5) im Gleichungssystem.

Boundary Conditions

Text

1+1

tmp

Die Lösung des verbleibenden Gleichungssystems mit

$\alpha =\frac{a}{{\ell }_0},{\ell }_0=a+b$

ist

$\begin{array}{ccl}{\mathrm{\Phi }}_0\cdot & \frac{EI}{{\ell }_0^2\cdot F}& =\frac{2\cdot \alpha -3\cdot {\alpha }^2+{\alpha }^3}{6}\\ W_1\cdot & \frac{EI}{{\ell }_0^3\cdot F}& =\frac{{\alpha }^2-2\cdot {\alpha }^3+{\alpha }^4}{3}\\ {\mathrm{\Phi }}_1\cdot & \frac{EI}{{\ell }_0^2\cdot F}& =\frac{\alpha -3\cdot {\alpha }^2+2\cdot {\alpha }^3}{3}\\ {\mathrm{\Phi }}_2\cdot & \frac{EI}{{\ell }_0^2\cdot F}& =\frac{{\alpha }^3-\alpha }{6}\end{array}$

Und die Ergebnisse können wir mit der analytischen Lösung vergleichen:

Solve

Text

1+1

tmp

Die Auslenkung des Kraft-Angriffspunktes w(a) tragen wir hier auf:

Absenkung w(a) des Kraftangriffspunktes

Post-Process: Results

Text

1+1

tmp

Der Vergleich der tabellierten Lösungen mit unseren Lösungen:

⚠ Ist das jetzt die analytische Lösung?:

Sind die FEM-Ergebnisse nicht exakt die der analytischen Lösung? Stimmt das? Und wenn ja:  wie kommt das?

Post-Process: Compare with Analytic Solution

Text

1+1

---

Links

• W8Zt, W8Zu

Literature

• ...

Lageplan