finite elemente-Startprobleme

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

ng83

Forum-Anfänger


	Beiträge: 43

	Anmeldedatum: 26.11.12

	Wohnort: Deutschland

	Version: ---

Verfasst am: 17.06.2015, 11:30 Titel: finite elemente-Startprobleme

Hallo ihr Lieben,
wir sollen das Lösen eines 2d finite elements programmieren. Als Vorarbeit haben wir den Algorithmus händisch am Bsp. eines lin. Dreieckelementes geübt.Bis dahin dachte ich, dass ich es halbwegs verstanden habe, nun macht sich aber zunehmend Verunsicherung und Verzweiflung breit...

gegeben:
1. lineares dreieckselemet
2. Gebiet ist das EInheitsquadrat [0,1]x[0,1]
3. -laplace u=0 mit
u(x_1,x_2)=x_2 auf dem rechten Rand und
u(x_1,x_2) =0 auf dem unteren Rand (Dirichletrand)
4. du/dn=-x_2 auf dem linken Rand
du/dn=x_1 auf dem oberen Rand
Neumannrand

die Aufgabe lautet: "Programmiert mal", wir haben außerdem noch die folg Dateien mitbekommen:

Code:

%function femstart( netdat )
%FEM-Hauptroutine

% Einlesen Netzdatenfile netdat

fileID = fopen(netdat);
% mfilename = febib.m1name(typ) - m-File fuer Elementmatrix Int(grad u * grad v)
typ = num2str(cell2mat(textscan(fileID,'%d','HeaderLines',1)));

%% Knotenkoordinaten
% n - Anzahl Knoten (Anzahl Freiheitsgrade)
n = cell2mat(textscan(fileID,'%d','HeaderLines',1));
% coord(n,2) - Koordinaten der Knoten
coord = cell2mat(textscan(fileID,'%f %f','HeaderLines',1,'Delimiter','\t'));

%% Elementzusammenhang
% nelem - Anzahl Elemente
nelem = cell2mat(textscan(fileID,'%d ','HeaderLines',1));
% elnode(nelem,kn+1) letzte Spalte: Materialbereichsnummer
elnode = cell2mat(textscan(fileID,'%d %d %d %d','HeaderLines',1,'Delimiter','\t'));

%% Randseiten
% nedges - Anzahl Randseiten
nedges = cell2mat(textscan(fileID,'%d ','HeaderLines',1));
% edges(nedge,kned+1) letzte Spalte: Randbereichsnummer
edges = cell2mat(textscan(fileID,'%d %d %d ','HeaderLines',1,'Delimiter','\t'));
fclose(fileID);

% Datentyp table mit (z.B.) Namen febib
febib = feminit;
% kn = febib.kn(typ) - Anzahl Knoten (Freiheitsgrade) pro Element
kn = febib.kn(typ);
% kned = febib.kned(typ) - Anzahl Knoten auf einer Elementseite
kned = febib.kned(typ);

% TODO
% Einlesen Materialdatenfile matdat
fileID = fopen(matdat);
% Koeffizienten

% rechte Seite

% Randbedingungen

%end

Funktion ohne Link?

sowie

Code:

function febib = feminit
%Initialisierung FEM-Rechnung

% Elementbibliothek

febib = table( {'lineares Dreieckselement';'quadratisches Dreieckselement';'bilineares Viereckselement'},...
[3; 6; 4], ... % Anzahl Knoten (Freiheitsgrade)
[2; 3; 2],... % Anzahl Knoten auf Randseite
{'elem01a';'elem02a';'elem11a'},... % Steifigkeitsmatrix
{'';'';''},... % Konvektionsmatrix
{'';'';''},... % Massematrix
{'symelem01a';'';''},... % Routine fuer symbolische Rechnung
'VariableNames',{'Elementtyp','kn','kned','m1name','m2name','m3nam3','symname'},...
'RowNames',{'101','102','201'} );

end

Funktion ohne Link?

Bitte, bitte helft mir. Ich weiß weder was ich mit dem Code anfangen, noch wie ich überhaupt beginnen soll.
Bin wirklich am Verzweifeln.

VG ng


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 |

RSS

Hosted by:

Copyright © 2007 - 2025 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.