MMC Loi

Majeure Mécanique
UP1 - Mécanique des Matériaux
Thermodynamique des milieux continus et lois de

comportement
J. Bruchon
École des Mines de Saint-Étienne
Centre SMS
Octobre 2016
1
Table des matières
1 Cinématique 4
1.1 Loi du mouvement . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2 Le gradient de déformation . . . . . . . . . . . . . . . . . . . . 8
1.2.1 Produit mixte . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.2 Transport d'un élément de volume . . . . . . . . . . . 9
1.2.3 Transport d'un élément de surface orienté . . . . . . . 10
1.2.4 Décomposition polaire du gradient de la transformation 11
1.3 Mesures de déformation . . . . . . . . . . . . . . . . . . . . . 14
1.4 Tenseur des contraintes de Cauchy . . . . . . . . . . . . . . . 19
1.5 Lois de conservation dans la conguration courante . . . . . . 20
1.5.1 Conservation de la masse . . . . . . . . . . . . . . . . . 20
1.5.2 Conservation de la quantité de mouvement . . . . . . . 21
1.5.3 Conservation du moment de la quantité de mouvement 23
1.6 Tenseurs des contraintes alternatifs . . . . . . . . . . . . . . . 24
2 Eorts intérieurs - Théorème de l'énergie cinétique 26

2.1 Puissance des eorts intérieurs en description Eulérienne . . . 26
3 Conservation de l'énergie 29
3.1 Formulation locale du premier principe . . . . . . . . . . . . . 30
3.2 Récapitulatif . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
4 Formulation des lois de comportement 32

4.1 Variété des comportements - Lois 1D . . . . . . . . . . . . . . 32
4.2 Le principe d'indiérence matérielle ou d'objectivité . . . . . . 36
4.2.1 Dénitions : repère et référentiel . . . . . . . . . . . . . 36
4.2.2 Loi de transformation des tenseurs . . . . . . . . . . . 37
4.2.3 Invariance des lois de la MMC vis-à-vis d'un change-
ment de référentiel . . . . . . . . . . . . . . . . . . . . 39
4.2.4 Caractère intrinsèque de certaines variables en MMC . 40
4.2.5 Résultats d'opérateurs sur les grandeurs objectives . . 40
4.3 Les principes de la théorie du comportement mécanique des
matériaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.3.1 Déterminisme et fonctionnelle mémoire . . . . . . . . . 45
4.3.2 Application : uides visqueux . . . . . . . . . . . . . . 46
4.4 Fonctions isotropes de tenseurs . . . . . . . . . . . . . . . . . 47
4.4.1 Invariants d'un tenseur . . . . . . . . . . . . . . . . . . 47
4.4.2 Théorème de Cayley - Hamilton . . . . . . . . . . . . . 48
4.4.3 Théorèmes de représentation . . . . . . . . . . . . . . . 48
2
4.4.4 Cas de R3 . . . . . . . . . . . . . . . . . . . . . . . . . 49
4.5 Thermodynamique des milieux continus . . . . . . . . . . . . . 50
4.5.1 Second principe de la thermodynamique des milieux
continus . . . . . . . . . . . . . . . . . . . . . . . . . . 50
4.5.2 Loi de Fourier . . . . . . . . . . . . . . . . . . . . . . . 52
5 Hyperélasticité 53
5.1 Réversibilité des lois hyperélastiques . . . . . . . . . . . . . . 54
5.2 Hyperélasticité isotrope . . . . . . . . . . . . . . . . . . . . . . 55
5.2.1 Description matérielle (S, C) . . . . . . . . . . . . . . . 55
5.2.2 Descriptions matérielles alternatives . . . . . . . . . . . 56
5.2.3 Description spatiale . . . . . . . . . . . . . . . . . . . . 57
5.3 Matériaux Néo-hookéens . . . . . . . . . . . . . . . . . . . . . 59
3
1 Cinématique
1.1 Loi du mouvement
On considère un corps B qui occupe initialement (à t = t0 ) une région D0
R3 , de frontière régulière notée ∂D0 . Les sollicitations
de l'espace euclidien
appliquées à B s'accompagnent d'un changement de conguration depuis celle
de référence, vers la conguration courante notée Dt . On supposera ici que la
conguration de référence coïncide avec la conguration initiale D0 . Au cours
de cette évolution, un point matériel P appartenant à B passe de la position
X , repérée par ses coordonnées cartésiennes XI dans le repère EI , I = 1, 2, 3
à une position x repérée par ses coordonnées cartésiennes xi dans le repère
ei , i = 1, 2, 3. Il est souvent commode d'identier ei et EI , c'est-à-dire de
travailler dans le même système de coordonnées. La loi du mouvement est la
donnée d'une fonction vectorielle f bijective (non pénétration de la matière)
telle que
x = f (X, t; t0 ) (1.1)
(voir gure 1.1).
Figure 1.1: Description du mouvement

On introduit le vecteur de déplacement u qui est le vecteur qui lie un point
matériel de la conguration de référence et ce même point dans la congura-
tion courante :
u(X, t) = f (X, t; t0 ) − X = x − X (1.2)
Si le mouvement du corps est décrit en utilisant les coordonnées matérielles

XI , I = 1, 2, 3, on dira que la description est matérielle ou Lagrangienne. Cela
revient à observer le mouvement du corps comme celui d'une quantité de ma-
tière qui se déplace dans l'espace à partir d'une référence connue. Ce mouve-
ment peut également être décrit en utilisant les coordonnées xi , i = 1, 2, 3. On
dira alors que la description est spatiale ou Eulérienne. Cela consiste à obser-
ver la quantité de matière qui passe dans un certain domaine xé de l'espace,
4
sans référence à la position au temps initial. Dans ce qui suit, nous adop-
terons la convention d'écriture suivante : les lettres capitales seront utilisées
pour désigner les composantes des coordonnées Lagrangiennes (se trouvant
dans la conguration de référence), et les lettres minuscules, pour désigner
les composantes des coordonnées Eulériennes (qui se trouvent donc dans la
conguration courante).
Notons que la loi du mouvement permet de passer d'une représentation à
l'autre.
An de visualiser la diérence entre description lagrangienne (ou maté-
rielle) et description eulérienne (ou spatiale) du mouvement, considérons la
traction uniaxiale d'une barre de longueur initiale L = 2. La loi du mouve-
ment de cette barre est dénie par
x = (1 + t)X, avec 0≤X≤2
De plus, on considère une certaine distribution de température sur la barre,

donnée en description lagrangienne par un champ
T (X, t) = Xt2 , avec 0≤X≤2
La loi du mouvement permet de passer à une description eulérienne, avec un

champ de température déni sur la conguration courante
T (x, t) = xt2 /(1 + t)
Figure 1.2: Traction uniaxiale : description eulérienne et lagrangienne

La gure 1.2 montre l'élongation de la barre et l'évolution de sa tempéra-
ture. Cette gure met en avant le fait qu'une particule matérielle reste, tout
5
au long de son mouvement, repérée par sa coordonnée X (xe) de la congu-
ration initiale, tandis que sa coordonnée x évolue avec le temps. La tempé-
rature peut être obtenue par deux approches diérentes. Au temps t = 3, la
2
température de la particule repérée par X = 2, est T = 2 × 3 = 18. Alter-
nativement, cette même particule au temps t = 3 est à la position x = 8. La
2
température en cette position et à cet instant est T (8, 3) = 8×3 /(1+3) = 18.
Autrement dit, T (X, t) est la température au temps t de la particule repérée
par X dans la conguration initiale, tandis que T (x, t) est la température au
temps t et au point x de la conguration courante.
Intéressons nous maintenant au calcul de la variation de la température

au cours du temps, en congurations lagrangienne et eulérienne. Simplions
d'abord le problème ci-dessus, en disant que la température de chaque par-
ticule matérielle X est constante, et prenons
T (X, t) = X
Évidemment, nous avons dans ce cas, puisque X ne dépend pas de t,

dT ∂T
(X, t) = (X, t) = 0
dt ∂t
lorsque l'on travaille dans la conguration initiale.
Dans la conguration courante, le champ de température s'écrit mainte-
nant
x
T (x, t) = ,
1+t
et sa dérivée par rapport au temps vaut
∂T x
(x, t) = −
∂t (1 + t)2
ce qui est diérent de zéro. Que peut-on faire de cette dérivée ? Rien du tout.
Elle est en eet la limite du rapport
T (x, t + ∆t) − T (x, t)

∆t
lorsque ∆t tend vers zéro. Or, la particule se trouvant en x
t n'est plus la
à
même que celle se trouvant en x à t + ∆t, et ainsi la diérence T (x, t + ∆t) −
T (x, t) est non nulle du fait du mouvement du milieu continu. Les équations
de la mécanique résultent de bilans eectués sur des quantités (masse, éner-
gie, quantité de mouvement, ...) associées à une particule ou un groupe de
particules que l'on suit dans leur mouvement. Calculer la variation en temps
6
d'une quantité T associée à une particule, ne nécessite, en description la-
grangienne, que de calculer la dérivée partielle ∂T /∂t(X, t), mais demande,
en conguration eulérienne, de considérer la limite de
T (x + ∆x, t + ∆t) − T (x, t)

(1.3)
∆t
lorsque ∆t tend vers zéro, où x + ∆x est la position à t + ∆t de la particule
se trouvant en x à t. Ainsi, ∆x est donné par :
dx(t)
∆x = ∆t = v(x, t)∆t
dt
x
Dans notre exemple, la vitesse v vaut v(x, t) = 1+t
, et l'on peut vérier
que la limite de l'expression (1.3) est bien zéro.
Plus généralement, pour calculer cette limite on considère le développe-
ment
∂T
T (x + ∆x, t + ∆t) = T (x, t + ∆t) + (x, t + ∆t)∆x + o(∆x)
∂x
Ainsi, la quantité dT /dt
dT T (x + ∆x, t + ∆t) − T (x, t)

(x, t) = lim
dt ∆t→0 ∆t
T (x, t + ∆t) − T (x, t) ∂T
= lim + (x, t)v(x, t)
∆t→0 ∆t ∂x
∂T ∂T
= (x, t) + (x, t)v(x, t)
∂t ∂x
s'appelle dérivée particulaire de la quantité T exprimée en variables eulé-
rienne. Elle désigne, par dénition, la variation en temps d'une quantité T
attachée à une particule que l'on suit dans son mouvement. Pour un mouve-
3
ment dans R , elle s'écrit
3
dT ∂T X ∂T (x, t) dxi (t) ∂T
(x, t) = (x, t) + = + v · ∇T (1.4)
dt ∂t i=1
∂xi dt ∂t
et est la dérivée totale en temps de la fonction T (x(t), t) où x(t) est la

trajectoire d'une particule.
Notons que la dérivée particulaire de la ième composante de la vitesse
s'écrit, sur la conguration actuelle
dvi ∂vi
(x, t) = + v · ∇vi
dt ∂t
7
ce qui est non linéaire par rapport à la vitesse. Nous utiliserons cette expres-
sion lors de l'expression de la conservation de la quantité de mouvement.
Notons également que lorsque les équations de la mécanique sont appro-
chées par une méthode numérique comme les éléments nis, le domaine de
dénition spatial est discrétisé par un maillage qui peut avoir une vitesse vm
quelconque ( i.e. diérente de zéro et de v ). Ce mouvement du maillage doit
être pris en compte dans les équations en remplaçant v·∇T par (v−vm )·∇T
dans (1.4).
Enn, la dérivée particulaire de T exprimée en variables de Lagrange est
confondue avec sa dérivée partielle par rapport au temps :
dT ∂T
(X, t) = (X, t)
dt ∂t
1.2 Le gradient de déformation
Figure 1.3: Transformation d'un vecteur matériel élémentaire de la con-

guration de référence sur la conguration matérielle.
Considérons comme sur la gure 1.3 deux points matériels de la congu-

ration de référence,X et X + dX . Ces points forment un vecteur (ou bre)
matériel(le)X + dX − X = dX . Dans la conguration courante, ces points
correspondent, par la loi du mouvement, à x et x + dx, et la bre matérielle
devient dx. On a les relations suivantes, pour chaque composante i :
xi = fi (X, t; t0 )
et
xi + dxi = fi (X + dX, t; t0 )
Cette dernière relation permet d'écrire, au premier ordre,
3
X ∂fi
xi + dxi = fi (X, t; t0 ) + (X, t; t0 )dXk
k=1
∂X k
Ainsi, on en conclut que
∂fi ∂fi ∂fi

dxi = dX1 + dX2 + dX3
∂X1 ∂X2 ∂X3
8
pour toute composante i. Donc :
∂f1 ∂f1 ∂f1

 
∂X1 ∂X2 ∂X3

    
dx1  dX1

 dx2  =  ∂f2 ∂f2 ∂f1
 dX2 


dx3 ∂X1 ∂X2 ∂X3
 dX3
 
 ∂f3 ∂f3 ∂f3
∂X1 ∂X2 ∂X3
que l'on écrit
dx = F dX
F est le gradient de la transformation f . Ainsi, cette expression fait cor-
respondre à tout vecteur innitésimal dX de la conguration initiale, son
vecteur transporté par le mouvement dans la conguration à l'instant t, dx.
Cette expression est équivalente à
dx
F = (1.5)
dX
ou encore, en utilisant la relation x=X +u :
∂u
F = +I (1.6)
∂X
où I est le tenseur identité.
Le déterminant de F est noté J et est appelé jacobien de la transformation.
1.2.1 Produit mixte

