Gelöste Aufgaben/StaF: Unterschied zwischen den Versionen

Aus numpedia
Zur Navigation springen Zur Suche springen
Keine Bearbeitungszusammenfassung
Keine Bearbeitungszusammenfassung
Zeile 1: Zeile 1:
[[index.php?title=Kategorie:Gelöste Aufgaben]]
[[Category:Gelöste Aufgaben]]
[[index.php?title=Kategorie:A*x=b]]
[[Category:A*x=b]]
[[index.php?title=Kategorie:Lineare Algebra]]
[[Category:Lineare Algebra]]
[[index.php?title=Kategorie:Numerische Lösung]]
[[Category:Numerische Lösung]]
[[index.php?title=Kategorie:Randwertproblem]]
[[Category:Randwertproblem]]
[[index.php?title=Kategorie:Prinzip der virtuellen Verrückungen]]
[[Category:Prinzip der virtuellen Verrückungen]]
[[index.php?title=Kategorie:Stab]][[index.php?title=Kategorie:Dehnstab]]
[[Category:Stab]]
[[index.php?title=Kategorie:Euler-Bernoulli-Balken]]
[[Category:Dehnstab]]
[[index.php?title=Kategorie:Finite-Elemente-Methode]]
[[Category:Euler-Bernoulli-Balken]]
[[index.php?title=Kategorie:Maxima]]
[[Category:Finite-Elemente-Methode]]
[[index.php?title=Kategorie:Stabwerk]]
[[Category:Maxima]]
[[Category:Stabwerk]]


==Aufgabenstellung==
==Aufgabenstellung==
Zeile 28: Zeile 29:


