Gelöste Aufgaben/Kw60

Aufgabenstellung

Der Euler-Bernoulli-Balken ABC (Länge ℓ₀, Biegesteifigkeikeit EI, Querschnittsflläche A) ist am rechten Rand fest eingespannt und wird durch das Gewicht

• aufgrund seiner eigenen Masse m₀ sowie
• aufgrund von zwei Personen der Massen m_A bzw. m_B

belastet.

Lageplan

Gesucht ist eine Näherungslösung mit der Methode der Finiten Elemente und zwei Elementen. Dabei sollen die Verläufe der Biegelinie und die Schnittlasten mit der analytischen Lösung verlglichen werden.

Es ist ℓ_B =  ℓ₀ /3, m_A =  m₀ /2 und m_B =  m₀ /2.

tmp

Das Modell erstellen wir mit der FEM-Formulierung für den Euler-Bernoulli-Balken mit zwei Elementen.

Header

Text

1+1

tmp

Zum späteren Plotten verwenden wir Referenzgrößen, die wir der analytischen Standard-Lösung eines Kragbalkens entnehmen.

Dies sind

${\begin{array}{ccl}{{W}_{\mathit {ref}}}&=&\displaystyle {\frac {{{q}_{0}}\cdot {{\ell }_{0}^{4}}}{8\cdot {{I}_{0}}\cdot E}},\\{{\Phi }_{\mathit {ref}}}&=&\displaystyle {\frac {{{q}_{0}}\cdot {{\ell }_{0}^{3}}}{6\cdot {{I}_{0}}\cdot E}},\\{{M}_{\mathit {ref}}}&=&\displaystyle {\frac {{{m}_{0}}\cdot {{\ell }_{0}}\cdot g}{2}},\\{{Q}_{\mathit {ref}}}&=&{{m}_{0}}\cdot g,\\{{q}_{0}}&=&\displaystyle {\frac {{{m}_{0}}\cdot g}{{\ell }_{0}}}\end{array}}$.

Wir verwenden dimensionslose Größen als Abkürzungen:

$\ell _{1}=\lambda _{1}\,\ell _{0},\ell _{B}=\lambda _{2}\,\ell _{0},\lambda _{1}=1-\lambda _{2},m_{A}=\theta _{A}\,m_{0},m_{B}=\theta _{B}\,m_{0}$.

Declarations

Text

1+1

tmp

Die analytische Lösung erzeugen wir aus Superpositionen der eingeprägten Lösungen zu den drei Gewichten auf die Struktur.

Dazu müssen wir ein bisschen umdenken: in den Standard-Lösungen ist der Balken jeweils links fest eingespannt - hier rechts.

Wir nehmen das komplementäre Koordinatensystem

${\bar {\xi }}=1-\xi$

und lesen die Lösung ab:

