Discussion :

[zaitIX]

bonsoir,
Je vais essayer d'expliquer mon problème clairement.
Je dois résoudre une EDO pour caractériser le comportement mécanique d'un métal subissant des chargement cyclique, on obtient un cycle d'hytérésis.
Mon système est définit comme ceci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
function  dy=f_iso(t,y)
global    Q b K n E sigma0 epoint
if abs(y(1))-Q*(1-exp(-b*y(2)))-sigma0<0
    dy(2)=0
else
    dy(2)=((abs(y(1))-Q*(1-exp(-b*y(2)))-sigma0)/K)^n
end
dy(3)=sign(y(1))*dy(2);
dy(1)=E*(epoint-dy(3));
dy=dy';
end

Puis ma fonction principale est:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
global Q  b   K n E sigma0 epoint
b=5;
Q=-150 ;
E=140000;
sigma0=200;
K=800;
n=6;
epoint=0.001;
% 1ere branche
[Tn,N]=ode45('f_iso',[0 10],[0  0  0]);
sigma1=N(:,1);
ep1=N(:,3);
e1=ep1+sigma1./E;
plot(e1,sigma1)

Ce qui me donne la première branche.
Ensuite je veux tracer 9 autres branches sachant que à partir de la deuxième ça devient cyclique.
Pour cela je fais une boucle for :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
L=length(N)
for w = 0:10
epoint=-epoint
[Tn,N]=ode45('f_iso',[10*(2*w+1) (30+20*w)],[N(L,:)]);
sigma2=N(:,1);
ep2=N(:,3);
e2=ep2+sigma2./E;
L=length(N)
end

A chaque itération Ti=Ti+20 et Tf=Tf+20 et comme condition initiale je prends la valeur du dernier vecteur précédent .
Mais je ne sais pas comment ployer pour faire en sorte d'avoir toutes les courbes sigma=f(e) sur un meme grave avec la boucle for...
De plus avec cette boucle, les valeurs vont être écraser à chaque itération donc ça va pas marcher...

merci

[zaitIX]

PLus de précisions.
La boucle marche très bien fait , mais lorsque je plot, je n'ai que le tracé de la dernière itération cependant moi je voudrais tracer toutes les courbes mais je ne sais pas comment faire...

[FLB]

Salut,
si j'ai bien compris tu cherches à sauvegarder e2 et sigma2 à chaque itération de ta boucle :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
sauvegardeE2 = zeros(L,11);
sauvegardeSigma2 = zeros(L,11)
for w = 0:10
%[...]
L=length(N);
sauvegardeE2(1:L,w+1) = e2;
sauvegardeSigma2(1:L,w+1) = sigma2;
end
plot(sauvegardeE2,sauvegardeSigma2);

[Merel]

Bonjour,

Pourquoi ne pas tout simplement utiliser un hold on après le premier plot et rajouter un plot de e2 et sigma 2 dans la boucle ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
global Q  b   K n E sigma0 epoint
b=5;
Q=-150 ;
E=140000;
sigma0=200;
K=800;
n=6;
epoint=0.001;
% 1ere branche
[Tn,N]=ode45('f_iso',[0 10],[0  0  0]);
sigma1=N(:,1);
ep1=N(:,3);
e1=ep1+sigma1./E;
plot(e1,sigma1)
hold on

L=length(N);

for w = 0:10
epoint=-epoint;
[Tn,N]=ode45('f_iso',[10*(2*w+1) (30+20*w)],[N(L,:)]);
sigma2=N(:,1);
ep2=N(:,3);
e2=ep2+sigma2./E;
L=length(N);

plot(e2,sigma2)

end

[tachmou]

Salut,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
L=length(N)
for w = 0:10
epoint=-epoint
[Tn,N]=ode45('f_iso',[10*(2*w+1) (30+20*w)],[N(L,:)]);
sigma2=N(:,1);
ep2=N(:,3);
e2=ep2+sigma2./E;
L=length(N)
figure % ouvre une nouvelle figure et sera donc appelé à chaque itération de ta boucle
% Ici tu peux donc tracer ce que tu veux dans ta nouvelle figure
end

cordialement,