Gelöste Aufgaben/StaF: Unterschied zwischen den Versionen

Aus numpedia
Zur Navigation springen Zur Suche springen
Keine Bearbeitungszusammenfassung
Keine Bearbeitungszusammenfassung
 
(5 dazwischenliegende Versionen desselben Benutzers werden nicht angezeigt)
Zeile 17: Zeile 17:


<onlyinclude>
<onlyinclude>
[[Datei:Screenshot 20210111-063733~2.png|100px|left|mini|Caption]]
[[Datei:StaB-01.PNG|150px|left|mini|Stabwerk mit in den Knoten fest verbundenen Stäben.]]
Gesucht ist ein Vergleich zwischen der klassischen Stabwerkstheorie und einer Herangehensweise, bei der wir eine feste Verbindung der Stäbe in den Knoten ansetzten. Grundlage des Modells ist die FEM-Lösung der Felddifferentialgleichung im Vergleich zur Lösung in Problemstellung „Stab“.
Gesucht ist ein Vergleich zwischen der klassischen Stabwerkstheorie und einer Herangehensweise, bei der wir eine feste Verbindung der Stäbe in den Knoten ansetzten. Grundlage des Modells ist die FEM-Lösung der Felddifferentialgleichung im Vergleich zur Lösung in Problemstellung „Stab“.
</onlyinclude>
</onlyinclude>
Zeile 40: Zeile 40:
Mit den Konventionen für die Knoten-Verschiebungen aus [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Gel%C3%B6ste_Aufgaben/Stab Stab] ist  
Mit den Konventionen für die Knoten-Verschiebungen aus [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Gel%C3%B6ste_Aufgaben/Stab Stab] ist  
::<math>
::<math>
\delta W^a = -\delta W_{4,0} \cdot F
\delta W^a = +\delta W_{4,0} \cdot F
</math>.
</math>.
Für <math>\delta \Pi</math> gilt  
Für <math>\delta \Pi</math> gilt  
Zeile 124: Zeile 124:
::<math>
::<math>
\left(\begin{array}{c}U_{I,0}\\W_{I,0}\\\Phi_{I,0}\end{array}\right) =  
\left(\begin{array}{c}U_{I,0}\\W_{I,0}\\\Phi_{I,0}\end{array}\right) =  
\underline{\underline{D}}_R(\alpha)\cdot
\underline{\underline{D}}_R(\alpha_k)\cdot
\left(\begin{array}{c}U_{I,k}\\W_{I,k}\\\Phi_{I,k}\end{array}\right)
\left(\begin{array}{c}U_{I,k}\\W_{I,k}\\\Phi_{I,k}\end{array}\right)
</math>
</math>
mit der Transformationsmatrix
mit der Transformationsmatrix
::<math>
::<math>
\underline{\underline{D}}_R(\alpha) =  
\underline{\underline{D}}_R(\alpha_k) =  
\left(\begin{array}{c} \cos(\alpha)& \sin(\alpha)& 0\\
\left(\begin{array}{c} \cos(\alpha_k)& \sin(\alpha_k)& 0\\
                       -\sin(\alpha)& \cos(\alpha)& 0\\
                       -\sin(\alpha_k)& \cos(\alpha_k)& 0\\
                             0    &        0  & 1\end{array}\right)
                             0    &        0  & 1\end{array}\right).