Soient u, v et w trois vecteurs de R3 . Le produit mixte u · (v ∧ w), noté
[u, v, w], est égal au volume du parallélépipède formé par ces trois vecteurs
(voir gure 1.4 : aire du parallélogramme v∧w multipliée par la hauteur).
Ce produit est également égal au déterminant de (u, v, w) :
[u, v, w] = dét(u, v, w)
1.2.2 Transport d'un élément de volume

Considérons trois vecteurs (dX1 , dX2 , dX3 ) dénis dans la congura-
tion initiale D0 . Ces vecteurs forment un parallélépipède de volume dV =
[dX1 , dX2 , dX3 ]. Au cours du mouvement, ces trois vecteurs se transforment
9
Figure 1.4: Parallélépipède
en (dx1 , dx2 , dx3 ) et le volume initial se transforme en un volume dv donné

par
dv = [dx1 , dx2 , dx3 ]

= [F dX1 , F dX2 , F dX3 ]
= (détF )dét(dX1 , dX2 , dX3 )
Autrement dit,
dv = JdV avec J = détF > 0 (1.7)
Une transformation est dite isochore au point X si le volume est préservé

autour de X sous cette transformation. Si cette propriété est vériée en tout
X, la transformation est dite isochore (J ≡ 1).
1.2.3 Transport d'un élément de surface orienté

Soit, maintenant, dS une surface élémentaire construite à partir de deux
vecteurs élémentaires dX1 et dX2 . On a : dX1 ∧ dX2 = N dA, où dA est
l'aire de dS et N désigne la normale unitaire à dS , voir Figure 1.5.
Nous nous intéressons à savoir comment est transportée cette surface par
le mouvement. Soient les vecteurs transportés de dX1 et dX2 , dx1 = F dX1
et dx2 = F dX2 .
Nous voulons donc exprimer la surface ds construite à partir de dx1 et
dx2 . On a :
dx1 ∧ dx2 = nda
où da est l'aire de ds et n la normale unitaire à ds. Évaluons l'action de nda
10
Figure 1.5: Transport convectif d'une surface élémentaire
sur un vecteur v quelconque :
(nda) · v = (dx1 ∧ dx2 ) · v = [dx1 , dx2 , v]

= [F dX1 , F dX2 , v] = (detF )[dX1 , dX2 , F −1 v]
= J dA N · (F −1 v) = J(F −T N dA) · v
d'où, par identication :
nda = JF −T N dA (1.8)
1.2.4 Décomposition polaire du gradient de la transformation

Le tenseur gradient de la transformation F peut être décomposé en une
composante de déformation et une composante de rotation. Cela se fait grâce
à ce que l'on appelle la décomposition polaire. Il s'agit d'une décomposition
unique du tenseur gradient de déformation :
F = RU = V R (1.9)
avec U le tenseur de déformations pures à droite qui se trouve dans la con-

guration de référence, et V, le tenseur des déformations pures à gauche, qui
se trouve dans la conguration courante. R
est quant à lui le tenseur de
T
rotation local qui relie les deux congurations (détR = 1, RR = I ). Les
tenseurs de déformations à gauche et droite sont reliés par la relation (voir
gure 1.6) :
V = R U RT (1.10)
Les tenseurs U et V peuvent également s'écrire en fonction du tenseur F :
q q
V = F FT, U = FTF (1.11)
11
Les tenseurs de déformation pure jouent un rôle primordial car ils permettent
d'isoler la déformation pure, dans le gradient de la transformation. Cependant
en pratique, ils sont peu utilisés sous cette forme car leur extraction à partir
de F demande un nombre important d'opérations. On considérera en général
leur carré qui sont bien plus simples à évaluer.
Le carré de la longueur d'un segment dans la conguration courante peut
s'écrire :
dl2 = dxT dx = dX T F T F dX = dX T CdX (1.12)
avec C = FTF, ce qui en utilisant la relation (1.11) revient à C = U 2.

Ce tenseur est appelé tenseur des déformations de Cauchy-Green droit et il
s'agit d'un tenseur Lagrangien.
De la même façon, on peut exprimer le carré de la longueur d'un segment
dans la conguration de référence :
dL2 = dX T dX = dxT F −T F −1 dx = dxT b−1 dx (1.13)
T
où l'on a introduit le tenseur b = F F , ce qui en utilisant la relation (1.11)
2
revient à b = V . Ce tenseur est appelé tenseur des déformations de Cauchy-
Green gauche ou tenseur de Finger, et il s'agit d'un tenseur Eulérien.
Figure 1.6: Décomposition polaire du gradient de la transformation. Les

ellipses ont pour grand et petit axes les vecteurs propres respectifs de U et
V.
Exercice 1.1 On étudie la transformation faisant passer de la conguration

de la gure 1.7a à la conguration 1.7b.
1. Sachant que l'allongement relatif suivant e1 est l−ll0 0 = λ, écrire la
transformation correspondante dans la base canonique (e1 , e2 , e3 ).
2. En déduire les expressions des tenseurs F , U et R.
12
Figure 1.7: Cinématique de l'extension simple pour un corps carré muni
d'une grille : état initial (gauche) et nal (droite).
Exercice 1.2 On étudie la transformation illustrée sur la gure 1.8 qui

s'écrit :
x1 = X1 + γX2 ; x2 = X2 ; x3 = X3
où γ est l'amplitude du glissement.
Figure 1.8: Cinématique du glissement simple pour un corps carré muni

d'une grille : état initial (gauche) et nal (droite).
1. Donner l'expression de F , vérier que la transformation est isochore.

2. En déduire l'expression du tenseur Cauchy - Green droit C .
Les valeurs propres de C sont
1 p 1 p
λ1 = ( (γ + γ 2 + 4))2 ; λ2 = ( (−γ + γ 2 + 4))2 ; λ3 = 1
2 2
Les racines positives de ces valeurs propres fournissent les valeurs
propres du tenseur des déformations pures U . Les vecteurs propres
de C et U sont, par construction, identiques.
3. Que signie λ1 > 1 et λ2 < 1 ?
4. Vériez que U et R sont donnés (dans la base canonique) par :
   
α αγ/2 0 α αγ/2 0
2
[U ] =  αγ/2 α(1 + γ /2) 0 , [R] =
  −αγ/2 α 0 
0 0 1 0 0 1
avec α = 1/ 1 + (γ/2)2
p
5. À quel type de rotation R correspond ? Donnez l'expression de tan θ,

où θ est l'angle de rotation.
13
1.3 Mesures de déformation
Allongements - Soit M un vecteur unitaire décrivant la direction d'une
bre matérielle dans la conguration de référence. On note
dX = |dX|M
Dans la conguration courante, la longueur de la bre devient
|dx|2 = (F dX)T (F dX) = |dX|2 M T CM
On dénit alors le rapport d'allongement dans la direction M par :
|dx| v
λ(M ) = = uM T C M = |F M | = |U M |
u
|dX| t |{z}
F T F =U T U
Si M est choisi comme le ième vecteur de la base canonique {Ei }, i =

1, 2, 3, alors
q
p
λ(Ei ) = Cii = F1i2 + F2i2 + F3i2 (1.14)
ce qui donne une interprétation simple des composantes diagonales du ten-

seur de Cauchy - Green droit C , à savoir le carré de l'allongement du vecteur
de base correspondant.
Angle de glissement - Le cosinus de l'angle formé par deux vecteurs

unitaires est égal au produit scalaire de ces vecteurs. En reprenant les nota-
tions précédentes (voir gure 1.9), dans la conguration initiale on a, pour
deux directions matérielles M1 et M2 formant un angle Θ, cos Θ = M1 · M2 .
Cet angle devient θ dans la conguration actuelle :
(F M1 ) · (F M2 ) M1T CM2
cosθ = =
|F M1 ||F M2 | λ(M1 )λ(M2 )
On appelle γ = Θ − θ l'angle de glissement des directions M1 et M2 dans

le plan de glissement (M1 , M2 ). Si M1 = Ei , M2 = Ej , i 6= j , alors Θ = π/2
et
Cij
sin γ = cos(θ) = p (1.15)
Cii Cjj
ce qui donne une interprétation des composantes extra-diagonales du tenseur
de Cauchy-Green droit.
14
Figure 1.9: Angle de glissement : γ = Θ − θ.
Nous avons vu que grâce à la décomposition polaire il est possible de

séparer les rotations pures des élongations pures. Sous l'action de rotations
pures, la distance entre des particules voisines reste inchangée. L'élongation
pure caractérisée par U et V quant à elle, change la distance entre des
particules matérielles. An de quantier la distension, c'est-à-dire d'évaluer
de combien le tenseur U ou V s'écarte de l'identité I, il est nécessaire de
dénir une mesure de déformation. Cependant, il faut souligner qu'ici, tout
n'est qu'aaire de convention, et qu'il existe une multitude de mesures de
déformation. Dans la pratique, la mesure choisie sera celle qui permettra de
formuler de façon la plus simple la loi de comportement du matériau étudié.
Toutefois, il est possible de donner la dénition suivante. On appelle mesure
de déformation un tenseur ayant les propriétés suivantes :
1. Il est symétrique et sans dimension physique ;
2. Il est nul pour un mouvement de corps rigide et en F =I;

