Gelöste Aufgaben/FEAD

Aufgabenstellung

Statt Formfunktionen über die ganze Stablänge anzusetzten wie in FEAB gehen wir jetzt nach der Methode der Finiten Elemente vor.

Caption

Lageplan

Gesucht ist die Näherungslösung für die statische Auslenkung der Stab-Querschnitte mit dem Prinzip der Virtuellen Verrückungen und dem Ansatz nach der Methode der Finiten Elemente. Dazu wir unterteilen die Struktur in Elemente und setzen Lineare Trialfunctions für die Verformung in den Elemeten an.

Lösung mit Matlab

Lorem Ipsum ....

tmp

Die Lösung basiert auf Maxima 16.04.2 - wir interessieren uns also vor allem für die Struktur der Lösung.===Header=== Text

📑 toggle code

1+1

tmp

Für Maxima brauchen wir einige Deklarationen.

Die Anzahl der Finiten-Elementen I wählen wir zu

$I=3$

Wir wählen Elemente gleicher Länge, also ist die Elementlänge l_i im Element i

$\begin{array}{ll}{\ell }_i& ={\ell }_e\\ & ={\ell }_0/I\end{array}$.===Declarations=== Text

📑 toggle code

1+1

tmp

Wir können die Trial-Functions als

$u(x)=\sum _{i=0}^I{U}_i\cdot {\varphi }_i(x)$

mit

${\varphi }_i(x)=\{\begin{array}{cl}{\xi }_{i-1}& \text{ wobei }x={\ell }_e\cdot ((i-1)+{\xi }_{i-1}),0<{\xi }_{i-1}<1\\ (1-{\xi }_i)& \text{ wobei }x={\ell }_e\cdot (i+{\xi }_i),0<{\xi }_i<1\end{array}$.

Für i=2 sieht die Trial-Function dann so aus:

Trial-Function für i=2.

plot2d([[parametric, 1+t, t, [t,0,1]],[parametric, 2+t, 1-t, [t,0,1]]],

             [x,0,3], [y,-0.2,1.2],

             [xlabel, "x/l ->"], [ylabel, "ϕ2/1"])

Für die Methode der Finiten Elemente ist allerdings die Sichtweise je Elemement anschaulicher, also

$u_i(x_i)=\{\begin{array}{c}{U}_{i-1}\cdot (1-\left(\frac{x_i}{{\ell }_i}\right))+U_i\cdot \left(\frac{x_i}{{\ell }_i}\right)\\ 0\text{ sonst }\end{array}$,

Und für Element i=2 wird der Verlauf der Querschnitts-Verschiebung also durch

$u_2(x_2)=\{\begin{array}{c}{U}_1\cdot (1-\left(\frac{x_2}{{\ell }_2}\right))+U_2\cdot \left(\frac{x_2}{{\ell }_2}\right)\\ 0\text{ sonst }\end{array}$.

entsprechend

$\delta u_2(x_2)=\{\begin{array}{c}\delta U_1\cdot (1-\left(\frac{x_2}{{\ell }_2}\right))+\delta U_2\cdot \left(\frac{x_2}{{\ell }_2}\right)\\ 0\text{ sonst }\end{array}$

beschreiben.

Bei 3 Finiten Elementen sind die System-Koordinaten dann

$\underset{\_}{Q}=\left(\begin{array}{l}{U}_0\\ U_1\\ U_2\\ U_3\end{array}\right)$===Formfunctions=== Text

📑 toggle code

1+1

tmp

Die gesamte virtuelle Arbeit am System können wir nun als Summe der virtuellen Einzelarbeiten je Element hinschreiben, also

$\begin{array}{ll}\delta W& =\delta W^a-\delta \mathrm{\Pi }\\ & =\sum _{i=1}^I(\delta W_i^a-\delta {\mathrm{\Pi }}_i)\end{array}$.

Für jedes Element erhalten wir die virtuelle Formänderungsenergie

$\begin{array}{ll}\delta {\mathrm{\Pi }}_i& =\underset{0}{\overset{{\ell }_i}{\int }}EA\cdot {u_i}^{\prime }\cdot \delta {u_i}^{\prime }d{x}_i\\ {\displaystyle }& =(\delta U_{i-1},\delta U_i)\cdot \frac{EA}{{\ell }_i}\cdot \left(\begin{array}{rr}1& -1\\ -1& 1\end{array}\right)\cdot \left(\begin{array}{l}{U}_{i-1}\\ U_i\end{array}\right)\end{array}$

und die virtuelle Arbeit der äußeren, eingeprägten Gewichtskraft

$\begin{array}{ll}\delta W_i^a& =\underset{0}{\overset{{\ell }_i}{\int }}\varrho \cdot A\cdot g\cdot \delta u_{i}d{x}_i\\ {\displaystyle }& =\frac{1}{2}\cdot (\delta U_{i-1},\delta U_i)\cdot \varrho \cdot A\cdot g\cdot {\ell }_i\cdot \left(\begin{array}{l}1\\ 1\end{array}\right)\end{array}$.

Element-wise contributions to δW

Text

📑 toggle code

1+1

tmp

Für jedes Element müssen wir nun die virtuellen Arbeiten zum Gesamt-Gleichungssystem zusammenaddieren.

Die Gleichgewichtsbedingung

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

liefern ein Gleichungssystem, in das wir nun die kinematische Randbedingung U₀=0, δU₀=0 durch Streichen der ersten Zeile des Gleichungssystems und der ersten Spalte der Gesamt-Steifigkeitsmatrix einarbeiten.

Wir wählen nun noch für jedes Element die gleiche Element-Länge ℓ_i = ℓ_e, mit ℓ_e = ℓ₀/3 und es bleibt

$\frac{EA}{{\ell }_i}\left(\begin{array}{ccc}2& -1& 0\\ -1& 2& -1\\ 0& -1& 1\end{array}\right)\cdot \left(\begin{array}{l}{U}_1\\ U_2\\ U_3\end{array}\right)=\varrho A{\ell }_{i}g\cdot \left(\begin{array}{l}1\\ 1\\ \frac{1}{2}\end{array}\right)$.===Equlilibrium Conditions=== Text

📑 toggle code

1+1

tmp

Die Gleichgewichtsbedingungen lösen wir nun und erhalten

$U_0=0,U_1=\frac{5g{\ell }_i^2\varrho }{2E},U_2=\frac{4g{\ell }_i^2\varrho }{E},U_3=\frac{9g{\ell }_i^2\varrho }{2E}$===Solving=== Text

📑 toggle code

1+1

tmp

Wir tragen die Ergebnisse für die numerische Näherungslösung gegen die exakte Lösung auf. Dabei setzen wir l_i=l₀/3. Die Verschiebungen sind elementweise für die Elemente e = 1,2,3 aufgetragen:

Verschiebung u(x) - analytische und FE-Lösung

Post-processing

Text

📑 toggle code

1+1

{{Template:MyTipp|title=Spannungen im Stab|text=Tragen Sie auch die Spannungen im Stab über die Stablänge an! Berechnen Sie die Spannungen auf Basis der Dehnung 

$\epsilon =\frac{du}{dx}.$

Links

• ...

Literature

• ...