University of Liège

Aerospace & Mechanical Engineering

La propagation des fissures : la quintessence de la mécanique des

matériaux expliquée simplement

Ludovic Noels

Computational & Multiscale Mechanics of Materials – CM3

http://www.ltas-cm3.ulg.ac.be/
Allée de la découverte 9, B4000 Liège
L.Noels@ulg.ac.be

Propagation des fissures

 Table des matières
• Tenseurs contraintes / déformations
• Critère de rupture sur base des contraintes
• Introduction à l’approche de durée de vie totale
• Limite de l’approche de durée de vie totale
• Origine de la mécanique de la rupture
• Rupture fragile/ductile
• Mécanique linéaire de la rupture
• Limites de la mécanique linéaire de la rupture
• Cuve de réacteur nucléaire
• Méthodologie d’analyse
• Propagation de fissures sous chargement cyclique
• Conception en résistance à la fatigue

Propagation des fissures

 Tenseurs Contraintes / Déformations
• Chargement uni-axial l
– La force dépend de la section A
F A0
– Le déplacement dépend de la longueur
X
– Définition de la contrainte L
• Contrainte d’ingénieur

• Contrainte réelle (Cauchy)

• Unités : N/m2 or Pa
– Définition de la déformation
• Unités : -

• Propriétés matériaux s
– Courbe contrainte-déformation sTS
• Module de Young E
• Limite élastique sp0 : Déformation au delà sp0
de laquelle le comportement est
irréversible E E
• Limite en tension sTS
ep ee e

Propagation des fissures

 Tenseurs Contraintes / Déformations
• Vecteurs de traction surfacique szz
– Tx, Ty, Tz
– 𝑻𝑖 = 𝑗 𝜎𝑖𝑗 𝒆𝑗 Tz

– Unités : N/m2 ou Pa
szx szy

• Tenseur des contraintes

syz
– 𝝈 = 𝜎𝑖𝑗 𝒆𝑖 ⊗ 𝒆𝑗
z sxz
– 𝝈T = 𝝈

– Unités: N/m2 ou Pa syx Ty

sxx
Tx sxy syy

x
y

Propagation des fissures

 Rupture Fragile / Ductile

True s
• 2 types de comportement sTS
– Fragile : sp0
• (Presque) pas de déformation
s
avant la rupture (macroscopique) TS
– Ductile : sp0
• Déformations plastiques avant la
rupture (macroscopique) True e

Propagation des fissures

 Avant la mécanique de la rupture
• Conception en limitant les contraintes
– Limite élastique (sp0) ou
– Limite de rupture (sTS)
• ~1860, Wöhler
– Technicien aux chemins de fer allemands
– A étudié la rupture des axes
• La rupture se produit
– Après différents temps de service
s – Pour des charges plus faibles que prédites

w/2 w/2

1 cycle = t

2 reversals
w/2 w/2

• La rupture vient du chargement cyclique

– Approche de durée de vie totale
• Approche empirique de la fatigue

Propagation des fissures

 Approche empirique de durée de vie totale

• Caractérisation du chargement

smax

sa = Ds/2
Ds = smax-smin
1 cycle =
2 reversals sm
t
smin

– Contraintes minimales & maximale : smin & smax

– Contrainte moyenne : sm=(smax + smin)/2
– Amplitude : sa = Ds/2 = (smax - smin)/2
– Rapport de chargement : R = smin / smax
– Pour certaines conditions environnementales particulières (humidité, hautes
températures, …) :
• Fréquence des cycles
• Formes des cycles (sinus, rampes, …)

Propagation des fissures

 Approche empirique de durée de vie totale

• Caractérisation du matériaux (faibles contraintes)

– Pour sm = 0 & Nf cycles identiques avant rupture
• Courbes de Wöhler
sa = Ds/2

se
Nf

103 104 105 106 107

• se : limite d’endurance (vie >107 cycles)

• 1910, Loi de Basquin

– sf’ coefficient de fatigue

– b exposant de fatigue

– Paramètres venant de tests expérimentaux

• Très couteux en temps

Propagation des fissures

 Conception par approche de durée de vie totale
• 1952, De Havilland 106 Comet 1, UK (1)
– Premier avion de transport à réaction
– Cabine pressurisée (0.58 atm)

– Conception du fuselage par approche de

durée de vie totale
• 1952, un fuselage a été testé en fatigue
– Chargement statique à 1.12 atm, puis
– 10 000 cycles à 0.7 atm

