2- Principe énergétique

Théorie de Griffith

Elle est développée pour des matériaux présentant un comportement élastique linéaire isotrope.
Elle se base sur le premier principe de la thermodynamique. Le travail d’une sollicitation externe se transforme en
énergie interne, cette énergie se décompose en une énergie de déformation élastique qui se manifeste par la
déformation autour la fissure et une énergie dissipée lors de la création d’une nouvelle surface de fissuration.

Griffith a adopté une méthode de calcul de la variation de l'énergie de

déformation due à l’apparition d’une fissure interne située symétriquement
dans une plaque infinie en utilisant la solution d'Inglis.
Par un calcul approché, en considère une fissure centrale dans une plaque
infinie d'épaisseur uniforme soumise à une tension uni-axiale dans la direction
perpendiculaire à la fissure.

BEST FOR You 9

O R GAN I C S C O M P AN Y

2- Principe énergétique

Sous un chargement appliqué, la fissure s'étend de 0 à 2a. On considère deux régions

assimilées à deux triangles (ABC et ABD) qui vont subir une relaxation complète du champ
de contrainte. Ainsi, la réduction de l'énergie de déformation (We) lorsque la géométrie
passe de la plaque sans fissure à une plaque avec une fissure est égale à l'énergie de
déformation élastique du matériau initialement emmagasinée par ACBD.
Ainsi et pour une épaisseur unitaire (e=1) :

où E est le module de Young. (γs) est l'énergie de surface par unité de surface. L'énergie de
surface Ws de la plaque due à l’apparition de la fissure est donnée par :

En utilisant la solution d'Inglis (β=π) et en exprimant la variation de (We) suite à une

variation de (a).

L'augmentation correspondante Ws de l'énergie de surface suite à une variation de (a) est

donnée par :
BEST FOR You 10
O R GAN I C S C O M P AN Y

4
 2- Principe énergétique

Pour que la fissure progresse, le taux de libération de l'énergie de déformation par unité d'extension de la fissure doit
être supérieur ou égal au taux d'augmentation de l'énergie de surface.
C'est-à-dire :

Par conséquent, la contrainte de rupture ou contrainte critique en présence d'une

fissure de taille 2a est donnée par

Où Gc=2γs est l’énergie requise pour l’extension d’une unité de langueur de la fissure.

On note G comme le taux de libération de l’énergie de déformation suite à la progression de la fissure par une
unité de longueur : 𝜎²𝜋𝑎
G=
𝐸
La condition de propagation de BEST
la fissure
FORpeut se mettre sous la forme ∶ G ≥ 𝐺𝐶
You
O R GAN I C S C O M P AN Y
11

2- Principe énergétique
Effet de la déformation plastique.
L’expression précédente est développée en supposant que l’énergie de déformation élastique se transforme en sa
totalité en une énergie de surface. Cependant, et on plus de la création d’une fissure, une zone de déformation
plastique sera également créée au voisinage, ainsi l’énergie de déformation plastique sera transformée en une énergie
de surface et une énergie de déformation plastique.

Dans ce cas :

BEST FOR You 12

O R GAN I C S C O M P AN Y

5
 3- Approche basée sur le facteur d’intensité de contrainte

BEST FOR You 13

O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Hypothèse : Comportement parfaitement linéaire, milieu infini (dimension du défaut sont très inférieures aux
dimensions de la pièce).

Etat de contrainte au voisinage d’un Etat de contrainte au voisinage d’un

défaut circulaire défaut elliptique

rayon de courbure (ρ)

BEST FOR You 14

O R GAN I C S C O M P AN Y

6
3- Approche basée sur le facteur d’intensité de contrainte

Une fissure peut être assimilée à un défaut elliptique de rayon d’entaille nulle.
Dans ce cas, la contrainte au fond de fissure est infinie.

BEST FOR You 15

