MMCexercices Corrigés

INTRODUCTION A LA MECANIQUE DES
M
MIILIEUX CONTINUS
Exercices corrigés
Guilhem MOLLON
Polytech Grenoble
Département Géotechnique, Troisième année
Edition 1, 2012-2013
V1.07
Mécanique des Milieux Continus Polytech Grenoble, Geo3
1.. Calcul tensoriel
Exercice A. Montrer que la symétrie est une propriété tensorielle, c'est-à-dire

c'est que si un
tenseur ̿ est symétrique ( ) dans une base orthonormée donné / 0 , , 1,
alors cette propriété est également vraie dans toute autre base orthonormée /′
′ , ′ , ′ .
Exercice B. Montrer que l'antisymétrie est également une propriété tensorielle.
Exercice C. Montrer qu'un tenseur quelconque peut toujours se décomposer

décompose en une
partie symétrique et une partie antisymétrique.
Exercice D. Soit ̿ un tenseur symétrique et ̿ un tenseur antisymétrique, montrer que

l'on a toujours ̿: ̿ 0.
Exercice E. Soit ̿ un tenseur symétrique et un tenseur quelconque, montrer que l'on

a toujours : ̿: ̿: , où est la partie symétrique de .
Exercice F. Soit un champ vectoriel. Montrer que l'on a toujours $%(&'

$% , 0
Exercice G. Soit un champ scalaire. Montrer que l'on a toujours &' (3&4$ , 0
Exercice H. Soit la base curviligne polaire () , * , + ,, telle que :
∙ ∙
! ∙ ∙ #
" ∙
1. Calculer la surface extérieure 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 ∙ ) sur cette sphère, sachant que le vecteur
élémentaire dans la direction radiale vaut ) - .
5 . 6
2
Corrigés
Soit un tenseur ̿ symétrique. Dans la base / 0 , , 1 on peut donc écrire

Exercice A.
= . 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) :
′ = 78 79 89
Or la matrice est symétrique, on a donc 89 = 98 . On écrit donc :
′ = 78 79 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 :
′ = 79 78 89 = 78 79 89
On en déduit que ′ = ′ , ce qu'il fallait démontrer.
Soit un tenseur ̿ antisymétrique. Dans la base / = 0 , , 1 on peut donc écrire

Exercice B.
= − . 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) :
′ = 78 79 89
Or la matrice est symétrique, on a donc 89 = −98 . On écrit donc :
′ = −78 79 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 :
′ = −79 78 89 = −78 79 89
On en déduit que ′ = −′ , ce qu'il fallait démontrer.