– 1953-1954 : accidents par dépressurisation

brutale
• Avril 1954, tests sur un autre fuselage
– Simulation des cycles de pressurisation
de la cabine
– Rupture après seulement 3057 cycles à
partir d’un coin de hublot

Propagation des fissures

 Conception par approche de durée de vie totale
• 1952, De Havilland 106 Comet 1, UK (2)
– Août 1954, toit d’un Comet repêché
• Origine de fissure à la fenêtre de communication
• Fenêtres carrées rivetées concentration de
contraintes
• Rivets perforants Défauts initiaux

• Approche de durée de vie totale

– Prédit l’initiation de fissure par fatigue
– Ne tient pas compte des défauts initiaux
– Test initial du fuselage
• Ecrouissage du métal autour des défauts a cause du
test statique à 1.12 atm fermeture des défauts
• Fuselages de production sans ce pré-chargement

• 1958, Comet 3 et 4
– Fenêtres rondes et collées

Propagation des fissures

 Origine de la mécanique de la rupture
y
• Limites de l’approche de durée de vie totale
s∞
– Ne prends pas en comptes les défauts
– Que se passe-t-il quand un défaut existe?
• Exemple : Le Comet syy syy
• Théorie de la concentration de contraintes
– Plaque infinie avec un trou ellipsoïdal (1913, Inglis) 2b x
2a

s∞
– b→ 0 smax → ∞ rupture pour s∞ → 0
P

• En contradiction avec les expériences de Griffith et Irwin

– Eprouvette en verre avec une fissure de surface (griffe)
– La limite en tension dépend W
» De la taille de la fissure a et a
» De l’énergie de surface gs L

t
– Soit

P
Propagation des fissures
 Rupture fragile/ductile
• Mécanisme de la rupture fragile
sTS

True s
– (Presque) Pas de déformation plastique sp0
avant la rupture (macroscopique)
– Clivage: séparation selon les plans
cristallographiques
• En général a l’intérieur des grains
• Directions préférentielles (à faibles
énergies de liaison) True e
• Parfois entre les grains : corrosion, H2, …
– Critère de rupture
• 1920, Griffith:

Propagation des fissures

 Rupture fragile/ductile

True s
• Mécanisme de la rupture ductile
– Déformations plastiques avant la rupture
(macroscopique) sTS
• Mouvement des dislocations
sp0
• Nucléation de vides autour des inclusions
• Coalescences de microcavités
• Initiation et croissance de fissures True e
– Critère de rupture
• Qu’en est-il du critère de Griffith?

• 1950, Irwin, le travail plastique en pointe de

fissure doit être ajouté

Propagation des fissures

 Rupture fragile/ductile
• “WWII liberty ships”
– Acier à basse T° : fragile
•
• gs ~ 3400 J m-2
– Acier à T° ambiante : ductile
•
• 2gs + Wpl ~ 200 kJ m-2
– Température de transition :
• DBTT (Ductile-Brittle Transition
Temperature)
– Utilisation d’un acier de faible
qualité
• Hiver : DBTT~ température de l’eau
• Lors de la mise à l’eau, les défauts
existants sont devenus critiques
• 30% des « liberty ships » ont
souffert de rupture

Propagation des fissures

 Rupture fragile/ductile
• Développement de la mécanique de la rupture
– Comment prédire la ruine quand une fissure est présente?
– Comment prédire l’évolution d’une fissure sous chargement cyclique?

• Différentes analyses
– Rupture fragile Mécanique linéaire de la rupture
(Linear Elastic Fracture Mechanics-LEFM)

– Rupture avec plasticité confinée en front de fissure

Corrections à la théorie de la y
mécanique linéaire de la rupture a rp
sp

– Rupture ductile avec zones plastiques non-confinées x

Mécanique non-linéaire de la rupture
(Non Linear Fracture Mechanics NLFM)

– Chargements cycliques en deçà du chargement limite

Propagation de fissures (loi de Paris)

Propagation des fissures

 Rupture fragile/ductile
• Développement de la mécanique de la rupture
– Comment prédire la ruine quand une fissure est présente?
– Comment prédire l’évolution d’une fissure sous chargement cyclique?

• Différentes analyses
– Rupture fragile Mécanique linéaire de la rupture
(Linear Elastic Fracture Mechanics-LEFM)

– Rupture avec plasticité confinée en front de fissure

Corrections à la théorie de la y
mécanique linéaire de la rupture a rp
sp

