Gelöste Aufgaben/TC12: Unterschied zwischen den Versionen

Aus numpedia
Zur Navigation springen Zur Suche springen
Zeile 106: Zeile 106:
====Aus Rand "A"====
====Aus Rand "A"====
<table class="wikitable" style="background-color:white;">
<table class="wikitable" style="background-color:white;">
<tr><td>[[Datei:TC12-11A.png|rahmenlos|alternativtext=|150px]]
<tr><td>[[Datei:TC12-11A.png|rahmenlos|alternativtext=|90x90px]]
</td><td>''Geometrische Randbedingungen''
</td><td>''Geometrische Randbedingungen''
 
# <math>w_1(0)=0</math>
* keine


''Kraft- und Momenten-Randbedingungen''
''Kraft- und Momenten-Randbedingungen''


# <math>m_A\;g - S \sin(\alpha_A) + Q_{A,+}= 0 \text{ mit } Q_{A,+} = - E I\cdot w'''(x)|_{x=0}</math>
#<math>M_{A} = 0 \text{ mit } M_{A} = - EI_1\cdot w''(x)|_{x=0}</math>
# <math>M_{A,+}= 0 \text{ mit } M_{A,+} = - E I\cdot w'''(x)|_{x=0}</math>
</td></tr>
</td></tr>
</table>
</table>
Zeile 120: Zeile 118:
====Aus Übergang "B"====
====Aus Übergang "B"====
<table class="wikitable" style="background-color:white;">
<table class="wikitable" style="background-color:white;">
<tr><td>[[Datei:TC12-11B.png|rahmenlos|alternativtext=|180px]]
<tr><td>[[Datei:TC12-11B.png|rahmenlos|alternativtext=|120x120px]]
</td><td>''Geometrische Randbedingungen''
</td><td>''Geometrische Randbedingungen''
# <math>w_1(\ell_1)=w_2(\ell_1)</math>
# <math>w_1(\ell_1)=w_2(\ell_1)</math>
Zeile 127: Zeile 125:
''Kraft- und Momenten-Randbedingungen''
''Kraft- und Momenten-Randbedingungen''


# <math>-M_{B,-} + M_{B,+} = 0</math>
#<math>-M_{B,-} - M_{B0} + M_{B,+} = 0</math>
# <math>-Q_{B,-} - S\; \sin(\alpha_B) +Q_{B,+} = 0</math>
#<math>-Q_{B,-}-k_B\cdot w_B + -Q_{B,+} = 0</math>
</td></tr>
</td></tr>
</table>
</table>
Zeile 134: Zeile 132:
====Aus Rand "C"====
====Aus Rand "C"====
<table class="wikitable" style="background-color:white;">
<table class="wikitable" style="background-color:white;">
<tr><td>[[Datei:TC12-11C1.png|rahmenlos|alternativtext=|170px]]<br/>
<tr><td>[[Datei:TC12-11C1.png|rahmenlos|alternativtext=|125x125px]]<br/>
[[Datei:TC12-11C2.png|rahmenlos|alternativtext=|170px]]
[[Datei:TC12-11C2.png|rahmenlos|alternativtext=|130x130px]]
</td><td>''Geometrische Randbedingungen''
</td><td>''Geometrische Randbedingungen''
# <math>w_2(\ell)=0</math>
#<math>w_C = r\cdot\Phi_C \text{ mit } \Phi_C = - \phi_2(l_2)</math>


''Kraft- und Momenten-Randbedingungen''
''Kraft- und Momenten-Randbedingungen''


# <math>K_C\cdot \Phi_C - M_{C,-} = 0 \text{ mit } M_{C,-} = - E I\cdot w''(x)|_{x=\ell}</math>
#<math>r\cdot Q_{C} + M_C = 0</math>
</td></tr>
</td></tr>
</table>
</table>
Zeile 147: Zeile 145:
Einsetzen liefert 8 Gleichungen
Einsetzen liefert 8 Gleichungen


