3.2.

3 Opérations sur le tenseur des contraintes

Déviateur des contraintes

Le déviateur des contraintes est défini comme:

1
S= − tr( )I (127)
3
Son expression indicielle est: Sij = σij − σ3kk δij
On note ainsi, quelque soit le tenseur , σH = 31 tr( ) la part hydrostatique du
tenseur, ou simplement la contrainte hydrostatique. Ce terme renvoie à la notion
de fluide. Dans un milieu fluide, la contrainte hydrostatique correspond à −p, p
étant la pression. σH est invariant par changement de repère. Quand σH = 0, on
parlera d’état de contrainte purement déviatorique, état de contrainte caractéris-
tique du cisaillement.

La trace du déviateur est par ailleurs toujours nulle.

Réciprocité des contraintes - théorème de Cauchy

soit T# (M,#n) un vecteur contrainte tel que T# (M,#n) = (M).#n,

∀(#n,#n! ) deux vecteurs unitaires, on montre que t#n! .T# (M,#n) = t#n.T# (M,#n! )

Démonstration: n!i .(σij .nj ) = (n!i .σij ).nj = nj .(σji .n!i ) - lié à la symétrie du
tenseur des contraintes.

Conditions limites - surface libre

On a défini le vecteur contrainte à partir de 2 volumes imbriqués l’un dans

l’autre. Quand on fait tendre le volume Φ vers Ω, les surfaces se confondent. Il y
a donc correspondance entre le vecteur contrainte T# (M,#n) et la force surfacique
appliquée au point M considéré de la surface. Si l’effort appliqué est un torseur
équivalent, il y a correspondance entre la résultante du torseur et l’intégrale de
surface du vecteur contrainte et entre le moment du torseur et l’intégrale de sur-
face du moment du vecteur contrainte par rapport au même point.
Enfin si la surface n’est soumise à aucun effort, elle est dite "libre", le vecteur
contrainte est donc nul au point M considéré.

41
4 CHAPITRE 4 - LOIS DE CONSERVATION
On présente dans ce chapitre les lois de conversation des grandeurs physiques
utiles au physicien et au mécanicien.
On établit ces équations en effectuant un bilan des différentes grandeurs sur un
volume Ω fixe (description Eulerienne) limité par une surface ∂Ω. L’intégrale de
volume est l’outil essentiel utilisé ici. L’intégrale de surface est utilisé plus parti-
culièrement en électromagnétisme.

4.1 Dérivation par rapport au temps des intégrales spatiales

En mécanique, lorsqu’on calcule une intégrale (de surface ou de volume), on
peut utiliser aussi bien les coordonnées d’Euler et de Lagrange. Cela ne pose pas
de problème particulier au niveau de l’intégration mais, en revanche, le compor-
tement est différent vis-à-vis de la dérivation temporelle.

Dérivation d’une intégrale de volume

Dans le cas des représentations d’Euler, on s’intéresse à des volumes et des

surfaces fixes. On souhaite pouvoir permuter dérivation et intégration pour une
−−→
quantité A(OM,t).
−−→
d −−→ dA(OM,t)
!!! !!! !!
A(OM,t)dv = dv + autre terme? ds (128)
dt Ω Ω dt ∂Ω

Le terme supplémentaire est un terme de flux (convection à travers la sur-

face). Sa définition ne peut s’intuiter facilement !
En revanche, si on utilise la représentation de Lagrange, le support d’intégra-
tion est mobile. Ainsi, à l’aide d’un changement de variable, on exprime l’inté-
grale de volume par rapport au volume initial Ω0 :
−−→ −−→
!!! !!!
A(OM,t)dv = A(φ̄(OM 0 ,t),t)Jdv0 (129)
Ω(t) Ω0
−−→ −−→
avec OM = φ̄(OM 0 ,t). J est le jacobien du changement de variable, ce qui
est classique pour n’importe quel changement de variable. Le Jacobien est la va-
leur absolue du déterminant de la matrice jacobienne J.

42
 & Construction du Jacobien

La matrice jacobienne permet un passage entre variables d’une configuration

à une autre. On sait que le changement de configuration s’effectue selon dx # =
# 0 , F correspondant au gradient de la fonction vectorielle φ̄, c’est également la
Fdx
matrice Jacobienne.
dxi
Fij = 0
dxj
Exprimons alors le passage de V0 à V :
volume: V = (dx∧# dy).
# dz # = dwi dzi = 'ijk dxj dyk dzi = 'ijk Fjm dx0 Fkn dy 0 Fil dz 0
m n l
0 0 0
# # # 0
volume: V0 = (dx ∧ dy ).dz = 'lmn dxm dyn dzl 0 0

On utilise la propriété: 'lmn .'lmn = 6 pour écrire: 'lmn V0 = 6.dx0m dyn0 dzl0 , soit
finalement:
1
V = ( 'ijk 'lmn Fjm Fkn Fil ).V0 (130)
6
On vérifie que la quantité 61 'ijk 'lmn Fjm Fkn Fil est le déterminant de F, et que
cette quantité est toujours positive (si produit mixte dans le sens direct). Puisque
F correspond à la matrice Jacobienne du changement de variable, le jacobien vaut
J = det(F). La dérivée par rapport au temps devient donc:

−−→
d −−→ dA(φ̄(OM 0 ,t),t) −−→ dJ
!!! !!!
A(OM,t)dv = J + A(φ̄(OM 0 ,t),t) dv0
dt Ω(t) Ω0 dt dt
(131)
La jacobien est construit à partir des des composantes de F, matrice exprimée
dans un cadre lagrangien. La dérivée par rapport au temps de J est donc une
dérivée ronde:
dJ ∂J 1 ∂Fil
= = 'ijk 'lmn Fjm Fkn (132)
dt ∂t 2 ∂t
On cherche ensuite à reconnaître le jacobien dans la dérivée par rapport au temps:
∂Fil ∂ ∂xi ∂ ∂xi ∂Vi ∂Vi ∂xp
= 0
= 0 = 0 = = Vi,p Fpl (133)
∂t ∂t ∂xl ∂xl ∂t ∂xl ∂xp ∂x0l

Soit:
dJ 1
= 'ijk 'lmn Fjm Fkn Fpl Vi,p (134)
dt 2

43
 det(F) = 16 'pjk 'lmn Fjm Fkn Fpl
'pjk det(F) = 'lmn Fjm Fkn Fpl

Soit:
dJ 1
= 'ijk 'pjk det(F)Vi,p (135)
dt 2
On utilise enfin la propriété: 'ijk .'pjk = 2δip pour écrire:

dJ
= δip det(F)Vi,p = det(F)Vi,i = det(F)div(E V# ) = J.div(E V
#) (136)
dt
Si on reprend maintenant l’expression de la dérivée de l’intégrale de volume,

d −−→ ∂A −−→
!!! !!!
A(OM,t)dv = ( + grad(A).E V# +A.div(E V
# ))Jdv0 (137)
dt Ω(t) Ω0 ∂t
−−→ # + A.div(E V# ) =
Cette expression se simplifie en observant que: grad(A).E V
div(A.E V# )
On constate alors que:
∂A ∂A
!!! !!!
E# # ))dv
( + div(A. V ))Jdv0 = ( + div(A.E V (138)
Ω0 ∂t Ω ∂t

L’intégrale de volume se décompose en 2 contributions puisque le théorème

de Gauss-Ostrogradski (ou théorème de la divergence) indique, #n définissant la
normale locale à la surface, que:
!!! !!
div(#u)dv = #u.#nds (139)
Ω ∂Ω

On a donc l’expression finale suivante:

d −−→ ∂A
!!! !!! !!
A(OM,t)dv = dv + # .#nds
A.E V (140)
dt Ω(t) Ω ∂t ∂Ω

On remarque que les variations sont dues d’une part aux variations de A sur
tout le volume et d’autre part à un effet de frontière correspondant au déplacement
de matière à travers la surface ∂Ω et n’impliquant que les valeurs de A sur ∂Ω.

44
 Dérivation d’une intégrale de surface

L’intégrale de surface est utilisée en particulier lorsqu’il s’agit de calculer des

circulations via la formule de Stockes. Sa dérivation par rapport au temps permet
de traiter le cas des circulations dans un milieu mobile (expressions généralement
non traitées en électromagnétisme).

On s’intéresse en particulier à la dérivation par rapport au temps de l’intégrale

suivante:
d # − −→
!!
A( OM,t).#nds (141)
dt ∂Ω
# −
A(
−→
OM,t) est une quantité vectorielle, #n la normale sortante à la surface ∂Ω.
A l’aide d’un changement de variable, on exprime l’intégrale de surface par
rapport à la surface initiale ∂Ω0 de normale sortante #n0 , #n étant le résultat de la
transformation de #n0 par F.
−−→ # φ̄(−−→
!! !!
#
A(OM,t).#nds = A( OM 0 ,t),t).F.#n0 Jds0 (142)
∂Ω(t) ∂Ω0
−−→ −−→
avec OM = φ̄(OM 0 ,t). J est le jacobien du changement de variable.

La dérivée par rapport au temps devient donc:

d #
dA # dF .#n0 .J + A.F.#
# n0 dJ ds0
!! !!
# nds =
A.# .F.#n0 J + A. (143)
dt ∂Ω(t) ∂Ω0 dt dt dt

On montre facilement que le jacobien déterminé pour l’intégrale de volume

reste valable.

d #
∂A # dF .#n0 .J+A.F#
!! !!
# nds =
A.# ( # V
+grad(A). # ).F#n0 J+A. # n0 div(V# )Jds0
dt ∂Ω(t) ∂Ω0 ∂t dt
(144)
La dérivée par rapport au temps de F a déjà été établie.
dF # )F
= grad(V
dt
45
 On obtient finalement:
d
!! !!
∂A#
#
A.#nds = ( # V
+ grad(A). # + A.grad(
# # ) + A.div(
V # # )).F#n0 Jds0
V
dt ∂Ω(t) ∂Ω0 ∂t
(145)
d
!! !!
∂A#
# nds =
A.# ( # V
+grad(A). # +A.grad(
# # )+A.div(
V # # )).#nds (146)
V
dt ∂Ω(t) ∂Ω ∂t
On montre également que cette relation est équivalente à:

d
!! !! # −
∂A → # #
# nds =
A.# ( + rot(A ∧V)+V # .div(A)).#
# nds (147)
dt ∂Ω(t) ∂Ω ∂t

4.2 Conservation des quantités mécaniques

4.2.1 Conservation de la masse
On regarde l’évolution temporelle de la masse m contenue dans le volume Ω
fixe (eulerien). Par définition:
!!!
m(t) = ρdv (148)
Ω
la variation de masse en fonction du temps vaut donc.
dm(t) d ∂ρ
!!! !!!
= ρdv = # )dv
+ div(ρ.E V (149)
dt dt Ω Ω ∂t
On retrouve deux contributions, l’une liée à "compressibilité" du milieu (varia-
tion de masse volumique), l’autre liée au déplacement des particules c’est-à-dire
au flux à travers la surface. L’hypothèse de conservation de la masse implique la
nullité de dm(t)
dt
, ce qui implique la nullité de l’intégrale mais également de l’inté-
grande, soit:
∂ρ #)=0
+ div(ρ.E V (150)
∂t
Elle peut se mettre sous la forme développée suivante :
∂ρ −−→ # + ρ.div(E V# ) = 0
+ grad(ρ).E V (151)
∂t
ou encore
Dρ
+ ρ.div(E V# ) = 0 (152)
Dt
46
- Si on traite d’un fluide incompressible, ρ est constant, ce qui se traduit par:
#)=0
div(E V

#)=0
- Dans le cas stationnaire (pas de transitoire), l’équation devient: div(ρ.E V

4.2.2 Conservation de la quantité de mouvement

On reprend le même principe que pour la conservation de la masse, et on étudie
les variations temporelles de la quantité de mouvement contenue dans le volume
Ω La quantité de mouvement s’écrit :
!!!
#
P (t) = # dv
ρE V (153)
Ω

C’est une quantité vectorielle. La relation générale 138 a été exprimée pour un
scalaire. Elle est donc valable pour chaque composante du vecteur (cela s’applique
aux tenseurs d’un ordre supérieur), donc pour le vecteur lui-même.
La variation de la quantité de mouvement en fonction du temps s’écrit:

dP# (t) d #
∂ρE V −→
!!! !!!
E# # ⊗E V
# ))dv
= ρ V dv = ( + div(ρE V (154)
dt dt Ω Ω ∂t
−→ # ⊗E V
# )) = (ρVi Vj ),i = ρ,i Vi Vj + ρVi,i Vj + ρVi Vj,i
On développe: div(ρE V

dP# (t) ∂ρ E # # −−→

∂E V
!!!
= V +ρ # ).E V# +(ρdiv(E V# )).E V
+(grad(ρ).E V # +ρgrad(E V# ).E V
# dv
dt Ω ∂t ∂t
(155)
dP# (t) ∂ρ −−→ ∂ E#
V
!!!
= ( + grad(ρ).E V# + ρdiv(E V# ))E V# + ρ( + grad(E V# ).E V
# )dv
dt Ω ∂t ∂t
(156)
On reconnaît l’expression de la conservation de la masse et l’accélération eu-
lerienne, soit:

dP# (t) #
∂E V
!!! !!!
= ρ( # ).E V# )dv =
+ grad(E V ρE #Γdv (157)
dt Ω ∂t Ω

47
 Il faut ensuite faire le bilan des éléments qui peuvent faire varier la quantité de
mouvement au sein de Ω :
-les efforts volumiques : f#v => f# dv
"""
Ω v
-les efforts exercés sur la surface ∂Ω (contraintes): T# (M,#n) = σ.#n => ∂Ω T# ds
""

La conservation de la quantité de mouvement s’exprime donc selon:

!!! !!! !!
E#
ρ Γdv = f#v dv + σ.#nds (158)
Ω Ω ∂Ω
Soit, par application du théorème de Gauss-Ostrogradski:
−→
!!! !!! !!!
E#
ρ Γdv = #
fv dv + div(σ)dv (159)
Ω Ω Ω

En regroupant et en constatant que cette équation est indépendante du volume:

−→
div(σ) + f#v = ρE #Γ (160)
& Différentes explicitations de la relation de conservation de la quantité de
mouvement (sans loi de comportement)
−→
– fluides: hydrostatique: σij = −pδij soit div(σ) = −(pδij ),j = −p,i =
−−→
−grad(p)
L’accélération est par ailleurs nulle. La seule force agissante est le poids,
−−→
soit: f#v = ρ#g . On obtient: grad(p) = ρ#g .
−→
– solides: statique-élastique (élasticité): accélération nulle: div(σ) + f#v = #0

48