– Rupture ductile avec zones plastiques non-confinées x

Mécanique non-linéaire de la rupture
(Non Linear Fracture Mechanics NLFM)

– Chargements cycliques en deçà du chargement limite

Propagation de fissures (loi de Paris)

Propagation des fissures

 Mécanique linéaire de la rupture
• Singularité en pointe fissure pour un matériau élastique linéaire
– 1957, Irwin, 3 modes de rupture
Mode I Mode II Mode III
(ouverture) (glissement) (Cisaillement)
y x
y x y x
z
z z

r
– Définis par les conditions limites q
x
Mode I Mode II Mode III
or or

Propagation des fissures

 Mécanique linéaire de la rupture
• Singularité en pointe fissure pour un matériau élastique linéaire (2)
– Solutions asymptotiques
Mode I Mode II Mode III

– Introduction des facteurs d’ intensité de contraintes (Pa m1/2) - Stress

Intensity Factors – SIF

– Les Ki dépendent
• De la géométrie &

• Du chargement
Propagation des fissures
 Mécanique linéaire de la rupture

• 1957, Irwin : Nouveau critère de rupture

– smax → ∞ s ne peut pas être utilisé
– Comparer le facteur d’intensité de contraintes Ki (dépendant du problème)
à une valeur propre au matériau : la ténacité KiC
• Si Ki ≥ KiC la fissure se propage
• Ténacité Kic (toughness)
– Métaux : voir figure
– Béton : 0.2 - 1.4 MPa m1/2

Toughness KIC [MPa 𝑚] Toughness KIC [MPa 𝑚]

Pure Al
Trip 200 Aisi 403 12 cr
Ductile
140 Mild steels
Low steel 160
steel
alloys
100
Ductile/brittle
120 transition regime
60 Composites 80
Ti alloys
40
20 Al alloys s0Y [MPa] Brittle
Yield Temperature T [C]
0 0
0 500 1000 1500 2000 -200 -100 0 100 200

Propagation des fissures

 Mécanique linéaire de la rupture

• Evaluation des facteurs d’intensité de contrainte (SIF)

– Analytique (fissure 2a dans un plan ∞)
y y y

s∞ t∞
t∞
s∞ x s∞ x x
2a 2a 2a
t∞ t∞

s∞

– Numérique

• bi dépendent
– De la géométrie
– De la taille de fissure

Propagation des fissures

 Mécanique linéaire de la rupture

• Mesure expérimentale de la ténacité KIc

P, d
– En suivant la norme ASTM E399
– Préparation Thickness t
• Par exemple en considérant un Single Edge W
Notch Bending (SENB) spécimen v/2 v/2
• Usiné avec une encoche
pour initier la fissure L=4W
• Chargement cyclique pour initier la fissure
– Mesure de la ténacité
• Mesure la courbe P - d
• Crack Mouth Opening Displacement (CMOD=v) CMOD, v
mesuré en tête de fissure
• Pc obtenu par les courbes P-v

• KIc est déduit de Pc à partir de formules P P

de référence (éléments-finis) Pc Pc

• f(a/W) dépend du test (SENB, …) tg 0.95 tg tg 0.95 tg

v
Propagation des fissures
 Limites de la mécanique linéaire de la rupture

• La contrainte prédite est singulière syy

Zone of asymptotic Asymptotic syy
– Il existe des non-linéarités
solution (in terms of
(plasticité en pointe de fissure) 1/r1/2) dominance
True syy
– Solution asymptotique non-valide
loin de la pointe de fissure x

• Hypothèse de plasticité confinée

– Valide si la zone de plasticité est petite par rapport à la fissure
– Zone cohésive rp en pointe de fissure
y
a rp
– 1960, Dugdale-Barenblatt sp
• Mode I x

Méthode iterative
• Contrainte nominale est supérieure à la moitie de la limite élastique
– Hypothèse ne tient pas mécanique non linéaire de la rupture (NLFM)
Propagation des fissures
 Reactor Pressure Vessels (RPV) ou cuve de réacteur

• Pressurized Water Reactor (PWR)

– Conditions d’utilisation typiques d’une cuve : 15 MPa, 300º C
Ø ~ 3-5 m

Image from US. Nuclear Regulatory Commision,

Image from US. Nuclear Regulatory Commision,
Quadrennial Technology Review 2015, U.S. Department of Energy
http://www.eia.doe.gov/cneaf/nuclear/page/nuc_reactors/pwr.html