<math>\begin{array}{rl} \displaystyle \frac{{{C}_{1,0}}}{{{\mathit{EI}}_{1}}}&=0\cr \displaystyle  -{{C}_{1,2}}&=0\cr \displaystyle  \frac{{{l}_{1}^{3}}\cdot {{C}_{1,3}}}{6\cdot {{\mathit{EI}}_{1}}}+\frac{{{l}_{1}^{2}}\cdot {{C}_{1,2}}}{2\cdot {{\mathit{EI}}_{1}}}+\frac{{{l}_{1}}\cdot {{C}_{1,1}}}{{{\mathit{EI}}_{1}}}+\frac{{{C}_{1,0}}}{{{\mathit{EI}}_{1}}}+\frac{{{q}_{0}}\cdot {{l}_{1}^{4}}}{24\cdot {{\mathit{EI}}_{1}}}&=\frac{{{C}_{2,0}}}{{{\mathit{EI}}_{2}}}\cr \displaystyle  \frac{{{l}_{1}^{2}}\cdot {{C}_{1,3}}}{2\cdot {{\mathit{EI}}_{1}}}+\frac{{{l}_{1}}\cdot {{C}_{1,2}}}{{{\mathit{EI}}_{1}}}+\frac{{{C}_{1,1}}}{{{\mathit{EI}}_{1}}}+\frac{{{q}_{0}}\cdot {{l}_{1}^{3}}}{6\cdot {{\mathit{EI}}_{1}}}&=\frac{{{C}_{2,1}}}{{{\mathit{EI}}_{2}}}\cr \displaystyle  -\frac{{{C}_{2,0}}\cdot {{k}_{B}}}{{{\mathit{EI}}_{2}}}-{{C}_{2,3}}+{{C}_{1,3}}+{{q}_{0}}\cdot {{l}_{1}}&=0\cr \displaystyle  -{{C}_{2,2}}+{{l}_{1}}\cdot {{C}_{1,3}}+{{C}_{1,2}}+\frac{{{q}_{0}}\cdot {{l}_{1}^{2}}}{2}&={{M}_{B}}\cr \displaystyle  -\frac{{{l}_{2}^{2}}\cdot {{C}_{2,3}}\cdot r}{2\cdot {{\mathit{EI}}_{2}}}-\frac{{{l}_{2}}\cdot {{C}_{2,2}}\cdot r}{{{\mathit{EI}}_{2}}}-\frac{{{C}_{2,1}}\cdot r}{{{\mathit{EI}}_{2}}}&=\displaystyle\frac{{{l}_{2}^{3}}\cdot {{C}_{2,3}}}{6\cdot {{\mathit{EI}}_{2}}}+\frac{{{l}_{2}^{2}}\cdot {{C}_{2,2}}}{2\cdot {{\mathit{EI}}_{2}}}+\frac{{{l}_{2}}\cdot {{C}_{2,1}}}{{{\mathit{EI}}_{2}}}+\frac{{{C}_{2,0}}}{{{\mathit{EI}}_{2}}}\cr \displaystyle  -{{C}_{2,3}}\cdot r-{{l}_{2}}\cdot {{C}_{2,3}}-{{C}_{2,2}}&=0 \end{array}</math>
<math>\begin{array}{rl} \displaystyle \frac{{{C}_{1,0}}}{{{\mathit{EI}}_{1}}}&=0\\ \displaystyle  -{{C}_{1,2}}&=0\\ \displaystyle  \frac{{{l}_{1}^{3}}\cdot {{C}_{1,3}}}{6\cdot {{\mathit{EI}}_{1}}}+\frac{{{l}_{1}^{2}}\cdot {{C}_{1,2}}}{2\cdot {{\mathit{EI}}_{1}}}+\frac{{{l}_{1}}\cdot {{C}_{1,1}}}{{{\mathit{EI}}_{1}}}+\frac{{{C}_{1,0}}}{{{\mathit{EI}}_{1}}}+\frac{{{q}_{0}}\cdot {{l}_{1}^{4}}}{24\cdot {{\mathit{EI}}_{1}}}&=\frac{{{C}_{2,0}}}{{{\mathit{EI}}_{2}}}\\ \displaystyle  \frac{{{l}_{1}^{2}}\cdot {{C}_{1,3}}}{2\cdot {{\mathit{EI}}_{1}}}+\frac{{{l}_{1}}\cdot {{C}_{1,2}}}{{{\mathit{EI}}_{1}}}+\frac{{{C}_{1,1}}}{{{\mathit{EI}}_{1}}}+\frac{{{q}_{0}}\cdot {{l}_{1}^{3}}}{6\cdot {{\mathit{EI}}_{1}}}&=\frac{{{C}_{2,1}}}{{{\mathit{EI}}_{2}}}\\ \displaystyle  -\frac{{{C}_{2,0}}\cdot {{k}_{B}}}{{{\mathit{EI}}_{2}}}-{{C}_{2,3}}+{{C}_{1,3}}+{{q}_{0}}\cdot {{l}_{1}}&=0\\ \displaystyle  -{{C}_{2,2}}+{{l}_{1}}\cdot {{C}_{1,3}}+{{C}_{1,2}}+\frac{{{q}_{0}}\cdot {{l}_{1}^{2}}}{2}&={{M}_{B}}\\ \displaystyle  -\frac{{{l}_{2}^{2}}\cdot {{C}_{2,3}}\cdot r}{2\cdot {{\mathit{EI}}_{2}}}-\frac{{{l}_{2}}\cdot {{C}_{2,2}}\cdot r}{{{\mathit{EI}}_{2}}}-\frac{{{C}_{2,1}}\cdot r}{{{\mathit{EI}}_{2}}}&=\displaystyle\frac{{{l}_{2}^{3}}\cdot {{C}_{2,3}}}{6\cdot {{\mathit{EI}}_{2}}}+\frac{{{l}_{2}^{2}}\cdot {{C}_{2,2}}}{2\cdot {{\mathit{EI}}_{2}}}+\frac{{{l}_{2}}\cdot {{C}_{2,1}}}{{{\mathit{EI}}_{2}}}+\frac{{{C}_{2,0}}}{{{\mathit{EI}}_{2}}}\\ \displaystyle  -{{C}_{2,3}}\cdot r-{{l}_{2}}\cdot {{C}_{2,3}}-{{C}_{2,2}}&=0 \end{array}
</math>


