Gelöste Aufgaben/UEBF: Unterschied zwischen den Versionen

Aus numpedia
Zur Navigation springen Zur Suche springen
Keine Bearbeitungszusammenfassung
Keine Bearbeitungszusammenfassung
Zeile 26: Zeile 26:
Um die Lösung dimensionslos zu machen, nutzen wir die [[Sources/Lexikon/Standard-Lösungen|analytische Lösung]] des einseitig fest eingespannten Balkens  unter einer Einzellast bei ''B''. Hier ist die maximale Auslenkung  
Um die Lösung dimensionslos zu machen, nutzen wir die [[Sources/Lexikon/Standard-Lösungen|analytische Lösung]] des einseitig fest eingespannten Balkens  unter einer Einzellast bei ''B''. Hier ist die maximale Auslenkung  


.
<math>\displaystyle \hat{w} := \frac{\ell^3}{3\; E\,I_1}\,F_0</math>.


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


<!-------------------------------------------------------------------------------->
Wir berechnen die Potentielle Energie '''''U''''' des Systems und erhalten aus
 
<math>\displaystyle \frac{d\,U}{d\,W} \stackrel{!}{=} 0 </math>
 
die Gleichung für den gesuchten Koeffizienten ''W'' der Trial-Funktion.<!-------------------------------------------------------------------------------->
 