</math>
</math>
Damit wir für die Elementsteifigkeitsmatrix - mit beiden Anfangs- und Endknoten des Elements - vom "0"-System ins "k"-System transformieren, brauchen wir die neue Matrix
Es ist praktisch, an dieser Stelle die Abkürzung
::<math>
::<math>
\underline{\underline{D}}_T(\alpha) =  
\underline{q}_{I,0} = \left(\begin{array}{c}U_{I,0}\\W_{I,0}\\\Phi_{I,0}\end{array}\right)
\left(\begin{array}{c} \underline{\underline{D}}^T_R(\alpha)& \underline{\underline{0}}\\
</math>
                       \underline{\underline{0}}            & \underline{\underline{D}}^T_R(\alpha)
für die Koordinaten eines Knoten im Referenzsystem einzuführen. Also ist
::<math>
\left(\begin{array}{c}U_{I,k}\\W_{I,k}\\\Phi_{I,k}\end{array}\right) = \underline{\underline{D}}_R^T(\alpha_k)\cdot \underline{q}_{I,0} .
</math>
 
Damit wir für die Elementsteifigkeitsmatrix - mit beiden Anfangs- und Endknoten des Elements - vom "0"-System ins "k"-System transformieren, brauchen wir die neue Transformations-Matrix
::<math>
\underline{\underline{T}}(\alpha_k) :=  
\left(\begin{array}{cc} \underline{\underline{D}}_R^T(\alpha_k)& \underline{\underline{0}}\\
                       \underline{\underline{0}}            & \underline{\underline{D}}_R^T(\alpha_k)
\end{array}\right)
\end{array}\right)
</math>
</math>
Mit diesen Transformationsmatrizen ist
Mit diesen ist
::<math>
::<math>
\delta\Pi = \sum_{k=1}^4 \delta \underline{Q}_{k,0}^T\cdot  
\delta\Pi = \sum_{k=1}^4 \delta \underline{Q}_{k,0}^T\cdot  
                         \underbrace{
                         \underbrace{
                         \underline{\underline{D}}_T^T(\alpha_k)\cdot  
                         \underline{\underline{T}}^T(\alpha_k)\cdot  
                         \underline{\underline{K}}_{k,k}\cdot  
                         \underline{\underline{K}}_{k,k}\cdot  
                         \underline{\underline{D}}_T(\alpha_k)}_{:=
                         \underline{\underline{T}}(\alpha_k)}_{:=
                         \underline{\underline{K}}_{k,0} }\cdot
                         \underline{\underline{K}}_{k,0} }\cdot
                         \underline{Q}_{k,0}
                         \underline{Q}_{k,0}
</math>
</math>
Wir sammeln alle Koordianten der Knoten im ''0''-System in  
Die resultierenden Element-Steifigkeitsmatrizen sind im folgenden aufgeschreiben:
<math>\underline{Q}_{0}</math> und schreiben die Gleichgewichtsbedingungen in der Form
 
<table class="wikitable mw-collapsible mw-collapsed" style="background-color:white; float: none; margin-right:14px;">
<tr><th>Element-Steigigkeitsmatrizen mit globalen Koordinaten</th></tr>
<tr><th>Element #1</th></tr>
<tr><td>
<math>
\underline{\underline{K}}_{1,0} = \begin{pmatrix}\frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}} & -\left( \frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\\
0 & \frac{2 A E}{{\ell_0}} & 0 & 0 & -\left( \frac{2 A E}{{\ell_0}}\right)  & 0\\
\frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}} & 0 & \frac{2 {{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\\
-\left( \frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right)  & \frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right) \\
0 & -\left( \frac{2 A E}{{\ell_0}}\right)  & 0 & 0 & \frac{2 A E}{{\ell_0}} & 0\\
\frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}} & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right)  & 0 & \frac{2 {{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
<tr><th>Element #2</th></tr>
<tr><td>
<math>
\underline{\underline{K}}_{2,0} = \begin{pmatrix}\frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
\frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
-\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
-\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
<tr><th>Element #3</th></tr>
<tr><td>
<math>
\underline{\underline{K}}_{3,0} = \begin{pmatrix}\frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
-\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
-\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
\frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
<tr><th>Element #4</th></tr>
<tr><td>
<math>
\underline{\underline{K}}_{4,0} = \begin{pmatrix}\frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}} & -\left( \frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\\
0 & \frac{A E}{{\ell_0}} & 0 & 0 & -\left( \frac{A E}{{\ell_0}}\right)  & 0\\
\frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}} & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
-\left( \frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right) \\
0 & -\left( \frac{A E}{{\ell_0}}\right)  & 0 & 0 & \frac{A E}{{\ell_0}} & 0\\
\frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}} & 0 & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
</table>
 
Wir sammeln nun alle Koordianten der Knoten im ''0''-System in  
::<math>\underline{Q}_{0}^T = \left(\underline{q}_{1,0}^T, \underline{q}_{2,0}^T, \underline{q}_{3,0}^T, \underline{q}_{4,0}^T \right)</math>
und schreiben die Gleichgewichtsbedingungen in der Form
::<math>
::<math>
\delta\underline{Q}_{0}^T\cdot\left(\underline{\underline{K}}_{0}\cdot\underline{Q}_{0}  
\delta W = \delta\underline{Q}_{0}^T\cdot\left(\underline{\underline{K}}_{0}\cdot\underline{Q}_{0}  
- \underline{P}\right) = 0
- \underline{P}\right) \stackrel{!}{=} 0
</math>
</math>
an.
Dabei kommen die Beiträge zur Gesamt-Steifigkeitsmatrix (hier noch in der Fassung ohne Berücksichtigung der Randbedingungen) aus den vier Beiträgen der virtuellen Formänderungsenergie - die wir hier farblich gekennzeichnet haben:
::[[Datei:StaF-11.png|350px|none|Einsortieren der Anteile der virtuellen Formänderungsenergie nach den Knoten.]]