für die Integrationskonstanten.
für die Integrationskonstanten.
<!-------------------------------------------------------------------------------->
<!-------------------------------------------------------------------------------->


Zeile 169: Zeile 162:
==tmp==
==tmp==


<!-------------------------------------------------------------------------------->
Das Gleichungssystem wollen wir als
 
<math>\underline{\underline{A}}\cdot\underline{x}= \underline{b}</math>
 
schreiben, also
 
<math>\begin{pmatrix}\frac{1}{{{\mathit{EI}}_{1}}} & 0 & 0 & 0 & 0 & 0 & 0 & 0\\ 0 & 0 & -1 & 0 & 0 & 0 & 0 & 0\\ \frac{1}{{{\mathit{EI}}_{1}}} & \frac{{{\ell}_{1}}}{{{\mathit{EI}}_{1}}} & \frac{{{\ell}_{1}^{2}}}{2\cdot {{\mathit{EI}}_{1}}} & \frac{{{\ell}_{1}^{3}}}{6\cdot {{\mathit{EI}}_{1}}} & -\frac{1}{{{\mathit{EI}}_{2}}} & 0 & 0 & 0\\ 0 & \frac{1}{{{\mathit{EI}}_{1}}} & \frac{{{\ell}_{1}}}{{{\mathit{EI}}_{1}}} & \frac{{{\ell}_{1}^{2}}}{2\cdot {{\mathit{EI}}_{1}}} & 0 & -\frac{1}{{{\mathit{EI}}_{2}}} & 0 & 0\\ 0 & 0 & 0 & 1 & -\frac{{{k}_{B}}}{{{\mathit{EI}}_{2}}} & 0 & 0 & -1\\ 0 & 0 & 1 & {{\ell}_{1}} & 0 & 0 & -1 & 0\\ 0 & 0 & 0 & 0 & -\frac{1}{{{\mathit{EI}}_{2}}} & -\frac{{{\ell}_{2}}+r}{{{\mathit{EI}}_{2}}} & -\frac{{{\ell}_{2}^{2}}+2\cdot {{\ell}_{2}}\cdot r}{2\cdot {{\mathit{EI}}_{2}}} & -\frac{{{\ell}_{2}^{3}}+3\cdot {{\ell}_{2}^{2}}\cdot r}{6\cdot {{\mathit{EI}}_{2}}}\\ 0 & 0 & 0 & 0 & 0 & 0 & -1 & -r-{{\ell}_{2}}\end{pmatrix}\cdot\begin{pmatrix}{{C}_{1,0}}\\ {{C}_{1,1}}\\ {{C}_{1,2}}\\ {{C}_{1,3}}\\ {{C}_{2,0}}\\ {{C}_{2,1}}\\ {{C}_{2,2}}\\ {{C}_{2,3}}\end{pmatrix}=\begin{pmatrix}0\\ 0\\ -\frac{{{q}_{0}}\cdot {{\ell}_{1}^{4}}}{24\cdot {{\mathit{EI}}_{1}}}\\ -\frac{{{q}_{0}}\cdot {{\ell}_{1}^{3}}}{6\cdot {{\mathit{EI}}_{1}}}\\ -{{q}_{0}}\cdot {{\ell}_{1}}\\ {{M}_{B}}-\frac{{{q}_{0}}\cdot {{\ell}_{1}^{2}}}{2}\\ 0\\ 0\end{pmatrix}</math>.
 
