Laboratoire de Simulation des Matériaux Cours Milieux Continus

Département des Matériaux 3ème Semestre

EPFL J.-M Drezet J. Valloton, K. Shahim & M. Salgado

Corrigé No 11  05.12.08
Equations Constitutives
Elasticité
1. Cylindre creux sous contrainte
On écrit d'abord le champ de déplacement :
u~r = r0 (h(r0 ) − αz0 )
α 2 µ
u~z = r + α(1 + )z02 + βz0
2 0 λ
1. Puisqu'il s'agit d'un solide dans le régime élastique, il ne peut que subir
de petites déformations. Le tenseur correspondant ε donc a la forme :
ε = 12 (Grad u + Grad u T )
 
r0 ∂h(r 0)
∂r + h(r0 ) − αz0 0 0
ε= 0 h(r0 ) − αz0 0
 

µ
0 0 2αz0 (1 + λ ) + β
Noter que l'on respecte les coordonnées cylindriques, donc il y a bien
un terme qui entre en jeu pour le gradient du champ de déplacement
dans la composante θθ. Le tenseur de déformations ayant pour com-
posantes non nulles seulement celles qui font la diagonale, les vecteurs
propres correspondent bien au repère orthonormé (e~ , e~ , e~ ). Désor-
mais, le tenseur de contraintes se calcule en utilisant la relation (9.8)
r θ z

du polycopié pour le cas isotrope. Cela permet d'écrire :

σrr = (λ + 2µ)εrr + λ(εθθ + εzz )
σθθ = (λ + 2µ)εθθ + λ(εrr + εzz )
σzz = (λ + 2µ)εzz + λ(εθθ + εrr )
en remplaçant avec les composantes du tenseur de déformations on
obtient, après simplication :
(λ + 2µ)(r0 ∂h(r0)
 
∂r + h(r0 )) + λh(r0 ) + βλ
λ(r0 ∂h(r0)
 
∂r + h(r0 )) + (λ + 2µ)h(r0 ) + βλ
 
 

σ= λ(r0 ∂h(r0) 3λ+2µ
∂r + h(r0 )) + λh(r0 ) + 2αµz0 λ + (λ + 2µ)β 

0
 
 
 
 0 
0
2. L'équation d'équilibre mécanique div σ + ρ~g = 0 nous permettra de
trouver le champ de forces volumiques. En réalité, le champ est donné
par div σ, qu'on résout en respectant encore les coordonnées cylin-
driques (voir relation 3.36 pour la divergence d'un champ tensoriel)
1
 Laboratoire de Simulation des Matériaux Cours Milieux Continus
Département des Matériaux 3ème Semestre
EPFL J.-M Drezet J. Valloton, K. Shahim & M. Salgado

σrr −σθθ
 
∂σrr
∂r + r
divσ =  0
 

∂σzz
∂z
Pour la composante en e~ on obtient :
r
!
∂ 2 h(r0 ) ∂h(r0 ) ∂h(r0 )
(λ + 2µ) r0 2
+2 +λ +
∂r ∂r ∂r
1 ∂h(r0 ) ∂h(r0 )
   
(λ + 2µ) r0 + h + λh + βλ − λ(r0 + h) − (λ + 2µ)h(r0 ) − βλ
r ∂r ∂r
qu'on simplie en : !
∂ 2 h(r0 ) ∂h(r0 )
(λ + 2µ) r0 2
+3
∂r ∂r
Pour la composante en e~ z
3λ + 2µ
2αµ
λ
3. Pour que les forces de volume soient eectivement les forces de pesan-
teur, on pose simplement l'équation d'équilibre mécanique. Dès lors, on
peut résoudre l'équation de deuxième dégrée dans la première égalité
et retrouver α à partir de la deuxième : !
∂ 2 h(r0 ) ∂h(r0 )
(λ + 2µ) r0 +3 =0
∂r2 ∂r

On pose y = ∂h(r
∂r
)
et on obtient après double intégration :
0

C
h(r0 ) = +k
2r02
Pour calculer α, on pose
3λ + 2µ
2αµ = −ρg
λ
et on obtient :
−ρgλ
α=
2µ(3λ + 2µ)
Le tenseur de contraintes s'écrit alors :
λ(2k + β) − rµ2
 

σ =  λ(2k + β) + µ( r12 + 2k) 

 
2λk − ρgz + (λ + 2µ)β

2. Contraintes planes
En contraintes planes, on a σ = σyz = σxz = 0 . Le tenseur des contraintes
est donc :
zz

2
 Laboratoire de Simulation des Matériaux Cours Milieux Continus
Département des Matériaux 3ème Semestre
EPFL J.-M Drezet J. Valloton, K. Shahim & M. Salgado
 
σxx σxy 0
σ =  σxy σyy 0 
 
0 0 0

Pour trouver le tenseur de déformation, on peut inverser la relation :

σij = λT r(ε)δij + 2µεij

A partir de cette relation, on résout d'abord pour ε ij

δij λ
εij = 2µ − 2µ T rε

Du même, il est aisé de montrer que :

T r(σ) = (3λ + 2µ)T r(ε)

et donc, on obtient :
δij λ T rσ
εij = 2µ − 2µ 3λ+2µ

Et, en remplaçant λ et µ par leur expression en fonction de E et ν

εij = − Eν T r(σ)δ − ij + 1+ν
E σij

Dans le cas des contraintes planes, on a donc :

 ε= 
− ν (σxx + σyy ) + 1+ν
E σxx
1+ν
E σxy 0
 E 1+ν
E σxy − E (σyy + σxx ) + 1+ν
ν
E σyy 0

 
0 0 − Eν (σxx + σyy )

On voit donc que, en contraintes planes, les déformations ne sont pas néces-
sairement planes.
3. Déformations planes
En déformations planes, on a ε = εyz = εxz = 0 . Le tenseur de déformation
est donc de la forme :
zz

 
εxx εxy 0
ε =  εxy εyy 0 
 
0 0 0

De plus, le cas de l'élasticité isotrope, on a :

σij = λT r(ε)δij + 2µεij

3
Laboratoire de Simulation des Matériaux Cours Milieux Continus
Département des Matériaux 3ème Semestre
EPFL J.-M Drezet J. Valloton, K. Shahim & M. Salgado

avec λ = C = (1 + ν)(1
12
Eν
et µ = C 44 =
E

Le tenseru des contraintes est donc :

− 2ν) 2(1 + ν)

 
(λ + 2µ)εxx + λεyy 2µεxy 0
σ= 2µεxy (λ + 2µ)εyy + λεxx 0
 

0 0 λ(εxx + εyy )

On voit donc que, en déformation plane, les contraintes ne sont pas néces-
sairement planes. Il existe une contrainte σ normale au plan.
zz