Mécanique des

Milieux Continus
Golay Frédéric

SEATECH
MMC

Golay MMC -2-

Ce cours de mécanique des milieux continus est à la base de l’enseignement de mécanique à SEATECH. Les
notions abordées ici, transport de champs, lois de conservation, ..., seront reprises ultérieurement en
mécanique des solides et mécanique des fluides. Dans une première partie, nous aborderons les notations
tensorielles et vectorielles indispensables à toute étude scientifique, puis dans une deuxième partie, nous
étudierons la cinématique des milieux continus. Après avoir introduit la modélisation des efforts et les lois de
conservation par le principe des puissances virtuelles, nous appliquerons ces lois de conservation aux lois de
comportement de l’élasticité linéaire (en mécanique des solides) et aux lois de comportement des fluides
newtoniens (en mécanique des fluides).

-3- Golay MMC

MMC

Golay MMC -4-

 Sommaire

TABLE DES MATIÈRES

Notations tensorielles ....................................................................................................... 9

1 Vecteurs et tenseurs ............................................................................................... 9

1.1 Notations ............................................................................................................................................... 9
1.2 Changement de repère ........................................................................................................................ 12

2 Permutations et déterminants............................................................................... 14
2.1 Les symboles de permutation .............................................................................................................. 14
2.2 Déterminant d’une matrice ................................................................................................................. 14
2.3 Polynôme caractéristique .................................................................................................................... 15
2.4 Adjoint d’un tenseur antisymétrique ................................................................................................... 15

3 Calcul vectoriel et analyse vectorielle .................................................................... 16

3.1 Calcul vectoriel ..................................................................................................................................... 16
3.2 Analyse vectorielle ............................................................................................................................... 16
3.3 Transformation d’intégrales ................................................................................................................ 17

4 Formules essentielles en Mécanique des Milieux Continus .................................... 18

4.1 Coordonnées cartésiennes orthonormées .......................................................................................... 18
4.2 Coordonnées cylindriques ................................................................................................................... 19
4.3 Coordonnées sphériques ..................................................................................................................... 20
4.4 Comment retrouver les formules ........................................................................................................ 21

5 A retenir ............................................................................................................... 23

CINEMATIQUE ................................................................................................................. 25

1 Le mouvement et ses représentations ................................................................... 25

1.1 Configuration ....................................................................................................................................... 25
1.2 Variables de Lagrange et variables d’Euler .......................................................................................... 26
1.3 Dérivées particulaires .......................................................................................................................... 26

2 Déformation d’un milieux continu ......................................................................... 27

2.1 Notion de déformation ........................................................................................................................ 27
2.2 Tenseur des déformations ................................................................................................................... 28
2.3 Conditions de compatibilité ................................................................................................................. 30

3 Transport, dérivées particulaires ........................................................................... 30

3.1 Transport d’un volume ........................................................................................................................ 30
3.2 Transport d’une surface orientée ........................................................................................................ 31
3.3 Dérivée particulaire d’une intégrale de volume .................................................................................. 32
3.4 Dérivée particulaire d’une intégrale de surface .................................................................................. 33

4 A retenir ............................................................................................................... 35

EFFORTS DANS LES MILIEUX CONTINUS ........................................................................... 37

-5- Golay MMC
MMC

1 Définitions ............................................................................................................ 37
1.1 Forces ................................................................................................................................................... 37
1.2 Vecteur-contrainte et tenseur des contraintes .................................................................................... 37

2 Equilibre ............................................................................................................... 39
2.1 Le Principe des Puissances Virtuelles (Germain 1972) ......................................................................... 39
2.2 Puissance virtuelle des efforts intérieurs ............................................................................................. 39
2.3 Puissance virtuelle des efforts extérieurs ............................................................................................ 40
2.4 Application du Principe des Puissances Virtuelles ............................................................................... 40
2.5 Equilibre ............................................................................................................................................... 41
2.6 Autre présentation: Principe fondamental de la dynamique............................................................... 42

3 Quelques propriétés du tenseur des contraintes ................................................... 43

3.1 Symétrie du tenseur des contraintes ................................................................................................... 43
3.2 Contrainte normale et contrainte tangentielle .................................................................................... 44
3.3 Directions principales, contraintes principales .................................................................................... 44
3.4 Invariants .............................................................................................................................................. 44
3.5 Cercles de Mohr ................................................................................................................................... 44

4 Exemples de tenseur des contraintes .................................................................... 47

4.1 Tenseur uniaxial ................................................................................................................................... 47
4.2 Tenseur sphérique................................................................................................................................ 47

5 A retenir ............................................................................................................... 48

ELASTICITE ...................................................................................................................... 49

1 Approche expérimentale: essai de traction............................................................ 49

2 Loi de comportement élastique linéaire (en HPP) .................................................. 50

2.1 Forme générale .................................................................................................................................... 50
2.2 Matériau élastique homogène isotrope............................................................................................... 50
2.3 Matériau élastique homogène orthotrope .......................................................................................... 50
2.4 Matériau élastique homogène isotrope transverse ............................................................................. 51
2.5 Caractéristiques de quelques matériaux .............................................................................................. 51
2.6 Critères de limite d’élasticité ............................................................................................................... 52

3 Le problème d’élasticité ........................................................................................ 53

3.1 Ecriture générale .................................................................................................................................. 53
3.2 Formulation en déplacement ............................................................................................................... 53
3.3 Formulation en contrainte ................................................................................................................... 53
3.4 Théorème de superposition ................................................................................................................. 53
3.5 Elasticité plane ..................................................................................................................................... 54
3.6 Thermoélasticité .................................................................................................................................. 55

4 A retenir ............................................................................................................... 58

INTRODUCTION A LA MECANIQUE DES FLUIDES............................................................... 59

1 Loi de comportement ............................................................................................ 59

1.1 Fluide Newtonien ................................................................................................................................. 59
1.2 Fluide incompressible........................................................................................................................... 60
1.3 Fluide non-visqueux ............................................................................................................................. 60
1.4 Fluide au repos ..................................................................................................................................... 60

Golay MMC -6-

 Sommaire

2 Conservation de la masse ...................................................................................... 60

3 Equation du mouvement ....................................................................................... 61

4 A retenir ............................................................................................................... 62

Bibliographie ................................................................................................................... 63

Annexes: Rappels de mécaniques des solides rigides ....................................................... 65

1 Cinématiques du solide ......................................................................................... 65

1.1 Description du mouvement ................................................................................................................. 65
1.2 Composition des mouvements ............................................................................................................ 66

2 Cinétique .............................................................................................................. 68
2.1 Définitions ............................................................................................................................................ 68
2.2 Eléments de cinétique ......................................................................................................................... 68
2.3 Cinétique du solide rigide .................................................................................................................... 69

3 Equations fondamentales de la mécanique des solides .......................................... 72

3.1 Torseur associé aux efforts externes ................................................................................................... 72
3.2 Loi fondamentale de la dynamique ..................................................................................................... 72

-7- Golay MMC

MMC

Golay MMC -8-

 Notations tensorielles

NOTATIONS TENSORIELLES

1 Vecteurs et tenseurs
Avertissement : L’objectif de ce chapitre, est de familiariser les étudiants avec les notations tensorielles. Afin
d’en simplifier le contenu, nous ne considérerons que des bases orthonormées.

1.1 Notations

1.1.1 Vecteur
   
Dans un espace euclidien  à trois dimensions, soit e1, e2 , e3 une base orthonormée. Un vecteur V est

représenté par ses composantes V1, V2 , V3

    3 
V  V1e1 V2e2 V3e3  Viei
i 1 (1.1)
En utilisant la convention de sommation, ou convention d’Einstein, on écrit
 
V  Viei
(1.2)

où, chaque fois qu’un indice est répété, il convient de faire varier cet indice de 1 à 3 et de faire la somme. Dans
l’expression (2) l’indice i est un "indice muet".

En notation matricielle on écrira parfois


 

 
V1 
 
   



 

 
V  V  V2 

 

(1.3)
 
V 

 
3 
 


et le vecteur transposé
 
T

T

V  V  V  V1 V2 V3
(1.4)

1.1.2 Application linéaire de  dans 

  
Soit A une application linéaire, dans la base e1, e2 , e3 . Cette application est représentée par une matrice 3x3

notée A :
 
A A A 
 11 12 13 
A A A 
 21 22 23 
 
A
 31 A32
A33 

   

Si W est un vecteur tel que W  AV , alors les composantes de W sont données par

W1  A11V1  A12V2  A13V3

W2  A21V1  A22V2  A23V3
W3  A31V1  A32V2  A33V3

et en utilisant les conventions de sommation où j est un indice muet

-9- Golay MMC

MMC

Wi  AijVj
(1.5)
et en notation vectorielle

W   A V 
On définit les symboles de Kronecker par

1 si ij
ij  


0 si ij (1.6)



En particulier l’application identité 1 est représentée par la matrice




12 13 

 1 0 0
 11
  


 22

23    0 1 0
 21
  
   


 31 32 33   0 0 1 
 

La composition de deux applications linéaires se traduit par le produit de leur matrice représentative, c’est-à-
dire

C  AB ou encore C   A B 

     
et en notation indicielle

C ij  Aik Bkj
(1.7)

1.1.3 Formes bilinéaires

Soit A une forme bilinéaire sur  , c’est-à-dire une application bilinéaire de    dans  . Dans la base
  
e1, e2 , e3 elle est représentée par une matrice Aij telle que
 
 
A V W  AijVW
i j (1.8)

ou en notation matricielle
 
 
A V W  V A W 
 

En particulier, la forme bilinéaire représentée dans toute base par les symboles de Kronecker est le produit
  
scalaire. Si ( e1, e2 , e3 ) est une base orthonormée, alors
 
ei  e j  ij

et le produit scalaire de deux vecteurs est donné par

     
V W  Viei Wje j  VW i
e  e j  ijVW
j i i j
 VW
i i

ou en notation matricielle
 
V W  V W 

1.1.4 Tenseurs

[Link] Tenseur du second ordre


Un tenseur du second ordre T est un opérateur linéaire qui fait correspondre à tout vecteur V de l’espace

euclidien un vecteur W de ce même espace.

Golay MMC - 10 -
 Notations tensorielles

 
W T V 
 
Cet opérateur peut être représenté par une matrice 3x3, notée T  ou T  ou T , telle que
   

Wi  TijVj

ou en notation matricielle
 
W   T  V 
ou
 
W  TV

* Un tenseur est dit symétrique si Tij  Tji

* Un tenseur est dit antisymétrique si Tij  Tji

* Un tenseur est dit isotrope si Tij  t ij

* On peut toujours décomposer un tenseur en une partie symétrique et antisymétrique

S A
T  T T Tij  TijS  TijA
ou

1 1
TijS 
2

Tij  Tji 
TijA  Tij  Tji
2
 
avec et

[Link] Tenseur d’ordre supérieur


On peut définir un vecteur V par ses composantes Vi , ou par les coefficients de la forme linéaire
  
X  X V  X iVi , car la base choisie est orthonormée (voir les notions de vecteurs covariants et
contravariants).
On peut alors considérer le vecteur comme un tenseur du premier ordre.
De même, une fonction scalaire peut être considérée comme un tenseur d’ordre zéro.

Un tenseur du troisième ordre S est un opérateur linéaire qui, à tout vecteur Z fait correspondre un tenseur
du second ordre T .

T  S (Z ) ou encore Tij  Sijk Z k

[Link] Produit tensoriel

   
On définit le produit tensoriel du vecteur U par le vecteur V , noté U V , comme le tenseur d’ordre deux,
     
  
défini par la forme bilinéaire qui aux vecteurs X et Y fait correspondre U  X V Y

 
Les 9 produits tensoriels ei  e j définissent une base de l’espace vectoriel des tenseurs d’ordre deux, si bien
que l’on peut écrire un tenseur T comme

 
T  Tijei  e j

ou encore, par exemple,

- 11 - Golay MMC
MMC

 
u1v1 u1v2 u1v3 

    
 
u  v  ui v jei  e j   u2v1 u2v2 u2v3 
 
u v u3v2 u3v3 
 3 1
 

[Link] Contraction et produit contracté

  
Soit le produit tensoriel A  B  C , on appelle contraction, l’opération qui lui fait correspondre le vecteur
  
A(B  C ) . Le produit contracté d’un tenseur d’ordre 4 R et d’un tenseur d’ordre 3 S est défini par le tenseur
d’ordre 5

           
  
R  S  Rijklei  e j  ek  el  S pqrep  eq  er  Rijkm Smqrei  e j  ek  eq  er

