Discussion :

[sanesc]

Bonjour,
J'ai un problème j'ai pluieurs système 'poutre verticale' donc j'ai les déplacement sur plusieurs hauteurs. Je suis arrviée à représenter un graphique dynamique qui me donne l'évolution de la déformée au cours du temps, Voici ma partie de programme :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
conteneur=[fliplr(displ(:,10:16)) displ(:,2) displ(:,1)]';
prof_conteneur=[2.61 -1.18 -4.06 -5.98 -6.94 -8.84 -9.79 -10.43 -10.43];
 
figure(8)
set(figure(8),'color',[1 1 1])
for i=1:1:500
plot(conteneur(:,i),prof_conteneur,'kd--');
axis([-25 25 -10 3]);
axis equal
M(i)=getframe;
end

Jusqu'ici tout va bien, sauf que malgré que j'essaye de fixer les xmin xmax ymin et ymax ils bougent !!...
Cependant le gros problème est que je souhaite sur la même figure mettre un second profile de déformation et suivre son évolution au cours du temps pour comparer la déformée des deux systèmes type 'poutre verticale en cisaillement'....et là je n'arrive pas à en mettre 2...

est possible et comment le programmer ?

Merci d'avance pour toute aide !

[Jerome Briot]

Comme ceci :

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
conteneur = [fliplr(displ(:,10:16)) displ(:,2) displ(:,1)]';
prof_conteneur = [2.61 -1.18 -4.06 -5.98 -6.94 -8.84 -9.79 -10.43 -10.43];
 
figure(8)
set(figure(8),'color',[1 1 1])
 
p = plot(conteneur(:,1),prof_conteneur,'kd--');
 
axis([-25 25 -10 3]);
axis equal
 
M(1) = getframe;
 
for i = 2:500
 
    set(p,'xdata',conteneur(:,i));
 
    M(i) = getframe;
 
    drawnow
 
end

Pour deux courbes (c'est un exemple) :

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
conteneur = [fliplr(displ(:,10:16)) displ(:,2) displ(:,1)]';
prof_conteneur = [2.61 -1.18 -4.06 -5.98 -6.94 -8.84 -9.79 -10.43 -10.43];
 
figure(8)
set(figure(8),'color',[1 1 1])
 
p(1) = plot(conteneur(:,1),prof_conteneur,'kd--');
hold on
p(2) = plot(conteneur(:,1)-5,prof_conteneur,'kd--');
 
axis([-25 25 -10 3]);
axis equal
 
M(1) = getframe;
 
for i = 2:500
 
    set(p(1),'xdata',conteneur(:,i));
    set(p(2),'xdata',conteneur(:,i)-5);
 
    M(i) = getframe;
 
    drawnow
 
end

[sanesc]

Merci beaucoup, je viens de tester et cela marche !!!
bonne journée !