6. Réglages de la convergence#

Le mot-clef facteur CONVERGENCE permet de contrôler les paramètres de convergence du calcul.

Le mot-clef ITER_GLOB_MAXI donne le nombre d’itérations maximum effectué pour résoudre le problème global à chaque instant. L’augmentation excessive de ce paramètre est généralement le signe d’un problème dans la modélisation ou d’une discrétisation temporelle inadéquate.

Il existe deux critères de convergence: RESI_GLOB_RELA en relatif et RESI_GLOB_MAXI en absolu. Si on ne renseigne rien, alors tout se passe comme si RESI_GLOB_RELA = 1.E-6.

Pour RESI_GLOB_MAXI=resmax, l’algorithme continue les itérations globales tant que

\[\underset{i=1, \dots , \mathrm{nbddl}}{\max} {F}^{n} > \text{resmax}\]

où \({F}^{n}\) est le résidu de l’itération n (Cf. [R5.03.01] pour plus de détails). Si cet opérande est absent, le test n’est pas effectué.

Pour RESI_GLOB_RELA=resrel, l’algorithme continue les itérations globales tant que:

\[\underset{i=1, \dots , \mathrm{nbddl}}{\max} {F}^{n} > \text{resrel} . {\max} L\]

où \(L\) est le vecteur du chargement imposé et des réactions d’appuis (Cf. [R5.03.01] pour plus de détails).

Lorsque le chargement et les réactions d’appui deviennent nuls, c’est-à-dire lorsque \(L\) est nul (par exemple dans le cas d’une décharge totale), RESI_GLOB_RELA est affiché à -1, et on essaie de passer du critère de convergence relatif RESI_GLOB_RELA au critère de convergence absolu RESI_GLOB_MAXI. Cette opération est transparente pour l’utilisateur (message d’alarme émis dans le fichier .mess). Lorsque le vecteur \(L\) redevient différent de zéro, on repasse automatiquement au critère de convergence relatif RESI_GLOB_RELA.

Toutefois, ce mécanisme de basculement ne peut pas fonctionner au premier pas de temps. En effet, pour trouver une valeur de RESI_GLOB_MAXI raisonnable de manière automatique (puisque l’utilisateur ne l’a pas renseigné), on a besoin d’avoir eu au moins un pas convergé sur un mode RESI_GLOB_RELA . Dès lors, si le chargement est nul dès le premier instant, le calcul s’arrête. L’utilisateur doit déjà alors vérifier que le chargement nul est normal du point de vue de la modélisation qu’il réalise, et si tel est le cas, trouver un autre critère de convergence (RESI_GLOB_MAXI par exemple).

Si RESI_GLOB_RELA et RESI_GLOB_MAXI sont présents tous les deux, les deux tests sont effectués.

7. Critères de convergence avancés#

Le critère RESI_REFE_RELA \(\varepsilon_{refe}\) conduit à estimer la convergence de l’algorithme de Newton de la manière suivante. A partir d’une référence, qui peut être:

  • Une contrainte SIGM_REFE

  • Une déformation EPSI_REFE pour des éléments incompressibles, des éléments de grille et de membrane;

  • Une variable interne VARI_REFE si l’on utilise des lois non locales à gradient de déformation;

  • Un flux thermique FLUX_THER_REFE dans un cas THM;

  • Deux flux hydriques FLUX_HYD1_REFE et FLUX_HYD2_REFE dans un cas HHM;

  • Un déplacement DEPL_REFE si on utilise des éléments de joint avec un comportement de type CZM;

  • Une force EFFORT_REFE et un moment MOMENT_REFE si on utilise des éléments de structure (discrets, barres, poutres ou câbles);

  • Une force linéique EFFORT_REFE et un moment linéique MOMENT_REFE si on utilise des éléments de structure 2D (plaques et coques);

On calcule une référence de résidu \({F}^{\mathit{ref}}\) (un vecteur de même longueur que le vecteur résidu). La convergence sera réalisée si et seulement si:

\[\forall i \in [1, \dots, \mathit{nbddl} ] \mid {F}_{i}^{n}\mid < \varepsilon_{refe} \, {F}_{i}^{\mathit{ref}}\]

EFFORT_REFE et MOMENT_REFE n’attendant pas des valeurs de même dimension pour les éléments de structure 1D et 2D, il est fortement déconseillé d’utiliser ce critère de convergence sur un modèle comportant ces deux types d’élément.