3
On reprend les formules du cours. Soient les tenseurs et = donnés par les formules
Exercice C.
suivantes :
1
= ( + @ ,
2
1
= = ( − @ ,
2
On va démontrer que est symétrique et que = est antisymétrique. Pour cela on se

place dans une base / = 0 , , 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 = :
1 1
= ( + , = ( + , =
2 2
1 1
= = ( − , = − ( − , = −=
2 2
On a donc démontré que est symétrique et = antisymétrique.
Dans une base orthonormée / = 0 , , 1, on peut écrire que = et également
Exercice D.
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
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é :
̿: = ̿: ( + = , = ̿: + ̿: = = ̿:
Posons le résultat intermédiaire / = &' . Dans une base orthonormée donnée

Exercice F.
/ = 0 , , 1, on a par définition :
D D
−
C D D G
D D F
/=B
B D − D F
B F
D D
−
A D D EH
Par ailleurs, on a par définition : $% / = + +

IHJ IHK IHL
I-J I-K I-L
.
On en déduit que $% &' vaut :
D D D D D D D D D
$% &' = M − N+ M − N+ M − N
D D D D D D D D D
D D D D D D
$% &' = − + − + −
D D D D D D D D D D D D
quelconque, on peut donc en déduire directement : $% &' = 0.

L'ordre des dérivations partielles successives d'une fonction de plusieurs variables est
On se place également dans une base orthonormée / = 0 , , 1. Soit le résultat

Exercice G.
intermédiaire / = 3&4$ donné par la formule du cours :
D
CD G
D F
/=B
BD F
B F
D
AD EH
5
Par ailleurs, on peut écrire, toujours dans la base / ∶
D/ D/
−
C D D G
D/ D/ F
&' / = B
B D − D F
B 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
&' 3&4$ = B
BD D −
B D D F
F
D D D D
−
AD D D D EH
formule classique &' 3&4$ = 0.

Pour la même raison que dans l'exercice précédent, on peut donc écrire directement la
1. Un élément de surface $ de la sphère de rayon délimité par deux secteurs

Exercice H.
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 \]YW $ = 2^ V $ = 2^ \− ]WY
*XY *XY
6
Finalement, on obtient le résultat classique = 4^ .
2. On doit calculer :
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
sortis de l'intégrale. On peut donc remplacer T par son expression, et écrire dans la
revanche, les trois vecteurs de la base cartésienne ont cette propriété, et peuvent être
base / = 0 , , 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 $ = 2^
+XY
7
Donc :
0
0
= 2^ ` W a
V $
*XY H
La primitive de la fonction est la fonction

bcQ L *
, donc :
0 0
C 0 G 0
= 2^ B F = `4^ a
W
d f
3 Y 3
A EH H
8
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 0g' 1 − i 0g' 1
h = i 0g' 1 + i 0g' 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′.
homogène mY à l'instant 'Y = 0. Calculer le jacobien de la transformation, ainsi que la

5. On considère un milieu animé de ce mouvement, muni d'une masse volumique
masse volumique du milieu à l'instant '.
6. Calculer le champ de vitesse n (i , ', et le champ d'accélération o (i , ', en

coordonnées lagrangiennes.
champ de vitesse n 0 , ' 1 et le champ d'accélération o 0 , ' 1 en coordonnées

7. Exprimer les coordonnées initiales à partir des coordonnées actuelles. Calculer le
eulériennes.
q 0 , ' 1 et des taux de

8. Calculer les tenseurs des taux de déformations eulériens p
q 0 , ' 1.
rotation Ω
9. On définit les coordonnées polaires lagrangiennes 0, Θ, X 1 par le changement de

variables 0i , i , X 1 = 0 ∙ Θ, ∙ Θ, X 1 et les coordonnées eulériennes
0&, , x 1 par le changement de variables 0 , , x 1 = 0& ∙ , & ∙ , x 1.
Expliciter les fonctions iT et i+ définissant une nouvelle représentation lagrangienne
du mouvement de la forme :
& = iT 0, Θ, X , t1 #
v
= i+ 0, Θ, X , t1
9
10. On définit une base curviligne polaire /x = 0T , + , 1, telle que :
T = + #
v
+ = − +
On note alors le champ de vitesse dans cette nouvelle base :
n 0 , ' 1 = nT 0&, , , ' 1T + n+ 0&, , , ' 1+ + n 0&, , , ' 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 0&, , ' 1 et l'accélération tangentielle o+ 0&, , ' 1
du mouvement étudié.
12. Définir les trajectoires associées à ce mouvement
défini dans la base orthonormée / = 0 , , 1 par la représentation eulérienne

Examen partiel : Etude cinématique d'un tourbillon. On considère un mouvement,
suivante :
n = −0 , 1
h n = 0 , 1 #
n =0
Dans cette représentation, est une fonction scalaire des deux coordonnées et ,
définie par l'expression 0 , 1 = 0 + 1z où 7 est une constante positive.
y
W
Ce mouvement est donc uniquement défini pour ‖ + ‖ > 0.
1. De quel type de mouvement s'agit-il ?
n dans la base

2. Montrer que la matrice du tenseur gradient de vitesse } = 3&4$
/ = 0 , , 1 s'exprime par :
0 , 1 2 − 0
} = ∙ ~ − −2 0
+
0 0 0H
q 0 , ' 1 et du
3. En déduire immédiatement les matrices du tenseur des taux de rotation Ω
q 0 , ' 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
masse volumique en tout point et à tout instant ' > 0.

animé de ce mouvement. Calculer la divergence du champ de vitesse, et en déduire la
10
5. Montrer que l'accélération o 0 , , , ' 1 en coordonnées eulériennes cartésiennes

(c'est-à-dire dans la base / = 0 , , 1) s'exprime par :
−
o 0, ' 1 = 0 , 1 ∙ d− f
0 H
6. On définit les coordonnées polaires lagrangiennes 0, Θ, X 1 par le changement de

variables 0i , i , X 1 = 0 ∙ Θ, ∙ Θ, X 1 et les coordonnées eulériennes
0&, , x 1 par le changement de variables 0 , , x 1 = 0& ∙ , & ∙ , x 1.
On définit également une base curviligne polaire /x = 0T , + , 1, telle que :
T = + #
v
+ = − +
Calculer l'accélération o 0&, , , ' 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+ =
2^&
mouvement sous la forme & = iT 0, Θ, , t1 et = i+ 0, Θ, , t1.

Donner sans calcul l'expression de la représentation polaire lagrangienne du
cartésienne sous la forme = (i , ', s'exprime par :

8. Démontrer que la représentation lagrangienne du mouvement dans la base
cos0g' 1 −sin0g' 1 0
= ~ sin0g' 1 cos0g' 1 0 i
0 0 1
Montrer que, dans cette expression, on a : g =

= 0 , 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
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 1H
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 1H
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 = 1
$i et $i′ vaut : o $i , $i′ =

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

=0
5. Le jacobien de la transformation est le déterminant de j . On a donc :
g' −g' 0

= $'j = g' g' 0 = g' + g' = 1
0 0 1
Par conséquent la masse volumique du milieu est constante dans le temps et en tout
point.
12
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 , ', =
- 0,1

. On a :
= i 0g' 1 − i 0g' 1
h = i 0g' 1 + i 0g' 1#
=i
p
= g ∙ (−i 0g' 1 − i 0g' 1,
Donc :
n =
p'
p
n = = g ∙ (i 0g' 1 − i 0g' 1, #
p'
p
n = =0
p'
De même, le champ d'accélération s'obtient par o (i , ', =

0,1

, donc :
pn
o = = g ∙ (−i 0g' 1 + i 0g' 1,
p'
pn
o = = g ∙ (−i 0g' 1 − i 0g' 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 1 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 1 $ −
=
$ 4$ − − 4
13
Donc :
g' −g' z 1 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 = 0g' 1 + 0g' 1
hi = − 0g' 1 + 0g' 1#
i =
question précédente, on peut alors remplacer les coordonnées initiales i par leur
Dans l'expression du champ de vitesse et du champ d'accélération calculés à la
expression en fonction des coordonnées actuelles . On obtient la formulation

eulérienne suivante :
n = g ∙ −( 0g' 1 + 0g' 1,0g' 1 − (− 0g' 1 + 0g' 1,0g' 1

#
0 1 0 1 0 1 0 1 0 1 0 1
n = g ∙ ( g' + g' , g' − (− g' + g' , g'
n =0
n = −g
n = g #
n =0
o = g ∙ −( 0g' 1 + 0g' 1,0g' 1 + (− 0g' 1 + 0g' 1,0g' 1

#
0 1 0 1 0 1 0 1 0 1 0 1
o = g ∙ −( g' + g' , g' − (− g' + g' , g'
o =0
o = −g ∙
ho = −g ∙ #
o =0
n ,n , avec I = 0 et :
directement l'accélération par o = + (3&4$
Ce dernier résultat pouvait aussi s'obtenir en passant par la formule donnant
I
I I
0 −g 0

3&4$ n = ~g 0 0
0 0 0 H
14
n, avec la formule de

8. On calcule d'abord le tenseur gradient de vitesse } = 3&4$
n , = I . On remarque qu'il s'agit d'un gradient eulérien (g minuscule),
cours (3&4$ I-
donc les dérivations sont effectuées par rapport aux coordonnées actuelles :
0 −g 0

} = 3&4$ n = ~g 0 0
0 0 0H
Ce tenseur est antisymétrique, on peut donc directement écrire p q = 0 et Ω

q = L. Ces
résultats peuvent être retrouvés par le calcul avec les formules de cours :
1 1
q=
p (} + }@ , ' q=
Ω (} − }@ ,
2 2
9. En coordonnées polaires, on a les correspondances suivantes :
i = ∙ Θ# = & ∙ θ# & = iT 0, Θ, X , t1 #

v v v
i = ∙ Θ = & ∙ θ = i+ 0, Θ, X , t1
; ;
On cherche les expressions des fonctions iT et i+ . On peut écrire à partir de ces

expressions :

& = + ; θ = ; θ =
r r
i i
= i + i ; Θ = ; Θ =
R R
Or on sait d'après l'énoncé que, dans le repère cartésien, le mouvement s'exprime par :
= i 0g' 1 − i 0g' 1
h = i 0g' 1 + i 0g' 1#
=i
On peut remplacer et par leurs expressions dans la formules donnant & :
& = + = (i 0g' 1 − i 0g' 1, + (i 0g' 1 + i 0g' 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
i 0g' 1 − i 0g' 1
θ =
r
Dans cette dernière expression, on remplace i , i et & par leur expression polaire :
∙ Θ0g' 1 − ∙ Θ0g' 1
θ = = cos0Θ + g' 1
R
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 /x = 0T , + , 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 :
n = −g ∙ + g ∙
Dans cette expression, on remplace et par leurs expressions en coordonnées
n = −g0& ∙ θ1 ∙ + g0& ∙ θ1 ∙

polaires :
n = g& 0−θ ∙ + θ ∙ 1 = g&+

donc :
Dans la base /x = 0T , + , 1, la vitesse s'exprime donc par une formule très simple :
0
n = ¢g&£
0 Hx
champ d'accélération dans la base /x . On a montré que le champ d'accélération

11. Dans le même ordre d'idée, on cherche à exprimer la formulation eulérienne du
s'exprimait par :
o = −g ∙ − g ∙
On remplace et par leurs expressions en coordonnées polaires :
16
o = −g 0& ∙ θ1 ∙ − g 0& ∙ θ1 ∙
o = −g &0θ ∙ + θ ∙ 1 = −g &T
Dans la base /x = 0T , + , 1, l'accélération est donc purement centripète :
−g &
o=¢ 0 £
0 Hx
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.
n s'exprime dans la base / = 0 , , 1

2. le tenseur gradient de vitesse } = 3&4$
par la matrice de terme général } = . On travaille d'abord sur la dérivation de la
I
I-
fonction 0 , 1 = 0 + 1z
y
W
:
D0 , 1 7 2 2

=− ∙ = −0 , 1 ∙
D 2^ 0 + 1 +
D0 , 1 7 2 2
=− ∙ = −0 , 1 ∙
D 2^ 0 + 1 +
On peut alors calculer les termes non triviaux de la matrice de }, à commencer par les
deux premiers termes diagonaux :
Dn D0 , 1 2

= − ∙ = 0 , 1 ∙
D D +
17
Dn D0 , 1 2
= ∙ = −0 , 1 ∙
D D +
On a également deux termes non-diagonaux assez complexes :
Dn D0 , 1 2
= − ∙ − 0 , 1 = 0 , 1 ∙ − 0 , 1
D D +
Dn −
= 0 , 1
D +
Dn D0 , 1 2
= ∙ + 0 , 1 = −0 , 1 ∙ + 0 , 1
D D +
Dn −
= 0 , 1
D +
Tous les autres termes de la matrice de } sont nuls car on est en mouvement plan.
base / = 0 , , 1 par :

Finalement, le tenseur gradient de vitesse s'exprime sous forme matricielle dans la
0 , 1 2 − 0
} = ∙ ~ − −2 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 : Ω
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
$% n = + +
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
Dn Dn
=−
D D
Par conséquent on a $% n = 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 0 , , , ' 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 :
pn Dn
o 0 , ' 1 = = n ,n
+ (3&4$
p' D'
On est en mouvement permanent, donc la dérivée eulérienne de la vitesse est nulle :
Dn
=0
D'
n a déjà été calculé à la question 2, il suffit donc d'en

Le gradient de vitesse } = 3&4$
faire le produit contracté avec le vecteur vitesse :
0 , 1 2 − 0 −
n ,n =
o 0, ' 1 = (3&4$ ∙ ~ − −2 0 ∙ 0 , 1 ∙ d f
+ 0 H
0 0 0H
0 , 1 − ∙ 02 1 + ∙ 0 − 1
o 0 , ' 1 = ∙ ~− ∙ 0 − 1 + ∙ 0−2 1
+
0 H
0 , 1 − − −

o 0 , ' 1 = ∙ ~− − = 0 , 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
On doit simplement remplacer et par leurs expressions polaires dans l'expression

de l'accélération. On calcule d'abord la fonction :
7 1 7 1 7
0 , 1 = = =
2^ + 2^ 0& ∙ θ1 + 0& ∙ θ1 2^&
Par ailleurs on a montré dans la question précédente que l'accélération est exprimée
par :
o 0 , ' 1 = − 0 , 1 ∙ 0 ∙ + ∙ 1
o 0 , ' 1 = − 0 , 1 ∙ 0& ∙ θ ∙ + & ∙ θ ∙ 1
Or on a la formule de changement de repère suivante : T = +

Finalement :
7
o 0, ' 1 = − 0 , 1 ∙ & ∙ T = − ∙ T
4^ &
L'accélération est donc purement centripète.
7. Le champ de vitesse est donné dans l'énoncé sous forme eulérienne par :
n = −0 , 1 ∙ + 0 , 1 ∙ = 0 , 1 ∙ 0− ∙ + ∙ 1
En coordonnées polaires, on a donc :
n = 0 , 1 ∙ 0−& ∙ θ ∙ + & ∙ θ ∙ 1 = 0 , 1 ∙ & ∙ +
7
n= ∙
2^& +
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 = Θ + '#
2^&
=i
8. Pour passer des coordonnées polaires aux coordonnées cartésiennes, on utilise les
formules suivantes :
20

& = + ; θ = ; θ =
r r
i i
= i + i ; Θ = ; Θ =
R R
On a par ailleurs = Θ + ', et donc :

y
WT K
7
θ = cos Θ + ' =
2^& r
et donc, en développant le cosinus :
7 7
= cosΘ ∙ cos ' − sinΘ ∙ sin '
r 2^& 2^&
7 7
= cosΘ ∙ cos Z '[ − sinΘ ∙ sin Z '[
r 2^0 + 1 2^0 + 1

= cosΘ ∙ cos00 , 1 ∙ ' 1 − sinΘ ∙ sin00 , 1 ∙ ' 1
r
i i
= ∙ cos00 , 1 ∙ ' 1 − ∙ sin00 , 1 ∙ ' 1
r R R
Finalement, on a :
= i ∙ cos0g' 1 − i ∙ sin0g' 1
avec g = 0 , 1. On peut développer le sinus de la même manière, et obtenir

l'expression lagrangienne de recherchée :
7
sinθ = sin Θ + ' = sin0Θ + g' 1 =
2^& r
i i
= sinΘ ∙ cos0g' 1 + cosΘ ∙ sin0g' 1 = ∙ cos0g' 1 + ∙ sin0g' 1
r R R
Finalement, on a :
= i ∙ cos0g' 1 + i ∙ sin0g' 1
On obtient finalement le système recherché :
21
= i ∙ cos0g' 1 − i ∙ sin0g' 1
h = i ∙ sin0g' 1 + i ∙ cos0g' 1#
=i
g constante pour chaque particule et qui décroit en relation hyperbolique avec la

9. Les trajectoires des particules sont planes et circulaires, avec une vitesse angulaire
distance à l'axe origine. L'accélération est purement centripète. On appelle ce

mouvement un tourbillon.
égal au tenseur gradient de vitesse, dont la matrice dans la base / = 0 , , 1 est :
10. On a montré à la question 3 que le tenseur des taux de déformation eulériens est
0 , 1 2 − 0
} = ∙ ~ − −2 q
0 = p
+
0 0 0H
On se place au point = & ∙ , et on a donc = 0 et = &. Par conséquent, la

q exprimée dans la base / en ce point est donnée par :
matrice de p
00, &1 0 & 0 7 0 1 0

q=
p ∙ ~& 0 0 = ~1 0 0
& 2^&
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 :
$ = √2 + √2
2 2
√2 √2 #
$ = −
2 2
$ =
q vaut :
Dans cette base, la matrice de p
7 1 0 0
q=
p ~0 −1 0
2^&
0 0 0 Hx
22
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 :
3 2 4

l = ~2 3 4
4 4 1H
1. Calculer la trace de l , ainsi que son déterminant. Calculer la matrice du tenseur l ,
l¥¥ = M('& l , − '& (l ,N. En déduire sans calcul le polynôme caractéristique du
et en déduire le deuxième invariant du tenseur donné par l'expression classique

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 − ,̿ = −0 − l 10 − l 10 − l 1. En déduire
les valeurs des déformations principales et la forme de la matrice de l dans sa base
principale.
propre correspondant à la déformation principale intermédiaire l .

4. Ordonner les déformations principales, puis calculer les coordonnées du vecteur
23
Corrigés
Exercice A.
1. La trace de l s'obtient immédiatement :
3 2 4

'& l '& ¢~2 3 4 £ 3 . 3 . 1 7 l¥
4 4 1H
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 :
3 2 4
2 3 4 = 3 ∗ 3 ∗ 1 + 2 ∗ 4 ∗ 4 + 4 ∗ 2 ∗ 4 − 4 ∗ 3 ∗ 4 − 4 ∗ 3 ∗ 4 − 1 ∗ 2 ∗ 2
4 4 1
$' l = −27 = l¥¥¥
Le tenseur l est le produit contracté de ª̿ par lui-même. Sous forme matricielle dans la
base / = 0 , , 1, on a :
3 2 4 3 2 4 29 28 24

l = ~2 3 4 ∙ ~2 3 4 = ~28 29 24
4 4 1H 4 4 1H 24 24 33 H
La trace de ce tenseur vaut '& l = 29 + 29 + 33 = 91
Le second invariant principal de l s'obtient directement par la formule de l'énoncé :
1 1
l¥¥ = M('& l , − '& (l ,N = 07 − 911 = −21
2 2
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¥¥¥ = − + 7 ∙ + 21 ∙ − 27
24
2. On cherche le déterminant du tenseur l − ̿ qui s'exprime sous forme matricielle

dans la base / = 0 , , 1 par :
3 2 4 1 0 0 3− 2 4
l − ̿ = ~2 3 4 − ∙ ~0 1 0 = ~ 2 3− 4
4 4 1 H 0 0 1H 4 4 1− H
Donc :
3− 2 4
$'(l − l,̿ = 2 3− 4
4 4 1−
$'(l − ,̿ = 03 − 1 ∗ 01 − 1 + 2 ∗ 4 ∗ 4 + 4 ∗ 2 ∗ 4 …
… − 4 ∗ 03 − 1 ∗ 4 − 4 ∗ 4 ∗ 03 − 1 − 01 − 1 ∗ 2 ∗ 2
$'(l − ,̿ = 0 − 6 + 9101 − 1 + 32 + 32 − 48 + 16 − 48 + 16 − 4 + 4
$'(l − ,̿ = − 6 + 9 − + 6 − 9 − 36 + 36
$'(l − ,̿ = − + 7 ∙ + 21 ∙ − 27
On retrouve donc par le calcul détaillé le résultat de la question précédente.
directement ses racines, sous la forme $'(l − ,̿ = −0 − l 10 − l 10 − l 1. Cette
3. On cherche à obtenir un développement du polynôme caractéristique faisant intervenir
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 :
3− 2 1
$'(l − ,̿ = 2 3− 3
4 4 1−
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 :
3− 2 4 03 − 1 − 2 2 4 1− 2 4
2 3− 4 = 2 − 03 − 1 3− 4 = −1 3− 4
4 4 1− 4−4 4 1− 0 4 1−
25
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 :
−1 2 4
̿
$'(l − , = 0 − 11 1 3− 4
0 4 1−
On remplace ensuite la ligne L1 par la combinaison linéaire L1+L2 :
0 5− 8
̿ 0 1
$'(l − , = − 1 1 3− 4
0 4 1−
On a le droit de permuter deux lignes (ou bien deux colonnes) en changeant le signe du
déterminant, on peut donc écrire :
1 3− 4
̿ 0 1
$'(l − , = − − 1 0 5− 8
0 4 1−
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 :
1 4
$
0 $ = ° °
§
0 §
On peut donc écrire :
5− 8
$'(l − ,̿ = −0 − 11 ± ±
4 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 − ,̿ = −0 − 11\05 − 101 − 1 − 4 ∗ 8]
$'(l − , = −0 − 110 − 6 − 271
Le discriminant du polynôme de second degré est :
∆= 0−61 − 4 ∗ 1 ∗ 0−271 = 144
26
Et les deux racines de ce polynôme sont :
6 ± √∆ √144
=3±
2 2
Finalement, le polynôme caractéristique du tenseur des déformations linéarisées est

donné par :
$'(l − ,̿ = −0 − 110 − 910 + 31
matrice de l dans sa base principale :

Les déformations principales sont les racines de ce polynôme, et on peut écrire la
1 0 0
ª̿ = ~0 9 0
0 0 −3 H´QU 8Tµb8´¶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´QU 8Tµb8´¶U
l = 9
avec :
l =1 #
l = −3
On pourrait calculer les coordonnées dans / = 0 , , 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 :
3 2 4 , , ,
~2 3 4 · , ¸ = l · , ¸ = · , ¸
4 4 1 H , H , H , H
27
Dans cette expression, on a fait apparaître les trois coordonnées (2,1 , 2,2 , 2,3 , du vecteur
dans la base / = 0 , , 1. En mettant ce système matriciel sous forme d'un
système d'équations linéaires, on a :
3 ∗ 2,1 + 2 ∗ 2,2 + 4 ∗ 2,3 = 2,1

h2 ∗ 2,1 + 3 ∗ 2,2 + 4 ∗ 2,3 = 2,2 #
4 ∗ 2,1 + 4 ∗ 2,2 + 2,3 = 2,3
2 ∗ 2,1 + 2 ∗ 2,2 + 4 ∗ 2,3 = 0

Soit :
h2 ∗ 2,1 + 2 ∗ 2,2 + 4 ∗ 2,3 = 0#

4 ∗ 2,1 + 4 ∗ 2,2 = 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 :
2,1 = √2
2
√ 2#
2,2 = −
2
2,3 = 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
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 :
= 0) est égale à j
-La composante normale de la force exercée sur la face arrière (correspondant au plan
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
-La force exercée sur la face inclinée du tétraèdre est égale à º = º + º + º
0 , , 1.
Déterminer les composantes de la matrice du tenseur de Cauchy dans la base
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). 0, 1 = 0 ; 0 , 1 = »Y ; 0 , 1 = »Y
(b). 0 , 1 = »Y ; 0 , 1 = 0 ; 0 , 1 = »Y
(c). 0, 1 = »Y ; 0 , 1 = »Y ; 0 , 1 = 0
29
1. Exprimer les trois tenseurs de contraintes »

,
0´1 »,
0¼1 et
»0b1 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 = + + .
4. Calculer le vecteur contrainte exercé sur une petite surface de direction ½ normale à
.
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 / = 0 , , 1 est donnée par
l'expression suivante, où 7 est une constante adimensionnelle :
0.77 3.67 0
» = ~3.67 2.87 0
0 0 7.6
1. Montrer que le calcul des trois contraintes 0À 1, = 1,2,3 et l'utilisation des
notera » , » , et » . Illustrer le cas 7 = 1.

propriétés des cercles de Mohr permettent de calculer les contraintes principales qu'on
2. Déterminer les valeurs de 7 correspondant à un état triaxial de révolution.
3. Pour la suite de l'exercice, on pose 7 = 1. 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 = 1, puis pour 7 = 2, la valeur de la contrainte de cisaillement

maximum ainsi que la direction de la normale correspondante.
30
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 » .
0). L'énoncé dit que la composante normale à cette force est égale à j . On peut
On va d'abord travailler sur la force appliquée sur la face arrière du tétraèdre (plan
normale sortante − .

également remarquer que la facette en question (face triangulaire arrière) est de
On va chercher à exprimer le vecteur contrainte appliqué sur cette facette. On peut

l'énoncer ainsi à partir de la formule du cours :
0− 1 = » ∙ 0− 1
1 » » »
Avec : = ~0 et » = ~» » »
0 H » » » H
Donc on a :
» » » −1 −»
0− 1 = ~» » » ∙ ~ 0 = ~−»
» » » H 0 H −» H
j . Si on note la surface de la facette triangulaire, on peut intégrer la contrainte

Par ailleurs, on sait que la composante normale de la force appliquée à la facette vaut
0− 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
−»
j = × 0− 1 ∙ 0− 1 = ~ ∙ ~ 0 = »
−» H 0 H
La surface de la facette est donnée par =

bK
, et on peut finalement écrire :
2j
» =

31
Par des raisonnements identiques sur les facettes de normales − et − , on obtient

de la même manière :
2j
» =

2j
» =

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.
1
Sa normale sortante est un vecteur unitaire, colinéaire à la direction ~1. Il s'agit donc
1
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 = 1
1
et donc :
4=
√3
1 ⁄ √3
La normale sortante à la facette est donc le vecteur = ·1⁄√3¸ .
1 ⁄ √3 H
On a également besoin de la surface ′ de cette facette. Il s'agit d'un triangle équilatéral

de côté √2, et on obtient après quelques calculs trigonométriques et/ou Pythagore :
√3
x =
2
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) :
01 = »
et donc :
32
» » » 1⁄√3
01 = ~» » » ∙ ·1⁄√3¸
» » » H
1⁄√3 H
» + » + »
1
0 1 = ~ +» +»
»
√3 » + » + »
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 :
» + » + »
′
º = º + º + º = ′ ∙ 01 = ~» + » + »
√3 » + » + »
H
Cette équation matricielle nous fournit trois équations scalaires :
′
º = 0» + » + » 1
√3

′
º = 0» + » + » 1#
√3

º = ′ 0» + » + » 1
√3

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 :
2j
º = M + » + » N
2

2j
º = M» + + » N#
2

2j
º = M» + » + N
2

On dispose alors d'un système de trois équations à trois inconnues : » , » , et » .
33
2º − 2j
» + » =

2º − 2j #
» + » =

2º − 2j
» + » =

On obtient finalement les résultats demandés suivants :
0º − j 1 + 0º − j 1 − 0º − j 1
» =

0º − j 1 − 0º − j 1 + 0º − j 1 #
» =

− º − j + º − j 1 + 0º − j 1
0 1 0
» =

Exercice B.
1. Considérons l'état de chargement (a). On nous donne trois équations vectorielles :
0 , 1 = 0 ; 0 , 1 = »Y ; 0 , 1 = »Y
D'après la première équation, on peut écrire :
» » » 1 »
0 , 1 = »0´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 »
0 , 1 = »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 »
»
0 , 1 = »0´1 ∙ = ~ » » ∙ ~0 = ~» = »Y
» » » 1 »
On en déduit que » = » = 0 et on confirme que » = »Y .
34
la base / = 0 , , 1, et on peut écrire cette matrice :

On a donc calculé les six termes indépendants de la matrice du tenseur de Cauchy dans
0 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
0 0 »Y
(c) :
»0¼1 = ~ 0 0 0
»Y 0 0 H
0 »Y 0
»0b1 = ~»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
tenseur de contraintes dont la matrice dans la base / = 0 , , 1 s'exprime par :

correspondants. Si on superpose les chargements (a), (b), et (c), on obtient donc le
0 »Y »Y
»0´1Ã0¼1Ã0b1 = »0´1 + »0¼1 + »0b1 = ~»Y 0 »Y
»Y »Y 0 H
direction . Il nous faut d'abord déterminer les coordonnées du vecteur unitaire normal
3. On cherche à exprimer le vecteur contrainte sur une facette dont on nous donne la
à cette facette. Ce vecteur est colinéaire à , mais de norme égale à 1. On peut le noter
′ et l'exprimer ainsi :
1 1
′ = = = ~ 1
‖‖ √3 √3
1 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
0 , 1 = , ′ = »0´1Ã0¼1Ã0b1 ∙ ′ = ~»Y 0 »Y ∙ ~1
»Y »Y 0 H √ 3 1 H
Et donc :
35
2»Y 1 2»Y
0 , 1 = ~1 = = 2»Y ′
√3 1 √3
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 0 , ½1, sachant que ½ ∙ = 0. Ecrivons ces deux données
sous forme mathématique :
½ ∙ = ½ + ½ + ½ = ½ + ½ + ½ = 0
0 »Y »Y ½ ½ +½
0 , ½1 = »0´1Ã0¼1Ã0b1 ∙ ½ = ~»Y 0 ½
»Y ∙ ~ = »Y ~½ + ½
»Y »Y 0 H ½ H ½ + ½ H
½
On en déduit :
0 , ½1 = −»Y ~½ = −»Y ½

½ H
5. D'après les formules précédentes, on peut formuler plusieurs observations. On a
, ′ = 2»Y ′
montré que :
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 à 2»Y .
Par ailleurs, on a montré que, pour tout vecteur ½ perpendiculaire à ′, on a :
0 , ½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 :
2»Y 0 0
»0´1Ã0¼1Ã0b1 =~ 0 −»Y 0
0 0 −»Y H´QU 8Tµb8´¶U
La base principale est formée du vecteur ′ et de tout couple de vecteurs

perpendiculaires appartenant au plan normal à ′.
36
6. En ordonnant les valeurs de contraintes (et en supposant que »Y > 0), on a :
» 0 0
»0´1Ã0¼1Ã0b1 = ~0 » 0
0 0 » H´QU 8Tµb8´¶U
avec : » = 2»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 :
» − » 3
ÄÅ´- = = »Y
2 2
1 » + » + »
»Å = '&» = =0
3 3
sphérique, et le déviateur de » est égal à ».

La contrainte moyenne est nulle, donc le tenseur de contraintes n'a pas de partie
1. On dispose de la matrice du tenseur des contraintes dans la base / = 0 , , 1 :

Exercice C.
0.77 3.67 0
» = ~3.67 2.87 0
0 0 7.6
37
Suivons l'énoncé, et calculons les trois vecteurs contraintes existant au point ¾, pour
des facettes de normales respectives , , et :
0.77 3.67 0 1 0.77

0 1 = » = ~3.67 2.87 0 ~0 = ~3.67 = 0.77 + 3.67
0 0 7.6 0 0
0.77 3.67 0 0 3.67
0 1 = » = ~3.67 2.87 0 ~1 = ~2.87 = 3.67 + 2.87
0 0 7.6 0 0
0.77 3.67 0 0 0
0 1 = » = ~3.67 2.87 0 ~0 = ~ 0 = 7.6
0 0 7.6 1 7.6
Chacun de ces vecteurs peut être décomposé en une composante normale et une
pour le vecteur contrainte appliqué à la facette de normale , on peut écrire :

composante tangentielle, afin d'être représentés dans le plan de Mohr. Par exemple,
0 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.77 et = 3.67
En appliquant la même opération à 0 1 et 0 1, on peut écrire de la même manière
µ = 2.87 et = 3.67
µ = 7.6 et = 0
Avec ces résultats, on peut tout à fait placer sur le plan de Mohr les points
coordonnées 0µ , 1, 0µ , 1, et 0µ , 1. Il faut noter pour l'instant que la
correspondant à chacun de ces trois vecteurs contraintes. Il s'agira des points de
connaissance de ces points n'est pas suffisante pour tracer le tricercle de Mohr et donc
pour déterminer les contraintes principales.
du cercle de Mohr. On remarque que la composante tangentielle de 0 1 est nulle, ce

On doit maintenant appliquer un certain nombre de raisonnements liés à des propriétés
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
aux deux vecteurs contraintes liés à ces facettes sont situées sur le même cercle de
0µ , 1 et 0µ , 1 appartiennent

Mohr (c'est la définition d'un cercle donné dans le cours). Par conséquent, les points
artiennent tous les deux au même cercle de Mohr. C'est
0µ , − 1 et 0µ , − 11.

aussi le cas de leurs symétriques par rapport à l'axe horizontal, c'est-à-dire
c'est des points
deux facettes de normales et sont orthogonales, et donc que les points qui leur
On tire enfin parti des règles d'orientation données dans le cours, et on observe
observ que les
correspondent dans le plan de Mohr sont nécessairement diamétralement opposés sur

le cercle de Mohr auquel ils appartiennent. On en déduit que les coordonné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.77 + 2.87

M , 0N = M , 0N = 01.757, 011
2 2
Le rayon de ce cercle se calcule ensuite par Pythagore :
71 + 03.671 = 702.8 − 1.751 . 03.6

= 02.87 − 1.757 61 3.757
On en déduit qu'une des contraintes principales est égale à 1.757 − 3.757 −27, et
qu'une autre contrainte principale vaut 1.757 . 3.757 5.57.. Quant à la troisième,
on l'a déjà déterminée implicitement au début de l'exercice, et elle vaut 7.6. La
situation est représentée dans le plan de Mohr pour 7 1 :
39
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 à :
» −27 ; » 5.57 ; » 7.6
On peut relever dans le cas présent trois situations pour lesquelles l'état de contraintes
est triaxial de révolution :
» » ⇒ −27 5.57 ⇒ 7 0
» = » ⇒ −27 = 7.6 ⇒ 7 = −3.8
» = » ⇒ 5.57 7.6 ⇒ 7 1.38
3. On pose 7 1, on ordonne les contraintes, et on a donc » 7.6, » 5.5, et

» −2.
On cherche les directions principales de contraintes 0 , , 1. Plus précisément, on
0 , , 1.
cherche à trouver les coordonnées de ces trois vecteur propres dans la base d'origine
On a déjà démontré que est une direction propre, et lorsque 7 = 1 cette direction
propre est celle de la contrainte principale maximale. On peut donc énoncer sans calcul
0
que :
= = ~0
1 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
0.7 3.6 0 , 0.7 , + 3.6 ,

or on a » = ~3.6 2.8
0 ~ , = ~3.6 , + 2.8 ,
0 0 7.6 H 0 H 0 H
Donc on obtient » = (0.7 , + 3.6 , , + (3.6 , + 2.8 , ,
Par ailleurs on a » = » = 5.5 ∙ Ç , + , È
On en déduit le système suivant :
5.5 = 0.7 + 3.6 , #

v
, ,
5.5 , = 3.6 , + 2.8 ,
−4.8 , + 3.6 , = 0#
v
3.6 , − 2.7 , = 0
−4.8 3.6 , 0
=
3.6 −2.7 , 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.
Il y a une information que l'on n'a pas encore utilisée : le fait que le vecteur doit être
unitaire. On a donc un nouveau système :
−4.8 , + 3.6 , = 0
v #
, + , = 1
De la première équation il vient , = , , que l'on peut intégrer dans la seconde

É
9
équation :
+ =1
,
16 ,
Finalement, on obtient : , = 0.8 et , = 0.6. On note enfin :
0.6
= ~0.8
0 H
On pourrait appliquer la même méthode pour obtenir , , et , , mais on peut aussi

tirer parti du fait que 0 , , 1 est une base orthonormée directe. On peut donc
utiliser le produit vectoriel et écrire :
41
= ∧
Donc on a :
0 0.6 −0.8
= 0 ∧ 0.8 = 0.6
~ ~ ~
1 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.7 3.6 0 √3⁄2 2.4

0 1 = » = ~3.6 2.8 0 ~ 1⁄2 = ~4.5
0 0 7.6 H 0 H 0 H
Pour tracer le point correspondant à ce vecteur contrainte dans le plan de Mohr, on a
normale s'obtient, comme d'habitude, par projection du vecteur 0 1 sur la direction
besoin d'une composante normale et d'une composante tangentielle. La composante
normale à la facette, c'est-à-dire par un produit scalaire :
2.4 √3⁄2
µ = 0 1 ∙ = ~4.5 ∙ ~ 1⁄2 = 4.3
0 H 0 H
Le vecteur contrainte 0 1 est la somme d'un vecteur µ orienté dans la direction
normale et d'un autre vecteur orienté dans une direction tangentielle ' qui n'est pas
retranchant la composante normale du vecteur 0 1 à lui-même :

connue. La composante tangentielle est donc la norme du vecteur obtenu en
= Ë0 1 − µ Ë
2.4 √3⁄2
= Ì~4.5 − 4.3 ~ 1⁄2 Ì = 2.7
0 H 0 H
On peut alors reporter le point correspondant sur le graphique (voir figure précédente).
facette de direction normale = + qui appartient au plan des deux

On constate qu'il appartient au cercle 1. C'est logique, car ce point est défini pour une
√
directions principales 0 , 1 (c'est le plan normal à ).
42
5. Pour 7 = 1, on a » = 7.6 et » = −2 . La contrainte de cisaillement maximale est

donc donnée par :
» − »
ÄÅ´- = = 4.8
2
L'abscisse du point correspondant sur le plan de Mohr correspond à celle du centre du

grand cercle de Mohr, soit:
» + »
= 2.8
2
donc de coordonnées 02.8,4.81 (voir figure).

Le point du cercle de Mohr correspondant à la facette de cisaillement maximal est
au plan des deux directions principales et (qu'on appelle "plan de cisaillement

Par définition, la facette qui reçoit cette contrainte possède une normale qui appartient
orientation situé selon la bissectrice des directions principales et , qui est donc
maximal"). En fait, la facette de cisaillement maximal correspond toujours à une
donnée par :
0 − 0.8
Ï Ò
Î 2 Ñ
+ 0 + 0.6Ñ −0.4
=Î = ~ 0.3
2 Î 2 Ñ 0.5
Î 1+0 Ñ
Í 2 Ð
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 2, 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
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 let Ó.
-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.
-Champ de déplacement donné par Ô = 0 ∙ i + k ∙ i 1 + 0−k ∙ i + / ∙ i 1 .

-On se place dans la zone centrale du massif pour se prémunir des effets de bord.
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).
paramètres et /. En utilisant les conditions aux limites en déplacement, calculer k.

4. En utilisant les conditions aux limites en contraintes, calculer les expressions des
5. On pose l = 210000 ¾Õ4 et Ó = 0.3 (Massif en acier). Les dimensions de la

section sont } = 0.2Ö et × = 0.15Ö, et le chargement vaut : = 500 ¾Õ4. Calculer
le vecteur déplacement du point ¾. Tracer la déformée de la section.
44
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
d'ancrage latéraux. Le barrage est soumis à son poids m¼ 3 et à la poussée mU 3 de

particulier à une "tranche" de ce barragebarrage située suffisamment loin des massifs
ci étant située à une hauteur ℎ égale à celle du barrage. A sa

l'eau, la surface de celle-ci
base, le barrage est parfaitement encastré.
paramètres let Ó.. On néglige la pression atmosphérique.

Le matériau composant le barrage est supposé élastique linéaire et isotrope, de
forme générale des inconnues Ô, ª̿, et ».

1. Expliquer pourquoi il s'agit d'un problème de déformations planes, et donner la
2. Décrire précisément les conditions aux limites d'une "tranche" du barrage. Ecrire les
équations d'équilibre mécanique, d'abord
d'abord à l'intérieur du domaine puis au niveau de ses
conditions limites en contraintes.
des coordonnées et . Donner l'expression générale des termes de la matrice d'un

3. On postule que tous les termes du tenseur de contrainte sont des fonctions linéaires
tel champ de contrainte. Montrer qu'il vérifie l'équation de Beltrami, et est donc la
solution du problème élastique.
45
Rappel : l'équation de Beltrami est donnée par :
1
∆» +
3&4$ 3&4$ '& » = 0
1+Ó
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 0∆»1 = ∆(» ,. Les
autres termes de l'équation sont calculables.
4. A l'aide de la loi de Hooke, montrer que » = Ó0» + » 1.
calculer » en tout point du système. A l'aide de la loi de Hooke, en déduire ª̿ en tout

5. A l'aide des équations d'équilibre et des conditions aux limites en contraintes,
point du système.
est obtenue au point /, et calculer cette contrainte.

6. Montrer que la contrainte tangentielle la plus intense (en valeur absolue) du système
46
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 :
1
ª̿ = Ô + (3&4$
3&4$ Ô,
@
2
IØJ IØJ IØJ
ÏI-J I-L Ò
ÎIØK
I-K
IØK Ñ
k 0
Or
3&4$ Ô = ÎI- = ~−k / 0
IØK
I-K I-L Ñ
ÎIØL IØL Ñ 0 0 0
J
IØL
ÍI-J I-K I-L Ð
47
0 0
Finalement, ª̿ = ~ 0 / 0
0 0 0
est la base d'origine 0 , , 1, et que les déformations principales sont Ù, /, 0Ú.
On constate qu'il s'agit d'une déformation plane, que la base principale de déformation
3. On applique la loi de Hooke formulée en fonction des coefficients de Lamé :
» = 2Û ∙ ª̿ + ∙ '&0ª̿1 ∙ ̿
avec Û = et = 0ÃÝ10z
Ü ÜÝ
0ÃÝ1 Ý1
.
On a '&0ª̿1 = + /, et il vient donc :
0 0 1 0 0
» = 2Û ∙ ~ 0 / 0 + ∙ 0 + / 1 ∙ ~0 1 0
0 0 0 0 0 1
2Û + ∙ 0 + /1 0 0
Et finalement :
» = · 0 2Û/ + ∙ 0 + /1 0 ¸
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 »0− 1 = Þ
Par conséquent, on peut écrire :−(2Û + ∙ 0 + /1, ∙ = −: ∙
Et finalement : 2Û + ∙ 0 + /1 = :
On applique la même méthode sur une des faces latérales, et on obtient :
2Û/ + ∙ 0 + /1 = 0
Il en résulte un système linéaire de deux équations à deux inconnues et / :
48
02Û + 1 + B = :#
v
A + 02Û + 1/ = 0
La deuxième équation donne / =

záâ
ãÃá
, que l'on réinjecte dans la première pour
obtenir 2Û + − = :.
áK
ãÃá
Après calcul, il vient :
−2Û −
=: ' / = :
4Û0Û + 1 4Û0Û + 1
Sur la face inférieure, le déplacement est bloqué verticalement. En posant i = 0, on
Ô 00, i , i 1 = 0
peut donc écrire :
Or Ô 0i , i , i 1 = ∙ i + k ∙ i .
On en déduit que k ∙ i = 0 ∀0i , 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 :
−2Û −
∙ i = : ∙i
4Û0Û + 1
Ô= #
/ ∙ i = : ∙i
4Û0Û + 1
0
Les coordonnées du point ¾ sont i = × = 0.15Ö et i = } = 0.2Ö.
Par ailleurs, les coefficients de Lamé valent :
l 210000
Û= = = 80770 ¾Õ4
201 + Ó1 201 + 0.31
210000 ∗ 0.3
= = 121150 ¾Õ4
01 + 0.3101 − 2 ∗ 0.31
Il vient finalement = −2.1666 ∙ 10z et / = 9.2855 ∙ 10zÉ , et enfin :
−0.325ÖÖ
Ô0¾1 = 0.186ÖÖ #
0
49
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 :
Ô 0i , i 1 ª ª 0 » » 0

Ô = Ô 0i , i 1# ; ª̿ = ~ª ª 0 ; » = ~» » 0
0 0 0 0 0 0 »
utiliser indifféremment les vecteurs position initial i et final .

Par ailleurs, en élasticité linéaire on applique systématiquement l'HPP, et on peut
-Sur sa base il est parfaitement encastré : Ô = 0 ½& å=H

2. Sur une "tranche" du barrage, les conditions aux limites sont les suivantes :
-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
$ç% » + m3 = 0 ½& p
du solide chargé :
-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 :
»0− 1 = mU 3 ∙ ½& æ=
-Sur la surface à l'air libre, la normale sortante est = 0 − 1, et la

√
contrainte subie par la surface est nulle :
1
» ¢ 0 − 1£ = 0 ½& æH
√2
linéaires des coordonnées et , on a donc proportionnalité entre chaque terme de

3. L'énoncé postule que tous les termes du tenseur de contrainte sont des fonctions
ce tenseur et ces coordonnées. On peut donc écrire les relations suivantes :
» = 4 +

» = 4 +
» = 4 +
» = 4 +
50
Dans ces expressions, les termes Ù4 , , 4 , , 4 , , 4 , Ú sont des

constantes de proportionnalité qui sont pour l'instant inconnues.
sont nulles car on est en déformations planes : » = » = 0.

Quant aux deux autres composantes du tenseur de Cauchy, on a déjà démontré qu'elles
Pour que ce champ de contrainte soit admissible, il doit vérifier l'équation de Beltrami:
1
∆» +
3&4$ 3&4$ '& » = 0
1+Ó
On écrit d'abord les différents termes de cette équation sous forme indicielle. Le
premier terme permet d'écrire :
D » D » D »

0∆»1 = ∆(» , = + +
D D D
le deuxième terme est légèrement plus long à dériver. On développe d'abord la trace :
3&4$ '& », = 3&4$

(3&4$ 3&4$ 0» + » + » 1

On développe ensuite le gradient du champ scalaire de cette trace, qui est donc un
champ vectoriel :
D
Ï 0» + » + » 1Ò
C Î D ÑG
Î D ÑF
3&4$ '& », = B
(3&4$ B3&4$ ÎD 0» + » + » 1ÑF
B ÑF

Î
Î D
0» + » + » 1Ñ
A ÍD ÐE
On calcule enfin le gradient de ce champ vectoriel, qui est donc un champ tensoriel :
D D D
Ï 0» + » + » 1 0» + » + » 1 0» + » + » 1Ò
Î D D D D D Ñ
Î D D D Ñ
3&4$ '& », =
(3&4$ ÎD D 0» + » + » 1 D
0» + » + » 1
D D
0» + » + » 1Ñ
Î Ñ
Î D 0» + » + » 1
D
0» + » + » 1
D
0» + » + » 1 Ñ
ÍD D D D D Ð
51
Il vient :
D D » D » D »
3&4$ '& », =
(3&4$ 0» + » + » 1 = + +
D D D D D D D D
Finalement, l'équation de Beltrami peut s'écrire en notation indicielle :
D » D » D » 1 D » D » D »

+ + + Z + + [=0
D D D 1 + Ó 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 » = Ó0» + » 1. On

l'utilise d'abord dans le sens contraintes-déformations :
» 1−Ó Ó Ó 0 0 0 ª
Ï» Ò Ï Ó 1−Ó Ó 0 0 0 Ò Ïª Ò
Î» Ñ l Î Ó 1−Ó Ñ Îª Ñ
Î» Ñ = Î Ó 0 0 0 Ñ∙Î Ñ
Î Ñ
01 + Ó101 − 2Ó1 Î 0 0 0 1 − 2Ó 0 0 Ñ Îª Ñ
Î» Ñ Î 0 0 0 0 1 − 2Ó 0 Ñ Îª Ñ
Í» Ð Í 0 0 0 0 0 1 − 2Ó Ð Íª Ð
Il en résulte que » = 0ÃÝ10z 0Óª + Óª + 01 − Ó 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 ª :
ª 1 −Ó −Ó 0 0 0 »
Ïª Ò Ï−Ó 1 −Ó 0 0 0 Ò Ï » Ò
Îª Ñ 1 Î −Ó Ñ Î» Ñ
Îª Ñ = Î−Ó 1 0 0 0 Ñ∙Î Ñ
Î Ñ lÎ 0 0 0 1+Ó 0 0 Ñ Î» Ñ
Îª Ñ Î0 0 0 0 1+Ó 0 Ñ Î» Ñ
Íª Ð Í0 0 0 0 0 1 + Ó Ð Í» Ð
On en déduit : ª = 0» − Ó» − Ó» 1 et ª = 0−Ó» +» − Ó» 1

Ü Ü
52
On réinjecte ces deux expressions dans la formule de » :
l Ó Ó
» = Z 0» − Ó» − Ó» 1 + 0−Ó» +» − Ó» 1[
01 + Ó101 − 2Ó1 l l
Ó
» = (» 01 − Ó1 + » 01 − Ó1 + » 0−2Ó1,
01 + Ó101 − 2Ó1
On isole ensuite » , et on obtient :
2Ó Ó 01 − Ó 1
» Z1 + [= 0» + » 1
01 + Ó101 − 2Ó1 01 + Ó101 − 2Ó1
Ó 01 − Ó 1
Et finalement :
01 + Ó101 − 2Ó1
» = 0» + » 1 = Ó 0» + » 1
2Ó
1+0
1 + Ó101 − 2Ó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
constantes Ù4 , , 4 , , 4 , , 4 , Ú, et donc de définir entièrement le

contraintes donnés à la question 3. L'objectif est de trouver les expressions des
tenseur de Cauchy dans le système.
La première équation d'équilibre est celle liée à la forme locale du PFD, valable pour
$ç% » + m3 = 0 ½& p
tout point du système :
On calcule d'abord la divergence du tenseur de Cauchy, donnée par :
D» D» D»

($ç% », = + +
D D D
Les vecteurs $ç% » et 3 s'expriment donc par :
D» D» D»

Ï + + Ò
Î D D D Ñ
0
ÎD» D» D» Ñ
$ç% » = Î + + ' 3 = ~3
D D D Ñ
Î Ñ 0
ÎD» + D» + D» Ñ
Í D D D Ð
53
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
surface immergée = , que l'on a écrite :

La deuxième équation d'équilibre est celle qui concerne la condition limite sur la
»0− 1 = mU 3 ∙ ½& æ=
Le premier terme de cette équation vaut :
» » » −1 »
~» » » ∙ ~ 0 = − ~»
» » » 0 »
−» = mU 3
On peut donc écrire :
−» = 0 # ½& =

−» = 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 ½& æ=
−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
libre notée H , que l'on a écrite :

La dernière équation d'équilibre dont on dispose est celle qui concerne la surface à l'air
1
» ¢ 0 − 1£ = 0 ½& æH
√2
Le premier terme de cette équation s'écrit :
» » » 1⁄√2 » − »
1 1
» ¢ 0 − 1£ = ~» » » ∙ ·−1⁄ ¸ =
√2 ~» − »
√2 » » » √2 » − »
0
Par conséquent, l'équation d'équilibre se traduit par
» − » = 0
− » = 0# ½& æ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 ½& æ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
du fait que l'on a démontré plus haut que » = Ó0» + » 1, on peut définir
En travaillant dans ce système pour calculer toutes les variables, et en tenant compte
entièrement notre champ de contraintes en précisant les huit coefficients inconnus :
55
4 = 0 ' = −mU 3

4 = −mU 3 ' = 0 #
4 = 0 m − 2mU 1 3 ' = 0 mU − m 1 3
4 = Ó0m − 2mU 13 ' = −Óm3
et donc :
» = 4 + = −mU 3
» = 4 + = −mU 3
» = 0
» = 4 + = 0m − 2mU 13 + 0mU − m13
» =0
» = 4 + = Ó 0m − 2mU 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
ª = 0» − Ó» − Ó» 1 = 04 − Ó4 − Ó4 1 + 0 − Ó − Ó 1
l l l
1+Ó 1+Ó 1+Ó
ª = » = 4 +
l l l
1
ª = 0−Ó» +» − Ó» 1
l

= 0−Ó4 + 4 − Ó4 1 + 0−Ó + − Ó 1
l l
Tous calculs faits, on obtient :
−Ó3 3
ª = 01 + Ó10m − 2mU 1 + 01 + Ó10Óm − mU 1
l l
1+Ó
ª = − mU 3
l
3 3
ª = 01 − Ó 10m − 2mU 1 + \m0Ó − 11 + mU 01 + Ó1]
l l
Le champ de déformation est donc lui aussi entièrement défini.
56
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

MMCexercices Corrigés

Transféré par

Droits d'auteur :

Formats disponibles

MMCexercices Corrigés

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

MMCexercices Corrigés

Transféré par

Droits d'auteur :

Formats disponibles

INTRODUCTION A LA MECANIQUE DES

1.. Calcul tensoriel

Exercice A. Montrer que la symétrie est une propriété tensorielle, c'est-à-dire

Exercice B. Montrer que l'antisymétrie est également une propriété tensorielle.

Exercice C. Montrer qu'un tenseur  quelconque peut toujours se décomposer

Exercice D. Soit ̿ un tenseur symétrique et ̿ un tenseur antisymétrique, montrer que

Exercice E. Soit ̿ un tenseur symétrique et  un tenseur quelconque, montrer que l'on

Exercice F. Soit  un champ vectoriel. Montrer que l'on a toujours $%(&'

Exercice H. Soit la base curviligne polaire () , * , + ,, telle que :

Soit un tenseur ̿ symétrique. Dans la base /  0 ,  ,  1 on peut donc écrire

= . Soit une base orthonormée quelconque /′ = ′ , ′ , ′  différente de /.

′ = 78 79 89

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

′ = 78 79 98

′ = 79 78 89 = 78 79 89

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

Soit un tenseur ̿ antisymétrique. Dans la base / = 0 ,  ,  1 on peut donc écrire

= − . Soit une base orthonormée quelconque /′ = ′ , ′ , ′  différente de /.

′ = 78 79 89

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

′ = −78 79 98

′ = −79 78 89 = −78 79 89

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

On va démontrer que   est symétrique et que  = est antisymétrique. Pour cela on se

On a donc démontré que   est symétrique et  = antisymétrique.

̿: ̿ = 

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.

̿:  = ̿: (  +  = , = ̿:   + ̿:  = = ̿:  

Posons le résultat intermédiaire / = &' . Dans une base orthonormée donnée

/ = 0 ,  ,  1, on a par définition :

Par ailleurs, on a par définition : $% / = + +

On en déduit que $% &'  vaut :

quelconque, on peut donc en déduire directement : $% &'  = 0.

On se place également dans une base orthonormée / = 0 ,  ,  1. Soit le résultat

intermédiaire / = 3&4$  donné par la formule du cours :

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

formule classique &' 3&4$  = 0.

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 ,

Le rayon  est indépendant de  et , et le terme  est indépendant de , on peut

Finalement, on obtient le résultat classique = 4^ .

Avec : ) =  - +  5 +  6

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

La primitive de la fonction   est la fonction

Problème. On considère un mouvement défini dans la base /  0 ,  ,  1 par sa

1. Calculer le tenseur gradient j , le tenseur des dilatations k̿ , et le tenseur des

2. A quelle classe particulière ce mouvement appartient-il ?

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 n (i , ', et le champ d'accélération o (i , ', en

champ de vitesse n 0 , ' 1 et le champ d'accélération o 0 , ' 1 en coordonnées

q 0 , ' 1 et des taux de

9. On définit les coordonnées polaires lagrangiennes 0, Θ, X 1 par le changement de

10. On définit une base curviligne polaire /x = 0T , + ,  1, telle que :

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

Exercice C. Montrer qu'un tenseur quelconque peut toujours se décomposer

Exercice D. Soit ̿ un tenseur symétrique et ̿ un tenseur antisymétrique, montrer que

Exercice E. Soit ̿ un tenseur symétrique et un tenseur quelconque, montrer que l'on

Exercice F. Soit un champ vectoriel. Montrer que l'on a toujours $%(&'

Exercice H. Soit la base curviligne polaire () , * , + ,, telle que :

Soit un tenseur ̿ symétrique. Dans la base / 0 , , 1 on peut donc écrire

= . Soit une base orthonormée quelconque /′ = ′ , ′ , ′ différente de /.

Soit un tenseur ̿ antisymétrique. Dans la base / = 0 , , 1 on peut donc écrire

= − . Soit une base orthonormée quelconque /′ = ′ , ′ , ′ différente de /.

On va démontrer que est symétrique et que = est antisymétrique. Pour cela on se

On a donc démontré que est symétrique et = antisymétrique.

̿: ̿ =

̿: ̿ = + +

On sait que l'on a = = = 0, et par ailleurs :

On en déduit directement ̿: ̿ = 0.

̿: = ̿: ( + = , = ̿: + ̿: = = ̿:

Posons le résultat intermédiaire / = &' . Dans une base orthonormée donnée

/ = 0 , , 1, on a par définition :

Par ailleurs, on a par définition : $% / = + +

On en déduit que $% &' vaut :

quelconque, on peut donc en déduire directement : $% &' = 0.

On se place également dans une base orthonormée / = 0 , , 1. Soit le résultat

intermédiaire / = 3&4$ donné par la formule du cours :

formule classique &' 3&4$ = 0.

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 ,

Le rayon est indépendant de et , et le terme est indépendant de , on peut

Finalement, on obtient le résultat classique = 4^ .

Avec : ) = - + 5 + 6

La primitive de la fonction est la fonction

Problème. On considère un mouvement défini dans la base / 0 , , 1 par sa

1. Calculer le tenseur gradient j , le tenseur des dilatations k̿ , et le tenseur des

champ de vitesse n 0 , ' 1 et le champ d'accélération o 0 , ' 1 en coordonnées

q 0 , ' 1 et des taux de

9. On définit les coordonnées polaires lagrangiennes 0, Θ, X 1 par le changement de

10. On définit une base curviligne polaire /x = 0T , + , 1, telle que :

n 0 , ' 1 = nT 0&, , , ' 1T + n+ 0&, , , ' 1+ + n 0&, , , ' 1

défini dans la base orthonormée / = 0 , , 1 par la représentation eulérienne

n dans la base

5. Montrer que l'accélération o 0 , , , ' 1 en coordonnées eulériennes cartésiennes

6. On définit les coordonnées polaires lagrangiennes 0, Θ, X 1 par le changement de

On définit également une base curviligne polaire /x = 0T , + , 1, telle que :

Calculer l'accélération o 0&, , , ' 1 en coordonnées eulériennes polaires.

mouvement sous la forme & = iT 0, Θ, , t1 et = i+ 0, Θ, , t1.

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

Montrer que, dans cette expression, on a : g =

10. On considère le point de vecteur position = & ∙ . Calculer la matrice du tenseur

g' + g' g' ∙ g' − g' ∙ g' 0

Enfin, on sait que l = (k̿ − ,̿ , et on a donc l = 0.

$i et $i′ vaut : o $i , $i′ =

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

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

g' −g' z 1 g' g' g' −g' @

expression en fonction des coordonnées actuelles . On obtient la formulation

n = g ∙ −( 0g' 1 + 0g' 1,0g' 1 − (− 0g' 1 + 0g' 1,0g' 1

o = g ∙ −( 0g' 1 + 0g' 1,0g' 1 + (− 0g' 1 + 0g' 1,0g' 1

n, avec la formule de

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

i = ∙ Θ# = & ∙ θ# & = iT 0, Θ, X , t1 #

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

& = + = (i 0g' 1 − i 0g' 1, + (i 0g' 1 + i 0g' 1,

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

n = −g0& ∙ θ1 ∙ + g0& ∙ θ1 ∙

n = g& 0−θ ∙ + θ ∙ 1 = g&+

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

o = −g 0& ∙ θ1 ∙ − g 0& ∙ θ1 ∙

o = −g &0θ ∙ + θ ∙ 1 = −g &T

Dans la base /x = 0T , + , 1, l'accélération est donc purement centripète :

n s'exprime dans la base / = 0 , , 1

D0 , 1 7 2 2

Dn D0 , 1 2

base / = 0 , , 1 par :

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

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

5. L'accélération o 0 , , , ' 1 est la dérivée particulaire de la vitesse, et se calcule

n a déjà été calculé à la question 2, il suffit donc d'en

0 , 1 − − −

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

o 0 , ' 1 = − 0 , 1 ∙ 0 ∙ + ∙ 1