INBTP
COURS DE MECANIQUE DE LA
RUPTURE
CHAPITRE 1:
GENERALITES SUR LA RUPTURE
Rostand MOUTOU PITTI
HDR
RMP RUPTURE - Février 2023
PLAN DU CHAPITRE
• Introduction
1. Modes de fissuration élémentaires
2. Champs mécaniques en pointe de fissure
3. Comportement en zone d’élaboration
RMP RUPTURE - Février 2023
Différentes approchesINTRODUCTION
numériques pour la fissuration
Approche locale Approche globale
Basée sur l'étude physique des processus de Basée sur l'extension de la mécanique des
dégradation aux échelle fines qui influe sur milieux continus pour la prise en compte
le comportement du matériau des singularités des fissures macroscopiques
EF, X-FEM, BEM, Meshless
Mécanique de
Mécanique de la rupture
l'endommagement Calcul de FICs Griffith 1924, Irwin 1926
Lemaitre et Chaboche 74
Loi de propagation en fatigue Aspects numériques
Aspects numériques
- souvent linéaire
- fortes non-linéarités
maillage, levelset, - description de la discontinuité
- finesse des discrétisations
transfert de champs - singularité des champs, contact,
- phénomène de localisation
propagation géométrique
Méthodes non-locales
Modèles de zone cohésive Bazant 76, Peerlings 96
Approches variationnelles
Dugdale 60, Barenblatt 59 Francfort et Marigo 98
Transition endommagement/rupture
4RMP RUPTURE - Février 2023
ZONES DE FISSURATION
Elaboration
Zone 2
Zone 1
Fissure
Singulière Zone 3
Lointaine
RMP RUPTURE - Février 2023
MODES DE FISSURATION
x2
x1
x3
Mode I Mode II Mode III
Ouverture de Cisaillement dans le Cisaillement hors
fissure suivant x2 plan de la fissure plan de la fissure
suivant x1 suivant x3
RMP RUPTURE - Février 2023
CHAMPS MECANIQUES ISOTROPES EN POINTE DE
FISSURE
1. Notation en pointe de fissure
r
q
Fissure
Pointe de la fissure
RMP RUPTURE - Février 2023
2. Facteurs d’Intensité de Contrainte (FIC)
FIC en mode I
Contrainte en mode I FIC en mode II
K1(σ ) # θ 3⋅ θ & θ K 2(σ ) # θ 3⋅ θ & θ
σ 11 = ⋅ %1− sin ⋅ sin ( ⋅ cos + ⋅ %2 + cos ⋅ cos ( ⋅ sin
2⋅π ⋅r $ 2 2 ' 2 2⋅π ⋅r $ 2 2 ' 2
Contrainte en mode II
K1(σ ) # θ 3⋅ θ & θ K 2(σ ) θ 3⋅ θ θ
σ 22 = ⋅ %1+ sin ⋅ sin ( ⋅ cos − ⋅ cos ⋅ cos ⋅ sin
2⋅π ⋅r $ 2 2 ' 2 2⋅π ⋅r 2 2 2
Contrainte en mode de cisaillement
K1(σ ) θ 3⋅ θ θ K 2(σ ) # θ 3⋅ θ & θ
σ 12 = ⋅ cos ⋅ cos ⋅ sin − ⋅ %1− sin ⋅ sin ( ⋅ cos
2⋅π ⋅r 2 2 2 2⋅π ⋅r $ 2 2 ' 2
RMP RUPTURE - Février 2023
Déplacement en mode I
FIC en mode I
FIC en mode II
K1(σ ) r # θ 3⋅ θ & K 2(σ ) r # θ 3⋅ θ &
u1 = ⋅ ⋅ %( 2 ⋅ k −1) ⋅ cos − cos ( − ⋅ ⋅ %( 2 ⋅ k + 3) ⋅ sin + sin (
4⋅π 2⋅π $ 2 2 ' 4⋅π 2⋅π $ 2 2 '
Déplacement en mode II
K1(σ ) r # θ 3⋅ θ & K 2(σ ) r # θ 3⋅ θ &
u2 = ⋅ ⋅ %( 2 ⋅ k +1) ⋅ sin − sin ( + ⋅ ⋅ %( 2 ⋅ k + 3) ⋅ cos + cos (
4⋅π 2⋅π $ 2 2 ' 4⋅π 2⋅π $ 2 2 '
Expressions de k
Contrainte plane Déformation plane
3− υ
k = 3− 4 ⋅ υ k=
1+ υ
RMP RUPTURE - Février 2023
3. Facteurs d’Intensité de Déformation (FID)
Déplacements en milieu dissipatif dans l’espace LC
Fonctions angulaires
r "
uα ( r, θ , t ) = ⋅ #gαβ (θ ) ⋅ Cβ (t ) + hαβ (θ ) ⋅ Dβ (t )$%
2.π
FID dans l’espace temporel
" t
1 % ∂K (σ ) t "
k ( t − τ ) % ∂K (σ )
β β
Cβ ( t ) = ∫ $$ '
' ⋅ d τ et Dβ ( t ) = ∫ $
$ '
' ⋅ dτ
−∞ # 2.µ ( t − τ ) & ∂τ −∞ # 2.µ ( t − τ ) & ∂τ
Expressions
E CP DP
µ= 3− υ
2 ⋅ (1+ υ ) k = 3− 4 ⋅ υ k=
1+ υ
RMP RUPTURE - Février 2023
4. Facteurs d’Intensité d’Ouverture de fissure (FIO)
Schématisation des lèvres de la fissure
UMS
q
M UMI
Expression des déplacements
! ! !
[u ] (r, t ) = u (r, t, θ = +π ) − u (r, t, θ = −π )
RMP RUPTURE - Février 2023
Ouverture de la fissure en mode I et II
FIO en mode I FIO en mode II
1 1
" r % (ε )
2 " r % 2
(ε )
[ 1]( ) $
u r, t = ' ⋅ K 2 ( t ) et [ 2 ]( ) $
u r, t = ' ⋅ K1 (t )
#2⋅π & #2⋅π &
K β(ε ) (t ) = 2 ⋅ "#Cβ (t ) + Dβ (t )$%
Relation entre FIO et FIC dans l’espace temporel
(ε )
t
∂K β(σ ) & k +1 )
K ( t ) = ∫ Cβ ( t − τ ) ⋅
β dτ avec Cβ (t ) = ( + (t )
−∞ ∂τ ' µ *
RMP RUPTURE - Février 2023
CHAMPS MECANIQUES ORTHOTROPES EN POINTE
DE FISSURE
1. Notation en pointe de fissure
r
q
Fissure
Pointe de la fissure
RMP RUPTURE - Février 2023
2. Facteurs d’Intensité de Contrainte (FIC)
FIC en mode I
Contrainte en mode I FIC en mode II
K1(σ ) * s ⋅s $ s s '- K (σ ) * 1 $ s2 s 2 '-
σ 11 = ⋅ ℜe , 1 2 ⋅ && 2 − 1 ))/ + 2
⋅ ℜe , ⋅ && 2 − 1 ))/
2⋅π ⋅r ,+ s1 − s2 % ρ 2 ρ1 (/. 2 ⋅ π ⋅ r ,+ s1 − s2 % ρ 2 ρ1 (/.
Contrainte en mode II
K1(σ ) * 1 $ s s2 '- K 2(σ ) * 1 $ 1 1 '-
1
σ 22 = ⋅ ℜe , ⋅ && − ))/ + ⋅ ℜe , ⋅ && − ))/
2⋅π ⋅r ,+ s1 − s2 % ρ 2 ρ1 (/. 2 ⋅ π ⋅ r ,+ s1 − s2 % ρ 2 ρ1 (/.
Contrainte en mode de cisaillement
K1(σ ) * s ⋅s $ 1 1 '- K (σ ) * 1 $ s s '-
σ 12 = ⋅ ℜe , 1 2 ⋅ && − ))/ + 2
⋅ ℜe , ⋅ && 1 − 2 ))/
2⋅π ⋅r ,+ s1 − s2 % ρ1 ρ 2 (/. 2 ⋅ π ⋅ r ,+ s1 − s2 % ρ 2 ρ1 (/.
RMP RUPTURE - Février 2023
Déplacement en mode I
r $ 1 '
u1 = 2 ⋅ K ⋅ (σ )
1
2⋅π
⋅ ℜe &
% s1 − s2
(
⋅ p2 ⋅ s1 ⋅ ρ 2 − p1 ⋅ s2 ⋅ ρ1 )
(
)
r $ 1 '
(σ )
+ 2 ⋅ K2 ⋅
2⋅π
⋅ ℜe &
% s1 − s2
(
⋅ p2 ⋅ ρ 2 − p1 ⋅ ρ1 )
(
)
Déplacement en mode II
r $ 1 '
u2 = 2 ⋅ K ⋅ (σ )
1
2⋅π
⋅ ℜe &
% s1 − s2
(
⋅ q2 ⋅ s1 ⋅ ρ 2 − q1 ⋅ s2 ⋅ ρ1 )
(
)
r $ 1 '
(σ )
+ 2 ⋅ K2 ⋅
2⋅π
⋅ ℜe &
% s1 − s2
(
⋅ q2 ⋅ ρ 2 − q1 ⋅ ρ1 )
(
)
Expressions des polynômes caractéristiques
S22
ρ j = cos (θ ) + i ⋅ s j ⋅ sin (θ ) avec j ∈ {1;2} et p j = S11 ⋅ s 2j + S12 et q j = + S12
sj
Sij : racines du polynôme S11 ⋅ sβ4 + ( 2 ⋅ S12 + S33 ) ⋅ sβ2 + S22 = 0
RMP RUPTURE - Février 2023
4. Facteurs d’Intensité d’Ouverture de fissure (FIO)
Ouverture de la fissure en mode I et II dans l’espace LC
∗
" i ⋅ ( p2 − p1 ) $ r
"#u1 $%( r, p) = 4 ⋅ ℜe )
∗ (σ )
* ( p) ⋅ K 2 ( p) ⋅
# s1 − s2 % 2⋅π
∗
" i ⋅ ( q2 ⋅ s1 − q1 ⋅ s2 ) $ r
"#u2 $%( r, p) = 4 ⋅ ℜe )
∗
p ⋅ K (σ )
*( ) 1 ( ) p ⋅
# s1 − s2 % 2⋅π
FIO en mode I et II dans l’espace LC
∗
(ε )∗
% i ⋅ ( q2 ⋅ s1 − q1 ⋅ s2 ) ( ∗(σ )
K1 ( p) = 4 ⋅ ℜe ' * ( p) ⋅ K1 ( p)
& s1 − s2 )
∗
(ε )∗
% i ⋅ ( p2 − p1 ) ( ∗(σ )
K 2 ( p) = 4 ⋅ ℜe ' * ( ) 2 ( p)
p ⋅ K
& s1 − s2 )
RMP RUPTURE - Février 2023
FIO dans l’espace temporel
t (σ )
∂K
K β(ε ) = ∫ Cβ (t − τ ) ⋅ β dτ
0 ∂τ
Expressions des complaisances élastiques réduites
Mode I Mode II
$ i ⋅ ( q2 ⋅ s1 − q1 ⋅ s2 ) ' $ i ⋅ ( p2 − p1 ) '
C1 (t) = 4 ⋅ ℜe & )(t ) C2 (t) = 4 ⋅ ℜe & )(t )
% s1 − s2 ( % s1 − s2 (
RMP RUPTURE - Février 2023
COMPORTEMENT EN ZONE d’ELABORATION
1. Zone d’élaboration viscoélastique
Zone d’élaboration Majoration des contraintes
viscoélastique selon Dugdale
s ab s ab
s mb s mb
S b1 S b2
r r
a ab
1
a ab
2
1 1 " K (σ ) %
Taille de la zone d’élaboration α β = α β = ⋅ $$ β( β ) ''
π #σm &
RMP RUPTURE - Février 2023
2. Comportement mécanique en zone dégradée
s 22
s f =sm
a0 a
x1 a
Champ de contrainte transversale en pointe de fissure
FIC Longueur de la zone d’élaboration
1 2
π " Kβ %
1 (σ )
(σ ) !8$ 2
K β = # & ⋅ σ m( β ) ⋅ (α β ) 2 α β = ⋅ $$ ( β ) ''
"π % 8 #σm &
RMP RUPTURE - Février 2023
3. Zone d’élaboration en mode mixte
Hypothèses sur les FIC selon Dugdale
(σ ) (σ )
(m) K K
σ 22 = Α⋅ I et σ 12(m) = Α ⋅ II
π ⋅α π ⋅α
Longueur de la zone d’élaboration en mode mixte
(démontrer par chaque groupe)
2
# K (σ )
K (σ ) &
π
α = Α ' ⋅ % (m)
I II
+ (m) ( avec Α ' =
$ σ 12 σ 22 ' 8
Formuler toutes les hypothèses
nécessaires à la démonstration
RMP RUPTURE - Février 2023
UFR – MR2
COURS DE MECANIQUE DE LA
RUPTURE
CHAPITRE 1:
GENERALITES SUR LA RUPTURE
Rostand MOUTOU PITTI
MCF - HDR
RMP RUPTURE - Février 2023