O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Sous l’hypothèse d’un comportement élastique, linéaire et isotrope. Plusieurs auteurs (Irwin, Westergaard…) ont
proposé une solution au problème des équations d’équilibres en milieu continu. Si on considère un système de
repérage en coordonnées cylindriques (r,𝜃), il est possible de montrer que l’expression du champ de contrainte en
cas de fissure s’exprime sous la forme suivante :

(fij) et (gij) sont des fonctions de (𝜃) sans dimensions.

Sous les hypothèses précédentes, l’expression de (σij) dépendra de :

• La position du point de calcul ;
• La taille et la géométrie de la fissure ;
• Les paramètres géométriques de la pièce ;
• La natureBEST FOR You externe…
du chargement 16
O R GAN I C S C O M P AN Y

7
 3- Approche basée sur le facteur d’intensité de contrainte

Loin de la fissure, le premier terme tend vers (0), alors que la contrainte σyy doit tendre vers (σ).
Au voisinage de la fissure (r→0), le premier terme (à gauche) est le plus prépondérant et l’évolution de la contrainte
σyy dépend principalement de (r).
Pour toute fissure, on peut donc approcher l’évolution de la contrainte par une même expression au voisinage de la
fissure.

Pour plus de commodité on pose k= K/ 2𝜋

BEST FOR
K = Facteur d’intensité deYou
contrainte (FIC-SIF) 17
O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Modes de chargements élémentaires au niveau d’une fissure :

Le mode (II) est le mode par Le mode (III) est le mode par
Le mode (I) ou le mode par
glissement de translation, le glissement de rotation, le chargement
ouverture. C’est le mode le
chargement est appliqué est toujours parallèle dans le plan de
mode le plus dangereux. Le
parallèlement au plan de fissure et fissure mais perpendiculairement à la
chargement est appliqué
à la direction de sa propagation. direction de propagation .
perpendiculairement au plan
de propagation de la fissure

Il est possible d’analyser les modes de chargement en analysant le vecteur contrainte appliqué sur le plan de fissure.
Ce vecteur peut être décomposé en une composante normale au plan de fissure est une composante tangentielle.
La composante tangentielle se décompose également on deux composantes donnant les modes (2) et (3).
Un cas de chargement arbitraire induit l’ensemble des trois modes.
BEST FOR You 18
O R GAN I C S C O M P AN Y

8
3- Approche basée sur le facteur d’intensité de contrainte

Modes de chargements élémentaires au niveau d’une fissure :

Sous l’hypothèse de la linéarité de la loi de comportement, le tenseur de déformation résultant d’un chargement
aléatoire peut être décomposé en une somme des tenseurs élémentaires induits par chaque mode de chargement.
le même raisonnement est valable pour le tenseur de contrainte :

𝑇 𝐼 𝐼𝐼 𝐼𝐼𝐼
𝜀𝑖𝑗 = 𝜀𝑖𝑗 + 𝜀𝑖𝑗 + 𝜀𝑖𝑗 σ𝑇𝑖𝑗 = σ𝐼𝑖𝑗 + σ𝐼𝐼 𝐼𝐼𝐼
𝑖𝑗 + σ𝑖𝑗

BEST FOR You 19

O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Etat de contrainte pour chaque mode de chargement :

Il est possible d’identifier, pour chaque mode de chargement, une fonction (fij) et un facteur (k) associé. Ainsi, au
voisinage du fond de fissure, il est possible d’exprimer les composantes du tenseur de contrainte sous la forme
suivante :

Les facteurs (KI, KII, KIII) dépendent du chargement externe, de la géométrie de la pièce et de la taille de la
fissure.

BEST FOR You 20

O R GAN I C S C O M P AN Y

9
 3- Approche basée sur le facteur d’intensité de contrainte

Le tableau suivant résume l’expression des composantes des tenseurs de contraintes pour chaque
mode de chargement au voisinage de la fissure :

BEST FOR You 21

O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Dans le cas du mode (I) et lorsque (𝜃=0), la contrainte de cisaillement (τxy) est nulle.
A (𝜃=0), le plan contenant la fissure est soumis à une contrainte normale (σyy) sans cisaillement et (σxx =σyy).