Le produit doublement contracté d’un tenseur d’ordre 4 R et d’un tenseur d’ordre 3 S est défini par le
tenseur d’ordre 3

         
  
R  S  Rijklei  e j  ek  el  S pqrep  eq  er  Rijnm Smnrei  e j  er

Par exemple, le produit doublement contracté de deux tenseurs d’ordre 2 T et T  est le scalaire
   
  
T  T   Tijei  e j  T  pqep  ea  TijTji

1.2 Changement de repère

1.2.1 Matrice de passage

     
Soit e1, e2 , e3 une base orthonormée et e1, e2, e3 une autre base orthonormée.

On définit la matrice de passage Q telle que :

   
e1  Q11e1  Q12e2  Q13e3
   
e2  Q21e1  Q22e2  Q23e3
   
e3  Q31e1  Q32e2  Q33e3

ou encore, en notations indicielles

 
ei  Qije j

et en notation matricielle
 
e   Q  e 
Les deux bases étant orthonormées, on doit avoir
   
ij  ei  e j  Qikek  Q jlel  QikQ jl kl  QikQ jk

ce qui montre que la matrice inverse de Q est QT . En particulier on tire la relation inverse :
 
ei  Q jie j

1.2.2 Vecteurs
      
Soit V un vecteur de composantes Vi dans la base e1, e2 , e3 et Vi  dans la base e1, e2, e3 .
  
V  Viei  Vie i
Golay MMC - 12 -
 Notations tensorielles

En utilisant la matrice de passage

  
V  Viei  VQ e
i ki k

soit

Vk  VQ
i ki
et Vk  VQ
i ik

ou encore, en notation matricielle

   
 
V   Q  V
   et V   Q  V 
T

Remarque : le produit scalaire est un invariant, c’est à dire que cette fonction est indépendante du repère
choisi.
En notation indicielle
   
V W   VkW k  VQ WjQkj  ijVW
i ki i j
 VW
i i
 V W

et en notation matricielle
     T 


  


 
V W   V  W    Q  V  Q  W
 
 
  
 
 
T
 
 
 V Q  Q  W  V W  V W
     
1.2.3 Application linéaire
     
Soit A une application linéaire, de composantes Aij dans la base e1, e2 , e3 . et Aij dans la base e1, e2, e3 .

En notation indicielle

Wi   AikVk  QijWj  Qij AjmVm  Qij AjmQkmVk

d’où

Aik  Qij AjmQkm

et en notation matricielle
    
      
   
    
T
W   A V   Q  W  Q  A V  Q  A Q  V
 
soit
T
A  Q  A Q 
       

1.2.4 Forme bilinéaire

     
Soit A une application linéaire, de composantes Aij dans la base e1, e2 , e3 . et Aij dans la base e1, e2, e3 .
 
( W )  AijVW
AV i j
 AijVW
i
 j  AijQkiVkQ
 mjWm

soit
  AijQkiQmj
Akm

et en notation matricielle
     
 
A(V W )  V A W  V  A W  
     
  T      T   
     
T
T
 Q  V   A Q  W   V  Q  A Q  W 
           

- 13 - Golay MMC
MMC

soit
T
A  Q  A Q 
       

1.2.5 Tenseur d’ordre 2

Soit T un tenseur d’ordre 2, en notation indicielle
       
T  Tijei  e j  Tijei  e j  TijQkiek  Qmjem  TijQkiQmjek  em

puis

Tkm  TijQkiQmj

2 Permutations et déterminants

2.1 Les symboles de permutation

On introduit les symboles de permutation

1 si i j k est une permutation paire de 1 2 3


ijk 
 1 si i j k est une permutation impaire de 1 2 3



0 si deux indices sont répétés

Ces symboles représentent le produit mixte des vecteurs de base
  

ijk  ei , e j , ek 
ijk sont les composantes d’un tenseur du troisième ordre, qui représente, par exemple, la forme trilinéaire
produit mixte:
  
 
U V W  ijk U iVjWk

Avec un peu de patience on peut démontrer les résultats suivants

  
 il im in 

 
  
ijk lmn  Det  jl jm jn 
 

 
 km kn 

  kl
 
       
 ijk imn jm kn jn km
 ijk ijn  2km

 ijk ijk  6

2.2 Déterminant d’une matrice

Les symboles de permutation permettent le calcul du déterminant d’une matrice par

ijk Det(A)  mnp Aim Ajn Akp

(1.9)
ou encore

1
Det(A)    A A A
6 ijk mnp im jn kp
On peut également déterminer l’inverse d’une matrice

Golay MMC - 14 -
 Notations tensorielles

1
B  A1 et B ji    A A
2Det(A) imn jpq mp nq

2.3 Polynôme caractéristique

Les valeurs propres d’un tenseur du second ordre sont obtenues par la résolution de l’équation caractéristique
P ()  Det(A  I )

soit en développant

1
  (A  im )(Ajn  jn )(Akp  kp )  0
6 ijk mnp im
ou encore

P ()  I 3  I 2   2I 1   3

avec
 1
 I 3  ijk mnp Aim Ajn Akp  Det(A)
 6

I  A A  A A   1 (Tr A)2  Tr A2 
1 
ij ji  
 2 
2  ii jj 2  

 I1 Aii Tr A



I 1, I 2 , I 3 sont appelés les invariants fondamentaux du tenseur A.

2.4 Adjoint d’un tenseur antisymétrique

Soit  un tenseur antisymétrique
 12 31 
 0
  12 0 23 
 
 31 23 0 

on peut également lui associer le vecteur

 
 
 

 1  23 

 
 
 


   2   31 

   

      

 

   
 
 3  
 12  

 
 
 


soit
 0 3 2 


  3 0 1 
 
 2 1 0 


Le vecteur  est le vecteur adjoint du tenseur antisymétrique  . En notation indicielle on a :






ij  ijk k




i  21 ijk jk


 (1.10)

- 15 - Golay MMC
MMC

3 Calcul vectoriel et analyse vectorielle

3.1 Calcul vectoriel

Le produit vectoriel
  
c  a b
s’écrit en notation indicielle
 
ciei  ijk a jbkei

On peut montrer que

        
     
   
a b  c  a c b  b c a
       
         
a b  c  d  a c b d  a d b c

3.2 Analyse vectorielle

On note d’une virgule la dérivée partielle, soit  i   . Les opérateurs exposés dans cette partie seront
x i
exprimés dans un repère cartésien orthonormé.

* Soit f une fonction scalaire

Le gradient d’une fonction scalaire est un vecteur

 f 
 

 


  x 

 

1

   f 
grad f  f  fi ei  
 


  x 
2

 

  f 

 


 3
 x 

Le laplacien d’une fonction scalaire est un scalaire

2 f 2 f 2 f
 f  fii   
x 12 x 22 x 32

* Soit v un vecteur
La divergence d’un vecteur est un scalaire
 v v v
Div v  vii  1  2  3
x 1 x 2 x 3

Le rotationnel d’un vecteur est un vecteur

 v 3 v2 

  
 x 2 x 3 
    v v 

rot v    v  ijk vk  j ei   1  3 
 x x 1 
 3 
 2  v  v1 
  

 x 1 x 2 
Le gradient d’un vecteur est une matrice

Golay MMC - 16 -
 Notations tensorielles

 v v1 v1 
 1
 x x 2 x 3 
 1
    v v2 v2 
 v  vi j ei  e j   2 
 x 1 x 2 x 3 
 v v 3 v 3 
 3
 
 x 1 x 2 x 3 

Le laplacien d’un vecteur est un vecteur

  2v
 2v1 2v1  

 1
  


 2 2 2 
  
  x  x  x 3  v 
 

1 2

  1
   2
v  2
v  2
v 2 

v  
 v  vi jj ei   2  2 2
 
  2

 x 1 x 2 2 2 
x 3     


  2
v  2
v  2
v3 

  v 3 
 


3
 3
 

 2
x 22 x 32 
 x 1

 



* Soit T un tenseur du second ordre

La divergence d’un tenseur est un vecteur
 T 
 11  T12  T13 
 x x 2 x 3 
 1
  T T22 T23 
Div T  Tij  j ei   21   
 x x 2 x 3 
 1

 T31 T32 T33 
   
 x 1 x 2 x 3 

* Quelques formules utiles

   
 
Div f a  f Div a  a  grad f
       
 
Div a  b  b  rot a  a  rot b
 
 
Div rot a  0
 


 
rot grad f  0
 
 
grad f g  f grad g  g grad f
     
 
rot f a  f rot a  grad f  a


  

Div grad f   f

 
  
rot rot a  grad Div a  a 

3.3 Transformation d’intégrales


Soit  un domaine borné et  sa frontière, de normale n .
Soit  une fonction scalaire, alors
 
  n dS   grad  dV

Soit A un vecteur, alors

- 17 - Golay MMC
MMC

  
 A  n dS   Div(A) dV

Soit T un tenseur, alors


T  n dS   Div(T ) dV
  
Soit  un domaine plan de normale n , de frontière  . Soit U un vecteur défini sur ce domaine. Si  est le
vecteur unitaire tangent à  , alors
    
 rot(U )  n dS  U   dl

Tous ces résultats sont issus du théorème de la divergence

 t jkl nl dS   t jkl l dV

4 Formules essentielles en Mécanique des Milieux Continus

4.1 Coordonnées cartésiennes orthonormées

   
OM  xex  yey  zez
   
* Soit v  vxex  vyey  vzez un vecteur, alors

 v vx vx 
 x
 x y z 

 vi      v vy vy 
ei  e j  vi j ei  e j  
y
(v )  v 
x j   x y z 
 v v z vz 
 z
 x y z 
 
et
 v    vy
x i
 v
divv  i  vii  Tr grad (v )  v  I  x 
x y
 v
 z
z

 
 
v  div (v ) 
 2 vi 
x j x j
   
ei  vi jj ei  vxex  vyey  vzez

* Soit f une fonction scalaire, alors

 f 
 
 
 x 
f    f 

grad ( f )  f  ei  fi ei   y 
x i 
 
f 

 
 z 
 

et


f  div grad ( f )   2 f
x j x j
 f jj 
2 f

2 f

x 2 y 2 z 2
2 f

 
T Txy Txz 

 xx
   
* Soit T  Tij ei  e j  Tyx Tyy Tyz  un tenseur symétrique du deuxième ordre, alors:

 
 
T Tzy Tzz 
 zx
 

Golay MMC - 18 -
 Notations tensorielles

 T Txy Txz 

 xx   
 x y z 
Tij   T Tyy Tyz 
  yx
div(T )  e  Tij  j ei     
x j i  x y z 
 
 Tzx Tzy Tzz 
 x  y  z 
 
et
 Txz 
 2Tij
    Txx Txy
T  e  e j  Tij kk ei  e j  Tyx Tyy Tyz 
x k x k i  T Tzy Tzz 
 zx 

4.2 Coordonnées cylindriques

  
   OM  1 OM  OM 
OM  rer  zez et  er   e   ez
r r  z
   
d (OM )  er dr  rd e  ez dz
  
er e ez
0  0  0
r 
r r
er  e  ez
 e   er  0
  
er e ez
0  0  0
z z z
   
* Soit v  vrer  ve  vzez un vecteur, alors

 v   
 r 1  vr  v  vr 
 
 r r     z 
     
 v 1  v  v  v 
grad (v )  v    r
 r r   
 z 

 v v z 
 z 1 v z 
 r r  z 
 
et
  
  v
r
v
div v  Tr (v )  v  I  r  r 
r
1 v   v z
r 

z
   2 v vr    2 vr v  
 
v  div v  vr  2
 r 
 2 er  v  2
r   r  r 

 2 e  vzez

* Soit f une fonction scalaire, alors

 f  1 f  f 
grad ( f )  f  er  e  e
r r  z z
et

2 f 1 f 1 2 f 2 f
f  div f     2 
r 2 r r r  2 z 2

- 19 - Golay MMC
MMC

 
T

 rr

Tr  Trz 
 
* Soit T  T

 r T Tz  un tenseur symétrique du deuxième ordre, alors:

 
T

 zr Tz  Tzz 
 


 Trr 1 Tr  Trz Trr  T  

    


 r r  z r 


  T  T  T 2T 

 1 
div(T )   r
 
 z
 r


 r r  z r 


  T 1  T  T T 



zr
 z 
 zz
 zr 


  r r    z r 

 

4.3 Coordonnées sphériques

  
  OM  1 OM  1 OM 
OM  rer et  er   e   e
r r  rsin  
   
d (OM )  er dr  rd e  rsin  d  e
  