3. Son développement limité autour de F =I s'écrit
1 ∂u ∂u T ∂u
( +( ) ) + o( )
2 ∂X ∂X ∂X
Les mesures de déformation les plus utilisées appartiennent aux familles
suivantes (dénies par Hill, 1968 et 1978) :
- une famille de tenseurs des déformations lagrangiens, dénis relativement
à la conguration de référence et tels que
( 1
(U m − I), m 6= 0
Em = m (1.16)
ln U , m = 0
pour m=2 par exemple, on obtient le tenseur des déformations de Green-

Lagrange :
1
E = (F T F − I) (1.17)
2
15
- une famille de tenseurs des déformations eulériens qui sont dénis relative-
ment à la conguration courante et tels que :
( 1
(V m − I), m 6= 0
Am = m (1.18)
ln V , m = 0
À titre d'exemple, considérons un essai de compression ou extension sui-

vant la direction ex au cours duquel on suit la distance entre deux points
matériels séparés initialement par la distance l0 qui devient l à l'instant ac-
tuel. La déformation est supposée homogène (F ne dépend que de t) loca-
lement dans la zone où sont placés les deux points de mesure. Calculons les
expressions des diérentes mesures de déformations évoquées. Les scalaires
donnés ici sont les composantes (1,1) du tenseur utilisé. Remarquons d'abord
que l'information essentielle est donnée par le gradient de la transformation
F = l/l0 . Ensuite, nous avons
l l
C = ( )2 , b = ( )2 ,
l0 l0
et les déformations sont données par
l − l0 l − l0
E1 = , A−1 = ,
l0 l

1 l 2 1 l0 2
E2 = ( ) − 1 , A−2 = 1−( )
2 l0 2 l
l l
E0 = log , A0 = log
l0 l0
La mesure de déformation uniaxiale E1 qui se rapporte à la longueur ini-
tiale est souvent appelée déformation nominale ou déformation de l'ingénieur.
La déformation logarithmique, ou déformation naturelle, E0 est souvent uti-
lisée. En eet, lorsqu'on a à faire à des déformations importantes, il paraît
naturel de considérer la déformation totale comme la somme de petits incré-
ments de déformations. On a alors :
Z l
dl l
E0 = = log
l0 l l0
16
Figure 1.10: Évolution de diérentes mesures de déformation en fonction du
rapport l/l0 au cours d'une déformation pure homogène. Source : [FOREST].
Hypothèse des petites déformations - Enn, remarquons que le ten-

seur de Green-Lagrange (1.17), peut s'exprimer en fonction du vecteur des
déplacements u=x−X :
1
E= ∇X u + ∇TX u + ∇TX u∇X u (1.19)
2
Dans ce que l'on appelle la théorie géométriquement linéaire de la mécanique
des solides, les déformations du corps solide sont supposées être petites. Les
non-linéarités géométriques ne sont donc pas prises en compte. En négligeant
donc les contributions non-linéaires de la relation (1.19), celle-ci se réduit à :
1
E≈ ∇X u + ∇TX u = ε (1.20)
2
L'hypothèse des petites déformations ainsi faite équivaut à k∇X uk 1, ou
encore à F = O(I). Ainsi, le troisième point de la dénition d'une mesure de
déformation implique que toutes les mesures de déformation sont équivalentes
dans la théorie des petites déformations.
Exercice 1.3 En vous basant sur la relation (1.12), montrez la relation sui-
vante reliant la longueur dL d'une bre matérielle dans la conguration de
référence D0 , à sa longueur dl dans la conguration courante Dt :
dl2 − dL2
= M · EM (1.21)
2dL2
17
où M est le vecteur unitaire porté par la bre matérielle dans D0
Exercice 1.4 En vous basant sur la relation (1.13), et avec les mêmes no-
tations que précédemment, montrez que
dl2 − dL2
= m · A−2 m (1.22)
2dl2
où m est le vecteur unitaire porté par la bre matérielle dans Dt .
Le tenseur d'Almansi A−2 est, par dénition, égal à 21 (I − b−1 ).
Exercice 1.5 On considère la transformation suivante
1 1
x1 = (18 + 4X1 + 6X2 ) ; x2 = (14 + 6X2 )
4 4
qui agit comme représenté sur la gure 1.11
Figure 1.11
1. Calculez F et F −1 .
2. Calculez comment se transforment les vecteurs (matériels) de la base
initiale. Calculez les vecteurs réciproques de la base courante.
3. Calculez les tenseurs de Cauchy-Green droit et gauche, C et b, en
déduire les mesures de déformation de Green - Lagrange et d'Almansi.
4. En vous basant sur la relation (1.21), calculez la déformation de Green
- Lagrange associée au vecteur matériel dX2 = (0, 1)T .
5. En vous basant sur la relation (1.22), calculez la déformation d'Al-
mansi associée avec l'élément dx2 , transformé de dX2 .
18
Exercice 1.6 Un solide subit une déformation uni-axiale donnée par
x1 = X1 (1 + βt)
où β est un paramètre arbitraire.

Calculez :
1. Le gradient de la transformation F .
2. Le tenseur de Cauchy - Green droit C .
3. La dilatation selon les axes X1 et X2 .
4. L'angle entre ces axes au cours de la déformation.
5. Le tenseur des déformations de Green - Lagrange.
6. La déformation selon les deux axes.
7. Le tenseurs petites déformations ε.
1.4 Tenseur des contraintes de Cauchy

On considère un corps continu déformable B . On dénit le vecteur contrainte
t, dans la conguration courante, qui est tel que :
df
t = lim (1.23)
ds→0 ds
df , qui sont exercés sur une

Il caractérise les eorts intérieurs de cohésion
partie de B à travers un élément de surface ds, de normale extérieure n (voir
gure 1.12).
D'après le théorème de Cauchy, en tout point, et à chaque instant t, le vecteur
des contraintes t (x, t, n) dépend de la normale n de façon linéaire. Il existe
donc un tenseur de second ordre appelé tenseur des contraintes de Cauchy,
noté σ qui est tel que :
t (x, t, n) = σ (x, t) · n (1.24)
Figure 1.12: Contraintes de Cauchy
19
Il peut être utile, comme nous le verrons plus tard, de décomposer le
tenseur des contraintes en une partie dite déviatorique, qui est de trace nulle
et une partie dite sphérique, de sorte que :
σ = σ D − pI (1.25)
avec σD la partie déviatorique de σ, et pI sa partie sphérique. p est appelée

la pression hydrostatique et est dénie par :
1
p = − tr σ (1.26)
3
1.5 Lois de conservation dans la conguration courante

Cette section décrit les lois de conservation de la masse, de la quantité de
mouvement et du moment de la quantité de mouvement. Celles-ci sont des re-
lations fondamentales en mécanique des milieux continus que nous utiliserons
pour mettre en place notre formalisme.
1.5.1 Conservation de la masse

Ce principe exprime le fait que la masse d'un système matériel ω ⊂D
que l'on suit dans son mouvement reste constante. En introduisant ρ0 et ρ
la masse volumique du milieu dans sa conguration de référence et actuelle
respectivement, la conservation de la masse s'exprime par :
Z Z
m(ω) = ρ(x, t) dv = ρ0 dV = constante
ω ω0
Or nous avons déjà vu que dv = JdV , avec J = détF , d'où
ρ(x, t)J = ρ0 (X) (1.27)
ce qui correspond à la forme lagrangienne de la conservation de la masse.

En dérivant cette relation par rapport au temps le long de la trajectoire
x(t) d'une particule matérielle, on obtient l'expression eulérienne de la loi de
conservation de la masse locale :
∂ρ
+ ∇ · (ρv) = 0 (1.28)
∂t
∂J
où l'on a utilisé le fait que
∂t
= J∇ · v . L'équation (1.28) est aussi appelée
équation de continuité, et peut directement être vue comme une équation
de bilan exprimant la variation locale de la masse volumique sous le ux
convectif ρv .
20
1.5.2 Conservation de la quantité de mouvement
Nous donnons d'abord un résultat mathématique. Soient F (x, t) une fonc-
tion de l'espace et du temps, et ω(t) un domaine borné dépendant du temps.
Cela signie que sa frontière ∂ω évolue au cours du temps avec une certaine
vitesse v. Alors :
Z Z Z
d ∂F
F (x, t) dv = (x, t) dv + F (x, t)v · n ds (1.29)
dt ω(t) ω(t) ∂t ∂ω(t)
où n est le vecteur unitaire normale à ∂ω et pointant à l'extérieur de ω.

R
Ce résultat s'interprète ainsi : la variation en temps de la quantité
ω
F dv
(membre de gauche) est due, d'une part à la variation de F dans le domaine
ω, puisque F dépend de t, Fv
et d'autre part à un ux sur le bord ∂ω . Ce
ux traduit le fait que le mouvement de ∂ω à la vitesse v fait entrer ou sortir
du domaine une certaine quantité de F .
De plus, rappelons le théorème de la divergence. Soit q un vecteur déni

sur ω (dépendant ou non de t). Alors :
Z Z
q · n ds = ∇ · q dv (1.30)
∂ω ω
En appliquant le théorème de la divergence à l'intégrale surfacique de

l'équation (1.29) avec q = F v, cette expression se réécrit en :
Z Z
d ∂F
F (x, t) dv = + ∇ · (F v) dv (1.31)
dt ω(t) ω(t) ∂t
Revenons à la conservation de la quantité de mouvement. Soit P (t) la

quantité de mouvement du domaine matériel ω(t), dénie par :
Z
P (t) = ρv dv (1.32)
ω(t)
La loi de conservation de la quantité de mouvement postule que la dérivée

par rapport au temps de la quantité de mouvement est égale à l'ensemble
des forces qui s'exercent sur le volume matériel :
Z Z
dP (t)
= ρb dv + t ds (1.33)
dt ω(t) ∂ω(t)
où b représente les forces volumiques, t = σ·n représente les forces de contact

s'exerçant sur ∂ω , la surface de ω ⊂ D.
21
Développons la ième composant du membre de gauche en utilisant de
théorème (1.31) :
Z Z
d ∂ρvi
ρvi dv = + ∇ · (ρvi v) dv
dt ω(t) ω(t) ∂t
Or,
∂ρvi ∂vi ∂ρ dvi

+ ∇ · (ρvi v) = ρ + vi + vi ∇ · (ρv) +ρv · ∇vi = ρ
∂t ∂t | ∂t {z } dt
=0,éq.(1.28)
dvi
où
dt
= ∂v
∂t
i
+ v · ∇vi est ce que l'on a déni comme la dérivée particulaire
de vi . En notation tensorielle, on écrit
Z Z Z
d dv ∂v
ρv dv = ρ dv = ρ + (∇v) · v dv
dt ω(t) ω(t) dt ω(t) ∂t
avec [∇v]ij = ∂vi /∂xj .
De plus, une extension du théorème de la divergence (montrez-le en uti-

lisant la notation indicielle) dit que :
Z Z
σ · n ds = ∇ · σ dv
∂ω | {z } ω
t
où la divergence du tenseur du second ordre σ est un vecteur dont la ième

composante vaut :
3
X ∂σij
[∇ · σ]i =
j=1
∂xj
Ainsi, l'équation (1.33) se réécrit en
Z
dv
ρ − ρb − ∇ · σ dv = 0 (1.34)
ω(t) dt
Ceci devant être vérié pour tout volume matériel ω, on en déduit que l'in-
tégrande doit être nulle. Ainsi, la forme locale de la loi de conservation de la
quantité de mouvement s'écrit sur le domaine courant :
dv ∂v
ρ = ρ( + (∇v) · v) = ∇ · σ + ρb (1.35)
dt ∂t
22
1.5.3 Conservation du moment de la quantité de mouvement
On considère un point xe x0 et le vecteur position r = x − x0 . On note
H le moment de la quantité de mouvement par rapport à r qui est déni
par : Z
H= ρr ∧ v dv (1.36)
ω(t)
Le principe de conservation du moment de la quantité de mouvement exprime

le fait que la dérivée temporelle du moment de la quantité de mouvement
relative au point xe x0 est égale à la somme des moments agissant sur le
corps. L'équation s'écrit alors :
Z Z
dH
= r ∧ (ρb) dv + r ∧ t ds (1.37)
dt ω(t) ∂ω(t)
On peut montrer que la forme locale de la loi de conservation de la quantité

de mouvement implique :
σ = σT (1.38)
Ce qui signie que le tenseur des contraintes de Cauchy est symétrique. Cette
condition de symétrie a une interprétation physique très simple : regardant
un cube élémentaire autour du point M, la condition d'équilibre ∇·σ = 0
n'est pas susante pour garantir l'équilibre. Il faut également s'assurer de
l'absence de rotation. Si l'on regarde par exemple les composantes σ12 et σ21
sur toutes les faces verticales du cube élémentaire, nous constatons sur la
gure 1.13 qu'elles créent deux couples antagonistes qui doivent forcément
s'équilibrer pour éviter toute mise en rotation du cube autour de l'axe e3 .
D'où σ12 = σ21 et la symétrie de σ.
Figure 1.13: Équilibre en rotation autour de l'axe e3 d'un cube élémentaire

sous contrainte.
23
1.6 Tenseurs des contraintes alternatifs
Un matériau est dit élastique si il existe un état de référence sans contrainte,
et si après déformation le tenseur des contraintes ne dépend que du tenseur
des déformations calculé par rapport à cet état de référence. Cependant, le
1 T
tenseur des déformations de Green - Lagrange E = (F F − I) est déni en
2
termes de variables de Lagrange, tandis que σ , le tenseur des contraintes de
Cauchy est exprimé en termes de variables d'Euler. L'objectif de cette section
est donc d'introduire un tenseur des contraintes en variables de Lagrange.
Soit ω ⊂ D un domaine matériel occupant la position ω0 à l'instant

initial et la position ωt à l'instant t. Soit t(M, t, n) la contrainte au point M
à l'instant t pour la direction n. La force élémentaire df qui s'exerce sur un
élément de surface ds, d'aire da, passant par M et de normale n, est égale
à :
df = tda = σ · nda
Nous cherchons à exprimer df à l'aide de l'élément de surface dS d'aire
dA de la conguration de référence. Or, d'après (1.8), on a :
nda = JF −T N dA
D'où
df = Jσ(F −T N )dA
De plus, comme σ(F −T N ) = (σF −T )N , on est amené à dénir le tenseur
suivant :
Π = JσF −T (1.39)
appelé tenseur des contraintes de Boussinesq ou de Piola-Lagrange.

Il vient alors :
df = ΠN dA
Les composantes de Π s'écrivent :
3
X ∂Xj
Πij = J σip
p=1
∂xp
Ce tenseur nous permet d'écrire de façon simple les équations du mou-

vement en variables de Lagrange. Pour ce faire, nous partons du bilan de la
quantité de mouvement en description Eulérienne :
Z Z Z
d
ρv dv = t da + ρbdv
dt ωt ∂ωt ωt
24
Or, on vient de montrer que df = tda = ΠN dA. On peut donc réécrire
l'équation de bilan sur la conguration initiale :
Z Z Z
dv
ρ0 dV = ΠN dA + ρ0 b dV
ω0 dt ∂ω0 ω0
Par le théorème de la divergence et comme ceci est valable pour tout

domaine ω, on en déduit l'équation locale Lagrangienne du mouvement :
dv ∂v(X, t)
ρ0 = ρ0 = ∇ · Π + ρ0 b (1.40)
dt ∂t
Cependant, le tenseur de Piola-Lagrange présente deux inconvénients.

D'une part il n'est pas symétrique. D'autre part, on remarque, par la rela-
tion df = ΠN dA, que Π appliqué à un vecteur élémentaire d'aire dans la
conguration initiale donne la force s'exerçant sur l'élément de surface dans
la conguration actuelle. Ainsi, Π ne satisfait pas pleinement l'objectif initial
qui était d'exprimer les contraintes dans la conguration initiale.
Pour aller plus loin, on introduit une force ctive dF attachée à l'élément
de surface dS de la conguration initiale, dénie par :
df = F dF ou dF = F −1 df
Ceci signie que df se déduirait de dF par transport convectif. Dans la

réalité, la force élémentaire df n'a aucune raison de se comporter comme un
vecteur transporté dans le mouvement.
Nous partons de F dF = ΠN dA. En remplaçant Π par son expression :
F dF = JσF −T N dA
D'où :
dF = J(F −1 σF −T )N dA
On est amené à introduire le tenseur suivant, dit tenseur des contraintes
de Piola-Kirchho :
S = JF −1 σF −T (1.41)
Et nous avons
dF = SN dA
qui est l'expression symétrique de
df = σ · nda
25
La diérence essentielle est que dF n'a pas de réalité mécanique. En compa-
rant ces deux expressions à
df = ΠN dA
nous soulignons le caractère purement Eulérien de σ, purement Lagrangien

de S et hybride de Π. De plus, le tenseur S est symétrique. Cependant S n'a
aucune réalité physique. Son intérêt réside dans le fait qu'il est le tenseur des
contraintes conjugué du tenseur des déformations de Green-Lagrange comme
nous le verrons dans la section suivante. Seuls les tenseurs σ et Π possèdent
un sens physique et peuvent caractériser les eorts appliqués. De plus, la
composante (1,1) du tenseur de Boussinesq en traction uniaxiale dans la di-
rection e1 , correspond à la contrainte dite ingénieur c'est à dire : Π1 = F/S0
où F est la force appliquée sur la section d'aire S0 .
Remarquons aussi que puisque S = F −1 Π, alors Π = F S, et l'équation

(1.40) devient :
dv ∂v(X, t)
ρ0 = ρ0 = ∇ · (F S) + ρ0 b
dt ∂t
Enn, lorsque l'hypothèse des petites déformations est faite, c'est-à-dire
lorsque les termes en O(k∇uk) peuvent être négligés, les tenseurs σ, Π et
S sont équivalents.
2 Eorts intérieurs - Théorème de l'énergie ci-

nétique
2.1 Puissance des eorts intérieurs en description Eulé-
rienne
La forme locale de la conservation de la quantité de mouvement s'écrit :
dv
ρ = ∇ · σ + ρb
dt
En multipliant cette équation par la vitesse v il vient :
1 dv 2
ρ = v · ∇ · σ + ρv · b (2.1)
2 dt
D'autre part, Exercice :
∇ · (σv) = (∇ · σ) · v + σ : ∇v
26
On rappelle, que la divergence d'un tenseur du second ordre σ est un
vecteur, noté ∇ · σ, dont la i ième composante est dénie par
3
X ∂σij ∂σij
[∇ · (σ)]i = =
j=1
∂xj ∂xj
où on a utilisé la convention de sommation de l'indice répété j dans la dernière

expression.
De plus, le produit doublement contracté entre deux tenseurs du deuxième
ordre est déni par : A : B = Apq Bqp .
Le gradient des vitesses se décompose en la somme d'une partie symé-

trique D et une partie antisymétrique Ω :
∇v = D + Ω
avec
1 1
D = (∇v + (∇v)T ) et Ω = (∇v − (∇v)T )
2 2
On montre de plus ( exercice) que si A est symétrique et B antisymé-
trique, alors A : B = 0. Il s'en suit, σ étant symétrique, que
σ : ∇v = σ : D + σ : Ω = σ : D
Ainsi, on a montré la relation
∇ · (σv) = ∇ · (σ) · v + σ : D
1 dv 2
et ρ = v · ∇ · σ + ρv · b devient :
2 dt
1 dv 2
ρ = ∇ · (σv) − σ : D + ρv · b
2 dt
Considérons maintenant une partie ω d'un système matériel D, et inté-
grons la relation précédente sur ce domaine à un instant t :
1 dv 2
Z Z Z Z
ρ dv = ∇ · (σv)dv − σ : Ddv + ρv · bdv
ωt 2 dt ωt ωt ωt
Le premier terme peut être réécrit (les termes supplémentaires liés à la dé-
pendance temporelle de ω s'annulent grâce à l'équation de conservation de
la masse) en :
Z
d 1 2
ρv dv
dt ωt 2
27
Le deuxième terme peut être transformé en une intégrale de surface :
Z Z
∇ · (σv)dv = (σv) · nda
ωt ∂ωt
D'où :
Z Z Z Z
d 1 2
ρv dv = (σv) · nda − σ : Ddv + ρv · bdv
dt ωt 2 ∂ωt ωt ωt
Exercice : Montrez que (σv) · n = (σ · n) · v.

En introduisant t=σ·n le vecteur contrainte en M, on obtient :
Z Z Z Z
d 1 2
ρv dv = t · vda − σ : Ddv + ρv · bdv (2.2)
dt ωt 2 ∂ωt ωt ωt
Interprétons les termes de cette équation :

La grandeur Z
1 2
K= ρv dv
ωt 2
est l'énergie cinétique associée au domaine ω à l'instant t.
Posons Z Z
Pext = t · vda + ρv · bdv
∂ωt ωt
On reconnaît la puissance des eorts extérieurs au domaine ω à t. Ces
eorts sont décomposés en eorts de contact s'exerçant sur la surface
de ω et en eorts à distance s'exerçant au sein de ω.
Reste le terme Z
Pint = − σ : Ddv
ωt
qui est la puissance des eorts intérieurs. Le terme −σ : D =
−σij Dij est donc la puissance volumique (eulérienne) des eorts in-
térieurs. La puissance Pint est nulle dans tout mouvement rigidiant,
puisque on a alors D = 0.
Nous venons de montrer que la dérivée particulaire à l'instant t de l'énergie

cinétique associée à un domaine ω est égale à la somme de la puissance des
eorts extérieurs et de la puissance des eorts intérieurs. Ce que l'on peut
résumer par :
dK
= Pext + Pint (2.3)
dt
28
Remarquons que cette équation ne contient aucune loi physique nouvelle.
Elle dérive directement de la conservation de la quantité de mouvement.
Dans cette section, nous avons déni la puissance des eorts intérieurs
sur la conguration courante. Par la loi du mouvement, on peut se ramener
sur la conguration initiale, et l'on a, pour un domaine matériel ω (on ne
donne pas la démonstration) :
Z Z
−Pint = σ : D dv = Jσ : D dV
ωt ω0
Z
= Π : Ḟ dV
ω0
Z
= S : Ė dV
ω0
Ceci signie que la densité des eorts intérieurs rapportée à l'élément de

volume initial peut s'exprimer en fonction de plusieurs couples de variables
constitués d'un tenseur de contraintes et d'une certaine vitesse ou taux de
déformation :
Jσ : D = S : Ė = Π : Ḟ (2.4)
Pour chaque couple, la contrainte et la déformation associée sont dites conju-

guées : E S sont conjugués, ainsi que Π et F . Par contre, le gradient des
et
déformation D n'étant la dérivée d'aucune déformation, σ et D sont associés,
mais aucune déformation n'est conjuguée avec σ .
3 Conservation de l'énergie
La conservation de l'énergie constitue le premier principe de la thermo-
dynamique. Les diérentes contributions à l'énergie totale du domaine ω à
l'instant t sont : Z
1
L'énergie cinétique K= ρv 2 dv
ωt 2
La puissance extérieure des eorts appliqués :
Z Z
Pext (v) = t · vda + ρb · vdv
∂ωt ωt
L'énergie interne E du système, représentée par une densité massique

e d'énergie interne :
Z
E= ρe(x, t)dv
ωt
29
L'apport de chaleur Q au système sous la forme d'un apport surfacique
h(x, t, ∂ω) et d'un apport volumique ρr(x, t) d'origine non mécanique
(rayonnement, ...) :
Z Z
Q= hds + ρrdv
∂ωt ωt
Les grandeurs scalaires introduites e, h et r ne dépendent pas du référen-

tiel d'observation. On dit alors qu'elles sont objectives. L'énergie cinétique,
par contre, est liée au référentiel choisi, car elle fait intervenir la vitesse. La
première loi de la thermodynamique dans un mou-
stipule alors que,
vement observé dans un référentiel galiléen, la variation d'énergie interne
et cinétique d'un corps matériel dans la conguration courante est due à la
puissance mécanique des eorts extérieurs et à l'apport de chaleur :
Ė + K̇ = Pext + Q (3.1)
En utilisant le théorème de l'énergie cinétique (2.3), K est éliminée de

cette équation, et l'on fait apparaître les eorts intérieurs :
Ė = −Pint (v) + Q (3.2)
De même que le vecteur contrainte t ne dépend que de la normale n à la

surface, on postule ici que l'apport surfacique de chaleur h(x, t, ∂ω) ne dépend
de la surface que via sa normale sortante (principe de Fourier - Stokes) :
h(x, t, ∂ω) = h(x, t, n) = −q(x, t) · n
L'existence du ux de chaleur q est une conséquence directe de la dépen-

dance de h à la normale n. Le signe - est une convention : q représente le
ux entrant et n la normale sortante. De plus, l'invariance supposée de h par
changement de référentiel, implique que q est un vecteur objectif.
L'équation du bilan d'énergie (3.2) prend ainsi la forme suivante :
Z Z Z Z
ρėdv = σ : Ddv − q · nda + ρrdv (3.3)
ωt ωt ∂ωt ωt
3.1 Formulation locale du premier principe

An d'aboutir à une formulation locale de la conservation de l'énergie,
on part de (3.3), qui est valable pour tout sous-domaine matériel ω, et l'on
30
transforme, par le théorème de la divergence, l'intégrale de surface en une
intégrale de volume :
Z Z
q · nda = ∇ · qdv
∂ωt ωt
D'où : Z
(ρė − σ : D + ∇ · q − ρr)dv = 0
ωt
Cette égalité étant valable pour tout sous-ensemble matériel ω, et sous la

condition que les champs impliqués soient susamment réguliers (continuité),
on en déduit la forme locale du premier principe de la thermodyna-
mique :
ρė = σ : D − ∇ · q + ρr (3.4)
On voit donc apparaître un apport d'énergie mécanique, via le terme

σ : D. Soulignons que ė est la dérivée particulaire de l'énergie interne, que
l'on peut expliciter en :
∂e ∂ρe
ρė = ρ( + v · ∇e) = + ∇ · (ρev)
∂t ∂t
où l'on s'est servi de l'équation de conservation de la masse (1.28) ∂ρ/∂t +
∇ · (ρv) = 0 pour obtenir la dernière expression.
Par suite, l'équation (3.4) se réécrit en :
∂ρe
+ ∇ · (ρev + q) = σ : D + ρr
∂t
3.2 Récapitulatif
Si on récapitule les formes locales des équations de bilan de la thermo-
mécanique des milieux continus (en l'absence de discontinuité), on a :
Conservation de la masse : ρ̇ + ρ∇ · v = 0
Conservation de la quantité de mouvement : ∇ · σ + ρb = ργ
T
Conservation du moment cinétique : σ = σ
Conservation de l'énergie : ρė = σ : D − ∇ · q + ρr
Ce qui fait 15 inconnues (f (3), σ (6), ρ (1), e (1), r (1) et q (3)), pour
seulement 5 équations ! Il est donc nécessaire d'introduire des relations entre
ces grandeurs an de pouvoir clore le système. C'est le rôle des équations
d'état et des lois de comportement que nous allons maintenant introduire.
31
4 Formulation des lois de comportement
Concentrons nous dans un premier temps sur le problème mécanique.
Celui-ci comporte une équation scalaire (bilan local de la masse) et une équa-
tion vectorielle (bilan local de la quantité de mouvement), soit au total, 4
équations. Les inconnues du problème sont la masse volumique ρ, la trans-
formation f (X, t), et la contrainte σ en chaque point matériel, ce qui fait
10 inconnues. En l'état actuel le problème est donc mal posé, puisque il
manque 6 équations. Ces équations, à trouver, permettent d'exprimer les 6
σ en fonction des inconnues fon-
composantes du tenseur des contraintes
damentales du problème ρ et f . Ces équations, contrairement aux lois de
conservation, n'ont pas un caractère universel. Leur forme dépend de la na-
ture et de l'état du matériau étudié. Elles sont appelées lois de comporte-
ment, et caractérisent les propriétés mécaniques du matériau étudié. Elles
doivent donc être élaborées en forte relation avec l'expérience.
Nous allons voir cependant, que la forme de ces lois n'est pas arbitraire.
Certains principes, considérés comme (quasi-) universels imposent certaines
restrictions sur leur forme. De plus, nous verrons que ces lois mettent né-
cessairement en jeu la grandeur température. Nous aurons alors recours à
deux autres lois universelles, les première et seconde lois de la thermodyna-
mique, pour formuler les lois de comportement, ce qui exigera de connaître
le comportement thermique et entropique de la matière.
4.1 Variété des comportements - Lois 1D

En observant la réponse mécanique des matériaux à des sollicitations
diverses telles que les essais uniaxiaux (traction, compression), les essais de
cisaillements plans, etc. (voir gure 4.1), on distingue trois grandes classes
de comportements.
Élasticité. Le comportement élastique linéarisé uniaxial est modélisé
par le comportement d'un ressort dont l'allongement ∆l est propor-
tionnel à la charge appliquée F . Cette analogie n'est valable que dans
un certain domaine de charge, dit domaine d'élasticité. Le facteur de
proportionnalité est la souplesse 1/k , où k est la raideur du ressort
−1
(en N.m ou P a.m) :
F = k∆l
Viscosité. Le comportement visqueux uniaxial est celui d'un amor-
tisseur ou d'un piston dont la réponse dépend de la vitesse de défor-
mation, ici de manière linéaire pour simplier :
F = η∆l˙
32
Figure 4.1: Contrainte en fonction de la déformation dans un matériau
soumis à un essai de traction
où η est le coecient de viscosité de l'amortisseur, exprimé en P a.s.m

ou N.s.m−1 (la viscosité d'un uide s'exprime en P a.s).
Plasticité. Le comportement rigide parfaitement plastique uniaxial

se caractérise par l'existence d'une contrainte seuil en-deçà de laquelle
aucune déformation n'est observée et pour laquelle une déformation
permanente quelconque est possible. Ce comportement est modélisé
par analogie avec un patin posé sur un substrat rugueux, susceptible
de se déplacer à volonté pour une charge appliquée susante. La des-
cription de ce comportement fait appel à une fonction critère
f (F ) = |F | − F0
et le glissement (écoulement) du patin s'écrit
∆l˙ = 0 si f <0 ou si (f = 0 et f˙ < 0)
∆l˙ 6= 0 si (f = 0 et f˙ = 0)
Ces comportements élémentaires sont résumés sur la gure 4.2. Il est
possible de combiner ces comportements entre eux, ce qui conduit à des
comportements plus complexes rencontrés dans la pratique : viscoélasticité,
élastoplasticité, viscoplasticité, et élastoviscoplasticité.
33
Figure 4.2: Schématisation des trois comportements élémentaires dans le
cas uniaxial.
A titre d'exemple, on considère le modèle viscoélastique (pour les uides)

de Maxwell représenté sur la gure 4.3 et constitué d'un amortisseur et d'un
ressort montés en série. Dans cette conguration, l'eort F se transmet
d'un élément à l'autre du montage, tandis que les vitesses de déformation
s'ajoutent. On a alors :
˙ = ∆l˙amort (t) + ∆l˙ressort (t) = F (t) Ḟ (t)

∆l(t) +
η k
Figure 4.3: Viscoélasticité uniaxiale : modèle de Maxwell

En intégrant cette équation diérentielle, on obtient la réponse du modèle
à une histoire de sollicitation ˙ allant de s = 0 où le système était supposé
∆l(s)
au repos à l'instant actuel s = t :
Z t
k ˙
F (t) = k exp(− (t − s))∆l(s)ds
0 η
Ainsi, la réponse actuelle en force dépend en général de l'histoire com-
plète de déformation du matériau. La loi de comportement est donc ici une
34
fonctionnelle de l'histoire du matériau. Remarquons toutefois que cette his-
toire est pondérée par le terme exponentiel qui privilégie les instants proches
de l'instant présent t. Cette perte de la mémoire des sollicitations anciennes
est caractéristique de la viscoélasticité.
Un deuxième modèle de comportement viscoélastique est obtenu en consi-

dérant cette fois un ressort et un amortisseur montés en parallèle comme sur
la gure 4.4. Ce modèle est le modèle de Kelvin - Voigt décrivant le compor-
tement visco-élastique de certains solides.
Figure 4.4: Viscoélasticité uniaxiale : modèle de Kelvin - Voigt

Dans un montage en parallèle, la déformation ∆l est globale tandis que
les contraintes s'ajoutent. Ainsi,
F = k∆l + η∆l˙
Par exemple si à t=0 la déformation est nulle, et qu'une force F0 est

alors appliquée au matériau, on obtient la déformation suivante en fonction
du temps :
F0
∆l(t) = l(t) = (1 − e−λt )
k
où 1/λ = η/k est un temps de relaxation. Ainsi, ce comportement tend
vers un comportement purement élastique lorsque e−λt → 0. Si à présent la
contrainte exercée sur le matériau est annulée au temps t1 , alors le matériau
se décharge et tend à recouvrer son état initial à déformation nulle, suivant
la loi (voir gure 4.5)
l(t > t1 ) = l(t1 )e−λ(t−t1 )
35
Figure 4.5: Chargement - déchargement d'un matériau de type Kelvin -
Voigt
4.2 Le principe d'indiérence matérielle ou d'objecti-

vité
Les lois physiques sont soumises à un certain nombre de principes d'inva-
riance lors d'un changement de repère ou de référentiel. Ainsi, une loi issue
de la mécanique classique doit être invariante vis-à-vis d'une transformation
Galiléenne (translation rectiligne uniforme à vitesse constante), ou encore,
les équations de l'électromagnétisme et de la mécanique relativiste sont in-
variantes par transformation de Lorentz.
Il en va de même pour les lois constitutives de la mécanique des mi-
lieux continus (MMC). Celles-ci doivent être indépendantes, non seulement
de l'orientation des axes (ce qui est fait pas une écriture tensorielle), mais
aussi du mouvement du système de référence utilisé ( i.e. du mouvement de
invariance
l'observateur). Cette dernière indépendance est équivalente à l'
des lois constitutives par superposition d'un mouvement de corps
rigide (translation et rotation) au mouvement du matériau. On dit
alors que ces lois doivent obéir au principe d'objectivité.
Pour obéir à ce principe et être donc invariantes par changement de réfé-
rentiel, les équations constitutives de la MMC doivent être écrites en termes
de quantités indépendantes de l'observateur (indépendantes de l'orientation
des axes) et de son mouvement. De telles quantités sont dites objectives. À
titre d'exemples, la position d'un point dépend de l'observateur et n'est donc
pas une grandeur objective. Il en va de même pour la vitesse qui dépend du
mouvement de l'observateur. Par contre, les distances et les angles sont des
grandeurs objectives (dans le cadre de la mécanique classique).
4.2.1 Dénitions : repère et référentiel

An d'eectuer des observations d'un phénomène physique, il est néces-
saire de disposer d'un référentiel que l'on appellera également observateur.
36
Le référentiel est analogue à un corps rigide muni d'une horloge. Dans ce réfé-
rentiel, l'origine des temps, ainsi qu'un système de coordonnées, encore appelé
repère, sont choisis librement. Ce repère, généralement choisi orthonormé,
matérialise le référentiel et permet d'identier une particule P à un instant
t à l'aide de ses coordonnées. Il faut donc bien faire la distinction entre un
changement de repère, lié à un changement de système de coordonnées,
qui engendre une description mathématique diérente du phénomène étudié ;
et un changement de référentiel, lié à un changement d'observateur, qui
engendre une vision diérente du phénomène étudié.
Nous cherchons la transformation la plus général permettant de passer

d'un référentiel à un autre. Cette relation doit conserver les angles et les dis-
tances (mouvement de corps rigide) et s'exprime donc comme la composition
à chaque instant, d'une translation et d'une rotation. Ceci signie que si un
point matériel P est repéré par ses coordonnées x à un instant
dans un t
∗
certain référentiel, alors ce même point P se repéré par ses coordonnées x
∗ ∗ ∗
à un instant t dans un second référentiel, et la relation entre x et x , t et t
s'écrit :
x∗ (P , t∗ ) = Q(t) x(P , t) + c(t)
(4.1)
t∗ = t − τ
où Q(t) est une matrice de rotation caractérisée par
QT (t)Q(t) = I, et détQ(t) =1
Par exemple une rotation d'angle θ(t) autour de l'axe z s'exprime par
 
cos θ(t) − sin θ(t) 0
Q(t) =  sin θ(t) cos θ(t) 0 
0 0 1
Si les horloges des deux référentiels sont synchronisées, ce que l'on sup-
posera dans la suite, alors τ = 0, donc t = t∗ .
4.2.2 Loi de transformation des tenseurs

Examinons comment se transforment les grandeurs tensorielles sous l'ef-
fet du changement de référentiel (4.1).
Un scalaire a qui représente une quantité physique intrinsèque (masse,

énergie interne, température, ... en MMC) est invariant par la transformation
dénie en (4.1) :
a∗ = a (4.2)
37
Remarquons que la norme de la vitesse n'est pas un scalaire au sens précé-
dent puisque c dépend de t.
Soit un vecteur matériel b, déni par deux points matériels. Dans le

premier référentiel, les coordonnées de ces points sont notées x1 et x2 , et les
composantes du vecteur valent b = x2 − x1 . Dans un second référentiel, ces
∗ ∗
coordonnées deviennent x1 et x2 , et les composantes du vecteur deviennent
b∗ = x∗2 − x∗1 . La relation entre b et b∗ s'écrit :
b∗ = Q(t)b (4.3)
En eet,
b∗ = x∗2 − x∗1 = (Qx2 + c) − (Qx1 + c) = Q(x2 − x1 ) = Qb
Si nous considérons un tenseur du second ordre C , il devient par la

transformation (4.1) :
C ∗ = Q(t)C Q(t)T soit Cij∗ = Qia Qjb Cab (4.4)
En eet, considérons deux vecteurs matériels u et v reliés par u = Cv .

Lors du changement de référentiel, nous avons
u∗ = Qu et v ∗ = Qv
et la relation entre les bres matérielles représentées par u et v doit rester
satisfaite dans le second référentiel. On doit donc avoir
u∗ = C ∗ v ∗
soit
Qu = C ∗ Qv
ou encore
u = QT C ∗ Qv = Cv
d'où l'on déduit
C ∗ = Q C QT
Ceci peut se généraliser à des tenseurs d'ordre n, pour lesquels on
trouve
∗
Zijk··· = Qia Qjb Qkc · · · Zabc··· (4.5)
Les relations (4.2), (4.3), (4.4) décrivent la transformation d'un scalaire,

d'un vecteur reliant deux points matériels, et d'un tenseur reliant deux vec-
teur matériels, sous le changement de référentiel (4.1).
38
4.2.3 Invariance des lois de la MMC vis-à-vis d'un changement
de référentiel
En plus du principe d'invariance par changement de repère à t xé (Q =
cste, c = cste) assuré par l'écriture tensorielle, les équations constitutives
doivent être formulées de telle sorte que l'expression de la réponse du maté-
riau soit indépendante de l'observateur et de son mouvement, ou, de manière
équivalente, indépendante d'un changement de référentiel (4.1) pour la con-
guration courante. Pour la conguration de référence, les deux référentiels
entre lesquels on fait le changement, sont supposés se correspondre. Cette
condition d'indiérence matérielle est appelée principe d'objectivité.
Le principe d'objectivité conduit à n'introduire dans les équations de
comportement que des grandeurs objectives, i.e. se transformant comme
les coordonnées lors d'un changement de référentiel (4.1), et donc ayant une
signication physique intrinsèque indépendante du mouvement du référen-
tiel. En d'autres termes, une grandeur tensorielle d'entrées a, bi , Cij , Dijk ,
··· est dite objective si elle se transforme suivant les lois (4.2) - (4.5), non
seulement vis-à-vis d'un changement de repère (changement d'axes avec c et
Q constants), mais aussi vis-à-vis d'un changement de référentiel, où cette
fois c(t) et Q(t) dépendent du temps.
f , vectorielles g et tenso-
En résumé, si les fonctions constitutives scalaires
rielles H , dépendant par exemple des grandeurs objectives a, b et C , vérient
f (a∗ , b∗ , C ∗ ) = f (a, Qb, QCQT ) = f (a, b, C)
g(a∗ , b∗ , C ∗ ) = g(a, Qb, QCQT ) = Qg(a, b, C) (4.6)
H(a∗ , b∗ , C ∗ ) = H(a, Qb, QCQT ) = QH(a, b, C)QT
alors, l'équation constitutive correspondante respectera le principe d'objec-

