Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Exercices Corrigés - Guilhem Mollon PDF

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 57

INTRODUCTION A LA MECANIQUE DES

M
MIILIEUX CONTINUS

Exercices corrigés

Guilhem MOLLON

Polytech Grenoble
Départem
ment Géotechnique, Troisième année
e
Edition 1, 2012-2013
V1.07
Mécanique des Milieux Contin
nus Polyteech Grenoble, Geo3

1. Calcul tenssoriel

Exercice A. Montrer que la symétrie est une propriété tensorielle, c'eest-à-dire que si un
tenseur  est symétrique (   ) dans une base orthonormée donné /  0
 
  1,
alors cette propriété est éggalement vraie dans toute autre base orthonormée
o / 
  
   
 .

Exercice B. Montrer que l'aantisymétrie est également une propriété teensorielle.


Exercice C. Montrer qu'un tenseur  quelconque peut toujours se déécomposer en une
partie symétrique et une parrtie antisymétrique.

Exercice D. Soit  un tenseeur symétrique et  un tenseur antisymétrrique, montrer que


l'on a toujours    .

Exercice E. Soit  un tenseuur symétrique et  un tenseur quelconquee, montrer que l'on


a toujours :       , oùù   est la partie symétrique de .

Exercice F. Soit  un cham  2,  


mp vectoriel. Montrer que l'on a toujours $%(&'

Exercice G. Soit  un cham


mp scalaire. Montrer que l'on a toujours 
&  2,  
 2(3&4$
' 

Exercice H. Soit la base currviligne polaire (


)  
*  + , telle que :
,

      
 !      
"    

1. Calculer la surface extérrieure d'une sphère de rayon , sachant qu'un élément de


surface infinitésimal peut s'éécrire $    $$ .
2. Calculer l'intégrale du champ
c ) sur cette sphère, sachaant que le vecteur
  
élémentaire dans la directionn radiale vaut 
)   - .  
 5 .  
6

2
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Corrigés

Exercice A.
Soit un tenseur  symétrique. Dans la base /  0,  ,
   1 on peut donc écrire
 =  . Soit une base orthonormée quelconque /′ = ′  , ′
 , ′
  différente de /.
On appelle 7 la matrice de passage de / à /′. Si on appelle ′ les termes de la
matrice de  dans la base /′, on a d'après le cours (en notation d'Einstein) :
′ = 78 79 89

Or la matrice  est symétrique, on a donc 89 = 98 . On écrit donc :

′ = 78 79 98

Les indices  et ; du second membre sont muets, on pourrait donc les remplacer par
n'importe quelle lettre, et on peut aussi les intervertir :
′ = 79 78 89 = 78 79 89

On en déduit que ′ = ′ , ce qu'il fallait démontrer.

Exercice B.
Soit un tenseur  antisymétrique. Dans la base / = ,  
,
   1 on peut donc écrire
 = − . Soit une base orthonormée quelconque /′ = ′  , ′
 , ′
  différente de /.
On appelle 7 la matrice de passage de / à /′. Si on appelle ′ les termes de la
matrice de  dans la base /′, on a d'après le cours (en notation d'Einstein) :
′ = 78 79 89

Or la matrice  est symétrique, on a donc 89 = −98 . On écrit donc :

′ = −78 79 98

Les indices  et ;du second membre sont muets, on pourrait donc les remplacer par
n'importe quelle lettre, et on peut aussi les intervertir :
′ = −79 78 89 = −78 79 89

On en déduit que ′ = −′ , ce qu'il fallait démontrer.


