Discussion :

[comoliv02]

Bonjour à tous,

voilà mon problème:

je sollicite un element de volume avec une certaine contrainte qui augmente en fonction du temps. Je discritise cette contrainte par petit pas pour voir le comportement de mon element de volume.

Pour déterminer l'évolution de mon élément de volume, je définis une matrice masse M initiale. Cette matrice variera au cours du temps selon les sollicitations subies par l'element.
J'ai un problème de la forme M*y'=f(t,y) (f est la sollicitation).
Je considère que ma matrice masse est constante à chaque pas de contrainte considéré.

On m'a dit d'utiliser une fonction ODE pour calculer l'état final du comportement de mon element au dernier pas de sollicitation.

Quelqu'un peut il m'expliquer "simplement" comment fonctionne ODE (je ne comprends pas trop l'explication de l'aide de matlab) et comment l'appliquer à mon probleme


Merci d'avance.

Olivier

PS: je ne sais pas si mon explication de probleme est comprehensible par tous. N'hesitez pas à me poser des questions.

[rostomus]

Bonjour,

Je ne suis pas sûr d'avoir bien compris le help, mais au grosso modo:

• On crée une fonction qui definit la matrice M en fonction de "t" et "y", on la nomme par exemple "masse"
• On crée aussi une fonction qui definit le systeme ED, on la nomme "contrainte" par exemple.
• On utilise la fonction "odeset" pour poser la proprieté "Mass" à @masse.

un exemple:

Code :

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

1
2
3
4
function M=masse(t,y)
M=[1 0 0;...
   3 1 0;...
   0 0 t+1];

Code :

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

1
2
3
4
5
function dy=contrainte(t,y);
dy = zeros(3,1); 
dy(1)=-y(1)*t;
dy(2)=-y(1)*t;
dy(3)=-y(3)*y(1)*t;

Code :

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

1
2
3
4
 
option=odeset('Mass',@masse);
[t,y]=ode45(@contrainte,[0 10],[3 3 3],option);
plot(t,y)

[comoliv02]

Bonjour,

merci je pense avoir compris.

Bonne journée

Olivier