RK bei DGL 2.Ordnung < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) überfällig  | | Datum: | 23:36 Do 03.04.2008 | | Autor: | yuffie |
| Aufgabe | Aufgabe
DGL eines Balkens der Länge l lautet (a(x)y´´)´´ =f(x);
dabei beschreibt die Fkt a(x) die Materialeigenschaft- insbesondere die Querschnittsgestaltung- über der Längskoordinate x mit 0<x<l.
Es gelte a(x)=A-B(2x-l)² mit A>Bl² und f(x) identisch 1. Bestimmen sie dieses Problem für die Randbedingeungen:
y(0)=y(l)=0
y´(0)=y´´(l)=0
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt. |
Hallo, ich habe meine Aufgabe schon so verkürzt, dass ich aus der Differentialgleichung 4.Odnuung eine DGL 2.Ordnung gemacht habe.
y´´(x)=0.5(x²-l²)+(x-l)c1/ a(x)
mit a(x)=A- B(2x-l)²
Meine Materialeigenschaften a(x)= A- B(2x-l)². Und dieses A und B muss ja minimiert werden, je nachdem wie lang mein Balken ist, verhält sich A und B. Und dieses Problem mit Hilfe von Matlab unter Runge-Kutta hinzukriegen klappt nicht.
y´1= y2 y(0)=0
y´2= y´´(x) y´(0)=0
ist meine erste Transformation, um Runge-Kutta anzuwenden.
Hierbei habe ich nun das Problem, dass er A und B nicht als definierte Variablen anerkennt.
Desweiteren stört mich, dass meine zweite Umformung, nur noch von x abhäng und nicht mehr vom y. Mein Professor meint, das ist so korekkt, will mir aber nicht weiterhelfen, wie ich es am besten in Matlab eingebe.
m-Files
function function y=haus(x,c1)
y=[0; 0];
y(1)=y(2);
$ [mm] y(2)=0.5\cdot{}((x^2)-1)+(x-1)\cdot{}c1/A-B\cdot{}(2\cdot{}x-1)^2; [/mm] $
l-habe ich 1 gesetz(Prof meint kann ich als const ansehen)
nach weiten Überlegungen, bin ich der Meinung, dass das Ergebnis bei RK =0 ist.
Wer ist noch der Meinung? Denn y(2) ist nicht von y abhängig und kann als Wert nicht eingesetzt werden bei RK. y=f(yj+h*summe aji*kj)
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 00:20 Sa 12.04.2008 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
