Gelöste Aufgaben/UEBO

Aufgabenstellung

Diese Problemstellung liefert einen Näherungsansatz für eine Standardlösung zum Euler-Bernoulli-Balken.

Der Euler-Bernoulli-Balken AB wird durch ein Moment M zwischen den beiden gelenkigen Lagern belastet. 

Lageplan

Gesucht ist eine Lösung für die Biegelinie mit dem Ansatz von Ritz und zwei Trial-Funktionen.

(Weg "1" wie in UEBH beschrieben.)

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.

tmp

Wir berechnen die Potentielle Energie U des Systems in Abhängigkeit von den generalisierten Koordinaten W_i und erhalten aus

$\displaystyle {\frac {d\,U}{d\,W_{i}}}{\stackrel {!}{=}}0$

die Gleichung für den gesuchten Koeffizienten W_i der Trial-Funktionen.

Header

Text

1+1

tmp

Um die Lösung dimensionslos zu machen, nutzen wir die analytische Lösung des Problems , hier die Beträge der maximalen Auslenkung des Balkens für a = ℓ und der Verdrehung des Balkens am Momenten-Angriffspunkt für a = ℓ/2:

$W^{max}=\displaystyle {\frac {M\,\ell ^{2}}{9{\sqrt {3}}\,\cdot E\,I}}$

die maximale Auslenkung des Balkens für a=ℓ

$\Phi ^{M}=\displaystyle {\frac {M\,\ell }{12\cdot E\,I}}$ die Verdrehung des Balkens am Momenten-Angriffspunkt für a=ℓ/2

Dimensionslose Orts-Koordinaten sind

${\begin{array}{ll}x&=\xi \cdot \ell ,\\a&=\alpha \cdot \ell .\end{array}}$.

Declarations

Text

1+1

tmp

Bei der Suche nach passenden Trial-Functions ϕ lassen wir uns ebenfalls von der analytischen Lösung des Problems "inspirieren":

$w_{a}=\displaystyle {\frac {M\ell ^{2}}{3EI}}\left(\left(3\cdot {{\alpha }^{2}}-6\cdot \alpha +2\right)\cdot \xi +{{\xi }^{3}}-3<\xi -\alpha >^{2}\right)$

Der Funktionsverlauf von w^a hat zwei charakteristische Ausprägungen:

... für a=0	... für a=ℓ/2
Diese Lösung - mit dem angreifenden Moment in A - hat eine starke symmetrische Komponente bzgl der Stab-Mitte.	Diese Lösung - mit dem angreifenden Moment in der Stab-Mitte - ist punktsymmetrisch zum Momenten-Angriffspunkt.

Und so wählen wir unsere Trial-Functions als

$\phi _{1}=c_{1}\cdot \xi \cdot (1-\xi ){\text{ und }}\phi _{2}=c_{2}\cdot \xi \cdot (1-\xi )\cdot (\alpha -\xi )$.

Trial-Functions

Für α=7∙ℓ/10 sehen sie so aus;

Die Koeffizienten c₁ und c₂ haben wir dabei so gewählt, dass