${\begin{array}{lccl}EI\,w({\bar {\xi }})&=&&\displaystyle {\frac {{q_{0}}\,\ell _{0}^{4}\cdot \left(6\cdot {{\bar {\xi }}^{2}}-4\cdot {{\bar {\xi }}^{3}}+{{\bar {\xi }}^{4}}\right)}{24}}\\&&+&\displaystyle {\frac {m_{A}\,g\,\ell _{0}^{3}\cdot \left(3\cdot {{\bar {\xi }}^{2}}-{{\bar {\xi }}^{3}}\right)}{6}}\\&&+&\left\{{\begin{array}{ll}\displaystyle {\frac {m_{B}\,g\,\ell _{0}^{3}\cdot \left(3\cdot {{\lambda }_{2}}\cdot {{\bar {\xi }}^{2}}-{{\bar {\xi }}^{3}}\right)}{6}}&{\text{ für }}{\bar {\xi }}\leq \lambda _{2}\\\displaystyle {\frac {m_{B}\,g\,\ell _{0}^{3}\cdot \left(3\cdot {{\lambda }_{2}}\cdot {{\bar {\xi }}^{2}}-{{\bar {\xi }}^{3}}+{{\left({\bar {\xi }}-{{\lambda }_{2}}\right)}^{3}}\right)}{6}}&{\text{ sonst}}\end{array}}\right.\end{array}}$.

Analytic Solution

Text

1+1

tmp

Die System-Matrix setzen wir aus den Element-Steifigkeitsmatrizen aus FEM: Trial-Functions für kubische Ansatz-Polynome zusammen und arbeiten die geometrischen Randbedingungen ein, indem wir die Zeilen und Spalten bzgl der Verschiebung und Verdrehung in C (W₂, Φ₂) streichen.

Das Gleichungssystem lautet dann

${\frac {E\,I}{\ell _{0}^{3}}}\cdot {\begin{pmatrix}-{\frac {12}{{{\lambda }_{2}^{3}}-3\cdot {{\lambda }_{2}^{2}}+3\cdot {{\lambda }_{2}}-1}}&{\frac {6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{2}}-2\cdot {{\lambda }_{2}}+1}}&{\frac {12}{{{\lambda }_{2}^{3}}-3\cdot {{\lambda }_{2}^{2}}+3\cdot {{\lambda }_{2}}-1}}&{\frac {6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{2}}-2\cdot {{\lambda }_{2}}+1}}\\{\frac {6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{2}}-2\cdot {{\lambda }_{2}}+1}}&-{\frac {4\cdot {{\ell }_{0}^{2}}}{{{\lambda }_{2}}-1}}&-{\frac {6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{2}}-2\cdot {{\lambda }_{2}}+1}}&-{\frac {2\cdot {{\ell }_{0}^{2}}}{{{\lambda }_{2}}-1}}\\{\frac {12}{{{\lambda }_{2}^{3}}-3\cdot {{\lambda }_{2}^{2}}+3\cdot {{\lambda }_{2}}-1}}&-{\frac {6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{2}}-2\cdot {{\lambda }_{2}}+1}}&-{\frac {12-36\cdot {{\lambda }_{2}}+36\cdot {{\lambda }_{2}^{2}}}{{{\lambda }_{2}^{6}}-3\cdot {{\lambda }_{2}^{5}}+3\cdot {{\lambda }_{2}^{4}}-{{\lambda }_{2}^{3}}}}&-{\frac {12\cdot {{\ell }_{0}}\cdot {{\lambda }_{2}}-6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{4}}-2\cdot {{\lambda }_{2}^{3}}+{{\lambda }_{2}^{2}}}}\\{\frac {6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{2}}-2\cdot {{\lambda }_{2}}+1}}&-{\frac {2\cdot {{\ell }_{0}^{2}}}{{{\lambda }_{2}}-1}}&-{\frac {12\cdot {{\ell }_{0}}\cdot {{\lambda }_{2}}-6\cdot {{\ell }_{0}}}{{{\lambda }_{2}^{4}}-2\cdot {{\lambda }_{2}^{3}}+{{\lambda }_{2}^{2}}}}&-{\frac {4\cdot {{\ell }_{0}^{2}}}{{{\lambda }_{2}^{2}}-{{\lambda }_{2}}}}\end{pmatrix}}\cdot {\underline {Q}}=m_{0}\,g\cdot {\begin{pmatrix}{\frac {1-{{\lambda }_{2}}+2\cdot {{\theta }_{A}}}{2}}\\{\frac {{{\ell }_{0}}-2\cdot {{\ell }_{0}}\cdot {{\lambda }_{2}}+{{\ell }_{0}}\cdot {{\lambda }_{2}^{2}}}{12}}\\{\frac {1+2\cdot {{\theta }_{B}}}{2}}\\{\frac {2\cdot {{\ell }_{0}}\cdot {{\lambda }_{2}}-{{\ell }_{0}}}{12}}\end{pmatrix}}$

mit

${\underline {Q}}=\left({\begin{array}{c}W_{0}\\\Phi _{0}\\W_{1}\\\Phi _{1}\\\end{array}}\right)$.

FE Formulation

Text

1+1

tmp

Lösen des linearen Gleichungssystems liefert

${\begin{array}{ccl}\displaystyle {\frac {W_{0}}{W_{ref}}}&=&\displaystyle -{\frac {-3-8\cdot {{\theta }_{A}}+\left(4\cdot {{\lambda }_{2}^{3}}-12\cdot {{\lambda }_{2}^{2}}\right)\cdot {{\theta }_{B}}}{3}}\\\displaystyle {\frac {\Phi _{0}}{W_{ref}}}&=&\displaystyle -{\frac {4+12\cdot {{\theta }_{A}}+12\cdot {{\lambda }_{2}^{2}}\cdot {{\theta }_{B}}}{3\cdot {{\ell }_{0}}}}\\\displaystyle {\frac {W_{1}}{W_{ref}}}&=&\displaystyle {\frac {6\cdot {{\lambda }_{2}^{2}}-4\cdot {{\lambda }_{2}^{3}}+{{\lambda }_{2}^{4}}+\left(12\cdot {{\lambda }_{2}^{2}}-4\cdot {{\lambda }_{2}^{3}}\right)\cdot {{\theta }_{A}}+8\cdot {{\lambda }_{2}^{3}}\cdot {{\theta }_{B}}}{3}}\\\displaystyle {\frac {\Phi _{1}}{W_{ref}}}&=&\displaystyle -{\frac {12\cdot {{\lambda }_{2}}-12\cdot {{\lambda }_{2}^{2}}+4\cdot {{\lambda }_{2}^{3}}+\left(24\cdot {{\lambda }_{2}}-12\cdot {{\lambda }_{2}^{2}}\right)\cdot {{\theta }_{A}}+12\cdot {{\lambda }_{2}^{2}}\cdot {{\theta }_{B}}}{3\cdot {{\ell }_{0}}}}\\\displaystyle {\frac {W_{2}}{W_{ref}}}&=&0\\\displaystyle {\frac {\Phi _{2}}{W_{ref}}}&=&0\end{array}}$.

Solving

Text

1+1

tmp

Und die Ergebnisse der analytischen und FE-Rechnungen plotten wir ...

... für w(x):

Biegelinie w(x)

... für Φ(x):

Kippung w'(x)

... für M(x):

Biegemoment M(x)

... für Q(x):

Querkraft Q(x)

Post-Processing

Text

1+1

---

Links

• ...

Literature

• ...