29.03.2015 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!