.. _u4.34.03-adaptation: Mot-clé ADAPTATION ================== L'adaptation du pas de temps consiste à calculer un nouveau pas de temps :math:`\Delta {t}^{i+1}`, à partir des informations du pas de temps actuel :math:`\Delta {t}^{i}` (et des pas de temps précédents :math:`\Delta {t}^{i-1}`, :math:`\Delta {t}^{i-2}`, ...). Le but est donc de calculer :math:`\Delta {t}^{i+1} = c \cdot \Delta {t}^{i}`, où :math:`c` est un coefficient réel. Le procédé est le suivant: Pour l'occurrence :math:`k` du mot-clé ADAPTATION, on définit un événement et un mode de calcul de :math:`\Delta {t}^{i+1}`. Si l'événement est vérifié, alors on calcule :math:`{c}^{k}` avec le mode de calcul choisi. On obtient un tableau récapitulatif de ce type: .. csv-table:: "Occurrence n°1", ":math:`{c}^{1} = 4.2`" "Occurrence n°2", "événement non vérifié" "Occurrence n°3", ":math:`{c}^{3} = 1.7`" "Occurrence n°4", ":math:`{c}^{4} = 3.9`" Ensuite, on choisit :math:`c=min({c}^{k})` parmi les occurrences dont l'événement est vérifié. Dans cet exemple, le nouveau pas de temps vaudra :math:`\Delta {t}^{i+1} = 1.7 \cdot \Delta {t}^{i}` . Pour utiliser MODE_CALCUL_TPLUS='IMPLEX', il faut que ce mot clé ne soit utilisé qu'une fois. Opérande EVENEMENT ------------------ Cet opérande permet de spécifier le critère de déclenchement pour l'adaptation du pas de temps. Si EVENEMENT='SEUIL', alors l'événement n'est vérifié que si le seuil est franchi. Si EVENEMENT='TOUT_INST', alors l'événement est vérifié à tous les instants ; ce cas est obligatoire avec la méthode IMPLEX. Si EVENEMENT='AUCUN', alors l'événement n'est jamais vérifié. Opérandes NB_INCR_SEUIL, NOM_PARA, CRIT_COMP, VALE_I ---------------------------------------------------- Ces opérandes ne sont à renseigner qu'en cas d'événement de type seuil (EVENEMENT='SEUIL'). L'événement est vérifié si on a NB_INCR_SEUIL fois de suite un «feu vert», un feu vert étant caractérisé par la condition définie par NOM_PARA, CRIT_COMP et VALE_I. Par exemple si NB_INCR_SEUIL=2, NOM_PARA='NB_ITER_NEWTON', CRIT_COMP='LE' et VALE_I=5, alors l'événement est vérifié si on a fait 2 fois de suite moins de 5 itérations de Newton. Remarque Il n'y a pas de valeur par défaut pour VALE_I dans le catalogue de l'opérateur DEFI_LIST_INST. Cependant, si sa valeur n'est pas renseignée, on prendra pour VALE_I la moitié du nombre maximal d'itérations de Newton (déclaré dans les options de l'opérateur non-linéaire). .. _u4.34.03-mode_calc_tplus: Opérande MODE_CALCUL_TPLUS -------------------------- Cet opérande sert à spécifier le mode de calcul du nouveau pas de temps (en fait, du coefficient :math:`c`). Si MODE_CALCUL_TPLUS='FIXE', alors le coefficient :math:`c` est fixe. Sa valeur est déterminée grâce à l'opérande PCENT_AUGM. Le nouveau pas de temps vaudra alors :math:`\Delta {t}^{i+1} = (1+\frac{\text{PCENT_AUGM}}{100.})\Delta {t}^{i}`. Si MODE_CALCUL_TPLUS='DELTA_GRANDEUR', alors le coefficient :math:`c` est lié à la variation au cours du pas de temps d'une grandeur choisie. On utilise la formule suivante: .. math:: c = \underset{i \in N}{min}\left( \frac{{V}_{ref}}{∣\Delta {V}_{i}∣} \right) avec :math:`N` l'ensemble des noeuds ou des points de Gauss pour lesquels :math:`\Delta {V}_{i}` n'est pas nul. La valeur de référence :math:`{V}_{ref}` est choisie par l'utilisateur. Si le champ :math:`\Delta V` est identiquement nul, alors on considère que le critère (l'événement) n'est pas vérifié. Par défaut, toutes les valeurs de la composante sont examinées. Toutefois, on peut filtrer la position en employant les mots-clés GROUP_MA (pour les champs SIEF_ELGA et VARI_ELGA) ou GROUP_NO (pour le champ DEPL). Dans ce cas, l’évaluation ne porte que sur la liste de groupes de mailles ou la liste de groupes de noeuds. Si MODE_CALCUL_TPLUS='ITER_NEWTON', alors le coefficient :math:`c` est lié à la variation au cours du pas de temps du nombre d'itérations de Newton. On utilise la formule suivante: .. math:: c = \sqrt{\frac{{N}_{ref}}{{N}_{it}+1}} où :math:`{N}_{it}` est le nombre d'itérations de Newton effectuées et :math:`{N}_{ref}` une valeur de référence choisie par l'utilisateur. On rappelle que la prédiction est considérée comme l'itération de Newton n°0, ce qui signifie que si on converge dès la prédiction alors :math:`c=\sqrt{{N}_{ref}}`. Si MODE_CALCUL_TPLUS='IMPLEX', alors le coefficient :math:`c` est lié aux incréments de variables internes (déformation plastique cumulée ou endommagement) du pas précédent (voir [:ref:`R5.03.81 `]). Dans ce cas, le mot clé ADAPTATION ne doit être utilisé qu'une seule fois. Le coefficient :math:`c` sera forcément compris entre 0.5 et 1.2. De plus, par définition de la méthode IMPLEX il n'y aura aucune cause d'échec. Opérande PCENT_AUGM ------------------- Cet opérande sert à spécifier le pourcentage d'augmentation du pas de temps si MODE_CALCUL_TPLUS = 'FIXE'. Par défaut, on augmente de 100%, c'est-à-dire que l'on double la valeur du pas de temps. La valeur peut être négative et doit être strictement supérieure à -100. Opérandes VALE_REF, NOM_CHAM, NOM_CMP ------------------------------------- Ces opérandes servent à spécifier le calcul du coefficient :math:`c` d'augmentation du pas de temps si MODE_CALCUL_TPLUS='DELTA_GRANDEUR' (cf. :ref:`u4.34.03-delta_grandeur`) avec VALE_REF la valeur de référence :math:`{V}_{ref}` et NOM_CHAM le nom du champ de grandeur considéré. Opérande NB_ITER_NEWTON_REF --------------------------- Cet opérande sert à renseigner le nombre d'itérations de Newton de référence, noté :math:`{N}_{ref}` (cf. :ref:`u4.34.03-mode_calc_tplus`) dans le cas où MODE_CALCUL_TPLUS = 'ITER_NEWTON'. Opérande INFO ============= Pour le moment, cet opérande ne sert pas.