www.almohandiss.

com

Laboratoire de Mtallurgie Physique Cours Millieux Continus

Dpartement des Matriaux 3me Semestre

EPFL M. Rappaz M. Felberbaum M. Salgado

Srie No 5  17.10.07
Tenseur des contraintes II

1.
Soit dans un repre orthonorm ( . On considre le champ de contraintes
Contraintes Principales et Cercle de Mohr
e1 , e2 , e3 )
homgne dni par (en MPa)

10 0 0
= 0 5

15

(1)
0 15 5
1. Dterminer les contraintes principales ( , , ) dans l'ordre d-
croissant ainsi que les directions principales associes ( , , )
I II III
nP I nP II nP III

2. Vrier que les trois invariants du tenseur restent inchangs

3. Dterminer la matrice de passage [A] du repre ( , , ) au repre
propre ( , , ) telle que r = [A]r
e1 e2 e3
eI eII eIII
0

4. Montrer que [A]..[A] est eectivement diagonal.

T

5. Pour une direction quelconque, dterminer les composantes normale,

t , et tangentielle, t , de la densit de force . Dessiner le cercle de
n

Mohr.
n t

6. Dterminer l'ellipse des contraintes dans les trois plans ( , ), ( , eII eIII )
et ( , ).
eI eII
eI eIII

7. Dterminer les parties hydrostatique et dviatorique du tenseur dans

( , , ).
e1 e2 e3

2.
Un tat de contraintes inhomogne dni sur le cube unitaire x [0, 1],
Condition d'quilibre

y [0, 1] et z [0, 1] (en m) est caractris par le tenseur suivant :

1 3y 3z
= 3y

xy xyz

(2)
3z yz zz
est en MPa, si x, y et z en m.
1. Trouver les deux composantes manquantes ( ) et ( ) sur la base
des indications suivantes :
yz r zz r

A. Le tenseur satisfait la condition d'quilibre div = 0 en tout point.

B. La trace du tenseur est nulle en tout point.
C. Le long de l'axe x, le tenseur = x ). yz
2

2. Calculez les forces globales (en N) exerces sur chacune des faces de
ce cube. Vriez que la somme de ces forces est nulle. Pourquoi a-t-on
cela?
1

www.almohandiss.com
 www.almohandiss.com

Laboratoire de Mtallurgie Physique Cours Millieux Continus

Dpartement des Matriaux 3me Semestre
EPFL J.-M Drezet, J. Valloton, K. Shahim & M. Salgado

Corrig No 5  17.10.07
Tenseur des contraintes II

1. Contraintes Principales et Cercle de Mohr

Pour trouver les contraintes principales il sut de rsoudre le dterminant
de la matrice

10 p 5 3 0
~ = 5 3

0 p 0

(1)
0 0 10 p

 10 
 p 5 3 0 

DET =  5 3 0 p 0  = (p2 10p )(10p )(10p )(5 3)2
 
10 p 
 

Les valeurs des contraintes principales sont , et

 0 0
p1 = 15 pII = 5 pIII =
10.
De toute vidence, par la forme du tenseur des contraintes, le valeur
correspond au valeur propre align sur l'axe z. En eet, la rotation du tenseur
pIII

pour retrouver les valeurs propres se fait autour de l'axe z, et c'est pourquoi
on ne les met pas dans un ordre dcroissant.
Les directions principales associes ces vecteurs s'obtiennent en appli-
quant l'quation (4.27) du polycopier. Ainsi, on trouve que pour , n
( 3/2, 1/2, 0) ; pour , n (1/2, . 3/2, 0) ; et pour ,n
pI pI

(0, 0, 1).
pII pII pIII pIII

On vrie que les invariants du tenseur restent inchangs :

I () = 10 + 0 + 10 = 20 = I (p ) = 15 5+ 10 = 20
II () = (10)(0) + (10)(10) + (0)(10) 02 (5 3)2 02 = 25
II (p ) = (15)(5) + (15)(10) + (5)(10) 0 = 25


 10 
 5 3 0 

III () = Det() =  5 3 0 0  = 750
 
 
 0 0 10 

 15 0 0

III (p ) = Det(p ) =  0 5 0 = 750


 0 0 10
Dsormais il nous faut trouver la matrice de rotation pour passer de .
Pour ce faire, l'angle de rotation autour de z doit tre calcul. Puisque on
p

est sur a, on va dj s'assurer que les vecteurs propres sont orthonormales

avec la relation suivante :
|n ||n | cos = n~ n~
1 2 1 2 (3)
1

www.almohandiss.com
 www.almohandiss.com

Laboratoire de Mtallurgie Physique Cours Millieux Continus

