(9 dazwischenliegende Versionen desselben Benutzers werden nicht angezeigt) Zeile 1:
Zeile 1:
[[index.php?title=Kategorie:Gelöste Aufgaben]]
[[Category :Gelöste Aufgaben]]
[[index.php?title=Kategorie:A*x=b]]
[[Category :A*x=b]]
[[index.php?title=Kategorie:Lineare Algebra]]
[[Category :Lineare Algebra]]
[[index.php?title=Kategorie:Numerische Lösung]]
[[Category :Numerische Lösung]]
[[index.php?title=Kategorie:Randwertproblem]]
[[Category :Randwertproblem]]
[[index.php?title=Kategorie:Prinzip der virtuellen Verrückungen]]
[[Category :Prinzip der virtuellen Verrückungen]]
[[index.php?title=Kategorie:Stab]][[index.php?title=Kategorie:Dehnstab]]
[[Category :Stab]]
[[index.php?title=Kategorie:Euler-Bernoulli-Balken]]
[[Category :Dehnstab]]
[[index.php?title=Kategorie:Finite-Elemente-Methode]]
[[Category :Euler-Bernoulli-Balken]]
[[index.php?title=Kategorie:Maxima]]
[[Category :Finite-Elemente-Methode]]
[[index.php?title=Kategorie:Stabwerk]]
[[Category :Maxima]]
[[Category :Stabwerk]]
==Aufgabenstellung==
==Aufgabenstellung==
Zeile 16:
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>
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.
Wir stellen das Modell des Stabwerks mit dem Prinzip der virtuellen Verrückungen auf und vergleichen, wie sich diese von der Herangehensweise aus „[https://numpedia.rzbt.haw-hamburg.de/index.php?title=Gel%C3%B6ste_Aufgaben/Stab Stab]“ mit der analytischen Lösung unterscheidet.
== Lösung mit Maxima ==
== Lösung mit Maxima ==
Zeile 25:
Zeile 26:
===Declarations===
===Declarations===
Wir übernehmen alle Vereinbarungen und Parameter aus der Problemformulierung „Stab“.
Wir übernehmen alle Vereinbarungen und Parameter aus der Problemformulierung „[https://numpedia.rzbt.haw-hamburg.de/index.php?title=Gel%C3%B6ste_Aufgaben/Stab Stab]“ .
===Gleichgewichtsbedingungen===
===Gleichgewichtsbedingungen===
Für die Gleichgewichtsbedingung nach dem [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Werkzeuge/Gleichgewichtsbedingungen/Arbeitsprinzipe_der_Analytischen_Mechanik/Prinzip_der_virtuellen_Verr%C3%BCckungen Prinzip der virtuellen Verrückungen] ist
Für die Gleichgewichtsbedingung nach dem [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Werkzeuge/Gleichgewichtsbedingungen/Arbeitsprinzipe_der_Analytischen_Mechanik/Prinzip_der_virtuellen_Verr%C3%BCckungen Prinzip der virtuellen Verrückungen]
::<math>
\begin{array}{ccc}
\delta W & \stackrel{!}{=} & 0\\
&=& \delta\Pi - \delta W^a
\end{array}
</math>
benötigen wir die virtuelle Formänderungsenergie <math>\delta \Pi</math> und die virtuelle Arbeit der äußeren Kraft <math>\delta W^a</math> der äußeren Kräfte und Momente.
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>
\delta W^a = +\delta W_{4,0} \cdot F
</math>.
Für <math>\delta \Pi</math> gilt
::<math>
\delta \Pi = \sum_{i=0}^4 \delta \Pi_i
</math>
mit den virtuellen Formänderungsarbeiten der vier Stäbe.
Dabei haben wir Anteile der Arbeit aus der [[https://numpedia.rzbt.haw-hamburg.de/index.php?title=Sources/Anleitungen/FEM-Formulierung_f%C3%BCr_den_Euler-Bernoulli-Balken 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
::<math>U_{I,k},W_{I,k},\Phi_{I,k}</math> und <math>U_{J,k},W_{J,k},\Phi_{J,k}</math>.
Damit haben wir
::<math>
\delta \Pi_k =
\left(\delta W_{I,k},\delta \Phi_{I,k},\delta \Phi_{J,k},\delta W_{J,k}\right)
\cdot
\frac{EI}{\ell_i^3}
\cdot
\begin{pmatrix}12 & 6\, {\ell_i} & -12 & 6\, {\ell_i}\\ 6\, {\ell_i} & 4\, {\ell_i^{2}} & -6\, {\ell_i} & 2\, {\ell_i^{2}}\\ -12 & -6\, {\ell_i} & 12 & -6\, {\ell_i}\\ 6\, {\ell_i} & 2\, {\ell_i^{2}} & -6\, {\ell_i} & 4\, {\ell_i^{2}}\end{pmatrix} \cdot
\left(\begin{array}{c}W_{I,k}\\ \Phi_{I,k}\\W_{J,k}\\ \Phi_{J,k}\end{array}\right)
+
\left(\delta U_{I,k},\delta U_{J,k}\right)
\cdot
\frac{EA}{\ell_i}
\cdot
\begin{pmatrix}1&-1\\-1&1\end{pmatrix}
\cdot
\left(\begin{array}{c}U_{I,k}\\U_{J,k}\end{array}\right)
</math>
W = 0
Für den Stab ''k'' definieren wir
= delta \Pi - \delta Wa
::<math>\underline{Q}_{k,k} = \left(\begin{array}{c}
benötigen wir die virtuelle Formänderungsenergie \delta \Pi und die virtuelle Arbeit der äußeren Kraft F.
U_{I,k}\\
Mit allen Konventionen für die Knoten-Verschiebungen ist
W_{I,k}\\
Delta Wa = -delta W_{4,0} *F
\Phi_{I,k}\\
Für \delta \Pi gilt
U_{J,k}\\
W_{J,k}\\
\Phi_{J,k}\\
\end{array}\right)</math> sowie <math>\delta\underline{Q}_{k,k} = \left(\begin{array}{c}
\ delta U_{I,k} \\
\delta W_{I,k}\\
\delta \Phi_{I,k}\\
\delta U_{J,k}\\
\ delta W_{J ,k }\\
\delta \Phi_{J,k}\\
\end{array}\right)
</math>
und finden damit
::<math>
\delta \Pi_k = \delta \underline{Q}_{k,k}^T \cdot \underline{\underline{K}}_{k,k} \cdot \underline{Q}_{k,k}
</math>
mit der Element-Steifigkeitsmatrix des Elements ''k'' im ''k''-Koordinatensystem
::<math>
::<math>
\begin{pmatrix}\frac{2 {{A}^{2}} E \eta }{3 {\ell_0}} & 0 & -\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\right) & 0 & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & 0 & 0 & 0\\
\underline {\underline {K }}_{k,k } =
0 & \frac{2 {{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 }{2 {{\ell}_{0}^{2}}}\right) & 0 & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
\frac{EI }{\ell_i ^3}
-\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\right) & -\left( \frac{{{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & \frac{{{A}^{2}} E \eta +3 {{\ell}_{0}^{2}} A E}{2 {{\ell}_{0}^{3}}}+\frac{8 {{A}^{2}} E \eta }{{{\ell}_{0}^{3}}} & 0 & -\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\right) & -\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}}}\\
\cdot
0 & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & 0 & \frac{3 {{A}^{2}} E \eta +{{\ell}_{0}^{2}} A E}{2 {{\ell}_{0}^{3}}}+\frac{2 A E}{{\ell_0}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{2 {{\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) \\
\begin {pmatrix }0 & 0 & 0 & 0 & 0 & 0\\
\frac{{{A}^{2}} E \eta }{3 {\ell_0}} & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\right) & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{2 {{\ell}_{0}^{2}}}\right) & \frac{4 {{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}}\\
0 & 12 & 6 \, {\ell_i } & 0 &-12 & 6 \, {\ell_i }\\
0 & -\left( \frac{{{A}^{2}} E \eta }{2 {{\ell}_{0}^{2}}}\right) & -\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}}}+\frac{{{A}^{2}} E \eta }{{{\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{3 {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) \\
0 & 6 \, {\ell_i } & 4\, {\ell_i ^{2}} & 0 & -6 \, {\ell_i } & 2\, {\ell_i ^{2}}\\
0 & 0 & \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{A E}{{\ell_0}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\\
0 & 0 & 0 & 0 & 0 & 0\\
0 & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & \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{3 {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}} & \frac{2 {{A}^{2}} E \eta }{3 {\ell_0}}\end{pmatrix}
0 &-12 & -6 \, {\ell_i } & 0 & 12 & -6 \, {\ell_i }\\
\cdot
0 & 6 \, {\ell_i } & 2\, {\ell_i ^{2}} & -0 & 6\, {\ell_i } & 4\, {\ell_i ^{2}}
\begin{pmatrix}
{W_{1,0}}\\
{U_{1,0}}\\
{{\Phi }_{1,0}}\\
{W_{2,0}}\\
{U_{2,0}}\\
{{\Phi }_{2,0}}\\
{W_{3,0}}\\
{U_{3,0}}\\
{{\Phi }_{3,0}}\\
{W_{4,0}}\\
{U_{4,0}}\\
{{\Phi }_{4,0}}
\end{pmatrix}
\end{pmatrix}
=
+
\begin{pmatrix}
\frac{EA}{\ell_i}
0\\
\cdot
0\\
\begin{pmatrix}1 & 0 & 0 &-1 & 0 & 0\\
0\\
0 & 0 & 0 & 0 & 0 & 0\\
0\\
0 & 0 & 0 & 0 & 0 & 0\\
0\\
-1 & 0 & 0 & 1 & 0 & 0\\
F\\
0 & 0 & 0 & 0 & 0 & 0\\
0\\
0 & 0 & 0 & 0 & 0 & 0
0
\end{pmatrix}
\end{pmatrix}
</math>
</math>
Hier kommt jetzt irgendein Text.
===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 [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Gel%C3%B6ste_Aufgaben/Stab Stab] mit der Euler-Drehmatrix ineinander überführen.
Dafür haben wir
::<math>
\left(\begin{array}{c}U_{I,0}\\W_{I,0}\\\Phi_{I,0}\end{array}\right) =
\underline{\underline{D}}_R(\alpha_k)\cdot
\left(\begin{array}{c}U_{I,k}\\W_{I,k}\\\Phi_{I,k}\end{array}\right)
</math>
mit der Transformationsmatrix
::<math>
::<math>
Some Text
\underline{\underline{D}}_R(\alpha_k) =
\left(\begin{array}{c} \cos(\alpha_k)& \sin(\alpha_k)& 0\\
-\sin(\alpha_k)& \cos(\alpha_k)& 0\\
0 & 0 & 1\end{array}\right).
</math>
Es ist praktisch, an dieser Stelle die Abkürzung
::<math>
\underline{q}_{I,0} = \left(\begin{array}{c}U_{I,0}\\W_{I,0}\\\Phi_{I,0}\end{array}\right)
</math>
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>
</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
{{MyCodeBlock|title=Title
::<math >
|text=Text
\underline {\underline {T}}(\alpha_k) : =
|code=
\left(\begin{array}{cc} \underline{\underline{D}}_R^T(\alpha_k)& \underline{\underline{0}}\\
<syntaxhighlight lang="lisp" line start=1>
\underline{\underline{0}} & \underline{\underline{D}}_R^T(\alpha_k)
1+1
\end{array}\right)
</syntaxhighlight>
</math>
}}
Mit diesen ist
:: <math>
\delta\Pi = \sum_{k =1}^4 \delta \underline{Q}_{k,0}^T\cdot
\underbrace{
\underline{\underline{T}}^T(\alpha_k)\cdot
\underline{\underline{K}}_{k,k}\cdot
\underline{\underline{T}}(\alpha_k)}_{:=
\underline{\underline{K}}_{k,0} }\cdot
\underline{Q}_{k,0}
</math >
Die resultierenden Element-Steifigkeitsmatrizen sind im folgenden aufgeschreiben:
<table class="wikitable mw-collapsible mw-collapsed" style="background-color:white; float: none; margin-right:14px;">
<table class="wikitable mw-collapsible mw-collapsed" style="background-color:white; float: none; margin-right:14px;">
Zeile 95:
Zeile 167:
<tr><td>
<tr><td>
<math>
<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}}}\\
\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\\
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}}\\
\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}}\\
Zeile 106:
Zeile 178:
<tr><td>
<tr><td>
<math>
<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}}}\\
\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{\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}}\\
\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}}\\
Zeile 117:
Zeile 189:
<tr><td>
<tr><td>
<math>
<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}}}\\
\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) \\
-\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}}\\
\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}}\\
Zeile 128:
Zeile 200:
<tr><td>
<tr><td>
<math>
<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}}}\\
\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\\
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}}\\
\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}}\\
Zeile 137:
Zeile 209:
</td></tr>
</td></tr>
</table>
</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>
\delta W = \delta\underline{Q}_{0}^T\cdot\left(\underline{\underline{K}}_{0}\cdot\underline{Q}_{0}
- \underline{P}\right) \stackrel{!}{=} 0
</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.]]
===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:
::<math>
\begin{pmatrix}
\frac{2 {{A}^{2}} E \eta }{3 {\ell_0}} & 0 & 0 & -\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\right) & \frac{{{A}^{2}} E \eta }{3 {\ell_0}} & 0 & 0 & 0\\
0 & \frac{2 {{A}^{2}} E \eta }{3 {\ell_0}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & -\left( \frac{{{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & 0 & -\left( \frac{{{A}^{2}} E \eta }{2 {{\ell}_{0}^{2}}}\right) & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
0 & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & \frac{3 {{A}^{2}} E \eta +{{\ell}_{0}^{2}} A E}{2 {{\ell}_{0}^{3}}}+\frac{2 A E}{{\ell_0}} & 0 & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{2 {{\ell}_{0}^{2}}}\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 -\sqrt{3} {{\ell}_{0}^{2}} A E}{4 {{\ell}_{0}^{3}}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) \\
-\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\right) & -\left( \frac{{{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & 0 & \frac{{{A}^{2}} E \eta +3 {{\ell}_{0}^{2}} A E}{2 {{\ell}_{0}^{3}}}+\frac{8 {{A}^{2}} E \eta }{{{\ell}_{0}^{3}}} & -\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\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 +3 {{\ell}_{0}^{2}} A E}{4 {{\ell}_{0}^{3}}}\right) & \frac{{{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\\
\frac{{{A}^{2}} E \eta }{3 {\ell_0}} & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{2 {{\ell}_{0}^{2}}}\right) & -\left( \frac{2 {{A}^{2}} E \eta }{{{\ell}_{0}^{2}}}\right) & \frac{4 {{A}^{2}} E \eta }{3 {\ell_0}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}} & -\left( \frac{{{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & \frac{{{A}^{2}} E \eta }{6 {\ell_0}}\\
0 & 0 & -\left( \frac{3 {{A}^{2}} E \eta +{{\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{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}} & \frac{3 {{A}^{2}} E \eta +{{\ell}_{0}^{2}} A E}{4 {{\ell}_{0}^{3}}}+\frac{A E}{{\ell_0}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta -\sqrt{3} {{\ell}_{0}^{2}} A E}{4 {{\ell}_{0}^{3}}}\right) & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\\
0 & -\left( \frac{{{A}^{2}} E \eta }{2 {{\ell}_{0}^{2}}}\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 +3 {{\ell}_{0}^{2}} A E}{4 {{\ell}_{0}^{3}}}\right) & -\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) & \frac{{{A}^{2}} E \eta +3 {{\ell}_{0}^{2}} A E}{4 {{\ell}_{0}^{3}}}+\frac{{{A}^{2}} E \eta }{{{\ell}_{0}^{3}}} & -\left( \frac{3 {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) \\
0 & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & -\left( \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & \frac{{{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}} & \frac{{{A}^{2}} E \eta }{6 {\ell_0}} & \frac{\sqrt{3} {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}} & -\left( \frac{3 {{A}^{2}} E \eta }{4 {{\ell}_{0}^{2}}}\right) & \frac{2 {{A}^{2}} E \eta }{3 {\ell_0}}
\end{pmatrix}
\cdot
\begin{pmatrix}
{{\Phi }_{1,0}}\\
{{\Phi }_{2,0}}\\
{W_{3,0}}\\
{U_{3,0}}\\
{{\Phi }_{3,0}}\\
{W_{4,0}}\\
{U_{4,0}}\\
{{\Phi }_{4,0}}
\end{pmatrix}
=
\begin{pmatrix}
0\\
0\\
0\\
0\\
0\\
F\\
0\\
0
\end{pmatrix}
</math>
===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>
\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>
Und damit sind auch die Verläufe der Schnittlasten in den Stäben identisch - wir brauchen also die Ergebnisse nicht noch einmal auftragen.
<!-------------------------------------------------------------------------------->
{{MyCodeBlock|title=Title
|text=Text
|code=
<syntaxhighlight lang="lisp" line start=1>
1+1
</syntaxhighlight>
}}
<hr/>
<hr/>
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
δ
W
=
!
0
=
δ
Π
−
δ
W
a
{\displaystyle {\begin{array}{ccc}\delta W&{\stackrel {!}{=}}&0\\&=&\delta \Pi -\delta W^{a}\end{array}}}
benötigen wir die virtuelle Formänderungsenergie
δ
Π
{\displaystyle \delta \Pi }
und die virtuelle Arbeit der äußeren Kraft
δ
W
a
{\displaystyle \delta W^{a}}
der äußeren Kräfte und Momente.
Mit den Konventionen für die Knoten-Verschiebungen aus Stab ist
δ
W
a
=
+
δ
W
4
,
0
⋅
F
{\displaystyle \delta W^{a}=+\delta W_{4,0}\cdot F}
.
Für
δ
Π
{\displaystyle \delta \Pi }
gilt
δ
Π
=
∑
i
=
0
4
δ
Π
i
{\displaystyle \delta \Pi =\sum _{i=0}^{4}\delta \Pi _{i}}
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
U
I
,
k
,
W
I
,
k
,
Φ
I
,
k
{\displaystyle U_{I,k},W_{I,k},\Phi _{I,k}}
und
U
J
,
k
,
W
J
,
k
,
Φ
J
,
k
{\displaystyle U_{J,k},W_{J,k},\Phi _{J,k}}
.
Damit haben wir
δ
Π
k
=
(
δ
W
I
,
k
,
δ
Φ
I
,
k
,
δ
Φ
J
,
k
,
δ
W
J
,
k
)
⋅
E
I
ℓ
i
3
⋅
(
12
6
ℓ
i
−
12
6
ℓ
i
6
ℓ
i
4
ℓ
i
2
−
6
ℓ
i
2
ℓ
i
2
−
12
−
6
ℓ
i
12
−
6
ℓ
i
6
ℓ
i
2
ℓ
i
2
−
6
ℓ
i
4
ℓ
i
2
)
⋅
(
W
I
,
k
Φ
I
,
k
W
J
,
k
Φ
J
,
k
)
+
(
δ
U
I
,
k
,
δ
U
J
,
k
)
⋅
E
A
ℓ
i
⋅
(
1
−
1
−
1
1
)
⋅
(
U
I
,
k
U
J
,
k
)
{\displaystyle \delta \Pi _{k}=\left(\delta W_{I,k},\delta \Phi _{I,k},\delta \Phi _{J,k},\delta W_{J,k}\right)\cdot {\frac {EI}{\ell _{i}^{3}}}\cdot {\begin{pmatrix}12&6\,{\ell _{i}}&-12&6\,{\ell _{i}}\\6\,{\ell _{i}}&4\,{\ell _{i}^{2}}&-6\,{\ell _{i}}&2\,{\ell _{i}^{2}}\\-12&-6\,{\ell _{i}}&12&-6\,{\ell _{i}}\\6\,{\ell _{i}}&2\,{\ell _{i}^{2}}&-6\,{\ell _{i}}&4\,{\ell _{i}^{2}}\end{pmatrix}}\cdot \left({\begin{array}{c}W_{I,k}\\\Phi _{I,k}\\W_{J,k}\\\Phi _{J,k}\end{array}}\right)+\left(\delta U_{I,k},\delta U_{J,k}\right)\cdot {\frac {EA}{\ell _{i}}}\cdot {\begin{pmatrix}1&-1\\-1&1\end{pmatrix}}\cdot \left({\begin{array}{c}U_{I,k}\\U_{J,k}\end{array}}\right)}
Für den Stab k definieren wir
Q
_
k
,
k
=
(
U
I
,
k
W
I
,
k
Φ
I
,
k
U
J
,
k
W
J
,
k
Φ
J
,
k
)
{\displaystyle {\underline {Q}}_{k,k}=\left({\begin{array}{c}U_{I,k}\\W_{I,k}\\\Phi _{I,k}\\U_{J,k}\\W_{J,k}\\\Phi _{J,k}\\\end{array}}\right)}
sowie
δ
Q
_
k
,
k
=
(
δ
U
I
,
k
δ
W
I
,
k
δ
Φ
I
,
k
δ
U
J
,
k
δ
W
J
,
k
δ
Φ
J
,
k
)
{\displaystyle \delta {\underline {Q}}_{k,k}=\left({\begin{array}{c}\delta U_{I,k}\\\delta W_{I,k}\\\delta \Phi _{I,k}\\\delta U_{J,k}\\\delta W_{J,k}\\\delta \Phi _{J,k}\\\end{array}}\right)}
und finden damit
δ
Π
k
=
δ
Q
_
k
,
k
T
⋅
K
_
_
k
,
k
⋅
Q
_
k
,
k
{\displaystyle \delta \Pi _{k}=\delta {\underline {Q}}_{k,k}^{T}\cdot {\underline {\underline {K}}}_{k,k}\cdot {\underline {Q}}_{k,k}}
mit der Element-Steifigkeitsmatrix des Elements k im k -Koordinatensystem
K
_
_
k
,
k
=
E
I
ℓ
i
3
⋅
(
0
0
0
0
0
0
0
12
6
ℓ
i
0
−
12
6
ℓ
i
0
6
ℓ
i
4
ℓ
i
2
0
−
6
ℓ
i
2
ℓ
i
2
0
0
0
0
0
0
0
−
12
−
6
ℓ
i
0
12
−
6
ℓ
i
0
6
ℓ
i
2
ℓ
i
2
−
0
6
ℓ
i
4
ℓ
i
2
)
+
E
A
ℓ
i
⋅
(
1
0
0
−
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
−
1
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
)
{\displaystyle {\underline {\underline {K}}}_{k,k}={\frac {EI}{\ell _{i}^{3}}}\cdot {\begin{pmatrix}0&0&0&0&0&0\\0&12&6\,{\ell _{i}}&0&-12&6\,{\ell _{i}}\\0&6\,{\ell _{i}}&4\,{\ell _{i}^{2}}&0&-6\,{\ell _{i}}&2\,{\ell _{i}^{2}}\\0&0&0&0&0&0\\0&-12&-6\,{\ell _{i}}&0&12&-6\,{\ell _{i}}\\0&6\,{\ell _{i}}&2\,{\ell _{i}^{2}}&-0&6\,{\ell _{i}}&4\,{\ell _{i}^{2}}\end{pmatrix}}+{\frac {EA}{\ell _{i}}}\cdot {\begin{pmatrix}1&0&0&-1&0&0\\0&0&0&0&0&0\\0&0&0&0&0&0\\-1&0&0&1&0&0\\0&0&0&0&0&0\\0&0&0&0&0&0\end{pmatrix}}}
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
(
U
I
,
0
W
I
,
0
Φ
I
,
0
)
=
D
_
_
R
(
α
k
)
⋅
(
U
I
,
k
W
I
,
k
Φ
I
,
k
)
{\displaystyle \left({\begin{array}{c}U_{I,0}\\W_{I,0}\\\Phi _{I,0}\end{array}}\right)={\underline {\underline {D}}}_{R}(\alpha _{k})\cdot \left({\begin{array}{c}U_{I,k}\\W_{I,k}\\\Phi _{I,k}\end{array}}\right)}
mit der Transformationsmatrix
D
_
_
R
(
α
k
)
=
(
cos
(
α
k
)
sin
(
α
k
)
0
−
sin
(
α
k
)
cos
(
α
k
)
0
0
0
1
)
.
{\displaystyle {\underline {\underline {D}}}_{R}(\alpha _{k})=\left({\begin{array}{c}\cos(\alpha _{k})&\sin(\alpha _{k})&0\\-\sin(\alpha _{k})&\cos(\alpha _{k})&0\\0&0&1\end{array}}\right).}
Es ist praktisch, an dieser Stelle die Abkürzung
q
_
I
,
0
=
(
U
I
,
0
W
I
,
0
Φ
I
,
0
)
{\displaystyle {\underline {q}}_{I,0}=\left({\begin{array}{c}U_{I,0}\\W_{I,0}\\\Phi _{I,0}\end{array}}\right)}
für die Koordinaten eines Knoten im Referenzsystem einzuführen. Also ist
(
U
I
,
k
W
I
,
k
Φ
I
,
k
)
=
D
_
_
R
T
(
α
k
)
⋅
q
_
I
,
0
.
{\displaystyle \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}.}
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
T
_
_
(
α
k
)
:=
(
D
_
_
R
T
(
α
k
)
0
_
_
0
_
_
D
_
_
R
T
(
α
k
)
)
{\displaystyle {\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)}
Mit diesen ist
δ
Π
=
∑
k
=
1
4
δ
Q
_
k
,
0
T
⋅
T
_
_
T
(
α
k
)
⋅
K
_
_
k
,
k
⋅
T
_
_
(
α
k
)
⏟
:=
K
_
_
k
,
0
⋅
Q
_
k
,
0
{\displaystyle \delta \Pi =\sum _{k=1}^{4}\delta {\underline {Q}}_{k,0}^{T}\cdot \underbrace {{\underline {\underline {T}}}^{T}(\alpha _{k})\cdot {\underline {\underline {K}}}_{k,k}\cdot {\underline {\underline {T}}}(\alpha _{k})} _{:={\underline {\underline {K}}}_{k,0}}\cdot {\underline {Q}}_{k,0}}
Die resultierenden Element-Steifigkeitsmatrizen sind im folgenden aufgeschreiben:
Element-Steigigkeitsmatrizen mit globalen Koordinaten
Element #1
K
_
_
1
,
0
=
(
8
A
2
E
η
ℓ
0
3
0
2
A
2
E
η
ℓ
0
2
−
(
8
A
2
E
η
ℓ
0
3
)
0
2
A
2
E
η
ℓ
0
2
0
2
A
E
ℓ
0
0
0
−
(
2
A
E
ℓ
0
)
0
2
A
2
E
η
ℓ
0
2
0
2
A
2
E
η
3
ℓ
0
−
(
2
A
2
E
η
ℓ
0
2
)
0
A
2
E
η
3
ℓ
0
−
(
8
A
2
E
η
ℓ
0
3
)
0
−
(
2
A
2
E
η
ℓ
0
2
)
8
A
2
E
η
ℓ
0
3
0
−
(
2
A
2
E
η
ℓ
0
2
)
0
−
(
2
A
E
ℓ
0
)
0
0
2
A
E
ℓ
0
0
2
A
2
E
η
ℓ
0
2
0
A
2
E
η
3
ℓ
0
−
(
2
A
2
E
η
ℓ
0
2
)
0
2
A
2
E
η
3
ℓ
0
)
{\displaystyle {\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 {2AE}{\ell _{0}}}&0&0&-\left({\frac {2AE}{\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 {2AE}{\ell _{0}}}\right)&0&0&{\frac {2AE}{\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}}}
Element #2
K
_
_
2
,
0
=
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
A
2
E
η
4
ℓ
0
2
−
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
A
2
E
η
4
ℓ
0
2
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
3
A
2
E
η
4
ℓ
0
2
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
4
ℓ
0
2
A
2
E
η
4
ℓ
0
2
3
A
2
E
η
4
ℓ
0
2
A
2
E
η
3
ℓ
0
−
(
A
2
E
η
4
ℓ
0
2
)
−
(
3
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
6
ℓ
0
−
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
A
2
E
η
4
ℓ
0
2
)
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
3
A
2
E
η
4
ℓ
0
2
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
−
(
3
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
4
ℓ
0
2
3
A
2
E
η
4
ℓ
0
2
A
2
E
η
6
ℓ
0
−
(
A
2
E
η
4
ℓ
0
2
)
−
(
3
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
3
ℓ
0
)
{\displaystyle {\underline {\underline {K}}}_{2,0}={\begin{pmatrix}{\frac {{{A}^{2}}E\eta +3{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}&{\frac {{{A}^{2}}E\eta }{4{\ell _{0}^{2}}}}&-\left({\frac {{{A}^{2}}E\eta +3{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}\right)&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}&{\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}\right)&-\left({\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}\right)&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}\right)&-\left({\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}&{\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}}
Element #3
K
_
_
3
,
0
=
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
A
2
E
η
4
ℓ
0
2
−
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
A
2
E
η
4
ℓ
0
2
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
−
(
3
A
2
E
η
4
ℓ
0
2
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
3
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
4
ℓ
0
2
−
(
3
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
3
ℓ
0
−
(
A
2
E
η
4
ℓ
0
2
)
3
A
2
E
η
4
ℓ
0
2
A
2
E
η
6
ℓ
0
−
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
A
2
E
η
4
ℓ
0
2
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
4
ℓ
0
2
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
3
A
2
E
η
4
ℓ
0
2
A
2
E
η
4
ℓ
0
2
−
(
3
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
6
ℓ
0
−
(
A
2
E
η
4
ℓ
0
2
)
3
A
2
E
η
4
ℓ
0
2
A
2
E
η
3
ℓ
0
)
{\displaystyle {\underline {\underline {K}}}_{3,0}={\begin{pmatrix}{\frac {{{A}^{2}}E\eta +3{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}\right)&{\frac {{{A}^{2}}E\eta }{4{\ell _{0}^{2}}}}&-\left({\frac {{{A}^{2}}E\eta +3{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}\right)&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}\right)&{\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}&-\left({\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}\right)&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}&-\left({\frac {{{A}^{2}}E\eta }{4{\ell _{0}^{2}}}}\right)&{\frac {{{A}^{2}}E\eta +3{\ell _{0}^{2}}AE}{4{\ell _{0}^{3}}}}&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}&-\left({\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}AE}{4{\ell _{0}^{3}}}}\right)&{\frac {3{{A}^{2}}E\eta +{\ell _{0}^{2}}AE}{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}}}
Element #4
K
_
_
4
,
0
=
(
A
2
E
η
ℓ
0
3
0
A
2
E
η
2
ℓ
0
2
−
(
A
2
E
η
ℓ
0
3
)
0
A
2
E
η
2
ℓ
0
2
0
A
E
ℓ
0
0
0
−
(
A
E
ℓ
0
)
0
A
2
E
η
2
ℓ
0
2
0
A
2
E
η
3
ℓ
0
−
(
A
2
E
η
2
ℓ
0
2
)
0
A
2
E
η
6
ℓ
0
−
(
A
2
E
η
ℓ
0
3
)
0
−
(
A
2
E
η
2
ℓ
0
2
)
A
2
E
η
ℓ
0
3
0
−
(
A
2
E
η
2
ℓ
0
2
)
0
−
(
A
E
ℓ
0
)
0
0
A
E
ℓ
0
0
A
2
E
η
2
ℓ
0
2
0
A
2
E
η
6
ℓ
0
−
(
A
2
E
η
2
ℓ
0
2
)
0
A
2
E
η
3
ℓ
0
)
{\displaystyle {\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 {AE}{\ell _{0}}}&0&0&-\left({\frac {AE}{\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 {AE}{\ell _{0}}}\right)&0&0&{\frac {AE}{\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}}}
Wir sammeln nun alle Koordianten der Knoten im 0 -System in
Q
_
0
T
=
(
q
_
1
,
0
T
,
q
_
2
,
0
T
,
q
_
3
,
0
T
,
q
_
4
,
0
T
)
{\displaystyle {\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)}
und schreiben die Gleichgewichtsbedingungen in der Form
δ
W
=
δ
Q
_
0
T
⋅
(
K
_
_
0
⋅
Q
_
0
−
P
_
)
=
!
0
{\displaystyle \delta W=\delta {\underline {Q}}_{0}^{T}\cdot \left({\underline {\underline {K}}}_{0}\cdot {\underline {Q}}_{0}-{\underline {P}}\right){\stackrel {!}{=}}0}
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.
Einarbeitung der Randbedingungen
Die Randbedingungen arbeiten wir hier durch das Streichen der passenden Zeilen für
δ
U
1
,
0
,
δ
W
1
,
0
,
δ
U
2
,
0
,
δ
W
2
,
0
{\displaystyle \delta U_{1,0},\delta W_{1,0},\delta U_{2,0},\delta W_{2,0}}
sowie der passenden Spalten für
U
1
,
0
,
W
1
,
0
,
U
2
,
0
,
W
2
,
0
{\displaystyle U_{1,0},W_{1,0},U_{2,0},W_{2,0}}
ein.
Das resultierende Gleichungssystem ist dies:
(
2
A
2
E
η
3
ℓ
0
0
0
−
(
2
A
2
E
η
ℓ
0
2
)
A
2
E
η
3
ℓ
0
0
0
0
0
2
A
2
E
η
3
ℓ
0
−
(
3
A
2
E
η
4
ℓ
0
2
)
−
(
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
6
ℓ
0
0
−
(
A
2
E
η
2
ℓ
0
2
)
A
2
E
η
6
ℓ
0
0
−
(
3
A
2
E
η
4
ℓ
0
2
)
3
A
2
E
η
+
ℓ
0
2
A
E
2
ℓ
0
3
+
2
A
E
ℓ
0
0
−
(
3
A
2
E
η
2
ℓ
0
2
)
−
(
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
3
A
2
E
η
4
ℓ
0
2
)
−
(
2
A
2
E
η
ℓ
0
2
)
−
(
A
2
E
η
4
ℓ
0
2
)
0
A
2
E
η
+
3
ℓ
0
2
A
E
2
ℓ
0
3
+
8
A
2
E
η
ℓ
0
3
−
(
2
A
2
E
η
ℓ
0
2
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
)
A
2
E
η
4
ℓ
0
2
A
2
E
η
3
ℓ
0
A
2
E
η
6
ℓ
0
−
(
3
A
2
E
η
2
ℓ
0
2
)
−
(
2
A
2
E
η
ℓ
0
2
)
4
A
2
E
η
3
ℓ
0
3
A
2
E
η
4
ℓ
0
2
−
(
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
6
ℓ
0
0
0
−
(
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
3
A
2
E
η
4
ℓ
0
2
3
A
2
E
η
+
ℓ
0
2
A
E
4
ℓ
0
3
+
A
E
ℓ
0
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
3
A
2
E
η
4
ℓ
0
2
0
−
(
A
2
E
η
2
ℓ
0
2
)
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
−
(
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
)
−
(
A
2
E
η
4
ℓ
0
2
)
−
(
3
A
2
E
η
−
3
ℓ
0
2
A
E
4
ℓ
0
3
)
A
2
E
η
+
3
ℓ
0
2
A
E
4
ℓ
0
3
+
A
2
E
η
ℓ
0
3
−
(
3
A
2
E
η
4
ℓ
0
2
)
0
A
2
E
η
6
ℓ
0
−
(
3
A
2
E
η
4
ℓ
0
2
)
A
2
E
η
4
ℓ
0
2
A
2
E
η
6
ℓ
0
3
A
2
E
η
4
ℓ
0
2
−
(
3
A
2
E
η
4
ℓ
0
2
)
2
A
2
E
η
3
ℓ
0
)
⋅
(
Φ
1
,
0
Φ
2
,
0
W
3
,
0
U
3
,
0
Φ
3
,
0
W
4
,
0
U
4
,
0
Φ
4
,
0
)
=
(
0
0
0
0
0
F
0
0
)
{\displaystyle {\begin{pmatrix}{\frac {2{{A}^{2}}E\eta }{3{\ell _{0}}}}&0&0&-\left({\frac {2{{A}^{2}}E\eta }{{\ell }_{0}^{2}}}\right)&{\frac {{{A}^{2}}E\eta }{3{\ell _{0}}}}&0&0&0\\0&{\frac {2{{A}^{2}}E\eta }{3{\ell _{0}}}}&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&-\left({\frac {{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&{\frac {{{A}^{2}}E\eta }{6{\ell _{0}}}}&0&-\left({\frac {{{A}^{2}}E\eta }{2{{\ell }_{0}^{2}}}}\right)&{\frac {{{A}^{2}}E\eta }{6{\ell _{0}}}}\\0&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&{\frac {3{{A}^{2}}E\eta +{{\ell }_{0}^{2}}AE}{2{{\ell }_{0}^{3}}}}+{\frac {2AE}{\ell _{0}}}&0&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta }{2{{\ell }_{0}^{2}}}}\right)&-\left({\frac {3{{A}^{2}}E\eta +{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}\right)&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)\\-\left({\frac {2{{A}^{2}}E\eta }{{\ell }_{0}^{2}}}\right)&-\left({\frac {{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&0&{\frac {{{A}^{2}}E\eta +3{{\ell }_{0}^{2}}AE}{2{{\ell }_{0}^{3}}}}+{\frac {8{{A}^{2}}E\eta }{{\ell }_{0}^{3}}}&-\left({\frac {2{{A}^{2}}E\eta }{{\ell }_{0}^{2}}}\right)&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}&-\left({\frac {{{A}^{2}}E\eta +3{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}\right)&{\frac {{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\\{\frac {{{A}^{2}}E\eta }{3{\ell _{0}}}}&{\frac {{{A}^{2}}E\eta }{6{\ell _{0}}}}&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta }{2{{\ell }_{0}^{2}}}}\right)&-\left({\frac {2{{A}^{2}}E\eta }{{\ell }_{0}^{2}}}\right)&{\frac {4{{A}^{2}}E\eta }{3{\ell _{0}}}}&{\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}&-\left({\frac {{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&{\frac {{{A}^{2}}E\eta }{6{\ell _{0}}}}\\0&0&-\left({\frac {3{{A}^{2}}E\eta +{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}\right)&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}&{\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}&{\frac {3{{A}^{2}}E\eta +{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}+{\frac {AE}{\ell _{0}}}&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}\right)&{\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\\0&-\left({\frac {{{A}^{2}}E\eta }{2{{\ell }_{0}^{2}}}}\right)&{\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}&-\left({\frac {{{A}^{2}}E\eta +3{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}\right)&-\left({\frac {{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta -{\sqrt {3}}{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}\right)&{\frac {{{A}^{2}}E\eta +3{{\ell }_{0}^{2}}AE}{4{{\ell }_{0}^{3}}}}+{\frac {{{A}^{2}}E\eta }{{\ell }_{0}^{3}}}&-\left({\frac {3{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)\\0&{\frac {{{A}^{2}}E\eta }{6{\ell _{0}}}}&-\left({\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&{\frac {{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}&{\frac {{{A}^{2}}E\eta }{6{\ell _{0}}}}&{\frac {{\sqrt {3}}{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}&-\left({\frac {3{{A}^{2}}E\eta }{4{{\ell }_{0}^{2}}}}\right)&{\frac {2{{A}^{2}}E\eta }{3{\ell _{0}}}}\end{pmatrix}}\cdot {\begin{pmatrix}{{\Phi }_{1,0}}\\{{\Phi }_{2,0}}\\{W_{3,0}}\\{U_{3,0}}\\{{\Phi }_{3,0}}\\{W_{4,0}}\\{U_{4,0}}\\{{\Phi }_{4,0}}\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\\0\\F\\0\\0\end{pmatrix}}}
Solving
Das Ergebnis ist - für die gleichen Parameter wie in Stab -
U
1
,
0
=
0
W
1
,
0
=
0
Φ
1
,
0
=
3.97
F
E
a
2
U
2
,
0
=
0
W
2
,
0
=
0
Φ
2
,
0
=
2.45
F
E
a
2
U
3
,
0
=
57.7
F
E
a
W
3
,
0
=
1.67
⋅
10
2
F
E
a
Φ
3
,
0
=
2.06
F
E
a
2
U
4
,
0
=
−
57.7
F
E
a
W
4
,
0
=
3.67
⋅
10
2
F
E
a
Φ
4
,
0
=
3.12
F
E
a
2
{\displaystyle {\begin{array}{l}U_{1,0}=0\\W_{1,0}=0\\\Phi _{1,0}=3.97{\frac {F}{Ea^{2}}}\\U_{2,0}=0\\W_{2,0}=0\\\Phi _{2,0}=2.45{\frac {F}{Ea^{2}}}\\U_{3,0}=57.7{\frac {F}{Ea}}\\W_{3,0}=1.67\cdot 10^{2}{\frac {F}{Ea}}\\\Phi _{3,0}=2.06{\frac {F}{Ea^{2}}}\\U_{4,0}=-57.7{\frac {F}{Ea}}\\W_{4,0}=3.67\cdot 10^{2}{\frac {F}{Ea}}\\\Phi _{4,0}=3.12{\frac {F}{Ea^{2}}}\end{array}}}
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
Links
Literature