Discussion :

[membreComplexe12]

Salut tous,

j'ai une ODE à résoudre mais j'ai un soucis avec l'utilisation de cette fonction, voici le message d'erreur:

Code :

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

1
2
Warning: Failure at t=5.844315e-003.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.387779e-017) at time t. 
> In ode23 at 362

je ne comprends ce qu'il y a marqué mais je ne comprends pas pourquoi ni comment solutionner ceci...

Voici mon script principal (je n'ai mis que les éléments importants):

Code :

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

1
2
3
4
5
6
7
8
global theta0;
global beta0;
global cste;
K1=0.5; K2=0.4;
theta=(a*m*K1)/2;
beta=(K2*a)/2;
cste=(b*M)/(2*l);
[X{1} Y{1}] = ode23(@F10, [tempsI{1} tempsF{1}], initial{1});

=> j'utilise des CELL car ici ce que je fais je dois le faire pour plusieurs problemes différents {1}{2}{3}...

voici ma fonction:

Code :

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

1
2
3
4
5
function dYdE=F10(t,y)
global theta;
global beta;
global cste;
dYdE=theta+cste*(1/y)-beta*y;

pour info voici les ordre de grandeurs de mes variables:

Code :

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

1
2
3
4
5
6
initial{1}=1e-05 % choix que j'ai fais arbitrairement car en théorie ça doit être très proche de 0)
tempsI{1}=0.0058
tempsF{1}=0.0766
cste=2.81e+18; 
theta=3.879;  
beta=0.55;

j’espère que vous pourrez m'aider à comprendre et solutionner mon problème car je ne vois pas quoi faire du tout...

merci d'avance

[membreComplexe12]

j'ai essayé avec ODE45 et ça ne change rien, j'ai le meme probleme.

a votre avis ça vient du schema que j'utilise ou de mon probleme ?

[FR119492]

Salut!
Il y a deux explications possibles:

1. Une erreur dans la formulation de ton problème, par exemple si la fonction que tu veux intégrer présente une asymptote verticale à l'intérieur du domaine d'intégration.
2. Une erreur dans l'utilisation de Matlab.

Personnellement, je pencherais plutôt pour la première hypothèse. Pour trancher, je te suggère de reprogrammer ton problème avec un autre langage, comme le C ou le Fortran, et d'afficher ou d'imprimer les valeurs intermédiaires.
Jean-Marc Blanc

[membreComplexe12]

merci d'avoir pris le temps de répondre !

j'ai trouvé pourquoi ça plante mais pas comment résoudre ceci

=> en fait la valeur: cste=2.81e+18
pose problemes car si je mets une valeur plus faible (du meme ordre de grandeur que les autres variables) alors la methode numerique arrive à resoudre le probleme.

le probleme est que je ne sais pas comment résoudre si cste=2.18e18 (ce qui est mon cas réel )

[FR119492]

Salut!
J'ai une idée, mais pour pouvoir la concrétiser, il faudrait que tu m'indiques la signification physique et la valeur de toutes les constantes et toutes les variables qui interviennent dans ton problème, ainsi que les unités utilisées.
Jean-Marc Blanc

[membreComplexe12]

d'accord, voici:

1°) la constante "cste" est un mélange de plusieurs grandeurs caractéristiques dans un materiau:

cste=(z^2*y^2*o*j^3)/(2*p);

- z (0.3 sans unité)
- j (3 sans unité)
- y (2.6e10 Pascals)
- o (3e-10 mètre)
- p (5e-8 mètre)

2°) et voici les deux autres grandeurs

theta0=(z*y*o*j^2*K1)/2;
beta0=(K2*j)/2;

K1 et K2 n'ont pas d'unité c'est juste pour calibrer le model mais à priori je pense que c'est proche de 1

[FR119492]

Salut!
Essaie de changer d'unités, c'est-à-dire, par exemple, d'exprimer y en gigapascals, o et p en nanomètres.
Jean-Marc Blanc

[membreComplexe12]

en effet, je pense que c'est la seule solution...
merci JM

ssinon pourrais tu m'expliquer pourquoi les methodes numeriques buttes sur ce genre de probleme, j'aimerai comprendre...