|
|
| bensko |
Forum-Anfänger
|
 |
Beiträge: 43
|
|
|
|
Anmeldedatum: 01.10.13
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 12.04.2015, 16:42
Titel: Gleichungssystem lösen
|
 |
Hallo Zusammen, ich denke einer von Euch kann mir hier sicher schnell helfen.
Ich würde gerne folgendes mit Matlab Lösen können.
F = K * u
mit
F =
F1
3
2
K =
2 -2 0
-2 3 -1
0 -1 1
u =
u1
0
u3
Also ich möchte kleine FEM Berechnungen mittels Matlab durchführen.
Zur Aufgabe: 2 Stäbe die mit einander Verbunden sind --> 3 Knoten
Dann Habe ich die Steifigkeitsmatrix von beiden Aufgestellt und dann zu eine Globale Steifigkeitsmatrix zusammen geführt. Nun würde ich gerne die Verschiebungen und Kräfte an den Knoten herausbekommen.
Es müsste ja eigentlich u1 = -2 ; u3 = 2 und F1 = -4 herauskommen
Vielen Dank!!
|
|
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 18:32
Titel:
|
|
Hallo,
sollte mit fsolve (numerisch) oder solve (symbolisch) kein Problem sein.
Grüße,
Harald
|
|
|
|
| bensko |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 43
|
|
|
|
Anmeldedatum: 01.10.13
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 18:45
Titel:
|
|
Hallo Harald,
Ich habe unter der Hilfe von Matlab 'fsolve' gesucht. Bin da auf ein Beispiel gestoßen.
Daraus habe ich entnommen, dass ich als erstes meine Gleichungen umstellen soll. Sprich 3 Gleichungen :
eq1 = k(1:1 , 1:3)*u -F1;
eq2 = k(2:2 , 1:3)*u -F2;
eq3 = k(3:3 , 1:3)*u -F3;
Die weitere Vorgehensweise verstehe ich nicht. Wie muss ich nun vorgehen, sodass ich das Gl-System nach F1, u1 und u3 auflösen kann. Ich muss doch bestimmt die Variablen mit syms erstmal als Variable deklarieren oder so ?? 
|
|
|
|
| bensko |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 43
|
|
|
|
Anmeldedatum: 01.10.13
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 19:01
Titel: Re: Gleichungssystem lösen
|
|
Edit zum ersten Beitrag
| bensko hat Folgendes geschrieben: |
Es müsste ja eigentlich u1 = -2 ; u3 = 2 und F1 = -4 herauskommen
!! |
u1 = -2,5 ; u3 = 2 und F1 = -5
|
|
|
|
| bensko |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 43
|
|
|
|
Anmeldedatum: 01.10.13
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 19:45
Titel:
|
|
Hallo Leute ich habe das erstmal so gelöst. Evtl kann mir jm sagen wie man das schöner machen kann. Hier meine Lösung falls es jm interessiert
------------------------------------
% Verschiebung
ug = [u1;0;u3];
% Gleichungssystem lösen
eq1 = kg(1:1 , 1:3)*ug -F1;
eq2 = kg(2:2 , 1:3)*ug -F2;
eq3 = kg(3:3 , 1:3)*ug -F3;
GLS1= solve(eq1,eq2,eq3, F1 , u1, u3 );
GLS2 = [simplify(GLS1.F1); simplify(GLS1.u1);simplify(GLS1.u3)];
% Verschiebungsvektor
u = [GLS2(2:2);u2;GLS2(3:3)];
% Kraftvektor
F = [GLS2(1:1);F2;F3];
-------------------------------------------
Wenn jm eine bessere Lösung kennt, wäre es super wenn ich es auch wissen könnte. Schönen Abend noch zusammen
Bensko
|
|
|
|
| bensko |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 43
|
|
|
|
Anmeldedatum: 01.10.13
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 19:47
Titel:
|
|
Hallo Leute ich habe das erstmal so gelöst. Evtl kann mir jm sagen wie man das schöner machen kann. Hier meine Lösung falls es jm interessiert
Wenn jm eine bessere Lösung kennt, wäre es super wenn ich es auch wissen könnte. Schönen Abend noch zusammen
Bensko
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 19:50
Titel:
|
|
Hallo,
es sollte auch klappen, wenn du die Variablen definierst und dann die Gleichungen in Matrix-Vektor-Form angibst.
Grüße,
Harald
|
|
|
|
| Gast |
|
|
Beiträge: ---
|
|
|
|
Anmeldedatum: ---
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 20:08
Titel:
|
|
|
Wie würde das dann aussehen? Und mit welchem Befehl kann ich es dann lösen? Viele grüße
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 12.04.2015, 21:12
Titel:
|
|
Hallo,
mit solve. Wenn es ein anderer Befehl wäre, hätte ich das dazu gesagt.
Grüße,
Harald
|
|
|
|
| bensko |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 43
|
|
|
|
Anmeldedatum: 01.10.13
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 21.03.2016, 23:04
Titel:
|
|
Soeben bemerkt das ich mich gar nicht bedankt habe. Aber lieber spät als nie oder?
Danke vielmals Harald!
|
|
|
|
|
|
|
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
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
|  RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
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.
|
|
