Descriptions du mouvement

• A la différence de la mécanique des solides indéformables, la mécanique des milieux continus permet de prendre en compte
les déformations d’un corps.
.
Considérons un milieu continu occupant un volume V à l’instant initial (t =0), par exemple une
balle en caoutchouc avant son écrasement dans la paume d’une main. Cette balle peut être vue
comme l’assemblage d’une infinité de petits éléments de matière appelés “points matériels”.
Chaque point matériel va se déplacer et avoir sa propre trajectoire.

Considérons un domaine 2D qui se déforme selon un parallélogramme.

Supposons la transformation 𝑥Ԧ = 𝜙 𝑋, 𝑡 .
Position initiale Temps
Une fois la transformation du milieu continu définie, il est facile de définir les
notions de déplacement, de vitesse et d’accélération:
𝑢 𝑋, 𝑡 = 𝑥Ԧ 𝑋, 𝑡 - 𝑋
𝑑 𝑋, 𝑡 variables de Lagrange
𝑣Ԧ 𝑋, 𝑡 = 𝑢 𝑋, 𝑡
𝑑𝑡 𝑥,
Ԧ 𝑡 variables d'Euler
Glissement 𝑑
𝑎Ԧ 𝑋, 𝑡 = 𝑣Ԧ 𝑋, 𝑡
𝑑𝑡
Mécanique des milieux continus - Pr. Bouchra Kouddane 1
 Déformations
Gradient de la transformation
On considère le point matériel P dans la configuration de Lors de la transformation 𝜙, le vecteur élémentaire d𝑋 devient :
référence C0, après transformation, la position spatiale P’
dans la configuration Ct . Une particule matérielle M,
respectivement M’ , située au voisinage de P, respectivement
de P’, permet de définir les deux vecteurs élémentaires On définit alors le tenseur gradient de la transformation comme :
d𝑋 𝑒𝑡 𝑑𝑥:
Ԧ

𝑑𝑋 𝑑𝑥Ԧ 𝑑𝑥Ԧ 𝑑𝑋
Point de vue lagrangien Point de vue Eulérien

Mécanique des milieux continus - Pr. Bouchra Kouddane 2

 Déformations
Exemple : Etat de déformation homogène triaxiale
On considère une déformation homogène triaxiale définie par Par simple dérivation des coordonnées actuelles, on
les relations suivantes : peut obtenir les composantes du tenseur gradient. Ce
dernier prend une forme diagonale.

En prenant en compte l’hypothèse que le volume

dans la configuration de référence est unitaire, le Le déterminant du tenseur gradient donne la variation relative
volume dans la configuration actuelle est donné par de volume dans notre transformation.
le produit des valeurs du tenseur.

Il correspond aussi au déterminant de la matrice.

Mécanique des milieux continus - Pr. Bouchra Kouddane 3
Déformations
Tenseur des dilatations
Pour caractériser les déformations, il faut traduire le fait que les distances entre deux points et les angles entre deux directions vont
changer lors de la transformation. Nous allons dans un premier temps évaluer le produit scalaire entre deux vecteurs élémentaires de
la configuration initiale lorsqu’ils se transforment en 𝑑𝑥1 et 𝑑𝑥2 dans la configuration courante. Le produit scalaire de deux vecteurs
est fonction à la fois de la norme de ces vecteurs, donc des distances entre deux points, et de l’angle formé entre ces vecteurs.
Algèbre tensoriel
𝑩𝑉⋅𝑈 =𝑉⋅𝑩 ന 𝐓𝑈
𝑉⋅𝑨 ന𝑈 = 𝑨ന𝐓𝑉 ⋅ 𝑈
Produit contracté ന (𝑩 𝑪ന )=(𝑨
𝑨 ന 𝑩 ) 𝑪ന = 𝑨ന 𝑩 𝑪ന

Le tenseur 𝐂Ӗ = 𝐅ധ 𝐓 · 𝐅ധ est appelé tenseur des dilatations de Cauchy-

Green droit (Lagrangien).

Propriétés du tenseur 𝐂ധ

Le tenseur 𝐂ധ est donc

