Gelöste Aufgaben/UEBD

Aufgabenstellung

Der Euler-Bernoulli-Balken AB wird durch eine vertikale Einzelkraft F belastet. Er ist in A fest eingespannt, in B gelenkig gelagert.

UEBC ist eine Variante dieser Aufgabe.

Lageplan

Gesucht ist eine Lösung für die Biegelinie w(x) mit dem Ansatz von Rayleigh-Ritz und zweiTrial-Funktion.

Lösung mit Maxima

Beim Verfahren von Ritz arbeiten wir mit

• dem Prinzip vom Minimum der Potentiellen Energie und
• Ansatzfunktionen über die gesamte Länge des Balkens.

Um die Lösung dimensionslos zu machen, nutzen wir die analytische Lösung des beidseits gelenkig gelagerten Balkens unter einer mittigen Einzellast (vgl. W8Zu). Hier ist die Neigung des Querschnitts am rechten Rand und die maximale Auslenkung jeweils

${\begin{array}{lll}\displaystyle w'(\ell )&=-\displaystyle {\frac {F\;\ell ^{2}}{16\;EI}}&=:{\hat {\varphi }}\\\displaystyle w({\frac {\ell }{2}})&=\displaystyle {\frac {F\;\ell ^{3}}{48\;EI}}&=:{\hat {w}}\end{array}}$.

Damit können wir uns die Lösungen dieses Problems als Vielfaches der analytischen Lösung eines ähnlichen Problems denken.

tmp

Das Vorgehen in dieser Aufgabe ist identisch zu UEBC - hier allerdings mit zwei Trial-Functions.

Header

Text

1+1

tmp

Zum Dimensionslos-Machen der Gleichungen brauchen wir später die Größen W, Φ, die wir hier bereits einführen:

${\begin{array}{ll}W&=\displaystyle {\frac {F\ell ^{3}}{48EI}}\cdot q_{w},\\\Phi &=-\displaystyle {\frac {F\ell ^{2}}{(6EI}}\cdot ({\frac {1}{2}}-{\frac {1}{2}}^{3})\cdot q_{\varphi },\\x&=\xi \cdot \ell ,\\a&=\alpha \cdot \ell .\end{array}}$.

Daraus folgen später die dimensionslosen Koordinaten des Problems zu

${\underline {Q}}=\left({\begin{array}{c}q_{w}\\q_{\varphi }\end{array}}\right)$.

Declarations

Text

1+1

tmp

Gesucht ist hier ein Ansatz mit zwei Trial Functions, also

$w(x)=\displaystyle \sum _{i=1}^{2}W_{i}\cdot \phi _{i}(x)$

und dem gesuchten Koeffizienten W_i. Die ϕ_i müssen dabei alle geometrischen Randbedingungen erfüllen:

${\begin{array}{ll}w(0)&=0,\\w'(0)&=0,\\w(l)&=0.\end{array}}$

Diese Funktion zu konstruieren ist nicht ganz trivial. Wir multiplizieren dafür einzelne Terme der Form (x-x₀), die die gewünschte Nullstelle für x=x₀ erzeugen. An der Stelle x=0 benötigen wir eine doppelte Nullstelle, damit auch die Ableitung an dieser Stelle Null ist.

Für die beiden Trial-Funktionen wählen wir

• ϕ₁ so, dass sie am rechten Rand keine Steigung hat, aber eine Auslenkung am Kraftangriffspunkt, also $w_{1}(\xi )=C_{W}\cdot \left(\xi -0\right)^{2}\cdot \left(\xi -1\right)^{2}{\text{mit }}\xi =\displaystyle {\frac {x}{\ell }}$ und
• ϕ₂ so, dass sie am rechten Rand eine Steigung hat, aber keine Auslenkung am Kraftangriffspunkt, also $w_{2}(\xi )=C_{\Phi }\cdot \left(\xi -0\right)^{2}\cdot \left(\xi -1\right)\cdot \left(\xi -\alpha \right){\text{mit }}\xi =\displaystyle {\frac {x}{\ell }}$.

Die freien Koeffizienten C_i wählen wir so, dass wir sie anschaulich "denken" können - wir wählen

• die gesuchte Auslenkung an der Kraft-Einleitungsstelle x=a für 'ϕ₁', also

$\displaystyle w(a)={\hat {w}}\cdot q_{w}\;\;{\text{ mit der Referenz-Lösung }}\;\;{\hat {w}}={\frac {F\;\ell ^{3}}{48\;EI}}$ und

• die gesuchte Neigung am rechten Rand für ϕ₂ , also

$\displaystyle w'(\ell )={\hat {\varphi }}\cdot q_{\varphi }\;\;{\text{ mit der Referenz-Lösung }}\;\;{\hat {\varphi }}=-{\frac {F\;\ell ^{2}}{16\;EI}}$.

aus W8Zu. Mit den neuen, dimensionslosen Koordinaten Q ist

$\displaystyle C_{W}={\frac {\hat {w}}{{{\alpha }^{4}}-2\cdot {{\alpha }^{3}}+{{\alpha }^{2}}}}\cdot q_{w}$,

$\displaystyle C_{\Phi }=-{\frac {\hat {\varphi }}{\alpha -1}}\cdot q_{\varphi }\cdot \ell$

und damit

$\displaystyle w(x)=\underbrace {{\frac {F\ell ^{3}}{48EI}}\cdot q_{w}} _{\displaystyle =W}\cdot \underbrace {\frac {\xi ^{2}(\xi -1)^{2}}{\alpha ^{2}(1-\alpha )^{2}}} _{\displaystyle =\phi _{1}(\xi )}\underbrace {-{\frac {F\ell ^{2}}{16EI}}\cdot q_{\varphi }\ell } _{\displaystyle =\Phi }\cdot \underbrace {\frac {\xi ^{2}(\xi -1)(\xi -\alpha )}{(1-\alpha )}} _{\displaystyle =\phi _{2}(\xi )}$.

Trial-Functions

Für α=2/3 sehen die Trial-Funktionen so aus:

Formfunctions

Text

1+1

tmp

Für die Gleichgewichtsbedingungen setzten wir Π und A in U ein und schreiben die skalare Gleichung allgemein in Matrizenform an. Dabei müssen wir

$\displaystyle {\frac {d\phi }{x}}={\frac {d\phi }{\xi }}\cdot \underbrace {\displaystyle {\frac {d\xi }{x}}} _{\displaystyle ={\frac {1}{\ell }}}$

berücksichtigen und erhalten

$U=\displaystyle {\frac {1}{2}}\cdot \displaystyle {\underline {Q}}^{T}\cdot {\underline {\underline {A}}}\cdot {\underline {Q}}-{\underline {Q}}^{T}\cdot {\underline {b}}$.

Potential Energy

Text

1+1

tmp

Diese Gleichung erfüllt die Gleichgewichtsbedingungen

$\displaystyle {\frac {dU}{dW}}{\stackrel {!}{=}}0$,

wenn

${\underline {\underline {A}}}\cdot {\underline {Q}}={\underline {b}}$.

Das Gleichungssystem dazu ist

${\begin{pmatrix}{\frac {1}{60\cdot {{\alpha }^{8}}-240\cdot {{\alpha }^{7}}+360\cdot {{\alpha }^{6}}-240\cdot {{\alpha }^{5}}+60\cdot {{\alpha }^{4}}}}&{\frac {1}{20\cdot {{\alpha }^{5}}-60\cdot {{\alpha }^{4}}+60\cdot {{\alpha }^{3}}-20\cdot {{\alpha }^{2}}}}\\-{\frac {\ell }{60\cdot {{\alpha }^{5}}-180\cdot {{\alpha }^{4}}+180\cdot {{\alpha }^{3}}-60\cdot {{\alpha }^{2}}}}&-{\frac {\left(6-10\cdot \alpha +5\cdot {{\alpha }^{2}}\right)\cdot \ell }{20\cdot {{\alpha }^{2}}-40\cdot \alpha +20}}\end{pmatrix}}\cdot {\begin{pmatrix}q_{w}\\q_{\varphi }\end{pmatrix}}={\begin{pmatrix}-1\\0\end{pmatrix}}$

Equilibrium Conditions

Text

1+1

tmp

Diese Gleichgewichtsbedingung ist erfüllt , wenn

${\begin{array}{ll}{{q}_{w}}&=60\cdot {{\alpha }^{8}}-240\cdot {{\alpha }^{7}}+372\cdot {{\alpha }^{6}}-264\cdot {{\alpha }^{5}}+72\cdot {{\alpha }^{4}},\\{{q}_{\phi }}&=4\,{{\alpha }^{2}}\cdot (1-{\alpha })\end{array}}$.

Solving

Text

1+1

tmp

Auslenkung und Kippung w(x), ϕ(x) der Querschnitte

Wir schauen uns zunächst die Absenkung des Kraftangriffspunktes w(a) und die Neigung w'(l) für verschiedene Werte von a (α)  an:

Wie der Verlauf der approximierten Biegelinie aussieht, schauen wir uns jetzt an - und zwar für verschiedene Werte von α:

Parameterstudie für α

Sie sehen, dass für α=6/10 die Auslenkung dieses Systems nur knapp die Hälfte der maximalen Auslenkung eines beidseitig gelenkig gelagerten Balkens beträgt - hier in unserer Approximation mit zwei Trial-Funktion.

Aufgabe UEBC zeigt, welche Ergebnisse die Approximation mit einer Trial-Funktion liefert.

✔ {{{title}}}:

Probieren Sie aus, wie gut der Ritz-Ansatz den Verlauf des Schnittmoments approximiert!

Post-Processing

Text

1+1

---

Links

• ...

Literature

• ...