an. Die Randbedingungen arbeiten wir hier durch das Streichen der passenden Zeilen für <math>\delta U_{1,0},\delta W_{1,0},\delta U_{2,0},\delta W_{2,0}</math> sowie der passenden Spalten für <math>U_{1,0},W_{1,0},U_{2,0},W_{2,0}</math> ein.
===Einarbeitung der Randbedingungen===
Die Randbedingungen arbeiten wir hier durch das Streichen der passenden Zeilen für <math>\delta U_{1,0},\delta W_{1,0},\delta U_{2,0},\delta W_{2,0}</math> sowie der passenden Spalten für <math>U_{1,0},W_{1,0},U_{2,0},W_{2,0}</math> ein.


Das resultierende Gleichungssystem ist dies:
Das resultierende Gleichungssystem ist dies:
Zeile 195: Zeile 260:
\end{pmatrix}
\end{pmatrix}
</math>
</math>
Hier kommt jetzt irgendein Text.


===Solving===
Das Ergebnis ist - für die gleichen Parameter wie in [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Gel%C3%B6ste_Aufgaben/Stab Stab] -
::<math>
::<math>
Some Text
\begin{array}{l}
U_{1,0}=0\\
W_{1,0}=0\\
\Phi_{1,0}=3.97 \frac{F}{E a^2}\\
U_{2,0}=0\\
W_{2,0}=0\\
\Phi_{2,0}=2.45 \frac{F}{E a^2}\\
U_{3,0}=57.7 \frac{F}{E a}\\
W_{3,0}=1.67\cdot 10^2 \frac{F}{E a}\\
\Phi_{3,0}=2.06 \frac{F}{E a^2}\\
U_{4,0}=-57.7 \frac{F}{E a}\\
W_{4,0}= 3.67\cdot 10^2 \frac{F}{E a}\\
\Phi_{4,0}=3.12 \frac{F}{E a^2}
\end{array}
</math>
</math>
Und damit sind auch die Verläufe der Schnittlasten in den Stäben identisch - wir brauchen also die Ergebnisse nicht noch einmal auftragen.


<!-------------------------------------------------------------------------------->
<!-------------------------------------------------------------------------------->
Zeile 209: Zeile 289:
</syntaxhighlight>
</syntaxhighlight>
}}
}}
<table class="wikitable mw-collapsible mw-collapsed" style="background-color:white; float: none; margin-right:14px;">
<tr><th>Element-Steigigkeitsmatrizen mit globalen Koordinaten</th></tr>
<tr><th>Element #1</th></tr>
<tr><td>
<math>
{k_1} = \begin{pmatrix}\frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}} & -\left( \frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\\
0 & \frac{2 A E}{{\ell_0}} & 0 & 0 & -\left( \frac{2 A E}{{\ell_0}}\right)  & 0\\
\frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}} & 0 & \frac{2 {{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\\
-\left( \frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right)  & \frac{8 {{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right) \\
0 & -\left( \frac{2 A E}{{\ell_0}}\right)  & 0 & 0 & \frac{2 A E}{{\ell_0}} & 0\\
\frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}} & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{2 {{A}^{2}} E \eta }{{\ell_{0}^{2}}}\right)  & 0 & \frac{2 {{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
<tr><th>Element #2</th></tr>
<tr><td>
<math>
{k_2} = \begin{pmatrix}\frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
\frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
-\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
-\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
<tr><th>Element #3</th></tr>
<tr><td>
<math>
{k_3} = \begin{pmatrix}\frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
-\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
-\left( \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta +3 {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right) \\
\frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & -\left( \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}}\right)  & \frac{3 {{A}^{2}} E \eta +{\ell_{0}^{2}} A E}{4 {\ell_{0}^{3}}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\\
\frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{4 {\ell_{0}^{2}}}\right)  & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {\ell_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
<tr><th>Element #4</th></tr>
<tr><td>
<math>
{k_4} = \begin{pmatrix}\frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}} & -\left( \frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\\
0 & \frac{A E}{{\ell_0}} & 0 & 0 & -\left( \frac{A E}{{\ell_0}}\right)  & 0\\
\frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}} & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
-\left( \frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}}\right)  & 0 & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right)  & \frac{{{A}^{2}} E \eta }{{\ell_{0}^{3}}} & 0 & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right) \\
0 & -\left( \frac{A E}{{\ell_0}}\right)  & 0 & 0 & \frac{A E}{{\ell_0}} & 0\\
\frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}} & 0 & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{{{A}^{2}} E \eta }{2 {\ell_{0}^{2}}}\right)  & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
</math>
</td></tr>
</table>