Loin de la fissure, l’état de contrainte est régi par les conditions aux limites.
Au voisinage du fond de la fissure, la contrainte est régie par les équations du tableau précédent.
La connaissance de la valeur de (KI) permet de définir complètement l’état de contrainte au voisinage du fond de fissure.
En raison de cette dépendance, il est possible de juger la tenue mécanique d’une pièce fissurée en évaluant le facteur (K).
BEST FOR You 22
O R GAN I C S C O M P AN Y

10
 3- Approche basée sur le facteur d’intensité de contrainte

Le facteur d'intensité de contrainte définit complètement les conditions en fond de fissure.

Si K est connu, il est possible de résoudre toutes les composantes de contrainte, de déformation et
de déplacement en fonction de r et θ.
Cette description à paramètre unique de l’état en fond de fissure s'avère être l'un des concepts les
plus importants de la mécanique de la rupture.

BEST FOR You 23

O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Expression du facteur d’intensité de contrainte

Le facteur d'intensité de contrainte (K) est déterminé en fonction des sollicitations externes, de
la géométrie de la pièce et de la fissure.
Les expressions de (K) ont été dérivées pour un certain nombre de configurations simples.
Pour des situations plus complexes, le facteur d'intensité de contrainte peut être estimé par
analyse numérique.

l'intensité de contrainte a comme unités [contrainte]*[ 𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑒].

Une configuration pour laquelle il existe une solution est une fissure traversante dans une
plaque infinie soumise à une contrainte de traction à distance.

L’expression du facteur (KI) est :

BEST FOR You 24

O R GAN I C S C O M P AN Y

11
3- Approche basée sur le facteur d’intensité de contrainte

Dans le cas d’une fissure sur une plaque infinie, l’expression du facteur d’intensité de contrainte se met sous la forme :

Il est important de constater que les défauts en surface présentent un danger plus
important que les défauts internes.

BEST FOR You 25

O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Dans le cas ou une fissure présente une inclinaison par rapport au

chargement uni-axial externe appliqué, la fissure présentera une
combinaison entre le mode (I) et (II), (KIII=0).
Si on construit un nouveau repère, (x’: suivant la fissure) (y’: normal à la
fissure), il est possible de décomposer le vecteur contrainte appliqué en
une composante normale et tangentielle à la surface de la fissure.

BEST FOR You 26

O R GAN I C S C O M P AN Y

12
 3- Approche basée sur le facteur d’intensité de contrainte
Fissure en milieu 3D

Dans le cas d’une fissure circulaire de rayon (a) dans un volume infini
soumis à un chargement uni-axial, l’expression du facteur d’intensité de
contrainte s’exprime par :

Dans le cas d’une fissure elliptique interne ou débouchante dans un milieu infini. (Grand axe = 2c,
Petit axe =2a)
Fissure Interne Fissure Externe

Les solutions ne sont valables que si la fissure est petite

BEST FORpar You
rapport aux
27
dimensions de la plaque et a≤c.
O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Ténacité (Toughness)

C’est une propriété mécanique qui décrit la résistance d’un matériau à la propagation d’une fissure déjà existante.

Dans le cas d’une fissure initiée au niveau d’une pièce, l’état de contrainte peut être décrit au voisinage de la fissure
par le facteur d’intensité de contrainte (KI) dans le cas d’une sollicitation en mode (I).
Lorsque l’état de contrainte au voisinage de cette fissure dépasse un état limite, la fissure subira une propagation.
Lorsque la propagation de la fissure est initiée, il est possible de décrire cet état limite par son facteur d’intensité de
contrainte associé qui sera notée le facteur d’intensité de contrainte critique (KIC).

Plus ce facteur est important, plus la résistance du

BEST FOR
matériau face aux fissures existantes estYou
importante. 28
O R GAN I C S C O M P AN Y

13
 3- Approche basée sur le facteur d’intensité de contrainte