symétrique défini positif.
Mécanique des milieux continus - Pr. Bouchra Kouddane 4
Déformations
Tenseur des dilatations
Inversement, on peut exprimer le produit scalaire des vecteurs
élémentaires dans la configuration de référence à partir des
vecteurs dans la configuration actuelle :
ന =𝑭
Où 𝑩 ന·𝑭ന 𝑻 est appelé tenseur des dilatation de Cauchy-Green
gauche (Eulérien).

Exemple : Etat de déformation homogène triaxiale (suite)

Avec la formule 𝐂ധ = 𝐅ധ 𝐓 · 𝐅ധ , on peut facilement Contrainte
calculer les composantes du tenseur de Cauchy
Green droit qui est alors un tenseur diagonal.
Les tenseurs 𝑩 ന tout comme 𝑭
ന et 𝑪, ന , ne
sont pas des mesures de déformations
pertinentes car pour un mouvement de
ന=𝑩
corps rigide, on a 𝑪 ന = ധ𝑰.

Mécanique des milieux continus - Pr. Bouchra Kouddane 5

 Déformations
Tenseurs des déformations

Une manière d’appréhender la notion de déformation est de On peut également calculer la différence en fonction du
calculer la différence : tenseur des dilatations de Cauchy-Green gauche :

ന est le tenseur des déformations de Green-Lagrange.

𝑬 𝒆ധ est le tenseur des déformations de Euler-almansi. C’est
C’est un tenseur du deuxième ordre symétrique et un tenseur du deuxième ordre symétrique et eulérien. Il peut
ന:
lagrangien. Il peut s’exprimer à partir de 𝑭 ന:
s’exprimer à partir de 𝑭

Mécanique des milieux continus - Pr. Bouchra Kouddane 6

Déformations
Interprétation des résultats : Dilatation linéaire
La connaissance de 𝑪 ന ou 𝑬
ന en un point M de la configuration initiale C0 permettent d’évaluer la dilatation,
autrement-dit la variation relative de longueur, de n’importe quel vecteur élémentaire 𝑑 𝑥Ԧ en ce point au cours de
la transformation.
La dilatation λ dans la
direction de 𝑑𝑋 ∶

En introduisant les vecteurs unitaires 𝑁 et 𝑛

La dilatation en fct du tenseur des dilatations

de Cauchy-Green droit

Dans une base 𝐞𝐢 𝐞𝐣 𝐞𝐤

Attention, il s’agit du terme diagonal
ന et non d’un indice muet répété
de 𝑪

Mécanique des milieux continus - Pr. Bouchra Kouddane 7

Déformations
Interprétation des résultats : Déformation dans une direction

Grâce au tenseur de déformation introduits précédemment, il est maintenant possible de calculer la déformation
dans une direction donnée.
On considère deux vecteurs 𝑑X1 et 𝑑𝑋2 égaux, on note alors 𝑑X1 = 𝑑𝑋2 = 𝑑X. Au cours de la transformation,
ce vecteur sera transformé en 𝑑x1 = 𝑑x2 = 𝑑x .
En décomposant les vecteurs 𝑑X et 𝑑x selon leur norme et leur orientation, on obtient :

On peut écrire :

En se rappelant que l’expression de la dilatation dans la direction N secrit 𝜆 𝑁 = 𝑑𝑙/𝑑𝐿, on peut déduire
immédiatement que :
La déformation en fct du tenseur des
déformations de Green-Lagrange

Mécanique des milieux continus - Pr. Bouchra Kouddane 8

Déformations
Interprétation des résultats : allongement unitaire dans une direction

Toujours dans le but de comparer l’état du milieu entre les deux configurations, nous pouvons évaluer la variation
relative de longueur que subit un vecteur 𝑑X au cours de la transformation. Pour cela, on appelle allongement unitaire
dans la direction 𝑑X la quantité adimensionnelle :

En décomposant les vecteurs 𝑑X et 𝑑x selon leur norme et leur orientation, on obtient :

Pour résumer :

Mécanique des milieux continus - Pr. Bouchra Kouddane 9

 Déformations
Interprétation des résultats : Distorsion angulaire
Glissement de deux vecteurs On peut également se demander s’il est possible d’évaluer les changements éventuels de
initialement orthogonaux l’angle formé par deux vecteurs d𝑋1 et d𝑋2 orthogonaux dans la configuration initiale.

