Méthode XFEM pour la Modélisation de grandes ... - Cast3M - CEA
Méthode XFEM pour la Modélisation de grandes ... - Cast3M - CEA
Méthode XFEM pour la Modélisation de grandes ... - Cast3M - CEA
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Chapitre I. Rappels et Cadre <strong>de</strong> <strong>la</strong> Recherche<br />
3.1 Approches globales<br />
Les approches globales sont basées sur le calcul d’intégrales <strong>de</strong> contour et sont utilisées <strong>pour</strong><br />
modéliser l’amorçage puis <strong>la</strong> propagation d’une fissure. Pour les matériaux ductiles, <strong>la</strong><br />
p<strong>la</strong>sticité en pointe <strong>de</strong> fissure n’est plus confinée et l’application <strong>de</strong> ces approches n’est<br />
possible que dans un cadre restreint.<br />
3.1.1 Rappel : Mécanique linéaire <strong>de</strong> <strong>la</strong> rupture<br />
La mécanique linéaire <strong>de</strong> <strong>la</strong> rupture est basée sur l'hypothèse d'é<strong>la</strong>sticité <strong>de</strong> <strong>la</strong> totalité <strong>de</strong> <strong>la</strong><br />
structure fissurée dont on étudie le comportement. Dans ce cadre, Irwin [Irwin1957] a montré<br />
<strong>la</strong> re<strong>la</strong>tion <strong>de</strong> proportionnalité entre le facteur d'intensité <strong>de</strong> contraintes K I et le premier terme<br />
du développement limité <strong>de</strong>s champs <strong>de</strong> contraintes en pointe <strong>de</strong> fissure. Le champ <strong>de</strong><br />
contraintes en un point M du voisinage <strong>de</strong> <strong>la</strong> pointe <strong>de</strong> fissure O est donné par :<br />
Figure I-2 - Repère <strong>de</strong> pointe <strong>de</strong> fissure<br />
σ<br />
σ<br />
σ<br />
xx<br />
yy<br />
xy<br />
K<br />
θ ⎛ θ θ ⎞<br />
⎜<br />
⎟<br />
2 ⎝ 2 2 ⎠<br />
θ ⎛ θ θ ⎞<br />
⎜<br />
⎟<br />
2 ⎝ 2 2 ⎠<br />
θ ⎛ θ θ ⎞<br />
⎜<br />
⎟<br />
2 ⎝ 2 2 ⎠<br />
I<br />
= cos 1−<br />
sin sin 3<br />
2πr<br />
K<br />
I<br />
= cos 1+<br />
sin sin 3<br />
2πr<br />
K<br />
I<br />
= cos sin cos3<br />
2πr<br />
(I-21)<br />
Où r et θ sont respectivement le rayon et l’angle donnant <strong>la</strong> position du point M dans le repère<br />
<strong>de</strong> pointe <strong>de</strong> fissure défini sur <strong>la</strong> figure I-2.<br />
Griffith [Griffith1920] fut le premier a donné une interprétation énergétique <strong>de</strong> <strong>la</strong> rupture<br />
fragile. Il introduisit le taux <strong>de</strong> restitution d'énergie G qui correspond à l'énergie é<strong>la</strong>stique<br />
surfacique libérée <strong>pour</strong> une avancée infinitésimale <strong>de</strong> fissure da dans un soli<strong>de</strong> é<strong>la</strong>stique<br />
linéaire.<br />
Une autre manière <strong>de</strong> calculer <strong>la</strong> valeur non signée <strong>de</strong> l'énergie é<strong>la</strong>stique libérée par l'avancée<br />
<strong>de</strong> fissure est <strong>de</strong> calculer le travail ∆W nécessaire <strong>pour</strong> refermer <strong>la</strong> fissure sur <strong>la</strong> longueur da.<br />
Soit <strong>pour</strong> une structure p<strong>la</strong>ne d'épaisseur B :<br />
Cette expression permet <strong>de</strong> plus <strong>de</strong> relier G et K I :<br />
da<br />
∆ W = B∫ σ . u dx = G Bda<br />
(I-22)<br />
0<br />
yy<br />
y<br />
2<br />
K<br />
I<br />
G = (1 −ν<br />
²) en déformations p<strong>la</strong>nes ; G =<br />
E<br />
K<br />
E<br />
2<br />
I<br />
en contraintes p<strong>la</strong>nes<br />
(I-23)<br />
3.1.2 Bases <strong>de</strong> <strong>la</strong> mécanique <strong>de</strong> <strong>la</strong> rupture non-linéaire :<br />
l’intégrale J<br />
Le calcul <strong>de</strong> l'intégrale <strong>de</strong> contour J a été proposé par Rice [Rice1968], son résultat donne le<br />
travail d'ouverture <strong>de</strong>s lèvres lors d'une extension <strong>de</strong> fissure dans un matériau é<strong>la</strong>stique<br />
linéaire ou non. Elle se calcule sur un contour fermé Г entourant <strong>la</strong> pointe <strong>de</strong> fissure, d'où son<br />
écriture dans le cas d'une structure p<strong>la</strong>ne :<br />
⎛ ∂u<br />
⎞<br />
J = ∫ ⎜w dy −T<br />
⋅ ⋅ds⎟<br />
(I-24)<br />
Γ ∂x<br />
⎝<br />
⎠<br />
11