Effet de l’épaisseur de l’éprouvette

La ténacité est une propriété mécanique qui dépend de :

• La température d’essai.
• L’épaisseur de l’éprouvette (Etat de contraintes).
• La vitesse de l’essai.

BEST FOR You 31

O R GAN I C S C O M P AN Y

Application
Une plaque d’acier d’une épaisseur de 30 mm et de largeur 12 m est soumise à un chargement uni axial en traction de
50 MPa. Sa ténacité est de 28.3 MPa.
Si une fissure au milieu de la plaque présente une fissure de longueur de 65mm, estimer la contrainte nécessaire à la
propagation brutale de la fissure et comparer cette contrainte à la limite élastique du matériau qui est de 240 Mpa .

BEST FOR You 32

O R GAN I C S C O M P AN Y

15
 3- Approche basée sur le facteur d’intensité de contrainte

Effet du comportement plastique

La prise en considération de la plasticité au fond de l’entaille permet d‘obtenir la zone affectée par une déformation
permanente. Dans ce cas, la contrainte au voisinage du fond d‘entaille présente une valeur finie.
Une première approximation de la limite de la zone plastique (ry) sur le plan de la fissure est obtenue en remplaçant,
dans l‘expression de la contrainte normale d‘ouverture, la contrainte (σyy) par (Re).

Sur le plan de la fissure (𝜃=0)

L’intégration du comportement élasto-plastique dans le calcul de l’évolution du champ de contrainte permet

d’obtenir une approximation plus réalistique de l’étendue de la zone plastique :

BEST FOR You 35

O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Zone de plasticité au voisinage du fond de fissure

Pour estimer la géométrie de la zone de plasticité, il est

possible de se baser sur le critère de l’énergie de distorsion.

Pour un état plan de contrainte, les contraintes principales

s’expriment par :

Dans le cas d’un état plan de déformation, il faudra ajouter la composante (σ3):

Si on considère un mode 1 de chargement, et on remplace les

contraintes déjà obtenues pour identifier les contraintes principales :

BEST FOR You 36

O R GAN I C S C O M P AN Y

16
3- Approche basée sur le facteur d’intensité de contrainte

On trouve : Etat plan de contrainte

Etat plan de
déformation

En intégrant les expressions précédentes dans l’expression de la contrainte équivalente et si on remplace la

contrainte équivalente par la limite élastique, on peut déduire l’expression de la coordonnée (r) pour chaque
position en fonction (𝜃).

Etat plan de contrainte

Etat plan de déformation

Ces équations se basent sur l’approximation simplifiée de la zone de plasticité

BEST FOR You 37
O R GAN I C S C O M P AN Y

3- Approche basée sur le facteur d’intensité de contrainte

Plasticité au voisinage de la fissure estimée à partir de la solution élastique, et la contrainte équivalente de Von Mises

BEST FOR You 38

O R GAN I C S C O M P AN Y

17
 4- Propagation des fissures
sous chargement cyclique

BEST FOR You 41

O R GAN I C S C O M P AN Y

4- Propagation des fissures sous chargement cyclique

Sous chargement statique, l’endommagement brutal d’un matériau présentant

une fissure de langueur (a) et qui est soumis à un mode (I) est conditionné par la
condition :

KI≥ KIC

Sous chargement cyclique appliqué sur une fissure préexistante. Il est possible
d’amorcer la propagation de la fissure à des niveaux de chargement inférieurs au
seuil décrit par la formule précédente.
Cela est causé principalement par la formation des bandes de glissement au
voisinage de la fissure.

BEST FOR You 42

O R GAN I C S C O M P AN Y

18
 4- Propagation des fissures sous chargement cyclique
Loi de propagation de fissure par fatigue

En présence d’une fissure initiale sollicitée en mode (I), l’application d’un

chargement cyclique conduit à une fluctuation cyclique du facteur d’intensité
de contrainte.
𝐾𝐼 𝑡 = Y(géométrie) σ𝑂 (𝑡) 𝜋𝑎