<hr/>
<hr/>

Aktuelle Version vom 16. November 2024, 15:39 Uhr


Aufgabenstellung

Wir untersuchen die Belastung eines ebenen Stabwerks. Die Stäbe haben wie skizziert die Länge ℓ bzw. ℓ/2. Die Struktur wird mit der Kraft F belastet.


Stabwerk mit in den Knoten fest verbundenen Stäben.

Gesucht ist ein Vergleich zwischen der klassischen Stabwerkstheorie und einer Herangehensweise, bei der wir eine feste Verbindung der Stäbe in den Knoten ansetzten. Grundlage des Modells ist die FEM-Lösung der Felddifferentialgleichung im Vergleich zur Lösung in Problemstellung „Stab“.

Wir stellen das Modell des Stabwerks mit dem Prinzip der virtuellen Verrückungen auf und vergleichen, wie sich diese von der Herangehensweise aus „Stab“ mit der analytischen Lösung unterscheidet.

Lösung mit Maxima

Wir nutzen das Computer-Algebra-System Maxima zur Lösung. Das macht hier Sinn, weil wir die Herangehensweise mit der aus Stab vergleichen wollen – für die wir ebenfalls Maxima eingesetzt haben.

Declarations

Wir übernehmen alle Vereinbarungen und Parameter aus der Problemformulierung „Stab“.

Gleichgewichtsbedingungen

Für die Gleichgewichtsbedingung nach dem Prinzip der virtuellen Verrückungen

benötigen wir die virtuelle Formänderungsenergie und die virtuelle Arbeit der äußeren Kraft der äußeren Kräfte und Momente.

Mit den Konventionen für die Knoten-Verschiebungen aus Stab ist

.

Für gilt

mit den virtuellen Formänderungsarbeiten der vier Stäbe.

Dabei haben wir Anteile der Arbeit aus der [Biegung] und der Längs-Dehnung des Stabes.

Für den Stab k mit den Knoten I und J haben wir als Koodinaten der Knoten

und .

Damit haben wir

Für den Stab k definieren wir

sowie

und finden damit

mit der Element-Steifigkeitsmatrix des Elements k im k-Koordinatensystem

Transformation der Koordinaten in das globale System

In den Ausdrücken der virtuellen Formänderungsenergie stehen die Koordinaten des lokalen Koordinatensystems von k. Die müssen wir, wie in Stab mit der Euler-Drehmatrix ineinander überführen.

Dafür haben wir

mit der Transformationsmatrix

Es ist praktisch, an dieser Stelle die Abkürzung

für die Koordinaten eines Knoten im Referenzsystem einzuführen. Also ist

Damit wir für die Elementsteifigkeitsmatrix - mit beiden Anfangs- und Endknoten des Elements - vom "0"-System ins "k"-System transformieren, brauchen wir die neue Transformations-Matrix

Mit diesen ist

Die resultierenden Element-Steifigkeitsmatrizen sind im folgenden aufgeschreiben:

Element-Steigigkeitsmatrizen mit globalen Koordinaten
Element #1

Element #2

Element #3

Element #4

Wir sammeln nun alle Koordianten der Knoten im 0-System in

und schreiben die Gleichgewichtsbedingungen in der Form

an. Dabei kommen die Beiträge zur Gesamt-Steifigkeitsmatrix (hier noch in der Fassung ohne Berücksichtigung der Randbedingungen) aus den vier Beiträgen der virtuellen Formänderungsenergie - die wir hier farblich gekennzeichnet haben:

Einsortieren der Anteile der virtuellen Formänderungsenergie nach den Knoten.
Einsortieren der Anteile der virtuellen Formänderungsenergie nach den Knoten.

Einarbeitung der Randbedingungen

Die Randbedingungen arbeiten wir hier durch das Streichen der passenden Zeilen für sowie der passenden Spalten für ein.

Das resultierende Gleichungssystem ist dies:

Solving

Das Ergebnis ist - für die gleichen Parameter wie in Stab -

Und damit sind auch die Verläufe der Schnittlasten in den Stäben identisch - wir brauchen also die Ergebnisse nicht noch einmal auftragen.

Title

Text


1+1





Links

  • ...

Literature

  • ...