FEM Lastvektor ermittlung < Numerik < Hochschule < Mathe < Vorhilfe
|
| Status: |
(Frage) überfällig  | | Datum: | 10:41 Mi 09.04.2008 | | Autor: | skydave |
Hallo,
ich habe eine Frage bezüglich den Werten des Lastvektors bei der Lösung eines klassischen Elastizitätsproblems mittels der finite Elemente Methode.
kurz zum Kontext:
Zur Ermittlung der Verformung eines Gegenstandes bei äußerem Krafteinfluß wird die Erhaltungsgleichung "innere Arbeit = äußere Arbeit" mittels FEM gelöst. Dabei wird die innere Arbeit durch ein Integral der potentiellen Elastitizätsenergie über das Volumen des Körpers in Abhängigkeit von der (unbekannten) Verformung des Körpers beschrieben.
Auf der rechten Seite werden die extern einwirkenden Kräfte durch ein Oberflächenintegral aufgesammelt.
Diese Formulierung mündet dann bei der FEM in ein Gleichungssystem Ax=b
wobei hier A die Steifigkeitsmatrix ist und Ax die gesamte innere Arbeit an den Knoten beschreibt. Der Störvektor b ist der Lastvektor und entspricht den externen Kräften die an den jeweiligen Knoten anliegen.
Nun zur Frage:
Für Knoten an denen Dirichlet-Randbedingungen definiert sind lassen sich die Werte des Störvektors b ermitteln. Wie aber ermittel ich die Werte des Störvektors an den freien Knoten? Ich hab hier ein Verständnissproblem und würde mich freuen wenn mir jemand bei der Klärung helfen könnte.
Noch ein Beispiel um die Frage zu verdeutlichen:
Gegeben ist zum Beispiel eine diskretisierte Wand (Knoten und Elemente).
Diese Wand ist fest im Boden verankert, hat also Dirichlet-Randbedingungen an den Knoten der Wandunterseite. Zusätzlich ist eine Reihe von Knoten an einer oberen Kante verschoben (ebenfalls Dirichlet Randbedingungen). Mittels FEM kann nun ohne weiteres die Verschiebung aller freien Knoten berechnet werden. Für mich stellt sich aber die Frage, wie bei der Aufstellung des Gleichungssystem die Werte des Störvektors b für die freien Knoten ermittelt werden.
Vielen Dank schonmal im Voraus!
David
P.S.:
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
P.P.S.:
Falls die Frage in dieses Forum nicht hingehört würde ich mich über Hinweise freuen wo ich sie besser stellen sollte.
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 11:20 Sa 10.05.2008 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
