WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Gleichungssystem lösen

 

bensko
Forum-Anfänger

Forum-Anfänger


Beiträge: 43
Anmeldedatum: 01.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 12.04.2015, 16:42     Titel: Gleichungssystem lösen
  Antworten mit Zitat      
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!!
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister


Beiträge: 24.432
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 12.04.2015, 18:32     Titel:
  Antworten mit Zitat      
Hallo,

sollte mit fsolve (numerisch) oder solve (symbolisch) kein Problem sein.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
bensko
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 43
Anmeldedatum: 01.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 12.04.2015, 18:45     Titel:
  Antworten mit Zitat      
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 ?? Question
Private Nachricht senden Benutzer-Profile anzeigen
 
bensko
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 43
Anmeldedatum: 01.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 12.04.2015, 19:01     Titel: Re: Gleichungssystem lösen
  Antworten mit Zitat      
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
Private Nachricht senden Benutzer-Profile anzeigen
 
bensko
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 43
Anmeldedatum: 01.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 12.04.2015, 19:45     Titel:
  Antworten mit Zitat      
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
Private Nachricht senden Benutzer-Profile anzeigen
 
bensko
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 43
Anmeldedatum: 01.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 12.04.2015, 19:47     Titel:
  Antworten mit Zitat      
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


Code:
% 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
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.432
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 12.04.2015, 19:50     Titel:
  Antworten mit Zitat      
Hallo,

es sollte auch klappen, wenn du die Variablen definierst und dann die Gleichungen in Matrix-Vektor-Form angibst.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Gast



Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 12.04.2015, 20:08     Titel:
  Antworten mit Zitat      
Wie würde das dann aussehen? Und mit welchem Befehl kann ich es dann lösen? Viele grüße
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.432
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 12.04.2015, 21:12     Titel:
  Antworten mit Zitat      
Hallo,

mit solve. Wenn es ein anderer Befehl wäre, hätte ich das dazu gesagt.

Code:
syms F1 u1 u3
F = [F1; 3; 2];
K = [2 -2 0; -2 3 -1; 0 -1 1];
u = [u1; 0; u3];
s = solve(F - K * u)


Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
bensko
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 43
Anmeldedatum: 01.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 21.03.2016, 23:04     Titel:
  Antworten mit Zitat      
Soeben bemerkt das ich mich gar nicht bedankt habe. Aber lieber spät als nie oder?

Danke vielmals Harald!
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

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 | goMatlab RSS Button 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.