Die Matrix-Elemente sind für die Koeffizientenmatrix
 
<math>\begin{array}{l} a_{1,1} = 1/EI_{1}\\ a_{2,3} = -1\\ a_{3,1} = 1/EI_{1}\\ a_{3,2} = \ell_{1}/EI_{1}\\ a_{3,3} = \ell_{1}^2/(2\cdot EI_{1})\\ a_{3,4} = \ell_{1}^3/(6\cdot EI_{1})\\ a_{3,5} = -1/EI_{2}\\ a_{4,2} = 1/EI_{1}\\ a_{4,3} = \ell_{1}/EI_{1}\\ a_{4,4} = \ell_{1}^2/(2\cdot EI_{1})\\ a_{4,6} = -1/EI_{2}\\ a_{5,4} = 1\\ a_{5,5} = -k_B/EI_{2}\\ a_{5,8} = -1\\ a_{6,3} = 1\\ a_{6,4} = \ell_{1}\\ a_{6,7} = -1\\ a_{7,5} = -1/EI_{2}\\ a_{7,6} = -(r+\ell_{2})/EI_{2}\\ a_{7,7} = -(2\cdot \ell_{2}\cdot r+\ell_{2}^2)/(2\cdot EI_{2})\\ a_{7,8} = -(3\cdot \ell_{2}^2\cdot r+\ell_{2}^3)/(6\cdot EI_{2})\\ a_{8,7} = -1\\ a_{8,8} = -r-\ell_{2}\\ \end{array}</math>
 
und für die rechte Seite
 
