Gelöste Aufgaben/FEB4

Aufgabenstellung

Wie in FEB3 geht es hier um die Ausnutzung von Symmetrien mit einer skalaren Gleichgewichtsformulierung.

Lageplan

Das skizzierte System besteht aus Euler-Bernoulli-Balken (jeweils Länge a) der Biegesteifigkeit EI. Die Konstruktion wird in A durch ein gelenkiges Festlager gehalten und in B durch die Kraft F belastet.

Gesucht ist die Ersatzfeder-Steifigkeit des Systems, das Sie aus Finiten Elementen mit Euler-Bernoulli-Balken zusammensetzen.

Gesucht ist die Auslenkung des Punktes B. Das Modell soll aus wenigen Finiten Elementen bestehen, überlegen Sie, welche Modelle sinnvoll sind. Die neutralen Fasern der Balken seien in Längsrichtung undehnbar – die Querschnitts-Schwerpunkte verschieben sich also nicht in Balken-Längsrichtung.

Gegeben: a, E I, F

tmp

Für das FE-Modell wollen wir nicht alle Sektionen der Feder modellieren - wir schneiden die Feder entlang ihrer Längsachse durch und nehmen eins der so entstandenen Elemente exemplarisch heraus.

Federelement.

Von diesen Elementen haben wir fünf in der Feder - die Längung dieses Elements der Feder entspricht also 1/5 der Feder-Längung in B.

Im Schnitt der horizontalen Sektionen mit der Längsachse der Feder ist das Schnittmoment jeweils Null - wir dürfen uns also an den Rändern unseres Sub-Systems ein Gelenk "denken" und mit dem hier gezeichneten System arbeiten.

Aber Achtung: so wie gezeichnet ist es statisch unbestimmt: wir können im Gleichgewicht beliebige Verschiebungen des Systems in Feder-Längsrichtung angeben.

Geht's noch einfacher?

• Wie müsste ein System aussehen, das noch einfacher gestrickt ist?
• Welche neuen Randbedingungen ergeben sich?

Header

Text

1+1

tmp

Wir haben also ein Ersatz-System mit drei Finiten Elementen und vier Knoten:

Knoten und Stäbe.

Die Knoten 0 und 3 entstehen in unserem Ersatz-System, weil hier das Biegemoment in den horizontalen Sektionen (I und III) Null ist. Diese beiden Sektionen sind nur a/2 lang, die vertikale Sektion ist a lang.

Damit das System in Längsrichtung festgelegt ist, führen wir in Knoten "0" ein gelenkiges Festlager ein - die Verschiebung in Knoten "3" ist dann die Auslenkung von Knoten "3" relativ zu Knoten "0".

Dementsprechnd haben die drei Finiten Elementen die Längen

${\begin{array}{l}\ell _{1}=a/2\\\ell _{2}=a\\\ell _{3}=a/2\\\end{array}}$ .

Die Element-Steifigkeitsmatrix kopieren wir aus Abschnitt FEM-Formulierung für den Euler-Bernoulli-Balken zu:

$\displaystyle K_{i}(\ell _{i})={\frac {EI}{\ell _{i}^{3}}}\cdot {\begin{pmatrix}12&6\cdot \ell &-12&6\cdot \ell _{i}\\6\cdot \ell _{i}&4\cdot {{\ell _{i}}^{2}}&-6\cdot \ell _{i}&2\cdot {{\ell _{i}}^{2}}\\-12&-6\cdot \ell _{i}&12&-6\cdot \ell _{i}\\6\cdot \ell _{i}&2\cdot {{\ell _{i}}^{2}}&-6\cdot \ell _{i}&4\cdot {{\ell _{i}}^{2}}\end{pmatrix}}$ .===Declarations=== Text

1+1

tmp

Für jedes Finite Element führen wir zuerst das standardisierte lokale x/z-Koordinatensystem ein. Und wir definieren neue, anschauliche - globale - Koordinaten des Systems, um sie von den Standard-Koordinaten W_i , Φ_i unterscheiden zu können:

Lokale Orts-Koordinaten xi, zi	Globale Knoten-Koordinaten Ui, Ψi

Wir haben also je Element die eingezeichneten globalen Knoten-Koordinaten

$U_{i-1},\Psi _{i-1},U_{i},\Psi _{i}$

sowie die lokalen FE-Koordinaten  

$W_{i-1},\Phi _{i-1},W_{i},\Phi _{i}$.

die wir miteinander "verheiraten" müssen. Das machen wir beim Komponieren des Gesamt-Gleichungssystems.

Die virtuelle Formänderungsenergie eines Finiten Elements ist nun

$\delta \Pi _{i}=\left(W_{i-1},\Phi _{i-1},W_{i},\Phi _{i}\right)^{T}\cdot {\underline {\underline {K}}}_{i}(\ell _{i})\cdot \left({\begin{array}{c}W_{i-1}\\\Phi _{i-1}\\W_{i}\\\Phi _{i}\end{array}}\right)$,

mit

$\delta \Pi =\displaystyle \sum _{1}^{3}\delta \Pi _{i}$.

Die virtuelle Arbeit der Kraft ist dabei

$\delta W^{a}=F\cdot \delta U_{3}$

Aus den Bildern für die Knoten-Koordinaten oben lesen wir nun die (geometrischen) Rand- und Übergangsbedingungen ab und skizzieren darunter die Einarbeitung dieser Bedingungen in die virtuellen Formänderungsenergien.

Element I	Element II	Element III

Dabei sind folgende Koordinaten sind in den drei Elementen gleich oder unterscheiden sich nur durch ein Vorzeichen:

${\begin{array}{ccccc}I&&II&&III\\\Phi _{1}&=&\Phi _{1}&&\\&&\Phi _{2}&=&\Phi _{2}\\W_{1}&&=&&-W_{2}&{\text{(Sektion II wird nicht länger!)}}\end{array}}$

Was noch fehlt ist die Übersetzung der lokalen in die globalen Koordinaten

lokal		global
	↔

Text

1+1

tmp

Globale Koordinaten

===Assembly of System Matrices===

Text

1+1

tmp

Solving

Text

1+1

tmp

Verformung des Elements.

===Post-Processing===

Text

1+1

Links

• ...

Literature

• ...

g