Mécanique linéaire de la rupture

Introduction aux éléments finis étendus

Daniel Weisz-Patrault

Ecole des Ponts

18 Mars 2020

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 1 / 71

 Objectifs

Rappels des Elements Finis

Principes de base de la X-FEM
Culture générale sur les méthodes numériques
Critère énergétique de propagation
Méthodes locales (endommagement, modèles à gradient) pas
abordées

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 2 / 71

Plan de le séance

1 Rappels

2 Eléments finis étendus : X-FEM

3 Propagation

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 3 / 71

Plan de le séance

1 Rappels

2 Eléments finis étendus : X-FEM

3 Propagation

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 4 / 71

Rappels

Formulation variationnelle
Méthode des éléments finis

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 5 / 71

Formulation variationnelle

Ωt ⊂ R3 : domaine à l’instant t
∂Ω : frontière à l’instant t
∂ΩT : frontière où l’on impose le vecteur contrainte
∂Ωu : frontière où l’on impose les déplacement

∂Ω = ∂ΩT ∪ ∂Ωu and ∂ΩT ∩ ∂Ωu = ∅

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 6 / 71

Formulation variationnelle

Espace des vitesses généralisées cinématiquement admissibles

n o
C = V : (x, t) ∈ Ωt × R+ 7→ V (x, t) / ∀x ∈ ∂Ωu , u(x, t) = ud (x, t)

Espace des vitesses virtuelles cinématiquement admissibles

C ∗ = {V ∗ : x ∈ Ωt 7→ V ∗ (x) / ∀x ∈ ∂Ωu , u∗ (x, t) = 0}

Espace des vitesses rigidifiantes

∗
= V ∗R : x ∈ Ωt 7→ V T + ω.x, ∀V T ∈ R3 / ∀ω ∈ Mas

CR 3

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 7 / 71

Formulation variationnelle

Puissance des efforts extérieurs

Z Z
∗ ∗
P V E(V ) = ρf .V dΩ + T .V ∗ dS
Ωt ∂ΩT

Puissance des efforts intérieurs

Z
∗
F .V ∗ − σ : ∇ [V ∗ ] dΩ


P V I(V ) =
Ωt

Puissance des efforts d’accélération

Z
∗
P V A(V ) = ργ.V ∗ dΩ
Ωt

γ champ d’accélération réel.

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 8 / 71

Formulation variationnelle
Condition de cohérence

La puissance intérieure d’un mouvement rigide est nulle

∀V T ∈ R3 ∀ω ∈ Mas
3 / P V I(V T + ω.x) = 0
D’où
∀V T ∈ R3 ∀ω ∈ Mas 3
Z



F . V T + ω.x − σ : ω dΩ = 0
Ωt
D’où Z
F = 0 et σ : ωdΩ = 0
Ωt
Tenseur anti-symmétrique/symétrique

F = 0 et σ ∈ Ms3

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 9 / 71

Formulation variationnelle
Puisque σ est symmétrique on introduit le taux de déformation

1
d∗ (V ∗ ) = ∇ [V ∗ ] + t ∇ [V ∗ ]


2
Puissance des efforts intérieurs
Z
∗
P V I(V ) = − σ : d∗ (V ∗ ) dΩ
Ωt

Lien entre les contraintes et les variables cinématiques

Comportement : σ = f (d, · · · )

On cherche V ∈ C telle que pour σ vérifiant le comportement :

∀V ∗ ∈ C ∗ P V I(V ∗ ) + P V E(V ∗ ) = P V A(V ∗ )

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 10 / 71

Formulation variationnelle
Principe des puissances virtuelles : V ∈ C telle que ∀V ∗ ∈ C ∗
Z Z Z Z
σ : d∗ (V ∗ ) dΩ = ρf .V ∗ dΩ + T .V ∗ dS − ργ.V ∗ dΩ
Ωt Ωt ∂ΩT Ωt

Théorème de la divergence Ωt div σ.V ∗ dΩ = ΩT σ.n .V ∗ dS

  RR