er e e
0  0  0
r 
r r
er  e  e
 e   er  0
   
er  e  e  
 sin e   cos e   sin er  cos e
  
   
Soit v  vrer  ve  ve un vecteur, alors
 
     

 vr 1  vr 
 1  1 vr
 




   v    v 
 


r 
 r  


r  sin  

  

  
 
   
1  v
   
1  1 v
  
v    
grad (v )  v  


  v 
r 

  cotg v  





r 

r  


 r  sin  

 

  
   
1 v 1  1 v
 
 v  
 





 cotg v  v 
r  

 r r  r  sin   
 


et
  v v 1 v 1 v v
divv  v  I  r  2 r    cot g  
r r r  r sin   r

  

v  2 v  1
(sin v ) 1 v 
 
 r r 2  r sin 
  sin   

 
    v v   
 
v  div (v )  

 v   r 



2 

r 2  
v 

cos  

2 sin2  sin2    

 
 


 

  vr v v   
 v  2 
  

 cotg    
 r 2 sin     2 sin    

 


Golay MMC - 20 -
 Notations tensorielles

* Soit f une fonction scalaire, alors


 f 


 


  r 

 
 1  f 

grad ( f )  
 


 r   

 1 f  

 


 r sin   
 

et

 
2 f
f  div grad ( f )  2 
r
2 2 f

1
r  2 r 2
cotg 
f

1 2 f
 r 2 sin2  2
 
T
 rr

Tr  Tr  
 
* Soit T  T
 r T T  un tenseur symétrique du deuxième ordre, alors:
 
 
T
 r T T 
 


 


 Trr Tr  Tr  

1  1 
 1 2Trr  T  T  Tr  cot g  

 


 


 





r r  r sin   r 



 

 Tr  T  T 

1  1  
 1 (T  T )cotg   3Tr 

  

div(T )  


 






r r  r sin   r 




 


 Tr  T  T 








 r
1
r 

 1
r sin  


r
 
 1 2Tcotg   3Tr 










 



4.4 Comment retrouver les formules

Nous nous plaçons par exemple en coordonnées cylindriques. On note
    
v  vrer  ve  vzez  viei avec i  r   z et  i    1   
r r  z
Donc, avec cette convention
 
 e  er
er ,  et e,  
r r
Chercher le gradient d’un tenseur consiste à augmenter l’ordre de ce tenseur, soit

()  () j  e j

Si on applique cette remarque à un vecteur, on obtient :

  
(v )  (viei ) j  e j

En n’oubliant pas de dériver les vecteurs de base, car nous sommes dans un système de coordonnées
cylindrique,
        
v  vi j ei  e j  vi ei, j  e j  vi j ei  e j  vi ei,  e
     
 vi j ei  e j  vr er ,  e  v e,  e
  v   v  
 vi j ei  e j  r e  e   er  e
r r
Pour obtenir l’opérateur divergence, il suffit de contracter doublement avec le tenseur unité d’ordre 2,

div()  ()  1
soit dans le cas d’un vecteur :

- 21 - Golay MMC
MMC

  v v v 1 v vz
div(v )  (v )  1  vii  r  r  r  
r r r r  z
et donc l’opérateur Laplacien pour un scalaire

r  2 1  1  2  2
  div   ii     2  2
r r 2 r r r  2 z
Appliquons maintenant cette méthodologie à un tenseur d’ordre 2.
  
(T )  T e
ij i 
 e j  ek
k
        
 Tij k ei  e j  ek  Tij ei,k  e j  ek  Tij ei  e j ,k  ek
        
 Tij k ei  e j  ek  Tij ei,  e j  e  Tij ei  e j ,  e
   Trj    T j   
 Tij k ei  e j  ek  e  e j  e  er  e j  e
r r
Tir    Ti   
 e  e  e  e  er  e
r i r i
Pour obtenir la trace de ce tenseur d’ordre 3 on contracte les deux derniers indices :

   Tr   T  Tir 
div T   (T )  1  Tij  j ei  e  e  e
  r  r r r i
 T 1 Tr  Trz T Trr  
  rr     e
 r r  z r r  r
 T 1 T Tz Tr  Tr  
  r     e
 r r  z r r  
 T 1 Tz  Tzz Tzr  
  zr    e
 r r  z r  z

On peut donc maintenant retrouver l’opérateur Laplacien d’un vecteur :

 
v  div v  
v v
vr     v   r
 vr    v   r  r e  vir e
 vi jjei  e  er  e  r
r r r r r i
 2 v vr    2 vr v   
 vr  2  2 er  v  2  2  e  vzez
 r  r   r  
r 

Golay MMC - 22 -
 Notations tensorielles

5 A retenir

Convention de sommation :
 
V  Viei

Produits tensoriels :
 
u1v1 u1v2 u1v3 

    
 
u  v  ui v jei  e j   u2v1 u2v2 u2v3 
 
u v u3v2 u3v3 
 3 1
 

Symboles de permutation :

1 si i j  k est une permutation paire de 1 2 3


  
ijk   1 si i j  k est une permutation impaire de 1 2 3
 ei , e j , ek  

0
 si deux indices sont répétés


Produit vectoriel :
   
c  a  b  ijk a jbkei

Quelques opérateurs :
        
Div v  vii , rot v    v  ijk vk  j ei ,  v  vi j ei  e j , Div T  Tij  j ei

En systèmes de coordonnées cylindrique ou sphérique, mieux vaut utiliser un formulaire !

- 23 - Golay MMC
MMC

Golay MMC - 24 -
 Cinématique

CINEMATIQUE

1 Le mouvement et ses représentations

1.1 Configuration
  
L’espace physique est rapporté à un repère orthonormé direct (O e1 e2  e3 ) . L’ensemble des particules ou
points matériels constituant le milieu continu étudié, occupe à chaque instant t, un ensemble de positions dans
l’espace: c’est la configuration du système à l’instant t, noté (t ) (d’intérieur (t ) et de frontière (t ) ).

On introduit aussi la notion de configuration de référence: c’est la configuration particulière du système à un

instant t0 fixé. Souvent on prendra 0  (0) , et on parlera alors de configuration initiale.

Toute particule M 0 de 0 est repérée par son vecteur position X (t ) dans la configuration de référence. Toute

particule M de (t ) est repérée par son vecteur position x (t ) dans la configuration actuelle (à l’instant t).

0 
 
 X t (t )

e3 0

e2 (t ) M
x
 M0  
e1 X  
u Xt

Figure 1 : Configurations de référence et actuelle

La position de chaque particule M sera donc déterminée si on connaît sa position dans la configuration de
référence et une fonction  telle que:
 
 
x (t )   X  t
(2.1)
  
 définit le mouvement par rapport à (O e1 e2  e3 ) . On devra donc déterminer trois fonctions scalaires, telles
que:

 x 1  1 (X1 X 2  X 3  t )


x   (X  X  X  t )
 2 (2.2)


2 1 2 3

x   (X  X  X t)
 3 3 1 2 3

Dire que le milieu est continu, c’est dire que  est une fonction continue et biunivoque de X . On supposera
que  est différentiable. Le déplacement par rapport à la configuration 0 , à l’instant t, de la particule M 0 est
le vecteur
  
u (X  t )  x (X  t )  X
(2.3)

- 25 - Golay MMC
MMC

1.2 Variables de Lagrange et variables d’Euler

Une grandeur attachée à une particule (masse volumique, vitesse, ...) peut être définie,

- Soit en fonction de X et t : variables de Lagrange

- Soit en fonction de x et t : variables d’Euler
Le vecteur vitesse d’une particule M est défini par
 
  dOM (X  t )
V (X  t )  
dt t (2.4)

Le vecteur accélération d’une particule M est défini par

  
  dV (X  t ) 2(X  t )
(X  t )  
dt t 2 (2.5)

1.2.1 Trajectoire
On appelle trajectoire d’une particule, la courbe géométrique lieu des positions occupées par cette particule au
 
 
cours du temps. x (t )   X  t est une représentation paramétrée en temps de la trajectoire. Par définition
de la vitesse,

  dOM dx  dx  dx 
V (x  t )   1 e1  2 e2  3 e3
dt dt dt dt
les trajectoires peuvent être obtenues par la résolution des trois équations

dx 1 dx 2 dx 3
   dt
V1 (x 1 x 2  x 3  t ) V2 (x 1 x 2  x 3  t ) V3 (x 1 x 2  x 3  t ) (2.6)

1.2.2 Lignes de courant

A un instant donné, on appelle lignes de courant du mouvement, les lignes qui sont en tout point tangentes au
vecteur vitesse de la particule située en ce point. Soit pour t fixé, deux équations :

dx 1 dx 2 dx 3
 
V1 (x 1 x 2  x 3  t ) V2 (x 1 x 2  x 3  t ) V3 (x 1 x 2  x 3  t ) (2.7)
   
Remarque: Pour un mouvement stationnaire (ou permanent) V (x  t )  V (x ) . Les lignes de courant et les
trajectoires sont confondues.

1.3 Dérivées particulaires

1.3.1 Définition
Lorsque l’on suit une particule dans son mouvement, la grandeur A attachée à la particule ne dépend que de
t. Par définition, on appelle dérivée particulaire de A à l’instant t , la dérivée de A par rapport à la seule
variable t .

En variables de Lagrange: A  A(X  t )

dA  A 
(X  t )  (X  t )
dt t (2.8)

En variables d’Euler: A  A(x  t )

Golay MMC - 26 -
 Cinématique

 A  A 
dA(x  t )  (x  t )dt  (x  t )dx j
t x j
dA  A  A  dx j
(x  t )  (x  t )  (x  t )
dt t x j dt
dA  A  A 
(x  t )  (x  t )  (x  t )Vj
dt t x j

ou encore

dA A
dt

t
 V  A   (2.9)

1.3.2 Application à l’accélération

  
   
(x  t ) 
dV (x  t ) V
dt

t
 V  V   (2.10)

que l’on peut également écrire


  V 1  2   
(x  t )   V  rotV V
t 2

2 Déformation d’un milieu continu

2.1 Notion de déformation

On dira qu’un milieu continu en mouvement subit des déformations si les distances relatives des points
matériels varient au cours du temps.
En différenciant (2.1), on obtient:

   i 
dx (t )    dX dx iei  dX jei
X j

On note F l’application linéaire qui fait passer de l’espace vectoriel dans lequel peut varier dX dans l’espace

vectoriel où varie a priori dx . Cette application linéaire, appelée tenseur gradient ou application linéaire
tangente, permet donc le passage de la configuration 0 à la configuration (t ) .

0
F (t )
 0
e3
 
e2 dX (t ) 
dx
 M
e1 M0

Figure 2 : Application linéaire tangente

En notation indicielle,
 x x 1 x 1 
 1
 X X 2 X 3 
 1
i x i  x x 2 x 2 
Fij   soit F   2 
X j X j  X1 X 2 X 3  (2.11)
 x x 3 x 3 
 3
 
 X1 X 2 X 3 

- 27 - Golay MMC
MMC

2.2 Tenseur des déformations

2.2.1 Définition
Le tenseur gradient décrit la transformation locale au voisinage d’une particule donnée. Afin de rendre compte
des déformations, c’est à dire des changements de forme autour de cette particule, on s’intéresse à l’évolution
du produit scalaire de deux vecteurs matériels pris respectivement dans les deux configurations 0 et (t ) .
    
Considérons trois particules voisines X , X  dX , X  dX  . Après déformations, elles occupent dans (t )
    
les positions respectives x , x  dx , x  dx  .

0
(t )
 0
e3  
 dx  (t )
e2 dX 

 dX M 
e1 dx
M0

Figure 3 : Notion de déformation

   x   

        
dx  dx   F (X  t )dX   F (X  t )dX    k dX i    x k dX j 
     X i   
  X j


d’où sa variation autour de la transformation

 
    
 x x k 
  
dx  dx   dX  dX    k  ij  dX idX j  Fki Fkj  ij  dX idX j
 X X  
 i j 

soit
     
dx  dx   dX  dX   2 dX  dX 

en posant

1  T 
 F (X  t ) F (X  t )  1
2   (2.12)

L’application linéaire  est appelée tenseur des déformations. Cette application est symétrique mais dépend
  
bien sûr de la base (O e1 e2  e3 ) initialement choisie.

2.2.2 Remarques
* S’il n’y a pas de déformations, alors   0 (et inversement).
T
* C  F F est appelé le tenseur des dilatations. Ce tenseur est symétrique.
On peut démontrer :

Théorème 1: Les valeurs propres de C sont strictement positives.

 
Théorème 2: Det F   0 t
 