tivité :
F (a, bi , Cij , · · · ) = 0 ⇒ F (a∗ , b∗i , Cij∗ , · · · ) = 0 (4.7)
pour tout changement de référentiel.
Quelques remarques :
Le principe d'objectivité n'est pas une conséquence du principe d'in-
variance de Galilée de la mécanique classique (ou l'on ne considère
que le cas Q = I ).
Seules les équations constitutives de la MMC doivent vérier le prin-
cipe d'objectivité. Les équations de conservation ne doivent pas né-
39
cessairement y obéir. Par contre, leur écriture dans un référentiel non
Galiléen doit être correctement posée (composition des vitesses et des
accélérations : force de Coriolis, etc.)
4.2.4 Caractère intrinsèque de certaines variables en MMC

Il est admis que les grandeurs suivantes ont un caractère intrinsèque et
sont donc objectives.
la masse volumique, l'énergie interne, l'énergie libre, l'entropie, la tem-
pérature, l'enthalpie, la production de chaleur par unité de volume et
de temps, se comportent comme des scalaire objectifs.
Les forces de volume, les forces de contact et le vecteur ux de chaleur,
se comportent comme des vecteurs objectifs.
4.2.5 Résultats d'opérateurs sur les grandeurs objectives

1. Les combinaisons linéaires de grandeurs objectives sont objectives.
2. Le produit de deux grandeurs objectives est objectif.
3. La puissance entière d'une grandeur objective est objective.
4. Le double produit contracté de deux grandeurs objectives est objectif :
A∗ij Bji
∗
= Qia Qjb Aab Qjc Qid Bcd = δad δbc Aab Bcd = Aab Bba
car QT Q = I , i.e. Qia Qib = δab