Z Z Z
∗ ∗ ∗
div σ .V ∗ dΩ

 
σ : d (V ) dΩ = σ.n .V dS −
Ωt ∂ΩT Ωt

Equations fortes du problème

  

 div σ = ργ − ρf sur Ωt


σ.n = T sur ∂ΩT


 u = ud

sur ∂Ωu

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 11 / 71

Formulation variationnelle
Réciproque : on multiplie les équations fortes par V ∗ ∈ C ∗
Intégration sur Ωt
Z Z
div σ .V ∗ dΩ = ρ γ − f .V ∗ dΩ
 

Ωt Ωt

Théorème de la divergence (V ∗ = 0 sur ∂Ωu )

Z Z Z
div σ .V ∗ dΩ = σ.n .V ∗ dS − σ : d∗ (V ∗ ) dΩ
 

Ωt ∂ΩT | {z } Ωt
T

On cherche donc V ∈ C telle que ∀V ∗ ∈ C ∗

Z Z Z Z
∗ ∗ ∗ ∗
σ : d (V ) dΩ = ρf .V dΩ + T .V dS − ργ.V ∗ dΩ
Ωt Ωt ∂ΩT Ωt

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 12 / 71

Rappels

Formulation variationnelle
Méthode des éléments finis

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 13 / 71

Méthode des éléments finis

Recherche de solutions approchées

Méthode des éléments finis
(i) Formulation faible (ou variationnelle)
(ii) Méthode d’approximation (de Galerkin)
(iii) Interpolation sur un maillage

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 14 / 71

Méthode des éléments finis

Formulation faible
V ∈ C solution si et seulement si ∀V ∗ ∈ C ∗
Z Z Z Z
∗ ∗ ∗ ∗
σ : d (V ) dΩ = ρf .V dΩ + T .V dS − ργ.V ∗ dΩ
Ωt Ωt ∂ΩT Ωt

Méthode d’approximation de Galerkin

Problème : C de dimension infinie
Projection de C sur un espace de dimension finie

CA = p (C)

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 15 / 71

Méthode des éléments finis

Maillage
Espace particulier CA
Découpage de Ω = N
S e
e=1 Ωe
Fonctions continues et polynômiales par morceaux de degrée n
 
f ∈ C 0 Ωt , R3 / ∀(x, y, z) ∈ Ωe 



 

 

h i   

n n n
Pn (Ωe )1≤e≤Ne = XXX

 f (x, y, z) = fijk xi y j z k 

 


 i=0 j=0 k=0 


i+j+k≤n

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 16 / 71

Méthode des éléments finis