${\begin{array}{ll}\phi _{1}(\displaystyle {\frac {1}{2}})&=1\\\phi _{2}'(\alpha )&=1\end{array}}$.

Mit den neuen, gesuchten Wichtungsfaktoren q_w und q_ϕ ist die Ansatzfunktion zur Lösung mit dem Verfahren von Rayleigh-Ritz damit

$\displaystyle w(\xi )={\frac {M\ell ^{2}}{3EI}}\left(q_{w}\left(4{\frac {\xi -\xi ^{2}}{3{\sqrt {3}}}}\right)+q_{\phi }\left(-\xi ^{3}+(1+\alpha )\xi ^{2}-\alpha \xi \right)\right)$

Aufgrund der gewählten Skalierungsfaktoren erwarten wir als Ergebnis näherungsweise

• für α=½: q_w ≈ 0 und q_ϕ ≈ 1,
• für α= 0: q_w ≈ 1 und q_ϕ ≈ 0.

Formfunctions

Text

1+1

tmp

Für die Gleichgewichtsbedingungen setzten wir Π (aus Abschnitt Euler-Bernoulli-Balken) 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 mit der Arbeitsfunktion des Moments

$A=M\cdot w'|_{\displaystyle \xi =\alpha }$

das Potential in Matrix-Schreibweise:

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

wobei

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

Einsetzen der Ansatzfunktion in die Formänderungsenergie und die Arbeitsfunktion liefert für die Matrizen A und b:

${\underline {\underline {A}}}={\begin{pmatrix}\displaystyle {\frac {128}{81}}&\displaystyle {\frac {16}{{3}^{\frac {5}{2}}}}-{\frac {32\cdot \alpha }{{3}^{\frac {5}{2}}}}\\\displaystyle {\frac {16}{{3}^{\frac {5}{2}}}}-{\frac {32\cdot \alpha }{{3}^{\frac {5}{2}}}}&\displaystyle {\frac {8\cdot {{\alpha }^{2}}}{3}}-{\frac {8\cdot \alpha }{3}}+{\frac {8}{3}}\end{pmatrix}}$,

${\underline {b}}={\begin{pmatrix}\displaystyle {\frac {8}{{3}^{\frac {3}{2}}}}-{\frac {16\cdot \alpha }{{3}^{\frac {3}{2}}}}\\\displaystyle 2\cdot \alpha -2\cdot {{\alpha }^{2}}\end{pmatrix}}$.

Potential Energy

Text

1+1

tmp

Diese Gleichung erfüllt die Gleichgewichtsbedingungen

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

wenn

${\begin{pmatrix}\displaystyle {\frac {128}{81}}&\displaystyle {\frac {16}{{3}^{\frac {5}{2}}}}-{\frac {32\cdot \alpha }{{3}^{\frac {5}{2}}}}\\\displaystyle {\frac {16}{{3}^{\frac {5}{2}}}}-{\frac {32\cdot \alpha }{{3}^{\frac {5}{2}}}}&\displaystyle {\frac {8\cdot {{\alpha }^{2}}}{3}}-{\frac {8\cdot \alpha }{3}}+{\frac {8}{3}}\end{pmatrix}}\cdot {\begin{pmatrix}{{q}_{w}}\\{{q}_{p}}\end{pmatrix}}={\begin{pmatrix}\displaystyle {\frac {8}{{3}^{\frac {3}{2}}}}-{\frac {16\cdot \alpha }{{3}^{\frac {3}{2}}}}\\\displaystyle 2\cdot \alpha -2\cdot {{\alpha }^{2}}\end{pmatrix}}$.

Equilibrium Conditions

Text

1+1

tmp

Auflösen der Gleichungen nach den unbekannten Koordinaten q_w und q_ϕ liefert

${\begin{array}{ll}q_{w}&=-{\frac {3^{\frac {3}{2}}}{8}}\left(-2+7\alpha -9\alpha ^{2}+6\alpha ^{3}\right)\\q_{\phi }&=-{\frac {1}{2}}\left(1-6\alpha +6\alpha ^{2}\right)\end{array}}$.

Damit ist die gesuchte Näherungs-Lösung

$\displaystyle w(\xi )={\frac {M\,\ell ^{2}}{6\,EI}}\left(\left(3\cdot {{\alpha }^{2}}-6\cdot \alpha +2\right)\cdot \xi +\left(-9\cdot {{\alpha }^{2}}+12\cdot \alpha -3\right)\cdot {{\xi }^{2}}+\left(6\cdot {{\alpha }^{2}}-6\cdot \alpha +1\right)\cdot {{\xi }^{3}}\right)$.

Solving

Text

1+1

tmp

Verlauf der Koordinaten q_w, q_ϕ

Die gesuchten Koordinaten  q_w und q_Φ sind dimensionslos. Wir können sie direkt für verschiedene Werte von α auftragen.

Wir sehen:

• für α=½: die Lösung wird - wie erwartet - nur durch ϕ₂ beschreiben - also q_w ≈ 0 und q_ϕ ≈ 1; allerdings ist die Qualität der Lösung mit q_ϕ = 1/4 sehr schlecht - hier drückt der Sprung in der Momenten-Kennlinie der analytischen Lösung auf das Ergebnis (s.u.).
• für α= 0: die Lösung wird - wie erwartet - primär durch ϕ₁ beschreiben, also q_w ≈ 1 und q_ϕ ≈ 0. Hier zeigt die Lösung mit q_w = 1.3 und q_ϕ = -0.5 einen recht großen Lösungs-Anteil der punktsymmetrischen Trial-Function.

Parameterstudie Biegelinie

Und so sieht die normierte Biegelinie des Balkens im Vergleich von Ritz-Näherung zu analytischer Lösung für verschiedene Werte von a aus:

Die dicken Linien gehören zu Näherung nach dem Ritz-Ansatz, die dünnen zur analytischen Lösung. Je weiter der Momenten-Angriffspunkt in die Balken-Mitte rückt und besonders für α=1/2 liefert der Ritz-Ansatz kein überzeugendes Ergebnis. Hier müssten wir mehr Trial-Functions "spendieren".

Vergleich: analytische / numerische Lösung für den Biegemomenten-Verlauf.

Post-Processing

Text

1+1

Wieso die Näherungslösung - besonders für α=½ - so schlecht ist, erkennt man beim Auftragen der Biegemomente im Stab für

• die analytische Lösung

$\displaystyle {\frac {M(x)}{M}}=\left\{{\begin{array}{cl}-\xi &\ldots {\text{ für }}\xi <{\frac {1}{2}}\\1-\xi &\ldots {\text{ sonst. }}\end{array}}\right.$ und

• die numerische Lösung

$\displaystyle {\frac {M(x)}{M}}={\frac {1}{4}}(2\xi -1)$.

Post-Processing - Nachtrag

Text

1+1

---

Links

• ...

Literature

• ...