Similairement au calcul en fatigue, il est possible d’associer à un chargement

cyclique en présence d’une fissure un rapport (R), un FIC moyen….

σ𝑀𝑖𝑛 𝐾𝑀𝑖𝑛
𝑅= = ∆𝐾 = 𝐾𝑀𝑎𝑥 − 𝐾𝑀𝑖𝑛
σ𝑀𝑎𝑥 𝐾𝑀𝑎𝑥

BEST FOR You 45

O R GAN I C S C O M P AN Y

4- Propagation des fissures sous chargement cyclique

Observations expérimentales

Si on considère 3 éprouvettes (E1, E2, E3) munies de fissures initiales identiques (a0) et soumises à 3
chargements différents (S1>S2 >S3) :

• L’augmentation de la contrainte induit une augmentation de la vitesse de fissuration (da/dN).

• La taille finale de la fissure (af)avant la rupture finale et d’autant plus courte que le chargement est important.
• La durée de vie durant la phase de propagation de fissure se réduit en augmentant le chargement.
• Notant que la rupture finale dépend de KIc du matériau et de la contrainte appliquée.

R=0

BEST FOR You 46

O R GAN I C S C O M P AN Y

20
 4- Propagation des fissures sous chargement cyclique

Observations expérimentales

A partir de la courbe précédente, il est possible d’identifier à chaque cycle de

chargement :
Propagation
ΔK (t) : Variation du FIC en considérant la longueur de la fissure actuelle. de la fissure
(da/dN)(t) : La progression de la fissure après l’application du cycle de chargement.

Le tracé de (da/dN)(t) en fonction de ΔK(t) donne le résultat suivant

Interprétation :
- La progression de la fissure après un cycle de chargement dépend de
ΔK appliqué à ce cycle qui dépend de la taille de la fissure actuelle.

- A des valeurs spécifiques de (ΔK=2) très faibles, la vitesse de la fissure

est inférieure à 10−10 (détection quasi-impossible) .
- A des valeurs de ΔK inférieures (2), la fissure est stable.
BEST FOR You 47
O R GAN I C S C O M P AN Y

4- Propagation des fissures sous chargement cyclique

Zone 1 : Caractérisée par une propagation très lente de la fissure

(10−10m/cycle). Ainsi le facteur d’intensité de contrainte (ΔK) ne subit pas
de variation, et la fissure préserve une vitesse de progression casi-nulle.
Dans cette même zone on trouve (ΔKth) qui définit le seuil nécessaire pour
amorcée la propagation.

Zone 2 : Caractérisée par une dépendance linéaire entre (ΔK) et (da/dN) sur
une échelle doublement logarithmique.

𝑑𝑎
Modèle linéaire proposé par Paris = 𝐴(Δ𝐾)𝑛
𝑑𝑁
A, n : Deux paramètres déterminés expérimentalement en
traçant la courbe de Paris

Zone 3 : La vitesse de propagation change de comportement et

devient plus importante. Cette zone est délimitée par Kc.

BEST FOR You 48

O R GAN I C S C O M P AN Y

21
 4- Propagation des fissures sous chargement cyclique

Estimation de la durée de vie en fatigue

Lors de la conception d’une nouvelle pièce à base d’un matériau spécifique. Il est important d’évaluer l’étendue du
cycle de vie en cas d’apparition des microfissures suite aux procédés d’obtentions ou suite à l’exploitation et la
maintenance. Cela est obtenu en évaluant d’abord la ténacité et les propriétés relatives à la propagation de la
fissure.
Sachant la longueur initiale (a0) et la longueur finale (af) d’une fissure avant l’éclatement final, il est possible de se
baser sur la loi de Paris pour estimer la durée de vie.

BEST FOR You 49

O R GAN I C S C O M P AN Y

4- Propagation des fissures sous chargement cyclique

BEST FOR You 50

O R GAN I C S C O M P AN Y

22