5. Si a est un scalaire objectif (a∗ = a) alors son gradient ∇a est
objectif. En eet
∂a∗ ∂a ∂xj ∂a ∂[QT x∗ ]j ∂a
∗
= ∗
= ∗
= Qij
∂xi ∂xj ∂xi ∂xj ∂xi ∂xj
Les dérivées temporelles de grandeurs spatiales objectives ne

sont pas nécessairement des grandeurs objectives. Ainsi,
1. La dérivée d'un vecteur objectif déni sur la conguration courante
b∗ = Qb n'est pas objectif puisque l'on a
ḃ∗ = Q̇b + Qḃ 6= Qḃ
2. De même, la dérivée temporelle d'un tenseur spatial objectif déni sur

∗ T
la conguration courante C = Q C Q n'est pas objective puisque
∗
Ċ = Q̇ C QT + Q Ċ QT + Q C Q̇T 6= Q Ċ QT
40
La raison de cette non-objectivité est que la diérence entre deux
tenseurs relatifs à deux instants diérents dépend du mouvement du
système de référence (en particulier de sa rotation, via le terme Q̇)
entre ces deux instants. Le problème est que la variation en temps du
tenseur des contraintes de Cauchy σ apparaît dans de nombreuses lois
de comportement. On doit alors fabriquer une dérivée temporelle qui
est objective (voir exercice 4.2).
Grandeurs géométriques
1. Le temps n'est pas une grandeur objective. En eet, t = t−τ 6= t. Le
∗
temps ne peut donc pas intervenir explicitement dans l'écriture d'une

loi de comportement.
2. On vérie aisément que le vecteur position x, la vitesse v, ou encore

l'accélération ne sont pas des vecteurs objectifs.
3. Un vecteur reliant deux points est objectif. On en déduit, à l'aide

du produit scalaire, que la distance entre deux points et l'angle entre
deux vecteurs sont des scalaires objectifs. De même, la normale à une
surface est un vecteur objectif.
Exercice 4.1 Exemples de grandeurs objectives et non-objectives

Grandeurs cinématiques
1. On appelle tenseur hybride du second ordre un tenseur T qui, à un
vecteur A de la conguration de référence, fait correspondre un vec-
teur a de la conguration actuelle : a = T A. Montrez que si a est un
vecteur objectif, alors T ∗ = QT . Un tenseur hybride se transformant
de la sorte est dit objectif.
On remarque ici que l'objectivité de T ne s'exprime pas comme dans