Propagation des fissures

 Reactor Pressure Vessels (RPV) ou cuve de réacteur

• Matériaux de la cuve
– Matériau de base
• Alliage Manganèse–nickel–molybdène
Ø ~ 3-5 m
• E.g. SA508 Cl. 3
– Revêtement
• Inox

tbase ~ 0.2 m tclad ~ 0.01 m Image from US. Nuclear Regulatory Commision,
http://www.eia.doe.gov/cneaf/nuclear/page/nuc_reactors/pwr.html

Propagation des fissures 27

 Reactor Pressure Vessels (RPV) ou cuve de réacteur

• Matériaux de la cuve : présente une transition ductile/fragile

– Définition de T100 : Température de référence pour laquelle 50% des
spécimens normalisés rompent avec une ténacité de 100 [MPa 𝑚]

Conditions
d’utilisation
Nécessite
NLFM

Integrity of reactor pressure vessels in nuclear power plants : assessment of irradiation embrittlement effects in reactor pressure vessel steels. International Atomic Energy Agency, 2009.
(IAEA nuclear energy series), ISSN 1995-7807, ISBN 978-92-0-101709-3

Propagation des fissures 28

 Reactor Pressure Vessels (RPV) ou cuve de réacteur

• Pressurized Thermal Shock (PTS)

– Perte de liquide de refroidissement (LOCA)
– Injection d’eau froide
– Chute brutale de la température de paroi
– Chute de la ténacité
– N.B. chute de la pression, mais existence de
contraintes thermiques

Toughness KIC [MPa 𝑚]

200 Aisi 403 12 cr

Ductile
Christopher Boyd, Interactions of Thermal-Hydraulics with Fuel Behavior,
Structural Mechanics, and Computational Fluid Dynamics
160 Office of Nuclear Regulatory Research Nuclear Regulatory Commission

Ductile/brittle
120 transition regime
80

40
Brittle Temperature T [C]
0
-200 -100 0 100 200

Fictitious material

Propagation des fissures 29

 Reactor Pressure Vessels (RPV) ou cuve de réacteur

• Fragilisation par irradiation

– Irradiation durant l’utilisation du RPV
– Fragilisation augmente avec
• La fluence (nombre de neutrons/unité de surface) décalage température
• L’énergie des neutron (>1MeV) de transition avec l'âge

Toughness KIC [MPa 𝑚] Ductile Un-irradiated

250
Loss
200 Ductile/brittle Irradiated
transition regime
150
Shift
100 Brittle

50
Temperature T [C]
0
-200 -100 0 100 200
Fictitious material

Propagation des fissures 30

 Reactor Pressure Vessels (RPV) ou cuve de réacteur

• Fragilisation par irradiation PTS

– Définie par un décalage de T100
• T0 = T100 + DT100

– Evaluée en termes de
• La composition chimique
• La fluence (Fn)

DT100
– Mesurée
• Spécimens de surveillance
• Proches de la paroi

Integrity of reactor pressure vessels in nuclear power plants : assessment

of irradiation embrittlement effects in reactor pressure vessel steels.
International Atomic Energy Agency, 2009. (IAEA nuclear energy series),
ISSN 1995-7807, ISBN 978-92-0-101709-3

Propagation des fissures 31

 Evaluation

• Méthodologie

Caractérisation des Conditions de

Propriétés matérielles
défauts chargement

• Ténacité du matériau
• Conditions critiques
• Défauts supposés de base (expériences)
(suivant la norme)
• Inspections non-
destructives en
service (Ultrasonique
• Evaluation des • Conditions
…)
températures, environnementales
contraintes (fluence en fin de vie..)

• Evaluation des facteurs d’intensité de contrainte • Ténacité en service

• Méthodes analytiques • Modèles
• Méthodes numériques • Expériences

Propagation des fissures 32

 Références
• Notes
– Lecture Notes on Fracture Mechanics, Alan T. Zehnder, Cornell University,
Ithaca, http://hdl.handle.net/1813/3075
– Fracture Mechanics Online Class, L. Noels, ULg, http://www.ltas-
cm3.ulg.ac.be/FractureMechanics
– Fracture Mechanics, Piet Schreurs, TUe,
http://www.mate.tue.nl/~piet/edu/frm/sht/bmsht.html
• Livres
– Fracture Mechanics: Fundamentals and applications, D. T. Anderson. CRC
press, 1991.
• RPV
– Documents
• Integrity of reactor pressure vessels in nuclear power plants : assessment of
irradiation embrittlement effects in reactor pressure vessel steels. International
Atomic Energy Agency, 2009. (IAEA nuclear energy series), ISSN 1995-7807,
ISBN 978-92-0-101709-3
• Pressurized thermal shock in nuclear power plants: good practices for
assessment, International Atomic Energy Agency, 2010 (iaea-tecdoc-1627), ISSN
1011-4289, ISBN 978-92-0-111109-8

