Gelöste Aufgaben/StaB

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 analytische Lösung der Felddifferentialgleichung.

Caption

In der klassischen Stabwerkstheorie (rechts) werden nur die Stab-Längskräfte berücksichtigt – das Ersatzmodell bildet die Verbindungen der Stäbe durch Gelenke ab. Trotzdem wurden die Stäbe in der Ausführung mit Knotenblechen vernietet oder verschweißt – von Gelenk also keine Spur.

Wir untersuchen, welchen Einfluss diese Vereinfachung auf die berechneten Belastungen der Struktur hat. Hier sollen die beiden Modell-Varianten – also die Modellierung mit Gelenken einerseits und Knotenblechen andererseits – miteinander verglichen werden. Gesucht ist die analytische Lösung der beiden Randwertprobleme und ein Vergleich.

Lösung mit Maxima

Wir nutzen das Computer-Algebra-System Maxima zur Lösung. Das macht Sinn, weil wir eine große Anzahl von Gleichungen verwalten müssen.

Declarations

Caption

Wir beginnen mit der Einführung von passenden Bezeichnungen für unser System:  

In Knoten I und II sind gelenkige Lager. Stäbe 1,2 und 3 sind im Knoten III verbunden, Stäbe 2 und 4 im Knoten II und Stäbe 3 und 4 im Knoten IV.

Caption

Jeder Stab – hier Stab k mit den Knoten I und J – trägt sein eigenes x_k, z_k – Koordinatensystem und ist im Winkel α_k zur Horizontalen eingebaut. Hier sind α₂ = 60° und α₃ = -60° sowie α₁ = 0 und α₄ = 0.

Als Systemparameter je Stab definieren wir vier Stablängen und die jeweils identischen Steifigkeiten für Dehnung und Biegung. Der Querschnitt der Stäbe sei ein doppeltsymmetrisches Profil mit der Querschnittsfläche $A=a^{2}$ und dem Flächenmoment $I=\eta {\frac {A^{2}}{12}}$.

Für das Referenzmodell, das unser System als Stabwerk mit gelenkigen Verbindungen erfasst, wählen wir die Stabkräfte S_i als Zugkräfte und schreiben folgende Gleichgewichtsbedingungen an:

${\frac {S_{3}}{2}}\operatorname {-} {\frac {S_{2}}{2}}\operatorname {-} {S_{1}}\operatorname {=} 0\operatorname {,} \operatorname {-} \left({\frac {{\sqrt {3}}{S_{3}}}{2}}\right)\operatorname {-} {\frac {{\sqrt {3}}{S_{2}}}{2}}\operatorname {=} 0\operatorname {,} \operatorname {-} {S_{4}}\operatorname {-} {\frac {S_{3}}{2}}\operatorname {=} 0\operatorname {,} F\operatorname {-} {\frac {{\sqrt {3}}{S_{3}}}{2}}\operatorname {=} 0$

Die Stabkräfte dazu sind

${S_{1}}\operatorname {=} {\frac {2F}{\sqrt {3}}}\operatorname {,} {S_{2}}\operatorname {=} \operatorname {-} \left({\frac {2F}{\sqrt {3}}}\right)\operatorname {,} {S_{3}}\operatorname {=} {\frac {2F}{\sqrt {3}}}\operatorname {,} {S_{4}}\operatorname {=} \operatorname {-} \left({\frac {F}{\sqrt {3}}}\right)$.

1+1

Ein Randwertproblem besteht aus den Gleichgewichtsbedingung im Feld - der Felddifferentialgleichung - und den Randbedingungen. Diese schreiben wir im Folgenden auf:

Die Lösung der Felddifferentialgleichung

Wenn wir die Stäbe in den Knotenpunkten fest verbinden, wird das System statisch überbestimmt – wir müssen die Verformungen des Systems berücksichtigen. Durch die geometrische Anordnung der Stäbe sind die Längs- und Querverformung der vier Stäbe miteinander gekoppelt. Wir müssen also zwei Differentialbeziehungen berücksichtigen:

• in Stab-Längsrichtung mit $u_{k}(x_{k})$ und
• in Stab-Querrichtung mit $w_{k}(x_{k})$.

Für Stab 3 sieht dies folgendermaßen aus:

Caption

Die Differentialbeziehungen des Gleichgewichts in u und w sowie deren Lösungen für unsere Stäbe sind

$EIw^{IV}=0$

$EAu''=0$