l'expression (4.4). Ceci provient du caractère hybride de T qui est
composé d'une partie eulérienne (dénie sur la conguration courante)
et d'une partie lagrangienne (dénie sur la conguration de référence
et donc non aectée par le changement de référentiel). Les relations
(4.4) et (4.5) ne sont vraies que pour des tenseurs entièrement dénis
sur la conguration courante.
2. Montrez que le tenseur gradient de la transformation F est ob-
jectif.
3. En déduire l'objectivité du jacobien J .
41
Il en va de même pour la mesure de la surface matérielle dS . Ainsi,
outre les angles et les longueurs, les mesures de volumes et surfaces
sont préservées par changement de référentiel.
4. On dénit le tenseur gradient de vitesse (eulérien) par : L = ∇v .
Montrez d'abord que ce tenseur peut s'écrire sous la forme
L = Ḟ F −1
puis que l'on a
L∗ = Q L QT + Q̇QT
d'où la non objectivité de ce tenseur.
5. Montrer que le tenseur des vitesses de déformation D, partie sy-
métrique de L est objectif.
Ceci montre que la combinaison linéaire de grandeurs non objectives

peut conduire à une grandeur objective. En revanche, la partie antisy-
métrique de L, Ω, n'est pas objective.
Déformations
1. Montrez que le tenseur de Cauchy - Green droit C est objectif, ainsi
que le tenseur de Green - Lagrange E .
2. Montrez que le tenseur de Cauchy - Green gauche b = F F T et le
tenseur d'Almansi A−2 = 12 (I − b−1 ) sont objectifs. Montrez que ceci
implique que A−2 et A∗−2 expriment le même changement intrinsèque
de longueur.
Contraintes
1. En partant de la dénition du tenseur de Cauchy σ, montrez que
celui-ci est objectif, à savoir :
σ ∗ = Q σ QT
2. En déduire que me tenseur des contraintes de Piola - Lagrange

Π et le tenseur des contraintes de Piola - Kirchho S sont
objectifs.
Exercice 4.2 Dérivée temporelle de la contrainte de Cauchy
Un matériau hypoélastique est un matériau pour lequel l'évolution des
contraintes de Cauchy est décrite par :
σ̇ = f (σ)D
42
où D est le tenseur des vitesses de déformation. Dans l'hypothèse des petites
déformations, l'expression ci-dessus est viable, puisque les rotations sont né-
gligées. Nous allons voir qu'en grandes transformations, la non-objectivité de
σ̇ conduit à des aberrations.
(a) Élongation pure suivant e1 et e2 . (b) Élongation pure suivant e01 et e02 et
rotation autour de e3 .
Figure 4.6: Déformation plane d'une bre décrite dans deux référentiels
diérents. Source : [PONTHOT].
1. Rappelez pourquoi σ̇ n'est pas objectif.

On considère la structure représentée sur la gure 4.6a. Celle-ci est sou-
mise à une déformation plane. Les positions courantes xi sont données en
fonction des positions initiales XA par
x1 = h1 (t)X1 , x2 = h2 (t)X2 , x3 = X 3
où les facteurs hi sont indépendants des coordonnées. La déformation est

donc homogène et hi représente le rapport entre la longueur courante et la
longueur initiale d'une bre matérielle alignée avec l'axe ei .
2. Calculez le tenseur gradient de la transformation F . En déduire le
tenseur gradient des vitesses de déformation L, puis les tenseurs des
taux de déformation D et des taux de rotation Ω.
Nous considérons maintenant la même structure, mais soumise à un mou-
vement de rotation autour de e3 , d'une vitesse angulaire θ̇. L'analyse faite
ci-dessus est celle d'un observateur lié au référentiel de la bre.
3. Écrire maintenant dans le référentiel du laboratoire (voir Figure 4.6b)

les positions courantes x1 et x2 en fonction de X1 et X2 . Vous vous
aiderez de la matrice Q, matrice de rotation des axes d'une bre ma-
térielle entre l'instant initial et l'instant actuel.
43
4. En déduire F L , le gradient de la transformation dans le référentiel du
laboratoire.
5. Montrez que LL = Q L QT + ω . Vous expliciterez ω en fonction de θ̇.
6. En déduire la forme de DL et ΩL .
Avec la convention choisie pour Q, on a σ = QT σ L Q, et σ L = Q σ QT .
7. Montrez que : σ˙L = Q σ̇ QT + ω σ L − σ L ω
Supposons maintenant que la structure soit précontrainte (σ 6= 0 à t = 0)

et tourne sans se déformer (h1 = h2 = 1 et θ̇ 6= 0). Cette précontrainte peut
être induite par une force appliquée sur la barre et qui tourne avec celle-ci.
De manière évidente on a :
D = DL = 0
La loi hypoélastique donne

σ̇ = f (σ)D = 0, et σ̇ L = f (σ L )DL = 0
8. Montrez que ceci entre en contradiction avec l'expression trouvée à la

question précédente.
Conclusions de l'exercice : le terme parasite ω σL − σL ω représente
un accroissement ctif et donc indésirable de contraintes dû à la rotation
relative des deux référentiels. Il est donc impératif, dans la loi hypoélastique,
de n'utiliser que des grandeurs objectives. Nous avons utilisé le tenseur des
vitesses de déformation D au lieu du tenseur du gradient des vitesses L. Nous
devons également substituer à la dérivée temporelle σ̇ une dérivée temporelle
∇
objective notée σ. Les lois hypoélastiques sont donc reformulées en
∇
σ = f (σ)D
Soulignons encore une fois que dans le cadre des petites déformations
les rotations sont négligeables, donc Q=I et Q̇ = 0. Dans ce cas, ω =0
et l'inconsistance est levée. Dans le cas général, il faut dénir une dérivée
objective. Un moyen simple est de poser :
∇
σ = σ̇ − Ω σ + σ Ω
44
où Ω = ṘRT est le taux de rotation, où R est la rotation issue de la décom-
position F = RU .
Cette dérivée est dite dérivée de Jaumann, car introduite par Jaumann
en 1903. On peut montrer qu'elle est objective. Il est important de noter
que ce n'est pas la seule dérivée temporelle de tenseur objective qui existe.
L'addition à cette dérivée de produits de tenseurs objectifs (D et σ) fournit
une dérivée objective.
4.3 Les principes de la théorie du comportement méca-

nique des matériaux
4.3.1 Déterminisme et fonctionnelle mémoire
On cherche à déterminer le tenseur des contraintes de Cauchy σ en un
point matériel x = f (X, t) d'un domaine ω dans sa conguration courante
à l'instant t. Le principe du déterminisme dit que la contrainte σ est déter-
minée par l'histoire présente et passée du corps matériel, i.e., l'histoire de sa
déformation entre l'origine s = 0 et l'instant actuel s = t. Nous avons vu
que l'expression de la loi de comportement devait, pour respecter le principe
d'objectivité, ne faire intervenir que des grandeurs objectives, ce qui exclut
une dépendance directe à f , t ou encore v . On suppose alors l'existence d'une
fonctionnelle mémoire F telle que :
∂ 2f ∂ nf

∂f
σ(x, t) = F (X, s), (X, s), · · · , (X, s) (4.8)
0≤s≤t,n>0 ∂X ∂X 2 ∂X n
avec x = f (X, t). Remarquons que cette expression réduit la dépendance
spatiale de la contrainte à un voisinage arbitrairement petit autour du point
matériel X étudié.
Dans la pratique, on vérie en général qu'il est susant de se limiter au
cas n = 1. On parle alors de théorie du premier gradient et le milieu
continu ainsi modélisé est qualié de matériellement simple. On rappelle
∂f
que l'on a déjà introduit le premier gradient de f, noté F = . D'où :
∂X

σ(x, t) = F F (X, s) (4.9)
0≤s≤t
L'objectivité du tenseur des contraintes, σ ∗ = Q σ QT , se traduit par la

condition suivante sur la fonctionnelle mémoire :

F∗ F (X, s) QT

Q(s)F (X, s) = Q F
0≤s≤t 0≤s≤t
45
pour toute matrice orthogonale Q.
Classiquement, le principe d'invariance de forme est utilisé, stipulant que
la forme de la loi de comportement est invariante par changement de réfé-
∗
rentiel : F (·) = F(·). Ce qui entraîne,

F (X, s) QT

F Q(s)F (X, s) = Q F (4.10)
0≤s≤t 0≤s≤t
Remarquons que le principe d'invariance implique en fait que F ne doit

dépendre directement que de grandeurs objectives. Par exemple, si nous
∗
considérons un scalaire a non objectif, alors a = α + a avec α 6= 0 un
réel. D'où, nous devrions avoir :

F (X, s), a QT

F Q(s)F (X, s), α + a = Q F
0≤s≤t 0≤s≤t
En dérivant l'expression ci-dessus par rapport à α, nous en déduisons que

∂F
= 0.
∂α
Enn, notons qu'en plus de la condition (4.10), l'expression d'une loi de
comportement doit également respecter les symétries matérielles du milieu
étudié. Ceci se traduit par l'invariance de F sous certaines transformations
dépendant du matériau étudié.
4.3.2 Application : uides visqueux