Propagation des fissures

 Chargement cyclique

• Rupture par fatigue

– Tests pour différents DP = Pmax - Pmin P

– Nucléation : fissure initiée pour K < Kc

• Surface : déformations résultant
de la propagation des dislocations a

P
N cycles

Pmax Pmin

Persistent slip
band (PSB)

Propagation des fissures

 Chargement cyclique

• Rupture par fatigue (2) Perfect crystal

– Etape I de croissance de fissures
• Le long d’un plan cristallin
– Etape II de croissance de fissures
• A travers plusieurs grains
– Le long des plans cristallins
– Droit macroscopiquement
• Striation de la surface de rupture :
correspond aux cycles

Propagation des fissures

 Chargement cyclique

• Rupture par fatigue (3) P

– Hypothèse de plasticité confinée 5

• Paramètres d’influence

a (cm)
– DP & 4 a DP2 > DP1
Rupture
– Pmin / Pmax
3 DP1
• Fatigue décrite par P
– DK = Kmax- Kmin & 2
– R = Kmin/Kmax
Nf (105)
1 2 3 4 5 6 7

Inspection structurelle Croissance

possible (pour DP1) rapide

Intervalle de striation
• Il existe DKth tel que si DK ~< DKth : la fissure est dormante

Propagation des fissures

 Chargement cyclique

• Taux de croissance

da / d Nf
Zone I Zone II Zone III
– Zone I
• Etape I de croissance de fissure R
• DKth dépend de R
m
– Zone II
• Etape II de croissance de fissure 1
• 1963, Paris-Erdogan
DKth logDK

– Dépend de la géométrie, du chargement, de la fréquence

– Acier : DKth ~ 2-5 MPa m1/2 , C ~ 0.07-0.11 10-11 [m (MPa m1/2)-m], m ~ 4
– Acier en environnement salin : DKth ~ 1-1.5 MPa m1/2, C ~ 1.6 10-11 [idem],
m ~ 3.3
• Attention : K dépend de a intégration requise pour avoir a(Nf)
– Mode I :
– Zone III
• Croissance rapide jusque rupture
• Rupture statique (clivage) sous Kmax(a): il existe une taille critique af telle que
Kmax(af) = Kc
Propagation des fissures
 Conception
• « Infinite life design »
– sa < se : vie « infinie »
– Non économique
• « Safe life design »
– Pas de fissure avant un nombre défini de cycles
• Le composant est remplacé à la fin de la vie
prédite, même s’il n’y a pas de fissure apparue
• Emphase sur la prévention d’initiation
• Approche empirique par nature
– Suppose que la structure est vierge de
défauts
– Se base sur les courbes sa – Nf
• Ajouter un facteur de sécurité
– Application : aubes de rotors/stators
• Une fois qu’une fissure apparait, la durée de vie
en heures est très courte

Propagation des fissures

 Conception

• « Fail safe design »

– Même si un composant rompt, la structure doit garder assez d’intégrité pour
être utilisée en toute sécurité
– Prédéfinition de chemins de fissuration et d’arrêt de fissure
– Nécessite des inspections périodiques

– Exemple: 1988, B737, Aloha Airlines 243

• 2 coques de fuselage non collées
• Environnement marin rouille et augmentation du volume
• Fatigue des rivets
• La fissure a suivi le chemin prédéfini

Propagation des fissures

 Conception

• « Damage tolerant design »

– Suppose l’existence de défauts au début des opérations
– Caractériser l’effet des défauts sur la vie de la structure
• Evaluation des défauts initiaux (A partir
du procédé de fabrication et
d’inspections non-destructives
• Estimation de la croissance &
élaboration des campagnes d’inspection
(e.g. en années, en nombre de vol…)
• Vérifier la croissance durant les
inspections
• Prédire la fin de vie (af)
• Retirer ou réparer la structure avant

– Inspections non-destructives
• Optique
• Rayons-X
• Ultrasonique (réflexion sur les surfaces de fissure)

Propagation des fissures