Théorème 3:  est symétrique et possède les mêmes vecteurs propres que C .

* Variation de longueur
   
Soit dX   dX  dl 0 ex et dx  dl , alors

Golay MMC - 28 -
 Cinématique

     
dx  dx   dX  dX   dl 2  dl 02  2 dX  dX   2dl 02 xx

ou encore, si les déformations sont petites

dl dl  dl 0
 1  2xx  1  xx  xx 
dl 0 dl 0

xx représente au premier ordre la variation de longueur dans la direction x .

* Variation d’angle
   
Soit dX  dl 0 ex , dX  dl 0 ey , alors

     
dx  dx  dX  dX  cos dldl   2 dX  dX   2dl 02 xy

ou encore,

2xy  cos  1  2xx 1  2yy

donc xy représente au premier ordre la variation d’angle entre les directions x et y .

2.2.3 Autre écriture

D’après (2.3) et (2.1)

x u
F (X  t )  (X  t )  1  (X  t )
X X
soit
 T T 
1  u u u u 
   (X  t )  (X  t )  (X  t ) (X  t )
2 X X 
 (2.13)
 X X

ou encore en notation indicielle

1  ui u j uk uk 

ij     
2  X j X i X i X j 

2.2.4 Cas des petites perturbations

 u
Cette hypothèse correspond au cas où u (X  t ) et (X  t ) sont petits.
X

En reprenant (20) et en ne retenant que les termes d’ordre 1, on obtient:

 T 
1  u 
 HPP   (X  t )  u (X  t )
2  X  (2.14)
 X 

ou encore en notation indicielle

1  ui u j 

ijHPP   
2  X j X i 

- 29 - Golay MMC
MMC

2.3 Conditions de compatibilité


A tout déplacement u on fait correspondre une déformation  . On peut aussi se poser le problème inverse.
Ce problème est dit ’problème de compatibilité géométrique d’un champ de déformation’, ou encore ’problème
d’intégrabilité d’un champ de déformation’.
Les conditions de compatibilité peuvent être établies dans le cas général, cependant nous ne les établirons que
dans le cas des petites perturbations.

Décomposons maintenant le gradient des déplacements en une partie symétrique  et une partie
antisymétrique  .

u
(X  t )  (X  t )  (X  t )
X
  1  ui u j 
T
1  u 
  (X  t )  u (X  t ) ij   
2  X X

 2  X j X i 

On a

ij k  ki j  jk i

soit en dérivant une nouvelle fois

ij kl  ij lk i j k l dans { 1 2 3}

 i j k l ij kl  kl ij  ik  jl  jl ik  0

(2.15)
ou encore

 22323  3322  2233  permutation circulaire
Six équations 


  3231  1233  3321  permutation circulaire

 1323

Réciproquement, si  vérifie (2.15), alors les formes différentielles

 
d ij   ki j  jk i  dx k

sont exactes; elles permettent donc de construire le champ  de tenseur antisymétrique. On vérifie ensuite
que les formes différentielles
 
dui  ik  ik  dx k

sont exactes, d’où la possibilité de construire un champ de déplacement u (X  t ) défini dans 0 .

3 Transport, dérivées particulaires

3.1 Transport d’un volume

  
Soit d0 un élément de volume de la configuration de référence, défini par trois vecteurs dX1 dX 2  dX 3 . Par la
  
transformation, ces trois vecteurs se transportent en trois vecteurs dx 1 dx 2  dx 3 qui définissent dans la
configuration actuelle un volume d .

Golay MMC - 30 -
 Cinématique

 d
 dx 3
dX 3 d0

 dx 2
dX 2

dx 1

dX1

Figure 4 : Transport d’un élément de volume

  
Le volume d est représenté par le produit mixte des vecteurs dx 1 dx 2  dx 3 :
  
d   dx 1  dx 2   dx 3

donc

d   ijk dx 1 j dx 2k dx 3i

Or, d’après (2.11)

d   ijk Fjp Fkq Fir dX1p dX 2q dX 3r

et, d’après (1.9)

  

d   pqr det(F ) dX1p dX 2q dX 3r  det(F ) dX1  dX 2  dX 3 
donc en définitive

d   Det(F ) d 0
(2.16)

3.2 Transport d’une surface orientée


Soit dS un élément de surface de la configuration de référence de normale N . Par la transformation, cette

surface se transporte en une surface ds de normale n dans la configuration actuelle. En considérant un
 
vecteur V dans la configuration de référence qui se transporte en un vecteur v dans la configuration actuelle,
   
on peut définir l’élément de volume (dS N ) V qui se transporte en un élément de volume (ds n )  v .


N 
n
dS ds

Figure 5 : Transport d’un élément de surface

D’après (2.16)
   
ds n  v  det(F ) dS N V

et comme avec (2.11)

 
v FV

    T   
ds n  FV   ds F n  V  det F dS N V

   

- 31 - Golay MMC
MMC

on obtient finalement

 T 
ds n  det(F )F dS N
(2.17)

3.3 Dérivée particulaire d’une intégrale de volume


Soit K (t )  (t ) k (x  t ) d  , une intégrale de volume sur le domaine (t ) dans la configuration de référence.
Pour en déterminer la dérivée temporelle, nous devons au préalable exprimer K (t ) sur la configuration de
référence pour "passer" la dérivation sous l’intégrale. En effectuant le changement de variable (2.1), et en
utilisant (2.16)

d   Det(F ) d 0  J d 0

on obtient

K (t )   k ((X  t ) t ) J d 0
0

puis

dK dk dJ 
   J  k  d 
dt 0 
 dt dt  0

A ce stade nous devons expliciter dJ / dt . En utilisant les notations indicielles, et en particulier les symboles
de permutation, on a:

1
J  det F    F F F
6 ijk pqr ip jq kr
soit

dJ 1 Fip
 ijk pqr F F
dt 2 t jq kr
or

  
Fip 2i (X  t )   i    V (X  t )   v (x t )  vi xl  vi F
t

t X p

X p  

 t  X i  
X p i
 x X x lp
p l p l

donc

dJ 1 v
 ijk pqr i Flp Fjq Fkr
dt 2 xl

mais

pqr Flp Fjq Fkr  ljk det F

soit

dJ 1 v v v
 ijk ljk i det F  il i det F  i J
dt 2 x l xl x i

dJ 
 J div v
dt (2.18)

En reportant dans l’expression de dK / dt

Golay MMC - 32 -
 Cinématique

dK dk 
   J  k J divv  d 0
dt 0 
 dt 

puis en exprimant l’intégrale sur la configuration actuelle, on obtient finalement

dK dk 
 (t )   k divv  d 
dt  dt  (2.19)

En utilisant les égalités suivantes,

dk k 
  v  k
dt t
  
div(kv )  v  k  kdivv

on peut écrire (2.19) sous la forme

dK  k 
 (t )   div kv  d 
dt  t 

ou encore, en utilisant le théorème de la divergence

dK k  
 (t ) d   (t ) kv  n d 
dt t

Application fondamentale: conservation de la masse

La masse d’un système matériel qu’on suit dans son mouvement reste constante.

M  (t ) (x  t ) d  dM
0
et dt

où  est la masse volumique. On a alors:

d   
  divv  0  div v   0
dt ou t (2.20)

3.4 Dérivée particulaire d’une intégrale de surface

  
Soit K (t )  (t ) k (x  t )  n d  , une intégrale de surface sur le domaine (t ) dans la configuration de
référence. Pour en déterminer la dérivée temporelle, nous devons au préalable exprimer K (t ) sur la
configuration de référence pour "passer" la dérivation sous l’intégrale. En effectuant le changement de variable
(2.1), et en utilisant (2.17)

 T 
d  n  det(F )F d 0N

on obtient
  T 
 
K (t )   k (X  t ) t  J F d 0 N
0

puis

 T   d  T   
dK  dk
    J F N  k  J F  N d 
dt     0
0
 dt dt   
T
on doit donc calculer dF / dt

- 33 - Golay MMC
MMC

1 1
1 1 1
dF dF dF dF 1
F F I  F F 0   F F
dt dt dt dt

dF 1 vi X k   vi   
F  e  ej  e  e j  v
dt X k  x j i x j i

donc
T T
dF  1    T
 F v   T v F
dt  

et

  
dK  dk T 
 T    T  
    J F N  k  J divv F N  k  J T v F N d 0
0 
dt  dt 


dK  dk     T 
    divv k  k  T v J F Nd 0
0  
dt  dt

puis en exprimant l’intégrale sur la configuration actuelle, on obtient finalement


 
dK  dk     
 (t )   divv k  v k   nd 
 dt
dt  
(2.21)
en utilisant la dérivée particulaire, (2.21) s’écrit

 
dK  k        
 (t )   divv k   k v  v k   nd 
dt  t 


      
dK
dt
 (t ) 
k
 t

  
 rot k  v  v div k   nd 


car
         
 
rot k  v  k divv  vdiv k   k v   v k

Golay MMC - 34 -
 Cinématique

4 A retenir


On appelle Variables de Lagrange le temps et la position initiale : X et t

On appelle Variables d’Euler le temps et la position courante : x et t

Dérivée particulaire

dA A
dt

t
 V  A 
Application linéaire tangente

x i  
F e  ej
X j i

Tenseur des déformations

1  T 
 F (X  t ) F (X  t )  1
2  

Tenseur des déformations sous l’hypothèse des petites perturbations

1 T 

2

 u  u 
Transport d’un volume

d   Det(F ) d 0

Transport d’une surface

 T 
ds n  det(F )F dS N

Dérivée d’une intégrale de volume

dK  k 
 (t )   div kv  d 
dt  t 

Dérivée d’une intégrale de surface


   
dK  dk  
 (t )   divv k  v k   nd 
dt 
 dt 

- 35 - Golay MMC
MMC

Golay MMC - 36 -
 Equilibre

EFFORTS DANS LES MILIEUX CONTINUS

1 Définitions

1.1 Forces
Elles résument les effets mécaniques, autres que cinématiques, exercés sur le milieu continu considéré par le
reste du domaine physique. Leur schématisation à chaque instant repose sur la définition d’un champ de
   
vecteur (x  t ) et d’une mesure positive  , définis sur la configuration actuelle (t ) . (x  t ) est une densité
de force pour la mesure  .
 
* Si  est une mesure de volume, alors (x  t ) est une force volumique (densité volumique de force) définie
dans  (t) de la configuration actuelle, par la fonction
   
f  x  (t )  f (x  t )   3
 
* Si  est une mesure de surface, alors (x  t ) est une force surfacique (densité surfacique de force) définie
sur F (t ) de la configuration actuelle, par la fonction
   
F x  F (t )  F (x  t )   3

* ... etc ...

Remarques :
* Les forces sont définies sur la configuration actuelle.

* A un instant donné et en un point donné x de (t ) , on ne peut imposer à la fois le déplacement et la force!.
Mais l’un des deux doit être imposé. On note F (t ) la frontière où la force est imposée, et U (t ) la
frontière où le déplacement est imposé. Dans le cas des appuis mobiles, les composantes non imposées
cinématiquement le sont pour les forces

* Le monde extérieur au milieu considéré doit, pour imposer le déplacement U (t ) au bord U (t ) , exercer
 
des forces que nous noterons R(x  t ) . Comme elles sont à priori inconnues, nous les appellerons réactions
pour éviter de les confondre avec les autres forces qui, elles, sont données.

1.2 Vecteur-contrainte et tenseur des contraintes

1.2.1 Contrainte de Cauchy


dF
C Soit un corps (C) en équilibre par application d’un système d’actions mécaniques

2 extérieures. Imaginons qu’une surface  divise (C) en deux parties (1) et (2). La
 n
partie(1) est en équilibre sous les actions mécaniques extérieures qui lui sont
M
d appliquées et les actions mécaniques exercées par la partie (2). Nous admettrons que
1   
sur chaque élément de surface d de  , (2) exerce sur (1) une force dF (x  t n )12 de
  
densité superficielle T (x  t n ) .

- 37 - Golay MMC
MMC

     
dF (x  t n )12  T (x  t n ) d 
(3.1)
   
T (x  t n ) est le vecteur contrainte au point x , relativement à la facette d définie par son vecteur normal n .

La densité surfacique de forces exercées en x dépend de x, t et aussi de l’orientation de la surface  au


voisinage de x. Elle est linéairement dépendante de n . On introduit alors l’application  telle que:
    
T (x  t n )  (x  t ) n
(3.2)

L’application (x  t ) s’appelle le tenseur des contraintes de Cauchy en x à l’instant t ; il caractérise, dans la
configuration actuelle, les efforts intérieurs de cohésion exercés sur une partie du solide à travers l’élément de