Étendons le modèle 1D de comportement visqueux vu plus haut au cas
général tridimensionnel. La fonctionnelle mémoire ne peut dépendre du gra-
dient des vitesses L qui est non-objectif. Elle va donc dépendre de sa partie
symétrique, D, le tenseur des vitesses de déformation, qui est une quantité
objective. Pour un uide, nous écrivons donc :
σ(x, t) = F(D(x, t))
Et l'on doit avoir :
QF(D)QT = F(QDQT )
Il est alors possible de montrer (théorème de représentation d'une fonction

isotrope dépendant d'un tenseur) qu'une telle relation implique nécessaire-
ment que :
σ = α0 I + α1 D + α2 D 2
où les αi dépendent ici de la masse volumique ρ et des invariants de D
(i.e. trD , ((trD) − tr(D ))/2 et detD ). On obtient ainsi la loi des uides
2 2
46
visqueux de Reiner-Rivlin, établie en 1945. La linéarisation par rapport à D
de la loi précédente conduit à la loi des uides de Navier-Stokes :
σ = −p(ρ)I + η1 (trD)I + 2η2 D
où p est l'état de contrainte lorsque D = 0, et η1 , η2 sont les paramètres de

viscosité.
4.4 Fonctions isotropes de tenseurs

Dans cette section, nous nous intéressons aux fonctions scalaires, vecto-
rielles ou tensorielles, ayant comme arguments des tenseurs, et dites isotropes.
Nous avons déjà rencontré de telles fonctions, puisque ce sont les fonctions
vériant les propriétés (4.6) traduisant le principe d'objectivité. Ainsi, l'iso-
tropie d'une fonction est mathématiquement dénie comme suit :
Une fonction scalaire f (b, C) est dite isotrope vis-à-vis du groupe
d'isotropie Q si la relation
f (b, C) = f (Qb, QCQT )
est valable pour tout tenseur orthogonal Q ∈ Q. En particulier, le

groupe orthogonal complet de R ,
3
i.e
l'ensemble des tenseurs du se-
cond ordre dont la valeur absolue du déterminant est unitaire, est noté
SO(3) (Special Orthogonal Group).
Une fonction vectorielle g(b, C) est dite isotrope vis-à-vis du groupe

Qg(b, C) = g(Qb, QCQT )
est valable pour tout tenseur orthogonal Q ∈ Q.

Une fonction tensorielle H(b, C) est dite isotrope vis-à-vis du groupe
QH(b, C)QT = H(Qb, QCQT )
est valable pour tout tenseur orthogonal Q ∈ Q.
4.4.1 Invariants d'un tenseur

Une fonction scalaire isotrope d'un seul tenseur C est appelée un inva-
riant de C . La trace et le déterminant d'un tenseur sont des invariants de
47
ce tenseur. Une classe importante d'invariants d'un tenseur C sont les inva-
riants principaux, notés en général Ik (C), avec k = 1, · · · , n (n étant la
dimension de l'espace), et dénis comme les coecients du polynôme
det(C + θI) = θn + I1 (C)θn−1 + · · · + In−1 (C)θ + In (C)
avec en particulier
I1 (C) = tr C et In (C) = det(C) (4.11)
D'autres invariants importants sont les moments du tenseur C, dénis

par
I k = tr C k , k = 1, · · · , n (4.12)
On montre qu'il est possible d'exprimer les invariants principaux en fonction

des moments et inversement.
4.4.2 Théorème de Cayley - Hamilton

Soient C un tenseur de dimension n et d'invariants Ik , k = 1, · · · , n, et
P (λ) son polynôme caractéristique donné par
P (λ) = det(C − λI) = (−1)n λn + (−1)n−1 I1 λn−1 + · · · + In (4.13)
Alors, le tenseur C satisfait à son polynôme caractéristique dans le sens

où,
(−1)n C n + (−1)n−1 I1 C n−1 + · · · + In I = 0 (4.14)
De plus, les valeurs propres de C sont les racines de (4.13).
4.4.3 Théorèmes de représentation

Dans la suite C désigne un tenseur symétrique tel que det C > 0. On
peut alors démontrer les théorèmes suivants :
1. Fonction scalaire d'une variable tensorielle : toute fonction sca-

laire d'un seul tenseur symétrique C est isotrope si et seulement si
elle peut s'exprimer de manière unique en fonction des n invariants
principaux (ou alternativement des moments) de C :
f (C) = f (I1 , · · · , In ) = f (I 1 , · · · , I n )
48
2. Fonction tensorielle d'une variable tensorielle : toute fonction
tensorielle symétrique d'un seul tenseur symétrique est isotrope si et
seulement si elle peut s'écrire
H(C) = φ0 I + φ1 C + φ2 C 2 + · · · φn−1 C n−1
où les φi sont des fonctions scalaires des invariants de C : φi ≡

φi (I1 , · · · , In ).
3. Fonction vectorielle d'un tenseur et d'un vecteur : toute fonc-
tion vectorielle g(b, C) d'un vecteur b et d'un tenseur symétrique C,
est isotrope si et seulement si elle peut s'exprimer de manière unique
sous la forme
g(b, C) = φ0 I + φ1 C + · · · + φn−1 C n−1 b

où les φi sont des fonctions scalaires, non seulement des invariants de

C, mais aussi de b · b, b · Cb, ..., b · C n−1 b.
Ces théorèmes peuvent se généraliser (pour n ≤ 3) à des fonctions de
plusieurs tenseurs. Ainsi, pour qu'une fonction scalaire dépendant de deux
tenseurs symétriques soit isotrope, il faut et il sut qu'elle s'écrive de ma-
nière unique en fonction des invariants de ces deux tenseurs, mais aussi des
invariants de leurs produits :
f (A, B) = f (tr A, · · · , tr B, · · · , tr AB, tr AB 2 , tr A2 B, tr A2 B 2 )
4.4.4 Cas de R3
Toute matrice C, carrée symétrique réelle de dimension 3, possède trois
valeurs propres réelles, λ1 , λ2 et λ3 , racines du polynôme caractéristique
P (θ) = det(C − θI) = −θ3 + I1 θ2 − I2 θ + I3
Les invariants de C s'écrivent en fonction de ses valeurs propres :
I1 = tr C = C : I = λ1 + λ2 + λ3 (4.15)
1
I2 = det C(tr C −1 ) = (trC)2 − tr(C 2 ) = λ1 λ2 + λ2 λ3 + λ3 λ1

2
I3 = det C = λ1 λ2 λ3
On peut également calculer (par calcul indiciel), la variation des invariants

par rapport à C. On obtient ainsi :
dI1
=I (4.16)
dC
49
dI2
= I2 C −1 − I3 C −2 (4.17)
dC
dI3
= I3 C −1 (4.18)
dC
Le théorème de Cayley - Hamilton s'écrit
C 3 − I1 C 2 + I2 C − I3 I = 0
d'où l'on déduit
1 2
C −1 = (C − I1 C + I2 I)
I3
et
I2

−2 1 I2 2 I1 I2
C = C + (1 − )C + ( 1 − I2 )I
I3 I3 I3 I3
4.5 Thermodynamique des milieux continus

La température intervient à divers titres dans le comportement thermo-
mécanique des milieux continus : introduction de contraintes d'origine ther-
mique, dépendance des paramètres matériau à la température.
Le bilan d'énergie introduit 5 nouvelles variables : la température T, le
ux de chaleur q et l'énergie interne e. Nous devrons donc établir des lois
de comportement thermiques an de clore le système. Ces lois dépendent du
matériau. Au nal, le problème thermomécanique est constitué de 5 équations
de conservation, 10 équations de comportement, et 15 inconnues à trouver.
Suivant la démarche précédente, on pose :
σ(x, t) = F (F (X, s), T (X, s), G(X, s))

0≤s≤t
q(x, t) = G (F (X, s), T (X, s), G(X, s))

0≤s≤t
e(x, t) = H (F (X, s), T (X, s), G(X, s))

0≤s≤t
avec
∂T
G(X, t) = (X, t)
∂X
4.5.1 Second principe de la thermodynamique des milieux conti-

nus
L'entropie S d'un ensemble matériel ω est dénie à partir d'une densité
massique d'entropie s(x, t)
Z
S(ωt ) = ρsdv
ωt
50
La variation d'entropie est due pour partie à des apports surfaciques et
volumiques Z Z
Φ(ωt ) = − φ · nda + ρrs dv
∂ωt ωt
où φ est le vecteur ux d'entropie, et rs une densité de production locale
d'entropie. La production interne d'entropie est alors la diérence
Σ(ωt ) = Ṡ(ωt ) − Φ(ωt )
Le second principe de la thermodynamique stipule que cette production

est positive ou nulle :
Σ(ωt ) ≥ 0
soit, de manière explicite :
Z Z Z
d
ρsdv + φ · nda − ρrs dv ≥ 0
dt ωt ∂ωt ωt
Par application du théorème de la divergence pour convertir l'intégrale surfa-

cique en intégrale volumique, et par utilisation de l'équation de conservation
de la masse (1.28) pour permuter les opérateurs de dérivation en temps et
d'intégration en espace, on aboutit à la forme locale de l'inégalité d'entropie :
ρṡ + ∇ · φ − ρrs ≥ 0 (4.19)
Le ux d'entropie est généralement relié au ux de chaleur et à la température

par :
q
φ=
T
Tandis que la production locale d'entropie vaut
r
rs =
T
où r est la production de chaleur. Ce qui mène à :
q r
ρṡ + ∇ · −ρ ≥0 (4.20)
T T
On introduit à présent l'énergie libre Ψ de densité massique ψ , à titre d'exemple
(mécanique des solides). On dénit
ψ = e − Ts
et Z
Ψ(ωt ) = ρψdv
ωt
51
1
On peut alors exprimer la variation d'entropie comme ṡ = (ė − Ṫ s − ψ̇).
T
On calcule alors la dissipation volumique d :
q
d = ρT ṡ + T ∇ · − ρr
T
q
= ρ(ė − Ṫ s − ψ̇) + T ∇ · − ρr
T
q
= σ : D − ∇ · q − ρ(Ṫ s + ψ̇) + T ∇ ·
T
où l'on a utilisé le bilan d'énergie (3.4) ρė = σ : D − ∇ · q + ρr.
En remarquant que
q ∇T
T∇ · =∇·q−q· ,
T T
la dissipation se réduit nalement à :
∇T
d = σ : D − ρ(Ṫ s + ψ̇) − q ·
T
On obtient ainsi l'inégalité de Clausius-Duhem qui constitue une
condition d'admissibilité thermodynamique des lois de comportement ther-
momécanique :
∇T
d = σ : D − ρ(Ṫ s + ψ̇) − q · ≥0 (4.21)
T
On distingue deux contributions à la dissipation :
La dissipation intrinsèque dint = σ : D − ρ(Ṫ s + ψ̇)

∇T
La dissipation thermique dth = −q ·
T
De sorte que
d = dint + dth ≥ 0
Comme généralement les lois donnant la contrainte, l'énergie libre et l'entro-
pie ne dépendent pas explicitement de la température, la condition d'admis-
sibilité devient :
dint ≥ 0 et dth ≥ 0
4.5.2 Loi de Fourier

Remarquons que l'inégalité dth ≥ 0 signie que le ux de chaleur est
opposé au gradient de température. La chaleur va donc toujours du chaud
vers le froid.
Si l'on pose la loi de Fourier q = −k∇T , alors on doit avoir la conductivité
thermique k ≥ 0.
52
5 Hyperélasticité
Un milieu (matériellement simple) est dit élastique si l'état de contrainte
actuel est entièrement déterminé par le gradient de la transformation à l'ins-
tant actuel et non par son histoire passée :
σ(x, t) = FD0 (F (X, t)) (5.1)
En particulier, le comportement ne dépend pas de la vitesse de sollicitation

D . La forme de la loi dépend a priori de la conguration de référence choisie
D0 .
Considérons un jeu de variables pour représenter les matériaux (solides) élas-
1 T
tiques : {E(X, t), T (X, t), G(X, t)} où E= (F F − I) est la mesure de
2
∂T
déformation de Green-Lagrange, et
∂X
G=
est le gradient de température
par rapport à la conguration de référence.
L'énergie libre est aussi fonction de ces variables :
ψ0 (E(X, t), T (X, t), G(X, t))

On exprime alors le second principe de la thermodynamique sous forme locale
en formulation lagrangienne. Pour un domaine matériel ωt , on a :
Z Z Z
d q r
ρs(x, t)dv + · nda − ρ dv ≥ 0
dt ωt ∂ωt T ωt T
En utilisant le changement de variable x = f (X, t), on se ramène à des

intégrales sur le domaine de référence :
Z Z Z
d Q r0
ρ0 s0 (X, t)dV + · N dA − ρ0 dV ≥ 0
dt ω0 ∂ω0 T ω0 T
avec s0 (X, t) = s(x, t) et r0 (X, t) = r(x, t), ω0 = f −1 (ωt ), dv = JdV . Le
ux de chaleur Q est donné par :
q · nda = Q · N dA ⇒ Q = JF −1 q
D'où l'inégalité locale :
Q r0
ρ0 ṡ0 + ∇X · − ρ0 ≥ 0
T T
En suivant la même démarche que précédemment, on obtient l'inégalité de
Clausius-Duhem en description lagrangienne :
G
D = S : Ė − ρ0 (Ṫ s0 + ψ˙0 ) − Q · ≥0 (5.2)
T
53
La dérivée totale en temps de l'énergie libre est donnée par :
∂ψ0 ∂ψ0 ∂ψ0

ψ̇0 = : Ė + Ṫ + · Ġ
∂E ∂T ∂G
En substituant ce développement dans (5.2) on trouve :
∂ψ0 ∂ψ0 ∂ψ0 G

D = (S − ρ0 ) : Ė − ρ0 ( + s0 )Ṫ − ρ0 · Ġ − Q · ≥0
∂E ∂T ∂G T
Ceci étant valable pour toute évolution {Ė, Ġ, Ṫ }, les quantités qui mul-
tiplient ces dérivées doivent s'annuler an d'avoir une dissipation D non
négative. On en déduit alors les relations constitutives suivantes :
∂ψ0 ∂ψ0 ∂ψ0

S = ρ0 ; s0 = − ; =0 (5.3)
∂E ∂T ∂G
Ces relations sont des lois de comportement dites hyperélastiques pour les
contraintes et l'entropie. Elles montrent que l'énergie libre est un potentiel
d'élasticité dont découlent les relations contrainte - déformation et entropie
- température. L'énergie libre ne peut dépendre du gradient de température
G.
5.1 Réversibilité des lois hyperélastiques

Une loi hyperélastique est une loi d'élasticité au sens (5.1), particulière au
sens qu'il existe un potentiel élastique tel que nous l'avons déni. Cependant,
toutes les lois ne sont pas hyperélastiques : il existe des lois élastiques ne
dérivant pas d'un potentiel (voir exercice 5.1). La caractéristique principale
d'une loi hyperélastique est qu'avec une telle loi, le travail élastique fourni
ne dépend pas du trajet de déformation suivi. Ceci n'est pas le cas en toute
généralité pour les lois élastiques (5.1).
Considérons un trajet de déformation imposé allant d'un état de défor-
mation initial E = EA à un état nal E = EB à température constante. On
calcule la densité de travail mécanique fourni en un point matériel lors de
cette transformation en intégrant le long du chemin la densité de puissance :
Z tB Z tB
∂ψ0
W (E A → E B ) = S : Ė dt = ρ0
: Ė dt
tA tA ∂E
= ρ0 ψ0 (E B , T ) − ρ0 ψ0 (E A , T ) (5.4)
Ainsi, le travail W ne dépend que des valeurs initiales et nales des variables
d'état et non du chemin suivi entre ces deux extrémités. Si ce chemin est un
54
cycle, i.e.
si E = E , alors le travail fourni sur l'ensemble du cycle est nul :
A B
l'énergie élastique stockée au cours d'une partie du trajet, est intégralement
restituée durant le trajet retour, même si celui-ci s'eectue par un autre
chemin. C'est le caractère réversible du comportement hyperélastique.
5.2 Hyperélasticité isotrope

5.2.1 Description matérielle (S, C)
D'après (5.3), nous avons
ρ0 ψ0 = W (E) = W (C)
Un corps hyperélastique est dit isotrope si son énergie libre est une fonc-
tion isotrope de ses arguments. Le théorème de représentation d'une fonction
scalaire isotrope dépendant d'une variable tensorielle, vu à la section 4.4.3
permet d'écrire W comme une fonction de trois invariants de C (ici préféré
à E) :
W (C) = W (I1 , I2 , I3 )
On a donc, pour un tel matériau :
!
∂W ∂W ∂I1 ∂W ∂I2 ∂W ∂I3
S=2 =2 + +
∂C ∂I1 ∂C ∂I2 ∂C ∂I3 ∂C
où les invariants de C choisis sont les invariants principaux
I1 = trC,
1
(trC)2 − trC 2

I2 =
2
et
I3 = detC = J 2
La relation reliant les contraintes aux déformations est donc obtenue en
calculant les dérivées des invariants par rapport à C .
∂I1
Calcul de . La dérivée d'un scalaire par un tenseur est un tenseur.
∂C
Dans le cas présent, l'entrée (i, j) de ce tenseur est donnée par :
∂I1 ∂Ckk
= = δij
∂Cij ∂Cij
D'où
∂I1
=I (5.5)
∂C
55
∂I2 ∂tr2 C
Calcul de
∂C
. Remarquons d'abord que
∂C
= 2I1 I . Ensuite,
∂trC 2 ∂Ckl Clk

= = 2Cji = 2Cij
∂Cij ∂Cij
D'où
∂I2
= I1 I − C (5.6)
∂C
∂I3
Calcul de . Le théorème de Cayley-Hamilton (4.14) s'écrit
∂C
C 3 − I1 C 2 + I2 C − I3 I = 0
D'où l'on déduit par contraction à gauche et à droite avec I
3I3 = trC 3 − I1 trC 2 + I2 trC
Ainsi,
∂I3 ∂trC 3 2
∂trC 2
3 = − trC I − I1 + (I1 I − C)trC + I2 I
∂C ∂C ∂C
= 3C 2 − trC 2 I − 2I1 C + (trC)2 I − CI1 + I2 I
= 3C −1 (C 3 − I1 C 2 + I2 C) = 3C −1 I3 (5.7)
En dénitive,
∂I3
= I3 C −1 (5.8)
∂C
Les relations que nous venons d'établir conduisent à :
S = 2 (W1 + I1 W2 )I − W2 C + I3 W3 C −1

(5.9)
où l'on a posé
∂Wi
Wi =
∂Ii
5.2.2 Descriptions matérielles alternatives

Il est possible de choisir d'autres invariants de C. Par exemple, en consi-
dérant
W = W (IC , IIC , IIIC ),
56
avec
IC = trC = C : I (5.10)
1
IIC = trC 2 /2 = C : C (5.11)
2
3
IIIC = trC /3 (5.12)
on obtient
∂W ∂W ∂W
S=2 I+ C+ C2 (5.13)
∂IC ∂IIC ∂IIIC
De plus, si l'on choisit maintenant W comme étant une fonction des in-
variant du tenseur de Green-Lagrange E,
W = W (IE , IIE , IIIE ),
alors
∂W ∂W ∂W
S= I+ E+ E2 (5.14)
∂IE ∂IIE ∂IIIE
En comparant cette expression à la loi de comportement générale d'un corps
élastique isotrope (voir théorème de représentation),
S = α0 I + α1 E + α2 E 2
où les αi sont fonctions des invariants de E, nous constatons qu'un solide

élastique est hyperélastique si il existe un potentiel W tel que :
∂W ∂W ∂W
α0 = , α1 = , α2 =
∂IE ∂IIE ∂IIIE
Un tel potentiel existe si les fonctions αi vérient les conditions :
∂α0 ∂α1 ∂α1 ∂α2 ∂α2 ∂α0

= , = , =
∂IIE ∂IE ∂IIIE ∂IIE ∂IE ∂IIIE
5.2.3 Description spatiale

La seule contrainte ayant une signication physique est la contrainte de
Cauchy σ, dénie dans la conguration courante (ou spatiale). Rappelons
que l'on passe de S à σ par la formule
σ = J −1 F SF T
57
Ainsi, via (5.13), et en utilisant le tenseur de Cauchy - Green gauche b =
FFT,

∂W −1 ∂W ∂W
σ = 2J FFT + F CF T + F C 2F T
∂I ∂IIC ∂IIIC
C
∂W ∂W 2 ∂W 3
= 2J −1 b+ b + b
∂IC ∂IIC ∂IIIC
Cette expression fait toujours intervenir les dérivées partielles de W par rap-
port aux invariants du tenseur lagrangien C. Cependant, on montre que les
invariants de b sont identiques à ceux de C:
Ib = trb = tr(F F T ) = tr(F T F ) = trC = IC
2IIb = trb2 = tr(F F T F F T ) = tr(F T F F T F ) = trC 2 = 2IIC

3IIIb = trb3 = tr(F F T F F T F F T ) = tr(F T F F T F F T F ) = trC 3 = 3IIIC
où l'on se sert à chaque fait du fait que pour deux tenseurs quelconques
d'ordre 2 A et B, alors tr(AB) = tr(BA). Ainsi,

∂w ∂w ∂w 2
σ = 2 I+ b+ b b
∂Ib ∂IIb ∂IIIb
∂w
= 2 b
∂b
Ainsi, dans le cas isotrope, il est possible de dénir un potentiel élas-

tique w(b, T ) = ρψ(b, T ) = ρ0 ψ0 (E, T ) = W (E, T ) dont dérive le tenseur
des contraintes. Remarquons cependant qu'en toute généralité, il n'est pas
possible de travailler dans la conguration courante avec un couple de va-
riables contraintes de Cauchy - déformations, puisque aucune déformation
n'est conjuguée avec σ. On ne peut pas, alors, exprimer directement un po-
tentiel w(b).
Exercice 5.1 Exemple d'une loi élastique non hyperélastique. On

considère la loi élastique isotrope suivante :
α
S= (trE 2 )I
2
1. Montrez que cette loi n'est pas hyperélastique.
58
2. On considère deux chemins de déformation paramétrés par le réel λ,
0≤λ≤1:
[E 1 (λ)] = diag(λ, λ, 0) et [E 2 (λ)] = diag(λ, λ2 , 0)
exprimés dans la base canonique de R3 .

Ces deux chemins partent du même état de déformation nulle E A =
E 1 (λ = 0) = E 2 (λ = 0) = 0, pour arriver au même état nal E B =
E 1 (λ = 1) = E 2 (λ = 1) représenté par diag(1, 1, 0), mais en passant
par des états transitoires diérents.
Calculez le travail élastique W1 eectué le long du chemin 1, et W2
eectué le long du chemin 2. Conclure.
5.3 Matériaux Néo-hookéens

Jusqu'à présent nous n'avons pas considéré de loi hyperélastique isotrope
particulière. C'est ce que nous faisons ici en considérant une loi simple, celle
des matériaux néo-hookéens compressibles. Le potentiel élastique est donné
par :
µ λ
W = (IC − 3) − µ log J + (log J)2 (5.15)
2 2
où λ et µ sont des paramètres du matériau, IC = trC et J = detF . L'énergie
stockée s'annule bien en l'absence de déformation, i.e. pour F = I . Dans
2
l'expression de cette loi, un choix naturel d'invariants est donc I1 et I3 = J .
D'où, l'expression (5.9) donne :
S = 2 W1 I + J 2 W3 C −1

(5.16)
Le modèle néo-hookéen est simple, et rend compte du comportement d'un

élastomère pour des déformations allant jusqu'à 100%.
Exercice 5.2 Montrez que

S = µ(I − C −1 ) + λ(log J)C −1
et que
µ λ
σ= (b − I) + (log J)I
J J
Que deviennent S et σ pour un matériau incompressible ?
59
Figure 5.1: Contrainte réelle σ11 lors d'essais de compression / traction
uniaxiale. Source : Wikipédia. Voir aussi : L R G Treloar (1944). Stress-
strain data for vulcanised rubber under various types of deformation. Trans.
Faraday Soc. 40 : 5970.
Dans l'approximation des petites déformations, on dénit le tenseur des

1 T
déformations ε(u) = (∇u + ∇ u). On a alors J ≈ 1 + trε ≈ 1, log J ≈ trε,
2
et b ≈ 2ε + I . D'où, sous cette hypothèse :
σ ≈ 2µε + λtrεI
ce qui permet d'identier les coecients µ et λ comme étant les coecients

de Lamé en petites déformations.
Exercice 5.3 Dilatation pure - On considère une dilatation pure, caractéri-

sée par xi = ΛXi . Montrez que le tenseur de Cauchy représente alors un état
de contrainte hydrostatique (tenseur sphérique) avec une pression donnée par
µ λ
p = − (J 2/3 − 1) − log J
J J
Exercice 5.4 Cisaillement simple - On considère la transformation (cisaille-
ment simple) de la gure 5.2 pour un matériau néo-hookéen. Calculez le ten-
60
Figure 5.2: Cisaillement simple
seur des contraintes de Cauchy. Montrez que la pression est non nulle même
si la transformation est isochore (J = 1).
Exercice 5.5 Traction biaxiale sur un élastomère. On considère une plaque
d'un élastomère, de 200 mm de côté dans le plan (X1 , X2 ) et d'épaisseur 10
mm, comme indiqué sur la gure 5.3. Le matériau est supposé incompressible
sous les sollicitations envisagées, et présente un module néo-Hooke de µ =
0, 75M P a. Quelles forces F1 et F2 doit-on appliquer dans les directions X1
et X2 respectivement pour étirer de manière homogène la plaque, de 400 mm
selon X1 et 300 mm selon X2 ? (Indication : quel tenseur des contraintes
choisir ?).
Figure 5.3: Traction biaxiale.
Références
[COIRIER] Coirier, J. Mécanique des milieux continus - 4e édition : Cours
et exercices corrigés. Dunod, 2013.
[FOREST] Forest, S. et al. Mécanique des Milieux continus. École des Mines
de Paris. Cours disponible sur internet.
61
[FORTIN] Fortin A. et Garon A. Les éléments nis : de la théorie à la pra-
tique. Cours disponible sur internet.
[GERMAIN] Germain, P. Cours de Mécanique des Milieux Continus. Mas-

son, 1973.
[LE MAÎTRE] Le Maître, J. et Chaboche, J.-L. Mécanique des matériaux

solides. Dunod, 2009.
[PONTHOT] Ponthot J.-P. Cours de Mécanique des Milieux Continus, Uni-

versité de Liège.
62

MMC Loi

Transféré par

Droits d'auteur :

Formats disponibles

MMC Loi

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

MMC Loi

Transféré par

Droits d'auteur :

Formats disponibles

Majeure Mécanique

UP1 - Mécanique des Matériaux

Thermodynamique des milieux continus et lois de

2 Eorts intérieurs - Théorème de l'énergie cinétique 26

4 Formulation des lois de comportement 32

(voir gure 1.1).

Figure 1.1: Description du mouvement

Si le mouvement du corps est décrit en utilisant les coordonnées matérielles

x = (1 + t)X, avec 0≤X≤2

De plus, on considère une certaine distribution de température sur la barre,

T (X, t) = Xt2 , avec 0≤X≤2

La loi du mouvement permet de passer à une description eulérienne, avec un

T (x, t) = xt2 /(1 + t)

Figure 1.2: Traction uniaxiale : description eulérienne et lagrangienne

Intéressons nous maintenant au calcul de la variation de la température

Évidemment, nous avons dans ce cas, puisque X ne dépend pas de t,

T (x, t + ∆t) − T (x, t)

T (x + ∆x, t + ∆t) − T (x, t)

dT T (x + ∆x, t + ∆t) − T (x, t)

et est la dérivée totale en temps de la fonction T (x(t), t) où x(t) est la

1.2 Le gradient de déformation

Figure 1.3: Transformation d'un vecteur matériel élémentaire de la con-

Considérons comme sur la gure 1.3 deux points matériels de la congu-

Ainsi, on en conclut que

∂fi ∂fi ∂fi

∂f1 ∂f1 ∂f1

∂X1 ∂X2 ∂X3

1.2.1 Produit mixte

1.2.2 Transport d'un élément de volume

en (dx1 , dx2 , dx3 ) et le volume initial se transforme en un volume dv donné

dv = [dx1 , dx2 , dx3 ]

dv = JdV avec J = détF > 0 (1.7)

Une transformation est dite isochore au point X si le volume est préservé

1.2.3 Transport d'un élément de surface orienté

sur un vecteur v quelconque :

(nda) · v = (dx1 ∧ dx2 ) · v = [dx1 , dx2 , v]

d'où, par identication :

1.2.4 Décomposition polaire du gradient de la transformation

avec U le tenseur de déformations pures à droite qui se trouve dans la con-

Les tenseurs U et V peuvent également s'écrire en fonction du tenseur F :

avec C = FTF, ce qui en utilisant la relation (1.11) revient à C = U 2.

dL2 = dX T dX = dxT F −T F −1 dx = dxT b−1 dx (1.13)

Figure 1.6: Décomposition polaire du gradient de la transformation. Les

Exercice 1.1 On étudie la transformation faisant passer de la conguration

Exercice 1.2 On étudie la transformation illustrée sur la gure 1.8 qui

Figure 1.8: Cinématique du glissement simple pour un corps carré muni

1. Donner l'expression de F , vérier que la transformation est isochore.

5. À quel type de rotation R correspond ? Donnez l'expression de tan θ,

Dans la conguration courante, la longueur de la bre devient

|dx|2 = (F dX)T (F dX) = |dX|2 M T CM

On dénit alors le rapport d'allongement dans la direction M par :

Si M est choisi comme le ième vecteur de la base canonique {Ei }, i =

ce qui donne une interprétation simple des composantes diagonales du ten-

Angle de glissement - Le cosinus de l'angle formé par deux vecteurs

On appelle γ = Θ − θ l'angle de glissement des directions M1 et M2 dans

Nous avons vu que grâce à la décomposition polaire il est possible de

1. Il est symétrique et sans dimension physique ;

2. Il est nul pour un mouvement de corps rigide et en F =I;

2 Eorts intérieurs - Théorème de l'énergie cinétique 26

(voir gure 1.1).

Figure 1.3: Transformation d'un vecteur matériel élémentaire de la con-

Considérons comme sur la gure 1.3 deux points matériels de la congu-

d'où, par identication :

avec U le tenseur de déformations pures à droite qui se trouve dans la con-

Exercice 1.1 On étudie la transformation faisant passer de la conguration

Exercice 1.2 On étudie la transformation illustrée sur la gure 1.8 qui

1. Donner l'expression de F , vérier que la transformation est isochore.

Dans la conguration courante, la longueur de la bre devient

On dénit alors le rapport d'allongement dans la direction M par :

Hypothèse des petites déformations - Enn, remarquons que le ten-

1.5 Lois de conservation dans la conguration courante

De plus, rappelons le théorème de la divergence. Soit q un vecteur déni

Ceci signie que df se déduirait de dF par transport convectif. Dans la

2 Eorts intérieurs - Théorème de l'énergie ci-

Ceci signie que la densité des eorts intérieurs rapportée à l'élément de

L'énergie interne E du système, représentée par une densité massique

L'existence du ux de chaleur q est une conséquence directe de la dépen-

où η est le coecient de viscosité de l'amortisseur, exprimé en P a.s.m

Plasticité. Le comportement rigide parfaitement plastique uniaxial

et le glissement (écoulement) du patin s'écrit

A titre d'exemple, on considère le modèle viscoélastique (pour les uides)

4.2 Le principe d'indiérence matérielle ou d'objecti-

4.2.1 Dénitions : repère et référentiel