Bonjour,
svp qlq'un peut m'aider a résoudre un probleme d'element fini sur matlab voici le debut de programme et l'enoncé de probleme, merci d'avance.
%preparation des données
E=10000;
nu=0.3;
C=[1 nu 0; nu 1 0; 0 0 (1-nu)/2]*E/(1-nu*nu);

Nombre_node=18;
Nombre_elem=3;

coor=[0 0;
5 0;
10 0;
0 15;
10 15;
0 30;
5 30;
10 30;
0 50;
10 50;
0 70;
5 70;
10 70;
0 85;
10 85;
0 100;
5 100;
10 100];
coor(:,1)=coor(:,1)*4;

elem = [1 3 8 6 2 5 7 4;
6 8 13 11 7 10 12 9;
11 13 18 16 12 15 17 14];

hold on;
% //plot
for ne =1:Nombre_elem
x(1)=coor(elem(ne,1),1);
y(1)=coor(elem(ne,1),2);
x(2)=coor(elem(ne,5),1);
y(2)=coor(elem(ne,5),2);
x(3)=coor(elem(ne,2),1);
y(3)=coor(elem(ne,2),2);
x(4)=coor(elem(ne,6),1);
y(4)=coor(elem(ne,6),2);
x(5)=coor(elem(ne,3),1);
y(5)=coor(elem(ne,3),2);
x(6)=coor(elem(ne,7),1);
y(6)=coor(elem(ne,7),2);
x(7)=coor(elem(ne,4),1);
y(7)=coor(elem(ne,4),2);
x(8)=coor(elem(ne,8),1);
y(8)=coor(elem(ne,8),2);
x(9)=coor(elem(ne,1),1);
y(9)=coor(elem(ne,1),2);
plot(x,y)
end

%construction matrice de rigidité globale

%conditions aux limites

%resoudre le systeme

%calcul deformations et contraintes
% calcul du Jacobien J

Configuration: Windows Vista Internet Explorer 7.0