surface n d

1.2.2 Autre écriture du tenseur des contraintes

En utilisant (2.17), (3.1) devient:
      
 
dF x (X  t ) t n(N  t )   N (X ) dS

où  est le tenseur
      
(X  t )  N  R 3  (X  t N )  (X  t )N  3

défini par
 T
(X  t )  (det F )  F
(3.3)

Cette application linéaire (X  t ) , définie pour X  0 , s’appelle le premier tenseur des contraintes de Piola-
Kirchoff en X à l’instant t; la composante ij est la i ème composante du vecteur contrainte exercée sur la

déformée d’une surface unité, normale à e j , de la configuration de référence. On prendra garde au fait que le

tenseur  n’est pas symétrique.

Si maintenant on cherche le vecteur "force de cohésion" dans la configuration de référence
   1       
   
dF 0 X  t N  F (X  t ) dF x (X  t ) t n(N  t )  S N (X ) dS

où S est le tenseur défini par

1
S F 
(3.4)

Cette application linéaire S (X  t ) , définie pour X  0 , s’appelle le second tenseur des contraintes de Piola-
Kirchoff en X à l’instant t. Attention, sa composante Sij n’est pas la i ème composante du vecteur contrainte

exercée sur la déformée d’une surface unité, normale à e j , de la configuration de référence, mais seulement la
i eme composante de son transporté dans la configuration de référence.
Selon le jeu d’écriture adopté, on a donc trois descriptions des contraintes :
1 1
  T
  T
  det F   F  det F  F S F
    (3.5)

Golay MMC - 38 -
 Equilibre

2 Equilibre

2.1 Le Principe des Puissances Virtuelles (Germain 1972)

Pour schématiser les efforts mis en jeu, il est commode d’imaginer des mouvements fictifs (ou virtuels) et
d’analyser le travail ou la puissance qui en résulte. Par exemple, pour évaluer les forces de gravité agissant sur
un objet, on peut imaginer de le soulever (mouvement virtuel de bas en haut).
Un milieu matériel étant isolé, on peut distinguer les actions extérieures qui agissent sur le milieu, des actions
intérieures qui représentent les liaisons existant entre toutes les parties du milieu.
Axiome d’objectivité
La puissance virtuelle des efforts intérieurs associée à tout mouvement rigidifiant est nulle.
Axiome d’équilibre
Pour tout milieu matériel repéré dans un référentiel absolu, à chaque instant et pour tout mouvement virtuel, la
puissance virtuelle des quantités d’accélération a est égale à la somme des puissances virtuelles des efforts
intérieurs i et des efforts extérieurs e .

2.2 Puissance virtuelle des efforts intérieurs


F


e3
 
e2 n
 (t )
e1 O 
f
(t )

Soit un milieu continu (t ) d’intérieur (t ) et de frontière  (t). Isolons maintenant un domaine (t ) de

frontière  (t) intérieur à  (t), et soit n la normale en un point de (t ) . A un instant t fixé, un mouvement

virtuel défini par une vitesse virtuelle v est appliqué à (t ) . Cette vitesse est supposée continue et
continûment dérivable sur (t ) .

Pour déterminer la puissance virtuelle des efforts intérieurs nous ferons les hypothèses suivantes:

* i admet une densité volumique p i 

i   pi dx

* i est en chaque point une forme linéaire des valeurs en ce point de dv et de ses dérivées premières:

En décomposant le gradient des vitesses virtuelles en une partie symétrique D et une partie antisymétrique
W ,

v
 D  W
x
 T

1  v  v  1  vi v j 
W   
Wij    
2  x x  2  x j x i 

- 39 - Golay MMC
MMC

 T
  v v j 
1  v 1 
D    v  Dij  i
 

2  x x  2  x j x i 


la densité volumique des efforts intérieurs devient:

pi  Ai Vi  Bij Wji  ij D ji

Le premier axiome du principe des puissances virtuelles impose que pour tout mouvement de solide rigide la
puissance des efforts intérieurs soit nulle. D’où:
 
- Soit un mouvement de translation: v  0 , W  0 et D  0

alors
 
i   pi dx   A  v dx  0  
dans
    
soit A  v  0  v , ou encore A  0
 
- Soit un mouvement de rotation: v  0 , W  0 et D  0

alors

i   pi dx   B  W dx  0  
dans

soit B  W  0  W , ou encore B  0 .

Donc en définitive:

i      D dx
(3.6)

On peut montrer que le tenseur  introduit ici correspond bien au tenseur des contraintes de Cauchy.

2.3 Puissance virtuelle des efforts extérieurs

Les efforts extérieurs comprennent
- des efforts exercés à distance par des systèmes extérieurs à  , supposés définis par une densité volumique

de forces f ,

- des efforts de cohésion schématisés par une densité surfacique de force T sur 
   
e   f  v dx  T  v dx
(3.7)

2.4 Application du Principe des Puissances Virtuelles


Si  est l’accélération et  la masse volumique de chacun des points de  , alors

      
d  v
a    v  v dx     div(v  v )  v dx
dt  t 

 
 v       
a    v  v  v  div(v )v   v dx
 t t 


et en utilisant la conservation de la masse (2.20) et la définition de l'accélération (2.10)

Golay MMC - 40 -
 Equilibre

 
a      v dx
(3.8)
En application du Principe des Puissances Virtuelles on obtient :
     
   D dx   f  v dx  T  v dx      v dx
(3.9)

Pour exploiter le fait que (3.9) est vérifié pour tout mouvement virtuel, nous allons faire apparaître v dans
chacun des termes.
En appliquant le théorème de la divergence, le premier terme devient :

v   
   D dx     dx      v  n dx   divx   v dx
x
Soit :
         
 T    n  v dx    f  divx      v dx  v
   

Ou encore
  


 f  div x
   dans 
   (3.10)

 T   n sur 




2.5 Equilibre
En considérant les développements du paragraphe précédent et en se ramenant au domaine (t ) , nous

pouvons donc écrire les équations d’équilibre d’un solide soumis à un champ de forces extérieures f dans (t )
 
, à un champ de forces extérieures F e sur F (t ) et à un déplacement imposé U i sur U (t ) .

Dans la configuration actuelle :

 
f (x  t )  divx (x  t )  0  x  (t )
(3.11)

 


F e (x  t )  x  F (t )
(x  t )  n(x  t )   

R(x  t )  x  U (t ) (3.12)


Dans la configuration de référence :
  
De même, si on note f 0 , R 0 et F 0 les densités volumiques et surfaciques de forces mesurées dans la
configuration de référence:
 
f 0(X  t )  divX (X  t )  0  x  0
(3.13)

 

F 0 (x  t )  x  x 1 (F (t ) t )
(X  t )  N (X  t )  
 

 (x  t )  x  0U
R 0
 (3.14)
Cas des petites perturbations

Reprenons (3.10), en l’exprimant en fonction de X

ij
fi x (X  t ) t   (x (X  t ) t )  0  x (X  t )  (t )
x j

- 41 - Golay MMC
MMC

ij X k
fi x (X  t ) t   (X  t ) (X  t )  0  X  0
X k x j

  
Or x (X  t )  X  u(X  t ) soit x (X  t )  1  u (X  t )
X X
On peut donc écrire l’équation d’équilibre sous la forme
1
ij  
 u 
fi x (X  t ) t   (X  t ) 1  (X  t )  0  X  0
X k  X 
  kj

Sous l’hypothèse des petites perturbations, on peut alors écrire:

1
 
 u  u
1  (X  t ) 1 (X  t )
 X  X
 

soit

ij  uk 
fi x (X  t ) t   (X  t ) jk  (X  t )  0  X  0
X k  X j 

Enfin, en ne retenant que les termes d’ordre 0, et après avoir effectué un développement de fi au voisinage de
X, on obtient:

ij
fi X  t   (X  t )  0  X  0
X j

soit
 
f (X  t )  divX (X  t )  0  x  0
(3.15)
  
Le raisonnement qui a permis de remplacer f x (X , t ) par f X , t  , permet aussi de remplacer F e (x (X , t )) par
  
F e (X , t ) et R(x (X , t )) par R(X , t ) . Donc, comme condition sur la frontière on obtient:

 
 

F e (X , t )
  X  0F
(X  t )  N (X  t )   


 R(X t )

 X   (3.16)

 0U



2.6 Autre présentation: Principe fondamental de la dynamique

(3.10) revient à écrire le Principe Fondamental de la dynamique. Dans un repère galiléen, pour tout système 
, le torseur dynamique (dérivée par rapport au temps du torseur cinématique) est égal à la somme des torseurs
des actions intérieures. Soit:

d  d 
 v dm   v d 
dt  dt  
   
dv 
  v divv d 
  dt 


  
 dv  d  
   v  v divv d 
 
  dt dt
   d   

    v   divv d 
 
 dt 



Golay MMC - 42 -
 Equilibre

donc avec la conservation de la masse

d  
 v dm   d 
dt  
 
  fd    nd 
 

et le théorème de la divergence
 
 d    fd    div d 
  

on retrouve le bilan de la quantité de mouvement

 
div   f  
L’équation de bilan sur les moments du principe fondamental de la dynamique s’écrit:

 dv    
 OM  dm   OM    n dx   OM  f dx (3.17)
dt

3 Quelques propriétés du tenseur des contraintes

Dans tous les développements à venir, nous nous placerons dans le cas des petites perturbations pour un
solide en équilibre. En conséquence, nous omettrons les variables x et t.

3.1 Symétrie du tenseur des contraintes

On sait que
    
 
OM    f dx  OM  n dx

soit en notation indicielle

 
 ijk x j (k  fk ) e i dx   ijk x j kl nl e i dx

puis, par application du théorème de la divergence

  
  ijk x j (k  fk )  (ijk x j kl ) e i dx  0
 xl 
  x (  f   )      dx  0
  ijk j k k kl l ijk kj  e i

et par application de l’équation du mouvement


 ijk kjei dx  0  (t )

c’est à dire

ijk kj  0 i

ce qui implique

123 23  132 32  0 213 13  23131  0 312 12  32121  0
23  32  0  13  31  0  12  21  0

donc en définitive

pq  qp

Le tenseur des contraintes est symétrique

- 43 - Golay MMC
MMC

3.2 Contrainte normale et contrainte tangentielle

  
T n  Considérons une facette de normale n . Tout naturellement, le vecteur contrainte
 
n T (n ) peut être décomposé en une composante normale n et une composante
 tangentielle  .
 n

    
n  T (n )  n  n    n
(3.18)
et

      2
2

    n   n    n 

    (3.19)

On dira que n est positive en traction et négative en compression.

3.3 Directions principales, contraintes principales

La matrice représentant le tenseur des contraintes est symétrique, elle est donc diagonalisable. Les valeurs
propres sont réelles et appelées contraintes principales I , II , III  . Les vecteurs propres, orthogonaux deux

à deux, sont les directions principales n I  n II  n III  . On a donc:

        
I  T (nI )  nI  II  T (nII )  nII  III  T (nIII )  nIII

3.4 Invariants
Le tenseur des contraintes possède trois invariants définis mathématiquement comme les coefficients de
 
l’équation caractéristique det    1 . C’est à dire les quantités scalaires :
 

I  Tr ()
(3.20)

1  
II  Tr ()2  Tr ( 2 )
2   (3.21)

III  Det()
(3.22)
Exprimées en fonction des contraintes principales, on obtient

I  I  II  III

II  I II  II III  III I

III  I II III

3.5 Cercles de Mohr

Connaissant le tenseur des contraintes  , on se propose de déterminer le domaine engendré par l’extrémité

du vecteur contrainte quand n varie. Par commodité, nous nous plaçons dans une base orthonormée dirigée
suivant les directions principales de  . Soit

Golay MMC - 44 -
 Equilibre


 
   
 


 n1 
 
 I 


0 0    n1 I 

 

     
n  n2    0 


II 0  et T   n2 II 
 

 
   
 




 3
n  0 
 0 III  


 3
n III 

 
   
 


avec n12  n22  n 32  1

on trouve aisément

n  I n12  II n22  III n 32

et

 2  n2  I2 n12  II2 n22  III

2
n 32

Dans l’hypothèse où les contraintes principales sont distinctes, on obtient alors après résolution du système:

 2  (n  II )(n  III )

n12 
(I  II )(I  III )

 2  (n  I )(n  III )

n22 
(II  I )(II  III )

2
 2  (n  I )(n  II )
n 3

(III  I )(III  II )

Si on ordonne les contraintes principales de telle sorte que I  II  III , alors

 2  (n  II )(n  III )  0

 2  (n  I )(n  III )  0

 2  (n  I )(n  II )  0

