FEM Steifigkeitsmatrix < DGL < Numerik < Hochschule < Mathe < Vorhilfe
|
Zur Berechnung der Steifigkeitmatrix bei der FEM-Methode habe ich im Skript folgende Formel gegeben (Siehe Anhang).
[Dateianhang nicht öffentlich]
Mit Hilfe der Formel wollte ich nun die Aufgabe [Dateianhang nicht öffentlich] lösen.Kriegs aber einfach nicht hin
Ausserdem bekannt ist:
[mm] v_1'=\frac{1}{h} [/mm] in [0,0.2]
[mm] v_1'=\frac{-1}{h} [/mm] in [0.2,0.4]
[mm] v_2'=\frac{1}{h} [/mm] in [0.2,0.4]
[mm] v_2'=\frac{-1}{h} [/mm] in [0.4,0.6]
[mm] v_3'=\frac{1}{h} [/mm] in [0.4,0.6]
usw.
also
[mm] a_{12}=\int_{0.2}^{0.4} \! [/mm] (4x+3) [mm] \, [/mm] dx * [mm] v_1'v_2' [/mm] = [mm] \int_{0.2}^{0.4} \! -\frac{4x+3}{h^2} \, [/mm] dx=-4 [mm] \int_{0.2}^{0.4} \! \frac{x}{h^2} \, [/mm] dx -3 [mm] \int_{0.2}^{0.4} \! \frac{1}{h^2} \, [/mm] dx
das wiederum ist:
[mm] -\frac{4}{h} \left[\frac{1}{2}x^2\right]_{0.2}^{0.4}-\frac{3}{h}\left[x\right]_{0.2}^{0.4} =-\frac{0.84}{h}
[/mm]
Leider ist das aber totaler Blödsinn,rauskommen müsste nämlich -21/h.
Was stimmt hier nicht? Ich vermute ich habe irgendwas beim integrieren verkehrt gemacht....
HELP PLS
Dateianhänge: Anhang Nr. 1 (Typ: PNG) [nicht öffentlich] Anhang Nr. 2 (Typ: PNG) [nicht öffentlich]
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 19:20 Do 06.09.2012 | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|