Verfasst am: 24.11.2014, 14:19
Titel: u=k\f ist der "problemlöser" für FEM wie läuft
Hallo neztwelt!
Meine Frage bezieht sich auf u=k\f Wie genau läuft die Berechnung ab.
Wird dort mit Navier-Stokes gearbeitet oder Euler oder oder?
Und dann noch etwas:
In der MAtLAb Hilfe wird von einer stiffness matrix gesprochen, wozu brauche ich dieses beim Berechnen von nichtlinearen partiellen DGL bei newtonschen Fluiden?
Bei u = k\f geht es ausschließlich um die Lösung eines linearen Gleichungssystems. Das hat also mit Navier-Stokes und dergleichen überhaupt nichts mehr zu tun. Lediglich das Aufstellen der Gleichungen (also: was sind u, k, f) kann damit zu tun haben.
zur zweiten Frage:
Bitte verlinken, wo du das genau gelesen hast.
Grüße,
Harald
navmet
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 05.12.2014, 10:45
Titel: u=k\f ist der "problemlöser" für FEM wie läuft
Hallo Harald,
vielen Dank für deine Antwort auf meine erste Frage. Es werden also lineare Gleichungssysteme gelöst. Es wird also Impuls, Energie und Kontinuität völlig aussen vor gelassen. Das heißt das auch die innere Reibung nicht beachtet wird, sehe ich das richtig?
zur zweiten Frage:
Ich habe in der Hilfe im Programm MatLab einfach mal u=k\f eingegeben und dannn findet sich unter assempde ein Absatz mit der Aussage das u=?k\f der Problemlöser sei und eben auch die stiffness matrix.
Es wird also Impuls, Energie und Kontinuität völlig aussen vor gelassen. Das heißt das auch die innere Reibung nicht beachtet wird, sehe ich das richtig?
Das wird davon abhängen, wie das Gleichungssystem konstruiert ist.
Zitat:
Ich habe in der Hilfe im Programm MatLab einfach mal u=k\f eingegeben
Dann bekomme ich keine Treffer.
Grüße,
Harald
navmet
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 08.12.2014, 15:13
Titel: u=k\f ist der "problemlöser" für FEM wie läuft
Hallo Harald,
habe MatLab Version 2011a kann vielleicht daran liegen.
Jedenfalls kommt bei mir wenn ich u=k\f eingebe:
[K,F] = assempde(b,p,e,t,c,a,f) assembles the PDE problem by approximating the dirichlet boundary condition with stiff springs ( see the eliptic System for details). K and F are the stiffness Matrix and right-hand side, respctively. The solution to the FEM formulation of the PDE problem is u=k\f.
Deshalb kam ich darauf, dass bei u=k\f eine partielle differentialgleichung (PDE) gelöst wird.
es wird definitiv ein lineares Gleichungssystem gelöst.
Allerdings ist es natürlich möglich, dass die Gleichungen eine partielle DGL beschreiben - und genau das wird die Idee von assempde sein.
Grüße,
Harald
navmet
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 11.12.2014, 16:02
Titel: u=k\f ist der "problemlöser" für FEM wie läuft
Hallo Harald,
falls du dich mit dem folgenden Code auskennst, kannst du mir dazu eine Hilfe geben?
Code:
for e=1:T
Knoten=t(e,:);
Pe=[ones(3,1),p(Knoten,:)];
Flaeche=abs(det(Pe))/2;
C=inv(Pe);
grad=C(2:3,:);
Ke=Flaeche*(grad')*grad;
Fe=Flaeche/3;
K(Knoten,Knoten)=K(Knoten,Knoten)+Ke;
F(Knoten)=F(Knoten)+Fe;
end
Schön, wenn Sie mir helfen könnten.
Wenn diesmal nicht, möchte ich mich trotzdem bedanken, denn Sie haben mir bisher gut geholfen.
Danke für Ihre Mühe und Zeit !!
ich kann da leider auch nur sagen, was die Operationen im einzelnen machen. Warum man dies bei der Anwendung so macht? Keine Ahnung.
Leider scheint der Autor des Codes auf das Kommentieren verzichtet zu haben.
Grüße,
Harald
navmet
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 15.12.2014, 09:55
Titel: u=k\f ist der "problemlöser" für FEM wie läuft
Hallo Harald,
vielen Dank trotzdem.
Ich schließe dann dieses Thema.
Freundliche Grüße und eine angenehme Weihnachtszeit,
Dirk!
Einstellungen und Berechtigungen
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
RSS