3
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Exercice C.
On reprend les formules du cours. Soient les tenseurs   et   donnés par les formules
suivantes :


  = (   @ 


  = ( −  @ 

On va démontrer que   est symétrique et que   est antisymétrique. Pour cela on se


place dans une base / = ,
 ,
  1, pour laquelle le tenseur  s'exprime sous la forme
d'une matrice de terme générale  . Explicitions les termes des matrices de   et   :

 
 = (    = (    = 
 
 
 = ( −   = − ( −   = −
 

On a donc démontré que   est symétrique et   antisymétrique.

Exercice D.
Dans une base orthonormée / = ,
 
,
  1, on peut écrire que  =  et également
que  = − . On en déduit que, comme pour tout tenseur antisymétrique, les termes
diagonaux de la matrice de  sont nuls dans toute base. D'après le cours, le produit
doublement contracté de ces deux tenseurs est un scalaire égal à :

:  =  

Développons cette notation d'Einstein sous forme explicite :

:  =        


       
       

On sait que l'on a  =  =  = 0, et par ailleurs :

 =  ;  =  ;  = 


 = − ;  = − ;  = −

On en déduit directement :  = 0.

4
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Exercice E.
Le résultat découle directement de celui de l'exercice précédent et de la distributivité
de l'opérateur produit doublement contracté :

:  = : (      = :    :   = :  

Exercice F.
Posons le résultat intermédiaire / = &'  2. Dans une base orthonormée donnée
/ = ,
 ,
  1, on a par définition :

D D

C D D G
D D
 = B  −  F
/ B D F
B  D F
D D

A D D EH

Par ailleurs, on a par définition : $%2/ = IH


I-
J IH
I-
IH
 K  L.
I- J K L

On en déduit que $%2&'


 2 vaut :

D D D D D D D D D


 2 =
$%2&' M − N M − N M − N
D D D D D D D D D

D   D   D   D   D   D  
 2 =
$%2&' −  −  −
D D D D D D D D D D D D

L'ordre des dérivations partielles successives d'une fonction de plusieurs variables est
quelconque, on peut donc en déduire directement :2$%2&'
 2 = 0.

Exercice G.
On se place également dans une base orthonormée / = ,  ,
  1. Soit le résultat
 2 donné par la formule du cours :
intermédiaire / = 3&4$
D
CD G
D F
 = B
/ BD F
B F
D
AD EH

5
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Par ailleurs, on peut écrire, toujours dans la base / ∶


D/ D/

C D D G
D/ D/
 = B  −  F
 2/
&' B D F
B  D F
D/ D/

A D D EH

On en déduit :
D D D D

CD D D D G
D D D D F
 2 = B
 23&4$
&' BD D − D D F
B    F
D D D D

AD D D D EH

Pour la même raison que dans l'exercice précédent, on peut donc écrire directement la
formule classique &'  2 = 0
 23&4$ .

Exercice H.
1. Un élément de surface $ de la sphère de rayon  délimité par deux secteurs
d'angles infinitésimaux $ et $ s'écrit :
$ =   $$

On cherche à calculer  = PQ8RèTU $. En paramétrant la surface en fonction de  et ,


on peut expliciter cette intégrale :
W W
=V V   $$
*XY +XY

Le rayon  est indépendant de  et , et le terme  est indépendant de , on peut


donc écrire :
W W

= V  ZV $ [ $
*XY +XY

W W
= V 
 \ ]W
Y $ = ^ V 
 $ = ^ \− ]WY
*XY *XY

6
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Finalement, on obtient le résultat classique  = ^.

2. On doit calculer :
2
W W
 = V  ∙ $
T =V V  ∙ T ∙   $$
Q8RèTU *XY +XY

Avec : 
) =  -   5   
6

On est en coordonnées curvilignes, donc on ne peut pas sortir le vecteur 


) car il n'est
pas indépendant du point d'intégration (du point de la sphère de surface $). En
revanche, les trois vecteurs de la base cartésienne ont cette propriété, et peuvent être
sortis de l'intégrale. On peut donc remplacer T par son expression, et écrire dans la
base / = , , 1 :
W W
V V  ∙   ∙   $$
C *XY +XY G
B W W F

= B V V  ∙   ∙   $$ F
B *XY +XY F
B W W F
V V  ∙  ∙   $$
A *XY +XY EH

W W W W
V V   $$ V   ZV
$[ $
C *XY +XY G C *XY +XY G
B W W F B W W F
 =   B V V   $$ F =   BV   ZV  $ [ $ F
B *XY +XY F B *XY +XY F
B W W F B W W F
V V    $$ V    ZV $ [ $
A *XY +XY EH A *XY +XY EH

Du fait de la périodicité des fonctions trigonométriques, on a :


W W
V $ = V $ = 0
+XY +XY

Par ailleurs, on a :
W
V $ = ^
+XY

7
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Donc :
0
0
 = ^ ` W a
V    $
*XY H

L
La primitive de la fonction    est la fonction bcQ *, donc :

0 0
C 0 G 0
 = ^ B   W F = `^ a
d f
 Y  H
A EH

8
Mécanique des Milieux Continus Polytech Grenoble, Geo3

2. Cinématique et déformations

Problème. On considère un mouvement défini dans la base /  0,


 ,
  1 par sa
représentation lagrangienne (g est une constante positive) :
 = i g' 1 − i g' 1
h = i g' 1  i g' 1#
 = i

1. Calculer le tenseur gradient j, le tenseur des dilatations k, et le tenseur des
déformations l de ce mouvement au point i et à l'instant '.
2. A quelle classe particulière ce mouvement appartient-il ?
3. Pour un instant ' donné, calculer la dilatation en un point i et dans une direction 
$i

4. Pour un instant ' donné, calculer le glissement en un point i et pour deux directions
orthogonales 
$i et 
$i′.

5. On considère un milieu animé de ce mouvement, muni d'une masse volumique


homogène mY à l'instant 'Y = 0. Calculer le jacobien de la transformation, ainsi que la
masse volumique du milieu à l'instant '.
6. Calculer le champ de vitesse  (i , '
n et le champ d'accélération o(i , ' en
coordonnées lagrangiennes.
7. Exprimer les coordonnées initiales à partir des coordonnées actuelles. Calculer le
champ de vitesse n, '1 et le champ d'accélération o, '1 en coordonnées
eulériennes.
8. Calculer les tenseurs des taux de déformations eulériens q , ' 1
p et des taux de
rotation Ωq, '1.
9. On définit les coordonnées polaires lagrangiennes , Θ, 1 par le changement de
variables i, i, 1 =  ∙ Θ,  ∙ Θ, 1 et les coordonnées eulériennes
&, , x 1 par le changement de variables  ,  , x 1 = & ∙ , & ∙ , x 1.
Expliciter les fonctions iT et i+ définissant une nouvelle représentation lagrangienne
du mouvement de la forme :
& = iT , Θ,   , t1 #
v
 = i+ , Θ,   , t1
9
Mécanique des Milieux Continus Polytech Grenoble, Geo3

10. On définit une base curviligne polaire / = ,


T ,
+   1, telle que :

T = 
  
 #
v
+ = −
  
 

On note alors le champ de vitesse dans cette nouvelle base :

 , ' 1 = nT &, ,  , ' 1


n T  n+ &, ,  , ' 1+  n &, ,  , ' 1


Indiquer l'expression des composantes polaires nT , n+ et n du champ de vitesse


eulérien.

11. Calculer l'accélération centrifuge oT &, , ' 1 et l'accélération tangentielle o+ &, , ' 1
du mouvement étudié.

12. Définir les trajectoires associées à ce mouvement

Examen partiel : Etude cinématique d'un tourbillon. On considère un mouvement,


défini dans la base orthonormée / = ,
 ,
  1 par la représentation eulérienne
suivante :

n = − ,  1
h n =  ,  1 #
n = 0

Dans cette représentation,  est une fonction scalaire des deux coordonnées  et  ,
y
définie par l'expression  ,  1 =    1z où 7 est une constante positive.
W
Ce mouvement est donc uniquement défini pour {   { > 0.

1. De quel type de mouvement s'agit-il ?

 2n
2. Montrer que la matrice du tenseur gradient de vitesse } = 3&4$  dans la base
/ = ,
 ,
  1 s'exprime par :

 ,  1    −  0


} =  ∙~  
0
    −  − 
0 0 0H

q , ' 1 et du
3. En déduire immédiatement les matrices du tenseur des taux de rotation Ω
q , ' 1.
tenseur des taux de déformation eulériens p

4. On considère un milieu continu de masse volumique mY à l'état initial (' = 0) et


animé de ce mouvement. Calculer la divergence du champ de vitesse, et en déduire la
masse volumique en tout point et à tout instant ' > 0.
10
Mécanique des Milieux Continus Polytech Grenoble, Geo3

5. Montrer que l'accélération o  ,  ,  , ' 1 en coordonnées eulériennes cartésiennes


(c'est-à-dire dans la base / = ,
 ,
  1) s'exprime par :

−

o , ' 1 =   ,  1 ∙ d− f
0 H

6. On définit les coordonnées polaires lagrangiennes , Θ,   1 par le changement de


variables i , i ,   1 =  ∙ Θ,  ∙ Θ,   1 et les coordonnées eulériennes
&, , x 1 par le changement de variables  ,  , x 1 = & ∙ , & ∙ , x 1.

On définit également une base curviligne polaire / = ,


T ,
+   1, telle que :

T = 
  
 #
v
+ = −
  
 

Calculer l'accélération o &, ,  , ' 1 en coordonnées eulériennes polaires.

7. Exprimer le champ de vitesse en coordonnées eulériennes polaires. En déduire que


les particules ont des trajectoires circulaires autour de l'origine, de vitesse angulaire :

7
n+ =
^&

Donner sans calcul l'expression de la représentation polaire lagrangienne du


mouvement sous la forme & = iT , Θ,  , t1 et  = i+ , Θ,  , t1.

8. Démontrer que la représentation lagrangienne du mouvement dans la base


cartésienne sous la forme  = (i , ' s'exprime par :

sg' 1 −sg' 1 0
 = ~ sg' 1 sg' 1 0 i
0 0 

Montrer que, dans cette expression, on a : g = =  ,  1.
)

9. Définir les trajectoires associées à ce mouvement, ainsi que les champs de vitesse et
d'accélération

10. On considère le point de vecteur position  = & ∙  . Calculer la matrice du tenseur
des taux de déformations eulériens en ce point. De quel type de déformation s'agit-il ?
En déduire les valeurs propres et les vecteurs propres du tenseur des taux de
déformations eulériens en ce point.

11
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Corrigés

Problème.
I-
1. Le terme général de la matrice du tenseur gradient dans la base / est j  I , donc

on a :
g' −g' 0

j  ~ g' g' 0
0 0 H

Pour le tenseur des dilatations, on a k = j8 j8 en notation d'Einstein, donc :

  g'   g' g' ∙ g' − g' ∙ g' 0


k = ~g' ∙ g' − g' ∙ g'   g'   g' 0 =  
0 0 H

Enfin, on sait que l =  (k −  , et on a donc l = 0.

2. Le tenseur des déformations est nul, on est donc en présence d'un mouvement
rigidifiant.

3. Puisque k =  , on peut dire que toute direction est direction principale. La dilatation
dans une direction quelconque 
$i vaut donc ($i  = k = 

4. Pour la même raison, le glissement entre deux directions orthogonales quelconques


 et 
$i  =  =0
 , $i′
$i′ vaut : o $i
 

5. Le jacobien de la transformation est le déterminant de j . On a donc :


g' −g' 0
 = $'j =  g' g' 0 =   g'   g' = 
0 0 

Par conséquent la masse volumique du milieu est constante dans le temps et en tout
point.

12
Mécanique des Milieux Continus Polytech Grenoble, Geo3

6. Le champ de vitesse n(i, ' s'obtient par dérivation matérielle du champ de


vecteurs exprimant les positions actuelles des particules : n(i, ' = -,1

. On a :

 = i g' 1 − i g' 1
h = i g' 1  i g' 1#
 = i

Donc :
p
n = = g ∙ (−i g' 1 − i g' 1
 p'
p
n = = g ∙ (i g' 1 − i g' 1 #
 p'
 p
 n = =0
p'
 ,1
De même, le champ d'accélération s'obtient par o(i, ' =  , donc :

pn
o = = g ∙ (−i g' 1  i g' 1
 p'
pn
o = = g ∙ (−i g' 1 − i g' 1#
 p'
 pn
 o  = =0
p'
Ces champs s'expriment en fonction de i, et sont donc bien en coordonnées
lagrangiennes.

7. D'après l'énoncé, les coordonnées actuelles s'obtiennent à partir des coordonnées


initiales par le système suivant :
 g' −g' 0 i
~  = ~ g' g' 0 ∙ ~i 
 0 0  i

On isole en particulier les deux premières équations de ce système :


 g' −g' i
  =  ∙ 
 g' g' i

On cherche à inverser ce système pour exprimer i en fonction de . On utilise la


formule d'inversion d'une matrice 2*2 :

4  z  $ −
  =  
 $ 4$ −  − 4

13
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Donc :

g' −g' z  g' g' g' −g' @


  =  = 
g' g'  g'   g' −g'
  g' g' g'

Cette matrice est orthogonale, car son inverse est égale à sa transposée. Elle traduit
donc une rotation. On en déduit :

i =  g' 1   g' 1
hi = − g' 1   g' 1#
i = 

Dans l'expression du champ de vitesse et du champ d'accélération calculés à la


question précédente, on peut alors remplacer les coordonnées initiales i par leur
expression en fonction des coordonnées actuelles  . On obtient la formulation
eulérienne suivante :

n = g ∙ −( g' 1   g' 1g' 1 − (− g' 1   g' 1g' 1


#
 1  1  1  1  1  1
 n = g ∙ (  g'    g'  g' − (−  g'    g'  g' 
 n = 0

n = −g
 n = g #
n = 0


o = g ∙ −( g' 1   g' 1g' 1  (− g' 1   g' 1g' 1
 #
 1  1  1  1  1  1
o = g ∙ −(  g'    g'  g' − (−  g'    g'  g' 
 o = 0

o = −g ∙ 
ho = −g ∙  #
o = 0

Ce dernier résultat pouvait aussi s'obtenir en passant par la formule donnant



I 
 2n
directement l'accélération par o =  (3&4$  , avec I = 0
 n  et :
I I

0 −g 0
 2n
3&4$  = ~g 0 0
0 0 0 H

14
Mécanique des Milieux Continus Polytech Grenoble, Geo3

8. On calcule d'abord le tenseur gradient de vitesse } = 3&4$ 2n


, avec la formule de
  = I . On remarque qu'il s'agit d'un gradient eulérien (g minuscule),
 2n
cours (3&4$  I-
donc les dérivations sont effectuées par rapport aux coordonnées actuelles  :

0 −g 0
  
} = 3&4$ 2n = ~g 0 0
0 0 0H

Ce tenseur est antisymétrique, on peut donc directement écrire p q = 0 et Ω


q = . Ces
résultats peuvent être retrouvés par le calcul avec les formules de cours :

 
q=
p (}  }@ 22222'22222Ω
q = (} − }@ 
 

9. En coordonnées polaires, on a les correspondances suivantes :

i =  ∙ Θ#  = & ∙ θ# & = iT , Θ,   , t1 #


v ; v ; v
i =  ∙ Θ  = & ∙ θ  = i+ , Θ,   , t1

On cherche les expressions des fonctions iT et i+ . On peut écrire à partir de ces


expressions :

 
& =    2222222 2222222θ = 2222222 2222222θ =
r r
i  i
 = i  i 2222222 2222222Θ = 2222222 2222222Θ =
 

Or on sait d'après l'énoncé que, dans le repère cartésien, le mouvement s'exprime par :

 = i g' 1 − i g' 1
h = i g' 1  i g' 1#
 = i

On peut remplacer  et  par leurs expressions dans la formules donnant & :

 
& =    = (i g' 1 − i g' 1  (i g' 1  i g' 1

& = i  i = 

On en déduit l'expression de la fonction iT . On effectue la même opération avec


l'expression de θ :

15
Mécanique des Milieux Continus Polytech Grenoble, Geo3

i g' 1 − i g' 1
θ =
r

Dans cette dernière expression, on remplace i , i et & par leur expression polaire :

 ∙ Θg' 1 −  ∙ Θg' 1
θ = = sΘ  g' 1

On en déduit θ = Θ  g' , ce qui nous donne la forme de la fonction i+ . Finalement,


on a :

&= #
¡
θ = Θ  g'

10. On travaille maintenant dans la base curviligne polaire / = ,


T ,
+   1 telle que :

T = 
  
 #
v
+ = −
  


On cherche à exprimer le champ de vitesse eulérien dans cette base. On a démontré


dans la question 7 que le champ de vitesse pouvait s'exprimer par :

 = −g ∙   g ∙ 


n

Dans cette expression, on remplace  et  par leurs expressions en coordonnées


polaires :
n = −g& ∙ θ1 ∙   g& ∙ θ1 ∙ 
donc :
n = g& −θ ∙   1 = g&
  θ ∙  +

Dans la base / = ,


T ,
+   1, la vitesse s'exprime donc par une formule très simple :

0
n = ¢g&£
0 H

11. Dans le même ordre d'idée, on cherche à exprimer la formulation eulérienne du


champ d'accélération dans la base / . On a montré que le champ d'accélération
s'exprimait par :

o = −g  ∙  − g  ∙ 

On remplace  et  par leurs expressions en coordonnées polaires :

16
Mécanique des Milieux Continus Polytech Grenoble, Geo3

o = −g & ∙ θ1 ∙  − g & ∙ θ1 ∙ 

o = −g &θ ∙   1 = −g &T


  θ ∙ 

Dans la base / = ,


T ,
+   1, l'accélération est donc purement centripète :

−g &
o = ¢ 0 £
0 H

12. Les trajectoires des particules matérielles sont des cercles centrés sur l'origine. La
vitesse angulaire est constante, la vitesse est purement orthoradiale (perpendiculaire à
un rayon) et augmente proportionnellement à la distance à l'origine. L'accélération est
purement centripète (dirigée vers l'origine) et est également proportionnelle en norme
à la distance à l'origine.

Examen Partiel.

1. Il s'agit d'un mouvement plan, car la vitesse dans la direction 


 est nulle et les deux
autre composantes n et n sont invariantes par translation dans la direction   . Par
ailleurs la représentation eulérienne du champ de vitesse est constante dans le temps,
donc on est par définition en présence d'un mouvement permanent.
 2n
2. le tenseur gradient de vitesse } = 3&4$  s'exprime dans la base / = ,
 ,
  1
I
par la matrice de terme général } = I- . On travaille d'abord sur la dérivation de la

y
fonction  ,  1 = W
   1z :

D ,  1 7  


=− ∙  = −   ,
  1 ∙
D ^    1   

D ,  1 7  


=− ∙   = − ,  1 ∙ 
D ^    1    

On peut alors calculer les termes non triviaux de la matrice de }, à commencer par les
deux premiers termes diagonaux :
Dn D ,  1  
= − ∙ =  ,  1 ∙ 
D D   

17
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Dn D ,  1  


=  ∙ = − ,  1 ∙ 
D D   

On a également deux termes non-diagonaux assez complexes :

Dn D ,  1 


= − ∙ −  ,  1 =  ,  1 ∙  −  ,  1
D D   

Dn  − 


 1
=   ,  
D   

Dn D ,  1 


=  ∙   ,  1 = − ,  1 ∙    ,  1
D D   

Dn  − 


=  ,  1 
D   

Tous les autres termes de la matrice de } sont nuls car on est en mouvement plan.
Finalement, le tenseur gradient de vitesse s'exprime sous forme matricielle dans la
base / = ,
 ,
  1 par :

 ,  1    −  0


} =  ∙~  
0
    −  − 
0 0 0H

3. On constate que la matrice de } est symétrique. Par conséquent, sa décomposition


q et un tenseur symétrique p
en un tenseur asymétrique Ω q est immédiate :

q = 0
-tenseur des taux de rotation :2Ω
q = }
-tenseur des taux de déformation eulériens : p

4. la masse volumique du milieu à l'état initial (' = 0) est uniforme et vaut mY . Pour
calculer l'évolution de cette masse volumique, on peut évaluer le jacobien de la
transformation en tout point et à tout instant. Dans un problème posé en notation
eulérienne, il est plus simple de passer par la divergence du champ de vitesse :

Dn Dn Dn


 =
$%2n  
D D D

On sait que n est constant et nul, et on a par ailleurs démontré à la question


précédente que :

18
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Dn Dn
=−
D D

Par conséquent on a $%2n = 0. Un mouvement pour lequel la divergence de la vitesse


est nulle en tout point et à tout instant est un mouvement isochore, pour lequel la
masse volumique est constante. Par conséquent, en tout point et à tout instant, on a :
m = mY

Le mouvement est plan, permanent, et isochore.

5. L'accélération o ,  ,  , '1 est la dérivée particulaire de la vitesse, et se calcule


par la formule de cours suivante, qui fait apparaître la dérivée eulérienne :
 Dn
pn 
o , ' 1 = =  2n
 (3&4$ 
 n
p' D'

On est en mouvement permanent, donc la dérivée eulérienne de la vitesse est nulle :



Dn

=0
D'

 2n
Le gradient de vitesse } = 3&4$  a déjà été calculé à la question 2, il suffit donc d'en
faire le produit contracté avec le vecteur vitesse :

 ,  1    −  0 −


 2n
o , ' 1 = (3&4$  =    ∙ ~  −  
 n −  0 ∙  ,  1 ∙ d  f
    
0 H
0 0 0H
 
  ,  1 − ∙   1   ∙  −  1
o , ' 1 =  ∙ ~− ∙  −  1   ∙ −  1
  
0 H

  −
  ,  1 −  −  
 1
o , ' =  1 −
∙ ~−  −   =   ,  ∙ d  f
 
   0 H
0 H

6. On cherche à se placer en coordonnées polaires pour simplifier les expressions. Pour


cela, on utilise les relations :
 = & ∙ θ#
v
 = & ∙ θ

19
Mécanique des Milieux Continus Polytech Grenoble, Geo3

On doit simplement remplacer  et  par leurs expressions polaires dans l'expression


de l'accélération. On calcule d'abord la fonction  :
7  7  7
 ,  1 = = =
^    ^ & ∙ θ1  & ∙ θ1 ^& 

Par ailleurs on a montré dans la question précédente que l'accélération est exprimée
par :
o , ' 1 = −  ,  1 ∙  ∙    ∙  1

o , ' 1 = −  ,  1 ∙ & ∙ θ ∙   1


  & ∙ θ ∙ 

Or on a la formule de changement de repère suivante : T = 


  

Finalement :
7
o , ' 1 = −  ,  1 ∙ & ∙ T = − ∙ 
^  &  T

L'accélération est donc purement centripète.

7. Le champ de vitesse est donné dans l'énoncé sous forme eulérienne par :
 = − ,  1 ∙    ,  1 ∙  =  ,  1 ∙ − ∙    ∙  1
n

En coordonnées polaires, on a donc :


 =  ,  1 ∙ −& ∙ θ ∙ 
n  1 =  ,  1 ∙ & ∙ 
  & ∙ θ ∙  +

7
 =
n ∙ 
^& +

La vitesse est purement orthoradiale (normale au rayon), et on en déduit que les


particules ont une trajectoire circulaire autour de l'origine. Sans calcul, on obtient la
représentation lagrangienne :
&=
7
h = Θ  '#
^& 
 = i

8. Pour passer des coordonnées polaires aux coordonnées cartésiennes, on utilise les
formules suivantes :

20
Mécanique des Milieux Continus Polytech Grenoble, Geo3

 
& =    2222222 2222222θ = 2222222 2222222θ =
r r

i i
 = i  i 2222222 2222222Θ = 2222222 2222222Θ =
 
y
On a par ailleurs  = Θ  ' , et donc :
WT K

7 
θ = s Θ  ' =
^&  r

et donc, en développant le cosinus :

 7 7
= sΘ ∙ s  ' − sΘ ∙ s  '
r ^&  ^& 

 7 7
= sΘ ∙ s Z   '[ − sΘ ∙ s Z '[
r ^   1 ^   1


= sΘ ∙ s ,  1 ∙ ' 1 − sΘ ∙ s ,  1 ∙ ' 1
r
 i i
= ∙ s ,  1 ∙ ' 1 − ∙ s ,  1 ∙ ' 1
r  

Finalement, on a :

 = i ∙ sg' 1 − i ∙ sg' 1

avec g =  ,  1. On peut développer le sinus de la même manière, et obtenir


l'expression lagrangienne de  recherchée :

7 
sθ = s Θ  
' = sΘ  g' 1 =
^& r
 i i
= sΘ ∙ sg' 1  sΘ ∙ sg' 1 = ∙ sg' 1  ∙ sg' 1
r  

Finalement, on a :

 = i ∙ sg' 1  i ∙ sg' 1

On obtient finalement le système recherché :

21
Mécanique des Milieux Continus Polytech Grenoble, Geo3

 = i ∙ sg' 1 − i ∙ sg' 1
h = i ∙ sg' 1  i ∙ sg' 1#
 = i

9. Les trajectoires des particules sont planes et circulaires, avec une vitesse angulaire
g2constante pour chaque particule et qui décroit en relation hyperbolique avec la
distance à l'axe origine. L'accélération est purement centripète. On appelle ce
mouvement un tourbillon.

10. On a montré à la question 3 que le tenseur des taux de déformation eulériens est
égal au tenseur gradient de vitesse, dont la matrice dans la base / = ,
 ,
  1 est :

 ,  1    −  0


} =  ∙ ~ −  −  q
0 = p
  
0 0 0H

On se place au point  = & ∙  , et on a donc  = 0 et  = &. Par conséquent, la


matrice de pq exprimée dans la base / en ce point est donnée par :

0, &1 0 & 0 7 0  0


q=
p ∙ ~&  0 0  = ~ 0 0
& ^& 
0 0 0H 0 0 0H

On reconnaît un glissement simple tel que défini dans le cours.

La base propre d'une telle déformation est nommée /′ = ($ , $ , $ , avec :

 ¤ ¤
$ =  
 

   
¤ ¤ #

$  =  −
 

    
 
$ = 

Dans cette base, la matrice de pq vaut :

7  0 0
q=
p ~0 − 0
^& 
0 0 0 H

22
Mécanique des Milieux Continus Polytech Grenoble, Geo3

3. Déformations

Exercice A. Soit un milieu soumis à un tenseur de déformation l, dont la matrice dans
une base /  0,
 ,
  1 donnée est :

  
l = ~  
  H

1. Calculer la trace de l , ainsi que son déterminant. Calculer la matrice du tenseur l  ,


et en déduire le deuxième invariant du tenseur donné par l'expression classique
 
l¥¥ = M('&2l  − '&2(l  N. En déduire sans calcul le polynôme caractéristique du

tenseur l .

2. Développer le déterminant du tenseur l −   et retrouver l'expression du polynôme


caractéristique.

3. Développer différemment $'(l −  pour l'exprimer sous la forme suivante,


beaucoup plus intéressante : $'(l −  = − − l 1 − l 1 − l 1. En déduire
les valeurs des déformations principales et la forme de la matrice de l dans sa base
principale.

4. Ordonner les déformations principales, puis calculer les coordonnées du vecteur


propre 
 correspondant à la déformation principale intermédiaire l .

23
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Corrigés

Exercice A.
1. La trace de l s'obtient immédiatement :
  

'&2l  '& ¢~   £   .  .     l¥
  H

Le déterminant d'une matrice 3*3 se calcule par la formule classique : somme des produits
des termes selon les diagonales descendantes moins somme des produits des termes selon les
diagonales montantes :
4  
$  §  4 ∙  ∙    ∙ § ∙ 3   ∙ $ ∙ ℎ − 3 ∙  ∙  − ℎ ∙ § ∙ 4 −  ∙ $ ∙ 
3 ℎ 

Appliqué à l , on a :
  
   =  ∗  ∗    ∗  ∗    ∗  ∗  −  ∗  ∗  −  ∗  ∗  −  ∗  ∗ 
  

$'2l = − = l¥¥¥

Le tenseur l  est le produit contracté de ª par lui-même. Sous forme matricielle dans la
base / = ,
 ,
  1, on a :

        
l  = ~    ∙ ~    = ~  
  H   H    H

La trace de ce tenseur vaut '&2l =      = 

Le second invariant principal de l s'obtient directement par la formule de l'énoncé :


  
l¥¥ = M('&2l  − '&2(l  N =  − 1 = −
 

Par application de la formule du cours, on en déduit le polynôme caractéristique du


tenseur des déformations linéarisées :

$'(l −   = −  l¥ ∙  − l¥¥ ∙   l¥¥¥ = −   ∙    ∙  − 

24
Mécanique des Milieux Continus Polytech Grenoble, Geo3

2. On cherche le déterminant du tenseur l −   qui s'exprime sous forme matricielle


dans la base / = ,
 ,
  1 par :

    0 0 −  
l −   = ~   −  ∙ ~0  0 = ~  −  
   H 0 0 H   − H

Donc :

−  
$'(l − l =   −  
  −

$'(l −  =  − 1 ∗  − 1   ∗  ∗    ∗  ∗  …


… −  ∗  −  1 ∗  −  ∗  ∗  −  1 −  −  1 ∗  ∗ 

$'(l −  =  −   1 − 1     −    −    −   

$'(l −   =  −    −    −  −   

$'(l −   = −   ∙    ∙  − 

On retrouve donc par le calcul détaillé le résultat de la question précédente.

3. On cherche à obtenir un développement du polynôme caractéristique faisant intervenir


directement ses racines, sous la forme $'(l −   = − − l 1 − l 1 − l 1. Cette
opération est peu aisée si on démarre de l'expression développée du polynôme, et il est
plus judicieux de chercher cette expression dès le calcul du déterminant. On prend
donc comme point de départ :

−  
 
$'(l −  =   −  
  −

Selon les formules classiques de calcul de déterminant, il est possible de remplacer une
colonne (respectivement une ligne) par une combinaison linéaire des différentes
colonnes (respectivement des différentes lignes). On remplace donc la colonne C1 par
la combinaison C1-C2 :

−    −  1 −    −  
  −   =  −  −  1 −   =  −  −  
  − −  − 0  −

25
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Il est également autorisé de "sortir" une constante multiplicative d'une ligne ou d'une colonne
du déterminant, donc on peut simplifier la première colonne :
−  
 
$'(l −   =  − 1   −  
0  −
On remplace ensuite la ligne L1 par la combinaison linéaire L1+L2 :
0 − 
$'(l −  =  − 1  −  
0  −

On a le droit de permuter deux lignes (ou bien deux colonnes) en changeant le signe du
déterminant, on peut donc écrire :
 − 
$'(l −   = − − 1 0 −  
0  −

Enfin, on sait qu'un déterminant 3*3 de la forme obtenu peut être simplifié en un
déterminant 2*2 par la formule suivante :
 4 
 $
0  $ = ° °
 §
0  §

On peut donc écrire :


− 
$'(l −  = − − 1 ± ±
 −

Il reste quelques calculs, mais on sait qu'on a atteint le but principal, car un
déterminant 2*2 de cette forme est un polynôme de degré 2 dont on sait très bien
retrouver les racines. On a :
$'(l −   = − − 1\ − 1 − 1 −  ∗ ]

$'(l −  = − − 1 −  − 1

Le discriminant du polynôme de second degré est :


∆= −1 −  ∗  ∗ −1 = 

26
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Et les deux racines de ce polynôme sont :

  ¤∆ ¤
=
 

Finalement, le polynôme caractéristique du tenseur des déformations linéarisées est


donné par :

$'(l −  = − − 1 − 1  1

Les déformations principales sont les racines de ce polynôme, et on peut écrire la


matrice de l dans sa base principale :
 0 0
ª = ~0  0
0 0 − H´QU28Tµb8´¶U

4. On ordonne les déformations principales dans l'ordre décroissant, et on a :


l 0 0
l = ~ 0 l 0
0 0 l H´QU28Tµb8´¶U
avec :
l = 
 l =  #
l = −

On pourrait calculer les coordonnées dans / = ,  ,


  1 des trois vecteurs de la base
principale de déformation ( ,  ,  . On va le faire uniquement pour   car la
méthode est toujours la même.

D'après le cours, la direction principale de déformation 


 est définie par :

l 
 = l 


En développant de système, on obtient :

   , , ,


~   ·, ¸ = l ·, ¸ = ·, ¸
   H , H , H , H

27
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Dans cette expression, on a fait apparaître les trois coordonnées (, , , , ,  du vecteur

 dans la base / = ,
 ,
  1. En mettant ce système matriciel sous forme d'un
système d'équations linéaires, on a :

 ∗ ,   ∗ ,   ∗ , = ,


h ∗ ,   ∗ ,   ∗ , = , #
 ∗ ,   ∗ ,  , = ,
Soit :
 ∗ ,   ∗ ,   ∗ , = 0
h ∗ ,   ∗ ,   ∗ , = 0#
 ∗ ,   ∗ , = 0

On remarque immédiatement que , = −, et que , = 0. Or on sait que le vecteur

 est unitaire (de norme égale à 1), donc on peut écrire :

 , = ¤
 
¤ #
, = −
 
 , = 0

On note une indétermination dans le signe de , et , . Tous les signes sont possibles, à
condition que l'on vérifie , = −, .

Ce système était trivial, mais ce n'est pas toujours vrai. Dans le cas contraire, on peut
utiliser le pivot de Gauss, ou passer par une inversion de matrice.

28
Mécanique des Milieux Continus Polytech Grenoble, Geo3

4. Contraintes

Exercice A. Soit un milieu continu dont l'état de contrainte est homogène et


stationnaire (le tenseur de Cauchy est constant dans l'espace et le temps). On considère
au sein de ce milieu un domaine matériel p, en forme de tétraèdre tel que représenté
dans la figure ci-dessous dans un repère cartésien orthonormé 0¹, ,
  ,
  1. Les arêtes
parallèles aux axes de la base ont une longueur donnée .
Un expérimentateur réalise quelques mesures, et constate que :
-La composante normale de la force exercée sur la face arrière (correspondant au plan
 = 0) est égale à j
-La composante normale de la force exercée sur la face de gauche (correspondant au
plan  = 0) est égale à j
-La composante normale de la force exercée sur la face inférieure (correspondant au
plan  = 0) est égale à j
-La force exercée sur la face inclinée du tétraèdre est égale à º = º 
  º    º 


Déterminer les composantes de la matrice du tenseur de Cauchy dans la base


,
 ,
  1.

Exercice B. On considère un petit cube de centre  dans un milieu continu soumis à


des contraintes. On effectue trois expériences de chargement (a), (b), et (c),
respectivement caractérisées par les vecteurs contraintes suivants (»Y étant une
constante) :

(a).  
,  1 = 0 ; ,  1 = »Y 
  ; ,  1 = »Y 
 

(b). ,  1 = »Y 


  ;  
,  1 = 0 ; ,  1 = »Y 
 

(c). ,  1 = »Y 


  ; ,  1 = »Y 
  ; 
,  1 = 0

29
Mécanique des Milieux Continus Polytech Grenoble, Geo3

1. Exprimer les trois tenseurs de contraintes » ´1 , »


 ¼1 , et 
 »b1 correspondant à ces trois
cas de chargement.

2. On effectue les trois expériences simultanément en superposant les trois systèmes de


forces. Exprimer le tenseur de contraintes correspondant à cette nouvelle expérience.
Pour la suite de l'exercice on considérera uniquement ce nouvel état de contrainte.

3. Calculer le vecteur contrainte exercé sur une petite surface normale à la direction 
définie par l'expression  =      .

 normale à
4. Calculer le vecteur contrainte exercé sur une petite surface de direction ½
.

5. Donner sans calcul supplémentaire les contraintes principales et les directions


principales de contraintes.

6. Tracer le tricercle de Mohr correspondant à cet état de contraintes et donner la


valeur de la contrainte de cisaillement maximale. Donner également la valeur de la
contrainte moyenne et la matrice du déviateur de contraintes.

Exercice C. En un point ¾ d'un milieu continu et à l'instant ', la matrice du tenseur de


Cauchy dans une base cartésienne orthonormée fixe / = ,  
,
  1 est donnée par
l'expression suivante, où 7 est une constante adimensionnelle :

0.7 .7 0
~
» = .7 .7 0
0 0 .

1. Montrer que le calcul des trois contraintes  À 1,  = ,, et l'utilisation des
propriétés des cercles de Mohr permettent de calculer les contraintes principales qu'on
notera » , » , et » . Illustrer le cas 7 = .

2. Déterminer les valeurs de 7 correspondant à un état triaxial de révolution.

3. Pour la suite de l'exercice, on pose 7 = . Déterminer les directions principales de »


en ¾.

¤ 
 =
4. Calculer la contrainte pour une facette de normale    
  . Tracer le point
 
correspondant sur le cercle de Mohr.

5. Déterminer pour 7 = , puis pour 7 = , la valeur de la contrainte de cisaillement


maximum ainsi que la direction de la normale correspondante.

30
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Corrigés

Exercice A.
On recherche les termes de la matrice de » à partir d'un certain nombre de mesures de
forces effectuées sur les faces d'un domaine matériel. L'état de contrainte est supposé
homogène et stationnaire, et le tenseur de Cauchy est toujours symétrique, on cherche
donc en tout et pour tout six inconnues scalaires : » , » , » , » , » , et » .
On va d'abord travailler sur la force appliquée sur la face arrière du tétraèdre (plan
  0). L'énoncé dit que la composante normale à cette force est égale à j . On peut
également remarquer que la facette en question (face triangulaire arrière) est de
normale sortante −.
On va chercher à exprimer le vecteur contrainte appliqué sur cette facette. On peut
l'énoncer ainsi à partir de la formule du cours :
− 1 = » ∙ − 1

 » » »


Avec :  = ~0 et » = ~» » » 
0 H » » » H

Donc on a :
» » » − −»
− 1 = ~»
 » »  ∙ ~ 0  = ~−» 
» » » H 0 H −» H

Par ailleurs, on sait que la composante normale de la force appliquée à la facette vaut
j . Si on note  la surface de la facette triangulaire, on peut intégrer la contrainte
− 1 sur toute la surface de la facette (pour obtenir une force), et considérer sa

composante normale (par projection sur la normale, c'est-à-dire par produit scalaire
avec le vecteur normal sortant de la facette). On a donc :
−» −
− 1 ∙ − 1 =  ~−»  ∙ ~ 0  = »
j =  × 
−» H 0 H
K
La surface  de la facette est donnée par  = b , et on peut finalement écrire :

j
» =


31
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Par des raisonnements identiques sur les facettes de normales − et −, on obtient
de la même manière :
j
» =

j
» =


On a déterminé les trois termes diagonaux de la matrice de ». Il nous reste à


déterminer les termes non-diagonaux. On va donc s'intéresser à la facette restante.

Sa normale sortante est un vecteur unitaire, colinéaire à la direction ~. Il s'agit donc

4
d'un vecteur d'allure générale d4f, où 4 est un réel à déterminer. Par ailleurs, ce vecteur
4
doit être unitaire, donc sa norme doit être égale à 1. On va donc avoir :

4  4  4 = 

et donc :

4=
¤

 Â ¤
La normale sortante à la facette est donc le vecteur  = ·¤¸ .
 Â ¤ H

On a également besoin de la surface ′ de cette facette. Il s'agit d'un triangle équilatéral


de côté ¤, et on obtient après quelques calculs trigonométriques et/ou Pythagore :

  ¤
 =

On peut maintenant écrire le vecteur contrainte subi par cette facette en n'importe
lequel de ses points (puisque l'état de contrainte est homogène) :
1 = »

et donc :

32
Mécanique des Milieux Continus Polytech Grenoble, Geo3

» » » ¤


1 = ~» » »  ∙ ·¤¸
» » » H
¤ H

»  »  »



1 =
 ~»  »  » 
¤ »  »  »
H

L'énoncé précise que la force subie par la facette est mesurée comme étant égale à :

º = º 
  º 
  º 


Or cette force est également le produit du vecteur contrainte (homogène) par la surface
de la facette, on a donc :
»  »  »
′
º = º 
  º 
  º  1 =
 = ′ ∙  ~»  »  » 
¤ »  »  »
H

Cette équation matricielle nous fournit trois équations scalaires :


′
º = »  »  » 1
 ¤ 

′
º = »  »  » 1#
 ¤

º = ′ »  »  » 1
  ¤   

On peut alors remplacer les termes diagonaux de la matrice par leurs valeurs trouvées
en début d'exercice, et faire de même avec la surface de la facette :
  j
º = M   »  » N
  

 j
º = M»    » N#
  

  
j
º = M»  »   N
 
 

On dispose alors d'un système de trois équations à trois inconnues : » , » , et » .

33
Mécanique des Milieux Continus Polytech Grenoble, Geo3

º − j
»  » =
 
º − j #
»  » =
 
 º − j
»
   » =


On obtient finalement les résultats demandés suivants :

º − j 1  º − j 1 − º − j 1
 » = 
 
º − j 1 − º − j 1  º − j 1 #
» =
 
 −º − j 1  º − j 1  º − j 1
»
  =


Exercice B.
1. Considérons l'état de chargement (a). On nous donne trois équations vectorielles :

 
,  1 = 0 ; ,  1 = »Y 
  ; ,  1 = »Y 
 

D'après la première équation, on peut écrire :

» » »  »


,  1 = »´1 ∙  = ~»
 » 
»  ∙ ~0 = ~»  = 0
» » » 0 »

On en déduit directement que » = » = » = 0

On considère maintenant la deuxième équation, qui nous permet d'écrire :

» » » 0 »


,  1 = »´1 ∙  = ~»
 » »  ∙ ~ = ~»  = »Y 

» » » 0 »

On en déduit que » = » = 0 mais aussi que » = »Y .

Enfin, on considère la troisième équation, qui donne directement :

» » » 0 »


,  1 = »´1 ∙  = ~»
 » »  ∙ ~0 = ~»  = »Y 

» » »  »

On en déduit que » = » = 0 et on confirme que » = »Y .

34
Mécanique des Milieux Continus Polytech Grenoble, Geo3

On a donc calculé les six termes indépendants de la matrice du tenseur de Cauchy dans
la base / = ,
 ,
  1, et on peut écrire cette matrice :

0 0 0
»´1 = ~0 0 »Y 
0 »Y 0 H

On rappelle que cette matrice est uniquement valable pour l'état de chargement noté
(a) dans l'énoncé. En appliquant des raisonnements identiques, on obtient également
les matrices des tenseurs de contraintes correspondant aux états de chargement (b) et
(c) :
0 0 »Y
»¼1 = ~0 0 0
»Y 0 0 H

0 »Y 0
»b1 = ~»Y 0 0
0 0 0H

2. Sans calcul supplémentaire, on peut tirer parti du fait que la superposition de


plusieurs états de contraintes correspond à la somme des tenseurs de contraintes
correspondants. Si on superpose les chargements (a), (b), et (c), on obtient donc le
tenseur de contraintes dont la matrice dans la base / = ,
 ,
  1 s'exprime par :

0 »Y »Y
»´1Ã¼1Ãb1 = »´1  »¼1  »b1 = ~»Y 0 »Y 
»Y »Y 0 H

3. On cherche à exprimer le vecteur contrainte sur une facette dont on nous donne la
direction . Il nous faut d'abord déterminer les coordonnées du vecteur unitaire normal
à cette facette. Ce vecteur est colinéaire à , mais de norme égale à 1. On peut le noter

′ et l'exprimer ainsi :

    


′ = = = ~ 
{{ ¤ ¤
 H

Le vecteur contrainte appliqué à la facette définie par cette normale unitaire est ensuite
donné directement par la formule du cours :
0 »Y »Y
 
, 1 = 
  
 , ′ = »´1Ã¼1Ãb1 ∙ ′ = ~»Y 0 »Y  ∙ ~
»Y »Y 0 H ¤   H
Et donc :

35
Mécanique des Milieux Continus Polytech Grenoble, Geo3

»Y  »Y
, 1 =
 ~ =  = »Y 
′
¤  ¤
H

4. On cherche à déterminer le vecteur contrainte appliqué à une facette de normale ½,


dont on sait seulement que cette normale est également normale au vecteur  défini
plus haut. On cherche donc , ½1, sachant que ½ ∙  = 0. Ecrivons ces deux données
sous forme mathématique :

½
 ∙  = ½   ½   ½  = ½  ½  ½ = 0

0 »Y »Y ½ ½  ½
, ½
 1 = »´1Ã¼1Ãb1 ∙ ½
 = ~»Y 0 »Y  ∙ ~½  = »Y ~½  ½ 
»Y »Y 0 H ½ H ½  ½ H

On en déduit :
½
, ½
 1 = −»Y ~½  = −»Y ½

½ H

5. D'après les formules précédentes, on peut formuler plusieurs observations. On a


montré que :
  = »Y 
 , ′ ′

Par conséquent la direction 


′ est une des directions principales de contrainte. Le
vecteur unitaire 
′ est donc l'un des vecteurs de la base principale de contraintes, et la
contrainte principale correspondante est égale à »Y .

Par ailleurs, on a montré que, pour tout vecteur ½ perpendiculaire à 
′, on a :

, ½
 1 = −»Y ½


On en déduit que −»Y est une racine double et que le tenseur de Cauchy est
cylindrique. Dans sa base principale, sa matrice s'exprime par :
»Y 0 0
»´1Ã¼1Ãb1 =~ 0 −»Y 0 
0 0 −»Y H´QU28Tµb8´¶U

La base principale est formée du vecteur  ′ et de tout couple de vecteurs



perpendiculaires appartenant au plan normal à ′.
36
Mécanique des Milieux Continus Polytech Grenoble, Geo3

6. En ordonnant les valeurs de contraintes (et en supposant que »Y > 0), on a :


» 0 0
»´1Ã¼1Ãb1 = ~0 » 0
0 0 » H´QU28Tµb8´¶U

avec : » = »Y et » = » = −»Y . Le tricercle de Mohr a l'allure suivante :

Les valeurs propres » et » sont confondues, donc le tricercle se résume à un cercle


unique. Par ailleurs, la contrainte de cisaillement maximale et la contrainte moyenne
sont données par les formules de cours :
» − » 
ÄÅ´- = = »Y
 
 »  »  »
»Å = '&» = =0
 

La contrainte moyenne est nulle, donc le tenseur de contraintes n'a pas de partie
sphérique, et le déviateur de » est égal à ».

Exercice C.
1. On dispose de la matrice du tenseur des contraintes dans la base / = ,
 ,
  1 :

0.7 .7 0
» = ~.7 .7 0
0 0 .

37
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Suivons l'énoncé, et calculons les trois vecteurs contraintes existant au point ¾, pour
des facettes de normales respectives , , et 
 :

0.7 .7 0  0.7



 1 = » = ~.7 .7 0  ~0 = ~.7 = 0.7  .7
0 0 . 0 0
0.7 .7 0 0 .7

 1 = » = ~.7 .7 0  ~ = ~.7 = .7  .7
0 0 . 0 0
0.7 .7 0 0 0

 1 = » = ~.7 .7 0  ~0 = ~ 0  = .
0 0 .  .

Chacun de ces vecteurs peut être décomposé en une composante normale et une
composante tangentielle, afin d'être représentés dans le plan de Mohr. Par exemple,
pour le vecteur contrainte appliqué à la facette de normale , on peut écrire :

  1 = µ ∙    ∙ '

Dans cette équation très générale, ' est un vecteur unitaire appartenant au plan de la
facette, et µ et  sont respectivement la composante normale et la composante
tangentielle du vecteur contrainte. Par identification avec l'expression trouvée plus
haut, on peut écrire :
' = 

µ = 0.7 et  = .7

En appliquant la même opération à 


 1 et   1, on peut écrire de la même manière


µ = .7 et  = .7

µ = . et  = 0

Avec ces résultats, on peut tout à fait placer sur le plan de Mohr les points
correspondant à chacun de ces trois vecteurs contraintes. Il s'agira des points de
coordonnées µ,  1, µ ,  1, et µ, 1. Il faut noter pour l'instant que la
connaissance de ces points n'est pas suffisante pour tracer le tricercle de Mohr et donc
pour déterminer les contraintes principales.
On doit maintenant appliquer un certain nombre de raisonnements liés à des propriétés
du cercle de Mohr. On remarque que la composante tangentielle de   1 est nulle, ce
qui signifie par définition que  est une direction principal de contrainte. On en déduit
donc que les vecteurs  et 
 sont situés dans le plan des deux autres directions
principales (pour l'instant inconnues), et que les points du plan de Mohr correspondant
38
Mécanique des Milieux Contin
nus Polyteech Grenoble, Geo3

aux deux vecteurs contrainntes liés à ces facettes sont situées sur lee même cercle de
Mohr (c'est la définition d'uun cercle donné dans le cours). Par consséquent, les points
µ ,  1 et µ ,  1 appartiennent tous les deux au même cerclle de Mohr. C'est
aussi le cas de leurs symétrriques par rapport à l'axe horizontal, c'estt-à-dire des points
µ , − 1 et µ , − 1.

On tire enfin parti des règlees d'orientation données dans le cours, et ono observe que les
deux facettes de normales 
 et  sont orthogonales, et donc que les
l points qui leur
correspondent dans le plan de Mohr sont nécessairement diamétraleement opposés sur
le cercle de Mohr auquel ilss appartiennent. On en déduit que les coorrdonnées du centre
du cercle de Mohr que l'on recherche sont au milieu d'un segment reliant deux points
diamétralement opposés :
µ  µ
µ 0.7  .7
M , 0N = M , 0N = .¯7, 01
 

Le rayon de ce cercle se calccule ensuite par Pythagore :

7 1  .71 = 7. − .¯1 . 0


 = .7 − .¯7 1  ¯7

On en déduit qu'une des contraintes principales est égale à ¯7 − ¯7  −7, et
qu'une autre contrainte prinncipale vaut ¯7 . ¯7  ¯¯7. Quaant à la troisième,
on l'a déjà déterminée im mplicitement au début de l'exercice, et elle
e vaut . La
situation est représentée danns le plan de Mohr pour 7   :

39
Mécanique des Milieux Continus Polytech Grenoble, Geo3

2. Un état triaxial de révolution correspond à n'importe quel état de contraintes pour


lequel deux des contraintes principales sont égales.

Dans notre cas, sans accorder d'attention à l'ordre des contraintes principales (puisque
celui-ci dépend de 7 ), on peut énoncer que les contraintes principales sont égales à :

»  −72;2»  ¯¯72;2»  

On peut relever dans le cas présent trois situations pour lesquelles l'état de contraintes
est triaxial de révolution :

»  » 2 ⇒ 2 −7  ¯¯72 ⇒ 27  0

» = » 2 ⇒ 2 −7 = .2 ⇒ 27 = −.

» = » 2 ⇒ 2¯¯7  2 ⇒ 27  8

3. On pose 7  , on ordonne les contraintes, et on a donc »  , »  ¯¯, et


»  −.

On cherche les directions principales de contraintes 0,


 ,
  1. Plus précisément, on
cherche à trouver les coordonnées de ces trois vecteur propres dans la base d'origine
,
 ,
  1.

On a déjà démontré que  est une direction propre, et lorsque 7 =  cette direction
propre est celle de la contrainte principale maximale. On peut donc énoncer sans calcul
que :
0
 =  = ~0

 H

On cherche alors les coordonnées des deux autres vecteurs principaux, tels que :

 = , ∙   , ∙   , ∙ 

 = , ∙   , ∙   , ∙ 

Puisque 
 =  , les vecteurs 
 et  sont orthogonaux à 
 , et on a , = , = 0.

Il nous reste quatre scalaires à trouver : , , , , , , et , .

On applique la formule classique définissant une direction principale de contrainte :

»À = » 
À

40
Mécanique des Milieux Continus Polytech Grenoble, Geo3

0. . 0 , 0.,  .,


or on a » = ~. . 
0  ~ ,  = ~.,  ., 
0 0 . H 0 H 0 H

Donc on obtient » = (0.,  .,   (.,  ., 

Par ailleurs on a » = » 


 = ¯¯ ∙ Ç, 
  , È


On en déduit le système suivant :


¯¯, = 0.,  ., #
v
¯¯, = .,  .,

−.,  ., = 0#
v
., − ., = 0

−. . , 0


    =  
. −. , 0

La solution de ce système est triviale : , = 0 et , = 0. Cette solution n'est


évidemment pas celle qui nous intéresse. En se penchant un peu plus sur le système,
on observe que ses deux équations sont en fait identiques. C'est le genre de piège que
les calculs peuvent faire apparaître.
 doit être
Il y a une information que l'on n'a pas encore utilisée : le fait que le vecteur 
unitaire. On a donc un nouveau système :
−.,  ., = 0
v #
 
,  , =

De la première équation il vient , = É , , que l'on peut intégrer dans la seconde
équation :

 
,   =
 ,

Finalement, on obtient : , = 0. et , = 0.. On note enfin :


0.
 = ~0.
0 H

On pourrait appliquer la même méthode pour obtenir , , et , , mais on peut aussi
tirer parti du fait que ,
 ,
  1 est une base orthonormée directe. On peut donc
utiliser le produit vectoriel et écrire :
41
Mécanique des Milieux Continus Polytech Grenoble, Geo3

 =  ∧ 



Donc on a :
0 0. −0.
 = ~0 ∧ ~0. = ~ 0. 
 H 0 H 0 H

4. On définit la facette de normale  = ¤  


 . Le vecteur contrainte subi par
  
 
cette facette s'obtient directement par la formule du cours :
0. . 0 ¤Â .

 1 = »
 = ~. . 0  ~ Â  = ~.¯
0 0 . H 0 H 0 H

Pour tracer le point correspondant à ce vecteur contrainte dans le plan de Mohr, on a


besoin d'une composante normale et d'une composante tangentielle. La composante
 1 sur la direction
normale s'obtient, comme d'habitude, par projection du vecteur 
normale à la facette, c'est-à-dire par un produit scalaire :
. ¤Â
 
µ =   1  = ~.¯ ∙ ~ Â  = .
 ∙ 
0 H 0 H

 1 est la somme d'un vecteur µ 


Le vecteur contrainte   orienté dans la direction
normale et d'un autre vecteur orienté dans une direction tangentielle 
' qui n'est pas
connue. La composante tangentielle est donc la norme du vecteur obtenu en
retranchant la composante normale du vecteur 
 1 à lui-même :


 = Ë  1 − µ Ë


. ¤Â
 = Ì~.¯ − . ~ Â  Ì = .
0 H 0 H

On peut alors reporter le point correspondant sur le graphique (voir figure précédente).
On constate qu'il appartient au cercle 1. C'est logique, car ce point est défini pour une
¤ 
facette de direction normale   = 

 qui appartient au plan des deux
  
 
  1 (c'est le plan normal à  ).
directions principales ,

42
Mécanique des Milieux Continus Polytech Grenoble, Geo3

5. Pour 7 = , on a » = . et » = − . La contrainte de cisaillement maximale est


donc donnée par :
» − »
ÄÅ´- = = .

L'abscisse du point correspondant sur le plan de Mohr correspond à celle du centre du


grand cercle de Mohr, soit:
»  »
= .

Le point du cercle de Mohr correspondant à la facette de cisaillement maximal est


donc de coordonnées .,.1 (voir figure).
Par définition, la facette qui reçoit cette contrainte possède une normale qui appartient
au plan des deux directions principales  et 
 (qu'on appelle "plan de cisaillement
maximal"). En fait, la facette de cisaillement maximal correspond toujours à une
orientation situé selon la bissectrice des directions principales 
 et  , qui est donc
donnée par :
0 − 0.
Ï Ò

   Î0  0.Ñ −0.
= Î Ñ = ~ 0. 
 Î  Ñ 0.¯
Î 0 Ñ
Í  Ð

Il faut noter que ce vecteur n'est pas unitaire (même s'il suffit à définir la direction de
la normale à la facette).
Pour 7  , la méthode est strictement la même, excepté que le vecteur  a changé, et
il faut donc le recalculer pour obtenir la bonne orientation de la facette de cisaillement
maximal.

43
Mécanique des Milieux Continus Polytech Grenoble, Geo3

5. Elasticité

Exercice A. Soit le massif rectangulaire de grande longueur représenté sur la figure


suivante :

On pose plusieurs hypothèses :


-Comportement élastique linéaire, de paramètres l et Ó.
-Face supérieure soumise à un pression uniforme :.
-Face inférieure en appui glissant et non-frottant, fixée seulement sur l'axe 0¹, 1
 .
-Forces de volume négligeables.
-Hypothèse des petites perturbations.
-On se place dans la zone centrale du massif pour se prémunir des effets de bord.
-Champ de déplacement donné par Ô =  ∙ i  k ∙ i 1  −k ∙ i  / ∙ i 1 .

1. Représenter la section du massif contenant le point ¾, ainsi que les conditions aux
limites sur cette section.
2. Calculer le tenseur des déformations ª, et donner ses valeurs principales et ses
directions principales. De quel type de déformation s'agit-il ?
3. Exprimer analytiquement le tenseur de Cauchy » (utiliser les coefficients de Lamé
pour simplifier les notations).
4. En utilisant les conditions aux limites en contraintes, calculer les expressions des
paramètres2 et /. En utilisant les conditions aux limites en déplacement, calculer k .
5. On pose l = 00002¾Õ4 et Ó = 0. (Massif en acier). Les dimensions de la
section sont } = 0.Ö et  = 0.Ö, et le chargement vaut  = 002¾Õ4. Calculer
le vecteur déplacement du point ¾. Tracer la déformée de la section.
44
Mécanique des Milieux Contin
nus Polyteech Grenoble, Geo3

Exercice B. On s'intéresse à un barrage poids de section triangulaire, dont la géométrie


est proposée sur la figure suivante. Pour rendre le problème plan, on s'intéresse en
particulier à une "tranchee" de ce barrage située suffisamment loin des massifs
d'ancrage latéraux. Le barraage est soumis à son poids m¼ 3 et à la poussée
p mU 3 de
l'eau, la surface de celle-ci étant située à une hauteur ℎ égale à cellee du barrage. A sa
base, le barrage est parfaitem
ment encastré.
Le matériau composant lee barrage est supposé élastique linéairee et isotrope, de
paramètres l et Ó. On négligge la pression atmosphérique.

1. Expliquer pourquoi il s''agit d'un problème de déformations plaanes, et donner la


forme générale des inconnuees Ô, ª, et ».
2. Décrire précisément les conditions
c aux limites d'une "tranche" du barrage.
b Ecrire les
équations d'équilibre mécannique, d'abord à l'intérieur du domaine puiis au niveau de ses
conditions limites en contraiintes.
3. On postule que tous les termes
t du tenseur de contrainte sont des fonctions
f linéaires
des coordonnées  et  . Donner
D l'expression générale des termes de
d la matrice d'un
tel champ de contrainte. Montrer
M qu'il vérifie l'équation de Beltram
mi, et est donc la
solution du problème élastiqque.

45
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Rappel : l'équation de Beltrami est donnée par :



∆»    2'&2» = 0
3&4$ 23&4$
Ó

Le laplacien d'un tenseur (par exemple ∆») est le tenseur dont chaque terme est le
laplacien des termes du tenseur d'origine (ici »). On a donc ∆»1 = ∆(» . Les
autres termes de l'équation sont calculables.

4. A l'aide de la loi de Hooke, montrer que » = Ó»  »1.


5. A l'aide des équations d'équilibre et des conditions aux limites en contraintes,
calculer » en tout point du système. A l'aide de la loi de Hooke, en déduire ª en tout
point du système.
6. Montrer que la contrainte tangentielle la plus intense (en valeur absolue) du système
est obtenue au point /, et calculer cette contrainte.

46
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Corrigés

Exercice A.
1. La section du massif contenant le point ¾ est un rectangle. Il est soumis aux
conditions aux limites suivantes :
-face supérieure : pression uniforme égale à :.
-faces latérales : pression uniforme égale à 0.
-face inférieure : déplacement vertical bloqué, déplacement horizontal bloqué
au point central.

2. Le champ de déplacement est donné par l'expression :


 ∙ i  k ∙ i
  −k ∙ i  / ∙ i #
Ô
0

Il s'agit manifestement d'un mouvement plan. Le tenseur des déformations linéarisées


se calcule par :
 @
ª =  2Ô
3&4$  2Ô
  2 (3&4$  

IØJ IØJ IØJ
ÏI-J I-K I-L Ò
Î IØK Ñ
 k 0
Or 
3&4$ 2Ô  = ÎIØK IØK
= ~−k / 0
I-J I-K I-L Ñ
ÎIØL IØL IØL Ñ 0 0 0
ÍI-J I-K I-L Ð

47
Mécanique des Milieux Continus Polytech Grenoble, Geo3

 0 0
Finalement, ª = ~ 0 / 0
0 0 0

On constate qu'il s'agit d'une déformation plane, que la base principale de déformation
est la base d'origine ,
 ,
  1, et que les déformations principales sont Ù, /, 0Ú.

3. On applique la loi de Hooke formulée en fonction des coefficients de Lamé :

» = Û ∙ ª   ∙ '&ª1 ∙  
Ü ÜÝ
avec Û = ÃÝ1 et  = ÃÝ1zÝ1.

On a '&ª1 =   /, et il vient donc :


 0 0  0 0
» = Û ∙ ~ 0 / 0   ∙   /1 ∙ ~0  0
0 0 0 0 0 
Et finalement :
Û   ∙   /1 0 0
» = · 0 Û/   ∙   /1 0 ¸
0 0  ∙   / 1

Les tenseurs » et ª sont homogènes, c'est-à-dire qu'ils sont égaux en tout point du
système.
4. Sur la face supérieure, la pression appliquée peut se noter :

Þ = − ∙ 


Par ailleurs, la condition aux limites d'équilibre sur cette face supérieure s'écrit :

» = 
Þ

La normale unitaire sortante à cette surface est −, donc on a »− 1 = 
Þ

Par conséquent, on peut écrire :−(Û   ∙   /1 ∙  = − ∙ 




Et finalement :2Û   ∙   /1 = 
On applique la même méthode sur une des faces latérales, et on obtient :
Û/   ∙   /1 = 0

Il en résulte un système linéaire de deux équations à deux inconnues  et / :


48
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Û  1   = #
v
  Û  1/ = 0
záâ
La deuxième équation donne / = ãÃá , que l'on réinjecte dans la première pour
á K
obtenir  Û   − ãÃá  = .

Après calcul, il vient :


−Û −  
= 22'22/ = 
ÛÛ  1 ÛÛ  1

Sur la face inférieure, le déplacement est bloqué verticalement. En posant i = 0, on


peut donc écrire :
Ô 0, i , i 1 = 0

Or Ô i , i , i 1 =  ∙ i  k ∙ i .
On en déduit que k ∙ i = 022∀i , i 1, et donc que k = 0.
5. En utilisant les résultats précédents, on peut réécrire le champ de déplacement :
−Û − 
 ∙ i =  ∙i
 ÛÛ  1 
 =
Ô  #
/ ∙ i =  ∙i
 ÛÛ  1 
 0

Les coordonnées du point ¾ sont i =  = 0.Ö et i = } = 0.Ö.


Par ailleurs, les coefficients de Lamé valent :
l 0000
Û= = = 002¾Õ4
  Ó1   0.1

0000 ∗ 0.
= = 02¾Õ4
  0.1 −  ∗ 0.1

Il vient finalement  = −. ∙ 0z et / = . ∙ 0zÉ , et enfin :


−0.ÖÖ
Ô¾1 =  0.ÖÖ #
0

49
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Exercice B.
1. Il s'agit d'un problème de déformations planes car la section du barrage est identique
sur toute sa longueur, de même que le chargement qu'il subit. On en déduit que son
champ de déplacement est entièrement situé dans le plan perpendiculaire à sa longueur
(plan d'une "tranche"), et que sa déformation dans la direction de sa longueur est
empêchée par les tranches voisines. Les champs inconnus ont la forme suivante :
Ô i , i 1 ª ª 0 » » 0
Ô = Ô i , i 1# 2222222222222222ª = ~ª ª 02222222222222222» = ~» » 0 
  
0 0 0 0 0 0 »

Par ailleurs, en élasticité linéaire on applique systématiquement l'HPP, et on peut


utiliser indifféremment les vecteurs position initial i et final .
2. Sur une "tranche" du barrage, les conditions aux limites sont les suivantes :
-Sur sa base il est parfaitement encastré : Ô = 022½&2åH
-Sur la face immergée (notée æ ), il subit la pression de l'eau
-Sur la face à l'air libre (notée æH ), il subit la pression atmosphérique qui est
jugée négligeable.
Les équations d'équilibre du milieu continu sont :
-A l'intérieur du domaine, forme locale du PFD appliquée à l'équilibre statique
du solide chargé :
 2½&2p
$ç%2»  m3 = 0

-Sur la surface immergée, la normale sortante est  = −, et la contrainte


subie par la surface est proportionnelle à la profondeur  et orientée dans la
direction de 
 :

»− 1 = mU 3 ∙ 22½&2


 æ

-Sur la surface à l'air libre, la normale sortante est  = ¤ 
 −  1, et la
contrainte subie par la surface est nulle :


» ¢ 22½&2æH
 −  1£ = 0
¤

3. L'énoncé postule que tous les termes du tenseur de contrainte sont des fonctions
linéaires des coordonnées  et  , on a donc proportionnalité entre chaque terme de
ce tenseur et ces coordonnées. On peut donc écrire les relations suivantes :
» = 4    
» = 4    
» = 4    
» = 4    

50
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Dans ces expressions, les termes Ù4 ,  , 4 ,  , 4 ,  , 4 ,  Ú sont des
constantes de proportionnalité qui sont pour l'instant inconnues.
Quant aux deux autres composantes du tenseur de Cauchy, on a déjà démontré qu'elles
sont nulles car on est en déformations planes : » = » = 0.
Pour que ce champ de contrainte soit admissible, il doit vérifier l'équation de Beltrami:

∆»    2'&2» = 0
3&4$ 23&4$
Ó

On écrit d'abord les différents termes de cette équation sous forme indicielle. Le
premier terme permet d'écrire :
D  » D  » D  »
∆»1 = ∆(»  =  
D D D

le deuxième terme est légèrement plus long à dériver. On développe d'abord la trace :
 23&4$
(3&4$  2'&2» = 3&4$
 23&4$
 »  »  » 1
 

On développe ensuite le gradient du champ scalaire de cette trace, qui est donc un
champ vectoriel :
D
Ï »  »  » 1Ò
C D
Î  ÑG
Î D Ñ
 2'&2» = B3&4$
 23&4$
(3&4$  »  »  » 1ÑF
 B ÎD F
B Î  ÑF
Î D »  »  » 1Ñ
A ÍD   
ÐE

On calcule enfin le gradient de ce champ vectoriel, qui est donc un champ tensoriel :
D D D
Ï »  »  » 1 »  »  » 1 »  »  » 1Ò
Î D
 D D  D D  Ñ
Î D D D Ñ
 23&4$
 2'&2» =
(3&4$  ÎD D »  »  » 1 »  »  » 1
D  D D
»  »  » 1Ñ
Î   Ñ
Î D D D Ñ
»  »  » 1 »  »  » 1 »  »  » 1
ÍD D D D  D   
Ð

51
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Il vient :

D D  » D  » D  »


 23&4$
(3&4$  2'&2» = »  »  » 1 =  
 D D D D D D D D

Finalement, l'équation de Beltrami peut s'écrire en notation indicielle :

D  » D  » D  »  D  » D  » D  »


   Z   [=0
D D D   Ó D D D D D D

Il faut signaler que la notation d'Einstein aurait permis de réduire sensiblement la


longueur de cette démonstration.

On voit que cette dernière équation ne fait apparaître que des dérivées secondes des
termes du tenseur de contrainte par rapport aux coordonnées spatiales. Or ces termes
suivent une variation affine par rapport aux coordonnées spatiales, donc toutes ces
dérivées secondes sont donc nulles. Il en résulte que l'équation de Beltrami est vérifiée,
et que l'hypothèse d'un champ de contrainte affine est satisfaisante.

4. On cherche à utiliser la loi de Hooke pour démontrer que » = Ó»  » 1. On
l'utilise d'abord dans le sens contraintes-déformations :
» −Ó Ó Ó 0 0 0 ª
Ï» Ò Ï Ó −Ó Ó 0 0 0 Ò Ïª 
Ò
λ Ñ l Î Ó −Ó Ñ Îª Ñ
λ Ñ = Î Ó 0 0 0 Ñ ∙ Î  Ñ
Î  Ñ   Ó1 − Ó1 Î 0 0 0  − Ó 0 0 Ñ Îª Ñ
λ Ñ Î 0 0 0 0  − Ó 0 Ñ Îª Ñ
Í» Ð Í 0 0 0 0 0  − Ó Ð Íª Ð

Ü
Il en résulte que » = ÃÝ1zÝ1 Óª  Óª   − Ó1ª 1

On sait que ª = 0 en déformations planes. On utilise alors la loi de Hooke dans le
sens déformations-contraintes pour exprimer ª et ª :

ª  −Ó −Ó 0 0 0 »
Ϫ Ò Ï−Ó  −Ó 0 0 0 Ò Ï »
Ò
Ϊ Ñ  Î −Ó Ñ Î» Ñ
Ϊ Ñ = Î−Ó  0 0 0 Ñ ∙ Î  Ñ
Î  Ñ l Î 0 0 0 Ó 0 0 Ñ Î» Ñ
Ϊ Ñ Î0 0 0 0 Ó 0 Ñ Î» Ñ
ͪ Ð Í0 0 0 0 0   ÓÐ Í» Ð
 
On en déduit : ª = » − Ó» − Ó» 1 et ª = −Ó» » − Ó» 1
Ü Ü

52
Mécanique des Milieux Continus Polytech Grenoble, Geo3

On réinjecte ces deux expressions dans la formule de » :


l Ó Ó
» = Z » − Ó» − Ó» 1  −Ó» » − Ó» 1[
  Ó1 − Ó1 l l

Ó
» = (»  − Ó1  »  − Ó1  » −Ó1
  Ó1 − Ó1 

On isole ensuite » , et on obtient :


Ó  Ó  − Ó 1
» Z  [= »  » 1
  Ó1 − Ó1   Ó1 − Ó1 

Et finalement :
Ó  − Ó 1
  Ó1 − Ó1
» = »  » 1 = Ó »  » 1
Ó 

  Ó1 − Ó1

La dernière simplification est laissée au lecteur, mais elle a été vérifiée !


5. Les équations d'équilibre ont déjà été données à la question 2, on va donc
maintenant les utiliser en tirant parti de l'expression générale des termes du tenseur de
contraintes donnés à la question 3. L'objectif est de trouver les expressions des
constantes Ù4 ,  , 4 ,  , 4 ,  , 4 ,  Ú, et donc de définir entièrement le
tenseur de Cauchy dans le système.
La première équation d'équilibre est celle liée à la forme locale du PFD, valable pour
tout point du système :
 2»  m3 = 0
$ç% 2½&2p

On calcule d'abord la divergence du tenseur de Cauchy, donnée par :


D» D» D»
 2» =
($ç%  
 D D D

Les vecteurs 


$ç%2» et 3 s'expriment donc par :

D» D» D»


Ï   Ò
Î D  D D Ñ
Î D» D» D» Ñ 0

$ç%2» = Î   222222'2222223 = ~3
D D D Ñ
Î Ñ 0
ÎD»  D»  D» Ñ
Í D D D Ð

53
Mécanique des Milieux Continus Polytech Grenoble, Geo3

Puisque » = » = 0, on peut donc écrire le PFD sous la forme :


D» D»
  =0
 D D

D» D»
 = −m3#
 D D

 D»
=0
 D

La troisième équation est triviale (0 = 0), on peut donc tirer de ce système deux
relations pertinentes :
4   = 0 #
v
4   = −m3

La deuxième équation d'équilibre est celle qui concerne la condition limite sur la
surface immergée  , que l'on a écrite :
»− 1 = mU 3 ∙ 22½&2
 æ

Le premier terme de cette équation vaut :


» » » − »
~» » »  ∙ ~ 0  = − ~» 
» » » 0 »

On peut donc écrire :


−» = mU 3
 −» = 0 # 2½&2
−» = 0

La troisième équation est triviale (0 = 01. On remplace les termes de » par leurs
expressions dans les deux premières équations :
−4  −   = mU 3 #
v 2½&2æ
−4  −   = 0

Or sur la surface æ on peut écrire  = 0, et il vient donc :


−  = mU 3 #
v
−  = 0

Et finalement, on a les deux relations suivantes :


 = −mU 3#
v
 = 0

54
Mécanique des Milieux Continus Polytech Grenoble, Geo3

La dernière équation d'équilibre dont on dispose est celle qui concerne la surface à l'air
libre notée H , que l'on a écrite :


» ¢ 22½&2æH
 −  1£ = 0
¤

Le premier terme de cette équation s'écrit :


» » » ¤ » − »
  
» ¢  −  1£ = ~» » »  ∙ ·−Â ¸ =
¤ ~» − » 
¤ » » » ¤ » − »
0

Par conséquent, l'équation d'équilibre se traduit par


» − » = 0
» − » = 0# 22½&2æH
» − » = 0

La troisième équation de ce système est triviale une fois de plus, et les deux autres
nous permettent d'écrire :
4     − 4  −   = 0 #
v 22½&2æH
4     − 4  −   = 0

Or sur la surface æH on peut écrire  =  , et ces termes peuvent donc se simplifier
dans le système précédent. Il vient donc :
4   − 4 −  = 0 #
v
4   − 4 −  = 0

En récapitulant les trois équations d'équilibre, on obtient un système avec six équations
et six inconnues :
4   = 0
 4   = −m3

 = −mU 3 #
  = 0
4   − 4 −  = 0
4   − 4 −  = 0

En travaillant dans ce système pour calculer toutes les variables, et en tenant compte
du fait que l'on a démontré plus haut que » = Ó»  » 1, on peut définir
entièrement notre champ de contraintes en précisant les huit coefficients inconnus :
55
Mécanique des Milieux Continus Polytech Grenoble, Geo3

4 = 0 '  = −mU 3



4 = −mU 3 '  = 0 #
 4 =  m − mU 1 3 '  =  mU − m 1 3
4 = Óm − mU 13 '  = −Óm3

et donc :
» = 4     = −mU 3

» = 4     = −mU 3

» = 0

» = 4     = m − mU 13  mU − m13

» = 0

» = 4     = Ó m − mU 13 − Óm3

Le tenseur de contrainte est donc entièrement connu. Les trois valeurs non-nulles du
champ de déformation se déduisent ensuite de la loi de Hooke :
  
ª = » − Ó» − Ó» 1 =  4 − Ó4 − Ó4 1    − Ó − Ó 1
l l l
Ó Ó Ó
ª = » = 4    
l l l

ª = −Ó» » − Ó» 1
l
 
= −Ó4  4 − Ó4 1  −Ó   − Ó 1
l l

Tous calculs faits, on obtient :


−Ó3 3
ª =    Ó1m − mU 1     Ó1Óm − mU 1
l l
Ó
ª = − mU 3
l
3 3
ª =   − Ó  1m − mU 1   \mÓ  − 1  mU   Ó1]
l l

Le champ de déformation est donc lui aussi entièrement défini.

56
Mécanique des Milieux Continus Polytech Grenoble, Geo3

6. La seule contrainte tangentielle non nulle est la contrainte » , et cette contrainte
augmente (en valeur absolue) proportionnellement avec la coordonnée  . Elle est
donc maximale (en valeur absolue) au point B, où elle vaut :
» = −mU 3ℎ

57

Vous aimerez peut-être aussi