Si on note α l’angle formé par ces deux vecteurs transformés, on nomme

glissement des deux vecteurs initialement orthogonaux d𝑋1 et d𝑋2 l’angle
défini par :
𝜋
cos(𝛼)= cos( − 𝛾)=sin(𝛾)
2
Le produit scalaire des vecteurs
dans la configuration actuelle : d𝑥1 ⋅ d𝑥2 = d𝑥1 d𝑥2 cos(𝛼)

cos(𝛼) =

En introduisant les En appliquant cette formule à

vecteurs unitaires : deux vecteurs de base 𝑒𝑖 et 𝑒𝑗
NB: C’est là un des avantages des
tenseurs 𝐂ധ et 𝑬 ന , qui permettent de
calculer à la fois les variations de
L’angle 𝛾 peut aussi longueur (composantes diagonales) et
être exprimé par : d’angles (composantes hors-diagonales)
Mécanique des milieux continus - Pr. Bouchra Kouddane au cours de la transformation. 10
 Déformations
Exemple : Etat de déformation homogène triaxiale (suite)
Avec la formule 𝐂ധ = 𝐅ധ 𝐓 · 𝐅ധ , on peut facilement
calculer les composantes du tenseur de Cauchy Les dilatations définies par les vecteurs 𝑒1 𝑒2 𝑒3 :
Green droit qui est alors un tenseur diagonal.
𝜆 𝑒1 = 𝐶11 = 𝜆1
𝜆 𝑒2 = 𝐶22 = 𝜆2
𝜆 𝑒3 = 𝐶33 = 𝜆3

ന = 𝟏 (𝑪
Ensuite, la formule 𝑬 ന − ധ𝑰) nous permet On remarque que si 𝜆𝑖𝑖 >1, cette déformation est
𝟐
d’avoir les composantes du tenseur de Green effectivement une dilatation, mais peut également être
Lagrange. appelé une « contraction » dans le cas contraire.
Dans ce cas, nous n'avons pas de glissement.

Mécanique des milieux continus - Pr. Bouchra Kouddane 11

Déformations
Déplacements Vecteur position actuelle

Vecteur position initiale

Le déplacement du point M à l'instant t : 𝑈 𝑋, t = 𝑀𝑀′ = 𝑂𝑀′- 𝑂𝑀 = 𝑥(𝑋,t)−𝑋
Ԧ
𝜕𝑈𝑖 𝜕𝑥𝑖
Après dérivation on a : = − 𝛿𝑖𝑗
𝜕𝑋𝑗 𝜕𝑋𝑗
𝐼Ӗ
ന
𝐻 𝐹ധ
𝜕𝑈𝑖
ന = 𝐺𝑟𝑎𝑑 𝑈 =
Donc on a le tenseur gradient du déplacement : 𝐻 = 𝐹ധ − 𝐼 Ӗ
𝜕𝑋𝑗

On peut donc décomposer le tenseur gradient du déplacement en deux tenseurs symétrique et antisymétrique comme suit:

ന = 𝜺ധ + 𝝎
𝑯 ന

ന𝑠 = 1 𝐻
𝜀Ӗ = 𝐻 ന+𝐻
ന𝑇 𝜔 ന𝑎 = 1 𝐻
ന=𝐻 ന−𝐻
ന𝑇
2 2

1 𝜕𝑈𝑖 𝜕𝑈𝑗 1 𝜕𝑈𝑖 𝜕𝑈𝑗

𝜀𝑖𝑗 = ( + ) 𝜔𝑖𝑗 = ( − )
2 𝜕𝑋𝑗 𝜕𝑋𝑖 2 𝜕𝑋𝑗 𝜕𝑋𝑖
Tenseur de Tenseur de rotation
déformation linéarisé
Mécanique des milieux continus - Pr. Bouchra Kouddane 12
Déformations
Déplacements
On peut donc écrire les tenseurs 𝐶Ӗ et 𝐸ധ à l’aide de 𝐻
ന:
ന = 𝐹ധ − 𝐼ന , donc : 𝐹ധ = 𝐼 Ӗ + 𝐻
Comme : 𝐻 ന Car 𝜀 Ӗ =
1
ന+𝐻
𝐻 ന𝑇
2
On sait que : 𝐶Ӗ = 𝐹ധ 𝑇 · 𝐹ധ
Donc : 𝐶Ӗ = (𝐼ന + 𝐻
ന 𝑇 ) 𝐼ന + 𝐻
ന = 𝐼ന + 𝐻
ന+𝐻
ന𝑇 + 𝐻
ന𝑇 𝐻
ന = 𝐼ന + 2𝜀 Ӗ + 𝐻
ന𝑇 𝐻
ന