Dpartement des Matriaux 3me Semestre
EPFL J.-M Drezet, J. Valloton, K. Shahim & M. Salgado

Ayant dmontr que c'est le cas, on applique la mme quation aux vecteurs
n et n pour trouver = 30 degrs. Noter que pour la transformation de
repre r~ = [A]~r, [A] n'est rien d'autre que le tenseur qui dcrit l'orientation
1 pI
0

des vecteurs propres.

Pour une direction ~n quelconque ou ~n = n(x, y, z), on peut trouver la com-
posante normale et tangentielle de la densit de force en oprant ainsi :
tn = ~t ~n = ~n ~n ~n = 20x2 5y 2 + 10z 2 + 15xy + 5yz + 30xz;

t = t tn
NOTE : On a intrt avoir une normale sur laquelle on aimerait calculer
ces composants avant de passer au dveloppement de t dans les quations
prcdents.

Le cercle de Mohr et les ellipses des contraintes dans les trois plans sont
dessins ci-dessous :

Finalement, la partie hydrostatique d'un tenseur ( ) est dnie par 1/3T r()I
alors que
la partie dviatorique
du tenseur est dnie comme = :
M
D M
25/3 0 0
M = 0 25/3 0

et
0 0 25/3

5/3 5 3 0

D = 5 3 25/3 0

0 0 5/3

www.almohandiss.com
 www.almohandiss.com
Laboratoire de simulation des matriaux Cours Milieux Continus
Dpartement des matriaux 3me semestre
EPFL M.Rappaz/M. Felberbaum, M. Salgado
Exercice 1 : Equilibre mcanique sur un cube unitaire

 1 3y 3z


 =  3y xy  yz

3z   zz

 yz

On sait que div() = 0, le flux est donc nul !

 


 0
r

div( ) =  x + yz
= 0
 z

  yz  zz

 +

 y z

On sait que la trace du tenseur est nulle, on a donc zz = -1 - xy

On trouve alors que yz = -z x + A(x, y) (daprs la 2me quation ci-dessus)
La 3me quation donne la constante A(x, y) qui ne doit pas tre fonction de y (sinon la drive
de yz par rapport y est non nulle, ce qui est impossible vu que la drive de zz par rapport
z est nulle) ; en consquence, la constante A (x, y) = x2 vu que sur laxe x, yz = x2
On a alors le tenseur suivant :

 1 3y 3z 


 =  3y xy x ( x  z)

3z x ( x  z) (1 xy)



Pour calculer les forces globales sur lensemble du cube unitaire, il faut appliquer la relation
t =  n et intgrer sur chaque surface du cube entre 0 et 1. On a donc 6 intgrales faire.

Normale // Ox :
 1 3y 3z  1  1 3y 3z   1
1 1
 
1 1
 
   3y 0 0
 0 dydz +    3y y (1 z)
0 dydz =
0 0   0  
3z 0 1   0  3z 1(1 z) (11y)  0
0

 1   1 
 
1 1 1 1
=   3y dydz +    3y dydz = 0
0  0 
3z
0 0
 3z 

www.almohandiss.com
 www.almohandiss.com
Laboratoire de simulation des matriaux Cours Milieux Continus
Dpartement des matriaux 3me semestre
EPFL M.Rappaz/M. Felberbaum, M. Salgado
Normale // Oy :
 1 0 3z   0   1 3 3z   0
1 1
  1 1
 
  0 0 x ( ) 
x  z
1 dxdz +   3 x x ( x  z)
 1 dxdz =
0 0   0  
3z x ( x  z) 1   0  0
 3z x ( x  z) (1 x)   0
 0   3 
1 1
 1 1

=   0 dxdz +    x dxdz
0 0  0 0 
x ( x  z)  x ( x  z)

Normale // Oz :
 1 3y 0  0 1 3y 3   0
1 1
 
1 1
 
   3y xy x 2
 0 dxdy +    3y xy x ( x 1)
 0 dxdy =
0 0   0  
 0 x (1 xy) 1  3 x ( x 1) (1 xy)  1
2 0

 0   3 
 
1 1 1 1
=    x 2 dxdy +    x ( x 1) dxdy
0  0 
 (1 xy)
0 0
(1 xy)

On se retrouve donc avec :

 0   3   0   3 
1 1

1 1

1 1


1 1


   0
dxdz +    x
dxdz +    x 2
dxdy +    x ( x 1)
dxdy = 0
0 0 
0 0 
0 
0 

x ( x  z)  x ( x  z) (1 xy)

0 0
 (1 xy)

On voit en effet bien que les intgrales sannulent, et ceci est vident par le thorme de la
divergence !

www.almohandiss.com