${\begin{array}{ccl}w_{k}(x_{k})&=&{\frac {\displaystyle a}{6}}C_{k,3}\left({\frac {x_{k}}{\ell _{k}}}\right)^{3}+{\frac {\displaystyle a}{2}}C_{k,2}\left({\frac {x_{k}}{\ell _{k}}}\right)^{2}+aC_{k,1}\left({\frac {x_{k}}{\ell _{k}}}\right)+aC_{k,0}\\w'_{k}(x_{k})&=&{\frac {\displaystyle a}{2}}C_{k,3}{\frac {x_{k}^{2}}{\ell _{k}^{3}}}+aC_{k,2}{\frac {x_{k}}{\ell _{k}^{2}}}+aC_{k,1}{\frac {1}{\ell _{k}^{1}}}\\w''_{k}(x_{k})&=&aC_{k,3}{\frac {x_{k}}{\ell _{k}^{3}}}+aC_{k,2}{\frac {1}{\ell _{k}^{2}}}\\w'''_{k}(x_{k})&=&aC_{k,3}{\frac {1}{\ell _{k}^{3}}}\\\end{array}}$

${\begin{array}{ccl}u_{k}(x_{k})&=&aD_{k,1}\left({\frac {x_{k}}{\ell _{k}}}\right)+aD_{k,0}\\u'_{k}(x_{k})&=&aD_{k,1}{\frac {1}{\ell _{k}}}\end{array}}$

Für den Kippwinkel des Balkens nutzen wir zusätzlich

$\varphi _{k}(x_{k})=w'_{k}(x_{k})$.

Die Schnittkräfte im Stab ergeben sich dann aus

$Q_{k}(x_{k})=-EIw'''_{k}(x_{k})$ $M_{k}(x_{k})=-EIw''_{k}(x_{k})$

$N_{k}(x_{k})=+EAu'_{k}(x_{k})$

Diese Schnittgrößen jeweils am positiven und negativen Schnittufer sind hier dargestellt:

Caption

1+1

Randbedingungen: Kinematische Verträglichkeit der Knotenverschiebungen

In den Knoten hat jetzt jeder Stab seine eigenen, lokalen Koordinaten der Verschiebung

$U_{i,k},W_{i,k}$,

wobei i den Knoten und k den Stab indexiert. Die Verdrehung des Knotens nennen wir

$\Phi _{i,k}$.

Die Verdrehungen der Knoten an den Enden sind in unserem 2D-System natürlich für jeden Stab gleich - wir führen den Index k hier nur ein, damit wir eine gleiche Nomenklatur für alle Koordinaten haben.

Caption

Damit diese Koordinaten miteinander verträglich sind, führen wir ein Referenzsystem „0“ mit den entsprechenden Koordinaten der Auslenkung in horizontale und vertikale Richtung ein:

Caption

Diese verschiedenen Koordinatensysteme müssen wir ineinander umrechnen. Das können wir anschaulich tun – wie hier im Beispiel für Knoten III und Stab 3 und finden

${U_{3,3}}={\frac {{\sqrt {3}}{W_{3,0}}}{2}}+{\frac {U_{3,0}}{2}},{W_{3,3}}={\frac {W_{3,0}}{2}}-{\frac {{\sqrt {3}}{U_{3,0}}}{2}}$

Oder wir wählen einen formalen Ansatz mit der Euer-Transformationsmatrix

${\underline {\underline {D}}}_{R}(\alpha _{k})=\left({\begin{array}{ccc}\cos(\alpha _{k})&\sin(\alpha _{k})&0\\-\sin(\alpha _{k})&\cos(\alpha _{k})&0\\0&0&1\end{array}}\right)$,

mit der wir die obige kinematische Beziehung aus

$\left({\begin{array}{c}U_{i0}\\W_{i0}\\\Phi _{i0}\\\end{array}}\right)={\underline {\underline {D}}}_{R}(\alpha _{k})\cdot \left({\begin{array}{c}U_{ik}\\W_{ik}\\\Phi _{ik}\\\end{array}}\right)$

gewinnen. Die gesuchten kinematischen Beziehungen erhalten wir für 4 Stäbe an ihren jeweiligen Enden zu

${\begin{array}{l}{U_{1,1}}={U_{1,0}},{W_{1,1}}={W_{1,0}},{{\Phi }_{1,1}}={{\Phi }_{1,0}},\\{U_{3,1}}={U_{3,0}},{W_{3,1}}={W_{3,0}},{{\Phi }_{3,1}}={{\Phi }_{3,0}},\\{U_{2,2}}={\frac {U_{2,0}}{2}}-{\frac {{\sqrt {3}}{W_{2,0}}}{2}},{W_{2,2}}={\frac {W_{2,0}}{2}}+{\frac {{\sqrt {3}}{U_{2,0}}}{2}},{{\Phi }_{2,2}}={{\Phi }_{2,0}},\\{U_{3,2}}={\frac {U_{3,0}}{2}}-{\frac {{\sqrt {3}}{W_{3,0}}}{2}},{W_{3,2}}={\frac {W_{3,0}}{2}}+{\frac {{\sqrt {3}}{U_{3,0}}}{2}},{{\Phi }_{3,2}}={{\Phi }_{3,0}},\\{U_{3,3}}={\frac {{\sqrt {3}}{W_{3,0}}}{2}}+{\frac {U_{3,0}}{2}},{W_{3,3}}={\frac {W_{3,0}}{2}}-{\frac {{\sqrt {3}}{U_{3,0}}}{2}},{{\Phi }_{3,3}}={{\Phi }_{3,0}},\\{U_{4,3}}={\frac {{\sqrt {3}}{W_{4,0}}}{2}}+{\frac {U_{4,0}}{2}},{W_{4,3}}={\frac {W_{4,0}}{2}}-{\frac {{\sqrt {3}}{U_{4,0}}}{2}},{{\Phi }_{4,3}}={{\Phi }_{4,0}},\\{U_{2,4}}={U_{2,0}},{W_{2,4}}={W_{2,0}},{{\Phi }_{2,4}}={{\Phi }_{2,0}},\\{U_{4,4}}={U_{4,0}},{W_{4,4}}={W_{4,0}},{{\Phi }_{4,4}}={{\Phi }_{4,0}}\end{array}}$

Alle Gleichungen zu den Drehwinkeln sind - im besten Sinne des Wortes - trivial. Wir nehmen sie hier trotzdem mit, um die Systematik nicht zu brechen.

1+1

                

Randbedingungen: die kinematischen Zwangsbedingungen an den Lagern

Zusätzlich zur kinematischen Verträglichkeit müssen wir auch noch die Randbedingungen in Knoten I und II anschrieben:

${\begin{array}{l}U_{1,0}=0\\w_{1,0}=0\\U_{2,0}=0\\w_{2,0}=0\\\end{array}}$

1+1

Randbedingungen: Kräfte und Momente in den Knoten

Es fehlen noch die Gleichgewichtsbedingungen für die Knotenpunkte. Wir erhalten für jeden Knoten – außer denen der Lager – zwei Kraft- und ein Momentengleichgewicht. Für die Lager-Knoten können wir nur die Momentengleichgewichte anschreiben – ohne neue Unbekannte für die Lager-Reaktionskräfte einzufügen.

Caption

Die Schnittbilder je Knoten erhalten wir, indem wir die Schnittgrößen an den entsprechenden Schnittufern anzeichnen und dann die Gleichgewichtsbedingungen je Knoten ansetzen.

Auch hier arbeiten wir der Eulerschen Drehmatrix, um diese Arbeit zu formalisieren.

Für Knoten III:

$-\left({\begin{array}{c}N_{1}(\ell _{1})\\Q_{1}(\ell _{1})\\M_{1}(\ell _{1})\end{array}}\right)-{\underline {\underline {D}}}_{R}(\alpha _{2})\cdot \left({\begin{array}{c}N_{2}(\ell _{2})\\Q_{2}(\ell _{2})\\M_{2}(\ell _{2})\end{array}}\right)+{\underline {\underline {D}}}_{R}(\alpha _{3})\cdot \left({\begin{array}{c}N_{3}(\ell _{3})\\Q_{3}(\ell _{3})\\M_{3}(\ell _{3})\end{array}}\right)=\left({\begin{array}{c}0\\0\\0\end{array}}\right)$

Für Knoten IV:

$-{\underline {\underline {D}}}_{R}(\alpha _{3})\cdot \left({\begin{array}{c}N_{3}(\ell _{3})\\Q_{3}(\ell _{3})\\M_{3}(\ell _{3})\end{array}}\right)-\left({\begin{array}{c}N_{4}(\ell _{4})\\Q_{4}(\ell _{4})\\M_{4}(\ell _{4})\end{array}}\right)=\left({\begin{array}{c}0\\-F\\0\end{array}}\right)$

Für Knoten I und II gewinnen wir jeweils nur die Momentenbilanzen, da wir die Lagerreaktionskräfte nicht als Unbekannte einführen wollen:

$M_{1}(0)=0,$

$M_{2}(0)+M_{4}(0)=0,$

Der gesuchte Satz an Gleichgewichtsbedingungen für die vier Knoten lautet dann

${\begin{array}{l}{M_{1}}(0)=0\\{M_{4}}(0)+{M_{2}}(0)=0\\-\left({\frac {\displaystyle {\sqrt {3}}{Q_{2}}\left({\ell _{0}}\right)}{\displaystyle 2}}\right)-{\frac {{\displaystyle N_{2}}\left({\ell _{0}}\right)}{\displaystyle 2}}-{N_{1}}\left({\frac {\displaystyle \ell _{0}}{\displaystyle 2}}\right)-{\frac {\displaystyle {\sqrt {3}}{Q_{3}}(0)}{\displaystyle 2}}+{\frac {\displaystyle {N_{3}}(0)}{\displaystyle 2}}=0\\-\left({\frac {\displaystyle {Q_{2}}\left({\ell _{0}}\right)}{\displaystyle 2}}\right)+{\frac {\displaystyle {\sqrt {3}}{N_{2}}\left({\ell _{0}}\right)}{\displaystyle 2}}-{Q_{1}}\left({\frac {\displaystyle \ell _{0}}{\displaystyle 2}}\right)+{\frac {\displaystyle {Q_{3}}(0)}{\displaystyle 2}}+{\frac {\displaystyle {\sqrt {3}}{N_{3}}(0)}{\displaystyle 2}}=0\\-{M_{2}}\left({\ell _{0}}\right)-{M_{1}}\left({\frac {\ell _{0}}{2}}\right)+{M_{3}}(0)=0\\-{N_{4}}\left({\ell _{0}}\right)+{\frac {\displaystyle {\sqrt {3}}{Q_{3}}\left({\ell _{0}}\right)}{\displaystyle 2}}-{\frac {\displaystyle {N_{3}}\left({\ell _{0}}\right)}{\displaystyle 2}}=0\\-{Q_{4}}\left({\ell _{0}}\right)-{\frac {\displaystyle {Q_{3}}\left({\ell _{0}}\right)}{\displaystyle 2}}-{\frac {\displaystyle {\sqrt {3}}{N_{3}}\left({\ell _{0}}\right)}{\displaystyle 2}}=-F\\-{M_{4}}\left({\ell _{0}}\right)-{M_{3}}\left({\ell _{0}}\right)=0\end{array}}$

1+1

Solving

Wir wählen hier zur Lösung der linearen Gleichungen einen Ansatz, bei dem wir die lokalen Ansatzfunktionen für die Schnittkräfte und die Verschiebungen schon in die Gleichungen der kinematischen Verträglichkeit und des Gleichgewichts einsetzten.

So erhalten wir z.B. aus

${U_{3,3}}={\frac {{\sqrt {3}}{W_{3,0}}}{2}}+{\frac {U_{3,0}}{2}}$

die Gleichung

$aD_{3,0}={\frac {{\sqrt {3}}{W_{3,0}}}{2}}+{\frac {U_{3,0}}{2}}$.

Und aus der Momentenbilanz für Knoten IV erhalten wir entsprechend

$-aC_{4,2}{\frac {EI}{\ell _{4}^{2}}}-aC_{2,2}{\frac {EI}{\ell _{2}^{2}}}=0$

                 Durch das Einsetzen der Lösungen für u_k(x_k) und w_k(x_k) haben wir als Unbekannte nur noch je Stab die 6 Integrationskonstanten C_ij, D_ij sowie je Knoten die Knotenverschiebungen Koordinaten U_k0, W_k0 und Φ_k0, also die 36 Unbekannten

Für das lineare Gleichungssystem

${\underline {\underline {A}}}\cdot {\underline {X}}={\underline {b}}$

ist also

${\underline {X}}^{T}=\left({U_{1,0}},{W_{1,0}},{{\Phi }_{1,0}},{U_{2,0}},{W_{2,0}},{{\Phi }_{2,0}},{U_{3,0}},{W_{3,0}},{{\Phi }_{3,0}},{U_{4,0}},{W_{4,0}},{{\Phi }_{4,0}},{C_{1,0}},{C_{1,1}},{C_{1,2}},{C_{1,3}},{C_{2,0}},{C_{2,1}},{C_{2,2}},{C_{2,3}},{C_{3,0}},{C_{3,1}},{C_{3,2}},{C_{3,3}},{C_{4,0}},{C_{4,1}}\operatorname {,} {C_{4,2}},{C_{4,3}},{D_{1,0}},{D_{1,1}},{D_{2,0}},{D_{2,1}},{D_{3,0}},{D_{3,1}},{D_{4,0}},{D_{4,1}}\right)$

Für deren Lösung stehen 4 Randbedingungen, 8 Gleichgewichtsbedingungen und 24 kinematische Verträglichkeitsbeziehungen – also 36 Gleichungen - zur Verfügung. Diese Gleichungen müssen sukzessive in die Matrix A und die rechte Seite des Gleichungssystems eingearbeitet werden.

Maxima liefert dann die Lösungen für unsere 36 Unbekannten.

Wir erhalten z.B. für die Normalkraft in Stab 3 den Ausdruck

${N_{3}}(0)={N_{3}}(\ell _{3})={\frac {11{{3}^{\frac {3}{2}}}{{\ell }_{0}^{2}}{{A}^{2}}{{\eta }^{2}}+518{\sqrt {3}}{{\ell }_{0}^{4}}A\eta +88{{3}^{\frac {3}{2}}}{{\ell }_{0}^{6}}}{27{{A}^{3}}{{\eta }^{3}}+459{{\ell }_{0}^{2}}{{A}^{2}}{{\eta }^{2}}+1008{{\ell }_{0}^{4}}A\eta +396{{\ell }_{0}^{6}}}}\cdot F$

der für η=0 (die Biegesteifigkeit ist Null)

$N_{3}={\frac {2F}{\sqrt {3}}}=S_{3}$.

liefert – also die Lösung für das vereinfachte Stabwerk.

1+1

Postprocessing

Wir tragen für alle Stäbe das Verhältnis N_i/S_i auf:

Caption

Wenn die Biegesteifigkeit der Stäbe klein wird, konvergiert unsere Lösung der Normalkräfte also gegen die einfache Stabwerkslösung.

Bleibt noch zu klären, ob die Spannungen durch die Biegung im Querschnitt ebenfalls vernachlässigbar klein bleiben. Dazu benötigen wir die Ausdrücke für

$\eta ={\frac {\displaystyle (H^{4}-h^{4})}{(H^{2}-h^{2})^{2}}},H=(a{\sqrt {\eta +1}})/{\sqrt {2}},h=(a{\sqrt {\eta -1}})/{\sqrt {2}}$,

die wir in die Gleichungen für die Normalspannungen im Querschnitt durch Normalkraft und Biegemoment einsetzten. Für die grafische Darstellung nutzen wir die Referenzspannung

$\sigma _{ref,i}={\frac {\displaystyle S_{i}}{A}}$

und berechnen damit die Spannungsverhältnisse

$s_{i,N}={\frac {\displaystyle \sigma _{N,i}}{\sigma _{ref,i}}},s_{i,M,0}={\frac {\displaystyle \sigma _{M,i,0}}{\sigma _{ref,i}}},s_{i,M,\ell }={\frac {\displaystyle \sigma _{M,i,\ell }}{\sigma _{ref,i}}}$

mit

$\sigma _{N,i}={\frac {\displaystyle N_{i}}{A}},\sigma _{M,i,0}={\frac {\displaystyle M_{i}(0)}{W}},\sigma _{M,i,\ell }={\frac {\displaystyle M_{i}(\ell _{i})}{W}}$

mit dem Widerstandsmoment

$W={\frac {\displaystyle I}{\left({\frac {\displaystyle H}{\displaystyle 2}}\right)}}$.

Die Auslenkungen der Knoten für η=1 sind

${\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.74{\frac {F}{Ea}}\\{W_{3,0}}=166.64{\frac {F}{Ea}}\\{{\Phi }_{3,0}}=2.06{\frac {F}{E{{a}^{2}}}}\\{U_{4,0}}=-57.74{\frac {F}{Ea}}\\{W_{4,0}}=366.64{\frac {F}{Ea}}\\{{\Phi }_{4,0}}=3.12{\frac {F}{E{{a}^{2}}}}\end{array}}$.

1+1

---

Links

• ...

Literature

• ...