𝐶Ӗ = 𝐼ന + 2𝜀 Ӗ + 𝐻
ന𝑇 𝐻
ന

1
D'autre part : 𝐸ധ = 𝐶Ӗ − 𝐼 Ӗ
2
1
ന𝑇 ) + 1 𝐻 𝑇 𝐻
Donc : 𝐸ധ = (𝐻 + 𝐻 ന = 𝜀Ӗ + 1 𝐻
ന𝑇 𝐻
ന
2 2 2

1 𝑇
𝐸ധ = 𝜀 Ӗ + 𝐻
ന 𝐻ന
2

Mécanique des milieux continus - Pr. Bouchra Kouddane 13

Déformations
Déplacements
Décomposition polaire ou décomposition du déplacement
Le tenseur gradient de la transformation 𝐹ധ permet de relier localement les configurations initiales et actuelles. Cette
transformation peut toujours être décomposée en une étape de déformation pure et une étape de rotation. Cette décomposition
unique du tenseur gradient de transformation, appelée décomposition polaire, s’écrit sous la forme :
Où 𝑅ധ est le tenseur de rotation. On rappelle que 𝑅ധ est un tenseur orthogonal c’est-à-dire
que sa transposée coïncide avec son inverse. Les tenseurs 𝑈 ന et 𝑉ധ sont des tenseurs
symétriques définis positifs, appelés respectivement tenseur des déformations pures à
droite et tenseur des déformations pures à gauche. Les tenseurs U et V sont liés par la
relation :

Dans le cas particulier d’un mouvement de corps rigide, on a :

Le tenseur de Cauchy-Green droit 𝐶Ӗ en fonction du tenseur 𝑈

ന:
Le tenseur 𝑈ന est donc la
racine carrée du tenseur 𝐶Ӗ
Le tenseur de Cauchy-Green gauche 𝐵ധ en fonction du tenseur 𝑉ധ :
Le tenseur 𝑉ധ est donc la
racine carrée du tenseur 𝐵ധ

Mécanique des milieux continus - Pr. Bouchra Kouddane 14

Déformations
Déplacements
Le tenseur des déformations de Green-Lagrange 𝐸ധ en fonction du tenseur 𝑈
ന:

Le tenseur des déformations de Euler-almansi 𝑒Ӗ en fonction du tenseur 𝑉ധ :

NB : Bien que 𝑈 ന et 𝑉ധ soient appelés tenseurs de déformation pure, ils sont principalement associés à la dilatation
(changement de volume) et à la distorsion (changement de forme sans modification de volume) plutôt qu'à une déformation
pure au sens strict (sans changement de volume).

Mécanique des milieux continus - Pr. Bouchra Kouddane 15

Déformations
Hypothèse des Petites Perturbations (HPP)
Déformations en petites perturbations :
Cette hypothèse est appliquée lorsque les déplacements et déformations d'une structure sont petits par rapport à ses
dimensions initiales. Cela permet de simplifier les calculs des efforts et des contraintes.

Exemple dans le cadre du bâtiment :

Lorsqu'un immeuble est soumis à une charge (comme le vent ou le poids propre), son
déplacement par rapport à sa position initiale non chargée est faible. Il se déplace peu,
ce qui rend valide l'hypothèse de petites perturbations.

Exemple dans les essais de traction :

Lorsqu'on effectue un essai de traction sur une éprouvette métallique, les déplacements
et les déformations sont généralement faibles dans le domaine élastique (et même au
début de la zone plastique). Ici aussi, les déformations étant faibles par rapport à la
taille de l'éprouvette, l'HPP reste valable.

Mécanique des milieux continus - Pr. Bouchra Kouddane 16

 Déformations
Hypothèse des Petites Perturbations (HPP)
Les déformations de Green-Lagrange et Euler-Almansi dépendent des déplacements de manière non-linéaire :
1 𝑇
𝐸ധ = 𝜀 Ӗ + 𝐻
ന 𝐻ന
2 expression non-linéaire (quadratique) des déplacements

Peut-on alors linéariser ? Oui, on utilisant l'hypothèse de petites perturbations (HPP)

H.P.P : les déplacements entre la configuration de référence et la configuration actuelle sont très petits et le gradient
des déplacements est également petit.
petits déplacements = petites transformations = petites déformations

Puisque 𝐹ധ = 𝐼 Ӗ + 𝐻,
ന dans le cas des petites perturbations le tenseur gradient de la transformation 𝐹ധ est proche du
tenseur unité 𝐼,Ӗ et 𝑑𝑥 est proche de 𝑑𝑋
Configuration actuelle = Configuration de référence

Ainsi, on va pouvoir négliger les termes non linéaires dans l'expression de 𝐸ധ et 𝐶.Ӗ Ainsi on linéarise la cinématique,
ce qui simplifie grandement la résolution finale du problème.
On peut donc déduire que : ന ≃ 𝒆ധ ≃ 𝜺ധ
𝑬 et ന ≃ 𝑰ന + 2ധ𝜺
𝑪
Mécanique des milieux continus - Pr. Bouchra Kouddane 17
Déformations
Hypothèse des Petites Perturbations (HPP)
Dilatation et allongement unitaire dans une direction
Rappel: La relation entre
allongement et dilatation : =

1
On a : 𝜆 𝑢0 = 1 + δ 𝑢0 ധ 0 , 𝑢0 )]1/2 ≃ [1 + 2𝜀(𝑢
= [1 + 2𝐸(𝑢 Ӗ 0 , 𝑢0 )]2 ≃ 1 + 𝜀(𝑢
Ӗ 0 , 𝑢0 ) = 1+ 𝑢0 ⋅ 𝜀 Ӗ 𝑢0

Ӗ 0 , 𝑢0 ) = 𝑢0 ⋅ 𝜀 Ӗ 𝑢0
Donc l'allongement : δ(𝑢0 ) ≃ 𝜀(𝑢 Développement limité

Ӗ 𝑖 ,𝑒𝑖 ) = 𝑒𝑖 ⋅ 𝜀 Ӗ 𝑒𝑖 = 𝜀𝑖𝑖Ӗ
Si 𝑢0 = 𝑒𝑖 on obtient : δ(𝑒𝑖 ) ≃ 𝜀(𝑒 𝜀11 = δ 𝑒1
𝜀22 = δ 𝑒2
𝜀33 = δ 𝑒3
Donc : 𝛿𝑖 = 𝜆𝑖 − 1 ≃ 𝜀𝑖 ≃ 𝐸𝑖

Mécanique des milieux continus - Pr. Bouchra Kouddane 18

Déformations
Hypothèse des Petites Perturbations (HPP)
Glissement de deux directions orthogonales
On sait que : 𝑑𝑥 ⋅ 𝑑𝑥′- 𝑑𝑋 ⋅ 𝑑𝑋 ′ = 𝑑𝑋 ⋅ (𝐶Ӗ − 𝐼)Ӗ ⋅ 𝑑𝑋 ′ = 2𝑑𝑋 ⋅ 𝐸ന 𝑑𝑋 ′
Ӗ 0 ,𝑢0 ′) = 2𝐸(𝑢
Soit 𝑢0 et 𝑢0 ′ deux vecteurs unitaires tels que 𝑢0 ⋅ 𝑢0 ′ = 0 => 𝐶(𝑢 ധ 0 ,𝑢0 ′)

Ӗ 0 ,𝑢0 ′)
𝐶(𝑢 ധ 0 ,𝑢0 ′)
2𝐸(𝑢
sin 𝛾 𝑢0 , 𝑢0 ′ = 𝜆(𝑢0 )𝜆(𝑢0 ′)
= 𝜆(𝑢0 )𝜆(𝑢0 ′)
sin 𝛾 = 𝛾 dans HPP
𝜀 Ӗ ≪ 1 dans HPP
2ധ𝜀 (𝑢0 ,𝑢0 ′)
𝛾 𝑢0 , 𝑢0 ′ ≃ Ӗ 0 ,𝑢0 ′) = 2 𝑢0 ′ ⋅ 𝜀 Ӗ 𝑢0
≃ 2𝜀(𝑢
[1+ധ𝜀(𝑢0 ,𝑢0 )][1+ 𝜀ധ (𝑢0 ′,𝑢0 ′)]

Si 𝑢0 = 𝑒𝑖 et 𝑢′0 = 𝑒𝑗 (i ≠ 𝑗), on obtient : 𝛾(𝑒𝑖 ,𝑒𝑗 ) ≃ 2𝜀(𝑒

Ӗ 𝑖 ,𝑒𝑗 ) = 2𝜀𝑖𝑗

Mécanique des milieux continus - Pr. Bouchra Kouddane 19

Déformations
Hypothèse des Petites Perturbations (HPP)
Dilatation volumique
Dans la même logique que pour l'allongement unitaire, on peut exprimer la dilatation volumique de la manière
suivante :
𝑑𝑣 − 𝑑𝑣0 𝑑𝑣
= −1=𝐽−1
𝑑𝑣0 𝑑𝑣0

𝐽 = 𝜆1 𝜆2 𝜆3 ≃ (1 + 𝜀1 ) (1 + 𝜀2 ) (1 + 𝜀3 ) ≃ 1 + 𝜀1 + 𝜀2 + 𝜀3

𝑑𝑣−𝑑𝑣0
= 𝐽 − 1 = 𝜀1 + 𝜀2 + 𝜀3 = 𝑡𝑟𝜀 Ӗ
𝑑𝑣0

NB : La trace de 𝜀 Ӗ représente la variation de volume d'un élément élémentaire représentatif

Mécanique des milieux continus - Pr. Bouchra Kouddane 20

Déformations
Hypothèse des Petites Perturbations (HPP)
Quelques définitions relatives au tenseur 𝜀 Ӗ
• 𝜀 Ӗ est un tenseur symétrique.
Ԧ 𝑖=1,2,3 auxquelles sont associées des valeurs
• 𝜀 Ӗ possède un système orthonormé de vecteurs propres (𝐴)
propres.
• Un tenseur 𝜀 Ӗ est de révolution autour de la direction principale 𝐴3 si 𝜀1 = 𝜀2 ≠ 𝜀3 => Toute direction
orthogonale à 𝐴3 est direction principales associée à 𝜀1 = 𝜀2
• Un tenseur 𝜀 Ӗ est sphérique si 𝜀1 = 𝜀2 = 𝜀3 => toute direction de l'espace est principale
• Invariants élémentaires du tenseur 𝜀 Ӗ : Indépendants de la base choisie
𝜀𝐼 = 𝜀1 + 𝜀2 + 𝜀3 = 𝑡𝑟𝜀 Ӗ => dilatation volumique
1 2
𝜀𝐼𝐼 = 𝜀1 𝜀2 + 𝜀2 𝜀3 + 𝜀3 𝜀1 = 2 𝑡𝑟𝜀 Ӗ − 𝑡𝑟 𝜀 Ӗ2
𝜀𝐼𝐼𝐼 = 𝜀1 𝜀2 𝜀3 = det 𝜀 Ӗ

Mécanique des milieux continus - Pr. Bouchra Kouddane 21

Déformations
Hypothèse des Petites Perturbations (HPP)
Dilatation isotrope
𝑥Ԧ = 𝜆X
Soit 𝑥1 = 𝜆𝑋1 , 𝑥2 = 𝜆𝑋2 , 𝑥3 = 𝜆𝑋3 avec 𝜆 > 0
La dilatation dans la direction 𝑢0

Approximation en HPP (Hypothèse de Petites

Perturbations) : Pour de petites perturbations
(HPP), où λ∼1, on peut approximer le tenseur de
Green-Lagrange 𝐸ധ :

Mécanique des milieux continus - Pr. Bouchra Kouddane 22

Déformations
Etats de déformation particuliers
Extension simple
L'image présente une extension simple dans la direction de 𝑒1

Hypothèse de petites perturbations (HPP) :

Lorsque λ∼1, le tenseur de Green-Lagrange peut être approximé :

Mécanique des milieux continus - Pr. Bouchra Kouddane 23

Déformations
Etats de déformation particuliers
Glissement simple
𝑥1 = 𝑋1 + 2𝜈𝑋2 𝑥1 = 𝑋2 𝑥3 = 𝑋3

Hypothèse de petites perturbations (HPP) :

𝐸ധ est exprimé en fonction de 𝜈 et contient des termes quadratiques, tels
que 𝜈 2 . En utilisant l'hypothèse de petites perturbations, on considère que
ces termes quadratiques sont très petits comparés aux termes linéaires et
qu'ils peuvent donc être négligés.

Autrement dit :
Mécanique des milieux continus - Pr. Bouchra Kouddane 24
Déformations
Etats de déformation particuliers
Déformation plane
La déformation plane se produit lorsque les déplacements dans une structure ou un matériau se produisent principalement
dans un plan, et que les déformations dans la direction perpendiculaire à ce plan sont négligeables. Cela signifie que l'on peut
analyser le problème en considérant uniquement deux dimensions (2D).
𝑥1 = 𝑥1 (𝑋1 , 𝑋2 , 𝑡) 𝑥2 = 𝑥2 𝑋1 , 𝑋2 , 𝑡 𝑥3 = 𝑋3

× × 0
𝐹ധ et 𝐶Ӗ sont de la forme : × × 0
0 0 1
× × 0
ധ 𝐻,
𝐸, ന 𝜀 Ӗ et 𝜔
ന sont alors de la forme : × × 0 => La déformation est plane
0 0 0
𝜀1 0 0
Dans la base principale (𝐴1 , 𝐴2 , 𝐴3 ): 0 𝜀2 0
0 0 0
Les éléments diagonaux dans cette base correspondent aux valeurs propres du tenseur de déformation et sont appelés
déformations principales.
Mécanique des milieux continus - Pr. Bouchra Kouddane 25
Déformations
Conditions de compatibilité
Dans les sections précédentes, nous avons montré qu’il est assez simple de calculer le tenseur des déformations linéarisé par
dérivation du champs de déplacement. En revanche, l’inverse n’est pas vrai en général. En effet, pour qu’un tenseur
symétrique de second ordre corresponde à un tenseur de déformation, autrement dit qu’il soit intégrable et associé à un
champs de déplacement continu, le tenseur doit respecter les relations de compatibilité pour garantir la continuité du milieu
continu lors de sa déformation.

Champ de déformations compatible

Champ de déformations non compatible

Champ de déplacement Champ de déformation

Mécanique des milieux continus - Pr. Bouchra Kouddane 26

Déformations
Conditions de compatibilité
ന = 𝑮𝒓𝒂𝒅 𝑼 ?
Problématique : Existe-t-il un champ de déplacement 𝑼(𝑴𝟎 ,t) tel que 𝜺ധ soit la partie symétrique de 𝑯
𝜕𝑈𝑖
ന = 𝜀Ӗ + 𝜔
le tenseur gradient du déplacement est exprimé par : 𝐻 ന 𝐻𝑖𝑗 = = 𝑈𝑖,𝑗
𝜕𝑋𝑗
1 𝜕𝑈𝑖 𝜕𝑈𝑗 1 1 𝜕𝑈𝑖 𝜕𝑈𝑗 1
Avec: 𝜀𝑖𝑗 = + = (𝑈𝑖,𝑗 + 𝑈𝑗,𝑖 ) 𝜔𝑖𝑗 = − = (𝑈𝑖,𝑗 − 𝑈𝑗,𝑖 )
2 𝜕𝑋𝑗 𝜕𝑋𝑖 2 2 𝜕𝑋𝑗 𝜕𝑋𝑖 2

ന a partir de 𝜀 Ӗ il faut déjà calculer 𝜔

Pour trouver 𝐻 ന et 𝜀 Ӗ , on dérive alors 𝜔𝑖𝑗 par rapport a k :
ന , l'idée est de trouver un lien entre 𝜔
1 1 1
𝜔𝑖𝑗,𝑘 = (𝑈𝑖,𝑗𝑘 − 𝑈𝑗,𝑖𝑘 ) = (𝑈𝑖,𝑗𝑘 + 𝑈𝑘,𝑖𝑗 − 𝑈𝑘,𝑖𝑗 − 𝑈𝑗,𝑖𝑘 ) = 𝑈𝑖,𝑘 ,𝑗
+ 𝑈𝑘,𝑖 ,𝑗
− 𝑈𝑘,𝑗 − 𝑈𝑗,𝑘 = 𝜀𝑖𝑘,𝑗 − 𝜀𝑘𝑗,𝑖
2 2 2 ,𝑖 ,𝑖

𝜔𝑖𝑗,𝑘 = 𝜀𝑖𝑘,𝑗 − 𝜀𝑘𝑗,𝑖 𝜔𝑖𝑗,𝑘𝑙 = 𝜔𝑖𝑗,𝑙𝑘 pour 𝑙 ≠ 𝑘 Conditions d'intégrabilité au sens de Cauchy
Condition nécessaire et suffisante
pour système intégrable
On remplace et on obtient : 𝜀𝑖𝑘,𝑗𝑙 − 𝜀𝑘𝑗,𝑖𝑙 = 𝜀𝑖𝑙,𝑗𝑘 − 𝜀𝑙𝑗,𝑖𝑘 6 relations de compatibilité
Si les 6 relations de compatibilité sont vérifiées donc on peut intégrer 𝜔𝑖𝑗,𝑘 = 𝜀𝑖𝑘,𝑗 − 𝜀𝑘𝑗,𝑖 ന = 𝜀Ӗ + 𝜔
𝐻 ന
ന
Pour trouver 𝑈 on doit intégrer 𝐻 𝑈𝑖,𝑗 = 𝜀𝑖𝑗 + 𝜔𝑖𝑗 est-il intégrable ? Oui si 𝑈𝑖,𝑗𝑙 = 𝑈𝑖,𝑙𝑗 soit 𝜀𝑖𝑗,𝑙 + 𝜔𝑖𝑗,𝑙 = 𝜀𝑖𝑙,𝑗 + 𝜔𝑖𝑙,𝑗
Déjà satisfaite puisque : 𝜔𝑖𝑗,𝑙 = 𝜀𝑖𝑙,𝑗 − 𝜀𝑙𝑗,𝑖 et 𝜔𝑖𝑙,𝑗 = 𝜀𝑖𝑗,𝑙 − 𝜀𝑗𝑙,𝑖

Mécanique des milieux continus - Pr. Bouchra Kouddane 27

Déformations
Conditions de compatibilité
𝜕 2 𝜀11 𝜕 𝜕𝜀13 𝜕𝜀21 𝜕𝜀32
= + −
𝜕𝑋2 𝜕𝑋3 𝜕𝑋1 𝜕𝑋2 𝜕𝑋3 𝜕𝑋1
𝜕 2 𝜀22 𝜕 𝜕𝜀21 𝜕𝜀32 𝜕𝜀13
= + −
𝜕𝑋1 𝜕𝑋3 𝜕𝑋2 𝜕𝑋3 𝜕𝑋1 𝜕𝑋2
𝜕 2 𝜀33 𝜕 𝜕𝜀32 𝜕𝜀13 𝜕𝜀21
= + −
𝜕𝑋1 𝜕𝑋2 𝜕𝑋3 𝜕𝑋1 𝜕𝑋1 𝜕𝑋3
𝜕2 𝜀12 𝜕2 𝜀11 𝜕2 𝜀22
2 = 2 +
𝜕𝑋1 𝜕𝑋2 𝜕𝑋2 𝜕𝑋1 2
𝜕2 𝜀23 𝜕2 𝜀22 𝜕2 𝜀33
2 = +
𝜕𝑋2 𝜕𝑋3 𝜕𝑋3 2 𝜕𝑋2 2
𝜕2 𝜀31 𝜕2 𝜀33 𝜕2 𝜀11
2 = +
𝜕𝑋3 𝜕𝑋1 𝜕𝑋1 2 𝜕𝑋3 2

Mécanique des milieux continus - Pr. Bouchra Kouddane 28