Herezh++

Bonjour Julien,
concernant la sensibilité de la convergence :
- cela provient du calcul de l'hystérésis.
- en regardant les sorties de commentaire et les algo, j'ai l'impression que cela provient de certain tests qui sont trop draconiens. Mais compte tenu de l'imbrication de l'ensemble des calculs... ou plus simplement de la complexité de l'ensemble, ce n'est pas facile de fixer une stratégie d'où beaucoup de paramètres de réglage possibles que j'ai mis à disposition.
a) Dans le cas de l'utilisation de R-K la résolution s'effectue en explicite et l'algo contrôle la prévision d'erreur.
b) au départ
(cf. void Hysteresis3D::CalculContrainte_tdt(double & phi_2_dt,Tableau<double>& indicateurs_resolution),
...
case 5: case 6:// cas d'une résolution par intégration explicite par du kutta
...)
on commence par calculer la dérivée initiale, puis on appelle alg_edp.Pilotage_kutta qui appelle Runge_Kutta_step45 qui fait une prédiction explicite d'une valeur de contrainte à t+2./7.*delta t et ensuite appelle Hysteresis3D::Sigma_point qui a pour objectif de calculer la nouvelle dérivée au point prédit.
Or, avant de calculer la dérivée, Hysteresis3D::Sigma_point vérifie que la valeur de la nouvelle contrainte est bien inférieur à la limite de saturation. sinon -> erreur , et c'est cette vérif qui ne passe pas.
Lorsque l'on met un niveau de commentaire de 6, on voit que qu'il s'agit d'un pouillem ...
Mais on peut relâcher cette vérif:
cf. le code:
//modif : 1 juin 2014
// on vérifie que l'amplitude de la contrainte transmise, n'est pas supérieure à la saturation
// on fait la vérification sur le delta, car si on est sur une branche non initiale, on peut très bien
// dépasser le Qzero, localement, mais ce sera invalidé ensuite avec l'algo de coincidence
if ( QdeltaSigma > (Qzero*depassement_Q0*wBase)) {erreur = 1;
if (Permet_affichage() > 5)
cout << "\n QdeltaSigma =" << QdeltaSigma << " au lieu de "<< (Qzero*depassement_Q0*wBase)
<< " Hysteresis3D::Sigma_point " << endl ;
}
else // sinon on continue { erreur = 0;

On voit qu'il y a le paramètre : depassement_Q0 qui par défaut == 1

Donc j'ai fait un essai avec dans les paramètres de réglage:
depassement_Q0_ 1.0001

en me disant que c'était cohérent avec le calcul qui est effectué avec une erreur absolue de 1.e-4, cf. la mise en données:
cas_kutta_ 5 erreurAbsolue_ 1.e-4 erreurRelative_ 1.e-3#5

Dans ce cas, le calcul s'effectue (sur tes fichiers déposés) directement sans pb.

Donc j'ai l'impression que le paramètre depassement_Q0_ est intéressant à tester quand il y a un pb.
Il faut aussi augmenter le niveau de commentaire (en débug de préférance) pour essayer de comprendre ...
Enfin, le code est bien documenté au niveau des erreurs, ça peut-être utile de s'y référer.

Si tu peux essayer sur tes simul et me faire un retour, je suis preneur

concernant l'affichage des infos, c'est cohérent avec ce que je voulais afficher, car le fait de continuer l'algo via suite info, correspond finalement à un nouveau calcul.
Cependant j'ai rajouté la possibilité de supprimer l'affichage via l'ajout aux paramètres spécifiques d'algorithme d'un niveau d'affichage de -1
Cela fonctionne également avec un paramètre > 0 mais dans ce cas il y a affichage du maxi de ddl pour l'itération (donc une info de plus)