Maillage
Espace des vitesses généralisées approchées
Ne
( )
[
d
CA = V ∈ Pn : (x, t) ∈ Ωe × R+ 7→ V (x, t) / u(x, t) = u (x, t)
e=1

Espace des vitesses virtuelles approchées

Ne
( )
[
∗
CA = V ∗ ∈ Pn : x ∈ 7→ V ∗ (x) / ∀x ∈ ∂Ωu , u∗ (x, t) = 0
e=1

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 17 / 71

Méthode des éléments finis

Maillage

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 18 / 71

Méthode des éléments finis
Maillage
Dimension : dim [CA ] = N

N = Nombre de noeuds total − Nombre de noeuds bloqués

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 19 / 71

Méthode des éléments finis

Explicitation d’une base

Fonctions de Pn valant 1 sur un noeud et 0 sur tous les autres


∀1 ≤ i ≤ N, ϕi ∈ Pn et ϕi xj = δij

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 20 / 71

Méthode des éléments finis
Explicitation d’une base
Exemple 1D, n = 1 et n = 2

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 21 / 71

Méthode des éléments finis
Fonction de forme
Restriction : Nie = ϕi |Ωe

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 22 / 71

Méthode des éléments finis
Notion d’élément de référence
Coordonnées réelles (x, y, z), de calcul (ζ, η, ξ)

Interpolation de la géométrie (sur un élément)

3

 X
x(ζ, η) = xi N
ei (ζ, η)




i=1
 3
X
 y(ζ, η) = yi N


 ei (ζ, η)
i=1

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 23 / 71

Méthode des éléments finis
Notion d’élément de référence
Coordonnées réelles (x, y, z), de calcul (ζ, η, ξ)

Interpolation des déplacements (sur un élément)

3
X 3
X
u(ζ, η) = ei (ζ, η) u(x, y) =
ui N ui Nie (x, y)
i=1 i=1
3
X 3
X
v(ζ, η) = vi N
ei (ζ, η) v(x, y) = vi Nie (x, y)
i=1 i=1

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 24 / 71

Méthode des éléments finis

Ecriture matricielle par élément (triangulaire linéaire)

Déplacements u|Ωe (x, y, z)

u1
 

 v1 

 w1 
N1e 0 0 N2e 0 0 N3e 0
   
0 
 u2 

 0 N1 0e 0 N2 0e 0 N3e 0 . v2 
 
0 0 N1e 0 0 N2e 0 0 N3e  w2 
| {z } 
 u3


Ne (x, y, z)  
 v3 
w3
| {z }
ue

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 25 / 71

Méthode des éléments finis

Ecriture matricielle par élément (triangulaire linéaire)

Champ de vitesse virtuel V ∗|Ω (x, y, z)
e

U1∗
 
 V1∗ 
W1∗
 
 
N1e 0 0 N2e 0 0 N3e 0 U2∗
   
0  
V2∗
 
 0 N1e 0 0 N2 0e 0 N3e 0  . 
 
0 e
0 N1 0 0 N2 0 e 0 N3e  W2∗ 
| {z }  U3∗


Ne (x, y, z)  
 V3∗ 
W3∗
| {z }
V ∗e

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 26 / 71

Méthode des éléments finis

Ecriture matricielle par élément (triangulaire linéaire)

Taux de déformation virtuel d∗| (x, y, z)
Ωe

1  h i t h i ∗
d∗| (x, y, z) = ∇ N e + ∇ N e .V e
|2
Ωe
{z }
Be (x, y, z)

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 27 / 71

Méthode des éléments finis
Discrétisation du PPV en quasi-statique
V ∈ CA solution si et seulement si pour σ = f (d, · · · ), ∀V ∗ ∈ CA
∗

Z Z Z
σ : d∗ (V ∗ ) dΩ = ρf .V ∗ dΩ + T .V ∗ dS
Ωt Ωt ∂ΩT

Linéarité de l’intégrale
Ne Z
X Ne Z
X Z
∗ ∗ ∗
σ : d (V ) dΩ = ρf .V dΩ + T .V ∗ dS
e=1 Ωe e=1 Ωe Ωe ∩∂ΩT

Projection sur l’espace d’approximation de dimension finie

Ne Z
X Ne Z
X Z
∗ ∗
σ : Be .V dΩ =e ρf .Ne .V e dΩ + T .Ne .V ∗e dS
e=1 Ωe e=1 Ωe Ωe ∩∂ΩT

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 28 / 71

Méthode des éléments finis

Discrétisation du PPV en quasi-statique

Projection sur l’espace d’approximation de dimension finie
Ne
X Z Ne
X Z Z 
t ∗ t t ∗ t t
V e. Be : σ dΩ = V e. ρ Ne .f dΩ + Ne .T dS
e=1 Ωe e=1 Ωe
| {z }
fe

Cas linéaire σ = C : ε = C : Be .ue

Ne
X Z  XNe
t
V ∗e . t
Be : C : Be dΩ .ue = t ∗
V e .f e
e=1 Ωe e=1
| {z }
Ke

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 29 / 71

Méthode des éléments finis

Assemblage
Somme de produits matriciels
Ne
X Ne
X
t
V ∗e .Ke .ue = t
V ∗e .f e
e=1 e=1

Assemblage : [V ∗ ] valeurs aux noeuds sans répétition

t
[V ∗ ] . K . [u] = t [V ∗ ] . f
   

 
Problème : K pas inversible.

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 30 / 71

Méthode des éléments finis

Réduction du système
 
On cherche K R inversible
Lignes : équations 0 0 = 00 car [V ∗ ] contient des 0 (sur ∂Ωu )
Colonnes : [u] connu sur les composantes où [V ∗ ] est nul
Ces composantes connues sont mises dans le vecteur force

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 31 / 71

Méthode des éléments finis

Inversion du système
Système après réduction ∀ [V ∗ ]R
t
[V ∗ ]R . K R . [u]R = t [V ∗ ]R . f R
   

Simplification (car vrai pour tout [V ∗ ]R ne s’annulant pas)

 −1  
[u]R = K R . f R

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 32 / 71

Méthode des éléments finis
Démarche comportement non-linéaire
Procédure incrémentale (incréments n)
Incrément de chargement ∆f e imposé
On cherche l’incrément solution correspondant ∆uek
Construction et assemblage d’un résidu
Ne
X Z   
t ∗ t
V e. Be : σ (uek + ∆uek ) − f ek + ∆f e dΩ = 0
e=1 Ωe
| {z }
Rek (∆uek )

Assemblage du résidu
Annulation du résidu par Newton-Raphson
∂σ
Fait apparaître un module de rigidité tangeante ∂ε
A déterminer par un algorithme local
Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 33 / 71
Méthode des éléments finis

Newton-Raphson
Comment passer de l’incrément k − 1 à k ?
On connaît ∆f e
On cherche l’incrément solution correspondant [∆uk ]
Assemblage du résidu
t
[V ∗ ] . [Rk ] ([∆u]) = 0

Réduction pour enlever les equations 0 = 0

[Rk ]R ([∆uk ]) = 0

Chercher un point d’annulation d’une fonction vectorielle

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 34 / 71

Méthode des éléments finis

Newton-Raphson
[Rk] R [Rk] R

Converge rapidement Ne congerve jamais

[Δuk] (1) [Δuk] (1)

[Δuk] [Δuk]
[Δuk] (0) [Δuk] (2) [Δuk] (0)

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 35 / 71

Méthode des éléments finis
Newton-Raphson

hInitialisation
i : [∆uk ](0) , puis annulation de la tangente en
(i−1)
∆uk

∂ [Rk ]R (∆u)  
[Rk ]R ([∆uk ](i−1) ) + . [∆uk ](i) − [∆uk ](i−1) = 0
∂∆u
 −1
(i) ∂ [Rk ]R (∆u)
[∆uk ] =− . [Rk ]R ([∆uk ](i−1) ) + [∆uk ](i−1)
∂∆u
| {z }
K −1
T

Résidu élémentaire
Ne
X Z   
t ∗ t
V e. Be : σ (uek + ∆uek ) − f ek + ∆f e dΩ = 0
e=1 Ωe
| {z }
e e
Daniel Weisz-Patrault (Ecole des Ponts)
R k
RUPTU
(∆uk) 18 Mars 2020 36 / 71
Méthode des éléments finis
Newton-Raphson
Rigidité tangente élémentaire

∂σ
Z  
e t e e
KT = Be : (u + ∆uk ) dΩ
Ωe ∂∆uek k
Z " #
e t
∂σ ∂ε
KT = Be : : dΩ
Ωe ∂ε ∂∆uek
Approximation EF
ε = Be . (uek + ∆uek )

Rigidité tangente élémentaire

" #
∂σ
Z
K eT = t
Be : : Be dΩ
Ωe ∂ε

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 37 / 71

Méthode des éléments finis

Module tangent equivalent

Cf cours MPM élasto-plasticité HPP
   
∂f ∂f
 C : ∂σ ⊗ C : ∂σ 
σ̇ = C − 
 ∂f : C : ∂f − ∂f
 : ε̇
 (1)
∂σ ∂σ ∂pcum
| {z }
C eq

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 38 / 71

Méthode des éléments finis

Algorithme local d’intégration de la plasticité (ex simple)

Objectif savoir dans quelles mailles on plastifie
Comment mettre à jour les variables plastiques
Incrément k − 1 on connait tout
Incrément k on connaît ∆f e et [∆uk ](i) (itération de NP)
Loi d’écoulement de von Mises écrouissage isotrope

3 s
ε̇p = λ
2 σY (pcum )

Différences finis implicites

εpk − εpk−1 3 sk
= λ
∆t 2 σY (pcum,k )

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 39 / 71

Méthode des éléments finis
Algorithme local d’intégration de la plasticité (ex simple)
Notation ∆λ = λ ∆t
Déformation plastique

3 sk
εpk = εpk−1 + ∆λ
2 σY (pcum,k )

Relation (cf cours elasto-plasticité HPP) : ṗcum λ

pcum,k = pcum,k−1 + ∆λ

Evolution elasto-plastique : ∆λ ≥ 0, fk ≤ 0, ∆λ fk = 0

⇒ ∆λ = 0 et ∆εp = 0

 fk < 0 alors élastique
 ( ∂fk p
∂σ : σ̇k ≤ 0 ⇒ ∆λ = 0 et ∆ε = 0
f
 k = 0 alors ∂fk p ∂fk
∂σ : σ̇k > 0 ⇒ ∆λ > 0 et ∆ε = λ ∂σ


Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 40 / 71

Méthode des éléments finis

Algorithme local d’intégration de la plasticité (ex simple)

(i)
On pose ∆ε = B : ∆ue,k (connu itération de NP)
e

εk = εk−1 + ∆ε

Notation ek déviateur de εk
Test élastique (∆λ = 0, εpk = εpk−1 )

str
k
= 2µ(ek − εpk−1 )

r
3 tr tr
fktr = s : sk − σY (pcum,k−1 )
2 k

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 41 / 71

Méthode des éléments finis
Algorithme local d’intégration de la plasticité (ex simple)
Premier cas : fktr ≤ 0 on garde le test élastique
Deuxième cas : fktr > 0 correction plastique ∆λ > 0
r
3 tr tr
fk = s : sk − σY (pcum,k ) = 0
2 k

Réécriture
sk = 2µ(ek − εpk )
= 2µ(ek − εpk−1 ) − 2µ(εpk − εpk−1 )
| {z }
∆εp

D’où
3 sk
sk = str
k
− 2µ∆λ
2 σY (pcum,k )
Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 42 / 71
Méthode des éléments finis

Algorithme local d’intégration de la plasticité (ex simple)

Cas général fk = 0 se résout par Newton-Raphson
Cas écrouissage linéaire

σY (pcum,k ) = σ0 (1 + α pcum,k )

D’où  
3 ∆λ
str
k
= sk 1 + 2µ
2 σY (pcum,k )
r r
3 tr tr 3 3 ∆λ
sk : sk = sk : sk 1 + 2µ
2 | 2 {z } 2 σY (pcum,k )
σY (pcum,k )

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 43 / 71

Méthode des éléments finis
Algorithme local d’intégration de la plasticité (ex simple)
Cas écrouissage linéaire
str s
 
q k
3 tr tr
= k
σY (pcum,k )+2µ 23 ∆λ
1 + 2µ 32 σY (p∆λ
cum,k )
s
2 k
:s
k
sk
=
σY (pcum,k )

Ecriture de fk = 0
q
fk = 32 sk : sk − σY (pcum,k ) = 0
q
= 32 strk
: str
k
− 2µ 32 ∆λ − σ0 (1 + α pcum,k−1 ) − σ0 ∆λ = 0
fktr
fktr − (3µ + σ0 )∆λ = 0 ⇒ ∆λ =
3µ + σ0

Mise à jour des variables plastique et des contraintes

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 44 / 71
Méthode des éléments finis
Résumé
Incrément k-1
p ,u ext
σk-1 ,εk-1 , ε k-1 k-1 ,pcum,k-1
, , fk-1
Connus

k=k+1
f kext=f k-1
ext
+Δf ext
Newton-Raphson
Tangente: (i)=(i+1)
Δuk(i) itération NP
si >tol si <tol
Algorithme local Calcul
σk ,εk , ε kp , uk ,pcum,k, , fkext Résidu en k,(i)

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 45 / 71

Plan de le séance

1 Rappels

2 Eléments finis étendus : X-FEM

3 Propagation

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 46 / 71

Eléments finis étendus : X-FEM

Cas 1D : fissuration aux noeuds

Cas 1D : fissuration entre les noeuds
Extension en 2D/3D : fissuration aux noeuds
Extension en 2D/3D : cas général
Introduction aux fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 47 / 71

Cas 1D : fissuration aux noeuds

Déplacement dans les 2 éléments

u(x) = u1 ϕ1 (x) + u− − + +
2 ϕ2 (x) + u2 ϕ2 (x) + u3 ϕ3 (x)

Degrés de liberté : déplacement moyen et saut

u− +
2 + u2 u− +
2 − u2
hui = [u] =
2 2
Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 48 / 71
Cas 1D : fissuration aux noeuds

Fonction d’Heaviside généralisée

Déplacement dans les 2 éléments

u(x) = u1 ϕ1 (x) + hui ϕ2 (x) + u3 ϕ3 (x) + [u] ϕ2 (x)H(x)

Avec
ϕ2 (x) = ϕ− +
2 (x) + ϕ2 (x)

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 49 / 71

Eléments finis étendus : X-FEM

Cas 1D : fissuration aux noeuds

Cas 1D : fissuration entre les noeuds
Extension en 2D/3D : fissuration aux noeuds
Extension en 2D/3D : cas général
Introduction aux fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 50 / 71

Cas 1D : fissuration entre les noeuds

Déplacement avec enrichissement

u = u1 ϕ1 + u2 ϕ2 + u3 ϕ3 + u4 ϕ4 + u
e2 ϕ2 H + u
e3 ϕ3 H

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 51 / 71

Eléments finis étendus : X-FEM

Cas 1D : fissuration aux noeuds

Cas 1D : fissuration entre les noeuds
Extension en 2D/3D : fissuration aux noeuds
Extension en 2D/3D : cas général
Introduction aux fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 52 / 71

Extension en 2D/3D : fissuration aux noeuds

Déplacements
  10  
u X ui
= ϕi (x)
v vi
i=1

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 53 / 71

Extension en 2D/3D : fissuration aux noeuds

Degrés de liberté : déplacement moyen et saut

u9 + u10 u9 − u10
hui = [u] =
2 2
v9 + v10 v9 − v10
hvi = [v] =
2 2
Enrichissement
  8      
u X ui hui [u]
= ϕi (x) + ϕ11 + ϕ11 H(x, y)
v vi hvi [v]
i=1

Avec
ϕ11 (x, y) = ϕ9 (x, y) + ϕ10 (x, y)

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 54 / 71

Extension en 2D/3D : fissuration aux noeuds

Equivalence entre dédoubler les noeuds et enrichir les fonctions

de forme

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 55 / 71

Eléments finis étendus : X-FEM

Cas 1D : fissuration aux noeuds

Cas 1D : fissuration entre les noeuds
Extension en 2D/3D : fissuration aux noeuds
Extension en 2D/3D : cas général
Introduction aux fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 56 / 71

Extension en 2D/3D : cas général

I indices non encerclés

J indices encerclés
  X  X u 
u ui ej
= ϕi + ϕj H(x, y)
v vi vej
i∈I j∈J

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 57 / 71

Extension en 2D/3D : cas général

Enrichissement avec la fonction H sur les noeuds encerclés

Problème de la longueur de fissure
Fissure s’arrête soit en A soit en C

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 58 / 71

Extension en 2D/3D : cas général
Solution la plus singulière en plan : Westergaard
           
KI θ θ 3θ KII θ θ 3θ
σxx = √ cos 1 − sin sin −√ sin 2 + cos cos
2πr 2 2 2 2πr 2 2 2
           
KI θ θ 3θ KII θ θ 3θ
σxy = √ sin cos cos +√ cos 1 − sin sin
2πr 2 2 2 2πr 2 2 2
           
KI θ θ 3θ KII θ θ 3θ
σyy = √ cos 1 + sin sin +√ sin cos cos
2πr 2 2 2 2πr 2 2 2
r   r  
KI r θ KII r θ
ux = cos (κ − cos (θ)) + sin (2 + κ + cos (θ))
2µ 2π 2 2µ 2π 2
r   r  
KI r θ KII r θ
uy = sin (κ − cos (θ)) + cos (2 − κ − cos (θ))
2µ 2π 2 2µ 2π 2

Fonctions de forme additionnelles : repère local en tête de fissure

√
√
√ √
r sin 2θ , r cos 2θ , r sin 2θ cos(θ), r cos 2θ cos(θ)



(Fl )1≤l≤4 =

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 59 / 71

Extension en 2D/3D : cas général

I tous les noeuds, J encerclés, K encadrés

Enrichissement
  X  ui  X u  4  
u ej X X u
bk,l
= ϕi + ϕj H + ϕk Fl
v vi vej vbk,l
i∈I j∈J k∈K l=1

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 60 / 71

Extension en 2D/3D : cas général

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 61 / 71

Eléments finis étendus : X-FEM

Cas 1D : fissuration aux noeuds

Cas 1D : fissuration entre les noeuds
Extension en 2D/3D : fissuration aux noeuds
Extension en 2D/3D : cas général
Introduction aux fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 62 / 71

Introduction aux fonctions de niveaux

Description géométrique de la fissure

Essentiel pour évaluer H
Décorrélé du maillage
Description explicite difficile
Automatisation du maillage de la fissure
Forme de fissure complexe : maillage distordu
Branchement de fissures
Description implicite : fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 63 / 71

Introduction aux fonctions de niveaux
Projection sur Γ(t)

xp (t) = argmin kx∗ − xk

x∗ ∈Γ(t)

Fonction de niveau : distance (signée) entre Γ(t) et x

 
ψ(x, t) = sign (x − xp (t)).n x − xp (t)

Par définition
Γ(t) = {x ∈ Ωt / ψ(x, t) = 0}
φ fonction de niveau orthogonale à ψ
Description implicite
Interpolation sur la base élément finis
X
ψ(x, y, t) = ψi (t)ϕi (x, y)
i∈I

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 64 / 71

Introduction aux fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 65 / 71

Introduction aux fonctions de niveaux

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 66 / 71

Plan de le séance

1 Rappels

2 Eléments finis étendus : X-FEM

3 Propagation

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 67 / 71

Propagation

Bifurcation
Critère énergétique de Francfort et Marigo

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 68 / 71

Bifurcation
Branche bifurquée peut avoir une longueur infinitésimale
Calcul différentiel
Directions courantes de bifurcation
1) Contrainte circonférentielle maximale
∂σθθ ∂ 2 σθθ
=0 <0
∂θ ∂2θ
2) Taux de restitution d’énergie maximal
∂G ∂2G
=0 <0
∂θ ∂2θ
3) Critère mode I. α angle de bifurcation, l longueur de bifurcation

(KI∗ (α, l), KII

∗ e ∗ (α), K
e ∗ (α))
(
(α, l)) → (K I II
l→0
∗
KII (α) = 0
e

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 69 / 71

Propagation

Bifurcation
Critère énergétique de Francfort et Marigo

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 70 / 71

Critère énergétique de Francfort et Marigo
Instabilité de la propagation
On ne modélise pas les effets dynamiques entre deux instants
Propagation de taille finie

F [α, l] = E − W + D
Z
1
E= σ : ε dΩ : énergie élastique
2 Ωt
Z Z
W= ρf .u dΩ + T .u dS : potentiel des efforts extérieurs
Ωt ∂ΩT
D = l Gc énergie dissipée en fissurant
Minimisation
(α, l) = argmin F [α∗ , l∗ ]
(α∗ ,l∗ )

D : coût pour aller vers un état énergétique plus faible

Daniel Weisz-Patrault (Ecole des Ponts) RUPTU 18 Mars 2020 71 / 71