Code_Aster
Titre :
Réalisation d'une étude de génie civil sous charge[...]
Responsable :
Alexandre FOUCAULT
Version default
Date :
12/11/2014
Page :
38/48
Clé :
U2.06.10
Révision :
12711
7 Résolution du calcul dynamique transitoire non linéaire
7.1
Introduction
Dans ce chapitre, on se place dans le cadre d’un calcul en petites perturbations, sans choc. Les seules non linéarités sont d’origine matériau. Ceci correspond aux études de génie civil les plus fréquentes.
On pourrait cependant modéliser des non linéarités géométriques (grands déplacements, grandes déformations) ou du contact-frottement.
L’ensemble des mots-clés cités dans ce chapitre sont attachés à l’opérateur DYNA_NON_LINE. On se réfèrera à [U4.53.01] et [U4.51.03] pour plus de détails.
7.2
Aide à la convergence du calcul
7.2.1 Évolution temporelle du chargement
La commande DEFI_LIST_INST permet de réaliser un redécoupage automatique du pas de temps lorsque l’algorithme de Newton ne converge pas.
Par ailleurs, on rappelle qu’il est nécessaire, pour un calcul donné, de tester plusieurs pas de temps
afin d’analyser la stabilité des résultats en fonction du pas de temps choisi (cf.§2.3).
7.2.2 Critères de convergence
Le mot-clé CONVERGENCE permet de définir la valeur du résidu relatif sur l’équilibre
(RESI_GLOB_RELA). Par défaut, celui-ci est à
10
−
6
. Afin de faciliter la convergence du calcul, on peut relâcher ce critère de convergence. Il faut cependant le faire avec précaution, l’utilisation d’un critère de convergence trop grand pouvant entraîner des résultats éloignés de la solution réelle. On conseille de ne pas utiliser un RESI_GLOB_RELA supérieur à
10
−
4
.
Il existe de plus un résidu absolu sur l’équilibre (RESI_GLOB_MAXI). Initialement, ce critère est utilisé lorsque le chargement et les réactions d’appui deviennent nuls (par exemple dans le cas d’une décharge totale). Dans ce cas, on passe automatiquement du critère relatif au critère absolu.
Lorsqu’on utilise un modèle de comportement endommageant, on observe dans certains cas une divergence du résidu absolu alors que le résidu relatif reste faible. C’est pourquoi, il est conseillé lors des calculs de vérifier que RESI_GLOB_MAXI reste faible. Si celui-ci atteint des valeurs élevées, on préconise de reprendre le calcul en imposant une valeur maximale au résidu global en définissant un
RESI_GLOB_MAXI pour le calcul. On conseille d’utiliser la valeur par défaut à
10
−
6
.
Sous ITER_GLOB_MAXI, on peut modifier le nombre d'itérations maximum effectué pour résoudre le problème global à chaque instant (10 par défaut). Si on voit qu’au cours des itérations de Newton, la convergence du modèle de comportement est lente, on peut augmenter la valeur de ce paramètre.
Sous ITER_GLOB_ELAS, on peut modifier le nombre d'itérations maximum effectué avec la matrice de
décharge lorsqu’on utilise le mot clé PAS_MINI_ELAS du mot clé facteur NEWTON (voir §7.2.3) pour
résoudre le problème global à chaque instant (25 par défaut).
7.2.3 Algorithme de Newton
Le mot-clé NEWTON permet de préciser la matrice utilisée pour les itérations globales de la méthode de
Newton (Figure 7.2.3-a). On peut utiliser soit la matrice élastique (MATRICE=’ELASTIQUE’), soit la
matrice tangente (MATRICE=’TANGENTE’). De plus, dans ce dernier cas on peut passer automatiquement de la matrice tangente à la matrice de décharge lorsque le pas de temps est ou devient (par le redécoupage) inférieur à un pas minimal (PAS_MINI_ELAS). Pour les modèles
d’endommagement la matrice de décharge s’identifie à la matrice sécante (Figure 7.2.3-a). Cette
option peut être utile lorsque le redécoupage automatique du pas de temps ne suffit pas à faire
Manuel d'utilisation
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.06 : Dynamique
Code_Aster
Titre :
Réalisation d'une étude de génie civil sous charge[...]
Responsable :
Alexandre FOUCAULT
Version default
Date :
12/11/2014
Page :
39/48
Clé :
U2.06.10
Révision :
12711
converger un calcul. Pour les modèles d’endommagement (béton), la matrice tangente peut devenir singulière et il peut être préférable d’utiliser la matrice de décharge pour converger.
(sécante pour les modèles d’endommagement)
Figure 7.2.3-a : description schématique des opérateurs de résolution de la méthode de Newton dans Code_Aster (pour un modèle de comportement d’endommagement).
Dans un premier temps, on préconise d’utiliser la matrice tangente. Pour optimiser la convergence, il est fortement conseillé de réactualiser cette matrice tangente le plus souvent possible, la meilleure stratégie étant, si la taille du problème le permet, de la réactualiser à toutes les itérations
(MATRICE=’TANGENTE’, REAC_ITER = 1, [U4.51.03]). Ceci est d’autant plus conseillé si on a une
modélisation utilisant la méthode de Deborst (cf. 4.2).
Si des problèmes de convergence apparaissent (matrice tangente singulière), on peut choisir de basculer sur la matrice de décharge (activation du mot-clé PAS_MINI_ELAS). Comme la convergence avec la matrice de décharge est plus lente que celle avec la matrice tangente, le mot clé
ITER_GLOB_ELAS (cf. §7.2.2) permet de définir un nombre d’itérations maximal spécifique à
l’utilisation de la matrice de décharge. Pour les modèles de comportement endommageants, la matrice de décharge (qui correspond à la matrice sécante) dépend de l’état de déformation atteint. Il est donc nécessaire de réactualiser cette matrice par l’intermédiaire du mot-clé REAC_ITER_ELAS.
Dans le cas où les problèmes de convergence subsistent, il peut être intéressant de réaliser un calcul avec la matrice élastique pour voir comment se déroulent les itérations de Newton dans ce cas.
Remarque
La recherche linéaire permet théoriquement d’améliorer la convergence de la méthode de Newton
[R5.03.01]. Cependant celle-ci n’est pas opérationnelle avec l’opérateur DYNA_NON_LINE .
7.2.4 Non linéarités matériaux
Le mot-clé COMPORTEMENT permet de définir les modèles de comportements non linéaires mis en
œuvre dans DYNA_NON_LINE. En cas de non convergence du calcul, on peut être amené à modifier
les paramètres associés à l’algorithme de De Borst (cf. §4.2, pour un calcul coque multicouches).
Manuel d'utilisation
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.06 : Dynamique