<math>\begin{array}{l} b_{1} = 0\\ b_{2} = 0\\ b_{3} = -(q_0\cdot\ell_{1}^4)/(24\cdot EI_{1})\\ b_{4} = -(q_0\cdot\ell_{1}^3)/(6\cdot EI_{1})\\ b_{5} = -q_0\cdot\ell_{1}\\ b_{6} = M_B-(q_0\cdot\ell_{1}^2)/2\\ b_{7} = 0\\ b_{8} = 0 \end{array}</math><!-------------------------------------------------------------------------------->
 
 
{{MyCodeBlock|title=Prepare for Solver
{{MyCodeBlock|title=Prepare for Solver
|text=Text
|text=Text
Zeile 180: Zeile 189:
==tmp==
==tmp==


<!-------------------------------------------------------------------------------->
Das Lösen des Gleichungssystems liefert
 
<math>\begin{array}{ll} \displaystyle {{C}_{1,0}}&\displaystyle =0\\ \displaystyle {{C}_{1,1}}&\displaystyle =-\frac{{{q}_{0}}\cdot {{\ell}_{1}^{3}}}{56}\\ \displaystyle {{C}_{1,2}}&\displaystyle =0\\ \displaystyle {{C}_{1,3}}&\displaystyle =-\frac{{{q}_{0}}\cdot {{\ell}_{1}}}{7}\\ \displaystyle {{C}_{2,0}}&\displaystyle =0\\ \displaystyle {{C}_{2,1}}&\displaystyle =\frac{13\cdot {{q}_{0}}\cdot {{\ell}_{1}^{3}}}{84}\\ \displaystyle {{C}_{2,2}}&\displaystyle =-\frac{9\cdot {{q}_{0}}\cdot {{\ell}_{1}^{2}}}{14}\\ \displaystyle {{C}_{2,3}}&\displaystyle =\frac{6\cdot {{q}_{0}}\cdot {{\ell}_{1}}}{7} \end{array}</math><!-------------------------------------------------------------------------------->
 
{{MyCodeBlock|title=Solving
{{MyCodeBlock|title=Solving
|text=Text
|text=Text
Zeile 191: Zeile 203:
==tmp==
==tmp==


<!-------------------------------------------------------------------------------->
Zum Auftragen der Ergebnisse nutzen wir die  Standard-Lösungen - Lastfall 3 mit der maximalen Auslenkung
 
<math>W^* = \displaystyle \frac{5\;\ell_1^4\cdot q_0}{384 \; EI_1}</math>
 
und dem Winkel
 
<math>\displaystyle \Phi^* = \frac{q_0\;\ell_1^3}{24\;EI_1}</math>
 
Die Ergebnisse sind dann ...<!-------------------------------------------------------------------------------->
 
 
 
 
{{MyCodeBlock|title=Post-Processing
{{MyCodeBlock|title=Post-Processing
|text=Text
|text=Text

Version vom 7. April 2021, 14:42 Uhr


Aufgabenstellung

Ein Stab ABC (E-Modul: E) besteht aus zwei Sektionen mit den Längen 1 bzw. 2 sowie den Flächenmomenten I1 bzw. I2. Die Sektionen haben jeweils einen quadratischen Querschnitt, Sektion AB ist durch eine konstante Streckenlast q0 belastet, in B wirkt das Moment MB0. Der Stab ist in A durch ein gelenkiges Festlager gelagert. In C ist das Stabende fest mit dem Umfang einer Rolle vom Radius r verbunden, die in D frei drehbar gelagert ist. In B sind die beiden Sektionen fest miteinander verbunden. Die Feder in B ist eine Translationsfeder mit der Steifigkeit  kB.

Interessant ist die kinematische Randbedingung aus der Rolle.


Lageplan

Gesucht ist die Analytische Lösung für den Euler-Bernoulli-Balken und die Verläufe der Schnittgrößen.

Parameter

Ermitteln Sie für ein Euler-Bernoulli-Modell die analytischen Verläufe der Schnittgrößen und Verschiebungen im Balken für diese Parameter:

Erstellen Sie dazu ein Programm, mit einem Euler-Bernoulli-Modell für die Berechnung der analytischen Verläufe der Schnittgrößen und Verschiebungen im Balken. Bestimmen Sie Querschnitts-Abmessungen der Sektionen so, dass die maximale Auslenkung des Systems 10 mm beträgt.

Lösung mit Maxima

Die Aufgabe ist ein klassisches Randwertproblem:

  1. zwei Gebiete, in denen ein Euler-Bernoulli-Balken in AB und BC durch eine Streckenlast q belastet ist (in Bereich II ist die Streckenlast allerdings Null) und somit durch die Differentialbeziehung

berschrieben wird.

  1. Rand- und Übergangsbedingungen in den Punkten A, B, C

Wir verwenden xi und ξi als Koordinaten je Bereich, in der Übersicht sieht das Randwertproblem so aus:

Rand
A
Bereich IÜbergang
B
Bereich IIRand
C


tmp

Wir lösen die Feld-Differentialgleichung exakt und passen die Integrationskonstanten an die Rand- und Übergangsbedingungen an.

Header

Text


1+1




tmp

Wir definieren zunächst die bekannten Parameter zu

.

Weitere brauchen nicht.

Declarations

Text


1+1




tmp

In Bereich I und II gilt dieselbe Bewegungs-Differentialgleichung

,

die wir durch Integration lösen und dann bereichsweise anpassen.

So gilt für Bereich II: q0 = 0.

Die allgemeine Lösung ist mit

... für Bereich I:

... für Bereich II:

.


Generic Solutions for Euler-Bernoulli-Beam

Text


1+1




tmp

Für die 2*4 = 8 Integrationskonstanten

suchen wir jetzt die passenden Gleichungen aus Rand- und Übergangsbedingungen.

Zur besseren Übersicht nennen wir die Schnitt-Momente und -Kräfte nach den jeweiligen Knotenpunkten A, B, C und fügen als Index ein + / - hinzu, um die Seite (+: rechts vom Knoten, -: links vom Knoten) zu kennzeichnen, wenn diese unterschiedlich sind.

Aus Rand "A"

Geometrische Randbedingungen

Kraft- und Momenten-Randbedingungen

Aus Übergang "B"

Geometrische Randbedingungen

Kraft- und Momenten-Randbedingungen

Aus Rand "C"


Geometrische Randbedingungen

Kraft- und Momenten-Randbedingungen

Einsetzen liefert 8 Gleichungen

für die Integrationskonstanten.


Boundary Conditions

Text


1+1




tmp

Das Gleichungssystem wollen wir als

schreiben, also

.

Die Matrix-Elemente sind für die Koeffizientenmatrix

und für die rechte Seite


Prepare for Solver

Text


1+1




tmp

Das Lösen des Gleichungssystems liefert

Solving

Text


1+1




tmp

Zum Auftragen der Ergebnisse nutzen wir die  Standard-Lösungen - Lastfall 3 mit der maximalen Auslenkung

und dem Winkel

Die Ergebnisse sind dann ...



Post-Processing

Text


1+1






Links

  • ...

Literature

  • ...



Biegelinie w(x)
Biegemoment M(x)
Querkraft Q(x)
Lageplan
Kippung w'(x)