{{MyCodeBlock|title=Header
{{MyCodeBlock|title=Header
|text=Text
|text=Text
Zeile 43: Zeile 48:
==tmp==
==tmp==


<!-------------------------------------------------------------------------------->
Zum Dimensionslos-Machen der Gleichungen wählen wir
 
<math>\begin{array}{ll}\hat{w} &=\displaystyle  \frac{\ell^3}{3 E\,I_1}\cdot F,\\  x&=\xi\cdot \ell,\\ a&=\alpha\cdot \ell.\end{array}</math>.<!-------------------------------------------------------------------------------->
 
{{MyCodeBlock|title=Declarations
{{MyCodeBlock|title=Declarations
|text=Text
|text=Text
Zeile 54: Zeile 62:
==tmp==
==tmp==


<!-------------------------------------------------------------------------------->
Gesucht ist hier ein Ansatz mit einer einzigen Trial Function, also
 
<math>w(x) = W\cdot\phi(x)</math>
 
und dem gesuchten Koeffizienten ''W''. ''ϕ'' muss dabei alle geometrischen Randbedingungen erfüllen:
 
<math>\begin{array}{ll}w(0) &= 0,\\w'(0) &= 0.\end{array}</math>
 
Die Funktionen müssen nicht die Randbedingungen für Schnittlasten in ''B'' erfüllen - schließlich suchen wir nach einer Näherungslösung!
 
Diese Funktion ''ϕ'' zu konstruieren ist einfach: wir brauchen lediglich eine doppelte Nullstelle in ''x=0'', also z.B.
 
<math>\phi(\xi) = \xi^2 \text{ mit } \xi = \displaystyle\frac{x}{\ell}</math>.
 
Der freie Koeffizient ''W'' ist nun die Auslenkung in ''B'' - deshalb taufen wir ihn hier ''W<sub>B</sub>''. Wir wählen ihn so, dass Sie die exakte Lösung bestmöglich approximiert.<!-------------------------------------------------------------------------------->
 
{{MyCodeBlock|title=Formfunctions
{{MyCodeBlock|title=Formfunctions
|text=Text
|text=Text
Zeile 65: Zeile 88:
==tmp==
==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
 
<math>\displaystyle \frac{d\phi}{x} = \frac{d\phi}{\xi}\cdot\underbrace{\displaystyle\frac{d\xi}{x}}_{\displaystyle = \frac{1}{\ell}}</math>
 
berücksichtigen und erhalten
 
<math>U = \displaystyle \frac{1}{2} \cdot \displaystyle \underline{Q}^T \cdot \underline{\underline{A}}\cdot \underline{Q} - \underline{Q}^T\cdot \underline{b} </math>.
 
Weil wir in der Spaltenmatrix
 
<math>\underline{Q} = \left(W_B\right)</math>
 
nur eine Unbekannte ''W<sub>B</sub>'' haben, können wir es uns hier etwas einfacher machen, ''U'' lautet:
 
<math>\displaystyle U=\frac{E\cdot \left( {{I}_{0}}+{{I}_{1}}\right)}{{{\ell}^{3}}} \cdot {{W}_{B}^{2}}-\frac{{{a}^{2}}\cdot F}{{{\ell}^{2}}}\cdot {{W}_{B}}</math>.<!-------------------------------------------------------------------------------->
 
 
 
{{MyCodeBlock|title=Potential Energy
{{MyCodeBlock|title=Potential Energy
|text=Text
|text=Text
Zeile 76: Zeile 116:
==tmp==
==tmp==


<!-------------------------------------------------------------------------------->
Diese Gleichung erfüllt die Gleichgewichtsbedingung(en)
 
<math>\displaystyle \frac{dU}{dW_B} \stackrel{!}{=} 0</math>,
 
wenn
 
<math>\displaystyle \frac{2\cdot \left( {{I}_{0}}+{{I}_{1}}\right) \cdot {{W}_{B}}\cdot E}{{{\ell}^{3}}}-\frac{{{a}^{2}}\cdot F}{{{\ell}^{2}}}=0</math>.<!-------------------------------------------------------------------------------->
 
 
 
{{MyCodeBlock|title=Equilibrium Conditions
{{MyCodeBlock|title=Equilibrium Conditions
|text=Text
|text=Text
Zeile 87: Zeile 136:
==tmp==
==tmp==


<!-------------------------------------------------------------------------------->
Auflösen nach ''W<sub>B</sub>'' liefert
 
<math>\displaystyle W_B = \frac{{{a}^{2}}\cdot \ell}{2\cdot \left( {{I}_{1}}+{{I}_{0}}\right) \cdot E}\cdot F</math>.<!-------------------------------------------------------------------------------->
 
 
{{MyCodeBlock|title=Solving
{{MyCodeBlock|title=Solving
|text=Text
|text=Text
Zeile 97: Zeile 150:


==tmp==
==tmp==
[[Datei:UEBF-21.png|mini|Lageplan]]Und so sieht die normierte Absenkung des Punktes ''B'' als Funktion von ''a'' aus:
Dies ist nicht die Biegelinie des Balkens!
Aufgabe [[Gelöste Aufgaben/UEBD|UEBD]] zeigt, wie wir zu einer besseren Approximation kommen.
{{MyTip|title=Schnittmomenten-Verlauf|text=Probieren Sie aus, wie gut der Ritz-Ansatz den Verlauf des Schnittmoments approximiert!}}<!-------------------------------------------------------------------------------->


<!-------------------------------------------------------------------------------->
{{MyCodeBlock|title=Post-Processing
{{MyCodeBlock|title=Post-Processing
|text=Text
|text=Text
Zeile 114: Zeile 173:




<hr/>
<hr />
'''Links'''
'''Links'''
* ...
* ...
Zeile 120: Zeile 179:
'''Literature'''
'''Literature'''
* ...
* ...
[[Datei:UEBF-21.png|mini|Lageplan]]

Version vom 14. April 2021, 11:24 Uhr


Aufgabenstellung

Der Euler-Bernoulli-Balken AB wird durch eine vertikale Einzelkraft F0 belastet. Er ist in A fest eingespannt und hat ein linear veränderliches Flächenmoment I(x).


Lageplan

Gesucht ist eine Lösung für die Auslenkung w(x) mit dem Ansatz von Ritz und einer Trial-Funktion.

Damit ähnelt diese Aufgabe UEBC.

Lösung mit Maxima

Beim Verfahren von Ritz arbeiten wir mit

Um die Lösung dimensionslos zu machen, nutzen wir die analytische Lösung des einseitig fest eingespannten Balkens  unter einer Einzellast bei B. Hier ist die maximale Auslenkung

.

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

tmp

Wir berechnen die Potentielle Energie U des Systems und erhalten aus

die Gleichung für den gesuchten Koeffizienten W der Trial-Funktion.

Header

Text


1+1




tmp

Zum Dimensionslos-Machen der Gleichungen wählen wir

.

Declarations

Text


1+1




tmp

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

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

Die Funktionen müssen nicht die Randbedingungen für Schnittlasten in B erfüllen - schließlich suchen wir nach einer Näherungslösung!

Diese Funktion ϕ zu konstruieren ist einfach: wir brauchen lediglich eine doppelte Nullstelle in x=0, also z.B.

.

Der freie Koeffizient W ist nun die Auslenkung in B - deshalb taufen wir ihn hier WB. Wir wählen ihn so, dass Sie die exakte Lösung bestmöglich approximiert.

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

berücksichtigen und erhalten

.

Weil wir in der Spaltenmatrix

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

.


Potential Energy

Text


1+1




tmp

Diese Gleichung erfüllt die Gleichgewichtsbedingung(en)

,

wenn

.


Equilibrium Conditions

Text


1+1




tmp

Auflösen nach WB liefert

.


Solving

Text


1+1




tmp

Lageplan

Und so sieht die normierte Absenkung des Punktes B als Funktion von a aus:

Dies ist nicht die Biegelinie des Balkens!

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

  • ...