Prabel
Prabel
Prabel
THSE
MODLISATION AVEC LA MTHODE X-FEM DE LA
PROPAGATION DYNAMIQUE ET DE L'ARRT DE FISSURE
DE CLIVAGE DANS UN ACIER DE CUVE REP
Prsente devant
l'Institut National des Sciences Appliques de Lyon
pour obtenir
le GRADE DE DOCTEUR
Ecole doctorale :
Mcanique, Energtique, Gnie Civil, Acoustique
Spcialit :
MECANIQUE - GENIE MECANIQUE - GENIE CIVIL
par
Benoit PRABEL
Ingnieur
Jury
ALLIX Olivier Professeur (ENS Cachan) Rapporteur
BOMPARD Philippe Professeur (Centrale Paris) Rapporteur
SMITH David Professeur (University of Bristol) Prsident du jury
GILLES Philippe Docteur (AREVA) Examinateur
COMBESCURE Alain Professeur (INSA Lyon) Directeur de thse
MARIE Stphane HDR (CEA) Examinateur
Je tiens tout d'abord remercier Alain Combescure et Stphane Marie, qui m'ont encadr durant
ces trois annes de thse. Leur passion pour la mcanique est communicative et trs stimulante.
Alain Combescure a toujours su tre disponible et son implication dans les aspects
scientifiques, humains et matriels de la thse est trs apprciable. Je remercie Stphane Marie
pour toute l'aide qu'il m'a apport au quotidien sans jamais compter son temps, pour ses
prcieux conseils, pour les nombreuses discussions que nous avons eu et les rflexions qu'elles
amenrent.
Je suis reconnaissant envers Olivier Allix et Philippe Bompard pour avoir accepter la lourde
tche de rapporteur, ainsi qu'envers David Smith, et Philippe Gilles pour leur participation au
jury.
Je remercie galement Anthony Gravouil pour m'avoir transmis sa connaissance des X-FEM,
Lionel Ferry, pour m'avoir form sur le plan exprimental, Stphane Chapuliot pour m'avoir
accueilli dans son laboratoire (et dans son club de foot).
Je pense aussi tous les membres du Lamcos, et en particulier aux doctorants : commencer
par Ludo, Vincent, Guillaume, Bertrand et Nahiene, toujours prts prendre un caf; n'oublions
pas Julien, Thomas, Tomtom ("collgue"), Stphane, Tong, David, Johann, Pierre et tous les
autres...
Je souhaite remercier les membres du SEMT avec lesquels j'ai pu discut, et l'ensemble du
LISN, pour leur bonne humeur et leur aide prcieuse : Catherine (l'efficacit administrative en
personne), Michel, Grard, Cdric, Yann, Marie-Hlne, Thierry, Tamara, Olivier, Vincent,
Tha, et les autres...
Enfin, je suis heureux de ddier ce mmoire Amlie, ma femme.
Rsum
Remerciements 3
Rsum 5
Principales notations 17
Introduction gnrale 19
8
III.5.3. tude de la propagation....................................................................
.....................121
III.6. Analyse du facis de rupture....................................................................................
..123
III.6.1. Objectifs des observations...................................................................
..................123
III.6.2. Observations........................................................................................................
...123
III.7. Conclusions.............................................................................................
.....................125
Bibliographie 179
9
Tables des figures
Fig. 1 : (a) Modes de rupture. (b) Repre local en pointe de fissure et contrainte asymptotique
dans le cas du mode I.............................................................................................................
........23
Fig. 2 : Contour pour l'intgrale J...............................................................................................
...23
Fig. 3 : Solide en traction et courbe force-ouverture correspondante...........................................24
Fig. 4 : Dominance des champs en pointe de fissure (tir de [9])................................................25
Fig. 5 : Master Curves pour l'acier 16MND5 et essais isothermes sur CT25 (tir de [17])........26
Fig. 6 : Calage partir d'essais sur CT25 du modle de Beremin pour l'acier 16MND5 la
temprature de -125C (tir de [21])................................................................................
.............27
Fig. 7 : Observations exprimentales de Kalthoff sur des prouvettes en Araldite B.................29
Fig. 8 : Application soudaine d'une pression uniforme................................................................30
Fig. 9 : Milieu fissur...................................................................................................
.................31
Fig. 10 : Partition des domaines : = 0 ..................................................................
.......34
Fig. 11 : Dfinition des contours 1 et 2 ..................................................................
.................35
Fig. 12 : Fonctions universelles de l'acier 16MND5....................................................................36
Fig. 13 : (a) Cas de Broberg. (b) Sa gnralisation par Freund....................................................39
Fig. 14 : Fonctions caractrisant l'quation de propagation..........................................................40
Fig. 15 : Corrlation entre vitesse de propagation et facteur d'intensit des contraintes Kdyn.. .42
Fig. 16 : Schmatisation des deux rgimes viscoplastiques (tir de [10])...................................42
Fig. 17 : volution de la limite d'lasticit avec la vitesse de dformation pour diffrentes
tempratures, comparaison modle exprience (tir de [55]).......................................................43
Fig. 18 : Comparaison entre le modle de Symonds-Cowper et l'exprience (tir de [56])........44
Fig. 19 : Zone active et sillage plastique laisse par la propagation d'une fissure dans un
matriau lasto-viscoplastique.....................................................................................................
..44
Fig. 20 : Ratio Gext/Gtip en fonction de la vitesse de propagation et du paramtre matriau PC
(tir de [10])....................................................................................................................
...............45
Fig. 21 : Description des essais Robertson (a) et Pellini (b).........................................................47
Fig. 22 : prouvette Compact Crack Arrest (tir de [65])............................................................47
Fig. 23 : prouvette Double Cantilever Beam.......................................................................
.......48
Fig. 24 : prouvette Single Edge Notch Bending impacte par un projectile.............................48
Fig. 25 : Exprience de Kalthoff...................................................................................
................49
Fig. 26 : prouvette Compact Tension Specimen.........................................................................
50
Fig. 27 : Schma de l'prouvette annulaire pour choc thermique utilise au LISN [21].............51
Fig. 28 : Rsultats obtenus par Hajjaj [56]...................................................................................52
Fig. 29 : Mthode du relchement de noeud..........................................................................
.......52
Fig. 30 : Mthode des lments finis espace-temps.....................................................................54
Fig. 31 : Modle de zones cohsives...........................................................................................
..55
Fig. 32 : Mthodes sans maillage (tir de [35])............................................................................56
Fig. 33 : Dfinition de la gomtrie de la fissure.........................................................................58
Fig. 34 : Stratgie d'enrichissement des noeuds...........................................................................59
Fig. 35 : Comparaison de plusieurs mthodes numriques sur le cas test de Kalthoff................60
Fig. 36 : Paramtre d'endommagement pour une prouvette flexion 4 points (tir de [116]).....61
Fig. 37 : Interface soumise au champs de vitesse V..................................................................64
Fig. 38 : (a) Mthode des "marqueurs". (b) Mthode des volumes finis.....................................65
Fig. 39 : Dfinition d'une interface depuis une fonction de niveau..............................................65
Fig. 40 : Utilisation de la mthode des marqueurs (tir de [106])................................................66
Fig. 41 : Comparaison avec la solution thorique (a), d'un schma "upwind" (b), et d'un schma
bas sur les diffrences centres (c) (tir de [106])......................................................................67
Fig. 42 : Dfinition de la gomtrie de la fissure depuis les fonctions de niveau (,).............69
Fig. 43 : Les fonctions de niveau dfinissent un repre local en pointe de fissure......................69
Fig. 44 : Stratgie d'enrichissement des noeuds lors de la propagation de fissure......................71
Fig. 45 : Utilisation de la mthode X-FEM en dynamique explicite...........................................73
Fig. 46 : Intgration par sous dcoupage, changement de position des points de Gauss.............74
Fig. 47 : Dfinition d'une interface pour reprsenter implicitement une inclusion ou un trou (tir
de [133])...............................................................................................................................
..........76
Fig. 48 : (a) Maillage de la structure. (b) Grille auxiliaire...........................................................77
Fig. 49 : (a) R-initialisation ne maintenant pas la position de l'interface. (b) R-initialisation
maintenant la position de l'interface................................................................................
.............78
Fig. 50 : Schma de propagation de fissure...............................................................
...................79
Fig. 51 : Comparaison avec la thorie (en rouge) des isovaleurs de et obtenus avec
diffrents algorithmes utilisant la fonction signe de Peng............................................................80
Fig. 52 : Comparaison avec la thorie (en rouge) des isovaleurs de et obtenus avec
diffrents algorithmes utilisant une fonction signe proportionnelle.............................................81
Fig. 53 : Comparaison avec la thorie (en rouge) des isovaleurs de et obtenus avec
l'algorithme de Gravouil vectoriel utilisant une fonction signe proportionnelle..........................82
Fig. 54 : Normes des gradients de et obtenus par l'algorithme de Gravouil vectoriel en
utilisant une fonction signe proportionnelle.......................................................................
...........82
Fig. 55 : Sous dcoupage non conforme....................................................................
...................85
Fig. 56 : Courbe de comportement et discrtisation de l'nergie de dformation........................87
Fig. 57 : volution de l'iso-zro de la fonction de niveau reprsentant une bulle de gaz............88
Fig. 58 : Dfinition du champ q...............................................................................
.....................90
12
Fig. 59 : volution du facteur d'intensit des contraintes dynamique..........................................91
Fig. 60 : Comparaison avec la thorie du facteur d'intensit des contraintes dynamique calcul
depuis l'intgrale J...............................................................................................
..........................92
Fig. 61 : (a) Gomtrie de l'exprience de Kalthoff. (b) Trajet de fissure en fonction de la
vitesse d'impact et du rayon de fond d'entaille............................................................................
..92
Fig. 62 : Rsultats de la modlisation de l'exprience de Kalthoff..............................................95
Fig. 63 : Exprience de Grgoire [76].......................................................................................
....96
Fig. 64 : Vitesse sur les bords gauche (a) et droit (b) de l'prouvette..........................................96
Fig. 65 : Longueur de fissure en fonction du temps, comparaison calcul-exprience.................97
Fig. 66 : Trajet de fissure, comparaison calcul-exprience..........................................................98
Fig. 67 : Courbes de tractions lisses (trait continu) compares aux courbes exprimentales
(symboles) pour diffrentes tempratures (tir de [17]).............................................................103
Fig. 68 : Comparaison entre la tnacit exprimentale (sur CT25) et les courbes de transition
"Master Curve" pour l'acier 16MND5................................................................
........................104
Fig. 69 : volution du niveau de contrainte avec la temprature (a), et avec la vitesse de
dformation (b).....................................................................................................
.......................105
Fig. 70 : Comparaison avec l'exprience du modle de Symonds-Cowper modifi -150C (a)
et +25C (b)...............................................................................................................................
106
Fig. 71 : Exemple du montage d'une prouvette CT..................................................................108
Fig. 72 : Jauges brins Vishay de type CPD (a) et CPA (b)......................................................109
Fig. 73 : Schma lectrique du montage.........................................................................
............109
Fig. 74 : Exemple de rsultats exprimentaux obtenus pour la ct10ov......................................110
Fig. 75 : Exemple de longueur de fissure en fonction du temps (ct10ov).................................111
Fig. 76 : Plan de l'prouvette CT25............................................................................................
.111
Fig. 77 : (a) Trajet de fissure rectiligne. (b) Trajet de fissure branche.....................................112
Fig. 78 : Longueurs de fissure l'amorage et l'arrt..............................................................113
Fig. 79 : Courbes force - ouverture du clip........................................................................
.........114
Fig. 80 : Facteur d'intensit des contraintes l'amorage pour les prouvettes CT...................115
Fig. 81 : volution de la tnacit avec l'paisseur de l'prouvette.............................................116
Fig. 82 : Dfinition de la longueur de fissure projete sur le ligament......................................116
Fig. 83 : Rsultats d'essais pour les CT ayant un trajet de fissure rectiligne.............................117
Fig. 84 : Rsultats d'essais pour les CT ayant un trajet de fissure branche..............................118
Fig. 85 : Plans des anneaux....................................................................................
.....................119
Fig. 86 : Longueur de fissure initiale a0 et finale af pour les anneaux......................................120
Fig. 87 : Courbe force-ouverture pour les anneaux en mode I pur (a), et pour les anneaux
chargs en mode mixte (b)..........................................................................................................
.120
Fig. 88 : Longueur de fissure en fonction du temps pour les anneaux en mode I pur...............121
Fig. 89 : Trajet de fissure exprimental pour les anneaux en mode mixte.................................122
13
Fig. 90 : Longueur de fissure en fonction du temps dans le cas des anneaux en mode mixte...123
Fig. 91 : Facis de rupture de l'prouvette ct10oy mi-parcours de la propagation.................124
Fig. 92 : Facis de rupture de l'prouvette ct10oy proximit du front d'arrt.........................124
Fig. 93 : Facis de rupture de l'prouvette ct10pc l'amorage.................................................125
Fig. 94 : Dfinition du solide prsentant une discontinuit (partie grise)............................129
Fig. 95 : Inconnues du modle de propagation...........................................................................131
Fig. 96 : Maillage de l'prouvette CT....................................................................
.....................131
Fig. 97 : Comparaison entre la modlisation et l'exprience de la courbe force-ouverture
pendant la mise en charge (a), et imposition des donnes exprimentales concernant l'volution
de longueur de fissure avec le temps (b)..............................................................................
.......132
Fig. 98 : volution des termes nergtiques durant la propagation pour la ct5pm....................132
Fig. 99 : Fonction de pondration w..........................................................................................
..136
Fig. 100 : volution de diffrentes mesures de la contraintes pour la ct5pm............................136
Fig. 101 : volution de diffrentes mesures de la contraintes pour la ct5pm............................137
Fig. 102 : Identification de critres dpendant de la vitesse de propagation..............................137
Fig. 103 : Identification de critre dpendant de la vitesse de dformation totale quivalente
(RKR 100m).............................................................................................................
.................138
Fig. 104 : Identification de critre dpendant de la vitesse de dformation inlastique
quivalente (Half-Disc 200m)................................................................................................
...139
Fig. 105 : Comparaison entre l'volution de la contrainte critique Ic et celle de la limite
d'lasticit y avec la vitesse de dformation inlastique quivalente.......................................139
Fig. 106 : Comparaison entre une analyse quasi-statique et une analyse dynamique concernant
l'volution de la contrainte pour la ct5pm....................................................................
...............140
Fig. 107 : Longueur de fissure en fonction du temps pour la ct5pm (a), et la ct10oy (b).
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de propagation......................................................................
.................142
Fig. 108 : Longueur de fissure en fonction du temps pour la ct5pm (a), et la ct10oy (b).
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de dformation....................................................................................
...143
Fig. 109 : Longueur de fissure en fonction du temps pour l'anneau 520uo. Comparaison avec
l'exprience des rsultats de simulation utilisant une contrainte critique dpendant de la vitesse
de propagation.........................................................................................................................
.....144
Fig. 110 : Longueur de fissure en fonction du temps pour l'anneau 520uo. Comparaison avec
l'exprience des rsultats de simulation utilisant une contrainte critique dpendant de la vitesse
de dformation........................................................................................................
.....................144
Fig. 111 : volution angulaire des contraintes d'ouverture et principale interpoles sur un cercle
de rayon de 100m obtenue pour l'anneau 520uv......................................................................146
Fig. 112 : Trajet de fissure pour l'anneau en mode mixte 520uv. Comparaison avec l'exprience
des rsultats de simulation obtenus en utilisant comme direction de propagation celle du
maximum de la contrainte circonfrentielle........................................................................
........147
Fig. 113 : Trajet de fissure pour l'anneau en mode mixte 520uv. Comparaison avec l'exprience
14
des rsultats de simulation utilisant une contrainte critique dpendant de la vitesse de
propagation.................................................................................................................................
..147
Fig. 114 : Longueur de fissure en fonction du temps pour l'anneau en mode mixte 520uv.
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de propagation......................................................................
.................148
Fig. 115 : Trajet de fissure pour l'anneau en mode mixte 520uv. Comparaison avec l'exprience
des rsultats de simulation utilisant une contrainte critique dpendant de la vitesse de
dformation................................................................................................................
..................148
Fig. 116 : Longueur de fissure en fonction du temps pour l'anneau en mode mixte 520uv.
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de dformation....................................................................................
...149
Fig. 117 : volution de la norme de la vitesse en fonction du temps.........................................151
Fig. 118 : volution de la contrainte hydrostatique en fonction du temps (dforme x10).......153
Fig. 119 : Comparaison avec l'exprience des trajets de fissure obtenus par la simulation pour
les deux anneaux chargs en mode mixte 520uv et 520uw........................................................154
Fig. 120 : Courbes force-ouverture exprimentales et obtenues par la simulation....................154
Fig. 121 : Comparaison avec l'exprience des trajets de fissure obtenus par la simulation pour
les deux anneaux chargs en mode mixte 520uv et 520uw........................................................155
Fig. 122 : Numrotation des points d'intgration pour l'lment standard QUA4 (a), et pour les
lments enrichis QU4H ou QU4F (b)......................................................................................
..165
Fig. 123 : Principe des barres de Hopkinson...............................................................................
169
Fig. 124 : (a) Effet de la temprature. (b) Effet de la vitesse de dformation...........................170
Fig. 125 : Comparaison avec l'exprience du modle de Symonds-Cowper (D=1.E8 ; p=12). 172
Fig. 126 : Comparaison avec l'exprience du modle de Symonds-Cowper modifi pour
diffrentes tempratures...................................................................................
...........................175
Fig. 127 : Illustration du phnomne d'oscillations numriques entre deux valeurs du couple
(force, dplacement)............................................................................................
........................176
Fig. 128 : Passage du temps n au temps n+1 par prdiction lastique, puis retour la surface de
charge............................................................................................................................
...............177
15
Principales notations
oprateur gradient
= oprateur de drivation temporelle
t
{ }
Ni
Ni = N i H
Ni Fk
p dformation inlastique quivalente
ensemble des inconnues au noeud i
{}
ui
qi = ai
bi , k
rp rayon plastique estim autour de la pointe de fissure
s dviateur des contraintes
T temprature
t temps
u , u , u dplacement, vitesse, acclration
w densit d'nergie
X scalaire
X vecteur
X matrice, ou tenseur du second ordre
X tenseur d'ordre 4
xi {i=1, 2, 3} coordonnes de l'espace
18
Introduction gnrale
La prsence de dfauts dans les structures affaiblit leur rsistance mcanique, peut entraner
leur ruine, et tre l'origine de catastrophes graves. Afin d'viter celles-ci, ingnieurs et
scientifiques doivent donc tre capable de prdire le comportement des structures fissures, et
d'en garantir l'intgrit.
Selon les cas, la prsence de fissures peut tre constate lors d'inspection si le composant est
accessible au contrle, ou suppose lorsque l'on souhaite valuer la sret d'un composant.
La mcanique de la rupture est la science qui tudie le comportement des structures fissures.
Elle fournit les outils ncessaires pour l'analyse de l'intgrit de structures. Il s'agit d'un
domaine scientifique relativement rcent puisqu'il est apparu aprs la premire guerre mondiale,
notamment avec les tudes de Griffith en 1920 [1] et d'Irwin en 1957 [2].
La mcanique de la rupture couvre des phnomnes physiques varis. La rupture fragile est une
rupture brutale, qui se produit alors que la structure est globalement lastique. A l'oppos, la
dchirure ductile caractrise une progression stable de la fissure pour des niveaux de
chargement et de plasticit levs. La fissuration par fatigue se produit sous des conditions de
chargement cyclique. Enfin, le fluage dcrit les phnomnes de fissuration pour lesquels les
dformations viscoplastiques interviennent.
La mcanique de la rupture trouve notamment ses applications dans les domaines
aronautiques, de construction naval, de gnie civil, et dans les activits nuclaires...
Dans le domaine nuclaire, la sret des installations est un enjeu de premier plan. C'est
pourquoi, les tudes ralises envisagent les scnarii les plus pessimistes de dgradations des
structures, et les sollicitations les plus pnalisantes.
Par exemple, le cas de l'Accident par Perte de Rfrigrant Primaire envisage la rupture du
circuit primaire de refroidissement. Dans ce cas, il est prvu d'injecter de l'eau froide ( une
temprature de l'ordre d'une dizaine de degrs) dans le coeur du racteur. La paroi interne de la
cuve du racteur chute alors brutalement de 300C environ 10C. Ce violent gradient
thermique produit des contraintes de traction sur la paroi interne de la cuve. Dans ce cadre
accidentel, l'analyse de sret prvoit l'tude de la stabilit d'un hypothtique dfaut.
Cependant, la problmatique d'extension de dure de vie des centrales REP incite
perfectionner les tudes menes, afin de prciser de nouvelles marges garantissant l'intgrit des
structures. La prise en compte d'un possible arrt de fissure aprs son amorage fragile est l'un
des lments cls. Alors qu'en gnral, on se contente de vrifier le non amorage d'un
hypothtique dfaut, on cherche alors montrer que, mme dans le cas improbable o ce dfaut
amorcerait, la fissure initie s'arrterait avant d'avoir mis en pril la fonctionnalit de la
structure.
Une telle dmonstration pousse approfondir les modles utiliss, et amliorer leur
prdictivit. En effet, la codification [3] intgre un critre d'arrt de fissure, mais celui-ci est
bas sur une analyse statique du problme, qui, bien que juge pessimiste, est discutable. De
nombreux auteurs (comme Kalthoff [4] ou Hahn et al [5]) ont d'ailleurs montr que seule une
analyse dynamique permettait une interprtation physique des phnomnes de propagation et
d'arrt de fissure.
L'objectif du prsent mmoire est ainsi de proposer un modle de propagation permettant de
prdire le comportement d'une fissure de clivage, en terme de vitesse de propagation et d'arrt
de fissure.
Le chapitre bibliographique expose les concepts inhrents la mcanique de la rupture en
statique et en dynamique. Il dcrit galement les principaux essais de propagation dynamique et
d'arrt de fissure fragile, et donne un aperu des principales mthodes numriques utilises pour
modliser ces phnomnes.
Le second chapitre entre dans le dtail de la mthode des lments finis tendus (ou eXtended
Finite Element Method). C'est en effet, cette mthode moderne que nous avons choisi de
dvelopper dans le logiciel de calcul Cast3m, au vu des avantages qu'elle prsente pour la
simulation numrique de propagation arbitraire de fissure. Deux techniques pour amliorer la
mthode sont prsentes et prouves dans ce chapitre.
Le troisime chapitre regroupe les donnes exprimentales de l'tude. Une brve prsentation
du matriau y est propose. Les conditions exprimentales et les rsultats d'essais sont dcrits
pour chacune des configurations tudies. Une analyse fractographique des facis de rupture
identifie les mcanismes de ruine impliqus pendant la propagation de fissure.
Finalement, le quatrime chapitre dcrit les modlisations effectues pour l'analyse de la
propagation et de l'arrt de fissure fragile. En particulier, un modle de propagation bas sur la
contrainte principale maximale en pointe de fissure est propos et discut.
20
Chapitre I.
tude bibliographique
Sommaire
I.1. Introduction..................................................................................................
....................22
I.2. Mcanique de la rupture en statique..........................................................................
....22
I.2.1. Rappel de mcanique de la rupture en statique.........................................................22
I.2.2. Modles d'amorage de fissure en clivage...............................................................
..25
I.2.3. Remarque sur les analyses statiques pour traiter de l'arrt de fissure.......................28
I.3. Propagation dynamique de fissure.........................................................................
........30
I.3.1. Champs asymptotiques en dynamique..............................................................
.........30
I.3.2. Bilan nergtique.................................................................................................
.......33
I.3.3. Solutions analytiques de problmes de propagation de fissure.................................38
I.3.4. Analyse lasto-dynamique d'essai de rupture fragile................................................40
I.3.5. Analyse lasto-viscoplastique d'essais de rupture fragile.........................................42
I.4. Essais de propagation dynamique et d'arrt de fissure...............................................46
I.4.1. Essais dfinissant une temprature d'arrt.................................................................46
I.4.2. prouvettes drives de la CT (Compact Tension)...................................................47
I.4.3. prouvettes charges dynamiquement..............................................................
.........48
I.4.4. prouvettes annulaires et choc thermique.................................................................50
I.5. Mthodes numriques pour simuler la propagation de fissure..................................52
I.5.1. Mthodes "classiques" bases sur les lments finis.................................................52
I.5.2. Mthodes non bases sur les lments finis..............................................................55
I.5.3. Mthodes bases sur la partition de l'unit................................................................56
I.6. Conclusion de l'tude bibliographique..........................................................................61
Chapitre I. tude bibliographique
I.1. Introduction
Ce chapitre bibliographique est consacr l'analyse de la propagation dynamique de fissure
fragile. L'tat des lieux qui y est dress n'est pas exhaustif, mais permet de situer les
dveloppements raliss au cours de cette thse. Le lecteur intress par certains sujets est
invit se reporter des ouvrages plus complets. Le livre de Lemaitre et Chaboche [6] donne
quelques lments de comprhension concernant la fissuration dans le contexte gnral de la
mcanique des milieux continus. Les ouvrages de Bui [7], de Kanninen et Popelar [8], et
d'Anderson [9] approfondissent les notions de la mcanique de la rupture. Enfin, l'ouvrage de
Freund [10] constitue une rfrence incontournable dans le domaine de la mcanique de la
rupture en dynamique.
Le chapitre bibliographique s'organise de la manire suivante. Les notions essentielles de
mcanique de la rupture en statique sont exposes dans un premier paragraphe. Le second
paragraphe rcapitule les principaux modles de propagation dynamique de fissure rencontrs
dans la littrature. Le troisime paragraphe prsente des essais de rupture caractristiques de ce
domaine. Enfin, le quatrime paragraphe donne un aperu des mthodes numriques
frquemment utilises pour modliser la propagation de fissure.
KM
ij r ,
M
ij avec M ={ I , II , III } (2)
r0
2 r
o ij dsigne le tenseur des contraintes, ijM sont des fonctions angulaires de . Les trois
modes de rupture indics par la lettre M sont reprsents sur la Fig. 1 (a).
22
Chapitre I. tude bibliographique
(a) (b)
Fig. 1 : (a) Modes de rupture. (b) Repre local en pointe de fissure et contrainte asymptotique
dans le cas du mode I.
Les facteurs d'intensit des contraintes sont par ailleurs directement relis au taux de restitution
d'nergie par la relation (3).
{
E
K 2I K 2II K 2III *
en dformations planes
G = avec E = 1 (3)
*
E
E en contraintes planes
L'intgrale J donne par Rice [11] permet d'valuer le taux de restitution d'nergie dans le cas
d'un matriau lastique (linaire ou non). Cette intgrale donne par l'quation (4) est exprime
dans le repre local de la fissure (cf. Fig. 1 (b)).
t
u
J = w def
n 1 ij
x1
n j d o w def = ij ij dt (4)
0
L'intgrale J utilise un contour ferm qui entoure la pointe de fissure (cf. Fig. 2).
x2
x1
23
Chapitre I. tude bibliographique
Champs HRR : ij r
0
0
J
0 0 I n r n1
ij n ,
L'intensit des champs mcaniques se caractrise donc par l'intgrale J (dfinie par l'quation
(4)) dont Rice [11] tend le concept au cas non linaire en utilisant l'nergie de dformation
totale (incluant les dformations plastiques irrversibles). Dans cette formulation, on apparente
l'nergie de dformation totale un comportement lastique non linaire, ce qui limite
l'application de cette intgrale au cas d'un matriau lasto-plastique soumis un chargement
monotone proportionnel. C'est le cas rencontr lorsque l'on traite de l'amorage et c'est pourquoi
les analyses exprimentales utilisent J comme paramtre de rfrence pour dfinir l'tat de
sollicitation. Cette interprtation est illustre par la Fig. 3.
P , a
F , Jda = d
F 0 a
dF
a
ada
a
24
Chapitre I. tude bibliographique
Approches globales
L'approche globale dterministe utilise la notion de facteur d'intensit des contrainte critique.
Lorsque la plasticit est confine la pointe de fissure comme sur la Fig. 4(a), le facteur
d'intensit des contraintes constitue un paramtre permettant de caractriser lui seul l'tat de
sollicitation en pointe de fissure. Ce premier critre de rupture caractrise donc les matriaux
fragiles, et s'nonce simplement par la relation (6).
K I = K Ic (6)
Les codifications, comme l'ASTM [14] ou le RCCM [3], dfinissent pour les aciers de cuve une
courbe enveloppe minimale de la tnacit K Ic en fonction de la temprature T en pointe de
fissure (7).
K Ic = min { 220, 29.431.355exp 0.0261T RT NDT 88.9 } en MPa m (7)
25
Chapitre I. tude bibliographique
Fig. 5 : Master Curves pour l'acier 16MND5 et essais isothermes sur CT25 (tir de [17]).
Dans cette approche, les coefficients A et B dpendent de la probabilit de rupture p, mais
galement de la largeur du front de fissure. De nombreux travaux de R&D ont t mens afin
d'identifier la temprature de rfrence T 0 , notamment concernant le nombre minimum
d'essais ncessaires, ou la prise en compte de l'amorage ductile avant la rupture fragile [18].
Approches locales
Plutt que de considrer des paramtres globaux comme le facteur d'intensit des contraintes,
les approches locales considrent les champs mcaniques proximit du front de fissure.
Introduit historiquement par Ritchie, Knott et Rice [19], le modle RKR suppose que la rupture
par clivage se produit lorsque la contrainte 22 atteint une valeur critique une distance
donne de la pointe de fissure (cf. quation (9)). Cette distance est relie par ces auteurs une
caractristique microstructurale du matriau qui est gnralement prise comme tant environ
26
Chapitre I. tude bibliographique
P R = 1exp ln
[ V
1 dV
1 pR V i i
] (10)
Ce type de modle suppose que la distribution de dfauts prsents dans le matriau (caractrise
par le paramtre m) sont susceptibles de cliver avec une probabilit value depuis la relation
empirique (11).
m
ln
1
1 p R i
=
0
u (11)
[ ]
m 1/ m
P R = 1exp W
u
avec la contrainte de Weibull W
[
=
V
dV
V0
p
m
I
] (12)
Fig. 6 : Calage partir d'essais sur CT25 du modle de Beremin pour l'acier 16MND5 la
temprature de -125C (tir de [21]).
27
Chapitre I. tude bibliographique
Une application du modle de Beremin l'acier 16MND5 pour une temprature de -125C est
donne sur la Fig. 6.
Ce type de modle permet de prendre en compte l'aspect stochastique de la rupture fragile et de
traiter certains phnomnes locaux agissant sur les contraintes qu'un paramtre global (comme
le facteur d'intensit des contraintes) ne peut dceler.
Diverse amliorations du modle de Beremin ont t proposes pour rendre compte de
phnomnes particuliers.
Lorsque la plasticit devient importante, et que l'moussement ou la contraction en direction
transverse diminue la nocivit de la fissure, la correction de la contrainte principale propose
dans [20] peut tre adopte (13).
corrige
I = I exp
Ip
k
avec 2k4 (13)
[ ]
m
W Wmin
P R = 1exp (14)
u Wmin
P R t = 1exp max
[ Vp 0t
p 0
{ }
mI dV
u T V 0
] (15)
Cependant, le modle de Beremin s'avre aujourd'hui encore incomplet. Par exemple, Bordet et
al [24] et [25] montre que la nuclation des dfauts n'apparat pas ds le dbut de la
plastification, mais que l'on peut dfinir une probabilit de nuclation en fonction de la
dformation plastique et de la limite d'lasticit. L'auteur exprime ainsi la probabilit de rupture
comme le produit de la probabilit de nuclation par la probabilit de propagation (fonction de
la contrainte principale). Cette phase de nuclation des dfauts est galement traite par
Margolin et al [26].
L'application de ce type de modle ncessite de calculer trs prcisment les champs au
voisinage de la pointe de fissure. Cela suppose une connaissance approfondie du comportement
du matriau. De plus, l'identification des paramtres du modle demande de vastes campagnes
exprimentales.
I.2.3. Remarque sur les analyses statiques pour traiter de l'arrt de fissure
Lorsque l'on traite de l'arrt de fissure conscutif une rupture fragile d'une structure sollicite
en quasi-statique, on peut dans un premier temps penser raliser une analyse statique du
problme. D'un point de vue pratique, ngliger l'aspect dynamique de la propagation de fissure
28
Chapitre I. tude bibliographique
permet de simplifier les modles et de considrer une dfinition statique du facteur d'intensit
des contraintes l'arrt (not ici K a par opposition avec sa dfinition dynamique K A ).
(a)
K A dynamic analysis
K a static analysis
KA
Ka (b)
29
Chapitre I. tude bibliographique
est particulirement intressante car elle dmontre la ncessit de prendre en compte l'aspect
dynamique de la propagation pour prdire l'arrt. En effet les interprtations statiques de
diverses prouvettes aboutissent des valeurs du facteur d'intensit des contraintes l'arrt trs
variables contrairement l'analyse dynamique. La Fig. 7 (b) explique pourquoi les valeurs de
K a et de K A diffrent. Le mouvement de la structure ne se fige pas instantanment aprs
l'arrt de la fissure. Pendant un certain temps, le facteur d'intensit des contraintes oscille autour
de K a des valeurs infrieures K A .
Au vu de ces observations, il semble ncessaire d'aborder en profondeur l'aspect dynamique de
la propagation de fissure de clivage.
crack plane x2
x1
22 = sign x2 H t
{
where H t = 0 if t 0
1 if t 0
30
Chapitre I. tude bibliographique
dans les paragraphes suivants. Nous poserons d'abord les hypothses gnrales de l'analyse,
avant de donner l'expression des champs asymptotiques, solutions du problme de la
propagation dynamique d'une fissure dans un milieu lastique.
La masse volumique est dsigne par , tandis que et sont les coefficients de Lam,
relis au module d'Young et au coefficient de Poisson par les relations :
E E
= , = (19)
112 21
Les termes T , K , , et C reprsentent respectivement la variation de temprature, la
conductivit thermique, le coefficient de dilatation thermique, et la chaleur spcifique
dformation constante.
T =T T 0 , K =3 2 (20)
La contrainte et l'entropie drivant de ce potentiel, on peut crire leur loi d'volution (21).
{
wdef
{
ij = = =u k ,k ij 2 u i , j K T ij = tr 12
i , j u i , j
soit C (21)
C s= K tr p T
s= =Ku k , k T T0
T T0
31
Chapitre I. tude bibliographique
t
2
(24)
c rot u = 2 rotu
2
s
t
o cd =
2
et c s=
reprsentent respectivement la vitesse des ondes de dilatation et
de cisaillement en dformations planes. Nous nous plaons dsormais dans un contexte
isotherme. Freund [10] donne dans son ouvrage la forme des champs asymptotiques en lasto-
dynamique isotherme. Le champ de dplacement prend la forme donne par l'quation (25).
K dyn
M t
u i = ui r 0 r 1/ 2 U iM , a ui r 1 ... (25)
2
La partie correspondant aux contraintes asymptotiques utilisent les fonctions suivantes :
{ { }
I 1 2 d s
Ux= 21 s d cos 4 d s s cos
D 2 2
pour le mode M=I,
U Iy =
1
D { }
2 d 12s d sin d 4 d s sin s
2 2
{ = {4 sin 2 1 sin }
II 1 d 2 s
Ux s d s s s
D 2 2
pour le mode M=II,
= {4 cos 21 cos }
II 1 d 2 s
U y d s d s s
D 2 2
et U z =
III 1 2
D s{ s cos
s
2 } pour le mode M=III.
2 2
avec : i = 1
a
ci , i = 1
a sin
ci , et tan i = i tan pour i=d , s .
32
Chapitre I. tude bibliographique
{
{ }
2 2 2
I 1 1 s 12 d s d 4 s d s
xx = cos cos
D d 2 s 2
2
2 d 1s
I
xy =
D { 1 d
sin
1
sin
d 2 s 2
s
} pour le mode M=I,
{ }
2 2
I 1 1 s 4 s d
yy = cos d sin s
D d 2 s 2
{
{ }
2 2 2
II 2 s 12 d s d 1 s s
xx = sin sin
D d 2 s 2
II
xy =
{
1 4 s d
D d
d 12s 2
cos
2 s
cos s
2 } pour le mode M=II,
2 s 12s
II
yy =
D { 1
sin d
1
d 2 s 2
sin s
}
III 1 s III 1 s
et xz = sin , yz = cos pour le mode M=III.
s s 2 s 2
I.3.2. Bilan nergtique
Dfinition du taux de restitution d'nergie en lasto-dynamique
Ajoutons au premier principe de la thermodynamique un terme dissipatif correspondant la
dcohsion de la matire pour obtenir le bilan (27).
P int P acc = P ext Q P diss
{
d
{
ext ext
int
P = def
w Q P = F u
dt (27)
avec et Q = rdiv q
acc d cin
P = w P diss = G a
dt
La puissance intrieure et cintique est gale la somme de la puissance des efforts extrieurs,
du taux de chaleur et de la puissance dissipe. La densit d'nergie de dformation w def a t
cin 1
dfinie par l'quation (4). La densit d'nergie cintique vaut w = 2 u i ui . r et q
reprsentent respectivement la densit de production interne de chaleur et le vecteur courant de
chaleur.
On introduit le domaine =0 reprsent sur la Fig. 10.
33
Chapitre I. tude bibliographique
+ V
-
ext
De plus, on a le dveloppement :
d
Grce au thorme du transport : ...= t ... ...Vn , on peut r-crire les
dt
premiers termes de droite de l'quation (28):
d
dt
w def w cin d =
(31)
w def w cin d w def w cin V n d dtd w def w cin d
0
Finalement, on obtient :
34
Chapitre I. tude bibliographique
d
P diss = w def wcin V u nd dt w def wcin d
(32)
F
On retrouve l'intgrale flux d'nergie F introduite par Freund dans [28]. Celle-ci traduit la
somme de l'nergie mcanique qui traverse et de l'nergie rsultant des actions extrieures
agissant sur ce contour.
Si la puissance dissipe est indpendante du contour d'intgration, cela n'est pas le cas de
F . En effet, on dduit directement de l'quation (31) :
d d
F 2 F 1 =
dt w def wcin d
2
dt
w def w cin d
1
(33)
=
C
w def w cin V w def w cin d
t k
xk
Dans le cas d'un matriau lastique linaire, on trouve finalement :
Freund dfinit le taux de restitution d'nergie dynamique comme tant le flux d'nergie
absorbe par la pointe de fissure par unit d'avance de la fissure. Si on identifie la vitesse du
contour V celle de la fissure a , on obtient :
F 1
G = lim
0 a
= lim
0 a
w def w cin a j ij ui n j d (35)
+
+ 1
- 1
35
Chapitre I. tude bibliographique
aG = w def w cin a j ij ui n j d
2
(36)
lim ij ui , j a k u i , kj ui ui ak ui ,k d
1 0
C
On constate qu'il subsiste dans l'expression (36) une intgrale de surface dont l'une des
frontires ( 1 ) est vanouissante. Supposons que le premier contour 1 englobe une rgion
o l'hypothse de stationnarit (dcrite par l'quation (37) pour le champ f i ) s'applique pour
les champs de dplacement u i et de vitesse u i .
d fi fi
0 soit f i ak (37)
dt xk
Alors, l'quation de G se simplifie en une expression indpendante du contour :
G w w 1j ij ui ,1 n j d
def cin
(38)
Si on injecte dans l'expression prcdente les champs asymptotiques dcrits par les quations
(25) et (26), on obtient la relation (39).
12 1
G=
E
[ AI a K 2I A II a K 2II ] A III a K 2III
(39)
36
Chapitre I. tude bibliographique
J = lim w w kj ij u i , k a k / a n j d
def cin
1 0
(41)
1
J = lim w w 1j ij u i ,1 n j d
def cin
1 0
(42)
1
J = lim
1 0
H 1j n j d = lim H kj qk n j d
1 0
1 1
donc J = lim
1 0
2 + -
H kj q k n j d H kj , j q k H kj qk , j d
C (43)
w def
avec H kj , j = u u ui ,k ij , j ui ,k ij u i , kj = u u i , k u i ui ,k
u i , j i , jk
Si on rutilise les hypothses prcdemment faites de lvres libres d'effort et suffisamment
droites, on obtient finalement l'expression de l'intgrale J en lasto-dynamique :
J
dyn
= lim
1 0
[ w def cin
]
w kj ij ui ,k q k , j u i ui ,k ui ui ,k q k d (44)
C
La mise en oeuvre numrique de cette intgrale n'est pas aussi directe qu'elle y parat. La
prsence d'une limite demande une certaine prudence. Afin que les hypothses de stationnarit
locale et par consquent l'indpendance du contour s'appliquent, il est ncessaire de conserver
cette limite.
37
Chapitre I. tude bibliographique
ij u i , k , j ij u i ,kj ui u i ,k ui u i , k =0 (45)
t
Ensuite, l'identification des densits d'nergie de dformation rversible et d'nergie cintique
conduisent une loi de conservation proche de celles donnes par Eshelby [40] (en statique) et
par Fletcher [41] (en dynamique). Il y a seulement le terme concernant la thermique en plus.
On remarque qu'on peut arranger les diffrents termes de l'expression prcdente pour obtenir :
Si on intgre l'quation (46) (en prenant k gal 1) sur un domaine C qui suit la propagation
rectiligne de la fissure la vitesse V = a x1 , et si on estime que le contour 1 tends
suffisamment vers 0 pour que la couronne se confondent avec le domaine total 2=C 1 ,
alors on obtient l'intgrale indpendante du contour donne par Bui [7] :
d
J dyn
Bui = w def w cin ui ui ,1 a j1 ij ui ,1 n j d d t ui ui ,1 d (48)
2
Pour l'implantation numrique dans Cast3M, le dernier terme est valu en faisant la diffrence
des quantits aux instants t et tt .
Si l'on multiplie directement l'quation de conservation (46) par le vecteur q, et que l'on intgre
sur un domaine contenant la pointe de fissure, alors on peut en dduire de manire quasi-
identique l'expression de J donne par Rthor [34] :
dyn
J Rthor = w def w cin jk ij ui ,k qk , j d ui ui ,k ui ui , k qk d (50)
Enfin, si l'on multiplie l'quation de conservation dans sa forme modifie (45) par le vecteur q,
alors on trouve l'expression (51), dj donne par l'quation (44).
J
dyn
= w def wcin jk ij ui , k qk , j d ui ui , k ui ui , k qk d (51)
C C
Ces diffrentes expressions de l'intgrale Jdyn coexistent cause des hypothses de stationnarit
locale et des passages la limite des domaines d'intgration. Leur valuation numrique donne
en gnral des valeurs trs proches.
38
Chapitre I. tude bibliographique
restreint le nombre de solutions analytiques disponibles. On dispose des solutions tablies par
Mott [42], Broberg [43], Rose [44] ou Freund [10] qui correspondent des cas relativement
simples. Par exemple, le cas acadmique de Broberg, gnralis par la suite par Freund,
considre une propagation vitesse constante dans un milieu lastique infini soumis un
chargement constant de traction sur les lvres (cf. Fig. 13).
x2 x2
22 x1 at = 22 x1 at = px 1
u 2 x1 at = 0 x1 u 2 x1 at = 0 x1
a a
(a) (b)
Fig. 13 : (a) Cas de Broberg. (b) Sa gnralisation par Freund.
Broberg [43] tablit la relation (52).
K dyn t , a I a / cs R1/ a 2 2
stat = a cs (52)
K a 1 a/ cd 2
o K stat a= a , R et I sont des fonctions de la vitesse de propagation et des vitesses
caractristiques du matriau. Freund [10] donne la relation (53) qui est assez proche de la
prcdente.
K dyn a , a = k a K stat a ,0 (53)
o K stat a ,0 est le facteur d'intensit des contraintes obtenus pour le mme chargement et la
longueur de fissure l'instant considr en quasi-statique (54).
a t
K stat a ,0=
2
px
0 at x
dx (54)
39
Chapitre I. tude bibliographique
40
Chapitre I. tude bibliographique
K uI r=
0
2
1 r
[ u y uy ]
On remarque que cette dfinition vite l'introduction des fonctions universelles de la Fig. 12
(58)
F K I ,KI = 0
u
(60)
Maigre [47] dcrit la complexit qu'il existe exprimer un tel critre pour du verre qui est
pourtant un matriau trs fragile.
Il existe d'autre forme pour exprimer la dpendance de la tnacit avec la vitesse de propagation
comme par exemple celle employe dans [50] ou [51]. Cette dernire s'crit :
K dyn
Ic = K A T R T NDT A T R T NDT a
2
(62)
41
Chapitre I. tude bibliographique
(a) (b)
Fig. 15 : Corrlation entre vitesse de propagation et facteur d'intensit des contraintes Kdyn.
(a) tir de [49], et (b) tir de [51].
Dans les cas o la plasticit peut tre nglige, ce type d'approche est la fois simple et
efficace. Il ncessite une identification pralable de plusieurs paramtres qui n'est pas toujours
ais (cf. Fig. 15(b)). Par ailleurs, il est parfois ncessaire de diffrencier la tnacit l'arrt
K A de la tnacit l'amorage K I c pour certains matriaux.
42
Chapitre I. tude bibliographique
calage repose sur des essais dynamiques raliss la plupart du temps l'aide de barres de
Hopkinson (Split Hopkinson Pressure Bar) qui permettent d'obtenir des vitesses de
dformations de l'ordre de 102s-1 103s-1 (cf. [52]).
Rossol [53] ralise une identification viscoplastique d'un acier de type 16MND5 avec une loi de
type Cowper-Symonds [54] (cf . quation (63)) pour les trois tempratures (-100C, -40C,
20C) et les trois vitesses de dformations (10-3s-1, 0.56s-1, 5.52s-1).
[ ]
1
in stat in
, ,T = ,T 1 p
(63)
D
Il trouve le couple de paramtres : p=12 et D=108s-1 .
Afin de modliser l'essai Charpy, Tanguy [55] identifie galement le comportement de l'acier
16MND5 pour des vitesses de dformations s'talant de 5.10-4s-1 4.103s-1. L'auteur propose
l'utilisation d'une double loi de Norton correspondant aux deux mcanismes de plastification
dcrits par la Fig. 16. Il s'agit donc d'une loi lasto-viscoplastique additive avec quatre
paramtres dpendant de la temprature. La comparaison du modle propos avec les rsultats
exprimentaux issus d'essais sur barres de Hopkinson est prsente sur la Fig. 17.
43
Chapitre I. tude bibliographique
Modle de Freund
Une difficult majeure lorsque l'on traite de mcanique de la rupture en non-linaire, demeure
la dfinition d'un paramtre pertinent mesurant l'intensit des champs en pointe de fissure, car
avec la propagation se crent des zones locales de dcharge. L'expression analytique d'une
solution asymptotique dans un matriau dpendant de l'histoire n'est pas une chose simple et le
cas d'une fissure se propageant dans un matriau lasto-viscoplastique une vitesse pouvant
varier ne connat pas de solution analytique.
C'est pourquoi Freund et Hutchinson [57], puis Xu et al [58] cherchent exprimer l'intensit
des champs en pointe de fissure depuis l'intensit du chargement extrieur, moyennant quelques
hypothses. Ils supposent d'abord l'existence d'une zone trs proche de la pointe de fissure
l'intrieur de laquelle les vitesses de dformation sont si grandes que la limite d'lasticit est
trs leve et que les champs asymptotiques dominants (quantifis par Gtip) sont de la forme des
champs lasto-dynamiques (cf. Fig. 19).
plastic wake :
in0 act
-h
Fig. 19 : Zone active et sillage plastique laisse par la propagation d'une fissure dans un
matriau lasto-viscoplastique.
Ensuite, un bilan nergtique est ralis (cf. quation (64)).
44
Chapitre I. tude bibliographique
h
1
Gtip = Gext
a
ij ijp d U el * dh
act
(64)
h
Ce bilan nonce simplement que l'nergie consomme par la pointe de fissure (Gtip) est gale
l'nergie extrieure (Gext) laquelle on te l'nergie de dformation plastique et l'nergie
lastique (Uel*) bloque dans le sillage plastique. Le dernier terme est jug par Freund comme
ngligeable devant les autres, ce qui a t numriquement vrifi par lui-mme et par Xu pour
les cas qu'ils ont tudi.
Le terme correspondant la dissipation plastique est valu en supposant les contraintes
singulires en r-1/2 dans toute la zone plastique et les dformations inlastiques sont dduites du
comportement viscoplastique du matriau. Selon le choix de l'intensit de la singularit,
l'approximation varie quelque peu. Freund et al [57] choisissent la premire expression de
l'quation (65), surestimant ainsi les contraintes et sous-estimant la dissipation plastique.
Mataga et al [59] utilisent la seconde expression de l'quation (65).
K tip
Approximation de Freund et al : ij ij , a
2 r
(65)
Approximation de Mataga et al : ij
K tip K ext
2r
ij , a
Xu montre numriquement que cette dernire surestime la dissipation plastique. Quelque soit
l'expression utilise, l'introduction d'un critre du type Gtip =Gc amne l'expression du
rapport du taux de restitution d'nergie en pointe par le taux de restitution d'nergie extrieur.
Par exemple, la seconde approximation conduit l'quation (66) illustre par la Fig. 20.
G ext 1
= (66)
G tip 1 D a / cR P c
45
Chapitre I. tude bibliographique
jusqu' atteindre la partie ascendante de la courbe. Uniquement cette partie est stable par rapport
la sollicitation. En effet, si l'on perturbe la vitesse de propagation de telle sorte qu'elle
augmente un peu, alors la rponse mcanique du systme est la diminution de G (cf. paragraphe
I.3.3). Puisque l'on considre la branche ascendante de la Fig. 20, l'application du modle
produit alors une diminution de la vitesse de propagation, rgulant ainsi cette dernire. Comme
la perturbation est attnue par le modle, on en conclut que seule la partie ascendante de la
courbe est mcaniquement stable. Finalement l'arrt de fissure doit ncessairement se produire
pour le minimum de la courbe (point M de la Fig. 20).
Malheureusement, les vitesses de propagation prdites par ce modle sont bien suprieures
celles observes exprimentalement (ces dernires n'excdent gnralement pas 0.4 cR ).
Critres locaux
Plutt que de dfinir un paramtre global, une alternative consiste utiliser un critre local en
pointe de fissuree reprsentatif des mcanismes de ruine impliqus dans la rupture fragile. C'est
le cas du modle RKR [19] que nous avons prsent au paragraphe I.2.2. Ce critre a t utilis
dans le contexte d'une analyse lasto-viscoplastique de choc thermique sur anneau par Hajjaj
[56]. L'auteur propose une dpendance de la contrainte critique avec la temprature similaire
celle de la tnacit. Cette dmarche lui permet de retrouver des valeurs de longueur de fissure
l'arrt trs proches des rsultats exprimentaux, mme si la vitesse de propagation est moins
bien reproduite. Puisque ce type de critre utilise une valeur ponctuelle de la contrainte, sa mise
en oeuvre numrique introduit des imprcisions.
Il existe d'autres manires pour mettre en oeuvre un critre bas sur la contrainte principale
maximale afin de se soustraire aux possibles erreurs numriques. Par exemple, Wells et al [60]
et Remmers et al [61] s'inspirent des mthodes de dlocalisation utilises en mcanique de
l'endommagement, pour moyenner le tenseur des contraintes en pointe de fissure. Cet aspect
sera abord plus en dtails dans la partie consacre la modlisation des essais raliss.
Les concepts probabilistes (comme ceux utiliss par le modle de Beremin) sont difficilement
applicables la propagation de fissure, car la longueur de fissure un instant donne est la
fois la donne et le rsultat du calcul. Dans ce contexte, il semble difficile de transposer
l'hypothse de rupture de l'ensemble sujette la rupture du maillon le plus faible. De plus, les
mcanismes en jeu sont moins bien prcisment connus. C'est pourquoi, la connaissance de
l'auteur, l'analyse de la propagation de fissure est toujours dterministe.
Essai Pellini
Un cordon de soudure est dpos sur une plaque rectangulaire, puis est entaill (cf. Fig. 21(b)).
Un projectile vient impacter la plaque du cot oppos la soudure la mettant en flexion, et
ouvrant par consquent l'entaille. On dfinit la Nil Ductility Transition Temperature (NDTT)
46
Chapitre I. tude bibliographique
comme tant la temprature maximale pour laquelle le cordon est compltement travers par la
fissure amorce dans l'entaille. Ce type d'essai est recommand dans [63] pour caractriser
l'arrt de fissure.
cordon de soudure
impact entaille
entaille
impact
(a) (b)
Fig. 21 : Description des essais Robertson (a) et Pellini (b).
47
Chapitre I. tude bibliographique
a t
wx,t
projectile
SENB
entaille
48
Chapitre I. tude bibliographique
d'un projectile qui met en flexion l'prouvette (cf. Fig. 24). L'aspect dynamique de la
sollicitation ncessite une prise en compte trs fine de la dimension temporelle du problme, car
le facteur d'intensit des contraintes varie beaucoup pendant et aprs l'impact.
Expriences de Kalthoff
L'influence de la dynamique du chargement est particulirement approfondie par Kalthoff [27],
qui met en vidence l'influence de la vitesse de monte en charge dans l'une de ces expriences.
Cette dernire est prsente sur la Fig. 25. L'prouvette prsente deux entailles symtriques et
est impacte par un projectile dont on peut faire varier la vitesse. Dans [72] et [73], Kalthoff
montre une transition entre fissuration par ouverture et fissuration par bande de cisaillement. Le
mme phnomne est prsent pour d'autres matriaux comme le marbre [74].
projectile
prouvette
avec une double entaille
(a) (b)
Fig. 25 : Exprience de Kalthoff.
(a) Schma de principe. (b) Demi-prouvettes mettant en avant la fissuration par ouverture (
gauche) et par bande de cisaillement ( droite) sur de l'acier 42CrMo4 (tir de [73]).
La fissuration par ouverture correspond au cas prdit par les outils classiques de la mcanique
de la rupture (critre en contrainte d'ouverture maximale par exemple). Lorsque la vitesse
d'impact dpasse une certaine limite, des bandes de cisaillement apparaissent. Ce processus de
rupture consomme une quantit d'nergie suprieure la fissuration par ouverture et l'aspect de
la surface de rupture est dans ce cas beaucoup moins lisse.
49
Chapitre I. tude bibliographique
entaille
CCS
impact
Premires interprtations
Bouyne [65] a men une tude sur l'acier 21/4Cr-1 Mo (en comparaison avec l'acier 16MND5)
en utilisant les prouvettes CCA, annulaires isothermes et annulaires sous choc thermique. Il
compare l'volution avec la temprature de la tnacit l'amorage et l'arrt, selon le modle
utilis (modles de Ritchie-Knott-Rice [19], et de Beremin [20]). Ainsi, concernant l'amorage,
l'auteur met en vidence que les deux modles sous-estiment les valeurs exprimentales de
contraintes critiques au del de 100C.
L'auteur utilise ces modles d'amorage pour la propagation et l'arrt de fissure. Seuls les
anneaux isothermes sont tudis en dynamique, mme si l'auteur recommande cette voie
d'tude. Il dmontre en effet que le ralentissement la fissure s'accompagne de l'lvation du
facteur d'intensit des contraintes. Il estime galement la vitesse de dformation aux environs de
50
Chapitre I. tude bibliographique
103 104 s-1 et conseille ainsi sa prise en compte dans la loi de comportement pour des
modlisations futures. Pour ces analyses, l'auteur se contente de prendre en compte la viscosit
du matriau par le biais de la variation de la limite d'lasticit. D'autre part, la plasticit n'est
intgre que pour l'analyse des anneaux sous choc thermique.
Les modles d'amorage appliqus l'arrt conduisent des valeurs de contraintes critiques
survalues. L'auteur conclut que l'nergie consomme durant la propagation se divise en trois
termes : la premire partie sert au clivage, la seconde partie est absorbe par les dformations
inlastiques, et la dernire reprsente l'nergie ncessaire pour dchirer les ligaments non
rompus par le clivage et observs l'arrire du front.
Les essais de choc thermique ont galement t tudis au LISN (Cea Saclay). Les prouvettes
annulaires utilises permettent de combiner les sollicitations d'origine mcanique et thermique.
L'anneau entaill et pr-fissur sur le rayon extrieur, est initialement refroidi et mis en
compression. Une douzaine de trous travers lesquels on injecte de l'eau chauffe 90C
permet d'obtenir le choc thermique. Cette technologie permet de solliciter mcaniquement
l'prouvette et d'obtenir un gradient de temprature d'environ 100C, ce qui est nettement moins
abrupt qu'avec l'utilisation d'un inducteur (gradient thermique compris entre 300C et 600C).
Pendant le temps de la conduction thermique (de l'ordre de la seconde), la temprature
augmente et le gradient thermique dans le ligament diminue alors que l'effort mcanique est
maintenu constant. Cela se traduit par un trajet original dans le plan facteur d'intensit des
contraintes (K) - temprature (T). Finalement, l'interprtation de l'amorage permet d'prouver
les modles analyss [21].
Fig. 27 : Schma de l'prouvette annulaire pour choc thermique utilise au LISN [21].
51
Chapitre I. tude bibliographique
(a) (b)
Fig. 28 : Rsultats obtenus par Hajjaj [56].
(a) Contrainte critique en fonction de la temprature pour deux essais.
(b) Longueur de fissure en fonction du temps pour un calcul 3D.
Ces rsultats sont trs encourageants car ils montrent que l'utilisation d'un critre en contraintes
relativement simple permet d'obtenir une rponse satisfaisante, particulirement en ce qui
concerne la longueur de fissure l'arrt prdite.
52
Chapitre I. tude bibliographique
Enfin, il existe d'autres variantes du schma prcdent. Par exemple, Malluck et al [82] propose
un schma o on impose la force de dcrotre de manire maintenir le facteur d'intensit des
contraintes K constant durant le pas de temps.
Cette mthode ncessite la connaissance pralable du trajet de fissure et impose la pointe de
fissure de se situer sur un noeud. Elle permet de simuler des ruptures en 2D, mais est beaucoup
plus difficile mettre en oeuvre en 3D lorsque le front de fissure n'est pas rectiligne. Elle a
notamment t utilise par Kobayashi [48], Bouyne [65], et Yuritzinn [83] pour l'analyse
d'essais de rupture fragile. En dynamique, le relchement de noeuds provoque des oscillations
peu ralistes en pointe de fissure.
Mthodes de remaillage
Pour avoir une bonne ide des champs proximit du front de fissure, il faut conserver un
maillage suffisamment fin lors de la propagation. C'est pourquoi, partant de cette ide, il semble
ncessaire de remailler la fissure lors de sa propagation. Les algorithmes de remaillage sont
souvent complexes et dlicats mettre en oeuvre. On peut extraire deux principales sources
d'erreur.
D'abord, aprs une opration de remaillage, il est ncessaire de projeter les champs de la
solution du pas prcdent sur le nouveau maillage. Comme les solutions discrtises sur
diffrents maillages ne sont pas les mmes, une simple projection ne conduit pas des champs
en quilibre. Une tape supplmentaire de r-quilibrage est donc souvent introduite en prenant
garde conserver ferme la surface de fissure frachement introduite par le remaillage.
Malheureusement, cette tape de r-quilibrage ne garantit pas la conservation de l'nergie
mcanique du systme pendant la projection.
Ensuite, tout comme pour la mthode de relchement de noeud, il est ncessaire de librer la
nouvelle surface de fissure.
Finalement ces oprations de projection et de relchement modifient l'nergie du systme
discret de manire purement numrique. L'influence de ces mouvements d'nergie sur
l'application d'un critre de propagation est tudi en lastodynamique par Rthor dans [84].
A mi-chemin entre le relchement de noeud et le remaillage, le "Moving element method"
prsent par Nishioka [85] consiste repositionner les noeuds des lments proches de la pointe
de fissure pour tenir compte de sa nouvelle position.
L'utilisation d'lments finis espace-temps peut aussi tre considre comme une variante des
mthodes de remaillage puisque la discrtisation de la gomtrie varie au cours du temps. Le
principe de ces lments est de considrer le temps de la mme manire que les variables
d'espace, et de remplacer les traditionnels schmas temporels bas sur les diffrences finies par
des approximations de type lments finis, ou mme par des formulations de type Galerkin
discontinu comme dans [86] ou [87]. Comme on le voit sur la Fig. 30, le temps est la troisime
dimension d'un problme plan. L'avantage principal de ce type de formulation rside dans la
capacit raffiner le maillage o et quand on juge ncessaire de le faire grce un indicateur de
qualit de la solution. Comme cela est illustr par la Fig. 30, le front d'une onde lastique et la
pointe de fissure requirent une finesse de maillage localement plus leve que dans le reste de
la structure.
53
Chapitre I. tude bibliographique
(a) (b)
Fig. 30 : Mthode des lments finis espace-temps.
(a) Exemple de calcul d'une plaque fissure fixe soumise un chargement brutal de traction.
(b) Le raffinement suit la propagation des ondes lastiques (tir de [87]).
54
Chapitre I. tude bibliographique
[u]
G
[u]
(a) (b)
Fig. 31 : Modle de zones cohsives.
(a) Maillage. (b) Loi de comportement.
Intellectuellement trs intressante, cette mthode possde nanmoins quelques inconvnients.
D'abord, elle introduit un lment d'paisseur nul, mais de rigidit fini, qui modifie la rigidit
globale de la structure, ce qui, en dynamique, perturbe le trajet des ondes mcaniques. De plus,
il faut connatre le trajet de la fissure priori. L'argument le plus souvent reproch ce modle
est le fait qu'il introduit une dpendance au maillage. Ceci a t dmontr par Xu et Needleman
[89] en introduisant des lments cohsifs entre tous les lments, et plus rcemment par Zhou
et Molinari [90]. La direction de propagation montre alors une forte dpendance la
construction du maillage en terme d'orientation et de taille des lments.
Une premire volution de cette mthode est celle que propose Camacho, Pandolfi et Ortiz [91],
[92] et Zhou et Molinari [90]. Les auteurs n'introduisent de zones cohsives qu'aux bords des
lments qui vrifient un critre de rupture ou qui sont situs dans le prolongement immdiat de
la fissure prexistante.
Une seconde volution des zones cohsives est leur introduction via la mthode de partition de
l'unit. Nous la dveloppons donc dans le paragraphe I.5.3 de ce chapitre.
55
Chapitre I. tude bibliographique
non nulles sur une zone d'influence qui peut contenir d'autres noeuds comme cela est reprsent
sur la Fig. 32(a).
(a) (b)
Fig. 32 : Mthodes sans maillage (tir de [35]).
(a) Reprsentation des noeuds et de leur zone d'influence.
(b) Modification de la fonction de forme en prsence d'une fissure.
De par leur nature, ces mthodes sont trs efficaces lorsque l'on traite de trs grandes
dformations ou de fragmentation. Belytschko [94] et Duflot [95] donne un aperu des
possibilits de ces mthodes appliques la mcanique de la rupture. La prsence d'une fissure
peut tre modlise par l'ajout de critre(s). Le critre de "visibilit" modifie les fonctions poids
afin que deux noeuds de part et d'autre de la fissure ne s'influencent pas l'un l'autre (cf. Fig.
32(b)). Organ [35] et Krysl [96] ont dmontr la possibilit de traiter des problme de
fissuration dynamique en utilisant ce type de formulation pour la zone de propagation de la
fissure d'une part, et une mthode d'lments finis standard pour le reste de la structure d'autre
part. Les problmes de fragmentation se prtent galement trs bien aux mthodes particulaires.
Rabczuck [97] traite par exemple le cas d'une plaque en bton impacte par un projectile.
Parmi les mthodes sans maillage, la dynamique molculaire a rcemment fourni des rsultats
intressants d'un point de vue thorique sur la mcanique de la rupture dynamique. Son principe
est simple puisqu'il consiste rsoudre l'quilibre dynamique d'un grand nombre d'atomes
(plusieurs dizaines de millions) qui interagissent via un potentiel d'attraction. De la forme de
celui-ci dcoule le comportement du milieu (par exemple, on obtient un comportement
lastique linaire en utilisant un potentiel quadratique). En utilisant une loi de comportement
hyperlastique, Buehler [98] montre l'existence d'une vitesse seuil observe exprimentalement
au-del de laquelle le trajet de la fissure devient beaucoup plus chahut, pouvant mme devenir
branch. L'auteur s'est galement pench sur le problme de fissuration supersonique.
56
Chapitre I. tude bibliographique
N i x = 1 x (67)
i
Si {i , p } constitue une base de fonctions dont on espre tre une bonne approximation locale,
on peut construire et ajouter l'approximation l'ensemble de fonctions {N i i , p } . Il faut
toutefois tre prudent car rien ne garantit que cet ensemble de fonctions soit une base, et
l'existence de deux fonctions trop proches peut conduire un mauvais conditionnement de la
matrice rsoudre. Si cet ensemble constitue une base de fonctions indpendantes, alors on peut
enrichir l'approximation de la manire suivante :
u x = N i xui N i xi , p x ai , p (68)
i i, p
L'quation (68) permet de se rendre compte que la fonction x peut tre exactement
reprsente grce la prsence de la base de fonctions enrichis.
En pratique et pour des raisons de cot de calculs, on limite l'introduction de fonctions de la
base d'enrichissement aux sous-domaines i , i I o cela est ncessaire. I reprsente
par exemple l'ensemble des noeuds des lments appartenant la zone que l'on souhaite
enrichir. Ainsi on cre une couche d'lments de transition pour lesquels une partie seulement
des noeuds est enrichie et qui ne vrifie donc plus vraiment la proprit de partition de l'unit.
Chessa [100] dmontre que le traitement de cette zone de transition a des consquences sur
l'ordre de convergence de la mthode. L'auteur propose de masquer l'influence des fonctions
d'enrichissement en liminant les termes indsirables introduit par les fonctions
d'enrichissement dans cette couche d'lments de transition. Farhat [101] propose galement un
moyen de s'affranchir de cette zone de transition. L'auteur spare l'espace enrichi du reste du
domaine et impose la continuit des champs aux frontires entre ces deux sous-domaines avec
des multiplicateurs de Lagrange.
{ [
1 r
u1 x r=
2 2
0
2
K I cos cos K II sin 2cos
2 ]
u2 x r=
1
0 2
r
2 [
K I sin cos K II cos 2cos
2 2 ] (69)
u3 x r=
2 r
0
K sin
2 III
2
57
Chapitre I. tude bibliographique
{
34 en dformation plane
La constante de Kolosov est donne par : = 3
en contrainte plane
1
{
sin / 2
u i x r=
0 F j xui j o F j x = r sin / 2 sin (70)
j cos / 2
cos / 2 sin
Grce la mthode de partition de l'unit, ces enrichissements sont ajouts au champ de
dplacement.
Ensuite, Mos et al [103] liminent compltement les oprations de remaillage en ajoutant une
fonction discontinue dans la base d'enrichissement (cf. quation (71)). Ainsi, le (re-)maillage ne
reprsente plus une contrainte lors de la modlisation par lments fini de propagation de
fissures.
{
H x = 1 si 0
1 si 0
(71)
crack surface:
=0
crack front:
=0
58
Chapitre I. tude bibliographique
La stratgie d'enrichissement peut tre rsume par l'quation (72) et la Fig. 34.
u x
i N i x ui H x i a i
k =1,..4
F k xi b i ,k
(72)
soit u x N i xqi o Ni = { N i , N i H , N i F k }
i
59
Chapitre I. tude bibliographique
60
Chapitre I. tude bibliographique
Le premier cas est prsent par Simone dans [116] avec une loi non-locale d'endommagement
pour la partie continue. L'auteur montre que l'ajout de fonctions "saut" permet de matriser le
comportement du matriau une fois le dommage critique atteint (la Fig. 36 montre que le
paramtre d'endommagement cesse de crotre aprs avoir atteint sa valeur critique).
Fig. 36 : Paramtre d'endommagement pour une prouvette flexion 4 points (tir de [116]).
Comparaison entre un modle d'endommagement seul ("continuous failure"),
et un modle d'endommagement accompagn d'un dplacement enrichi par une fonction
discontinue ("continuous -discontinuous failure").
Le second cas reprsente une avance importante dans l'utilisation des modles cohsifs.
L'utilisation de la partition de l'unit permet de rsoudre les problmes de trajet de fissure
dpendant du maillage et de raideurs numriques de l'interface. Cette mthode est en fait trs
proche de la mthode X-FEM originale, puisqu'au lieu d'utiliser des fonctions singulires en
pointe de fissure, on introduit des zones cohsives censes tre reprsentatives du processus de
rupture (cf. modle de Dugdale-Barenblatt). Cette mthode est trs prise (cf. [109], [111]) car
elle permet de s'affranchir des problmes d'lments de transition ("blending elements") et
d'intgration numrique des fonctions singulires.
61
Chapitre I. tude bibliographique
norme ASTM) et les prouvettes impactes (comme les SENB impactes) ne semblent pas
adquates dans ce contexte, moins de raliser une mesure prcise permettant de dterminer le
chargement et les conditions aux limites.
Il apparat galement qu'une prise en compte fine du comportement du matriau ncessite son
identification exprimentale de hautes vitesses de dformation. Les mcanismes de rupture en
jeu requirent galement une attention toute particulire, car certains auteurs attribuent la
prsence de ligaments non rompus un rle trs important dans l'arrt de fissure.
Les mthodes numriques utilises pour simuler la propagation de fissure ne sont pas toutes
quivalentes. Il semble que les mthodes classiques bases sur la mthode des lments finis
montrent certaines limites face ce type d'analyse. L'utilisation de zones cohsives est
relativement difficile mettre en oeuvre notamment afin de se prmunir de la dpendance au
maillage. Les mthodes sans maillage semblent tre dans ce contexte une alternative attrayante,
bien que leur usage reste peu rpandu, et que leur couplage avec la mthode des lments finis
traditionnels demande des dveloppements coteux. Enfin, les mthodes bases sur la partition
de l'unit permettent de raliser des simulations efficaces et bien adaptes ce type de
problme.
Cette tude bibliographique souligne le fait que l'analyse de la propagation et de l'arrt de
fissure de clivage initie sous un chargement quasi-statique ncessite la mise en oeuvre d'outils
de diverses natures.
Pour ne pas limiter les interprtations et tre capable de reproduire prcisment les essais de
rupture, il faut possder un outil numrique adapt la propagation arbitraire de fissure. C'est
pourquoi, nous avons choisi de dvelopper la mthode des lments finis tendus dans Cast3m.
Le comportement du matriau aux grandes vitesses de dformation doit galement tre
identifi. Des essais de propagation et d'arrt de fissure, ainsi que la mesure de la vitesse de
propagation, constituent les donnes indispensables l'analyse des phnomnes.
Des critres permettant de raliser des simulations prdictives de la propagation et de l'arrt de
fissure doivent tre proposs et prouvs. La comparaison des rsultats fournis par les
simulations prdictives aux donnes exprimentales permet de juger de la pertinence des
modles de propagation proposs.
Ces aspects sont dvelopps dans la suite du mmoire au cours de trois chapitres traitant
respectivement de la mthode des lments finis tendus, de la campagne exprimentale, et de
l'analyse de ces essais.
62
Chapitre II.
dynamique non-linaire
Sommaire
II.1. Introduction..................................................................................................................
...64
II.2. tat de l'art...............................................................................................................
.......64
II.2.1. Gnralits sur la mthode des fonctions de niveau (Level set method).................64
II.2.2. Mthode des lments finis tendus (X-FEM).........................................................70
II.2.3. Conclusions..........................................................................................
.....................75
II.3. Actualisation des fonctions de niveau sur une grille auxiliaire.................................75
II.3.1. Principe de la mthode...................................................................................
...........75
II.3.2. Algorithme pour la propagation de fissures.............................................................77
II.4. Prise en compte des non linarits matriau...............................................................83
II.4.1. Remarques gnrales..............................................................................................
...83
II.4.2. Technique d'intgration...................................................................................
..........83
II.4.3. Conservation de l'nergie en dynamique non linaire..............................................86
II.5. Exemples d'application......................................................................
............................88
II.5.1. Translation d'une bulle de gaz..................................................................................88
II.5.2. Cas de Broberg............................................................................................
..............88
II.5.3. Exprience de Kalthoff.....................................................................
........................92
II.5.4. Plaque troue impacte........................................................................
.....................95
II.6. Conclusion............................................................................................
...........................98
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
II.1. Introduction
La mthode des lments finis tendus (X-FEM) a t brivement introduite au cours du
chapitre bibliographique. Le second chapitre entre dans le dtail de la mthode. Un tat de l'art
approfondi sur le sujet est prsent dans la premire partie. Ensuite les deux parties suivantes
abordent des techniques spcifiques la mthode, mises en oeuvre durant cette thse. Il s'agit
de l'actualisation des fonctions de niveau sur une grille auxiliaire d'une part, et de la prise en
compte du comportement non-linaire du matriau d'autre part. Enfin, ce chapitre se conclut par
quelques exemples numriques dmontrant l'efficacit de la mthode et des techniques
proposes.
64
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
t t
t 0 0 0 0 0 0 0 0
0.2 1 1 1 0.4 0 0 0
(a) (b)
Fig. 38 : (a) Mthode des "marqueurs". (b) Mthode des volumes finis.
x
n
xp
interface
65
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
L'quation de propagation de l'interface est de type Hamilton-Jacobi car on peut l'crire sous la
forme gnrale donne par l'quation (77).
= H x , , ,... (77)
t
Un point essentiel de la mthode rside dans la construction de schmas numriques appropris
la rsolution de ce type d'quation. En effet, la solution numrique doit vrifier la condition
d'entropie, c'est dire qu'elle doit tre la limite des solutions du mme problme augment d'un
terme visqueux, lorsque ce terme visqueux tend vers zro (cf. chapitre 4 et 5 de [106]). C'est
cette condition que l'on peut s'assurer d'obtenir une solution physique.
Considrons l'exemple donn par Sethian dans [106] d'une interface initialement sinusodale,
laquelle on applique un champ de vitesse d'expansion unit (cf. Fig. 40).
66
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
segments) conduit une interface qui se recoupe, formant une "swallowtail" (=queue
d'hirondelle). Cela s'explique en considrant deux points proches du "creux" de la sinusode
dont la trajectoire se coupe. L'introduction d'un terme visqueux (via la courbure ) dans la
vitesse V =1 permet d'viter ce phnomne, mais ne permet toujours pas d'obtenir la
solution physique attendue.
La construction de schmas numriques vrifiant la condition d'entropie passe par l'utilisation
d'approximation dcentre amont, ou "upwind" (=regardant dans la direction d'o vient le vent)
des diffrentes quantits. L'approximation numrique utilise n'a pas besoin d'tre extrmement
prcise, mais elle doit transmettre l'information de l'interface vers les points les plus loigns.
Ainsi, pour approximer le gradient dans le cas d'une vitesse positive, plutt que d'utiliser
l'approximation base sur les diffrences centre (78), on prfre utiliser la relation (79).
i1 i1
Schma des diffrences centres : , x (78)
2x
Schma "upwind" : , x max Di x , 0 min Dx
i , 0
i1 i i i 1 (79)
o l'on dfinit les oprateurs : Dx
i = et Dx
i =
x x
Un exemple comparant les deux approximations est fournie sur la Fig. 41.
On remarque que selon le signe de la vitesse, on bascule d'une expression du gradient une
autre, afin de toujours regarder dans le sens de l'information.
L'intgration temporelle est explicite, et nous utilisons le schma de Runge-Kutta d'ordre 2
dcrit par les relations (81). La premire tape donne le prdicteur n1* , rsultat du schma
d'Euler. La deuxime tape utilise ce prdicteur dans le calcul de H n1 * , permettant ainsi de
corriger le second membre.
67
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
n1* n n
= t H
t (81)
n1 = n H n H n1 *
2
La condition de stabilit du schma gnral (ou condition de Courant-Friedrichs-Lewy) peut
s'crire comme le produit de la condition du schma temporel par celle du schma spatial.
Gottlieb [118] montre que l'utilisation du schma temporel de Runge-Kutta d'ordre 2 dcrit par
les quations (81) est optimum puisqu'il n'affecte pas la stabilit. Par contre, il subsiste la limite
du schma spatial qui dans notre cas est donne par l'quation (82) qu'on peut retrouver dans
[117] ou [119].
x
t (82)
2max { V }
Malheureusement, la proprit de distance signe n'est pas ncessairement conserve. Une
opration de r-initialisation dcrite par l'quation (83), et base sur la rsolution d'une quation
de type Hamilton-Jacobi permet de rajuster la norme du gradient de la fonction de niveau.
sign 1 = 0 (83)
Lorsque l'tat stationnaire est atteint, par dfinition, la drive temporelle s'annule, et la norme
du gradient de la fonction de niveau gale bien l'unit.
Il existe d'autres dtails relatifs la mthode que nous n'exposerons pas dans ce mmoire. Nous
invitons le lecteur intress consulter l'ouvrage de Sethian [106] ou l'article d'Osher [120].
La mthode des fonctions de niveau a t applique avec succs dans des domaines aussi varis
que la mcanique des fluide avec la modlisation de l'coulement de fluide bi-phasiques [119]
ou du phnomne de combustion, la modlisation de la croissance et de la germination de
cristaux [121], l'amlioration d'image [122], etc...
68
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
crack surface:
=0
crack front:
=0
x crack surface
{xx=0
0
(85)
Finalement, grce leurs proprits, ces fonctions de niveau dfinissent un repre local en
pointe de fissure (cf. Fig. 43).
69
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
{
dyn dyn
a= a K 1 , K 2 , ,...
= K dyn dyn
1 , K 2 , ,...
puis
{
V =a cos
V = asin
2. Extension de ces quantits tout le domaine (dans le cas tridimensionnel)
V V
sign V =0 sign V =0
idem pour V
3. Modification de la vitesse du plan de fissure
{}+ V
V = (86)
t V
4. Actualisation et r-initialisation du plan de fissure pendant t
V = 0 sign 1=0
t
5. Actualisation du front de fissure pendant t
V = 0
t
6. Orthogonalisation et r-initialisation du plan de fissure
sign = 0 sign 1 = 0
Duflot [125] propose une lgre modification de la formule d'actualisation de afin de
permettre le traitement des grands angles de propagation. La relation utiliser dans ce contexte
est donne par (87) et introduit 0 le gradient de avant l'actualisation.
{}+ V 0
V = 0 o V = (87)
t t V
Il propose galement une mthode baptise "smooth level set method", mlant quations
gomtriques et itratives de type Hamilton-Jacobi afin d'obtenir des champs continus et lisses.
70
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
{
sin / 2
{
H x = 1 si 0
1 si 0
F j x = r sin / 2 sin
cos / 2
(88)
cos / 2 sin
o reprsente la seconde coordonne du repre local attach la pointe de fissure (cf. Fig.
43). Les inconnues supplmentaires correspondant ces deux types de fonctions
(respectivement notes a et bk pour les fonctions H et Fk) viennent s'ajouter aux degrs de
libert traditionnels (nots u). L'approximation du champ de dplacement devient ainsi :
u x N i xui N i x H x ai N i x F k x b i , k
(89)
i iI H i I F k =1, ..4
Dans l'quation (89), IH et IF dsignent respectivement l'ensemble des noeuds appartenant des
lments ayant t coups, et l'ensemble des noeuds contenant ou ayant contenu la pointe de
fissure.
71
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
qnn
n
q n1 = q n1
n1 n1
qn = (90)
q n =0 q n1
Rthor [34] a dmontr que cette stratgie, associe au schma de Newmark de l'acclration
moyenne permet d'assurer la conservation de l'nergie. L'unique dissipation est celle introduite
par la valeurs des nouveaux degrs de libert l'instant t n1 , et correspond au travail
d'ouverture des lvres de la fissure. Ce rsultat thorique est dmontr en lasto-dynamique
dans [34]. Sa gnralisation au cas d'un matriau non-linaire est fournie au paragraphe II.4.3
de ce chapitre.
A partir des champs HRR ([12] et [13]), Elguedj [107] effectue une analyse de Fourier pour
extraire une base de fonctions reprsentative des champs asymptotiques en pointe de fissure
pour un matriau lasto-plastique. Il propose ainsi la base de 6 fonctions dcrite par l'quation
(91).
1
F
HRR
j x = r n1
{ sin 2 , sin 2 sin , cos 2 , cos 2 sin , (91)
sin sin 3 , cos sin 3
2 2 }
Le coefficient n est le paramtre matriau de la loi de Ramberg-Osgood. Mise part cette
modification de la singularit, deux fonctions de sont ajoutes par rapport au cas lastique.
Schma temporel
Il existe de nombreux schmas temporels pour traiter de la dynamique des structures (cf. [127]
pour une revue dtaille), mais le plus simple et le plus rpandu reste sans doute celui de
Newmark [128]. Ce dernier suppose les champs de vitesse et de dplacement suffisamment
rguliers dans le temps pour que l'on puisse crire le dveloppement limit (92), qui se r-crit
sous la forme (93).
2 3
t t u u
un1 = u n t u n un 6 n1 n
2 6 t
(92)
un1 = u n t u n
t2
2
2 n1 n
t
u u
72
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
2
un1 = un t u n t 2 1 u n t 2 un1
(93)
un1 = un t 1 un t u n 1
Ainsi, il est possible d'crire les diffrentes quantits cinmatiques en fonction d'une seule, et
de rsoudre ainsi l'quation d'quilibre dynamique en fonction de la quantit choisie. Les
conditions de stabilit de ce schma sont bien connues pour la mthode des lments finis
classique (cf. [127] par exemple). Elles se rsument par les relations (94).
1
Pour 2 , le schma est inconditionnellement stable.
2
1 (94)
1
Pour 2
et 2 , le schma est stable si t t c = .
w /2
o w reprsente la plus grande frquence propre de la structure. Celle-ci est dtermine en
pratique par l'estimation de la plus grande frquence propre du plus petit lment. En gnral,
on estime ainsi grossirement t c min { xe } / E / o x e est la taille d'un lment.
Les deux schmas les plus frquemment utiliss sont :
1 1
Avec = 2 , = 4 comme paramtres, le schma implicite de l'acclration moyenne est
inconditionnellement stable.
1
Avec = 2 , =0 , le schma explicite de la diffrence centre ncessite l'utilisation d'un
pas de temps infrieur un pas de temps critique dfini par la structure.
Le choix =1/ 2 s'explique par la prcision du schma qui est de l'ordre de ot 2 dans ce
cas (alors qu'elle est de l'ordre de ot lorsque 1/ 2 ). Dans la littrature, il est conseill
de choisir le schma de l'acclration moyenne lorsque l'on traite de phnomnes relativement
lents, et par opposition, d'utiliser le schma de la diffrence centre pour traiter des problmes
de dynamique rapide.
Dans le cadre des lments finis tendus, l'analyse dynamique explicite pose le problme du
pas de temps critique qui est proportionnel la dimension du plus petit lment du maillage.
Ainsi, un lment "coup" de manire arbitraire peut avoir une dimension tendant vers zro.
(a) (b)
Fig. 45 : Utilisation de la mthode X-FEM en dynamique explicite.
Problme du pas de temps critique rsolu par De Borst [129] (a),
et par Menouillard [110] (b).
73
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
De Borst [129] choisit d'viter la situation en interdisant la fissure de trop s'approcher des
noeuds des lments qu'elle traverse. Menouillard [110] montre que l'utilisation d'une matrice
de masse diagonalise permet d'viter cette question, et d'effectuer des calculs explicites avec la
mthode X-FEM.
Afin de garantir la stabilit du schma utilis sans aucune rserve sur le pas de temps, les
rsultats prsents dans ce mmoire utilisent toujours le schma implicite de l'acclration
moyenne, qui est galement inconditionnellement stable dans le cadre de la mthode des
lments finis tendus (dmonstration disponible dans [34]). De plus, comme nous l'avons dj
signal, ce schma garantit la conservation de l'nergie en lastodynamique.
Technique d'intgration
Il est clair que la prsence de fonction discontinue et/ou singulire dans un lment rend
imprcise l'intgration de cet lment si l'on conserve la mme mthode que si l'lment n'tait
pas enrichi. Mos et al [103] propose donc de sous dcouper les lments enrichis, et
d'appliquer ces triangles une technique d'intgration classique. Samaniego et al [130] utilise
cette mthode dans le cas d'un matriau au comportement non linaire.
(a)
(b)
Fig. 46 : Intgration par sous dcoupage, changement de position des points de Gauss.
(a) Avant la propagation. (b) Aprs la propagation.
Comme l'illustre la Fig. 46, cette technique d'intgration introduit une modification du support
d'intgration pour l'lment contenant la pointe de fissure. Dans un calcul non linaire, cela
signifie qu'il faut projeter les champs sur les nouveaux points de Gauss. Or, un des principaux
avantage de la mthode X-FEM provient du fait qu'elle permet d'viter les oprations de re-
maillage et de projection. Il semble donc ncessaire d'adopter une autre stratgie d'intgration.
Ventura [131] montre qu'une intgration exacte des fonctions "saut" H est possible sans sous-
dcoupage explicite des lments enrichis. L'approche dveloppe par l'auteur est attrayante,
cependant, notre connaissance, le cas des fonctions singulires n'est pas traite ce jour.
74
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
II.2.3. Conclusions
La mthode des lments finis tendus est une mthode moderne et efficace pour simuler la
propagation de fissure. La reprsentation de la fissure est implicite et s'appuie sur la mthode
des fonctions de niveau (ou "Level Set method"). Cette description est ensuite utilise pour
enrichir le champ de dplacement (mthode des lments finis tendus ou "X-FEM").
Cette mthode numrique permet de raliser avec succs des simulations efficaces de
propagation arbitraire de fissure dans un domaine d'applications varies : fatigue (cf. [107] et
[108]), rupture fragile en lasto-dynamique (cf. [34] et [110]), multi-fissuration (cf. [112] et
[113]), etc...
Bien qu'elle constitue une avance notable dans le domaine de la simulation numrique de
propagation de fissure, certains points de la mthode peuvent encore tre amliors. Pour deux
d'entre eux, nous essayons de donner une rponse dans les prochains paragraphes.
En premier lieu, l'actualisation des fonctions de niveau est issue des mthodes de diffrences
centres, et son application au cas de maillage non-structur n'est pas immdiate, notamment
lorsque l'on considre des lments distordus ou dont la taille dans l'espace varie rapidement.
Par exemple considrons un maillage contenant la fois des lments gros et petits. Le pas de
temps critique tant proportionnel la taille du plus petit des lments, l'actualisation ncessite
un grand nombre d'itrations. De plus, la reprsentation des fonctions de niveau dans les
lments de grosse dimension est imprcise.
Une rponse simple et efficace ce type de problme consiste introduire une grille auxiliaire
diffrente du maillage de la structure sur laquelle l'actualisation des fonctions de niveau est bien
matrise. Cette technique est dveloppe dans le paragraphe II.3 de ce chapitre et dans [132].
En second lieu, les fonctions enrichies n'tant pas polynomiales, leur intgration numrique
pose problme. Classiquement, un sous dcoupage en triangles munis d'un nombre lev de
points de Gauss est utilis. Cela introduit des changement de points de Gauss proximit de la
pointe de fissure, et par consquent des projections de champs dans une zone o leur gradient
est particulirement lev. Cette tape est fortement pnalisante sur la qualit d'un calcul non
linaire.
Ainsi nous exposons dans le paragraphe II.4 de ce chapitre une technique d'intgration non
conforme, dveloppe dans [107] et [132], permettant d'obtenir une intgration prcise des
matrices lmentaires malgr la prsence des fonctions d'enrichissement, tout en vitant les
projections de champs au voisinage de la pointe de fissure.
Enfin, le paragraphe II.5 de ce chapitre donne quelques exemples d'applications des techniques
exposes.
75
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
Fig. 47 : Dfinition d'une interface pour reprsenter implicitement une inclusion ou un trou
(tir de [133]).
D'autre part, l'actualisation des fonctions de niveau sur un maillage non structur ncessite
l'utilisation d'un schma numrique de type Petrov-Galerkin dvelopp dans [134], relativement
plus complexe que celui que nous avons prsent au paragraphe II.2.1.
Partant du constat que la mthode des fonctions de niveau a initialement t dveloppe dans un
contexte de diffrences finies, nous proposons d'introduire une grille auxiliaire structure et
rgulire, indpendante du maillage mcanique (cf. Fig. 48).
Cette grille reprsente le support pour la dfinition et l'actualisation des fonctions de niveau.
Elle couvre la zone o est susceptible de se propager la fissure et sert exclusivement sa
reprsentation. C'est pourquoi, nous qualifions cette grille d'auxiliaire (cf. [132]).
Le calcul mcanique donne le champ de vitesse (il s'agit des deux scalaires a et en deux
dimensions). Cette donne est utilise sur la grille auxiliaire pour actualiser les fonctions de
niveaux en suivant l'algorithme expos par les relations (86). Ensuite les nouvelles fonctions de
niveau sont projetes sur le maillage de la structure en utilisant une approximation de type
lment fini. Par exemple, la valeur de au point p de coordonne x p est dtermine par la
relation (95).
x p = N i x p i (95)
i
Afin de prendre en compte la courbure ventuelle de la fissure, on constate qu'il est prfrable
d'utiliser un raffinement suffisant pour la grille auxiliaire. Typiquement, on remarque qu'une
taille d'lment de la grille gale un dixime du rayon de courbure suffit pour obtenir une
erreur sur la reprsentation de la fonction de niveau de l'ordre de quelques pourcents.
Il semble ncessaire d'assurer une prcision suffisante concernant la reprsentativit des
fonctions de niveau sur la grille auxiliaire. En effet, ces donnes sont les sources d'information
de l'tape de projection (95). Pour tre clair, il ne semble pas raisonnable d'envisager sur la
grille auxiliaire une prcision infrieure celle du maillage mcanique. On utilise donc la rgle
pratique (96).
76
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
(a)
(b)
77
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
Des solutions ont t proposes dans la littrature. Sussman et al ([119] et [135]) contraignent
par exemple la prservation du volume lmentaire de fluide contenu dans une cellule pendant
la r-initialisation en ajoutant dans l'quation d'volution un terme correctif.
Afin de rester dans une dmarche simple, on propose d'utiliser une fonction signe diffrente de
celles dveloppes par Peng [136] ou Sussman [119]. La fonction signe utilise pendant l'tape
de rinitialisation (mais aussi pour l'orthogonalisation) indique de quel cot de l'interface on se
situe, et permet d'adapter le schma numrique en consquence. La dfinition mathmatique de
la fonction signe est trs simple (cf. quation (97)), mais aussi trs discontinue.
1
sign = 0
1
{ si
si
si
0
=0
0
(97)
C'est pourquoi Peng [136] et Sussman [119] introduisent diffrents lissages proximit de
l'interface. La fonction signe de Peng s'crit par exemple :
sign = (98)
2 x 2
2
=0 =0
=0 =0
A B A B
(a) (b)
Fig. 49 : (a) R-initialisation ne maintenant pas la position de l'interface. (b) R-initialisation
maintenant la position de l'interface.
Le maintien de la position de l'interface peut tre vrifi en assurant la relation (99) illustre par
la Fig. 49 (b).
A
= A (99)
B B
Pour vrifier cette condition lors de la r-initialisation, on introduit une fonction signe
78
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
{
1 si x
sign = si x x (100)
x
1 si x
Exemple d'application
On considre un carr de dimension [1:1][1 :1] avec une taille d'lment de x=0.05
soit une grille de 41x41 points, ce qui est relativement grossier.
Un premier cas de propagation du front de fissure est dfinies par la vitesse, l'incrment de
temps et l'angle suivant :
V t
1 0,2 45
Une reprsentation graphique de ces donnes est propose sur la Fig. 50.
a t
79
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
(a)
(b)
(c)
Fig. 51 : Comparaison avec la thorie (en rouge) des isovaleurs de et obtenus avec
diffrents algorithmes utilisant la fonction signe de Peng.
(a) Stolarska. (b) Gravouil scalaire. (c) Gravouil vectoriel.
80
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
La Fig. 52 prsente les rsultats obtenus avec l'utilisation d'une fonction signe proportionnelle.
(a)
(b)
Fig. 52 : Comparaison avec la thorie (en rouge) des isovaleurs de et obtenus avec
diffrents algorithmes utilisant une fonction signe proportionnelle.
(a) Gravouil scalaire. (b) Gravouil vectoriel.
Les rsultats obtenus sont nettement meilleurs et un second cas test dfini par le tableau ci-
dessous montre galement le bon comportement de l'algorithme mme pour de trs grands
angles de bifurcations (cf. Fig. 53).
V t
1 0,3 70
81
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
Fig. 53 : Comparaison avec la thorie (en rouge) des isovaleurs de et obtenus avec
l'algorithme de Gravouil vectoriel utilisant une fonction signe proportionnelle.
L'utilisation de la fonction signe proportionnelle permet d'obtenir un gradient proche de l'unit
au voisinage de la pointe de fissure (cf. Fig. 54).
Il est important de noter qu'tant donn une fonction , il est impossible de construire une
fonction qui lui soit orthogonale et dont le norme du gradient vale 1. On favorise donc
volontairement la premire proprit, la seconde tant vrifie du mieux possible. C'est
pourquoi il existe une zone tire (en bleu sur la Fig. 54) et une zone contracte (en rouge sur la
Fig. 54) pour la norme du gradient de . Cependant, les rsultats sont satisfaisants, puisque la
zone en avant du front de fissure vrifie parfaitement les deux proprits souhaites.
82
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
83
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
Pour satisfaire ces deux exigences, nous dveloppons ici une technique d'intgration
alternative, galement expose dans [107] et [132]. Elle consiste sous-dcouper de manire
non conforme les lments enrichis et proches de la pointe de fissure.
En premier lieu, il faut tre capable de slectionner l'ensemble des lments que l'on souhaite
sous-dcouper. Il s'agit du regroupement :
des lments dont au moins l'un des noeuds est enrichi,
et des lments susceptibles de bientt s'enrichir pour lesquels on souhaite viter une
projection de champs.
Cette deuxime catgorie est constitue des lments dont l'un des noeuds vrifie une condition
de proximit avec la pointe de fissure (cf. quation (101)).
2 2 1/ 2 r p (101)
En utilisant le critre de Von Mises dans la solution lasto-statique, on obtient une estimation
du rayon plastique [6]. La formule (102) donne cette estimation pour deux angles particuliers en
dformation plane.
2
r p =0
1 K Ic
2 Y
12 2
(102)
2
r p =45
1 K Ic
4 Y
Une fois ces lments slectionns, on les sous-dcoupe. Ainsi, un quadrangle 4 noeuds est
divis en 16 sous quadrangles, amenant le nombre total de points d'intgration 64. La Fig. 55
rsume cette transition.
L'opration de projection de champ est donc ralise au niveau lmentaire pendant que le
matriau est encore lastique. Ainsi, la conservation de la masse et de l'nergie dpend
uniquement des termes mij et kij des matrices de masse et de raideur qui correspondent aux
degrs de liberts standards. Si ces derniers ne sont pas modifis lorsque l'on passe d'un
lments standard ( 4 points d'intgration) un lment sous-dcoup ( 64 points de Gauss),
alors on vrifie mieux la conservation de la masse et de l'nergie. L'approximation introduite est
de l'ordre de grandeur de l'erreur d'intgration de ces termes de masse et de raideur, c'est dire
trs faible.
84
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
rp
(a)
rp
(b)
85
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
(104)
X A[ X ] = [ 12 X ]
T T
AX si A est symtrique
n1
qn =
qnn
q n =0
et q
n1
n1 = q nn1
q n1 (106)
Cette tape intermdiaire de projection fait apparatre une force de re-fermeture des lvres :
[
M Mn
] [ ]
n ext
q n Bn Fn
q =0
T
n
=
M M
n B F crack
n (107)
M n 1 q n1 Bn1
n
Il est de premire importance de noter que l'initialisation des nouveaux degrs de libert 0
permet de conserver l'quilibre exact l'instant n (premire ligne de l'quation (107)), et dfinit
explicitement la force de re-fermeture des lvres.
crack n1 T
Fn = [ M ] qn [ B ] n (108)
86
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
T T
[1/ 2 q M q] [ q B ] 1/ 2 [ ]
T ext ext
= [ q] F 1/ 2[ F ]
(111)
T
1/ 2 [ q]T [ M ] q n [ B] n 1 /2 [ q ]T [ M ] qn [ B ] n T
t / 22 [ q] M n1 [ q ]
Avec la dfinition du travail d'une force entre les instants t n et t n1 donne par (104),
n 1
W n n1 = qt F t dt [ q ]T F 1 / 2[ F ] (112)
n
w erreur
w def
87
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
88
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
{
dyn 1
K t , a 0.9386 si a=500ms
stat = (115)
K a 0.6532 si a=1500 ms1
4 12 u
o u y x = KI x
E 2
(116)
u y xi xi
conduit
u
KI
E 2 i
41 2
xi
i
89
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
J
dyn
= lim
1 0
[ w def cin
]
w kj ij ui ,k q k , j u i ui ,k ui u i ,k q k d (117)
C
d
J dyn
Attigui = w df w cin jk ij ui ,k qk , j d d t ui ui ,k qk d (118)
Pour l'implantation numrique dans Cast3M de J dynAttigui , le dernier terme est valu en faisant la
diffrence des quantits aux instants t et tt .
Rsultats et discussion
La simulation est ralise en considrant une plaque carre de grande dimension (demi cot
gal 100m). La taille des lments en pointe de fissure est de 0.25m. La longueur initiale de
la fissure est de 1m. La plaque est initialement mise en traction de manire quasi-statique avant
d'imposer la vitesse de propagation. L'avance s'effectue par saut discret d'un lment par pas
de temps.
L'expression de J dyn (cf. quation (117)) fait intervenir des termes multiplis par le champ q,
non nuls sur un domaine vanouissant 1 reprsent sur la Fig. 58. Certains auteurs estiment
que les champs en pointe de fissure sont de trop mauvaise qualit pour tre utiliss et
conseillent de n'intgrer que sur C . Afin de lever toute ambigut sur la prise en compte ou
non de ce domaine, on compare les rsultats fournis avec et sans ce domaine.
90
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
(a)
(b)
(c)
91
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
(a) (b)
Fig. 61 : (a) Gomtrie de l'exprience de Kalthoff. (b) Trajet de fissure en fonction de la
vitesse d'impact et du rayon de fond d'entaille.
92
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
Il s'agit de la rupture d'une plaque possdant deux entailles symtriques suite l'impact d'un
projectile lanc la vitesse V. Pour des vitesses d'impact relativement basses, la rupture
s'effectue par ouverture, tandis que des bandes de cisaillement s'observent pour des vitesses
plus leves. Nous nous restreignons au premier cas de figure en adoptant V=16.5m/s comme
dans [34] et [140] . Les rsultats exprimentaux fournissent un trajet de fissure inclin d'environ
70 avec l'horizontale. Kalthoff prcise que dans ce cas, la rupture de l'prouvette est complte.
Le matriau tudi est un acier haute rsistance X2 NiCoMo18 9 5 (18Ni1900), dont la limite
d'lasticit peut tre value 2.1GPa. Puisque celle-ci est trs haute, on peut modliser le
comportement de cet acier comme lastique linaire avec pour caractristiques :
E K Ic
190 GPa 0,3 8000 kg.m-3 68 MPa.m0.5
Modle de propagation
Le modle de propagation utilis est bas sur la mcanique de la rupture linaire. La direction
est donne comme tant celle de la contrainte d'ouverture maximale. En utilisant les champs
asymptotique, on trouve que l'angle de propagation peut s'exprimer en fonction de KI et KII (cf.
formule (119)).
[ [ ] ]
2 1/2
1 KI KI
= 2atan sign K II 8 (119)
4 K II K II
En ce qui concerne la vitesse de propagation, plutt que de l'imposer comme dans [34], on
choisit d'utiliser le critre empirique de Kanninen [8] exprim dans la formule (120), o l'on
choisit les paramtres Gc =22E3 N.m1 , m=1 et v lim=c R . La tnacit Gc est estime
depuis les donnes fournies dans [88].
1/ 2 1 /m
[ ]
a = vlim 1
Gc
J dyn
(120)
Rsultats de la simulation
La symtrie nous permet de ne discrtiser que la moiti de l'prouvette. On choisit un maillage
uniforme avec une taille de maille de 1mm, ce qui conduit un total de 10000 lments. Pour
ce cas extrmement simple, le maillage de la structure est galement le support pour
l'actualisation des fonctions de niveau. Le pas de temps du calcul est pris gal 1s pendant la
mise en charge, puis 2s pendant la propagation.
La Fig. 62 reprsente la norme de la vitesse sur le maillage dform (a), et la contrainte de Von
Mises (b) pour diffrents instants de la propagation.
93
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
t=5s t=5s
t=27s t=27s
t=55s t=55s
94
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
t=95s t=95s
(a) (b)
Fig. 62 : Rsultats de la modlisation de l'exprience de Kalthoff.
(a) Norme de la vitesse sur maillage dform (x10). (b) Contrainte de Von Mises.
On constate que la contrainte dpasse rarement la limite d'lasticit validant l'hypothse de
matriau lastique.
Le trajet de fissure est correctement prdit dans l'espace. L'angle de la fissure avec l'horizontale
est d'environ 75 au dbut de la propagation, et de 55 la fin de la propagation. Les rsultats
exprimentaux donne une valeur moyenne de 70.
Concernant l'volution temporelle, aucune information exprimentale sur la vitesse d'avance
de la fissure n'est disponible la connaissance de l'auteur.
95
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
(a) (b)
Fig. 64 : Vitesse sur les bords gauche (a) et droit (b) de l'prouvette.
Modle de propagation
Le modle de propagation adopt pour simuler la propagation de la fissure est le mme que
celui de l'exemple prcdent (cf. quations (119) et (120)).
96
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
Rsultats de la simulation
On compare dans un premier les rponses fournies par une modlisation du comportement du
matriau comme tant lastique, et une modlisation avec un comportement lasto-plastique
parfait. Comme on le constate sur la Fig. 56 (a) et la Fig. 45, les deux rponses sont quasiment
les mmes. Cela permet de conclure que l'hypothse de plasticit confine est particulirement
bien adapte ce type de configuration.
De plus, le critre de propagation utilis permet de reproduire l'amorage, l'arrt et le re-
dmarrage de la fissure dans un contexte fortement marque par la dynamique. L'volution dans
le temps et dans l'espace de la propagation de fissure obtenue par simulation numrique est trs
proche des rsultats exprimentaux.
(a)
(b)
97
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
fissure varie relativement peu, seul l'arrt semble moins marqu. Mais dans un second temps,
pour des valeurs de Y infrieures 15MPa, la vitesse chute brutalement et la fissure propage
peine (cf. Fig. 56(b) et Fig. 45). Ce changement de comportement traduit l'influence de la
plasticit. En effet, une si faible limite d'lasticit entrane une trs forte plastification dans les
zones cisailles par l'impact de la barre d'entre ( gauche) et autour du trou. La plasticit
devient une source si importante de dissipation d'nergie que la pointe de fissure ne parvient pas
avancer.
II.6. Conclusion
La mthode des lments finis tendus a t expose de manire approfondie dans ce chapitre.
Elle permet de simuler efficacement la propagation arbitraire de fissure sans avoir remailler la
structure. Elle repose sur deux ingrdients :
la description implicite de la fissure l'aide d'un couple de fonctions de niveau orthogonales,
l'enrichissement du champ de dplacement avec des fonctions discontinue et asymptotiques.
Afin d'amliorer encore cette mthode, quelques ides ont t avances dans ce chapitre.
Nous avons d'abord propos l'introduction d'une grille auxiliaire dcouple du maillage de la
structure permettant une actualisation efficace des fonctions de niveau. En effet, la mthode des
fonctions de niveau perd de son efficacit et de sa simplicit lorsque l'on cherche l'appliquer
un maillage non structur ou irrgulier. La grille auxiliaire constitue un support de dfinition
rgulier sur lequel l'actualisation des fonctions de niveau est simple et bien matrise. Nous
avons galement propos l'utilisation d'une fonction signe proportionnelle afin de maintenir la
position de l'iso-zro durant l'tape de r-initialisation.
L'intgration numrique des lments enrichis est classiquement ralise en les sous-dcoupant
en triangles. Dans un calcul non-linaire, cette technique introduit des projections de champs en
pointe de fissure. Nous avons donc expos une technique de sous-dcoupage non conforme qui
utilise un grand nombre de points d'intgration par lment. D'une part, cette technique permet
d'intgrer les fonctions non polynomiales utilises dans l'enrichissement. D'autre part, en sous-
dcoupant les lments proches de la pointe de fissure, on effectue la transition d'un lment
standard en lment sous-dcoup avant que la plasticit ne l'affecte, et on vite ainsi la
projection des champs inlastiques.
98
Chapitre II. La mthode des lments finis tendus (X-FEM) en dynamique non-linaire
L'implantation de la mthode des lments finis tendus et des techniques prsentes dans ce
chapitre a t ralise dans Cast3m. Les exemples d'application montrent l'efficacit de la
mthode pour simuler la propagation arbitraire et l'arrt de fissure en dynamique non linaire.
Ces dveloppements numriques ainsi que le dernier exemple ont fait l'objet d'une publication
dans IJNME [132].
99
Chapitre III.
acier 16MND5
Sommaire
III.1. Introduction...................................................................................................
..............102
III.2. Prsentation du matriau............................................................................
...............102
III.2.1. Caractristiques gnrales.....................................................................................102
III.2.2. Effet de la vitesse de dformation sur le comportement.......................................104
III.3. Protocole exprimental.....................................................................................
..........107
III.3.1. Gnralits..................................................................................
...........................107
III.3.2. Mesure de la vitesse de propagation de fissure.....................................................108
III.4. Essais sur CT.........................................................................................
......................111
III.4.1. Gomtrie de l'prouvette...........................................................................
...........111
III.4.2. Donnes de prfissuration et d'amorage..............................................................112
III.4.3. Vitesse de propagation....................................................................................
.......116
III.5. Essais sur anneau..............................................................................
..........................119
III.5.1. Gomtrie de l'prouvette............................................................................
..........119
III.5.2. Donnes de prfissuration et d'amorage..............................................................119
III.5.3. tude de la propagation....................................................................
.....................121
III.6. Analyse du facis de rupture....................................................................................
..123
III.6.1. Objectifs des observations...................................................................
..................123
III.6.2. Observations........................................................................................................
...123
III.7. Conclusions.............................................................................................
.....................125
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
III.1. Introduction
Ce chapitre prsente l'ensemble des travaux exprimentaux raliss au cours de cette thse
concernant l'tude des phnomnes de propagation dynamique de fissure de clivage.
Le matriau tudi est un acier 16MND5 dj largement caractris lors de prcdents travaux
au CEA [141] et [17]. Afin de complter la connaissance de son comportement, une
identification aux grandes vitesses de dformation est mene. La premire partie de ce chapitre
prsente les principales caractristiques mcaniques du matriau.
La seconde partie de ce chapitre dcrit le protocole exprimental de cette campagne d'essais, en
dtaillant notamment les moyens de mesure de vitesse de propagation de fissure mis en oeuvre.
Plusieurs gomtries sont considres : des prouvettes CT de type CT25 mais d'paisseur
rduite (5 et 10mm), des anneaux fissurs en compression dont certains sont adapts pour
amorcer et propager en mode mixte. Les rsultats exprimentaux de ces deux gomtries font
l'objet des troisime et quatrime parties de ce chapitre.
Enfin, on prsente, dans une dernire partie, les analyses fractographiques des facis de rupture
qui dcrivent les mcanismes de ruine en jeu.
Il provient d'une virole de cuve de racteur eau pressurise (REP) franais (correspondant au
code matriau LISN 520). Le prlvement des prouvettes est effectu de manire les
solliciter dans le sens circonfrentiel de la cuve d'origine, et faire propager la fissure dans le
sens de l'paisseur.
{
si Y
E
= 1over n o Y = ' Y E ' 0 (122)
E
0
' Y
E'
' Y
K si Y
Les paramtres E , ' Y , 0 , E ' , K ' , et n sont des fonctions de la temprature donnes par
102
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
{
E = a EbE cE 2 d E 3
2 3 4
' Y = a ' b ' c ' d ' e '
Y Y Y Y Y
0 = 100a b [ atanc d e ] T C
0 0 0 0 0
o = (123)
E ' = a E ' bE ' [ atan c E ' d E ' eE ' ] 273
K ' = aK ' b K ' [ atan c K ' d K ' e K ' ]
1/ n = an bn [ atan cn d n en ]
La comparaison entre ce lissage et les essais de traction est donne par la Fig. 67.
Fig. 67 : Courbes de tractions lisses (trait continu) compares aux courbes exprimentales
(symboles) pour diffrentes tempratures (tir de [17]).
La conductivit et la capacit calorifique du matriau ont galement t values pour cette
gamme de temprature dans [17].
Proprits rupture
Plusieurs campagnes d'essais sur diffrentes gomtries et pour diffrentes sollicitations ont t
ralises sur cet acier au laboratoire. Cela a permis d'tudier en dtail l'amorage d'une fissure
en clivage et de l'ensemble des phnomnes associs (comme l'effet de triaxialit, le Warm Pre-
Stressing, ou la spcificit d'une sollicitation thermique, ...). La courbe de transition de cet acier
103
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
a ainsi pu tre dtermine. La Fig. 68 compare les donnes exprimentales d'essais sur CT25
aux "Master Curves".
(a)
(b)
Fig. 68 : Comparaison entre la tnacit exprimentale (sur CT25) et les courbes de transition
"Master Curve" pour l'acier 16MND5.
On remarque que tous les essais effectus une temprature infrieure -100C amorcent en
clivage, contrairement aux essais effectus pour une temprature situe entre -75C et -50C
qui prsentent parfois un amorage ductile. Le clivage disparat entre -50C et 25C, laissant
place une rupture ductile caractrise par un plateau valu 320MPa.m-0.5.
104
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
pour la campagne d'essais, des essais sur barres de Hopkinson ont t raliss par le Laboratoire
de Mcanique du Solide de l'cole Polytechnique de Paris (cf. [142] et [143]).
Le principe des essais sur barres de Hopkinson, l'identification d'une loi de comportement
lasto-viscoplastique permettant de prendre en compte l'adoucissement du matriau d
l'chauffement adiabatique, ainsi que l'intgration de ce comportement dans Cast3m, sont
dcrits dans l'annexe B. Nous ne prsentons que les principaux rsultats dans ce paragraphe.
(a)
(b)
105
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
[ ]
1
,0 ,T 1H
in in stat in in p
, ,T =
(124)
{
a=0.07592 0.07182T / 273.
o { p=4.250.015T
in in
H T , = abexp / c
et b=0.06352 0.2397T / 273.
c=0.02528 0.04124T / 273.
Deux modifications ont t apportes la loi originale de Symond-Cowper. D'une part, les
coefficients H et p dpendent de la temprature afin de reproduire au mieux les essais sur toute
la gamme de temprature tudie. D'autre part, le coefficient H dpend galement de la
dformation inlastique, afin de reproduire l'effet d'adoucissement thermique. Le modle ainsi
identifi peut tre compar aux rsultats exprimentaux sur la Fig. 70 pour deux tempratures et
pour l'ensemble des essais dans l'annexe B.
(a)
(b)
106
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
, in , in =
3
2
sij sij y in , in
(125)
o sij = ij 13 kk ij
in int t int
= (126)
t
Le lecteur intress pourra se reporter l'annexe B, qui donne plus de dtails sur la mise en
oeuvre dans le calcul lment fini de cette loi de comportement lasto-viscoplastique.
Mise en temprature
Les prouvettes sont refroidies par injection d'azote dans une enceinte isolante. Il faut en
gnral environ 1 heure pour atteindre la temprature de consigne (-125C). Ensuite, afin
d'assurer une temprature homogne dans l'prouvette, une tape de stabilisation est observe.
La dure minimale de celle-ci est de 30 minutes pour les CT, et d'1 heure pour les anneaux. La
temprature est mesure l'aide de thermocouple de type K ou T, directement fix sur
l'prouvette. L'ensemble des essais est ralis -125C, afin de s'assurer d'un amorage
systmatiquement par clivage. Par comparaison, la temprature de rfrence T0 dfinissant la
master curve pour ce matriau est de -120.5C [144].
Mise en charge
Lessai de traction rupture est quant lui pilot partir du logiciel WAVEMAKER
dINSTRON. Il permet d'enregistrer en cours dessai :
le dplacement du vrin,
la force,
louverture ( laide dun capteur douverture au pralable talonn),
(et ventuellement la tension d'une jauge brins).
107
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
clip gage
crack gage
thermocouple
108
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
5.1mm
0.25mm
25.4mm
1.27mm
(a) (b)
Fig. 72 : Jauges brins Vishay de type CPD (a) et CPA (b).
L'avance de la fissure produit la rupture des brins de la jauge dispose sur son trajet. Lorsqu'un
brin se rompt, la rsistance globale de la jauge diminue et la tension mesure augmente. Cette
tension est enregistre au cours du temps l'aide d'un PC quip de cartes d'acquisition rapide.
Ces dernires peuvent enregistrer jusqu' environ 17 millions points une frquence allant
jusqu' 60MHz, soit une dure de 0,28 secondes. En divisant la frquence, on augmente la
dure d'acquisition du mme facteur. Un programme d'acquisition a t cr l'aide du logiciel
LABVIEW.
L'une des voies de mesure est utilise comme trigger (pour dclencher l'acquisition). Comme la
carte d'acquisition rapide fonctionne avec des entres simples (et pas diffrentielles), des
problmes de rfrence de masse ont surgi. Un re-cblage des botiers lectriques, ainsi que
l'utilisation d'un PC annexe ddi au dclenchement du trigger ont permis de les rsoudre.
D'autre part, les jauges de diffrents constructeurs ont t testes. Il s'est avr que les jauges
Vishay fournissent un signal moins bruit que les jauges HBN. Le schma lectrique gnral
est dcrit par la Fig. 73.
Rs
Vin=+15V
Gage
......Rbi......
+Vout
Rp
Rei
Pour les jauge de type CPD (Fig. 72), on ajoute la rsistance Rb de chaque brin une rsistance
109
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
monte en srie note Re, afin de linariser le signal de sortie en fonction du nombre de brins
rompus.
Pour les autres types de jauges (comme le type CPA de la Fig. 72), le schma lectrique est le
mme mais avec les rsistances Re nulles, puisqu'il est impossible de cbler indpendamment
les brins. Ainsi, le signal de ce type de jauges n'est pas linaire avec le nombre de brins rompus.
Nanmoins, ces jauges ont un pas relativement fin (20 brins sur 5 ou 10mm selon le modle).
On les a donc utilise lorsque la longueur d'avance de fissure tait trop petite pour tre
mesure prcisment par les jauges de type CPD (20 brins sur 20mm).
Les rsistances propres des brins Rb valent environ 110 ohms, tandis que celles notes Rp et Rs
pour "parallle" et "srie" ont respectivement les valeurs de 82 et 180 Ohms, dtermines afin
d'utiliser pleinement la plage d'entre admissible des cartes (jusqu' 5V).
Exemple de rsultats
L'analyse des donnes brutes issues de la mesure s'effectue deux chelles de temps.
Pendant la mise en charge, on observe la rupture des brins qui recouvrent la prfissuration ou en
sont trs proches (cf. Fig. 74(a)).
Lors de la rupture fragile, la rupture des brins se fait grande vitesse (cf. Fig. 74(b)). On note
une bonne simultanit des signaux fournis par les deux jauges places sur chaque face de
l'prouvette.
(a) (b)
Fig. 74 : Exemple de rsultats exprimentaux obtenus pour la ct10ov.
(a) Courbes force-temps et tension de la jauge-temps (basse frquence).
(b) Signaux enregistrs par deux jauges situes sur les deux faces de la mme prouvette CT
(haute frquence).
A partir de la frquence d'acquisition et des numros des points chantillonns, on en dduit
l'instant de rupture de chacun des brins. Connaissant la position de chaque brin, on obtient la
courbe longueur de fissure-temps (cf. Fig. 75).
110
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
12.5mm
60mm
10mm
35mm
60
50mm
60mm
111
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
fatigue fatigue
precracking precracking
crack
branching
straight crack
propagation
(a) (b)
Fig. 77 : (a) Trajet de fissure rectiligne. (b) Trajet de fissure branche.
Le tableau suivant et la Fig. 78 rcapitulent les donnes exprimentales concernant les
longueurs de fissure l'amorage et l'arrt, ainsi que le type de trajet de fissure (rectiligne ou
branch).
Mesure
n essai paisseur Fissure
prouvette a0 (mm) a0/w d'avance af (mm)
LISN BN (mm) branche?
de fissure?
05-217 CT10_520OY 10 25,25 0,51 NON OUI 48,5
06-064 CT10_520OZ 10 24,81 0,5 NON OUI 48
06-065 CT10_520PC 10 24,76 0,5 OUI OUI 50
05-216 CT10_520PA 10 24,93 0,5 OUI OUI 50
06-066 CT10_520OU 10 25,34 0,51 OUI OUI 50
05-214 CT10_520OV 10 25,05 0,5 OUI OUI 50
06-069 CT5_520PM 5 24,2 0,48 NON OUI 48,5
05-212 CT5_520PJ 5 24,72 0,49 NON OUI 49,5
05-215 CT5_520PE 5 25,22 0,5 OUI NON 50
05-211 CT5_520PF 5 24,49 0,49 OUI NON 50
06-067 CT5_520PK 5 24,25 0,49 OUI OUI 50
05-213 CT5_520PI 5 24,48 0,49 OUI NON 50
06-063 CT5_520PG 5 24,69 0,49 OUI OUI 50
112
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
La longueur de pr-fissuration est calcule comme tant le rapport entre l'aire pr-fissure et sa
largeur, constituant ainsi une longueur moyenne. Elle s'lve environ 2mm pour les
prouvettes de largeur B=10mm et 1.5mm pour celles de largeur B=5mm, l'exception de
l'prouvette CT5_520PD qui ne semble pas pr-fissure du tout.
CT10_520OY
CT10_520OZ } fissure rectiligne
CT10_520PC
CT10_520PA
CT10_520OU
CT10_520OV
} fissure branche
CT5_520PM
CT5_520PJ
} fissure rectiligne
CT5_520PE
CT5_520PF
CT5_520PK
CT5_520PI
CT5_520PG
} fissure branche
a0 (mm)
23 28 33 38 43 48
longueur de fissure (mm) af (mm)
Donnes d'amorage
Les enregistrements effectus lors de la mise en charge permettent de tracer les courbes force-
ouverture pour les diffrentes prouvettes (cf. Fig. 79).
L'ensemble des donnes relatives l'amorage sont rsumes dans le tableau suivant.
n essai paisseur F rupture K1 J total KJ
prouvette a0 (mm) a0/w
LISN BN (mm) (kN) (MPam) (kJ/m) (Mpam)
05-217 CT10_520OY 10 25,25 0,505 8,99 39,42 6,92 39,95
06-064 CT10_520OZ 10 24,81 0,496 11,62 49,62 11,04 50,44
06-065 CT10_520PC 10 24,76 0,495 20,34 86,58 37,22 92,63
05-216 CT10_520PA 10 24,93 0,499 25,75 110,75 68,52 125,69
06-066 CT10_520OU 10 25,34 0,507 26,04 114,89 78,67 134,68
05-214 CT10_520OV 10 25,05 0,501 26,39 114,36 85,49 140,39
06-069 CT5_520PM 5 24,20 0,484 6,14 50,55 11,31 51,06
05-212 CT5_520PJ 5 24,72 0,494 8,59 72,98 25,09 76,05
05-215 CT5_520PE 5 25,22 0,504 13,86 121,38 97,52 149,95
05-211 CT5_520PF 5 24,49 0,490 17,21 144,09 174,49 200,57
06-067 CT5_520PK 5 24,25 0,485 16,61 137,13 196,39 212,78
05-213 CT5_520PI 5 24,48 0,490 17,43 145,85 340,82 280,32
06-063 CT5_520PG 5 24,69 0,494 17,49 148,29 369,36 291,81
06-068 CT5_520PD 5 23,00 0,460 19,48 149,48 392,23 300,71
113
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
(a) (b)
(c) (d)
Fig. 79 : Courbes force - ouverture du clip.
(a) ct5 avec un trajet rectiligne, (b) ct5 avec un trajet branche,
(c) ct10 avec un trajet rectiligne, (d) ct10 avec un trajet branche.
On remarque dj une forte dispersion dans l'instant de rupture qui se confirme dans le calcul
du facteur d'intensit des contraintes gnralis en lasto-plastique (cf. Fig. 80).
Le calcul de ce dernier est ralis l'aide de la norme ASTM E 1820 [145]. Rappelons
nanmoins que cette dernire demande que la condition 2w / B4 soit vrifie, ce qui n'est
pas le cas ici, puisque w=50mm, et B=5 ou 10mm. Si, malgr cela, on applique cette norme, on
peut calculer le facteur d'intensit des contraintes partir de la charge rupture et de l'aire sous
la courbe force-ouverture plastique. Ce type de formule suppose une condition de dformation
plane qui peut tre remise en cause par la faible paisseur de nos prouvettes, mais qui reste
cependant une bonne approximation en gnral.
Un lien fort entre l'instant d'amorage et le trajet de la fissure est mis en vidence ici. En effet,
les prouvettes qui ont rompu de manire rectiligne ont toutes cass plus tt que celles qui
prsentent une fissure branche.
114
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
CT10_520OY
CT10_520OZ } fissure rectiligne K1 (MPa.m)
CT10_520PC KJ (MPa.m)
CT10_520PA
CT10_520OU
CT10_520OV
} fissure branche
CT5_520PM
CT5_520PJ } fissure rectiligne
CT5_520PE
CT5_520PF
CT5_520PK
CT5_520PI
CT5_520PG
} fissure branche
Fig. 80 : Facteur d'intensit des contraintes l'amorage pour les prouvettes CT.
Malgr des conditions de symtrie videntes, la majorit des prouvettes prsentent un trajet de
fissure branch. Pour expliquer ce phnomne, plusieurs lments de rponses peuvent tre
avancs. Tout d'abord, les prouvettes n'ont pas d'entailles latrales qui rduisent l'paisseur du
ligament dans le plan de symtrie, crant ainsi un plan de rupture privilgi. Ensuite, la faible
paisseur des CT conduit une plastification plus tendue en pointe de fissure. Cette
plastification se dveloppe notamment pour les angles avoisinant 45 crant une zone plastique
en forme d' "oreille". Lorsque l'amorage est tardif, le niveau de chargement lev accentue
l'influence de cette zone plastique, o va se diriger la fissure dans les premiers instants de
propagation.
Compte tenu du chargement, un critre en contrainte circonfrentielle maximale
(traditionnellement utilis en mcanique linaire de la rupture) donne toujours le trajet
rectiligne comme direction privilgie. Or, on observe bien que si la rupture ne se produit pas
des niveaux de chargement faible, un branchement se produit pour des niveaux de chargement
plus levs. Il semble donc ncessaire de recourir un autre type de critres (bas sur la
contrainte de cisaillement par exemple) pour prdire l'ventuel branchement de la fissure. Nous
nous contentons dans ce mmoire de dcrire ce type de trajet sans y investir plus de temps,
mme son tude permettrait sans doute d'apporter un clairage nouveau sur la physique du
clivage.
Le second facteur influent sur la tnacit l'amorage est l'paisseur de l'prouvette. Pour
chaque trajet de fissure respectifs, les CT d'paisseur 10mm ont casses plus tt que celles de
5mm d'paisseur. Les modles bass sur une approche probabiliste (Master Curve [16] ou
Beremin [20]) prdisent une dispersion dans l'amorage fonction de l'paisseur de l'prouvette
la puissance comme indique par l'quation (128).
1/ 4
K = K min B ref / B K ref K min (128)
On utilise les rsultats prcdemment obtenus pour des CT d'paisseur 25mm comme
rfrences, et on applique la formule prcdente avec K min =20MPa m pour prdire une
valeur maximale de K aux paisseurs 10 et 5mm. En estimant K ref 130MPa m comme
tant la valeur maximale, on construit la courbe enveloppe trace sur la Fig. 81. On n'observe
pas de relle corrlation avec les rsultats exprimentaux. Pour placer les fortes valeurs de
tnacit l'amorage observes sur les CT d'paisseur 5mm sous la courbe enveloppe, il
115
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
faudrait vrifier une valeur maximale de K ref 210MPa m pour les CT d'paisseur 25mm.
x2 pointe de fissure
x tip x1
116
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
lisse les ventuelles oscillations, mais on risque alors de perdre des informations. Puisque l'on
cherche liminer les erreurs de mesure et extraire une tendance globale l'volution de la
vitesse dans le temps, une interpolation sur 5 points semble raisonnable.
(a)
(b)
117
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
(a)
(b)
118
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
R 15mm
60
5,5mm
25
52,5mm
47mm
25mm
60mm
110mm
(a) (b)
Fig. 85 : Plans des anneaux.
(a) Anneau pour chargement mcanique en mode I pur.
(b) Usinage supplmentaire pour chargement mcanique en mode mixte.
Afin de raliser des essais en mode mixte, on usine un couple d'encoches supplmentaire pour
permettre dincliner le plan de fissure de 25 avec lhorizontale lors de lessai rupture. Cette
gomtrie est dcrite sur la Fig. 85 (b).
119
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
520SI a0
520SH af
520UL
520UN
520UM
520UO
520US
520UT
520UQ
520UV
520UW } modemixte
0 2 4 6 8 10
longueur de fissure (mm)
(a) (b)
Fig. 87 : Courbe force-ouverture pour les anneaux en mode I pur (a), et pour les anneaux
chargs en mode mixte (b).
Bien qu'ayant des longueurs de fissure variables, les anneaux fournissent des courbes force-
120
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
ouverture relativement proches. On note que ces courbes sont loin d'tre linaires, ce qui
signifie qu'une zone importante de la structure a plastifi durant la mise en charge.
(a)
(b)
Fig. 88 : Longueur de fissure en fonction du temps pour les anneaux en mode I pur.
(a) Rsultats exprimentaux de l'ensemble des anneaux.
(b) Mise en vidence du dphasage temporel entre les mesures des deux faces d'un anneau.
On constate parfois l'existence d'un lger dcalage temporel entre les mesures ralises de part
et d'autre de l'prouvette (cf. Fig. 88 (b)). Ce dphasage s'explique par l'importante largeur des
anneaux par rapport l'avance totale de la fissure. Si l'amorage a lieu en un point du front
plus proche d'un bord que d'un autre, alors il existe un dcalage temporel entre les mesures
enregistres sur ces deux bord, correspondant la diffrence de temps ncessaire au dfaut pour
121
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
atteindre chacun des bords. D'ailleurs, ce dphasage temporel est infrieur au temps de
propagation de la fissure sur une distance gale l'paisseur des anneaux (en se basant sur les
vitesses mesures aux premiers instants), ce qui conforte cette hypothse. Bien que non
ngligeable, ce dphasage ne modifie pas l'volution de la vitesse de propagation avec la
longueur de fissure. Les profils obtenus restent trs proches.
(a)
(b)
122
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
L'volution de la longueur de fissure dans le temps peut donc tre reprsente de deux manires
diffrentes comme on le voit sur la Fig. 90.
(a) (b)
Fig. 90 : Longueur de fissure en fonction du temps dans le cas des anneaux en mode mixte.
(a) Dans le repre de la jauge. (b) Dans le repre curviligne suivant la pointe de fissure.
III.6.2. Observations
L'initiation et la propagation
La Fig. 91 est reprsentative du facis observ sur l'ensemble des gomtries o la vitesse de
propagation est importante. On constante un facis de rupture caractristique de la rupture par
clivage. Surtout, aucun ligament ductile n'est visible dans cette zone.
123
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
L'arrt
La Fig. 92 prsente le mme facis mais prs du front de fissure. On constate le mme type de
mcanisme, sans ligament ductile visible mais prsentant ici des marches de cisaillement. Ces
marches sont certainement cres la jonction de deux plans parallles plus faciles rompre
pour la fissure qui cette endroit dispose de nettement moins d'nergie que lors de la phase
prcdente.
crack arrest
front
124
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
Le branchement des CT
La Fig. 71 prsente le facis de rupture l'amorage d'une CT ayant branche. Il y apparat que
la fissure s'amorce directement avec un angle important. On distingue mme l'amorce de la
second fissure qui se fait dans la continuit de la premire. Ceci confirme pour ces prouvettes
minces le rle important de la plasticit lors de l'amorage.
crack
propagation
direction fatigue
crack
crack
propagation
fatigue direction
crack
(b)
(a)
(b)
(a)
III.7. Conclusions
Aprs une prsentation du matriau de l'tude et de son comportement mcanique (notamment
aux grandes vitesses de dformation), la campagne exprimentale de propagation et d'arrt de
fissure a t prsente au cours de ce chapitre. Deux gomtries ont t considres : il s'agit
d'prouvettes CT, et d'anneaux entaills en compression. Ces derniers ont galement t adapts
pour permettre un chargement en mode mixte.
125
Chapitre III. Essais de propagation et d'arrt de fissure sur un acier 16MND5
Pour chacune des configurations, les donnes relatives la prfissuration et la mise en charge
des prouvettes, ainsi que les informations concernant l'avance de la fissure dans le temps (et
ventuellement dans l'espace pour les anneaux en mode mixte) ont t prsentes en dtail.
Il s'avre que les prouvettes CT prsentent deux types de trajet de fissure trs diffrents. Soit la
fissure amorce de faibles niveaux de chargement et propage dans le plan de symtrie; soit
l'prouvette rompt plus tardivement et prsente une fissure branche ds l'initiation. Par contre,
les anneaux, plus larges, ont tous rompus dans le plan de symtrie lorsqu'ils ont t sollicits en
mode I. Le cas d'un chargement en mode mixte a galement t trait.
Les mesures de vitesse de propagation sont trs cohrentes entre elles. On constate que pour les
prouvettes CT, la propagation de fissure s'effectue une vitesse moyenne comprise entre 200
et 500m/s, alors qu'elle avoisine 700m/s pour les anneaux en mode I. De plus, la vitesse de
propagation semble augmenter avec la quantit d'nergie emmagasine par la structure et
disponible pour le clivage.
Les observations fractographiques ont montr des surfaces caractristiques du clivage, avec des
marches de cisaillement plus marques vers le front d'arrt. Aucun ligament ductile n'a t
rencontr. Ainsi, l'analyse de ces essais devra utiliser un critre de propagation et d'arrt de
fissure dcrivant uniquement les mcanismes de clivage.
Les donnes exprimentales rcoltes vont faire l'objet d'une analyse mcanique dtaille dans
le prochain chapitre. Elles constituent de ce fait la base physique de l'identification d'un modle
de propagation de fissure de clivage.
126
Chapitre IV.
Sommaire
IV.1. Introduction...........................................................................................
......................128
IV.1.1. Contexte..................................................................................................
...............128
IV.1.2. Dmarche.................................................................................................
..............128
IV.2. Interprtation des essais sur prouvettes CT...........................................................129
IV.2.1. Modle mcanique................................................................................................
.129
IV.2.2. Analyse nergtique.........................................................................
.....................132
IV.2.3. Analyse locale......................................................................................................
..133
IV.2.4. Sur la ncessit de prendre en compte les aspects dynamiques...........................140
IV.3. Analyses prdictives..................................................................................
..................140
IV.3.1. Schmas numriques associs aux modles de propagation................................140
IV.3.2. Analyses prdictives sur prouvettes CT..............................................................141
IV.3.3. Analyses prdictives sur prouvettes annulaires en mode I.................................143
IV.3.4. Analyses prdictives en mode mixte................................................................
.....145
IV.4. Conclusion......................................................................................................
..............155
Chapitre IV. Modlisation de la propagation de fissure de clivage
IV.1. Introduction
IV.1.1. Contexte
Les modles d'amorage font l'objet d'un plus grand nombre d'tudes que les modles de
propagation et d'arrt de fissure. Nanmoins, ces derniers sont ncessaires lorsque l'on envisage
de dmontrer qu'une fissure ayant amorc en fragile s'arrte avant de mettre en pril l'intgrit
de la structure.
Dans le chapitre bibliographique, nous avons vu que la plupart des analyses de propagation et
d'arrt de fissure fragile sont effectues dans un contexte lastodynamique. L'extension au cas
non linaire complexifie les analyses, notamment parce qu'il est difficile de dfinir un
paramtre unique reprsentatif de la sollicitation en pointe de fissure. La plasticit et les effets
de dchargement locaux lors de la propagation, invalide le sens physique de l'intgrale Jdyn. Les
modles locaux (comme le critre RKR [19]) semblent dans ce contexte plus attrayants. Ceux-
ci postulent en gnral l'existence d'une contrainte critique de clivage. Ce type de critre semble
priori tre facilement gnralisable au cas de la propagation et de l'arrt de fissure.
Ce chapitre prsente une analyse avance de la propagation et de l'arrt d'une fissure fragile, et
propose un modle de propagation, qui sera valid ou non par l'tude des diffrentes
configurations exprimentales prsentes au cours du chapitre prcdent.
IV.1.2. Dmarche
L'identification du modle de propagation s'appuie sur deux macros-tapes :
L'identification du modle en ralisant des simulations numriques o la vitesse de
propagation est impose ("generation phase").
L'application de ce modle pour prdire par la simulation l'avance de la fissure dans le
temps ("application phase").
Cette dmarche propose par Kanninen [51] a galement t utilise par Nishioka et al [36].
Nous proposons d'effectuer l'identification du modle de propagation partir des expriences
sur prouvettes CT, et d'appliquer ensuite ce modle ces mmes prouvettes CT, ainsi qu'aux
prouvettes annulaires en mode I et en mode mixte.
La pertinence du modle est value en comparant aux rsultats exprimentaux, la vitesse de
propagation et la longueur de fissure l'arrt obtenues par application du modle (calculs
prdictifs).
Lors de la phase d'application du modle, l'utilisation d'une gomtrie diffrente de celle qui a
servi l'identification est un indicateur quant la transfrabilit du modle. Pour qu'un critre
ait un sens physique, il est ncessaire qu'il ne dpende pas de la gomtrie tudie, ni d'une
variable de configuration (comme la longueur de prfissuration, ou la mixit des modes par
exemple). Cette dmarche donne du poids la dmonstration si la confrontation avec
l'exprience des rsultats de simulation est rellement prdictive.
La premire partie de ce chapitre (paragraphe IV.2) correspond l'interprtation des essais sur
prouvette CT. Il expose d'abord le modle mcanique utilis dans l'ensemble des simulations,
avant de prsenter les rsultats des simulations et de proposer une identification d'un modle de
propagation.
La seconde partie (paragraphe IV.3) concerne l'application du modle de propagation identifi.
Des simulations prdictives de propagation de fissure en mode I sur prouvettes CT, en mode I
sur prouvettes annulaires, et en mode mixte sur prouvettes annulaires y sont dtailles.
128
Chapitre IV. Modlisation de la propagation de fissure de clivage
u
F
+ tip
quation d'quilibre
En l'absence de force volumique, l'quilibre dynamique, les conditions initiales et aux limites
sont dfinis par le jeu d'quations (132).
quation d'quilibre : div = u
{
n = F ext sur F
conditions aux limites : n = 0 sur + -
(132)
u = uext sur u
{
0
u x ,t =0 = u
conditions initiales :
u x ,t =0 = u0
129
Chapitre IV. Modlisation de la propagation de fissure de clivage
Les trois premires lignes du jeu d'quations (132) peuvent se mettre sous la forme
variationnelle (133) (formulation faible).
vV ,
0
:v d uv d = ext
F v dS
F
(133)
0
avec V = { v , v=0 sur u }
quations constitutives
Le comportement du matriau est lasto-viscoplastique crouissage isotrope, avec une loi
d'coulement associe. La partition de la dformation en une partie lastique et une partie
inlastique est rappele par les quations (134).
Partition de la dformation : = e in
{
0 si , p , p 0 (134)
e 1 in
avec =D et =
p si , p , p =0
in in
o D est la matrice de Hooke, p= 23 : reprsente la vitesse de dformation inlastique
quivalente, et est la fonction seuil dfinie par (135).
, p , p = eq y p , p (135)
La surface de charge utilise la contrainte de Von Mises (136).
eq = sij = ij 13 kk ij
3
2
s:s avec le dviateur (136)
On remarque que la viscosit est prise en compte dans la contrainte seuil y qui dpend de la
vitesse de dformation inlastique quivalente p . La loi y p , p est en fait la courbe
d'crouissage issue de l'identification exprimentale du comportement du matriau, qu'on a
prsent au chapitre III et dans l'annexe B. Rappelons qu'il s'agit d'une loi de Symonds-Cowper
modifie pour tenir compte de l'adoucissement d l'chauffement cr par la dissipation
plastique (cet chauffement tant considr comme adiabatique, la temprature T n'est pas
indpendante et peut tre directement relie p et p ).
Modle de propagation
Pour tre complet dans notre description du problme, il faut ajouter les inconnues dcrivant la
l'volution de la fissure dans le temps. Pour un problme bidimensionnel en mode I pur, il s'agit
simplement de a t . Dans un cadre plus gnral (tridimensionnel en mode mixte), il faut
dcrire l'volution du front de fissure tip t (cf. Fig. 94).
On considre deux types de modlisations :
Les modlisations vitesse de propagation impose, o tip t = imp t .
Les modlisations prdictives, o tip t =F ut , t ,...
La dmarche adopte pour identifier le modle de propagation s'appuie sur ces deux types de
modlisations.
Avec la mthode X-FEM, le mouvement du front de fissure est dcrit naturellement par
l'volution des fonctions de niveau dans le temps. En chaque point du front, il faut tre capable
130
Chapitre IV. Modlisation de la propagation de fissure de clivage
de dterminer si la fissure propage, et dans le cas affirmatif, dans quelle direction et quelle
vitesse, comme cela est schmatis sur la Fig. 95.
x2
crack a
x1
Maillage
Le maillage utilis pour l'tude de l'prouvette CT est donn sur la Fig. 96. Il compte environ
4200 lments (soit 4350 noeuds) avec une taille de maille en pointe de fissure de 200m. Les
lments dans la zone de propagation de fissure sont uniquement des quadrangles 4 noeuds
(QUA4). Les fonctions d'interpolation gomtriques sont donc du 1er ordre.
131
Chapitre IV. Modlisation de la propagation de fissure de clivage
(a) (b)
Fig. 97 : Comparaison entre la modlisation et l'exprience de la courbe force-ouverture
pendant la mise en charge (a), et imposition des donnes exprimentales concernant
l'volution de longueur de fissure avec le temps (b).
On procde d'abord une tape de mise en charge quasi-statique reprsente sur la Fig. 97(a),
et ralise jusqu' l'amorage dtermin partir de la tnacit exprimentale. Ensuite, la
propagation de la fissure s'effectue en dynamique avec la loi de comportement lasto-
viscoplastique prcdemment dcrite.
132
Chapitre IV. Modlisation de la propagation de fissure de clivage
reste relativement faible dans l'prouvette. L'nergie cintique augmente pendant la moiti de la
propagation avant de diminuer. Des tendances similaires ont dj t observes par Hahn et al
[45] sur une Double Cantilever Beam en acier et par Kobayashi [48] sur une prouvette CT en
polymre.
133
Chapitre IV. Modlisation de la propagation de fissure de clivage
matriau. Nous pensons qu'il s'agit probablement d'un lment essentiel pour la comprhension
des phnomnes (cf. Kanninen [51]). Parmi les travaux les plus rcents sur ce sujet, Hajjaj [56]
suit une dmarche similaire pour simuler la propagation de fissure d'un anneau sous choc
thermique. Il modlise la propagation et larrt de fissure en considrant un comportement
lasto-viscoplastique du matriau tudi, un acier de cuve REP. Lauteur s'appuie dans son
analyse sur un critre RKR avec une contrainte critique de clivage dpendant uniquement de la
temprature, partir dune relation dcoulant de celle entre la tnacit lamorage et la
temprature.
134
Chapitre IV. Modlisation de la propagation de fissure de clivage
C'est le calcul de l'angle de propagation chez ces mmes auteurs qui a inspir la troisime
mthode de calcul de la contrainte principale.
Mariani et al [147] cherche par exemple une approximation polynomiale en r du tenseur des
contraintes sur un demi-disque en avant de la pointe de fissure. Puis, il utilise la direction
perpendiculaire la contrainte principale maximale comme direction de propagation.
Wells et al [60] et Remmers et al [148] moyennent le tenseur des contraintes sur un disque en
pointe de fissure de rayon gal trois fois la taille d'un lment via une fonction de pondration
Gaussienne w dfinie par l'quation (138). Ils dcident galement de fixer la direction d'avance
de la fissure comme tant la direction perpendiculaire la contrainte principale maximale de ce
tenseur moyen.
w=
1
2 3/ 2 l 3
exp
r 2
2l 2
(138)
L'ide de ces mthodes n'est pas de dlocaliser le tenseur des contraintes (comme on le fait pour
les modles d'endommagement par exemple), mais de s'affranchir de possibles erreurs
numriques locales en intgrant les contraintes sur un petit volume. Rappelons que le modle de
Beremin [20] calcule une probabilit de rupture en intgrant sur le volume plastique en pointe
de fissure, une fonction dpendant de la contrainte principale maximale.
Au lieu de prendre un disque complet de rayon gal trois fois la taille d'un lment, on choisit
de fixer la dimension l sur laquelle est effectue la moyenne du tenseur des contraintes. Cette
technique introduit donc la notion d'une distance caractristique comme le critre RKR. Cette
distance doit tre suffisamment faible pour caractriser l'tat en pointe de fissure, mais
galement assez leve pour dfinir une zone de dimension raisonnable par rapport la taille
d'un lment. Nous la fixons 200m. Comme une fonction Gaussienne est utilise pour
pondrer cette moyenne, plus les points d'intgration seront proches de la pointe de fissure et
plus ils auront un poids important. Mis part cette proprit, le choix de la forme de cette
fonction de pondration demeure arbitraire. On dcide ne pas prendre en compte les points pour
lesquels w0.1 , et de normaliser la contrainte par l'intgrale de la fonction poids seule.
D'autre part comme l'tat de contrainte en arrire de la pointe de fissure est faible et ne joue pas
de rle prdominant dans le processus de rupture, on ne moyenne les contraintes que sur un
demi disque en avant du front de fissure. Cette mthode de calcul de la contrainte moyenne sur
un demi disque est rsume par le jeu d'quations (139).
{ {10
2
w 0 r H si w 0 0.1 r si 0
avec w= , w 0 =exp et H =
0 si w 0 0.1 2l 2 si 0
car dsigne la fonction de niveau associe au front de fissure, et [ ... ] I reprsente la
contrainte principale maximale en traction.
135
Chapitre IV. Modlisation de la propagation de fissure de clivage
136
Chapitre IV. Modlisation de la propagation de fissure de clivage
137
Chapitre IV. Modlisation de la propagation de fissure de clivage
A
Ic a= m (140)
1 a/ vlim
Par rapport la formule gnrale ci-dessus, on a fix les paramtres m=1 et v lim=c R . Ainsi,
on obtient respectivement pour les critres de type RKR100m et HD200m les grandeurs
suivantes : A =1500MPa et A =1150MPa .
La dpendance observe de la contrainte critique avec la vitesse de propagation indique
certainement une sensibilit de cette contrainte vis vis de la vitesse de dformation associe.
En effet, d'un point de vue mcanique, c'est la vitesse de propagation qui entrane la rapidit des
phnomnes en pointe de fissure. D'un point de vue plus phnomnologique, la rupture est une
rponse du matriau des sollicitations qui se caractrisent par les champs de contraintes,
dformation, vitesse de dformation, tempratures, etc... En fait, il s'agit d'un problme coupl,
dans lequel la propagation dynamique est la fois la cause et la consquence des champs
mcaniques en pointe de fissure. Aussi, il apparat que la dpendance de la contrainte critique
de clivage la rapidit des phnomnes en pointe de fissure, peut tre traduite aussi bien par la
vitesse de propagation que par la vitesse de dformation, l'une et l'autre dcrivant la mme
chose mais d'un point de vue diffrent.
Pour le vrifier, on reprsente sur la Fig. 103 lvolution de la contrainte calcule avec le critre
RKR100m avec la vitesse de la dformation totale quivalente calcule au mme point.
Une dpendance de la contrainte critique avec la vitesse de dformation nest pas un fait
actuellement reconnu. Une contrainte critique constante est trs souvent suppose dans la
littrature. Hajjaj [56] a par exemple considr une unique dpendance la temprature de la
contrainte critique. Nanmoins, les essais de choc thermique qu'il a modlis prsente un
gradient de temprature de plus de 400C. La variation de la contrainte due la temprature est
si importante qu'elle peut masquer la dpendance avec la vitesse de dformation prsente ici.
138
Chapitre IV. Modlisation de la propagation de fissure de clivage
Un rsultat identique est obtenu avec la seconde mthode de calcul de la contrainte sur un demi
disque (cf. Fig. 104 et quation (142)).
139
Chapitre IV. Modlisation de la propagation de fissure de clivage
Fig. 106 : Comparaison entre une analyse quasi-statique et une analyse dynamique
concernant l'volution de la contrainte pour la ct5pm.
On s'aperoit qu'une analyse statique conduit des valeurs de contrainte variables au cours de la
propagation, et ne permet donc pas d'identifier un critre de propagation et d'arrt de fissure.
140
Chapitre IV. Modlisation de la propagation de fissure de clivage
B n1 , n1
i i
Application du critre : (145)
n1 = t n1 [ 1 B ]
i
Correction de la vitesse de propagation avec par ex. : t i1 i
141
Chapitre IV. Modlisation de la propagation de fissure de clivage
(a)
(b)
Fig. 107 : Longueur de fissure en fonction du temps pour la ct5pm (a), et la ct10oy (b).
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de propagation.
142
Chapitre IV. Modlisation de la propagation de fissure de clivage
(a)
(b)
Fig. 108 : Longueur de fissure en fonction du temps pour la ct5pm (a), et la ct10oy (b).
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de dformation.
143
Chapitre IV. Modlisation de la propagation de fissure de clivage
Fig. 109 : Longueur de fissure en fonction du temps pour l'anneau 520uo. Comparaison avec
l'exprience des rsultats de simulation utilisant une contrainte critique dpendant de la
vitesse de propagation.
Fig. 110 : Longueur de fissure en fonction du temps pour l'anneau 520uo. Comparaison avec
l'exprience des rsultats de simulation utilisant une contrainte critique dpendant de la
vitesse de dformation.
Le transfert des critres au cas des anneaux en mode I pur s'avre trs satisfaisant, car la
prdiction de la vitesse de propagation et de la longueur de fissure larrt est trs prcise. Ceci
est dautant plus remarquable, que les vitesses de propagation et les longueurs de fissure sont
bien diffrentes de celles considres pour les prouvettes CT : jusqu' 1000 m/s sur les
anneaux comparer au maximum de 600 m/s observes sur les prouvettes CT. on note
nanmoins une vitesse un peu moindre lorsque le critre utilisant une contrainte critique
dpendant de la vitesse de dformation est utilis.
144
Chapitre IV. Modlisation de la propagation de fissure de clivage
[ ] (146)
2
145
Chapitre IV. Modlisation de la propagation de fissure de clivage
Fig. 111 : volution angulaire des contraintes d'ouverture et principale interpoles sur un
cercle de rayon de 100m obtenue pour l'anneau 520uv.
Dans un second temps, nous avons ainsi prfr utiliser le critre not c I , bas sur la
direction perpendiculaire la contrainte principale maximale du tenseur des contraintes
moyenn. En effet, cette mthode donne un trajet de fissure de bonne qualit, et limine les
oscillations dues des problmes numriques. Les dtails relatifs au calcul de ce tenseur
moyenn ont dj t prsents par l'quation (139) au paragraphe IV.2.3.
Rsultats prdictifs
Les rsultats des simulations prdictives en terme de trajet de fissure et d'volution dans le
temps sont respectivement fournis par les Fig. 113 et Fig. 114 pour les critres bass sur une
contrainte critique dpendant de la vitesse de propagation, et par les Fig. 115 et Fig. 116 pour
les critres bass sur une contrainte critique dpendant de la vitesse de dformation.
146
Chapitre IV. Modlisation de la propagation de fissure de clivage
Fig. 112 : Trajet de fissure pour l'anneau en mode mixte 520uv. Comparaison avec
l'exprience des rsultats de simulation obtenus en utilisant comme direction de propagation
celle du maximum de la contrainte circonfrentielle.
Fig. 113 : Trajet de fissure pour l'anneau en mode mixte 520uv. Comparaison avec
l'exprience des rsultats de simulation utilisant une contrainte critique dpendant de la
vitesse de propagation.
147
Chapitre IV. Modlisation de la propagation de fissure de clivage
Fig. 114 : Longueur de fissure en fonction du temps pour l'anneau en mode mixte 520uv.
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de propagation.
Fig. 115 : Trajet de fissure pour l'anneau en mode mixte 520uv. Comparaison avec
l'exprience des rsultats de simulation utilisant une contrainte critique dpendant de la
vitesse de dformation.
148
Chapitre IV. Modlisation de la propagation de fissure de clivage
Fig. 116 : Longueur de fissure en fonction du temps pour l'anneau en mode mixte 520uv.
Comparaison avec l'exprience des rsultats de simulation utilisant une contrainte critique
dpendant de la vitesse de dformation.
On observe un bon accord entre les rsultats exprimentaux et les simulations par lments finis
pour les deux types de critre. Le trajet de fissure prdit se superpose aux points exprimentaux
relevs. La vitesse de propagation est un peu moins bien estime par le critre utilisant une
contrainte critique fonction de la vitesse de dformation, mais cela est probablement la
consquence d'un effet tridimensionnel. En effet, l'paisseur des anneaux diffre de celle des
prouvettes CT sur lesquels a t identifi le modle de propagation. Les deux critres
surestiment lgrement la longueur de fissure l'arrt, mais les volutions prdites restent de
bonne qualit.
149
Chapitre IV. Modlisation de la propagation de fissure de clivage
9s
19s
150
Chapitre IV. Modlisation de la propagation de fissure de clivage
37s
87s
151
Chapitre IV. Modlisation de la propagation de fissure de clivage
9s
19s
152
Chapitre IV. Modlisation de la propagation de fissure de clivage
37s
87s
153
Chapitre IV. Modlisation de la propagation de fissure de clivage
Fig. 119 : Comparaison avec l'exprience des trajets de fissure obtenus par la simulation pour
les deux anneaux chargs en mode mixte 520uv et 520uw.
On constate que le modle de propagation permet de reproduire correctement la diffrence de
trajet de fissure qu'il existe entre les deux anneaux chargs en mode mixte, mais galement que
l'arrt tardif observ sur l'anneau 520uv est encore plus flagrant sur l'anneau 520uw. Cela peut
s'expliquer en considrant la Fig. 120.
154
Chapitre IV. Modlisation de la propagation de fissure de clivage
rapport au cas de l'anneau 520uv, ou mme des autres anneaux et des prouvettes CT.
Afin de confirmer cette hypothse, on ralise des simulations identiques, mais en imposant cette
fois ci un amorage lorsque le produit force-ouverture de la simulation atteint le produit
exprimental, esprant ainsi obtenir une rupture nergie de dformation lastique quivalente.
L'instant t d'amorage est ainsi dfini par l'galit (137).
F t t simulation = F t Rupture t Rupture experience (148)
Les rsultats des simulations obtenues avec cet amorage sont prsents sur la Fig. 105.
Fig. 121 : Comparaison avec l'exprience des trajets de fissure obtenus par la simulation pour
les deux anneaux chargs en mode mixte 520uv et 520uw.
On constate que les longueurs l'arrt prdites sont nettement plus proches des rsultats
exprimentaux, confortant ainsi l'explication avance sur les prcdents rsultats, mais le trajet
de fissure s'loigne lgrement de celui relev exprimentalement sur le bord des prouvettes.
Finalement, il s'avre qu'une modlisation tridimensionnel du problme serait souhaitable. Cela
permettrait notamment de reproduire plus fidlement la courbe force-ouverture exprimentale,
et d'viter ainsi les mises au point d la dfinition de l'instant d'amorage. De plus, il est fort
probable qu'une tude tridimensionnelle favoriserait la description fine des nuances de trajet de
fissure que nous avons prsentes dans ce paragraphe.
En outre, on peut vraisemblablement penser que les conditions aux limites jouent un rle
dissipatif (frottement du cylindre d'appui sur l'anneau) dans le mouvement d'ensemble de la
structure, ce qui limiterait l'ouverture de la fissure. De manire gnrale, des donnes
exprimentales plus pointues concernant les conditions aux limites et le chargement, sont
souhaitables afin d'affiner les modlisations.
IV.4. Conclusion
Afin de prdire le comportement d'une structure fissure, il est ncessaire d'introduire un
modle de propagation dcrivant l'volution de la fissure en fonction des paramtres
mcaniques.
Dans une premire tape, on a impos aux prouvettes CT l'avance de la fissure en fonction du
temps mesure exprimentalement. Deux manires d'valuer la contrainte principale en pointe
de fissure ont t mises en avant : la mthode ponctuelle de Ritchie-Knott-Rice [19], et une
155
Chapitre IV. Modlisation de la propagation de fissure de clivage
mthode qui moyenne le tenseur des contraintes [60]. Pour chacune de ces mthodes, nous
avons montr qu'il existait une dpendance de la contrainte critique de clivage avec la vitesse
des phnomnes en jeu, qui peut tre dcrite soit par la vitesse de propagation de la fissure
directement, soit par l'valuation de la vitesse de dformation.
Les critres ainsi identifis ont ensuite t prouvs lors de simulations prdictives, c'est dire
en considrant la position du front de fissure comme une inconnue supplmentaire au problme
mcanique. Les simulations prdictives sur prouvettes CT donnent de bons rsultats, mais il
s'agit de la mme gomtrie qui a servit l'identification. Ainsi, afin d'assurer une certaine
transfrabilit des critres, des simulations des essais sur anneaux en compression ont t
menes. Les vitesses de propagation et les longueurs de fissures l'arrt obtenues par la
simulation sont en bon accord avec les rsultats exprimentaux. La gnralisation des modles
de propagation proposs au cas d'une sollicitation en mode mixte, a permis de mettre une
nouvelle fois en vidence de bons rsultats prdictifs. Les vitesses mais galement les trajets de
fissure exprimentaux ont pu tre retrouvs par la simulation. Une surestimation de la longueur
de fissure est nanmoins observe dans les cas o l'amorage est tardif et introduit dans les
simulations bidimensionnelles une quantit d'nergie lastique suprieure l'exprience. Cette
lgre erreur pourrait sans doute tre gomme en adoptant une modlisation tridimensionnelle
du problme.
156
Conclusion gnrale et perspectives
Dans un premier temps, la mthode des lments finis tendus (X-FEM) a t implante dans le
code de calcul Cast3m. Cet outil numrique permet de simuler efficacement la propagation
arbitraire de fissure, puisque l'on vite les oprations de maillage et de remaillage.
Afin de garantir un bon comportement de la mthode dans le cas inlastique, une mthode
d'intgration non conforme des lments proches du front de fissure est prsente. Elle permet
d'intgrer de manire prcise les fonctions d'enrichissements, et vite la projection de champs
au voisinage de la pointe de fissure lors de la propagation.
Une seconde amlioration prsente dans ce mmoire concerne l'utilisation des fonctions de
niveau (Level Set function). Ces dernires permettent de reprsenter implicitement la position
et la gomtrie de la fissure. Mais leur actualisation souffre d'imprcisions numriques
lorsqu'elle est ralise sur un maillage irrgulier, ou si la taille de maille varie rapidement. Afin
d'amliorer cette technique, on propose d'introduire une grille auxiliaire et rgulire, sur
laquelle on dfinit et actualise les fonctions de niveau. Ainsi, on dcouple compltement le
maillage de la structure, du support dfinissant les fonctions de niveau et donc de la gomtrie
de la fissure. Cette technique permet la fois d'augmenter si ncessaire le raffinement de la
description de la fissure, et d'utiliser un schma aux diffrences finies, simple et efficace, pour
actualiser les fonctions de niveau. Ce point est important car la mthode des fonctions de niveau
a originalement t dveloppe dans un contexte de mthodes aux diffrences finies.
La mthode des lment finis tendus et les deux techniques particulires rappeles ci-dessus
ont fait l'objet de plusieurs tests dvelopps dans le chapitre II de ce mmoire. Les bons
rsultats rapports permettent d'envisager la finalisation de leur implantation dans le code de
calcul Cast3m.
Dans un second temps, une campagne d'essais s'appuyant sur un acier 16MND5 utilis pour les
cuves de racteurs nuclaires REP, a t mene. Une loi de comportement lasto-viscoplastique
base sur le modle de Symonds-Cowper est propose afin de reproduire la rponse de ce
matriau aux grandes vitesses de dformation.
Les essais de rupture fragile concernent deux gomtries : des prouvettes Compact Tension, et
des prouvettes annulaires en compression. Ces dernires sont galement adaptes pour traiter
le cas du mode mixte. La mesures de la vitesse de propagation est obtenue en utilisant des
jauges brins colles sur la surface des prouvettes. Il s'avre que la propagation de fissure
s'effectue sur les prouvettes CT une vitesse comprise entre 200 et 500m/s, tandis qu'elle
atteint environ 700m/s sur les prouvettes annulaires sollicites en mode I. On constate que
cette vitesse semble augmenter avec la quantit d'nergie disponible dans la structure au
moment de l'amorage.
Pour des niveaux d'amorage trs levs, un branchement de la fissure est mme observ sur les
CT : la fissure se spare ds l'initiation en deux branches voluant hors du plan de symtrie de
l'prouvette. L'paisseur rduite et par consquent la plasticit, peut probablement tre un
facteur retardant l'amorage et conduisant ce type de trajet particulier.
Les analyses fractographiques entreprises sur l'ensemble des prouvettes montrent clairement
que le mcanisme l'origine de la ruine du matriau est le clivage. Des marches de cisaillement
plus marques proximit du front d'arrt semblent traduire la faible quantit d'nergie alors
disponible pour sparer les plans cristallographiques. En effet, vers la fin de la propagation, la
fissure n'a plus d'autre choix que de se diriger l o les plans cristallins sont les plus favorables.
Par ailleurs, contrairement des tudes prcdentes, aucun ligament ductile n'est relev. Ainsi,
un critre de propagation et d'arrt de fissure pertinent doit uniquement dcrire les mcanismes
de clivage.
Les donnes exprimentales sont le support de l'analyse mcanique entreprise dans le chapitre
IV. On adopte une dmarche en deux temps. D'abord, une modlisation par lments finis des
essais sur prouvettes CT est ralise en imposant de manire forte la vitesse de propagation de
la fissure mesure exprimentalement. Pendant cette tape, certaines grandeurs caractristiques
sont releves et discutes. Ensuite, on construit un modle de propagation bas sur
l'identification prcdente, et on l'applique l'ensemble des prouvettes : CT, annulaires en
mode I et en mode mixte.
Afin de construire un modle de propagation robuste et direct dans son interprtation, on oriente
la discussion vers un critre local, bas sur une "mesure" de la contrainte principale en pointe de
fissure. L'tape d'identification sur CT montre que quelque soit la dfinition de la contrainte
principale, on observe une dpendance de la contrainte critique de clivage avec la vitesse de
sollicitation. Cette vitesse de sollicitation peut tre quantifie par la vitesse de propagation, ou
par la vitesse dformation en pointe de fissure.
L'application de cette classe de modles aux prouvettes CT, permet de retrouver naturellement
l'volution du front de fissure en fonction du temps. L'extension aux anneaux en mode I fournit
galement des rsultats trs satisfaisants. Ceci est d'autant plus remarquable que les vitesses de
propagation des anneaux sont suprieures celles des CT, et que le modle identifi extrapole
la contrainte critique ces vitesses. Pour traiter du mode mixte, on montre qu'un critre local
bas sur la contrainte d'ouverture maximale, ou sur la contrainte principale du tenseur moyenn
en pointe de fissure permet de retrouver de manire extrmement prcise le trajet de la fissure
dans l'espace. Les vitesses de propagation restent encore satisfaisantes.
158
volutions relatives de la contrainte critique de clivage et de la limite d'lasticit. L'existence
d'un lien entre ces deux grandeurs permettrait une meilleure comprhension de la physique des
mcanismes du clivage, et une extension immdiate du modle de rupture aux tempratures.
Cette extension du modle reste une voie explorer afin de pouvoir traiter le cas du choc
thermique sur anneau.
Enfin, le dveloppement du modle de propagation par clivage doit aussi tre capable de
prdire les fissurations plus surprenantes, comme le branchement observ sur les prouvettes
CT paisseur rduite.
159
Chapitre V. Annexes
Sommaire
Syntaxe
La syntaxe de l'oprateur HPROPA est donne par :
cas 1 : on ralise l'actualisation d'une fonction de niveau.
CHPO2 = HPRO CHPO1 | CHPO3 FLOT1 (MOT1) (MOT2) |
| MOT3 CHPO4
cas2 : on calcule une grandeur associe une fonction de niveau.
CHPO2 = HPRO CHPO1 MOT4
Commentaire
Il est ncessaire que le maillage support de tous les champs par points fournis dans cet
oprateur soit le mme rectangle rgulier orient selon le repre global.
Nom Description Grandeur
CHPO1 Fonction de niveau l'instant t. t
cas 1
CHPO3 Vitesse de l'interface, scalaire (terme V ou V
d'expansion) ou vectorielle (terme d'advection).
FLOT1 Pas de temps. t
MOT1 = "NOREINIT". vite l'opration de
rinitialisation.
Annexe A : Implantation de la mthode des lments finis tendus (X-FEM) dans Cast3M
164
Annexe A : Implantation de la mthode des lments finis tendus (X-FEM) dans Cast3M
4 3
17 20 21 24 etc ...
1 2
14 58 9 12 13 16
(a) (b)
Fig. 122 : Numrotation des points d'intgration pour l'lment standard QUA4 (a), et pour
les lments enrichis QU4H ou QU4F (b).
La plupart des lments ncessaires pour l'introduction d'lments enrichis sont prsents dans
Cast3m. Cependant, deux points sont traiter avec prudence.
Il faut vrifier que chaque procdure utilise la bonne variable (par exemple, il faut parfois
boucler sur les fonctions de forme et pas sur les noeuds).
L'enrichissement utilise des fonctions dfinies localement. Il faut donc disposer de cette
information lorsque l'on en a besoin au cours du calcul.
Dans reshpx.eso (jumeau de reshpt.eso pour les lments enrichis), on n'introduit pas les
fonctions enrichis car celles-ci varient selon l'lment considr, mais les fonctions de forme
standards. On cre dans les oprateurs RIGI, MASS, SIGM, EPSI, BSIGM, des branchements
165
Annexe A : Implantation de la mthode des lments finis tendus (X-FEM) dans Cast3M
vers les fichiers rigix.eso, massex.eso, sigmax.eso, epix.eso, bsigmx.eso, qui eux-mmes utilise
shapx.eso (jumeau de shape.eso) lorsque le type d'lment correspond un lment enrichi.
Cette stratgie a t adopte afin de ne pas modifier trop lourdement les fichiers en place
utiliss pour les calculs avec des lments standards. Toutes ces procdures reposent sur une
bonne gestion de l'enrichissement durant la propagation.
Oprateur TRIELE
cas 1
CHAMX2 MODX2 = TRIE MODS1 CHPO1 CHPO2 (FLOT1)
cas 2
CHAMX2 MODX2 (CHAM21 CHAM22 ...) = TRIE CHAMX1 MODX1 CHPO1 CHPO2
(CHAM11 CHAM12 ...) (FLOT1)
Nom Description
CHPO1 Fonction de niveau reprsentant le front de fissure l'instant n
CHPO2 Fonction de niveau reprsentant le plan de fissure l'instant n
FLOT1 Permet de tester la proximit d'un noeud la pointe de fissure. Les
lments dont un des noeuds est proche de la pointe de fissure sont
intgrs comme des lments enrichis (sous-dcoupage non conforme).
CHAMX2 Champ par lment contenant l'information sur tous les
enrichissements depuis le dbut du calcul jusqu' l'instant n.
MODX2 Modle avec trois zones correspondantes aux trois types d'lments
dfinis l'instant n.
cas 1
MOD1S Modle standard (QUA4).
cas 2
CHAMX1 Champ par lment contenant l'information sur tous les
enrichissements depuis le dbut du calcul jusqu' l'instant n-1.
MODX1 Modle avec trois zones correspondantes aux trois types d'lments
dfinis l'instant n-1.
CHAM11, CHAM12, Champs par lment dfinis sur le modle MODX1.
...
CHAM21, CHAM22, Champs par lment dfinis sur le modle MODX2, issus de la
... projection des champs CHAM11, CHAM12, ...
166
Annexe A : Implantation de la mthode des lments finis tendus (X-FEM) dans Cast3M
167
Annexe B :
i t t t
V
specimen
2 u E 2 u
t 2 x2
=0 conduit {
V i = C i t r t
V t = C t t
(149)
puis { F i = S E i t r t
F t = S E t t
Rsultats
La vitesse de dformation varie lgrement durant l'essai. C'est pourquoi la dfinition d'une
seule vitesse de dformation est conventionnelle.
La temprature initiale de l'essai est obtenue en refroidissant dans un premier temps
l'chantillon par injection d'azote liquide, et en le laissant se rchauffer dans un second temps.
La dure du choc est suffisamment courte pour considrer le chargement comme tant
adiabatique. On peut considrer que toute l'nergie dissipe plastiquement est convertie en
chaleur. En prenant une capacit calorifique de 3,53E6 J/C/m3 [17], on obtient une
augmentation de la temprature non ngligeable pendant la dure de l'essai (entre 20 et 30C
10% de dformation). En procdant des essais interrompus, il a t montr que ce
rchauffement conduit une baisse significative de l'crouissage (cf. [143]).
Les effets de la temprature initiale et de la vitesse de chargement sont illustrs sur la Fig. 124.
(a)
(b)
170
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
vitesse de sollicitation augmente. Il s'agit d'une tendance gnrale, et on observe une variabilit
assez marque dans ces essais. Nanmoins, on peut affirmer que, dans notre intervalle d'tude,
la sensibilit de la contrainte avec la temprature est nettement plus marque que la variation de
la contrainte avec la vitesse de dformation.
=
[ A B in n C D in m ln
in
0
]
E in k 1T
(150)
La seconde ligne de l'quation (150) montre que la mme dpendance linaire avec la
temprature est suppose pour toutes les composantes de la contrainte. La contrainte
d'activation thermique est prpondrante pour les faibles vitesses de dformations ( 1s1 ),
alors que pour les vitesses leves, c'est le terme en puissance qui prvaut. Ce modle permet
de prendre en compte un crouissage viscoplastique diffrent du cas statique, et est capable de
donner de trs bons rsultats (cf. [154]).
La loi de Johnson-Cook [155] est une loi multiplicative dveloppe initialement pour des essais
balistiques. L'influence de chaque variable est clairement visible dans la relation (151).
in , in ,T = [ A B n ] [ 1C ln in ] 1
[ T T room
T melt T room ] (151)
[ ]
in
peut se r-crire in , in ,T = stat in ,T 1C ln stat
On remarque qu' priori, seul le paramtre C doit tre identifi puisque l'on connat le
comportement quasi-statique du matriau tudi. Nanmoins, stat peut galement tre
interprt comme un seuil partir duquel les effets de viscosit ne sont plus ngligeable. Dans
ce cas, il ne reprsente plus la vitesse de dformation statique et il doit galement tre identifi.
La loi de Symonds-Cowper [54] est une loi multiplicative puissance 2 paramtres (152).
[ ]
in 1
in in stat in p
(152)
, ,T = ,T 1
D
Cette loi a t utilise par Rossol [53] pour identifier le comportement de l'acier 16MND5. Les
coefficients issus de cette identification sont : D=1.E8 ; p=12.
Ce modle a ensuite t repris directement par Hajjaj [56]. La comparaison entre le modle et
l'exprience est illustre par la Fig. 125.
Hajjaj explique que l'cart entre le modle et exprience est d l'adoucissement thermique.
Cependant, cet effet devrait pouvoir tre pris en compte dans la modlisation afin d'tre
171
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
reprsentatif du comportement rel du matriau. Cela suppose donc soit une intgration
explicite du champ de temprature et de l'chauffement d la dformation plastique, soit une
intgration indirecte de cet chauffement via la loi de comportement. Cette seconde voie est
explore dans le prochain paragraphe.
ln
stat in
,0 ,T
VP in , in ,T 1 in 1
1 = ln ln D
p p
(153)
172
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
[ ]
1
,0 ,T 1H
in in stat in in p
, ,T =
o { p=4.250.015T
in in
H T , = abexp / c (154)
{
a=0.07592 0.07182T / 273.
et b=0.06352 0.2397T / 273.
c=0.02528 0.04124T / 273.
173
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
174
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
in int int t
Dfinition explicite de la vitesse de dformation : =
t
(155)
in int t int
Dfinition implicite de la vitesse de dformation : =
t
On distingue 3 principaux type d'approches pour intgrer une loi de comportement de ce type.
L'approche explicite.
L'approche implicite itrative.
L'approche implicite directe.
L'approche explicite
On suppose que la vitesse de dformation varie suffisamment lentement pour permettre
d'utiliser la valeur du pas prcdent. Cette mthode est inadapte au cas de la propagation d'une
fissure car les dformations inlastiques (et donc leur variations) sont localises en pointe de
175
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
fissure, laquelle n'occupe pas la mme place entre deux pas de temps successifs.
F 2 ip
F 2ip
D
B 2 ip1
F sol.
F 2ip 1
C
u
du 2ip 1
du 2ip
Fig. 127 : Illustration du phnomne d'oscillations numriques entre deux valeurs du couple
(force, dplacement).
Afin de contourner ce problme, on pourrait insrer une seconde boucle dans les itrations
plastiques qui chercherait pour un incrment de dplacement donn, quelle valeur de vitesse de
dformation in permet de vrifier la relation implicite de l'quation (155). Ce problme de
point fixe peut tre rsolu par dichotomie, par la mthodes des scantes, ou d'autres encore,
mais toutes consomment beaucoup de temps de calcul (multipli d'un facteur 30 environ).
176
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
el elastic predictor
return path
n
n1
=constant
elastic domain
=0 yield surface
Fig. 128 : Passage du temps n au temps n+1 par prdiction lastique, puis retour la surface
de charge.
Le calcul du prdicteur lastique est ralis depuis la dformation totale de l'itration plastique :
el= D1 tot . On en dduit un tat de contrainte. Si celui-ci est lastique, alors c'est fini.
Dans le cas contraire, il faut effectuer un retour sur la surface de charge dformation totale
constante. Pour cela, on diffrencie la fonction surface de charge :
, in , in
i ini ini i in ini in ini
, , in in (156)
n h
Par rapport l'algorithme initial, seul le dernier terme a t ajout. On propose d'autre part les
relations de rcurrence suivantes :
i 1= i =D1 n p
ini1 ini in in
= avec = p/ t (157)
ini1=ini in in = p
En reportant les relations prcdentes dans l'quation (149), on obtient l'expression du
multiplicateur plastique qui annule :
i
p = (158)
n D1 n h / t
Dtaillons le calcul de pour notre cas particulier.
177
Annexe B : Identification et modlisation du comportement dynamique de l'acier 16MND5
Prenant , in , in =
3
s s y in , in
2 ij ij
[ 1 H ] / [1H 1E-4 ] ,
1 1
et in in
y , =
stat in in p p
y (159)
y
= in = in =
stat
y H
in 1
in p
1
on voit que :
[ ]
1
p
p 1 H 1E-4
Une fois les multiplicateur plastique connu, il ne reste qu' actualiser les variables qui
interviennent dans le calcul. On rpte cette procdure jusqu' obtenir une valeur de la surface
de charge suffisamment proche de 0. On choisit en pratique le critre suivant :
min 1E-6 , 0.1 t y0 .
178
Bibliographie
[1] Griffith AA, The Phenomena of Rupture and Flow in Solids, Philosophical Transactions
of the Royal Society of London, 1921, vol. 221, pp. 163-198.
[2] Irwin GR, Analysis of stresses and strains near the end of a crack traversing a plate,
Journal of Applied Mechanics, 1957, vol. 24, pp. 361-364.
[3] Rgles de Conception et de Construction relatives aux matriels Mcaniques, Annexe
A16 - Tome I - Vol Z, RCC-M, AFCEN, 2002.
[4] Kalthoff JF, On the measurement of dynamic fracture toughness: a review of recent work,
International Journal of Fracture, 1985, vol. 27, pp. 277-298.
[5] Hahn GT, Hoagland RG, Kanninen MF, Rosenfield AR, A Preliminary Study of Fast
Fracture and Arrest in DCB Test Specimen, Proceedings of the Conference on Dynamic
Crack Propagation - Leyden (Netherlands), 1973.
[6] Lemaitre J, Chaboche JL, Mcanique des matriaux solides, Paris : Dunod, 1985.
[7] Bui HD, Mcanique de la rupture fragile, Paris : Masson, 1978.
[8] Kanninen MF, Popelar CH, Advanced fracture mechanics, New York : Oxford University
Press, 1985.
[9] Anderson TL, Fracture mechanics, Cambridge UK : Cambridge University Press, 1995.
[10] Freund LB, Dynamic fracture mechanics, Cambridge UK : Cambridge University Press,
1989.
[11] Rice JR, A path independant integral and the approximate analysis of strain concentration
by notches and cracks, Providence US : Division of Engineering, Brown University,
1967.
[12] Hutchinson JW, Singular behavior at the end of a tensile carck in a hardening material,
Journal of the Mechanics and Physics of Solids, 1968, vol. 16, pp. 13-31.
[13] Rice JR, Rosengren GF, Plane strain deformation near a crack tip in a power law
hardening material, Journal of the Mechanics and Physics of Solids, 1968, vol. 16, pp. 1-
12.
[14] American Society for Testing and Material, Analysis of flaws - Section XI, Appendix A,
ASTM Code, US ASTM, 1986.
[15] Smith DJ, Hadidimoud S, Fowler H, The effects of warm pre-stressing on cleavage
fracture. Part 1: evaluation of experiments., Engineering Fracture Mechanics, 2004, vol.
71, pp. 2015-2032.
[16] Wallin K, The scatter in KIc-results, Engineering Fracture Mechanics, 1984, vol. 19, pp.
1085-1093.
[17] Reytier M, Criter/Fister : conclusions/synthse des travaux effectus sur les critres
d'amorage sous choc thermique dans la transition, Rapport SEMT/LISN/RT/04-037/A,
CEA, 2004.
[18] Wallin K, The effect of ductile tearing on cleavage fracture probability in fracture
toughness testing, Engineering Fracture Mechanics, 1989, vol. 32, pp. 523-531.
[19] Ritchie RO, Knott JF, Rice J, On the relationship between critical tensile stress and
fracture stress in mild steels, Journal of Mechanics and Physics of Solids, 1973, vol. 21,
pp. 395-410.
[20] Beremin FM, A local criterion for cleavage fracture of a nuclear pressure vessel,
Metallurgical transactions A, 1983, vol. 14, pp. 2277-2287.
[21] Reytier M, Chapuliot S, Marie S, Ferry L, Nedelec M, Study of cleavage initiation under
thermal shock by tests on cracked rings and thermomechanical calculations, Nuclear
Engineering and Design, 2006, vol. 236, pp. 1039-1050.
[22] Yin S, Bass R, Williams P, Ludwig M, Keim E, Effect of shallow flaws and biaxial
loading on transition temperatures using a Weibull stress model, ASME PVP Conference
- Cleveland (USA), 2003.
[23] Lefevre W, Barbier G, Masson R, Rousselier G, A modified Beremin model to simulate
the warm pre-stress effect, Nuclear Engineering and Design, 2002, vol. 216, pp. 27-42.
[24] Bordet SR, Karstensen AD, Knowles DM, Wiesner CS, A new statistical local criterion
for cleavage fracture in steel. Part I: model presentation, Engineering Fracture Mechanics,
2005, vol. 72, pp. 435-452.
[25] Bordet SR, Karstensen AD, Knowles DM, Wiesner CS, A new statistical local criterion
for cleavage fracture in steel. Part II: application to an offshore structural steel,
Engineering Fracture Mechanics, 2005, vol. 72, pp. 453-474.
[26] Margolin BZ, Gulenko AG, Nikolaev VA, Ryadkov LN, Prediction of the dependence of
KJC (T) on neutron fluence for RPV steels on the basis of the Unified Curve concept,
International Journal of Pressure Vessels and Piping, 2005, vol. 82, pp. 679-686.
[27] Kalthoff JF, Shockey DA, Instability of cracks under impulse loads, Journal of Applied
Physics, 1977, vol. 48, pp. 986-996.
[28] Freund LB, Energy flux into the tip of an extending crack in an elastic solid, Journal of
Elasticity, 1972, vol. 2, pp. 341-349.
[29] Destuynder P, Djaoua M, Lescure S, Quelques remarques sur la mcanique de la rupture
lastique, Journal de Mcanique Thorique et Applique, 1983, vol. 2, pp. 113-135.
[30] Shih CF, Moran B, Nakamura T, Energy release rate along a three dimensional crack
front in a thermally stressed body, International Journal of Fracture, 1986, vol. 30, pp. 79-
102.
[31] Nikishkov GP, Atluri SN, Calculation of fracture mechanics parameters for an arbitrary
three dimensional crack by the equivalent domain integral method, International Journal
for Numerical Methods in Engineering, 1987, vol. 24, pp. 1801-1821.
[32] Attigui M, Petit C, Valeta MP, Identification des paramtres de fissuration en lasto-
dynamique par la mthode G-Thta, Rapport DMT 95-675, CEA, 1995.
[33] Attigui M, Petit C, Mixed mode separation in dynamic fracture mechanics: New path
independent integrals, International Journal of Fracture, 1997, vol. 84, pp. 19-36.
[34] Rthor J, Mthode des lments finis tendus en espace et en temps: Application la
propagation dynamique des fissures, Thse, Villeurbanne : INSA de Lyon, 2005.
[35] Organ DJ, Numerical solutions to dynamic fracture problems using the element free
Galerkin method, PhD thesis, Chicago : Nothwestern University, 1996.
180
[36] Nishioka T, Computational dynamic fracture mechanics, International Journal of Fracture,
1997, vol. 86, pp. 125-159.
[37] Okada H, Atluri SN, Omori Y, Kobayashi AS, Direct evaluation of T* integral from
experimentally measured near tip displacement field, for a plate with stably propagating
crack, International Journal of plasticity, 1999, vol. 15, pp. 869-897.
[38] Okada H, Atluri SN, Further study of the T* integral; Plane stress stable crack
propagation in ductile materials, Computational Mechanics, 1999, vol. 23, pp. 339-352.
[39] Aoki S, Kishimoto K, Sakata M, Energy-release rate in elastic-plastic fracture problems,
Journal of Applied Mechanics, 1981, vol. 48, pp. 825-829.
[40] Eshelby JD, The elastic energy-momentum tensor, Journal of Elasticity, 1975, vol. 5, pp.
321-335.
[41] Fletcher DC, Conservation laws in linear elastodynamics, Archive for Rational
Mechanics and Analysis, 1976, vol. 60, pp. 329-353.
[42] Mott NF, Brittle fracture in mild steel plates, Engineering, 1947, vol. 165, pp. 16-18.
[43] Broberg KB, The propagation of a brittle crack, Arkiv fur Fysik, 1960, vol. 18, pp. 159-
192.
[44] Rose LRF, Recent theoretical and experimental results on fast brittle fracture,
International Journal of Fracture, 1976, vol. 12, pp. 799-813.
[45] Hahn GT, Hoagland RG, Kanninen MF, Rosenfield AR, Crack arrest in steels,
Engineering Fracture Mechanics, 1975, vol. 7, pp. 583-591.
[46] Kanninen MF, Popelar C, Gehlen PC, Dynamic analysis of crack propagation and arrest
in the double-cantilever-beam specimen, ASTM special technical publication, 1977, vol.
627, pp. 19-38.
[47] Maigre H, Contribution thorique l'identification des grandeurs carctristiques en
mcanique de la rupture dynamique, Thse, Palaiseau : Ecole Polytechnique, 1990.
[48] Kobayashi AS, Seo K, Jou JY, Urabe Y, A dynamic analysis of modified compact-
tension specimens using homalite-100 and polycarbonate plates, Experimental
Mechanics, 1980, vol. 20, pp. 73-79.
[49] Zehnder AT, Rosakis AJ, Dynamic fracture initiation and propagation in 4340 steel under
impact loading, International Journal of Fracture, 1990, vol. 43, pp. 271-285.
[50] Schwartz CW, Crack-speed relations inferrred from large single-edge notched specimens
of a 533B steel, Report ORNL/Sub/79-7778/9, Oak Ridge National Laboratory, 1994.
[51] Kanninen MF, Hudak SJ, Couque HR, Dexter RJ, O'Donoghue PE, Viscoplastic-dynamic
crack propagation: Experimental and analysis research for crack arrest applications in
engineering structures, International Journal of Fracture, 1990, vol. 42, pp. 239-260.
[52] Zhao H, Gary G, On the use of SHPB techniques to determine the dynamic behavior of
materials in the range of small strains,, International Journal of Solids and Structures,
1996, vol. 33, pp. 3363-3375.
[53] Rossoll A, Local approach of ductile cast iron fracture toughness measured by Charpy
test, Diploma Thesis, Centrale Paris - TVFA TU Wien, 1993.
[54] Cowper GR, Symonds PS, Strain-hardening and Strain-rate Effects in the Impact Loading
of Cantilever Beams, Providence US : Division of Applied Mathematics, Brown
University, 1957.
181
[55] Tanguy B, Modlisation de l'essai Charpy par l'approche locale de la rupture. Application
au cas de l'acier 16MND5 dans le domaine de transition, Thse, Paris : Ecole des Mines
de Paris, 2001.
[56] Hajjaj M, Propagation et arrt de fissure dans les cuves de racteurs eau prssurise,
Thse, Paris : Ecole Centrale de Paris, 2006.
[57] Freund LB, Hutchinson JW, Lam PS, Analysis of high strain rate elastic plastic crack
growth, Engineering Fracture Mechanics, 1986, vol. 23, pp. 119-129.
[58] Xu Y, Saigal S, An element free Galerkin analysis of steady dynamic growth of a mode I
crack in elastic plastic materials, International Journal of Solids and Structures, 1999, vol.
36, pp. 1045-1079.
[59] Mataga PA, Freund LB, Hutchinson JW, Crack tip plasticity in dynamic fracture, Journal
of Physics and Chemistry of Solids, 1987, vol. 48, pp. 985-1005.
[60] Wells GN, Sluys LJ, A new method formodelling cohesive cracks using finite elements,
International Journal of Numerical Methods in Engineering, 2001, vol. 50, pp. 2667-
2682.
[61] Remmers JC, De Borst R, Needleman A, The cohesive segment method for the
simulation ofdynamic fracture, IUTAM Symposium on Discretization for evolving
discontinuities - Lyon (France), 2006.
[62] Robertson TS, Propagation of brittle fracture in steel, Journal of the Iron and Steel
Institute, 1953, vol. 175, pp. 361-374.
[63] Wiesner C.S., Predicting structural crack arrest behaviour using small-scale material
characterisation tests, International Journal of Pressure Vessels and Piping, 1996, vol. 69,
pp. 185-196.
[64] American Society for Testing and Material, Standard Test Method for Determining Plane-
Strain Crack-Arrest Fracture Toughness, KIa, of Ferritic Steels, ASTM E1221, US
ASTM, 2006.
[65] Bouyne E, Propagation et arrt de fissure de clivage dans l'acier 2 1/4 Cr-1 Mo, Thse,
Paris : Ecole des Mines de Paris, 1999.
[66] Kanninen MF, An augmented double cantilever beam model for studying crack
propagation and arrest, International Journal of Fracture, 1973, vol. 9, pp. 83-92.
[67] Hoagland RG, Rosenfield AR, Hahn GT, Mechanisms of fast fracture and arrest in steels,
Metallurgical Transactions, 1972, vol. 3, pp. 123-136.
[68] Hahn GT, Hoagland RG, Rosenfield AR, Sejnoha R, Rapid crack propagation in a high
strength steel, Metallurgical Transactions, 1974, vol. 5, pp. 475-482.
[69] Hahn GT, Hoagland RG, Lereim J, Markworth AJ, Rosenfield AR, Fast fracture
toughness and crack arrest toughness of reactor pressure vessel, ASTM special technical
publication, 1980, vol. 711, pp. 289-320.
[70] Smith E, The reinitiation of fracture at the tip of an arrested crack in a Reactor Pressure
Vessel: the effect of ligaments on the reinitiation K value, Journal of Nuclear Materials,
1986, vol. 137, pp. 203-206.
[71] Rosenfield AR, Majumdar BS, A Micromechanical Model for Cleavage-Crack
Reinitiation, Metallurgical Transactions, 1987, vol. 18, pp. 1053-1059.
[72] Kalthoff JF, Modes of dynamic shear failure in solids, International Journal of Fracture,
2000, vol. 101, pp. 1-31.
182
[73] Kalthoff JF, Burgel A, Influence of loading rate on shear fracture toughness for failure
mode transition, International Journal of Impact, 2004, vol. 30, pp. 957-971.
[74] Ramsey JM, Chester FM, Hybrid fracture and the transition from extension to shear
fracture, Nature, 2004, vol. 428, pp. 63-66.
[75] Rittel D, Maigre H, An investigation of dynamic crack initiation in PMMA, Mechanics of
Materials, 1996, vol. 23, pp. 229-239.
[76] Gregoire D, Maigre H, Rthor J, Simulations of dynamic crack propagation experiments
under combined mode loading, 16th European Conference of Fracture - Alexandroupolis
(Grce), 2006.
[77] Gregoire D, Maigre H, Rthor J, Combescure A, Dynamic crack propagation under
mixed-mode loading : Comparison between experiments and X-FEM simulations,
International Journal of Solids and Structures, 2007, vol. 44, pp. 6517-6534.
[78] Genty A, Etude exprimentale et numrique de l'amorage et de l'arrt de fissure sous
choc thermique dans un acier faiblement alli (16MND5), Thse, Paris : Ecole des Mines
de Paris, 1989.
[79] Bryan RH, Bass BR, Merkle JG, The Heavy-Section Steel Technology Pressurized-
Thermal-Shock Experiment, PTSE-1, Engineering Fracture Mechanics, 1986, vol. 23, pp.
81-97.
[80] Iung T, Pineau A, Dynamic crack propagation and crack arrest investigated with a new
specimen geometry: Part I: Experimental and numerical calculations, Fatigue and fracture
of engineering materials and structures, 1996, vol. 19, pp. 1357-1367.
[81] Iung T, Pineau A, Dynamic crack propagation and crack arrest investigated with a new
specimen geometry: Part II: Experimental study on a low-alloy ferritic stee, Fatigue and
fracture of engineering materials and structures, 1996, vol. 19, pp. 1369-1381.
[82] Malluck JF, King WW, Fast fracture simulated by the conventinnal finite elements : A
comparison of two energy-release algorithms, Crack Arrest Methodology and
Applications, 1980, vol. 711, pp. 38-53.
[83] Yuritzinn T, Blanchet N, Rapport d'avancement 2004 de la fiche "Analyse du critre KIa
pour l'arrt de fissure", Rapport SEMT/LISN/RT/04-048/A, CEA, 2004.
[84] Rthor J, Gravouil A, Combescure A, A stable numerical scheme for finite element
simulation of dynamic crack propagation with remeshing, Computer Methods in Applied
Mechanics and Engineering, 2004, vol. 193, pp. 4493-4510.
[85] Nishioka T, Atluri SN, Numerical modeling of dynamic crack propagation infinite bodies,
by moving singular elements, Part 1: Formulation, Journal of Applied Mechanics, 1980,
vol. 47, pp. 570-582.
[86] Huang H, Costanzo F, On the use of space-time finite elements in the solution of elasto-
dynamic fracture problems, International Journal of Fracture, 2004, vol. 127, pp. 119-146.
[87] Abedi R, Haber R, Thite S, Erickson J, An h-adaptive space-time discontinuous Galerkin
method for linear elastodynamics, Revue Europenne de Mcanique Numrique, 2006,
vol. 15, pp. 609-758.
[88] Song JH, Wang H, Belytchko T, A comparative study on finite elemetn methods for
dynamic fracture, Private communication, 2007, vol. , pp. .
[89] Xu XP, Needleman A, Numerical simulations of fast crack growth in brittle solids,
Journal of the Mechanics and Physics of Solids, 1994, vol. 42, pp. 1397-1407.
183
[90] Zhou F, Molinari JF, Dynamic crack propagation with cohesive elements: a methodology
to address mesh dependency, International Journal of Numerical Methods in Engineering,
2004, vol. 59, pp. 1-24.
[91] Camacho G, Ortiz M, Computational modelling of impact damage in brittle materials,
International Journal of Solids and Structures, 1996, vol. 33, pp. 2899-2938.
[92] Pandolfi A, Guduru PR, Ortiz M, Rosakis AJ, Three dimensional cohesive-element
analysis and experiments of dynamic fracture in C300 steel, International Journal of
Solids and Structures, 2000, vol. 37, pp. 3733-3760.
[93] Fedelinski P, Boundary element method in dynamic analysis of structures with cracks,
Engineering Analysis with Boundary Elements, 2004, vol. 28, pp. 1135-1147.
[94] Belytschko T, Krongauz Y, Organ D, Fleming M, Krysl P, Meshless methods: An
overview and recent developments, Computer Methods in Applied Mechanics and
Engineering, 1996, vol. 139, pp. 3-47.
[95] Duflot M, Application des mthodes sans maillage la mcanique de la rupture, Thse,
Lige : Universit de Lige, 2004.
[96] Krysl P, Belytschko T, The element free Galerkin method for dynamic propagation of
arbitrary 3D cracks, International Journal for Numerical Methods in Engineering, 1999,
vol. 44, pp. 767-800.
[97] Rabczuk T, Belytschko T, Adaptivity for structured meshfree particle methods in 2D and
3D, International Journal for Numerical Methods in Engineering, 2005, vol. 63, pp. 1559-
1582.
[98] Buehler MJ, Gao H, Hyperelasticity governs dynamical crack tip instabilities, Nature,
2006, vol. 439, pp. 307-310.
[99] Babuska I, Melenk JM, The partition of unity method, International Journal of Numerical
Methods in Engineering, 1997, vol. 40, pp. 727758.
[100] Chessa J, Wang H, Belytschko T, On the construction of blending elements for local
partition of unity enriched finite element, International Journal of Numerical Methods in
Engineering, 2003, vol. 57, pp. 1015-1038.
[101] Farhat C, Harari I, Franca LP, The discontinuous enrichment method, Computer Methods
in Applied Mechanics and Engineering, 2001, vol. 190, pp. 6455-6479.
[102] Belytschko T, Black T, Elastic crack growth in finite element with minimal remeshing,
International Journal of Numerical Methods in Engineering, 1999, vol. 45, pp. 601620.
[103] Moes N, Dolbow J, Belytschko T, A finite element method without remeshing,
International Journal of Numerical Methods in Engineering, 1999, vol. 46, pp. 131-150.
[104] Stolarska M, Chopp DL, Moes N, Belytschko T, Modelling crack growth by level sets in
the extemended finite element method, International Journal of Numerical Methods in
Engineering, 2001, vol. 51, pp. 943-960.
[105] Gravouil A, Moes N, Belytschko T, Non planar 3d crack growth by the extended finite
element and level sets : Part ii : Level set update, International Journal of Numerical
Methods in Engineering, 2002, vol. 53, pp. 2569-2586.
[106] Sethian JA, Level set methods and fast marching methods, Cambridge UK : Cambridge
University Press, 1999.
[107] Elguedj T, Gravouil A, Combescure A, Appropriate ext ended functions for xfem
simulation of plastic fracture mechanics, Computer Methods in Applied Mechanics and
184
Engineering, 2006, vol. 195, pp. 501-515.
[108] Ribeaucourt R, Gestion du contact avec frottement le long des faces de fissures dans le
cadre de la mthode X-FEM. Application la fatigue tribologique, Thse, Villeurbanne :
INSA de Lyon, 2006.
[109] Zi G, Belytschko T, New crack tip elements for X-FEM and applications to cohesive
cracks, International Journal for Numerical Methods in Engineering, 2003, vol. 57, pp.
2221-2240.
[110] Menouillard T, Rthor J, Combescure A, Bung H, Efficient explicit time stepping for
the eXtended Finite Element Method, International Journal of Numerical Methods in
Engineering, 2006, vol. 68, pp. 911-939.
[111] De Borst R, Remmers J, Needleman, Mesh-independent discrete numerical
representations of cohesive-zone models, Engineering Fracture Mechanics, 2006, vol. 73,
pp. 160-177.
[112] Budyn E, Zi G, Moes N, Belytschko T, A method for multiple crack growth in brittle
materials without remeshing, International Journal of Numerical Methods in Engineering,
2004, vol. 61, pp. 1741-1770.
[113] Liang J, Huang R, Prevost JH, Suo Z, Evolving pattern in thin films with the extended
finite element method, International Journal of Solids and Structures, 2003, vol. 40, pp.
2343-2354.
[114] Guidault PA, Allix O, Champaney L, Cornuault C, Une approche micro-macro pour le
suivi de fissure avec enrichissement local, Revue Europenne de Mcanique Numrique,
2006, vol. 15, pp. 187-198.
[115] Rannou J, Gravouil A, Combescure A, Elments finis tendus multi-grilles avec
raffinement local pour la mcanique de la rupture, 8me colloque national en calcul de
structure - Giens (France), 2007.
[116] Simone A, Wells G, Sluys L, From continuous to discontinuous failure in a gradient-
enhanced continuum damage model, Computer Methods in Applied Mechanics and
Engineering, 2003, vol. 192, pp. 4581-4607.
[117] Osher S, Sethian JA, Fronts propagating with curvature -dependent speed: algorithms
based on hamilton-jacobi formulations, Journal of Computational Physics, 1988, vol. 79,
pp. 12-49.
[118] Gottlieb S, Shu CW, Total Variation Diminishing Runge-Kutta Schemes, Mathematics of
Computation, 1998, vol. 67, pp. 73-85.
[119] Sussman M, Fatemi E, Smereka P, Osher S, An improved level set method for
icompressible two-phase flows, Computer and Fluids, 1998, vol. 27, pp. 663-680.
[120] Osher S, Fedkiw RP, Level set methods- An overview and some recent results, Journal of
Computational Physics, 2001, vol. 169, pp. 463-502.
[121] Sethian JA, Strain J, Crystal growth and dendritic solidification, Journal of Computational
Physics, 1992, vol. 98, pp. 231-253.
[122] Malladi R, Sethian JA, Image Processing: Flows under Min/Max Curvature and Mean
Curvature, Graphical Models and Image Processing, 1996, vol. 58, pp. 127-141.
[123] Burchard P, Cheng LT, Merriman B, Osher S, Motion of curves in three dimensions using
a level set approach, Journal of Computational Physics, 2001, vol. 170, pp. 720-741.
[124] Cheng LT, Burchard P, Merriman B, Osher S, Motion of curves constrained on surfaces
185
using a level set approach, Journal of Computational Physics, 2002, vol. 175, pp. 604-
644.
[125] Duflot M, A study of the representation of cracks with level sets, International Journal for
Numerical Methods in Engineering, 2006, vol. 70, pp. 1261-1302.
[126] Ventura G, Xu JX, Belytschko T, A vector level set method and a new discontinuity
approximations for crack growth by efg, International Journal of Numerical Methods in
Engineering, 2002, vol. 54, pp. 923-944.
[127] Hughes TJR, The Finite Element Method: Linear Static and Dynamic Finite Element
Analysis, Englewoods Cliffs NJ : Prentice-Hall, 1987.
[128] Newmark NM, A method of computation for structural dynamics, Proc. A.S.C.E., 1959,
vol. 85, pp. 6794.
[129] De Borst R, Remmers J, Needleman A, Computational aspect of cohesive zone models,
Advanced Fracture Mechanics for Life and Safety Assesments - Stockholm (Sweden),
2004.
[130] Samaniego E, Belytschko T, Continuum discontinuum modelling of shear band,
International Journal of Numerical Methods in Engineering, 2005, vol. 62, pp. 1857-
1872.
[131] Ventura G, On the elimination of quadrature subcells for discontinuous functions in the
eXtended Finite-Element Method, International Journal of Numerical Methods in
Engineering, 2006, vol. 66, pp. 761-795.
[132] Prabel B, Combescure A, Gravouil A, Marie S, Level set non-matching meshes:
Application to dynamic crack propagation in elastic-plastic media, International Journal
of Numerical Methods in Engineering, 2007, vol. 69, pp. 1553-1569.
[133] Sukumar N, Chopp DL, Moes N, Belytschko T, Modeling holes and inclusions by level
sets in the xfem, Computer Methods in Applied Mechanics and Engineering, 2001, vol.
190, pp. 61836200.
[134] Barth TJ, Sethian JA, Numerical schemes for the hamilton-jacobi and level set equations
on triangulated domains, Journal of Computational Physics, 1998, vol. 145, pp. 1-40.
[135] Sussman M, Fatemi E, An Efficient, Interface-Preserving Level Set Redistancing
Algorithm and Its Application to Interfacial Incompressible Fluid Flow, SIAM Journal on
Scientific Computing, 1999, vol. 20, pp. 1165-1191.
[136] Peng D, Osher S, Merriman B, Zhao H, Kang M, A pde-based fast local level set method,
Journal of Computational Physics, 1999, vol. 155, pp. 410-438.
[137] Gravouil A, Mthode multi-chelles en temps et en espace avec dcomposition de
domaines pour la dynamique non-linaire des structures, Thse, Cachan : Ecole Normale
Suprieure de Cachan, 2000.
[138] Verpeaux P, Charras T, Dynamique du solide : modification du schma de Newmark aux
cas non-linaire, 6me colloque national en calcul de structure - Giens (France), 2003.
[139] Yuritzinn T, Synthse de la fiche "Analyse du critre KIa pour l'arrt de fissure", Rapport
SEMT/LISN/RT/06-003/A, CEA, 2006.
[140] Belytschko T, Chen H, Xu J, Zi G, Dynamic crack propagation based on loss of
hyperbolicity and a new discontinuous enrichment, International Journal of Numerical
Methods in Engineering, 2003, vol. 58, pp. 1873-1905.
[141] Chapuliot S, Marie S, Dubarle PE, Premire campagne de caractrisation de l'acier
186
16MND5 des essais FISTER, Rapport SEMT/LISN/00-001/A, CEA, 2001.
[142] Gary G, Essais de compression dynamique de l'acier FISTER 16MND5, Rapport de
contrat (proposition CEA 5077), LMS X, 2005.
[143] Gary G, Essais de compression dynamique de l'acier FISTER 16MND5, Rapport de
contrat (proposition CEA 6101), LMS X, 2006.
[144] Reytier M, Elements visant la proposition d'un critre de rupture dans la transition
fragile/ductile, Rapport SEMT/LISN/RT/03-035/A, CEA, 2003.
[145] American Society for Testing and Material, A2 - Special requirements for testing compact
specimen, ASTM E1820, US ASTM, 1986.
[146] Nguyen TD, Govindjee S, Klein PA, Gao H, A material force for inelastic fracture
mechanics, Journal of Mechanics and Physics of Solids, 2005, vol. 53, pp. 91-121.
[147] Mariani S, Perego U, Extended finite element method for quasi-brittle fracture,
International Journal of Numerical Methods in Engineering, 2003, vol. 58, pp. 103-126.
[148] Remmers JC, De Borst R, Needleman A, The cohesive segment method for the
simulation ofdynamic fracture, IUTAM Symposium on Discretization for evolving
discontinuities, Lyon (France), 2006.
[149] Rgles de Conception et de Construction relatives aux matriels Mcaniques, Annexe
A16 - Tome I - Vol Z, RCC-M, AFCEN, 2002.
[150] Oliver J, Huespe AE, Theoretical and computational issues in modelling material failure
in strong discontinuity scenarios, Computer Methods in Applied Mechanics and
Engineering, 2004, vol. 193, pp. 2987-3014.
[151] Oliver J, Huespe AE, Continuum approach to material failure in strong discontinuity
settings, Computer Methods in Applied Mechanics and Engineering, 2004, vol. 193, pp.
3195-3220.
[152] Dumstorff P, Meschke G, Crack propagation criteria in the framework of X-FEM-based
structural analyses, International Journal for Numerical and Analytical Methods in
Geomechanics, 2007, vol. 31, pp. 239-259.
[153] Combescure A, Implantation d'un nouvel lment fini dans le code de calcul Castem
2000, Rapport DMT/93 654, CEA, 1993.
[154] Zhao H, A constitutive model for metals over a large range of strain rates: Identification
for mild-steel and aluminium sheets, Materials Science and Engineering, 1997, vol. 230,
pp. 95-99.
[155] Johnson GR, Cook WH, A constitutive model and data for metals subjected to large
strains, high strain rates and high temperatures, Proceedings of the 7th International
Symposium on Ballistics - The Haghe (Netherlands), 1983.
[156] Ortiz M, Simo JC, An analysis of a new class of algorithms for elastoplastic constitutive
relations, International Journal of Numerical Methods in Engineering, 1986, vol. 23, pp.
353-366.
[157] Wang WM, Sluys LJ, De Borst R, Viscoplasticity for instabilities due to strain softening
and strain-rate softening, International Journal of Numerical Methods in Engineering,
1997, vol. 40, pp. 3839-3864.
187
FOLIO ADMINISTRATIF
Prnoms : Benoit
Spcialit : Mcanique
RESUME :
Ce mmoire de thse prsente l'tude de la propagation et de l'arrt d'une fissure de clivage dans un acier de cuve REP.
L'tude bibliographique rappelle quelques lments de mcanique de la rupture en dynamique. Elle montre que la comprhension des
phnomnes associs la propagation instable d'une fissure de clivage n'est pas complte ce jour. Une bonne connaissance exprimentale
des grandeurs en jeu, ainsi qu'un outil numrique capable de reproduire fidlement ces phnomnes, sont des lments permettant de faire
progresser les modles de propagation, notamment en ce qui concerne la dfinition d'un paramtre de fissuration pertinent.
La mthode des lments finis tendus (X-FEM) a t implante dans le code de calcul Cast3m. Cet outil numrique permet de simuler la
propagation arbitraire de fissure sans avoir la mailler. Afin d'amliorer cette mthode, on propose deux techniques. D'une part, on dfinit
les fonctions de niveau qui dcrivent la gomtrie de la fissure sur une grille auxiliaire dcouple du maillage de la structure, pour
permettre une excellente qualit de la reprsentation de la fissure. D'autre part, afin d'tendre l'utilisation des lments enrichis au cas d'un
matriau inlastique, on propose une intgration non conforme des lments proches du front de fissure. Cette technique permet d'intgrer
de manire prcise les fonctions d'enrichissements, et vite la projection de champs au voisinage de la pointe de fissure.
La campagne d'essais de propagation de fissure concerne deux gomtries : des prouvettes Compact Tension, et des prouvettes
annulaires en compression. Une attention particulire est porte aux mesures de vitesse de propagation, car elles constituent une donne
essentielle en mcanique de la rupture dynamique. On relve une vitesse comprise entre 200 et 500m/s pour les prouvettes CT, et de
l'ordre de 700m/s pour les prouvettes annulaires. Le cas d'un chargement en mode mixte a galement t examin. Les analyses de facis
de rupture montrent que l'unique mcanisme de ruine rencontr pour tous les essais est le clivage (absence de ligaments ductiles).
Un modle de propagation bas sur la contrainte principale en pointe de fissure a t identifi en imposant la vitesse de propagation des
prouvettes CT. Une dpendance de la contrainte critique de clivage avec la vitesse de sollicitation (vitesse de propagation ou de
dformation) a t conclue. L'application de ce modle aux prouvettes CT, annulaires en mode I et annulaires en mode mixte, montre qu'il
est possible de prdire par la simulation numrique, le trajet de la fissure, sa vitesse d'avance, ainsi que sa longueur l'arrt, avec une trs
bonne.
MOTS-CLES : rupture, clivage, lments finis tendus, propagation de fissure, arrt de fissure, dynamique, plasticit, 16MND5.