Code_Aster calcul d'endommagement Manuel utilisateur

Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
Date : 28/02/2013 Page : 1/9
Clé : U2.05.06
Révision : 10597
Réalisation de calculs d'endommagement en quasistatique (rupture fragile)
Résumé :
Ce document a pour but de conseiller l'utilisateur de Code_Aster pour la réalisation de simulation à l'aide de
modèles d'endommagement et notamment des techniques à mettre en œuvre pour avoir une solution de bonne
qualité et de manière robuste.
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
Date : 28/02/2013 Page : 2/9
Clé : U2.05.06
Révision : 10597
Table des matières
1 Introduction et mise en garde..............................................................................................................3
2 Choix du modèle d'endommagement..................................................................................................3
3 Le maillage..........................................................................................................................................5
4 La résolution du problème avec STAT_NON_LINE.............................................................................5
4.1 Newton-Raphson simple (METHODE='NEWTON')......................................................................5
4.2 Newton-Raphson + recherche linéaire..........................................................................................6
4.3 Newton-Raphson + recherche linéaire méthode mixte.................................................................6
4.4 Pilotage ........................................................................................................................................7
4.4.1 Préliminaire valable pour les 2 types de pilotages ..............................................................7
4.4.2 Conseils pour la mise en œuvre du pilotage PRED_ELAS..................................................8
4.4.3 Conseils pour la mise en œuvre du pilotage DEFORMATION............................................8
4.5 La dynamique...............................................................................................................................9
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
1
Date : 28/02/2013 Page : 3/9
Clé : U2.05.06
Révision : 10597
Introduction et mise en garde
Pour étudier la dégradation d'une structure, plusieurs outils sont disponibles et notamment (par niveau
de complexité et temps calcul croissants) :
•la mécanique de la rupture avec des critères d'amorçage basés sur des critères en contrainte ou
le taux de restitution d'énergie ;
•les modèles de zones cohésives qui permettent d'étudier l'amorçage et la propagation de
fissures lorsqu'on connaît a priori le trajet potentiel de fissuration ;
•les modèles d'endommagement qui permettent de modéliser la dégradation du matériau de
façon diffuse ou localisée.
Des guides de mise en œuvre dans Code_Aster existent pour les 2 premières familles de méthode. Il
s'agit des documents [U2.05.01] : « Domaine d'application des opérateurs de mécanique de la rupture
et conseils d'utilisation » et [U2.05.07] : « Notice d'utilisation des éléments cohésifs ». Ce document
vise à compléter l'offre des guides d'utilisation en se concentrant sur la mise en oeuvre des calculs
d'endommagement en quasi-statique.
Toutefois, il faut préciser que, par rapport aux autres techniques, les calculs d'endommagement sont
les plus compliqués à mettre en oeuvre, les plus coûteux en temps de calcul mais aussi ceux qui sont
le moins maîtrisés autant du point de vue de la robustesse que du point de vue de la qualité du
résultat, notamment parce que :
•les problèmes d'endommagement n'ont en général pas une solution unique à l'échelle de la
structure ;
•les solutions obtenues avec les modèles d'endommagement locaux dépendent fortement du
maillage ;
•on ne sait pas encore proposer de moyens d'identifier les longueurs caractéristiques des
matériaux, et ces longueurs ne sont pas transposables d'une méthode de régularisation à
l'autre ;
•la gestion de mailles cassées posent souvent problème notamment avec les modèles
régularisées. Selon les modélisations, on peut observer des élargissements non physiques
des bandes endommagés plutôt que l'apparition de multi-fissuration ou de très fortes
difficultés de convergence.
Il convient donc de toujours avoir un regard critique sur les résultats obtenus.
2
Choix du modèle d'endommagement
Pour représenter la dégradation d'un matériau voire sa rupture, l'une des méthodes possibles (notamment
lorsqu'on ne connaît pas le mode de ruine) est d'utiliser une loi de comportement « adoucissante », c'est-à-dire
telle que, une fois passé un seuil (en contrainte ou en déformation), la contrainte diminue lorsque la déformation
augmente. Ce type de comportement peut être obtenu :
•en introduisant une variable d'endommagement D comprise entre 0 et 1 (= mécanique de
l'endommagement telle qu'introduite par Kachanov ou Lemaitre) : lois ENDO_FRAGILE, ENDO_ISOT_BETON,
ENDO_ORTH_BETON, MAZARS dans Code_Aster ;
•en utilisant des modèles de plasticité avec un « écrouissage » négatif comme par exemple les lois
BETON_DOUBLE_DP, DRUCK_PRAGER, mais aussi les différentes lois de sol disponibles dans Code_Aster ;
•en introduisant une variable telle que la porosité couplée à de la plasticité pour la rupture ductile : loi de
ROUSSELIER dans Code_Aster.
Quelle que soit la méthode employée, toutes ces lois adoucissantes « locales » conduisent à la même
difficulté : il arrive un moment où le problème devient mal posé et la solution obtenue devient fortement
dépendante du maillage. En effet, l'endommagement se localise dans une bande ayant pour épaisseur un seul
élément (d'où une énergie de fissuration qui tend vers zéro quand on raffine le maillage) et il y a dépendance du
« chemin » de fissuration à la topologie du maillage. La pertinence de ces calculs (bien qu'encore assez
répandus dans la littérature) est donc très contestable, et le fait de faire dépendre l'énergie dissipée de la taille
de maille n'apporte qu'une solution très partielle.
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
Date : 28/02/2013 Page : 4/9
Clé : U2.05.06
Révision : 10597
Pour contourner ce problème, il est couramment admis qu'il faut introduire dans le problème à résoudre une
longueur caractéristique, qui va contrôler l'épaisseur de la zone endommagée indépendamment du maillage et
permettre à nouveau d'avoir convergence de la solution lorsqu'on raffine le maillage.
Plusieurs méthodes de régularisation ont été mises au point ces dernières années qui ont chacune leurs
avantages et leurs limites. Dans Code_Aster, les méthodes disponibles s'appuient toutes sur une notion de
gradient (par opposition aux méthodes de type intégrales, largement répandues dans la littérature également).
On distingue les méthodes qui introduisent :
•le gradient de la variable d'endommagement (modèle GRAD_VARI, confer documentation [R5.04.01]) ;
•une déformation régularisée calculée à partir de la déformation et de son laplacien (modèle GRAD_EPSI, cf.
doc [R5.04.02]) ;
•les modélisations second gradient et second gradient de dilatation qui introduisent une énergie dépendant
totalement, ou en partie, des composantes du gradient de la déformation (modèles 2DG et DIL, cf. doc
[R5.04.03]) ;
•les modélisations à gradient de gonflement (modèle INCO_GD, cf. doc [R3.06.08]).
Dans tous les cas, le principe est de pénaliser d'un point de vue énergétique la localisation de
l'endommagement.
Le tableau ci-dessous récapitule pour les différentes lois de comportement, quelle modélisation régularisée est
disponible (et valide) dans Code_Aster.
Loi de comportement
Modélisation
ENDO_FRAGILE
GRAD_EPSI /
GRAD_VARI /
2DG
ENDO_ISOT_BETON
GRAD_EPSI /
GRAD_VARI /
2DG
MAZARS
GRAD_EPSI /
2DG
ENDO_ORTH_BETON
GRAD_EPSI /
2DG
DRUCK_PRAGER
2DG/DIL
ROUSSELIER
INCO_GD
Tableau 1: correspondance loi d'endommagement/modélisation non-locale
Remarques :
1) Il n'existe pas de version non-locale du modèle BETON_DOUBLE_DP. La version implantée inclut
toutefois une régularisation de type Hillerborg pour éviter le problème de l'énergie qui tend vers zéro
quand la taille des éléments tend vers zéro.
2) Toutes les lois de comportement peuvent être utilisées avec les modélisations 2DG et DIL. Cependant,
la modélisation 2DG n'a pour l'instant été utilisée qu'avec les lois de sols et la modélisation DIL, n'a
de sens que pour régulariser l' « endommagement volumique » : elle est donc bien adaptée pour
traiter le cas des matériaux dilatants et donc des sols. Outre DRUCK_PRAGER, il est donc possible
d'utiliser les lois CAM_CLAY (qui est un cas particulier de la loi de Hujeux) et HUJEUX (cf. thèse
d'Alexandre Foucault et CR-AMA-09-154), VISC_DRUCK_PRAG (voir note H-T64-2009-03498) et
LETK, BARCELONE, CJS, HOEK-BROWN ; mais le retour d'expérience avec ces lois est encore
faible.
3) La modélisation INCO_GD est pour l'instant mal maîtrisée, et doit faire l'objet d'études et de
développement complémentaires. On déconseille son utilisation.
4) Le coût CPU des modélisations régularisées est important, d'une part parce qu'elles nécessitent de
mailler assez finement (voir §5) et d'autre parce qu'elles introduisent des degrés de liberté
supplémentaires rendant les matrices à inverser beaucoup plus grandes et moins creuses.
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
3
Date : 28/02/2013 Page : 5/9
Clé : U2.05.06
Révision : 10597
Le maillage
Pour un calcul mécanique avec endommagement local, on peut utiliser indifféremment des maillages linéaires
ou quadratiques. Par contre, il est conseillé d'avoir des tailles de mailles assez homogènes pour que la
dissipation le soit également (rappel : la dissipation en local est liée à la taille de la maille). Il est également
conseillé d'avoir procédé à l'identification des paramètres post-pic sur des essais utilisant à peu près les mêmes
tailles de maille.
La plupart des modèles non-locaux s'appuient sur des maillages quadratiques. Le tableau suivant récapitule les
éléments disponibles pour chaque modélisation.
2D
AXI
3D
GRAD_EPSI
TRIA6
TRIA3
QUAD8
-
TETRA10
TETRA4
HEXA20
PENTA15
PYRAM13
GRAD_VARI
TRIA6
QUAD8
TRIA6
QUAD8
TETRA10
HEXA20
PENTA15
PYRAM13
2DG
TRIA7
QUAD9
-
-
DIL
TRIA7
QUAD9
TRIA6
QUAD8
-
-*
INCO_GD
TRIA6
QUAD8
TRIA6
QUAD8
TETRA10
HEXA20
Tableau 2: mailles disponibles en fonction des modélisations
Concernant la finesse du maillage, on rappelle que pour que la régularisation du problème soit effective, il faut
au moins trois mailles voire dix mailles (pour la modélisation GRAD_VARI) pour une bande d'endommagement.
Ainsi, si la taille de la zone endommagée fait 1 cm , il faut que dans cette zone, la taille des mailles soient
inférieures à 3 mm voire 1 mm pour GRAD_VARI.
4
La résolution du problème avec STAT_NON_LINE
Le but de ce chapitre est de donner des conseils sur la manière de résoudre les problèmes adoucissants avec
STAT_NON_LINE. On propose ici une certaine gradation des outils en fonction de leur originalité, de leur
efficacité, de leurs difficultés de mise en œuvre,... permettant à un nouvel utilisateur de se faire sa propre
expérience. En fonction de la loi de comportement utilisée, de la modélisation, de la structure, il pourra être plus
judicieux d'utiliser directement l'une ou l'autre des méthodes (notamment le pilotage ou la recherche linéaire
mixte).
4.1
Newton-Raphson simple (METHODE='NEWTON')
Il s'agit de la méthode la plus simple à utiliser et qui garantit la qualité de la solution obtenue, à condition bien
sûr d'utiliser un critère de convergence suffisamment petit. Cependant :
(i)les problèmes étant adoucissants, il est recommandé de travailler en déplacement imposé plutôt qu'en effort
imposé. Sinon, dès que l'effort global à appliquer tant à diminuer, il n'y aura plus aucun moyen de converger.
*
La version 3D sera disponible en version 10
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
Date : 28/02/2013 Page : 6/9
Clé : U2.05.06
Révision : 10597
(ii)Dans le cas de matériau fragile, on observe souvent une instabilité de la solution, c'est-à-dire une rupture
brutale du matériau (snap-back dans la réponse globale force-déplacement). Dans ce cas, il y a peu de chance
que Newton seul franchisse cette instabilité (sauf parfois dans le cas de petites instabilités ou du passage à
l'état totalement cassé), et il convient de basculer vers la recherche linéaire mixte ou surtout le pilotage.
F
U
Figure 1: exemple de réponse instable
On recommande en général, d'activer le découpage automatique du pas de temps (plusieurs niveaux, voir la
commande DEFI_LIST_INST) et d'utiliser la matrice tangente (réactualisée pour les lois types ENDO_**, pas
forcément réactualisée pour Mazars) mais de permettre éventuellement le basculement vers la matrice sécante
lorsqu'on a redécoupé beaucoup le pas de temps (mot-clé PAS_MINI_ELAS). Une fois que la structure est
fortement endommagée, la matrice sécante est parfois plus efficace. Signalons que pour les lois de sols, il est
inutile de réactualiser la matrice sécante, car on utilise la matrice élastique (REAC_ITER_ELAS=0).
Avec la matrice tangente, on peut autoriser entre 10 et 30 itérations de Newton selon la loi de comportement
utilisé (ITER_GLOB_MAXI). En revanche, avec la matrice sécante, la convergence est beaucoup plus lente et il
faut autoriser beaucoup plus d'itérations (ITER_GLOB_ELAS≥ 100).
4.2
Newton-Raphson + recherche linéaire
Dans certains cas, la recherche linéaire classique (mot-clé RECH_LINEAIRE, METHODE='CORDE') peut aider à
la convergence. Mais la méthode MIXTE est nettement plus robuste.
4.3
Newton-Raphson + recherche linéaire méthode mixte
La recherche linéaire méthode mixte (mot-clé RECH_LINEAIRE, METHODE='MIXTE') est capable d'aller
chercher des solutions plus loin que le bassin d'attraction de Newton et va ainsi permettre de franchir des snapback. Pour éviter d'aller chercher des solutions trop éloignées (les problèmes d'endommagement n'admettent
pas une solution unique), on conseille :
•d'utiliser de petits pas de temps,
•de ne pas activer le redécoupage du pas de temps (ou très peu), car en présence d'une grande instabilité la
taille du pas de temps importe peu
•d'autoriser un grand nombre d'itérations de Newton (ITER_GLOB_MAXI=1000),
•d'utiliser la matrice tangente réactualisée (REAC_ITER=1),
•d'autoriser un grand nombre d'itérations de recherche linéaire pour laisser à la méthode l'opportunité de
trouver la solution dans la direction de descente (ITER_LINE_MAXI=50)
C'est en général la méthode la plus robuste, mais elle peut s'avérer relativement coûteuse, puisqu'on autorise
beaucoup d'itérations. Par contre, elle n'est en général plus très efficace quand la structure est totalement
cassée. Dans ce cas, elle a tendance à faire beaucoup d'itérations et peut parfois augmenter artificiellement la
taille de la bande de localisation.
En cas de grands snap-back, correspondant à une rupture brutale de l'éprouvette, identifiable sur la courbe
force-déplacement par une diminution brutale de l'effort appliqué, il est conseillé d'utiliser le pilotage, qui va
chercher une solution continue afin de s'assurer de la solution de la qualité.
4.4
Pilotage
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
Date : 28/02/2013 Page : 7/9
Clé : U2.05.06
Révision : 10597
Le pilotage ou méthode de continuation est une méthode qui permet de suivre la réponse d'une structure en cas
d'instabilité (mot-clé PILOTAGE, confer documentation [R5.03.80]). Pour ce faire, l'intensité du chargement
devient une nouvelle inconnue du problème. En conséquence, cette méthode n'est pas applicable lorsque le
problème dépend explicitement du temps (endommagement couplé avec du fluage, chargement thermique,
etc.).
Deux modes de pilotage adaptés aux problèmes d'endommagement existent dans Code_Aster. Il s'agit d'une
part du pilotage par incrément de déformation (TYPE='DEFORMATION'), et d'autre part, du pilotage par
prédiction élastique (TYPE='PRED_ELAS'). Le pilotage par prédiction élastique est en général plus performant,
mais il n'est pas générique (contrairement au pilotage par déformation) ; il est disponible pour les lois
ENDO_FRAGILE, ENDO_ISOT_BETON et ENDO_ORTH_BETON.
4.4.1
Préliminaire valable pour les 2 types de pilotages
Lorsque vous décidez de piloter un chargement, vous ne maîtrisez plus son évolution au cours du temps. En
effet, l'intensité du chargement devient une inconnue du problème : elle peut augmenter ou diminuer au cours
des incréments. Vous pouvez néanmoins fixer la rapidité avec laquelle il va évoluer, grâce aux incréments de
temps et au paramètre COEF_MULT (voir plus loin).
Vous ne pouvez piloter qu'un chargement ou alors plusieurs chargements qui sont proportionnels et donc
définis dans la même occurrence de AFFE_CHAR_MECA et multipliée par la même fonction multiplicatrice. (ex.
traction + cisaillement simultanés). Ainsi, si les efforts sont décomposés en deux termes : les efforts fixes F fixe
et les efforts pilotés F pilo , le chargement vaudra F =F fixe F pilo où  sera l'intensité de l'effort piloté.
Admettons que le chargement soit composé des conditions aux limites et d'un effort de traction.
TRAC = AFFE_CHAR_MECA(DDL_IMPO=(_F(GROUP_MA='HAUT', DY = 0.0002)
FCT = DEFI_FONCTION( VALE = 0. , 0., 1., 1.)
Si TRAC est le chargement à piloter, dans STAT_NON_LINE, il suffit de remplacer :
EXCIT=_F(CHARGE= TRAC, FONC_MULT=FCT)
par
EXCIT=_F(CHARGE= TRAC, TYPE_CHARGE='FIXE_PILO').
Au cours des pas de temps, vous pouvez suivre l'intensité du chargement grâce à la valeur ETA_PILOTAGE
indiquée dans les tableaux de convergence de STAT_NON_LINE (voir exemple ci-dessous). Dans notre cas, à
la fin du pas de temps considéré, le déplacement vaudra :
dy=0,0002×eta_pilotage=0,0002×8,66245×10
-3
-------------------------------------------------------------------------------------|
ITERATIONS
|
RESIDU
|
RESIDU
|
PARAMETRE
|
OPTION
|
|
NEWTON
|
RELATIF
|
ABSOLU
|
PILOTAGE
|
ASSEMBLAGE
|
|
| RESI_GLOB_RELA | RESI_GLOB_MAXI | ETA_PILOTAGE |
|
-------------------------------------------------------------------------------------|
0
X | 6.19721E-02 X | 1.61697E+02
| 9.58385E-03
|TANGENTE
|
|
1
X | 1.35762E-04 X | 3.20143E-01
| 8.66483E-03
|TANGENTE
|
|
2
| 2.15218E-07
| 5.07372E-04
| 8.66245E-03
|TANGENTE
|
--------------------------------------------------------------------------------------
Vous pouvez contrôler les bornes du chargement en fixant dans le mot-clé PILOTAGE, les paramètres
ETA_PILO_MAX et ETA_PILO_MIN. Quand ces valeurs maximales ou minimales sont atteintes, le calcul
s'arrêtera, quel que soit l'instant. ETA_PILO_MIN est surtout utile quand on travaille à effort imposé pour arrêter
le calcul quand l'effort devient très petit (parce que la structure est cassée) et ETA_PILO_MAX quand on
travaille à déplacement imposé, lorsqu'on a atteint le déplacement souhaité.
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
Date : 28/02/2013 Page : 8/9
Clé : U2.05.06
Révision : 10597
On peut également renseigner les valeurs ETA_PILO_R_MIN et ETA_PILO_R_MAX qui fixe les valeurs
physiques pour ETA. ETA_PILO_R_MIN est notamment indispensable pour ENDO_FRAGILE qui présente le
même comportement en traction et en compression sous peine de sauter d'une solution à l'autre.
Enfin, il peut être intéressant de préciser le GROUP_MA sur lequel effectuer le pilotage. En effet, si
l'endommagement est localisé, cela peut permettre de limiter les temps de calculs. En revanche, c'est
indispensable dans le cas où il y a plusieurs lois de comportement dans le modèle, dont certaines non
« pilotables » (pour pilotage type PRED_ELAS).
En général, le calcul avec pilotage demande plus de pas de temps qu'un calcul direct, car on limite l'évolution
de l'endommagement. Si le calcul s'est arrêté parce qu'il a atteint la fin de la liste d'instant, et pas le
chargement, il suffit de « prolonger » la liste d'instant initiale au moment de la POURSUITE du calcul.
4.4.2
Conseils pour la mise en œuvre du pilotage PRED_ELAS.
Pour les lois ENDO_FRAGILE et ENDO_ISOT_BETON, le coefficient COEF_MULT est relié au pas de temps et à la
variation d'endommagement minimale recherchée :
COEF_MULT=
t
D
Ainsi, si vous souhaitez qu'il y ait au moins un point qui voit son endommagement croître de 20%, et que vos
pas de temps valent
 t=1 , vous devez définir COEF_MULT=
1.
=5 . Ainsi plus COEF_MULT est petit,
0.2
plus l'endommagement progresse rapidement.
Pour les autres lois, le coefficient COEF_MULT assure que, au moins un point de Gauss sorte du seuil
d'élasticité
f
linéarisé
COEF_MULT=
pred _ elas
d'une
t
COEF_MULT
quantité
;
on
a
donc
t
MAX  f pred_elas . Dans ce cas, la valeur du paramètre à fournir n'est pas évidente et
pts de Gauss
dépend des lois. On conseille de mettre 1 pour commencer. En fonction de la progression de l'endommagement
observé, on pourra accélérer le phénomène en diminuant COEF_MULT ou le ralentir en augmentant
COEF_MULT. On peut aussi préférer de manière équivalente modifier la liste des pas temps.
On conseille en général d'opter pour SELECTION='RESIDU',
Lorsqu'on pilote le chargement, il ne faut pas hésiter à autoriser beaucoup de redécoupage (en utilisant la
commande DEFI_LIST_INST) car certains passages sont délicats.
4.4.3
Conseils pour la mise en œuvre du pilotage DEFORMATION.
Pour le pilotage en déformation, il est obligatoire de faire au moins un incrément de charge sans le pilotage. On
conseille en général de ne basculer avec le pilotage que lorsqu'on ne converge plus avec les autres méthodes .
L'idée ici est d'assurer qu'au moins un point de la structure voit sa déformation progresser de façon monotone :
COEF_MULT× MAX 
pts de Gauss
.
.

=
t
∥ -∥
La valeur de COEF_MULT n'est pas toujours facile à calibrer, mais c'est en général des valeurs importantes
puisqu'inversement proportionnelle à l'incrément de déformation (exemple pour des pas de 0,1 s et un
0,1
accroissement de déformation de 10-6 , COEF_MULT= 10-6 =100 000 ). Après un premier test, on ajustera
si besoin la valeur de ce coefficient en augmentant la valeur pour charger moins vite, et en la diminuant pour
endommager plus vite. On peut aussi préférer agir sur la liste des pas de temps.
4.5
La dynamique
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Réalisation de calculs d'endommagement en quasi-st[...]
Responsable : Kyrylo KAZYMYRENKO
Date : 28/02/2013 Page : 9/9
Clé : U2.05.06
Révision : 10597
En dernier recours, on peut également essayer de se lancer dans un calcul dynamique. Dans certains cas, cela
peut apporter des solutions mais c'est à utiliser avec beaucoup de précautions. Une documentation U2 a
spécialement été rédigée pour rassembler les conseils, il s'agit de la doc [U2.04.07], qu'il est indispensable de
consulter avant de se lancer dans une simulation en dynamique (pour résoudre un problème quasi-statique).
Remarque :
Depuis la version 10 de Code_Aster, deux méthodes seront également disponibles pour l'utilisateur : la gestion
automatique du pas de temps (qui peut se révéler une aide intéressante) et la méthode IMPLEX pour les lois de
comportement ENDO_FRAGILE et ENDO_ISOT_BETON mais uniquement en version locale.
Manuel d'utilisation
Fascicule u2.05 : Mécanique de la rupture et de l'endommagement
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
">