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.
Caption
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
_
=
(
U
I
,
k
W
I
,
k
Φ
I
,
k
U
J
,
k
W
J
,
k
Φ
J
,
k
)
{\displaystyle {\underline {Q}}=\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
_
=
(
δ
U
I
,
k
δ
W
I
,
k
δ
Φ
I
,
k
δ
U
J
,
k
δ
W
J
,
k
δ
Φ
J
,
k
)
{\displaystyle \delta {\underline {Q}}=\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
_
T
⋅
K
_
_
k
⋅
Q
_
{\displaystyle \delta \Pi _{k}=\delta {\underline {Q}}^{T}\cdot {\underline {\underline {K}}}_{k}\cdot {\underline {Q}}}
mit der Element-Steifigkeitsmatrix
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}={\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}}}
(
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
)
⋅
(
W
1
,
0
U
1
,
0
Φ
1
,
0
W
2
,
0
U
2
,
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}{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}}={\begin{pmatrix}0\\0\\0\\0\\0\\F\\0\\0\end{pmatrix}}}
Hier kommt jetzt irgendein Text.
S
o
m
e
T
e
x
t
{\displaystyle SomeText}
Title
Text
Element-Steigigkeitsmatrizen mit globalen Koordinaten
Element #1
k
1
=
(
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 {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 {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
=
(
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 {k_{2}}={\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
=
(
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 {k_{3}}={\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
=
(
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 {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 {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}}}
Links
Literature