ou encore
2 2
 II  III   
  II  III 
  n 
2
   
 
 2  2 
(3.23)
2 2
 I  III   
  I  III 
  n 
2
  

 2   2 
(3.24)
2 2
 I  II   
  I  II 
  n 
2
   
 
 2  2 
(3.25)

 Dans le plan de Mohr, l’extrémité du vecteur contrainte, d’après (3.24), est donc
 intérieure au cercle centré sur On d’abscisse (I  III ) / 2 et de rayon
T (I  III ) / 2 . Par contre, d’après (3.23) (res. (3.25)), l’extrémité du vecteur
III I contrainte est extérieure au cercle centré sur On d’abscisses (II  III ) / 2
II
n (resp.( (I  II ) / 2 ) et de rayon (II  III ) / 2 (resp. (I  II ) / 2 ).

- 45 - Golay MMC
MMC

Description des Cercles principaux:

Nous allons étudier la description du grand Cercle de Mohr. Les facettes concernées
III  sont parallèles à la direction associée à la contrainte principale II .
 n
t    
On constitue avec les directions I,III,II un trièdre direct (O eI  eIII  eII ) , la normale n

de la facette évoluant dans le plan I III.
I
   
Et on définit l’angle   I , n  , et le vecteur t tel que (n, t , II ) soit direct.

On a alors
  
n  Cos  eI  Sin  eIII

et
  
T  I Cos  eI  III Sin  eIII
    
En utilisant les formules de changement de base de (O eI  eIII  eII ) à (n, t , II ) , on a donc

I  III I  III
n   Cos 2
2 2
I  III
  Sin 2
2


Lorsque la facette tourne autour de la direction de la contrainte
principale II d’un angle donné, l’extrémité du vecteur-contrainte tourne
I  III
sur le cercle de Mohr d’un angle double dans le sens opposé (autour du
III 2 I
centre du cercle).
2 n

T

Golay MMC - 46 -
 Equilibre

4 Exemples de tenseur des contraintes

4.1 Tenseur uniaxial

  0 0
   
   e1  e1   0 0 0
 
 0 0 0

L’équilibre des forces sur la frontière du domaine nous donne:

     
Sur 0 : n  e1 donc n  F 0 et F 0  e1

     
Sur 1 : n  e1 donc n  F 1 et F 1  e1

Sur la frontière latérale les pressions sont nulles.

On se trouve en présence d’un chargement uniaxial de traction/compression.
Si   0 c’est un état de tension uniaxiale
Si   0 c’est un état de compression uniaxiale

La direction principale est e1

4.2 Tenseur sphérique

p 0 0 

  pI   0 p 0 

 
 0 0 p 


Dans ce cas, toute direction est direction principale. La contrainte normale principale est -p. p est appelé la
pression. Si p  0 on a un état de compression, et si p  0 on a un état de tension.
Par exemple pour un fluide au repos:
D’après l’équation d’équilibre
 
div   g  0
 
divpI  g  0
  
gradp  g  0

p p p
soit   0 et  g
x 1 x 2 x 3

et p  p0  gx 3

Donc, pour un fluide au repos p  gx 3  Cste .

- 47 - Golay MMC
MMC

5 A retenir

Vecteur contrainte et tenseur des contraintes de Cauchy

    
T (x  t n )  (x  t ) n

Le tenseur des contraintes est symétrique !

Equilibre
 
f (x  t )  divx (x  t )  0  x  (t )

 

F e (x  t )  x  F (t )
(x  t )  n(x  t )  
 

R(x  t )  x  U (t )


Contrainte normale
    
n  T (n )  n  n    n

Contrainte tangentielle

      2
2

    n   n    n 

   

Golay MMC - 48 -
 Elasticité

ELASTICITE

1 Approche expérimentale: essai de traction

Pour déterminer l’évolution d’un système déformable, nous avons déjà déterminé les équations de la
cinématique et de la sthénique. A ces équations, il est maintenant nécessaire d’adjoindre une relation
supplémentaire reliant les efforts internes et les grandeurs cinématiques. Cette relation, appelée Loi de
Comportement, dépend du matériau considéré. La construction d’une loi de comportement est basée sur des
observations expérimentales.
Dans ce chapitre nous exposerons le modèle de comportement des matériaux élastiques, sous l’hypothèse des
petites perturbations.

Pour effectuer un essai de traction simple sur un métal, on utilise une

F éprouvette cylindrique caractérisée par :
- des extrémités surdimensionnées
- des congés de raccordement (pour éviter les concentrations de contrainte)
S0 L L L - une partie médiane cylindrique dans laquelle le champ de contrainte est
supposé homogène, de traction simple parallèlement à l’axe de l’éprouvette.
L’essai de traction consiste à enregistrer l’évolution de l’allongement relatif de la
longueur initiale L0 en fonction de la force de traction F , ou du rapport F / S 0
F
, où S 0 représente l’aire initiale de la section de l’éprouvette.

La figure ci-contre représente un tel enregistrement pour un acier

Plasticité inox. On remarque alors les propriétés suivantes :
irréversible - Le diagramme est indépendant de la vitesse de chargement
F Élasticité
11 
S réversible - La partie OA du diagramme est réversible. Si on charge jusqu’à un
σe A niveau inférieur à 0 , alors la décharge décrit la même courbe OA.

- La partie réversible est linéaire

- Si on effectue un chargement au-delà du seuil 0 , puis une
décharge, l’éprouvette présente une déformation permanente.
O
La partie réversible du diagramme de traction est, par définition,
L représentative du comportement élastique du matériau. 0 est la
11 
Déformation L
permanente limite initiale d’élasticité du matériau. La linéarité du segment OA A
caractérise le comportement élastique linéaire du matériau.

- 49 - Golay MMC
MMC

2 Loi de comportement élastique linéaire (en HPP)

2.1 Forme générale

A partir des observations expérimentales on peut écrire que les contraintes dépendent linéairement des
déformations. En l’absence d’effets thermique et de contraintes initiales on a :

(x  t )  C (x )  (x  t )
(4.1)

C est un tenseur du quatrième ordre, dont les composantes sont les coefficients d’élasticité du matériau.

ij (x  t )  C ijkl kl (x  t )

En utilisant les propriétés des tenseurs de contrainte et de déformation, on peut montrer que :

C ijkl  C jikl  C ijlk  C jilk

Le tenseur C , dont la matrice représentative comporte 81 composantes, ne dépend donc plus que de 21
paramètres indépendants.

2.2 Matériau élastique homogène isotrope

Toutes les directions sont équivalentes, de telle sorte que la loi de comportement est invariante dans toute
rotation de la configuration de référence. Ce modèle s’applique à la plupart des matériaux : acier, béton, ...
Si la configuration est libre de contraintes, alors la loi de comportement s’écrit :

   Tr () 1  2  
(4.2)
ou encore en notation indicielle

ij  kk ij  2ij

Les coefficients matériel  et  , qui dépendent de la particule considérée, sont appelés les coefficients de
Lamé. Leur expression en fonction du module d’Young E et du coefficient de Poisson  , est

E E
 et 
2 (1   ) (1   ) (1  2 )

ou

(3  2) 
E et 
 2(  )

avec, en inversant (4.2)

 1
  Tr () 1  
E E (4.3)

2.3 Matériau élastique homogène orthotrope

Le matériau possède trois directions privilégiées deux à deux orthogonales. La loi de comportement est
invariante par les symétries par rapport aux plans orthogonaux construits à partir de ces directions. Dans ces
matériaux, on peut classer les tôles laminées, les composites tissés, le bois, certains bétons structurés, ...
Dans ce cas on montre que la matrice de comportement est définie par 9 paramètres indépendants. Dans le
repère principal d’orthotropie, la loi se met sous la forme:

Golay MMC - 50 -
 Elasticité

 1 12 13 
 0 0 0 
 E E1 E1
 1 
  1 23 
   21 0 0 0   
 11    11 



  E2 E2 E2  
 
      32   22 


22   31 1  
 
 

  0 0 0   33 
 33 
  E3 E3 E3  


2

   


12  (4.4)
 0 1
0 
  
 12  
   0 0 0 



2 23 
  G12  

23 

   

2 
 0 1 13 
0 
  
 13 
 0 0 0  
 G23 
 
 0 1 
 0 0 0 0 
 G13 

Avec les conditions de symétrie

12 21 13  31  32 23

  
E1 E2 E1 E3 E3 E2

2.4 Matériau élastique homogène isotrope transverse

Un matériau homogène isotrope transverse est tel que la matrice de comportement est invariante par toute
rotation autour d’un axe privilégié. En utilisant cette invariance, on montre que seuls 5 paramètres
indépendants caractérisent le comportement. Si l’axe est porté par la direction 3, on a alors:
 1 12 13 
 0 0 0 
 E E1 E1
 1 
  1 13 
   21 0 0 0   
 11    11 
   E1 E1 E1   

 22     31  
 22 
   31 1   

33   0 0 0  
33 


  E3 E3 E3  


212 

  
12 

(4.5)
 0 1
0 
 
 
   0 0 0  

 223   G12  
 23 
     

213   0 1 
13 

   0 0 0 0  
 
 G13 
 
 0 1 
 0 0 0 0 
 G13 

2.5 Caractéristiques de quelques matériaux

Matériaux isotropes usuels:
Matériau E en Gpa   en kg/l

acier 210 0.285 7.8

fonte grise 90 à 120 0.29 7.1
aluminium 71 0.34 2.6
béton 10 0.15 2.4
fibre de verre E 73 0.15 2.54
Graphite HM 350 0.4 1.92
résine époxy 3.8 0.31 1.15

- 51 - Golay MMC
MMC

Matériaux composites:

Unidirectionnel Tissu Unidirectionnel Unidirectionnel

Verre/Epoxy Verre/Epoxy CarboneHT/Epoxy Kevlar/Epoxy
50% 50% 50% 50%

 en g/cm 3 1,87 1,87 1,49 1,32

E 1 en Mpa 38000 21000 116000 65000

E 2 en Mpa 11500 21000 7500 4900

12 0,28 0,26 0,32 0,34

2.6 Critères de limite d’élasticité

Les critères de résistance que nous allons définir représentent des valeurs limites pour les contraintes
maximales, et permettent de ce fait de garder un caractère élastique aux déformations.

2.6.1 Critère de Tresca

Il consiste à considérer de manière indépendante les trois contraintes de cisaillement maximal du tricercle de
Mohr. Soit en fonction des contraintes principales
 
Sup  I  II  I  III  II  III 
  20


 (4.6)

2.6.2 Critère de Von-Mises

1  
2 2
(   )  (   )  (  III )2   0
2  I II I III II (4.7)

ou encore

1
2

(11  22 )2  (11  33 )2  (22  33 )2  6(122  132  23
2
)   0

2.6.3 Critère de Hill

Le critère de Hill s’applique dans le cas de matériaux élastique orthotropes

F (11  22 )2  H (11  33 )2  G (22  33 )2  2L23

2 2
 2M 13 2
 2N 12 1
(4.8)
où F , H ,G, L, M , N sont des constantes fonctions des contraintes à ruptures.

Golay MMC - 52 -
 Elasticité

3 Le problème d’élasticité

3.1 Ecriture générale

Cinématique : + Equations de compatibilité
 

1
2

 u  T u 
 
u  U 0 (X ) sur U

Equilibre :
 
div   f  0 dans 
 
 F sur 
 n    F

R

sur U


Loi de comportement :

  Tr ()I  2

3.2 Formulation en déplacement

 
div   f  0
 
div (Tr ()I  2)  f  0
 
(Tr ())  2div()  f  0
    
(div u )  div( u )  div(T u )  f  0

soit l’équation de Navier

   
(  )(div u )  div( u )  f  0
(4.9)
Remarque: Si on prend la divergence de l’équation de Navier
 
(  2)(div u )  div( f )  0
  
Donc, si le champ de forces volumiques est tel que div f  0 alors div u est une fonction harmonique.

3.3 Formulation en contrainte

En partant de l’écriture des équations de compatibilité, on peut démontrer les équations de Michell

1    
div( )  (Tr ())  div f I   f  T f  0
1 1 (4.10)

Soit, si le champ de force est uniforme, on obtient les équations de Beltrami.

(1   )div( )  (Tr ())  0

(4.11)

3.4 Théorème de superposition

- 53 - Golay MMC
MMC

       
Si (U  f  F ) et (V  gG ) sont deux jeux de données engendrant respectivement des solutions u et v , alors
       
u  v est solution du problème de données (U  V   f   g F  G ) (Le problème est
évidemment linéaire).

3.5 Élasticité plane

3.5.1 Contraintes planes

Dans le cas où le chargement est dans le plan 12, la structure mince dans la direction 3, on peut faire
l’hypothèse que le problème est plan et libre de contraintes dans la direction 3.
Dans ce cas
 
11 (x 1 x 2 ) 12 (x 1 x 2 ) 0
 
  21 (x 1 x 2 ) 22 (x 1 x 2 ) 0
 
 0 0 0

et d’après la loi de comportement
  (x  x )  (x  x ) 0 
 11 1 2 11 1 2 

   21 (x 1 x 2 ) 22 (x 1 x 2 ) 0 

 
 0 0 33 (x 1 x 2 )

On remarquera que la déformation suivant 3 est non nulle.

3.5.2 Déformations planes

Dans le cas où le chargement est dans le plan 12, la structure très élancée dans la direction 3, sans possibilités
de déplacement suivant 3, on peut faire l’hypothèse que le problème est plan sous l’hypothèse des
déformations planes.
Dans ce cas
  (x  x )  (x  x ) 0
 11 1 2 12 1 2 
   21 (x 1 x 2 ) 22 (x 1 x 2 ) 0
 
 0 0 0

et d’après la loi de comportement
 (x  x )  (x  x ) 0 
 11 1 2 12 1 2 

  21 (x 1 x 2 ) 22 (x 1 x 2 ) 0 

 
 0 0 33 (x 1 x 2 )

On remarquera que la contrainte suivant 3 est non nulle.
D’après (4.3)

1   
33  33  11  22  33   0
E E
donc
 
33   11  22 

Nous allons prouver que les contraintes peuvent être déterminées par une seule fonction scalaire.
En appliquant l’équation d’équilibre (3.11) on a :

Golay MMC - 54 -
 Elasticité



111  122  0


  222  0

 211
donc






(x 1 x 2 )  11  2 et 12  1






(x1x2 )  21 2 et 22 1



comme le tenseur des contraintes est symétrique, on a 2  1  0 , donc

(x 1 x 2 )    2 et   1

en définitive on a prouvé


  22




 11


(x 1 x 2 )   11



 22


 12



 12



 est appelée la fonction d’Airy.

Le tenseur des contraintes devant vérifier l’équation de Beltrami (4.11), on a

(1   )ij kk  kk ij  0

d’où
  0

 est donc une fonction biharmonique.

3.6 Thermoélasticité

3.6.1 Thermodynamique : équations de bilan

Jusqu’à présent nous avons utilisé les équations de bilan suivantes :
Conservation de la masse

d 
  divv  0
dt
Conservation de la quantité de mouvement
 
div   f   dans 

Conservation du moment cinétique (3.17)

Nous introduisons maintenant l’équation de bilan de conservation d’énergie, ou encore le premier principe de
la thermodynamique :

d
dt
E  K   Pext  Q
où

E représente l’énergie interne E   e d  (e densité d’ énergie interne)



  
K représente l’énergie cinétique K   21 v  v d  ( v la vitesse)


- 55 - Golay MMC
MMC

   
Pext représente la puissance des efforts extérieurs Pext   f  v d    F  v d 
 

 
Q représente le taux de chaleur reçu Q   r d    q  n d  (q vecteur de chaleur et r source de chaleur)
 

Par application du premier principe, en utilisant (5.7) on a:

de        
 d    v   d    f  v d    F  v d    r d    q  n d 
 dt     

en utilisant la conservation de la quantité de mouvement (3.9)

de  
 d      d   r d   q  n d
 dt   

Soit, par application du théorème de la divergence, la forme locale du premier principe


e      r  divq
(4.12)
Nous présentons également, sans plus de discussion le second principe de la thermodynamique :
 
dS r q n
  d   d
dt T  T

où T est la température et S l’entropie. Ce second principe s’écrit sous sa forme locale


q r
s  div   0
T T (4.13)

où s représente l’entropie massique

3.6.2 Équation de la chaleur

On peut exprimer l’énergie interne massique e en fonction de l’entropie massique s , de la température T et
l’énergie libre  .
e    Ts
(4.14)
En thermoélasticité, sous l’hypothèse des petites perturbations, pour un écart de température par rapport à la
température au repos T  T0 petit, on a:

  (T )
Grace au second principe on peut montrer que

 
 et s 
 T

donc, le premier principe peut s’écrire

e    Ts

  Ts

et comme

   
     T     sT
 T 

on a

    sT  sT
  sT      r  divq

or
Golay MMC - 56 -
 Elasticité

   2  2   1  s 
s         T     T
 T  T 2
 T T
 T
c’est à dire que le premier principe s’écrit

 s  
T    T T  r  divq
T T

En introduisant la chaleur spécifique C  T s

T

 
T    CT  r  divq
T

puis la loi de Fourier q  kT , où k représente la conductivité thermique,


T    CT  r  divkT
T
En général la contribution mécanique est négligeable par rapport aux autres contributions, si bien que
l’équation de bilan de l’énergie conduit à l’équation de la chaleur :
 T  
CT  C   v  T   r  divkT
 t  (4.15)

Dans le cas où le problème à traiter est stationnaire, sans source de chaleur, avec une conductivité constante,
on retrouve l’équation habituelle :
T  0

3.6.3 Loi de comportement thermo-élastique

Dans le cadre de la thermoélasticité , l’énergie libre spécifique s’écrit comme un développement limité au
second ordre en déformation et température, ou plutôt en déformation et écart de température   T  T0
(supposés “ petits ”) :

1 1
(T )    C    s   b  2     
2 2
Par définition

  
 (T )  C      C     
 


où  représente le tenseur des dilatations thermiques

Dans le cas isotrope la loi de comportement thermo-élastique s’écrit :

  Tr ()1  2  (3  2)

- 57 - Golay MMC
MMC

4 A retenir

Loi de comportement élastique linéaire isotrope

   Tr () 1  2  

 1
  Tr () 1  
E E
Critère de Tresca
 

Sup  I  II  I  III  II  III 
  20
 

Le problème d’élasticité

  

 1
   u   u
  2
 T

u  U 0 (X ) sur U

  
div   f  0 dans 
  
  F sur 
 n    F
 
sur U
R

 

  Tr ()I  2


Equation de Navier
   
(  )(div u )  div( u )  f  0

En élasticité plane sous l’hypothèse des deformations planes :

 
33   11  22  et 33  0

Conservation de l’énergie

d
dt
E  K   Pext  Q
Forme locale de la conservation de l’énergie

e      r  divq

Equation de la chaleur
 T  
C   v  T   r  divkT
 t 

Loi de comportement thermoélastique isotrope

  Tr ()1  2  (3  2)

Golay MMC - 58 -
 Mécanique des fluides

INTRODUCTION A LA MECANIQUE DES FLUIDES

1 Loi de comportement
En mécanique des fluides, nous travaillerons toujours en variables d’Euler
Comme pour les matériaux solides (qui sont des fluides qui s’ignorent ..) les lois de comportement fluide sont
élaborées à partir de l’expérience.

Fluide viscoplastique
Fluide à seuil
τ
Fluide fluidifiant

Newton

Fluide
épaississant

dU
dy

1.1 Fluide Newtonien

Pour un fluide Newtonien, les contraintes sont une fonction affine des vitesses de déformation. Soit,

  pI  Tr (D )I  2D

(5.1)
où

1   T 
D
2

grad v  grad v  (5.2)

soit, en notation indicielle

1  vi v j 
Dij    
2  x j x i 

 est la viscosité dynamique (dimension Poiseuille  M )

LT
 est le second coefficient de viscosité

 2
On introduit également la viscosité cinématique   (dimension Stokes  L )
 T

- 59 - Golay MMC
MMC

1.2 Fluide incompressible

 
Si le fluide est incompressible, alors on a vu que divv  0 ou TrD  0

Donc, (5.1) devient

  pI  2D
(5.3)

1.3 Fluide non-visqueux

Si le fluide est parfait, alors on a ne tient pas compte de la viscosité, donc (5.1) devient

  pI
(5.4)
Le tenseur des contraintes est alors sphérique.
En particulier, l’action d’un fluide non visqueux sur une paroi est normale à la paroi (d’après l’équation
d’équilibre).

1.4 Fluide au repos

 
Si le fluide est au repos, alors v  0 , donc (5.1) devient

  pI
(5.5)

2 Conservation de la masse
La masse d’un système matériel qu’on suit dans son mouvement reste constante.
 dM
M  (t ) (x  t ) d  et 0
dt
où  est la masse volumique. On a alors (2.20)
 
d
dt
  divv  0 ou

t
 div v  0  (5.6)

Si on considère une grandeur différentiable  quelconque, on a alors pour un fluide incompressible

d d
(t )  dm  (t ) dm
(5.7)
dt dt
Démonstration:

d d
(t )  dm  (t )  d 
dt dt
d  
 (t )   div v  d 
 dt 
 d d 
 (t )    div v  d 
 dt dt 
 d d  
 (t )      div v   d 
 dt  dt  
 d  
 (t )   d
 dt 
d
 (t ) dm
dt
Soit  un domaine géométrique fixe traversé par le fluide,

Golay MMC - 60 -
 Mécanique des fluides

  
 v  n d    div(v ) d  (d aprés le théorème de la divergence)

   d (d après la conservation de la masse)
t

    d  (car  est fixé)
t

   dm
t
Si le fluide est incompressible, alors la masse volumique est constante et

 d 
 0
t dt

Si on note qm le débit massique à travers une surface S et qv le débit volumique, alors

   
qm   S v  n d     S v  n d   qv

Donc, en définitive:
 
Pour un domaine  fixe traversé par un fluide incompressible  v  n d   0 : le débit volumique à
travers la frontière  est nul.

3 Equation du mouvement
D’après l’équation du mouvement(3.10),

 dv
f  div   
dt
D’où, pour un fluide newtonien

dv   
  f  div pI  Tr (D )I  2D 
dt  
      
 f  div pI   div Tr (D )I   2div D 
     
   
 
 f  p   divv  div v  T v  
Soit l’équation de Navier-Stokes compressible

  

dv
dt
 
 f  p  (  ) divv  v
(5.8)
 
* Pour un fluide incompressible, divv  0 , donc (5.8) devient:
 
     
dv  v
     v  v   f  p  v

dt  t 
(5.9)
* Pour un fluide non visqueux, (5.8) devient:
 
   
dv  v 
    v  v   f  p
dt  t 

(5.10)

- 61 - Golay MMC
MMC

4 A retenir

Loi de comportement pour un fluide newtonien

  pI  Tr (D )I  2D

Conservation de la masse pour un fluide incompressible

 
divv  0
Grace à la conservation de la masse pour un fluide incompressible

d d
(t )  dm  (t ) dm
dt dt
Equation de Navier Stokes compressible

  

dv
dt
 
 f  p  (  ) divv  v

Golay MMC - 62 -
 Bibliographie

BIBLIOGRAPHIE

[1] Mécanique des Milieux Continus, Cours ESIM 1984, Equipe IMST Marseille.
[2] G. Duvaut, Mécanique des Milieux Continus, ed. Masson 1990.
[3] P. Germain - P. Muller, Introduction à la Mécanique des Milieux Continus, ed. Masson 1995.
[4] J. Salençon, Mécanique des Milieux Continus, ed. ellipse 1988.
[5] P. Germain, Mécanique, ed. ellipse, ecole polytechnique, tomes I et II.
[6] G. Dhatt, J.L. Batoz, Modélisation des structures par éléments finis: Solides élastiques, ed. Hermes, tome I.
[7] A. Bazergui, T. Bui-Quoc,A. Biron, G. McIntyre, C. Laberge, Résistance des matériaux, ed. de l’école
polytechnique de Montréal 1993.
[8] J. Coirier, Mécanique des Milieux Continus, ed. Dunos 1997.
[9] J. Lemaitre, J.L. Chaboche, Mécanique des matériaux solides, ed. Dunos 1996.
[10] O. Débordes, Thermodynamique des milieux continus, ESM2, Cours du DEA de Mécanique 2001.

- 63 - Golay MMC
MMC

Golay MMC - 64 -
 Annexes

ANNEXES: RAPPELS DE MÉCANIQUES DES SOLIDES RIGIDES

1 Cinématiques du solide
Avertissement: L’objectif de ce chapitre, est de familiariser les étudiants avec les notations tensorielles. Afin
d’en simplifier le contenu, nous ne considérerons que des bases orthonormées.

1.1 Description du mouvement

Soit S un ensemble de particules tel que la distances entre deux particules quelconques reste pratiquement
constante au cours du mouvement. On étudie l'ensemble S en le considérant indéformable: solide rigide.

1.1.1 Système de référence

  
Dans un espace euclidien  à trois dimensions, soit e1, e2 , e3 une base orthonormée. On définit un référentiel

d'observation par cet espace euclidien et le temps:  ,t  . On définit la dérivée d'un vecteur U par rapport
au temps dans ce référentiel par:

dU dU i 
 e
dt dt i


1.1.2 Mouvement d'un solide

  
Soit S un solide rigide en mouvement par rapport à  . Soit S O, e1, e2 , e3  un espace euclidien lié à S.
 
Considérons un vecteur lié à S, dont les composantes sont représentées par X dans  et X S dans S .

On note Q l'opérateur linéaire définissant le passage de  à S .

   
X S  QX et X  QT X S

Comme X S est indépendant du temps puisque lié à S, on a:

dX dQT  dQT  
 XS  QX  LX
dt dt 
dt 


or

QTQ  I

C'est à dire

dQT dQ
Q  QT 0
dt dt
L  LT  0

L est un opérateur linéaire antisymétrique, on peut donc définir un vecteur  (d'après (1.10)) tel que:

dX   
 LX    X
dt
 (6.1)
avec

- 65 - Golay MMC
MMC

 1 
  ijk Ljiek
2

1.1.3 Torseur cinématique

Soient A et P deux particules du solide S.
  
OP  OA  AP
donc par dérivation
   
dOP dOA dAP dOA  
      AP
dt dt dt dt
   

soit
   
V P   V A    AP
(6.2)

On définit le torseur cinématique par le vecteur vitesse de A par rapport à  , V A et le vecteur de rotation

 

instantanée  :   


V 

A




 
  

 


1.1.4 Accélération
   
V P   V A    AP
   
dV P  dV A d   dAP
   AP   
dt dt dt dt
   

Soit
 
           
 P    A 
d
dt

 AP      AP   A  
d
dt
 
 AP    AP   2 AP
  (6.3)

1.2 Composition des mouvements

L'espace temps est commun à tous les référentiels d'observation. on considère deux référentiels a  a , t et 
 
b  , t .

1.2.1 Dérivation composée

     
 
Soit U un vecteur dans la base  a O, e1a , e2a , e3a , U  U ieia par dérivation:
.

dU dU i  a
 e
dt a dt i


     
 
Soit U un vecteur dans la base  b O, e1b , e2b , e3b , U  U ieib par dérivation:
.
 
dU dU i  b deib
 e U i
dt dt i dt
a

Golay MMC - 66 -
 Annexes

Car  b est en mouvement par rapport à  a , et d'après (6.1)

 
deib 
 q /b  eib
dt

d'où, avec q /b représentant le vecteur rotation de  b par rapport à  a :
 
dU dU  
   q / b  U
dt dt
a b (6.4)

1.2.2 Composition des vitesses

Soit P  S :
  
O a P  O aO b  O b P
  
dO a P dO aO b dO b P
 
dt a
dt a
dt
  a

Et donc d’après (6.4)

  
dO a P dO aO b dO b P  
   q / b  O b P
dt a
dt dt
 a b

soit
  
V a (P )  V b (P )  V e (P )
    (6.5)
Vitesse / a Vitesse / b Vitesse d ' entrainement

1.2.3 Composition des accélérations

On dérive (6.5) par rapport à a :
     
d 2O a P d 2O aO b d 2O b P  dO b
P d q / b   dO b
P
b
   q / b    O P  q / b 
dt 2 dt 2 dt 2 dt dt dt
a
 a
 b b a a

    
 a (P )   a (O b )   b (P )  q / b V b (P ) 
           
 d q / b   dO b P
 
q b q b 
  /    /    O P   /   
b q b
 q / b  O b P 
 dt   dt 
 b   b 

     
d   q / b     
 a (P )   b (P )   a (O b )   O b P  q / b  q / b  O b P   2
q / b V b (P)
dt  
Accélération de Coriolis

b
Accélération d ' entrainement

   
 a (P )   b (P )   e (P )   c (P ) (6.6)

- 67 - Golay MMC
MMC

2 Cinétique

La cinématique ne s'intéresse au mouvement des corps que du point de vue de l'espace et du temps: durée,
vitesse, distance, etc …; tandis qu'en cinétique on introduit, en plus, le concept de masse, c'est à dire qu'on
tient compte aussi de la masse

2.1 Définitions
On définit la masse d’un solide S par :

m  S dm(P )  S (P, t )dv

(6.7)
Où  représente la densité volumique de masse.

On définit G le centre de masse (ou d’inertie) du solide S par :

 
 O   S OPdm(P )  mOG
(6.8)

2.2 Eléments de cinétique

2.2.1 Torseur cinétique

On définit le Torseur Cinétique ou Torseur des quantités de mouvement par :
 



 R  S V (P )dm(P ) Résultante cinétique de S /
    
 
k A  S AP V (P )dm(P ) Moment cinétique en A / (6.9)


On peut remarquer que si le repère  est fixe, alors :
 d  
R S OP (P )dm(P )  mV (G )
dt

2.2.2 Torseur dynamique

On définit le Torseur dynamique par :
 



 d  S (P )dm(P ) Résultante dynamique de S /
    
 
 A  S AP  (P )dm(P ) Moment dynamique en A / (6.10)



2.2.3 Relation entre torseur cinématique et torseur dynamique

En dérivant par rapport au temps dans  on obtient :

 dR
d
dt (6.11)
 
dk A dAP   
 S V (P )dm(P )  S AP  (P )dm(P )
dt dt
 
dAO  dOP   
 S V (P )dm(P )  S V (P )dm(P )  S AP  (P )dm(P )
 dt   dt
 V (A)  S V (P )dm(P )   A

 dk A  
A  V (A)  mV (G )
dt (6.12)

Golay MMC - 68 -
 Annexes

      

 A  S AP  (P )dm(P )  S AG  (P )dm(P )  S GP  (P )dm(P )
   
 A  AG  d  G

 dk G  
A   AG  m (G )
dt (6.13)

2.2.4 Energie cinétique

On définit l’énergie cinétique du solide S par :

1  2
T (S )  S V (P )dm(P )
2 (6.14)

2.2.5 Théorème de Koenig

     
Soit  O, e1, e2 , e3  un espace euclidien et G G, e1, e2 , e3  un espace euclidien barycentrique lié au solide S.
   
k A  k G  AG  mV (G )
(6.15)
   
 A  G  AG  m (G )
(6.16)
 
dOP dOP
T (S )  S  dm
dt dt
     
1 dOG dOG 1 dGP dGP dGP dOG
T (S )  S  dm  S  dm  S  dm
2 dt dt 2 dt dt dt dt

1
T (S )  T (S )  mV 2 (G )
G
2 (6.17)

2.3 Cinétique du solide rigide

2.3.1 Opérateur d’inertie

On définit l’opérateur d’inertie par  A tel que :
  
A :
 

u     A (u )   AP  u  AP dm 
 
Si AP  x iei alors
  
      
AP  u  AP  ijk x j kpq u p xqei  pi qj x j u p xqei  qi pj x j u p xqei  x j ui x jei  x j u j x iei

  
 
 



  

 
 
  
AP  u  AP  x j2ek  ek  uiei   x k x iek  ei   u je j  x j2ek  ek  x k x jek  e j  uiei  
Et l’opérateur d’inertie est représenté par la matrice :
 I I 12 I 13 
 1
I A  I 12 I2 I 23 
 
I 13 I 23 I3 

Où

- 69 - Golay MMC
MMC


I 1   x 22  x 32 dm  I 12   x 1x 2dm
I2   x 2
1
 x dm2
3
I 13   x 1x 3dm
I3   x 2
1
 x dm2
2
I 23   x 2x 3dm

2.3.2 Influence des symétries matérielles

 
• Si le solide S possède un plan de symétrie A, e1, e2  , alors

 x 1x 3dm   x 1x 3dm   x 1x 3dm   x 1x 3dm   x 1x 3dm  0
x 3 0 x 3 0 x 3 0 x 3 0

Soit
 I I 12 0 
 1
I A  I 12 I2 0 
 
 0 0 I3 


• Si le solide S possède un axe de symétrie A, e3  , alors

  dx 1dx 2dx 3    rdrd dx 3

 

Et comme

 x 1x 3 dx 1dx 2dx 3   r cos x 3 rdrd dx 3  0

 

2
 x 1x 2 dx 1dx 2dx 3   r sin  cos x 3 rdrd dx 3  0
 

on a finalement
I 0 0 
 1
I A   0 I2 0 
 
 0 0 I3 


• Moment d’inertie par rapport à une droite ∆ (de vecteur unitaire  ) passant par A
Soit H le projeté orthogonal d’un point P du solide S, on a alors :
 2   2  2    2   2 
 
2
I    PH dm    AP  AH dm     AP  AP . dm
S S    S   

              

     dm   . [Link]   AP AP.dm
2
I     . AP .AP  AP .
S   S

        


I   .  AP .AP  AP AP . dm
S     
        
 
Et comme a  b  c  a  c b  a  b c ,  
       

I   . AP    AP dm     A ()    I A  
S


Golay MMC - 70 -
 Annexes

• Théorème de Huyggens généralisé

  


 A (u )   AP  u  AP dm


    
 
  AG  u  AP dm   GP  u  AP dm
  
 
    

 
 
 AG  u   APdm   GP  u  AG dm   GP  u  GP dm
 
  
    
 
   
 AG  u  mAG   GPdm  u  AG   GP  u  GP dm   
Soit
 
 
 
 A (u )  AG  u  mAG   G (u ) 
• Théorème de Huyggens appliqué au moment d’inertie par rapport à une droite ∆ (de vecteur unitaire
   
 ) passant par A (tel que AG ⊥  ) et ∆’ (de vecteur unitaire  )passant par G.
              2     
 
I      A ()  I A     AG    mAG   G ()    mAG   m   AG AG  IG  
     
soit
     2   2
I   I A     mAG  IG    mAG  IG
 

2.3.3 Moment cinétique du solide

  
k A  S AP V (P )dm(P )

Soit Q un point quelconque du solide

    

k A  S AP  V (Q )    QP dm(P ) 
      

k A  mAG V (Q )  S AQ  QP    QP dm(P )   
        

k A  mAG V (Q )  AQ    S QPdm(P )  S QP    QP dm(P )   
D’où
      
k A  mAG V (Q )  mAQ    QG   G ()  
Dans le cas particulier où Q=G, on obtient :
   
k A  mAG V (G )   G ()

C'est-à-dire
 
k G   G ()

2.3.4 Energie cinétique du solide

1  2
T (S )  S V (P )dm(P )
2
Soit Q un point quelconque du solide
- 71 - Golay MMC
MMC

   2
T (S ) 
1
2
 
S V (Q )    QP dm(P )

      2
     
2
1 1
T (S )  S V (Q ) dm(P )  S V (Q )    QP dm(P )  S   QP dm(P )
2 2

m  2        

T (S ) 
2
 1
 
V (Q )  mV (Q )    QG  S   QP    QP
2  dm(P )
Soit

m  2    1  
T (S ) 
2
 
V (Q )  mV (Q )    QG     Q ()
2
Dans le cas particulier où Q=G, on obtient :

m  2 1  
T (S )  V (G )     G ()
2 2

3 Equations fondamentales de la mécanique des solides

3.1 Torseur associé aux efforts externes


Soit f (P ) une densité volumique de force exercée sur le solide S .

Soit F (P ) une densité surfacique de force exercée sur la frontière du solide S .

Soit F(P ) une densité linéique de force exercée sur une courbe  .

Soit F i une force ponctuelle exercée en un point Pi de S .

Le torseur des efforts extérieurs est défini par :

    


 
R   f (P )   F (P )   F (P )   F i
Fe (S )  
  S  S 
   i    

C A   AP  f (P )   AP  F (P )   AP  F(P )   AP  F i


 S S  i

3.2 Loi fondamentale de la dynamique

Il existe au moins un référentiel Galiléen associé à une chronologie, tel que :
S , t Torseur dynamique Torseur des forces extérieures

Ou encore
    
S , t  
 d,  A  Fe (R,C A )

En conséquence, on peut énoncer :

Théorème de la résultante dynamique : dans un référentiel galiléen
 
R  m (G )

Théorème du moment dynamique : dans un référentiel galiléen, soit A un point fixe


 dk A 
A  CA
dt

Golay MMC - 72 -
 Annexes

- 73 - Golay MMC
MMC

Golay MMC - 74 -