Gelöste Aufgaben/UEBC

Aufgabenstellung

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

Lageplan

Gesucht ist eine Lösung für die Biegelinie w(x) mit dem Ansatz von Rayleigh-Ritz und einer Trial-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 maximale Auslenkung

$\displaystyle w^{*}({\frac {\ell }{2}})={\frac {F\;\ell ^{3}}{48\;EI}}=:{\hat {w}}$.

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

tmp

Wir lösen die Aufgabe mit Computer-Algebra - mit Maxima.

Header

Text

1+1

tmp

Zum Dimensionslos-Machen der Gleichungen wählen wir

${\begin{array}{ll}W&=\displaystyle {\frac {F\ell ^{3}}{48EI}}\cdot q,\\x&=\xi \cdot \ell ,\\a&=\alpha \cdot \ell .\end{array}}$.

Declarations

Text

1+1

tmp

Gesucht ist hier ein Ansatz mit einer einzigen Trial Function, also

$w(x)=W\cdot \phi (x)$

und dem gesuchten Koeffizienten W. ϕ muss dabei alle geometrischen Randbedingungen erfüllen:

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

Die Funktionen müssen nicht die Randbedingungen für das Moment an den Rändern A und B erfüllen!

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, also

$w(\xi )=C\cdot \left(\xi -0\right)^{2}\cdot \left(\xi -1\right){\text{mit }}\xi =\displaystyle {\frac {x}{\ell }}$.

Den freien Koeffizienten C wählen wir so, dass wir ihn anschaulich "denken" können - wir wählen die gesuchte Auslenkung an der Kraft-Einleitungsstelle x=a aus W8Zu, also

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

Mit der neuen, dimensionslosen Koordinate q ist

$\displaystyle C={\frac {q\cdot {\hat {w}}}{\alpha ^{3}-\alpha }}$,

und damit

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

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}}$.

Weil wir in

${\underline {Q}}=\left(q\right)$

nur eine Unbekannte q haben, können wir es uns hier etwas einfacher machen, U lautet:

$\displaystyle U={\frac {\left(\left(-24\cdot {{\alpha }^{6}}+48\cdot {{\alpha }^{5}}-24\cdot {{\alpha }^{4}}\right)\cdot {{\ell }^{3}}\cdot q+{{\ell }^{3}}\cdot {{q}^{2}}\right)\cdot {{F}^{2}}}{\left(1152\cdot {{\alpha }^{6}}-2304\cdot {{\alpha }^{5}}+1152\cdot {{\alpha }^{4}}\right)\cdot {\mathit {EI}}}}$.

Potential Energy

Text

1+1

tmp

Diese Gleichung erfüllt die Gleichgewichtsbedingungen

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

wenn

$\displaystyle {\frac {\left(-12\cdot {{\alpha }^{4}}+24\cdot {{\alpha }^{5}}-12\cdot {{\alpha }^{6}}+q\right)\cdot F}{12\cdot {{\alpha }^{6}}-24\cdot {{\alpha }^{5}}+12\cdot {{\alpha }^{4}}}}=0$.

Equilibrium Conditions

Text

1+1

tmp

Diese Gleichgewichtsbedingung ist erfüllt , wenn

$q=12\cdot {{\alpha }^{6}}-24\cdot {{\alpha }^{5}}+12\cdot {{\alpha }^{4}}$.

Solving

Text

1+1

tmp

Auslenkung w(x)

Wir schauen uns zunächst die Absenkung des Kraftangriffspunktes w(a) für verschiedene Werte von a (α)  an: Die maximale Auslenkung des Kraftangriffspunktes ist hier bei α=2/3 erreicht. Das ist nicht identisch mit der größten Absenkung der Querschnitte im Gebiet!

Parameter-Studie für α

Dafür schauen wir uns den Verlauf der approximierten Biegelinie an - und zwar für verschiedene Werte von α:

Sie sehen, dass für α=6/10 die Auslenkung dieses Systems nur ein Viertel der maximalen Auslenkung eines beidseitig gelenkig gelagerten Balkens beträgt - jedenfalls in unserer Approximation mit einer Trial-Funktion.

Aufgabe UEBD zeigt, wie wir zu einer besseren Approximation kommen.

✔ Schnittmomenten-Verlauf:

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

Post-Processing

Text

1+1

---

Links

• ...

Literature

• ...