===Gleichgewichtsbedingungen===
===Gleichgewichtsbedingungen===
Für die Gleichgewichtsbedingung nach dem [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Werkzeuge/Gleichgewichtsbedingungen/Arbeitsprinzipe_der_Analytischen_Mechanik/Prinzip_der_virtuellen_Verr%C3%BCckungen Prinzip der virtuellen Verrückungen] ist
Für die Gleichgewichtsbedingung nach dem [https://numpedia.rzbt.haw-hamburg.de/index.php?title=Werkzeuge/Gleichgewichtsbedingungen/Arbeitsprinzipe_der_Analytischen_Mechanik/Prinzip_der_virtuellen_Verr%C3%BCckungen Prinzip der virtuellen Verrückungen]
::<math>
\begin{array}{ccc}
\delta W & \stackrel{!}{=} & 0\\
&=& \delta\Pi - \delta W^a
\end{array}
</math>
benötigen wir die virtuelle Formänderungsenergie <math>\delta \Pi</math> und die virtuelle Arbeit der äußeren Kraft <math>\delta W^a</math> der äußeren Kräfte und Momente.
 
Mit den Konventionen für die Knoten-Verschiebungen aus [[https://numpedia.rzbt.haw-hamburg.de/index.php?title=Gel%C3%B6ste_Aufgaben/Stab Stab]] ist  
::<math>
\delta W^a = -\delta W_{4,0} \cdot F
</math>.
Für <math>\delta \Pi</math> gilt
::<math>
\delta \Pi = \sum_{i=0}^4 \delta \Pi_i
</math>
mit den virtuellen Formänderungsarbeiten der vier Stäbe.
 
Dabei haben wir Anteile der Arbeit aus der [[https://numpedia.rzbt.haw-hamburg.de/index.php?title=Sources/Anleitungen/FEM-Formulierung_f%C3%BCr_den_Euler-Bernoulli-Balken Biegung]] und der Längs-Dehnung des Stabes.
 
Für den Stab ''k'' mit den Knoten ''I'' und ''J'' haben wir als Koodinaten der Knoten
::<math>U_{I,k},W_{I,k},\Phi_{I,k}</math> und <math>U_{J,k},W_{J,k},\Phi_{J,k}</math>.
Damit haben wir
::<math>
\delta \Pi_k =
\left(\delta W_{I,k},\delta \Phi_{I,k},\delta \Phi_{J,k},\delta W_{J,k}\right)
\cdot
\frac{EI}{\ell_i^3}
\cdot
\begin{pmatrix}12 & 6\, {\ell_i} & -12 & 6\, {\ell_i}\\ 6\, {\ell_i} & 4\, {\ell_i^{2}} & -6\, {\ell_i} & 2\, {\ell_i^{2}}\\ -12 & -6\, {\ell_i} & 12 & -6\, {\ell_i}\\ 6\, {\ell_i} & 2\, {\ell_i^{2}} & -6\, {\ell_i} & 4\, {\ell_i^{2}}\end{pmatrix} \cdot
\left(\begin{array}{c}W_{I,k}\\ \Phi_{I,k}\\W_{J,k}\\ \Phi_{J,k}\end{array}\right)
+
\left(\delta U_{I,k},\delta U_{J,k}\right)
\cdot
\frac{EA}{\ell_i}
\cdot
\begin{pmatrix}1&-1\\-1&1\end{pmatrix}
\cdot
\left(\begin{array}{c}U_{I,k}\\U_{J,k}\end{array}\right)
</math>
 
 
 
Für den Stab ''k'' definieren wir
::<math>\underline{Q} = \left(\begin{array}{c}
U_{I,k}\\
W_{I,k}\\
\Phi_{I,k}\\
U_{J,k}\\
W_{J,k}\\
\Phi_{J,k}\\
\end{array}\right)</math> sowie <math>\delta\underline{Q} = \left(\begin{array}{c}
\delta U_{I,k}\\
\delta W_{I,k}\\
\delta \Phi_{I,k}\\
\delta U_{J,k}\\
\delta W_{J,k}\\
\delta \Phi_{J,k}\\
\end{array}\right)</math>
 
 
 
 
 
 
 
 
 
 




W = 0
        = delta \Pi - \delta Wa
benötigen wir die virtuelle Formänderungsenergie \delta \Pi und die virtuelle Arbeit der äußeren Kraft F.
Mit allen Konventionen für die Knoten-Verschiebungen ist
Delta Wa =  -delta W_{4,0} *F
Für \delta \Pi gilt





Version vom 21. Oktober 2024, 20:15 Uhr


Aufgabenstellung

Wir untersuchen die Belastung eines ebenen Stabwerks. Die Stäbe haben wie skizziert die Länge ℓ bzw. ℓ/2. Die Struktur wird mit der Kraft F belastet.


Caption

Gesucht ist ein Vergleich zwischen der klassischen Stabwerkstheorie und einer Herangehensweise, bei der wir eine feste Verbindung der Stäbe in den Knoten ansetzten. Grundlage des Modells ist die FEM-Lösung der Felddifferentialgleichung im Vergleich zur Lösung in Problemstellung „Stab“.

Wir stellen das Modell des Stabwerks mit dem Prinzip der virtuellen Verrückungen auf und vergleichen, wie sich diese von der Herangehensweise aus „Stab“ mit der analytischen Lösung unterscheidet.

Lösung mit Maxima

Wir nutzen das Computer-Algebra-System Maxima zur Lösung. Das macht hier Sinn, weil wir die Herangehensweise mit der aus Stab vergleichen wollen – für die wir ebenfalls Maxima eingesetzt haben.

Declarations

Wir übernehmen alle Vereinbarungen und Parameter aus der Problemformulierung „Stab“.

Gleichgewichtsbedingungen

Für die Gleichgewichtsbedingung nach dem Prinzip der virtuellen Verrückungen

δW=!0=δΠδWa

benötigen wir die virtuelle Formänderungsenergie δΠ und die virtuelle Arbeit der äußeren Kraft δWa der äußeren Kräfte und Momente.

Mit den Konventionen für die Knoten-Verschiebungen aus [Stab] ist

δWa=δW4,0F.

Für δΠ gilt

δΠ=i=04δΠi

mit den virtuellen Formänderungsarbeiten der vier Stäbe.

Dabei haben wir Anteile der Arbeit aus der [Biegung] und der Längs-Dehnung des Stabes.

Für den Stab k mit den Knoten I und J haben wir als Koodinaten der Knoten

UI,k,WI,k,ΦI,k und UJ,k,WJ,k,ΦJ,k.

Damit haben wir

δΠk=(δWI,k,δΦI,k,δΦJ,k,δWJ,k)EIi3(126i126i6i4i26i2i2126i126i6i2i26i4i2)(WI,kΦI,kWJ,kΦJ,k)+(δUI,k,δUJ,k)EAi(1111)(UI,kUJ,k)


Für den Stab k definieren wir

Q_=(UI,kWI,kΦI,kUJ,kWJ,kΦJ,k) sowie δQ_=(δUI,kδWI,kδΦI,kδUJ,kδWJ,kδΦJ,k)








(2A2Eη300(2A2Eη02)0A2Eη3000002A2Eη30(A2Eη402)(3A2Eη402)A2Eη60(A2Eη202)0A2Eη60(2A2Eη02)(A2Eη402)A2Eη+302AE203+8A2Eη030(2A2Eη02)(A2Eη+302AE403)3A2Eη302AE403A2Eη4020(3A2Eη402)03A2Eη+02AE203+2AE0(3A2Eη202)3A2Eη302AE403(3A2Eη+02AE403)(3A2Eη402)A2Eη30A2Eη60(2A2Eη02)(3A2Eη202)4A2Eη30(A2Eη402)3A2Eη402A2Eη600(A2Eη202)(A2Eη+302AE403)3A2Eη302AE403(A2Eη402)A2Eη+302AE403+A2Eη03(3A2Eη302AE403)(3A2Eη402)003A2Eη302AE403(3A2Eη+02AE403)3A2Eη402(3A2Eη302AE403)3A2Eη+02AE403+AE03A2Eη4020A2Eη60A2Eη402(3A2Eη402)A2Eη60(3A2Eη402)3A2Eη4022A2Eη30)(W1,0U1,0Φ1,0W2,0U2,0Φ2,0W3,0U3,0Φ3,0W4,0U4,0Φ4,0)=(00000F00)

Hier kommt jetzt irgendein Text.

SomeText

Title

Text


1+1




Element-Steigigkeitsmatrizen mit globalen Koordinaten
Element #1

k1=(8A2Eη0302A2Eη02(8A2Eη03)02A2Eη0202AE000(2AE0)02A2Eη0202A2Eη30(2A2Eη02)0A2Eη30(8A2Eη03)0(2A2Eη02)8A2Eη030(2A2Eη02)0(2AE0)002AE002A2Eη020A2Eη30(2A2Eη02)02A2Eη30)

Element #2

k2=(A2Eη+302AE4033A2Eη302AE403A2Eη402(A2Eη+302AE403)(3A2Eη302AE403)A2Eη4023A2Eη302AE4033A2Eη+02AE4033A2Eη402(3A2Eη302AE403)(3A2Eη+02AE403)3A2Eη402A2Eη4023A2Eη402A2Eη30(A2Eη402)(3A2Eη402)A2Eη60(A2Eη+302AE403)(3A2Eη302AE403)(A2Eη402)A2Eη+302AE4033A2Eη302AE403(A2Eη402)(3A2Eη302AE403)(3A2Eη+02AE403)(3A2Eη402)3A2Eη302AE4033A2Eη+02AE403(3A2Eη402)A2Eη4023A2Eη402A2Eη60(A2Eη402)(3A2Eη402)A2Eη30)

Element #3

k3=(A2Eη+302AE403(3A2Eη302AE403)A2Eη402(A2Eη+302AE403)3A2Eη302AE403A2Eη402(3A2Eη302AE403)3A2Eη+02AE403(3A2Eη402)3A2Eη302AE403(3A2Eη+02AE403)(3A2Eη402)A2Eη402(3A2Eη402)A2Eη30(A2Eη402)3A2Eη402A2Eη60(A2Eη+302AE403)3A2Eη302AE403(A2Eη402)A2Eη+302AE403(3A2Eη302AE403)(A2Eη402)3A2Eη302AE403(3A2Eη+02AE403)3A2Eη402(3A2Eη302AE403)3A2Eη+02AE4033A2Eη402A2Eη402(3A2Eη402)A2Eη60(A2Eη402)3A2Eη402A2Eη30)

Element #4

k4=(A2Eη030A2Eη202(A2Eη03)0A2Eη2020AE000(AE0)0A2Eη2020A2Eη30(A2Eη202)0A2Eη60(A2Eη03)0(A2Eη202)A2Eη030(A2Eη202)0(AE0)00AE00A2Eη2020A2Eη60(A2Eη202)0A2Eη30)


Links

  • ...

Literature

  • ...