Université Hassan II

École Nationale Supérieure d’Électricité et de Mécanique

Département Génie Mécanique

P O L Y C O P I E
Élément de module

Mécanique des Solides Déformables

Théorie, Applications et Simulation Numérique.

Cours enseigné par : Pr. JALAL Salah Eddine

Année universitaire : 2022 - 2023

Destiné aux étudiants de la 1ère année de l’ENSEMľ - Spécialité : Génie Mécanique.

Mécanique des Solides 2 1ère année du cycle d’Ingénieurs
Table des matières

Introduction 1

1 Notion du milieu continu 3

1.1 Concept du milieu continu . . . . . . . . . . . . . . . . . . . 3
1.1.1 Volume élémentaire représentatif . . . . . . . . . . . 3
1.1.2 Concept d’un milieu continu . . . . . . . . . . . . . 4
1.1.3 Forces . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.4 Lois de conservation . . . . . . . . . . . . . . . . . . 5
Conservation de la masse . . . . . . . . . . . . . . . 5
1.1.5 Lois de mouvement (Euler 1776) . . . . . . . . . . . 6
Énergie . . . . . . . . . . . . . . . . . . . . . . . . . 7
Remarques . . . . . . . . . . . . . . . . . . . . . . . 7

2 Calcul indiciel 9
2.1 Principe d’objectivité . . . . . . . . . . . . . . . . . . . . . . 9
2.2 Tenseurs et appareil mathématique . . . . . . . . . . . . . . 9
2.3 Calcul indiciel . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . 10
2.3.2 Règles importantes . . . . . . . . . . . . . . . . . . 10
2.3.3 Symbole de Kronecker . . . . . . . . . . . . . . . . . 11
2.3.4 Transformation des coordonnées . . . . . . . . . . . 12
Rotation des axes . . . . . . . . . . . . . . . . . . . 12
Transformation des coordonnées . . . . . . . . . . . 14
2.4 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

3 Champ gradient de déformation et ses applications 17

3.1 Champ gradient de déformation . . . . . . . . . . . . . . . . 17
3.1.1 Déformation d’un milieu continu . . . . . . . . . . . 17
3.1.2 Tenseur gradient de la déformation . . . . . . . . . . 18
3.2 Champ de tenseur de dilatation . . . . . . . . . . . . . . . . 19
3.2.1 Transport du produit scalaire . . . . . . . . . . . . . 19
3.2.2 Tenseur de dilatation . . . . . . . . . . . . . . . . . 20

i
 TABLE DES MATIÈRES

3.2.3 Dilatation relative des longueurs . . . . . . . . . . . 21

3.2.4 Glissement des vecteurs orthogonaux . . . . . . . . . 21
3.2.5 Transport de la base canonique . . . . . . . . . . . . 22
3.2.6 Transport des volumes . . . . . . . . . . . . . . . . 23

4 Tenseur des contraintes et équations d’équilibre 27

4.1 Milieu déformable . . . . . . . . . . . . . . . . . . . . . . . 27
4.1.1 Forces externes et équilibre mécanique . . . . . . . . 27
4.1.2 Etat mécanique local . . . . . . . . . . . . . . . . . 27
Contraintes . . . . . . . . . . . . . . . . . . . . . . 28
Traction de surface . . . . . . . . . . . . . . . . . . 28
4.1.3 Repère local, contraintes normales et de tangentielles 29
Contrainte normale ff . . . . . . . . . . . . . . . . . 29
Contrainte de cisaillement fi . . . . . . . . . . . . . 29
4.1.4 Tenseur de contraintes . . . . . . . . . . . . . . . . 30
Équilibre des forces . . . . . . . . . . . . . . . . . . 30
Équilibre des moments . . . . . . . . . . . . . . . . 31
4.1.5 Propriétés du tenseur contraintes . . . . . . . . . . . 32
4.2 Équations de mouvement . . . . . . . . . . . . . . . . . . . 33
4.2.1 Parallélépipède élémentaire . . . . . . . . . . . . . . 33
4.2.2 Équations de mouvement . . . . . . . . . . . . . . . 33
4.2.3 Symétrie du tenseur contrainte et principe de récipro-
cité des ffij . . . . . . . . . . . . . . . . . . . . . . . 35
4.2.4 Équilibre statique . . . . . . . . . . . . . . . . . . . 36
4.3 Tenseur déviateur et tenseur sphérique . . . . . . . . . . . . 36
4.4 Espace des contraintes principales . . . . . . . . . . . . . . . 37

5 Tenseur des déformations et équations de compatibilité 41

5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . 41
5.2 État local de déformation . . . . . . . . . . . . . . . . . . . 41
5.2.1 Petites et grandes déformations . . . . . . . . . . . . 41
5.2.2 Principe de superposition . . . . . . . . . . . . . . . 42
5.2.3 Extension et glissement . . . . . . . . . . . . . . . . 43
5.3 Les grandes déformations . . . . . . . . . . . . . . . . . . . 43
5.4 Les petites déformations . . . . . . . . . . . . . . . . . . . . 45
5.4.1 Le vecteur déplacement . . . . . . . . . . . . . . . . 45
5.4.2 Gradient de déplacement et petites déformations . . . 46
5.4.3 Translation, rotation et déformation . . . . . . . . . 46
Partie antisymétrique du tenseur de déplacement . . . 47
Partie symétrique du tenseur gradient de déplacement 48
5.5 Tenseur des déformations . . . . . . . . . . . . . . . . . . . 48

Mécanique des Solides ii 1ère année du cycle d’Ingénieurs

 TABLE DES MATIÈRES

5.5.1 Définition . . . . . . . . . . . . . . . . . . . . . . . 48
5.5.2 Repère local, extension, distorsion et glissement . . . 48
5.5.3 Propriétés des déformations . . . . . . . . . . . . . . 49
5.6 Interprétation physique . . . . . . . . . . . . . . . . . . . . . 49
5.7 Glissement . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

6 Loi de comportement 53
6.1 Essais de la traction simple . . . . . . . . . . . . . . . . . . 53
6.1.1 Cas de l’acier doux . . . . . . . . . . . . . . . . . . 53
6.1.2 Contraction latérale . . . . . . . . . . . . . . . . . . 55
6.1.3 Cas de l’alliage d’aluminium . . . . . . . . . . . . . 55
6.2 Hypothèses de linéarité . . . . . . . . . . . . . . . . . . . . 56
6.2.1 Linéarité géométrique . . . . . . . . . . . . . . . . . 56
6.2.2 Linéarité matérielle . . . . . . . . . . . . . . . . . . 57
6.2.3 Linéarité des conditions aux limites . . . . . . . . . . 57
6.3 Aspects phénoménologiques . . . . . . . . . . . . . . . . . . 58
6.4 Élasticité linéaire . . . . . . . . . . . . . . . . . . . . . . . 58
6.4.1 Domaine de validité et d’emploi . . . . . . . . . . . 58
6.4.2 Formulation . . . . . . . . . . . . . . . . . . . . . . 59
Potentiel d’état . . . . . . . . . . . . . . . . . . . . 59
Élasticité linéaire isotrope : loi de comportement . . 60

7 Forme intégrale de l’équilibre et de la cinématique 63

7.1 Forme intégrale de l’équilibre . . . . . . . . . . . . . . . . . 63
7.1.1 L’intégration par partie est la clef de la forme intégrale 63
Remarque . . . . . . . . . . . . . . . . . . . . . . . 63
7.2 Forme intégrale de la cinématique . . . . . . . . . . . . . . 64
7.3 Principe des déplacements virtuels (ou des travaux virtuels) . . 64
Remarques . . . . . . . . . . . . . . . . . . . . . . . 65

8 Introduction à la Plasticité 67
8.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . 67
8.2 Comportement plastique en sollicitation uni-axiale . . . . . . 67
8.2.1 Seuil de la plasticité . . . . . . . . . . . . . . . . . . 67
8.2.2 Influence des déformations plastiques sur les proprié-
tés mécaniques . . . . . . . . . . . . . . . . . . . . . 69
Déformation plastique modérées . . . . . . . . . . . 69
Déformations plastiques importantes . . . . . . . . . 70
8.2.3 Comportement plastique pur . . . . . . . . . . . . . 70
L’élasto-plasticité parfaite . . . . . . . . . . . . . . . 70
L’élasto-plasticité écrouissante . . . . . . . . . . . . 71
8.3 Domaine d’élasticité en sollicitation quelconque . . . . . . . 73

1ère année du cycle d’Ingénieurs iii Mécanique des Solides

 TABLE DES MATIÈRES

8.3.1 Généralités . . . . . . . . . . . . . . . . . . . . . . . 73
Matériau parfaitement plastique . . . . . . . . . . . . 74
8.3.2 Quelques résultats expérimentaux concernant les mé-
taux . . . . . . . . . . . . . . . . . . . . . . . . . . 74
8.4 Théorie de l’élasto-plasticité . . . . . . . . . . . . . . . . . 76
8.4.1 Hypothèse de stabilité des matériaux, postulat de
Drucker . . . . . . . . . . . . . . . . . . . . . . . . 78
8.4.2 Propriétés . . . . . . . . . . . . . . . . . . . . . . . 80
Loi de normalité . . . . . . . . . . . . . . . . . . . . 80
Convexité . . . . . . . . . . . . . . . . . . . . . . . 80
8.4.3 Puissance de dissipation plastique . . . . . . . . . . . 81
8.4.4 Règles de l’écrouissage . . . . . . . . . . . . . . . . 82
Paramètres d’écrouissage . . . . . . . . . . . . . . . 82
Idéalisation des règles de l’écrouissage . . . . . . . . 82

Mécanique des Solides iv 1ère année du cycle d’Ingénieurs

Introduction

La mécanique des solides que vous avez étudié jusqu’à maintenant, est
limitée aux problèmes de mouvement des corps rigides, i.e. si l’on choisi
deux points différent faisant partie du milieu étudié, la distance qui les sé-
pare ne change pas quelques soient les sollicitations extérieures. Ceci donne
naissance au principe de l’équiprojectivité qui est considéré comme l’un des
piliers de la mécanique des solides indéformables.
Mais dans la réalité, presque tous les éléments structuraux et toutes les
structures sont déformables. Commençant par les disques dentées d’une
montre jusqu’au tripodes de télécommunication en passant par les pignons
des appareils électropositifs, les courroies de transmission de puissance, les
ponts, les pylônes, les câbles, le fuselage des avions, ... etc. D’où, la né-
cessité d’une nouvelle mécanique appelé la mécanique des milieux continus
déformables. Plus particulièrement, la mécanique des solides déformables.
Cette nouvelle science a pour objectif l’étude, d’une manière mathématique
rigoureuse, le comportement des solides soumis à des sollicitations exté-
rieures diverses. Ce comportement se caractérise par la transmission des
efforts à l’intérieur de la matière, ce qui cause l’évolution des déplacements
de chaque point matériel.
Dans ce cours, on établit des règles et relations fondamentales intrinsèques,
certaines indépendamment des propriétés particulières du matériau, d’autres
valables pour des classes étendues de matériaux. Ce caractère général donne
une grande puissance à ces sciences pour les applications de l’ingénieur.
Lorsque le matériau est solide et obéit à un comportement linéaire réver-
sible (loi de Hooke, linéarisation matérielle) et que les déformations et les
déplacements sont très petits (linéarisation géométrique), on aboutit à la
théorie de l’élasticité, branche la plus développée de la mécanique des so-
lides.
La mécanique des solides établit un ensemble d’équations qui décrivent les
états de contrainte, l’état de déformation et le champ de déplacement en
tout point du corps. Ces équations différentielles, accompagnées des condi-
tions aux limites, sont généralement difficiles à résoudre ; elles présentent
néanmoins un intérêt considérable pour l’ingénieur, puisqu’elles offrent -dans

1
 INTRODUCTION

le cadre des hypothèses choisies- la solution exacte de tout problème tridi-

mensionnel.
Contrairement à la mécanique des structures, la mécanique des solides ne
réduit pas un corps à une géométrie simplifiée (une poutre à son axe par
exemple) et n’introduit aucune hypothèse restrictive quant à la cinématique
de la déformation (loi de Bernoulli par exemple).

Mécanique des Solides 2 1ère année du cycle d’Ingénieurs

Chapitre 1

Notion du milieu continu

1.1 Concept du milieu continu

1.1.1 Volume élémentaire représentatif

Soit un matériau solide poreux, on nomme porosité l’ensemble des vides

(pores) dans ce solide (Fig. 1.1). Ces vides sont remplis généralement par
des fluides (liquide ou gaz). C’est une grandeur physique qui conditionne
les capacités d’écoulement et de rétention d’un substrat. Dans une portion
d’un solide poreux, la porosité peut être définie comme le rapport entre le
volume des vides et le volume total.
Vpores
ffi = (1.1)
Vtotal

Quand le volume, objet de la mesure, est très petit (V1 ), il peut être

Figure 1.1 – Définition du volume élémentaire représentatif VER à l’aide de

la porosité (d’après Kutilek et Nielsen, 1994)

entièrement situé dans un pore et la porosité est égale à 1 ; il peut aussi être

3
 CHAPITRE 1. NOTION DU MILIEU CONTINU

situé dans la phase solide et la porosité soit égale à 0. Au fur et à mesure

que le volume utilisé pour la mesure croît, la porosité des échantillons
centrés dans les pores diminue ; inversement, la porosité des échantillons
centrés dans la phase solide augmente. Ces deux porosités tendent vers une
même valeur et sont identiques à partir d’un certain volume d’échantillon.
Ce volume est appelé volume élémentaire représentatif ou tout simplement
VER.

1.1.2 Concept d’un milieu continu

Lorsque l’on ignore la structure moléculaire de la matière au profit d’une

réponse globale moyenne, alors, on se place à l’échelle macroscopique. La
matière perd alors son caractère discontinu et est admise uniformément ré-
partie dans le volume étudié. Considérons par exemple, autour d’un point
P d’un corps ˙, un volume matériel V de masse m. Diminuons V ; si la
limite définie dans l’expression (1.2) est bien définie en tout point P de ˙,
la masse est dite continûment répartie dans V .

m
ȷ = lim < +1 (1.2)
V !0 V

On peut définir semblablement des quantités comme les forces, les contraintes,

Figure 1.2 – Conception du milieu continu

les déformations, l’énergie, etc., qui sont ainsi des fonctions continues des
coordonnées. On appelle milieu continu un tel corps (solide, liquide, gaz. . .)
et la mécanique des milieux continus la science qui l’étudie.
La notion du milieu continu conduit à des théories qui fournissent des résul-
tats en excellent accord avec l’observation expérimentale. Les corps étudiés
doivent avoir des dimensions largement plus grandes que leurs constituants
élémentaires, e.g. les grains pour un matériau granulaire.
Dans ce polycopie, le concept milieu continu est adopté.

Mécanique des Solides 4 1ère année du cycle d’Ingénieurs

 CHAPITRE 1. NOTION DU MILIEU CONTINU

1.1.3 Forces

En mécanique des milieux continus, on distingue deux types de forces

extérieures (actions) :
Les forces de volume agissent sur les éléments de volume dV (ou de masse
dm) du corps, telles que les forces gravitationnelles (Fig. 1.3), électroma-
gnétique ou d’inertie. On les définit par unité de volume et on désigne par
bx ; by ; bz ou bi (i = 1; 2; 3) leurs composantes.
Les forces de surfaces ou tractions de surface sont des forces agissant sur
la surface A délimitant le corps, e.g. la pression de l’eau ou la succion du
vent par exemple. On les définit par unité d’aire (1.3) et on désigne leurs
composantes par tx ; ty ; tz ou ti (i = 1; 2; 3).

Figure 1.3 – Forces extérieures d’un milieu continu

1.1.4 Lois de conservation

Les ingrédients essentiels de la mécanique des milieux continus sont les

lois de conservation de la masse, de la quantité de mouvement, du moment
cinétique et, dans une moindre mesure, de l’énergie.
On utilise ici le système d’axes cartésien direct (x; y; z) ou (x1 ; x2 ; x3 ).

Conservation de la masse

Soit ȷ la masse volumique (ou densité de masse), V le volume et t le

temps classique. La masse totale est donnée par l’équation (1.3).

m= ȷ dV (1.3)
V

et la loi postule qu’elle reste constante (Eq. 1.4).


dm d
= ȷ dV = 0 (1.4)
dt dt V

1ère année du cycle d’Ingénieurs 5 Mécanique des Solides

 CHAPITRE 1. NOTION DU MILIEU CONTINU

1.1.5 Lois de mouvement (Euler 1776)

En appelant ~v la vitesse, ȷ~
v la densité de la quantité de mouvement et
~R la résultante des forces extérieures, la loi de la quantité de mouvement
F
est donnée par l’équation (1.5).

d
ȷ~ ~R
v dV = F (1.5)
dt V

Si ~ x ˆ ȷ~
x est un vecteur issu de l’origine (Fig. 1.4), ~ v la densité du moment
~ R le moment résultant, la loi du moment cinétique est donnée
cinétique et M
par l’équation (1.6).

d
x ˆ ȷ~
~ v dV = M ~R (1.6)
dt V

Les lois d’Euler sont parfois considérées comme une extension des lois de
Newton valables pour un ensemble de particules, via l’hypothèse du milieu
continu ; en fait, il est préférable de les accepter comme un postulat. Elles
peuvent s’énoncer comme suit :
La dérivée par rapport au temps de la quantité de mouvement et du moment
cinétique est égale, respectivement, à la force résultante et au moment
cinétique.
On peut exprimer F ~R et M~ R à l’aide des forces de volume ~b et des tractions

Figure 1.4 – Lois de mouvement

de surface ~
t (Chp. 4) ; A désigne la surface délimitant le volume V (Fig. 1.4) ;
en tenant compte des équations (1.5 et 1.6) en considérant que les limites
dV
d’intégration sont fixes. En appelant ~a = l’accélération, on obtient les
dt
équations (1.7a et 1.7b).
  
ȷ~
a dV = ~b dV + ~
t dA (1.7a)
V V  A 
x ˆ ȷ~
~ a dV = x ˆ ~b dV +
~ xˆ~
~ t dA (1.7b)
V V A

Mécanique des Solides 6 1ère année du cycle d’Ingénieurs

 CHAPITRE 1. NOTION DU MILIEU CONTINU

Énergie

L’énergie cinétique Ecin , l’énergie interne Eint , la puissance des forces

extérieures P et les puissances d’autres sources extérieures : thermiques,
électriques, ...) pouvant contribuer à l’énergie, groupées dans la notation Q,
sont liées par la loi de conservation d’énergie définie dans l’équation (1.8).
dEcin dEint
+ =P +Q (1.8)
dt dt
Dans un milieu continu, l’énergie cinétique s’exprime par l’équation (1.9).

1
Ecin = v´~
ȷ~ v dV (1.9)
V 2

et l’énergie interne peut s’écrire suivant l’équation(1.10).


Eint = e0 dV (1.10)
V

A l’aide de l’énergie interne par unité de volume e0 ; enfin, la puissance des

forces extérieures vaut, avec ~b et ~
t
 
P = ~b ´ ~
v dV + ~
t´~
v dA (1.11)
V A

Remarques

Les lois de conservation s’appliquent tant au corps complet qu’un frag-

ment arbitraire de ce corps.
Un fragment est obtenu par des coupes, auxquelles s’étendent les intégrales
d’aires (Chp. 4).

1ère année du cycle d’Ingénieurs 7 Mécanique des Solides

 CHAPITRE 1. NOTION DU MILIEU CONTINU

Mécanique des Solides 8 1ère année du cycle d’Ingénieurs

Chapitre 2

Calcul indiciel

2.1 Principe d’objectivité

Être objectif, c’est être invariant lors d’un changement de repère, être
indépendant du système d’axes ou de l’observateur. e.g. la masse, l’énergie
et la loi de comportement d’un matériau doivent être objectives.

2.2 Tenseurs et appareil mathématique

Les lois qui décrivent les phénomènes physiques doivent obéir au prin-
cipe d’objectivité. A ce point de vue, la notion symbolique (ou vectorielle)
est attirante, parce qu’elle fournit une description naturelle des phénomènes
physiques qui évite toute référence à un système de coordonnées, et souligne,
par conséquent, le caractère invariant des scalaires et des vecteurs.
Cependant, en mécanique des milieux continus on est amené à introduire des
êtres invariants, appelés tenseurs, d’un degré d’extension supérieur à celui
des vecteurs. Le calcul tensoriel édicte les règles permettant d’écrire les
équations de façon objective.
On se rattache constamment à un système d’axes cartésien direct et on tra-
vaille avec les composantes de ces nouvelles grandeurs. Ces composantes
cartésiennes ont d’ailleurs l’avantage de représenter, le plus souvent, les
grandeurs physiques réelles avec lesquelles l’ingénieur travaille concrète-
ment. De plus, on utilise la notation indicielle ; elle permet d’écrire de
façon compacte, et de manipuler avec aisance, des grandeurs et équations
qui, explicitées, s’avèrent souvent fort encombrantes.
Le tenseur est, comme le vecteur, caractérisé par ses composantes. Dans
l’espace euclidien tridimensionnel, un tenseur d’ordre n possède 3n compo-
santes. Le tenseur d’ordre 2 a neuf composantes Aij (i = 1; 2; 3; j = 1; 2; 3).
Le vecteur est un tenseur d’ordre 1, à trois composantes ak (k = 1; 2; 3), e.g.
la force, le déplacement ou l’accélération. Le scalaire est un tenseur d’ordre

9
 CHAPITRE 2. CALCUL INDICIEL

0, a une composante, e.g. la masse, la température ou l’énergie.

2.3 Calcul indiciel

2.3.1 Introduction

En calcul indiciel, on applique le calcul algébrique ordinaire. Il n’a pas

de règles spéciales (contrairement au calcul vectoriel par exemple, qui dis-
tingue le produit scalaire du produit vectoriel). Le nombre des indices montre
clairement l’ordre du tenseur. Cette simplicité a un revers : on dépend des
axes, ce qui oblige à définir la règle de transformation des composantes
d’une grandeur physique objective, lors d’un changement d’axes.
En calcul indiciel, x, y et z sont remplacés par x1 , x2 et x3 . Le système
d’axes est (x1 ; x2 ; x3 ) (Fig. 2.1). Les composantes d’un vecteur ~
a sont a1 ,
a2 et a3 ou simplement ai (Eq. 2.1).
3
X
~
a = a1~
e1 + a2~
e2 + a3~
e3 = ai~
ei (2.1)
i=1

Les indices prennent les valeurs 1 à 3 ; une grandeur telle que bi a 3 com-
posantes, Aij 9 composantes, Mijk 27 composantes, Eijkl 81 composantes,
etc.

Figure 2.1 – Repère orthonormé (O; ~

e1 ; ~
e2 ; ~
e3 ) et notation indicielle (~
ei sont
les vecteurs unitaires)

2.3.2 Règles importantes

On peut compacter l’écriture indicielle en observant quelques règles. La

première, essentielle, est la suivante :

Mécanique des Solides 10 1ère année du cycle d’Ingénieurs

 CHAPITRE 2. CALCUL INDICIEL

Règle 1 ou convention de sommation (Eïnstein)

Chaque indice apparaît deux fois dans un même monôme, ce dernier re-
présente la somme de trois thermes, en donnant aux indices les valeurs 1,2
et 3 successivement. Voici quelques exemples :

ai bi = a1 b1 + a2 b2 + a3 b3 (2.2a)
akm bm = ak1 b1 + ak2 b2 + ak3 b3 (2.2b)
gss = g11 + g22 + g33 (2.2c)

A titre d’exemple, la somme définie dans l’équation (2.2b) devrait s’écrire

comme suit
3
X
akm bm
m=1

mais, vu que les indices sont toujours attachés aux coordonnées (x; y; z),
i.e. (x1 ; x2 ; x3 ), le signe ˚ est finalement inutile.
L’indice répété s’appelle indice de sommation ; il est dit muet, puisqu’il dis-
paraît une fois la somme effectuée. Il en résulte la règle suivante :
Un même indice muet ne peut apparaître que deux fois dans un même mo-
nôme. Ainsi
cij cik xk

est correct et indique une double somme sur les indices i et k, alors que :

cmj cmm xm

est incompréhensible.

Règle 2

un indice qui n’apparaît qu’une seule fois dans un monôme doit aussi ap-
paraître une fois dans tous les monômes de l’expression.
Cet indice s’appelle indice d’identification ou indice libre et représente ordi-
nairement le numéro de l’équation traitée, e.g.

xj = aij bi + dk (2.3)

2.3.3 Symbole de Kronecker

Le symbole de Kronecker est une grandeur à deux indices notées ‹ij et

définie par l’équation (2.4).

1 si i = j
‹ij = (2.4)
0 si i 6= j

1ère année du cycle d’Ingénieurs 11 Mécanique des Solides

 CHAPITRE 2. CALCUL INDICIEL

Ce symbole est symétrique (‹ij = ‹ji ) et permet d’écrire, par exemple

ei ´ ~
~ ej = ‹ij (2.5)

En groupant les neuf valeurs dans une matrice carrée 3 ˆ 3, on constate

que le symbole de Kronecker est l’image indicielle de la matrice identité (ou
unité) I (Eq. 2.6).
0 1
B 1 0 0 C
B C
I = [‹ij ] = B
B 0 1 0 C
C (2.6)
@ A
0 0 1
Le delta de Kronecker fait apparaitre les propriétés illustrées dans les équa-
tions (2.7a, 2.7b et 2.7c).

‹kk = 3 = trI (2.7a)

ai ‹ij = aj (2.7b)
Aim ‹mj = Aij (2.7c)

L’équation (2.7a) exprime la trace de (‹kk = ‹11 + ‹22 + ‹33 = 3). Les deux
équations (2.7b) et (2.7c) traduisent la règle de substitution d’un indice ;
dans l’équation (2.8b) par exemple, ai devient aj , ou inversement (il suffit
d’effectuer la somme pour prouver). On a semblablement

‹ik ‹jk = ‹ij (2.8a)

‹ji ai = aj (2.8b)
‹im Amj = Aij (2.8c)

où (2.8a) et (2.8b) découlent la symétrie du symbole de Kronecker (i.e.

‹ij ai = ai ‹ji = ai ‹ij = aj (Eq. 2.7b)). En groupant les neufs composantes
Aij dans la matrice carrée A de dimensions 3 ˆ 3, on observe que (2.7c) et
(2.8c) expriment le produit matriciel AI=A.

2.3.4 Transformation des coordonnées

Rotation des axes

Soit (x; y; z) et (x0 ; y 0 ; z 0 ) deux systèmes d’axes cartésiens directs ayant

le même origine O (Fig. 2.2). Si par exemple, cos(x0 ; y) désigne le cosinus
directeur de l’axe x0 par rapport à l’axe y, les coordonnées x; y; z et x0 ; y 0 ; z 0
d’un point quelconque P sont liées par les équations (2.9a), (2.9b) et (2.9c)
qui représentent la transformation de coordonnées, ici une rotation, entre
(x; y; z) et (x0 ; y 0 ; z 0 ).

x0 = x cos(x0 ; x) + y cos(x0 ; y) + z cos(x0 ; z) (2.9a)

y 0 = x cos(y 0 ; x) + y cos(y 0 ; y) + z cos(y 0 ; z) (2.9b)
0 0 0 0
z = x cos(z ; x) + y cos(z ; y) + z cos(z ; z) (2.9c)

Mécanique des Solides 12 1ère année du cycle d’Ingénieurs

 CHAPITRE 2. CALCUL INDICIEL

Introduisons la notation indicielle ; on désigne par cij le cosinus directeur

Figure 2.2 – Transformation de coordonnées (rotations)

de l’axe x0i par rapport à l’axe xi (Eq. 2.10).

cij = cos(x0i ; xj ) = e
~0 i ´ ~
ej (2.10)

Les relations (2.9a), (2.9b) et (2.9c) peuvent s’écrire selon l’équation (2.11a),
alors que la transformation inverse (passage des x0i au xj ) est exprimée par
l’équation (2.11b).

x0i = cij xj (2.11a)

xj = cij x0i (2.11b)

Matriciellement, ces équations s’écrivent conformément aux équation (2.12a)

et (2.12b).
0
X = CX (2.12a)
T 0
X=C X (2.12b)

où C = [cij ] est la matrice de rotation, composée des cosinus directeurs, et

où l’indice T désigne la transposition. On sait que la matrice de rotation est
orthogonale, donc
cik cjk = ‹ij (2.13)
Matriciellement, l’équation (2.13) peut s’écrire comme suit

CCT = I (2.14)

On note aussi que le déterminant d’une matrice orthogonale est égal à l’unité
(Eq. 2.15).
det C = 1 (2.15)

1ère année du cycle d’Ingénieurs 13 Mécanique des Solides

 CHAPITRE 2. CALCUL INDICIEL

Transformation des coordonnées

Le scalaire, ou le tenseur d’ordre 0, n’a qu’une composante, invariante

dans tous système d’axes. Pour un vecteur, ou tenseur d’ordre 1, les compo-
santes se transforment de manière identique aux équations (2.11a) et (2.11b),
il suffit de considérer le vecteur ~ ~ . On a
v = OP

vi0 = cir vr (2.16)

La règle est : si les composantes d’un vecteur se transforment selon l’équa-

tion (2.16) dans un changement d’axes, la grandeur est un tenseur d’ordre
1, soit un vecteur.
La règle de transformation des composantes d’un tenseur d’ordre 2 s’obtient
alors par une simple généralisation (Eq. 2.17).

0
Tij = cir cjs Trs (2.17)

et s’énonce, si les composantes Trs d’une grandeur se transforment selon

l’équation (Eq. 2.17) dans le changement d’axes, la grandeur est un tenseur
d’ordre 2.
Matriciellemnt, l’équation (2.17) peut s’écrire

T = CTCT (2.18)

où T et T0 sont des matrices carrées de dimensions 3 ˆ 3 dont les élé-

ments sont les composantes des tenseurs. C’est la manière la plus simple
de représenter ces neuf composantes (Eq. 2.19).
0 1
B T11 T12 T13 C
B C
T = B
B T21 T22 T23 C
C (2.19)
@ A
T31 T32 T33

En généralisant les équations (2.16) et (2.17), on définit la règle de trans-

formation des composantes d’un tenseur d’ordre n (Eq. 2.20).

0
Tijk:::n = cir cjs ckt :::cnw Trst:::w (2.20)

2.4 Exercices

Exercice 2.1 Les expressions indicielles suivantes sont-elles correctes ?

1. aij xj + bi = 0
2. aij cjj = ei
2
3. 52 u = ‹mn @x@m @x
u
n


4. W = Tij deij

Mécanique des Solides 14 1ère année du cycle d’Ingénieurs

 CHAPITRE 2. CALCUL INDICIEL

@xi @xj
5. eij = L
@xk kj
+ L
@xk ki
1
6. Eij = (Jki Jkj ` ‹ij ) + Gmij
2
7. f = aik (bjk + cij )dij
8. dai = J @X m
@xi
dAm
p
9. ds = dxi dxi

Exercice 2.2 Les expressions indicielles suivantes sont-elles correctes ?

1. wi = kemni wmn
2. eij = Askk ‹ij + Bsik
@ 2 Tij
3. @xi @xj
= R▽2 Tkk
4. eijk eirs = ‹jr ‹ks ` ‹js ‹kr
5. bsn wn = Aekmn brk bsm btn ur us
@Xi
6. N
@xr i
dS = krst ds dxt
7. Tkm = k h k @Xk
T hm @Xm
hr @xr rs hs @xs
n
X @fp @fp
8. eij = Cp Tkm
p=1 @Tij @Tkm

Exercice 2.3 On donne xi = aij yj , mais yi = bij zj . Éliminer les y pour

avoir x en fonction de z.

Exercice 2.4 Développer entièrement l’expression ci-après et donner l’ex-

pression matricielle.
a = Dij xi xj

Exercice 2.5 Ma gomme a la forme d’un parallélépipède rectangle de di-

mensions :
l1 ˆ l2 ˆ l3 = 50mm ˆ 20mm ˆ 10mm
Elle me sert souvent pour illustrer les notions de base de la mécanique des
milieux continus.
J’ai posé ma gomme debout sur la table, i.e. en contact avec une face l2 ˆl3 ,
et je l’ai comprimé en appliquant une force de 400N sur la face opposée.
Ses nouvelles dimensions h1 ˆ h2 ˆ h3 sont telles que h1 = 45mm et
h2 = 21mm.
1. A quoi est égal h3 ? Calculer le module d’Young E et le coefficient de
Poisson ‌.
2. J’ai ensuite pris ma gomme par ces mêmes faces et je l’ai étiré avec une
force de 200N de chaque côté. Quelles sont ses nouvelles dimensions ?
Finalement, j’ai plongé dans l’océan à 400m de fond avec ma gomme.
Quelles sont ses nouvelles dimensions ?

1ère année du cycle d’Ingénieurs 15 Mécanique des Solides

 CHAPITRE 2. CALCUL INDICIEL

Mécanique des Solides 16 1ère année du cycle d’Ingénieurs

Chapitre 3

Champ gradient de
déformation et ses
applications

3.1 Champ gradient de déformation

Le champ gradient d’une déformation associe à chaque point de la confi-

guration de référence une matrice 3 ˆ 3 qui transforme un petit vecteur
(infinitésimal) du voisinage proche de ce point en un petit vecteur du voisi-
nage du point image de la configuration déformée. Cette matrice contient les
composantes d’un tenseur qui obéissent aux règles de changement de base.
A partir de ce tenseur, on calcule le tenseur de dilatation qui permet de
comparer localement les longueurs et les volumes de ces petits voisinages
lors de la déformation.

3.1.1 Déformation d’un milieu continu

Les milieux continus qui seront étudiés se déforment et se déplacent

dans un espace euclidien V que l’on identifie à R3 par le choix d’une base
“ ”
orthonormée e(1) ; e(2) ; e(3) , on considère un milieu formé des particules
matérielles a d’un ensemble fermé ˙ 2 R3 , que l’on appellera configuration
de référence. Une déformation du milieu est définie pa l’application X qui
associe chaque point a une nouvelle position x = X(a). Cette application
est définie par ce qui suit
X : ˙ `! R3
(3.1)
a 7`! X(a)
A côté des notations a, x ou X, qui désignent les matrices colonnes de
coordonnées dans R3 , nous utiliserons la notation indicielle (Eq. 3.2).
xi = Xi (a) = Xi (a1 ; a2 ; a3 ) avec : i 2 f1; 2; 3g (3.2)

17
 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

Figure 3.1 – Déformation d’un milieu continu

On note ˙ l’image de ˙0 par l’application X et on l’appelle configuration

déformée ou configuration actuelle.

3.1.2 Tenseur gradient de la déformation

Une déformation quelconque va transformer un réseau de droites en un ré-

seau de courbes de forme quelconque. Néanmoins, on peut approcher l’image
de très petits segments de droites (infinitésimaux). Ainsi, dans les petits voi-
sinages d’un point a et de son image a = X(a), on peut approcher la
déformation X par une application linéaire associant à des petits vecteurs du
voisinage de a des petits vecteurs au voisinage de x (Fig. 3.1). D’un point de
vue mathématique, ces petits vecteurs sont en fait des vecteurs infinitésimaux
et l’application linéaire qui les relie est la jacobienne de la déformation au
point a. Néanmoins, on va continuer à utiliser le terme de petit vecteur pour
simplifier.
Explorons alors les images x0 = X(a0 ) des points a0 situés dans un petit
voisinage de a (Fig. 3.2). Soit da = a0 ` a petit vecteur traduisant cette
exploration autour de a et dx = x0 ` x = X(a0 ) ` X(a) le petit vecteur
transporté par la déformation.
On suppose que l’application X est différentiable (Eq. 3.3).

x0 = X(a0 ) = X(a + da) = X(a) + F (a) da + O[(da)2 ] (3.3)

Ce développement limité fait apparaitre une matrice 3 ˆ 3 que l’on note

Figure 3.2 – Linéarisation de l’application X au voisinage du point a

Mécanique des Solides 18 1ère année du cycle d’Ingénieurs

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

F (a) dont les composantes sont définies dans l’équation indicielle (Eq. 3.4).
@Xi
Fij (a) = (a) (3.4)
@aj
Cette matrice est souvent appelée jacobienne de l’application X au point a
ou différentielle de X au point a, elle sera appelée ici gradient de l’appli-
cation X au point a (Eq. 3.5).
F (a) = grad X(a) (3.5)
De même qu’un vecteur a pour composantes une matrice colonne, on peut
considérer que la matrice jacobienne est un tenseur d’ordre 2. Ceci signifie
que cette matrice définit un objet mathématique intrinsèque indépendant
de la base. Lors d’un changement de base, les composantes de ce tenseur
doivent changer suivant les règles précitées.
Le gradient de la déformation est donc un champ de tenseurs défini par la
relation (3.6).
dx = F (a) da + O[(da)2 ] (3.6)
En notation indicielle

Figure 3.3 – Exploration des voisinages de a et x

dxi = Fij (a) + O[(da)2 ] (3.7)

Où nous adoptons la convention de sommation d’Einsteïn (§2.3.2). D’autre
part, nous omettrons désormais de préciser le terme O[(da)2 ] dans les ex-
pressions, ce qui revient à considérer les petits vecteurs comme des vecteurs
infinitésimaux.

3.2 Champ de tenseur de dilatation

3.2.1 Transport du produit scalaire

Pour étudier les propriétés physiques d’un milieu continu il est important
de mesurer comment varient les longueurs et les volumes lors d’une défor-
mation. Nous allons voir que cette information est contenue dans le tenseur

1ère année du cycle d’Ingénieurs 19 Mécanique des Solides

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

des dilatations. On considère donc une déformation X qui transforme la

configuration de référence˙0 en une configuration ˙.

dx = x0 ` x = F (a) da (3.8a)
dx0 = x00 ` x = F (a) da0 (3.8b)

On se place au voisinage de a 2 ˙0 et l’on considère deux petit vecteurs

da = a0 ` a et da0 = a00 ` a reliant a aux deux points voisins a0 et a00
(Fig. 3.4). Soitx, x0 et x00 les images de a, a0 et a00 par l’application X.
On note dx et dx0 les vecteurs transportés par la déformation dont les
composantes vérifient les équations (3.8a) et (3.8b). On s’intéresse à la

Figure 3.4 – Transport du produit scalaire au voisinage de a

fonction C(a; da; da0 ) qui associe aux petits vecteurs da et da0 , au voisinage
de a, le produit scalaire de leurs images (Eq. 3.9).

C(a; da; da0 ) = dx ´ dx0 = dxT dx0 = dxi dx0i (3.9)

où "´" désigne le produit scalaire de deux vecteurs et dxT la matrice ligne

transposée de la matrice colonne dx. Ce produit scalaire est transporté par
la déformation (Eq.3.10).

C(a; da; da0 ) = dx ´ dx0 = (F da) ´ (F da0 )

= (F da)T ´ (F da0 )
(3.10)
= daT F T F da0
= dai Fni Fnj da0j

3.2.2 Tenseur de dilatation

Le champ de tenseur symétrique C (Eq. 3.11a et 3.11b) est appelé ten-

seur des dilatations de Cauchy, ou encore tenseur de dilations de Cauchy-
Green. On peut le voir comme définissant la forme bilinéaire symétrique
C(a; da; da0 ) dans l’espace des vecteurs da du voisinage de a, cette forme

Mécanique des Solides 20 1ère année du cycle d’Ingénieurs

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

étant caractérisée par la relation (Eq. 3.11c).

C(a) = F T (a)F (a) (3.11a)

Cij (a) = Cji (a) = Fni (a)Fnj (a) (3.11b)
0 0 T
dx ´ dx = C(a; da; da ) = da C(a)da T
= dai Cij da0j (3.11c)

Cette forme bilinéaire symétrique définit un nouveau produit scalaire qui

indique comment sont déformés les longueurs et les angles.

3.2.3 Dilatation relative des longueurs

Le tenseur des dilatations C(a) permet de mesurer ˜(a; da), la dilatation

relative des longueurs dans la direction da et au voisinage du point a. Cette
dilatation est définie par la relation (3.12).
q
kdxk C(a; da; da)
˜(a; da) = = (3.12)
kdak kdak

où kdxk est la norme de l’image dx du petit vecteur da pris dans le voisi-

nage de a. La dilatation relative des longueurs est donc une grandeur sans
dimensions, car elle exprime une longueur dilatée par unité de longueur. En
utilisant l’expression de C en fonction de C, on peut écrire

v v
T
u
u da C(a)da u dai Cij daj
u u
˜(a; da) = t = t
(3.13)
kdak2 dan dan

La fonction ˜ est donc le rapport de deux normes, la norme définie par le

produit scalaire C et la norme canonique.

3.2.4 Glissement des vecteurs orthogonaux

Le tenseur de dilatation C permet aussi de calculer l’angle de glisse-

ment ‚(a; da; da0 ) de deux petits vecteurs orthogonaux da et da0 défini par
l’équation (3.13).

dx ´ dx0 = kdxkkdx0 k sin ‚(a; da; da0 ) (3.14)

avec da´da0 = 0 Cette définition entraîne que les vecteurs déformés dx et dx0
ı
forment un angle „ = `‚. L’angle de glissement ‚ est encore appelé angle
2
de cisaillement. Il mesure la variation d’angle de deux vecteurs orthogonaux,
après la déformation. On remarque que l’on peut exprimer ‚ en fonction de

1ère année du cycle d’Ingénieurs 21 Mécanique des Solides

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

Figure 3.5 – Angle de glissement des petits vecteurs da et da0

C à l’aide de la relation (3.15).

C(a; da; da0 )
sin ‚(a; da; da0 ) = q q
C(a; da; da C(a; da0 ; da0
daT C(a)da0
= r r (3.15)
daT C(a)da da0T C(a)da0
dai Cij da0j
= p q
dam Cmn dan da0p Cpq da0q
avec da ´ da0 = da ´ da0 = 0

3.2.5 Transport de la base canonique

Pour chaque point a, les composantes diagonales C11 , C22 et C33 du ten-
seur des dilatations C(a) peuvent s’interpréter comme le carré des dilatations
(ou contractions) relatives des vecteurs de base. Les autres composantes C12 ,
C13 et C23 , sont reliées au glissement des couples des vecteurs de base.
En effet, soit da(1) = da e(1) , da(2) = da e(2) et da(3) = da e(3) trois pe-
tits vecteurs dans le voisinage a et proportionnels aux vecteurs de la base
canonique e(1) ; e(2) ; e(3) . Soit dx(1) , dx(2) et dx(3) les trois petits vecteurs
images transportés par la déformation, i.e. les petits vecteurs du voisinage
de x = X(a) qui vérifient l’équation (3.16).

dx(i) = F (a) da(i) (3.16)

Les coefficients diagonaux C11 , C22 et C33 sont égaux aux carrés des dila-
tations relatives dans les trois directions de la base canonique, au voisinage
de a, comme le montre la relation (3.17).

kdx(1) k2 (da e(1) )T C(da e(1) )

2 (1)
˜ [a; da e ]= = = C11 (3.17)
kda(1) k2 (da)2
ainsi que les relations similaires obtenues pour les indices 2 et 3.
Les coefficients non diagonaux Cij avec i 6= j vérifient des relations déduites
de l’exemple définit par l’équation (3.18).
T
dx(1) ´dx(2) = kdx(1) kkdx(2) k sin ‚12 = (da)2 e(1) C e(2) = (da)2 C12 (3.18)

Mécanique des Solides 22 1ère année du cycle d’Ingénieurs

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

Figure 3.6 – Transport de la base canonique au voisinage de a

où ‚ij = ‚[a; da e(i) ; da e(j) ] est le glissement du couple de vecteurs ortho-

gonaux [da e(i) ; da e(j) ], i.e. la variation d’angle qu’ils subissent lors de la
déformation. On peut alors écrire
C12
sin ‚12 = p (3.19)
C11 C22
ainsi que les relations similaires pour les couples 13 et 23.

3.2.6 Transport des volumes

Rappelons que le volume d’un petit parallélépipède engendré par trois

vecteurs quelconques da, da0 et da00 est égal à la valeur absolue de leur
produit mixte (3.20).
dV0 = j(da; da0 ; da00 )j (3.20)

Le produit mixte est défini par la relation (3.21).

Figure 3.7 – Transport d’un volume infinitésimal parallélépipédique

(da; da0 ; da00 ) = da ´ (da0 ˆ da00 ) = "ijk dai da0j da00k (3.21)

où "ijk est le pseudo-tenseur fondamental alterné qui vaut 0 si deux indices

sont égaux, 1 si les indices sont une permutation circulaire de f1; 2; 3g et `1
sinon ("123 = 1 mais "213 = `1, etc.). Un pseudo-tenseur diffère d’un tenseur
dans la mesure où le signe de ses composantes dépend de l’orientation de
l’espace. Si w = u ˆ v est le produit vectoriel de deux vecteurs u et v
quelconques, ses composantes s’écrivent wi = "ijk uj uk . Le produit mixte

1ère année du cycle d’Ingénieurs 23 Mécanique des Solides

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

n’est autre que le déterminant de la matrice carrée formée des matrices

colonnes da; da0 ; da00 . Après déformation, le volume élémentaire devient

dV =j(dx; dx0 ; dx00 )j

=j(F da; F da0 ; F da00 )j
(3.22)
=j det F jj(da; da0 ; da00 )j
=jdetF j dV0

En notant J(a) =j det F (a)j le Jacobien de la déformation X, on obtient la

relation (3.23).
dV
= J(a) (3.23)
dV0
Le Jacobien J(a) exprime donc la dilatation volumique relative induite par
la déformation. Cette relation s’étend à un petit volume dV0 de forme quel-
conque, pris dans un voisinage de a, qui est transformé par la déformation
X en petit volume dV .

Exercice 3.1 Ma pièce métallique a la forme d’un cube de côté l1 = l2 =

l3 = 20mm. Le choix des axes est tel que a 2 [0; 1]3 si a est un point de
la pièce. Je l’ai tordue à l’aide d’un étau et une clé à fourche. Sa grande
déformation x = X(a) est modélisée par les fonctions suivantes :

x1 = ka1
x2 = a2
x3 = a3 + ˛a21

avec ˛ = 0; 1mm`1 et k = 2.
Avant la déformation, j’ai dessiné sur une face l1 ˆ l1 de ma pièce une petite
croix composée de quatre petits segments de longueur 0; 1mm ayant une
extrémité au centre de la face et dirigés dans les directions 1 et 3.
1. Calculer sans aucune approximation les images des quatre vecteurs da
composant la croix. Comparer avec les quatre vecteurs dx obtenus en
considérant les da comme des petits vecteurs (infinitésimaux). Vérifier
que l’erreur est bien d’ordre O[(da)2 ].
2. Calculer de manière exacte les quatre longueurs des segments de la
croix après déformation. Calculer le tenseur des dilatations ainsi que
les longueurs des quatre petits segments déformés en les considérant
comme petits (infinitésimaux). Comparer avec le calcul exact.
3. Quels angles font entre eux les petits segments de la croix après le
déformation ?
Avant la déformation, j’ai percé un petit trou au centre de la croix en
enfonçant de 2mm la tête d’un pointeau de diamètre 0; 5mm.
4. Que devient le volume de ce trou après la déformation ?

Mécanique des Solides 24 1ère année du cycle d’Ingénieurs

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

5. Quelle forme a la pièce après la déformation ? (faire un tracé)

Exercice 3.2 On considère un solide élastique homogène et isotrope ayant

la forme d’un cylindre de révolution de rayon R et de longueur l. On choisi
l’axe (Oa3 ) sur l’axe du cylindre. On suppose que l’on exerce das forces sur
ce solide de nature à le déformer, et l’on mesure le champs de déplacement
correspondant. Dans cette expérience, la déformation x = X(a) du solide

Figure 3.8 – Solide cylindrique élastique

est décrite par

x1 = a1 + ka1 a2
“ ”
x2 = a2 + k a21 + a22
x3 = a3
1. Calculer le tenseur de dilatation C(a) en un point a quelconque du
solide.
2. Interpréter les composantes diagonales de ce tenseur en termes de di-
latation de petits vecteurs.
3. Interpréter les composantes non diagonales en termes de variation d’angle.
4. Calculer la dilatation relative des longueurs dans la direction (1; 1; 1)
pour des fibres prises autour d’un point quelconque a.

Exercice 3.3 On considère un parallélépipède rectangle de longueur l1 , l2

et l3 , composé d’un matériau élastique homogène et isotrope. On suppose
qu’il est non contraint et que sa densité et ȷ0 . On choisit les axes orthonormés
Oa1 , Oa2 et Oa3 parallèles aux arêtes.
On effectue une expérience de traction sur le solide en appliquant des forces
de contact à sa surface. On suppose que la déformation x = X(a) s’écrit

x1 = a1 + ka1
x2 = a2 ` ‌ka2
x3 = a3 ` ‌ka3

1ère année du cycle d’Ingénieurs 25 Mécanique des Solides

 CHAPITRE 3. CHAMP GRADIENT DE DÉFORMATION ET SES APPLICATIONS

1. Calculer le tenseur des dilatations C(a) et interpréter ses composantes.

2. Lors de cette expérience, on mesure le champ de température dont la
représentation eulérienne est

T (E) (x) = k1 x1 + k2 x2 + k3 x3

Calculer sa représentation lagrangienne T (L) (a).

3. Calculer les représentations eulérienne et lagrangienne du champ de
densité ȷ obtenu après le déformation.
4. Interpréter physiquement cette expérience.

Mécanique des Solides 26 1ère année du cycle d’Ingénieurs

Chapitre 4

Tenseur des contraintes et

équations d’équilibre

4.1 Milieu déformable

4.1.1 Forces externes et équilibre mécanique

La loi fondamentale de la dynamique appliquée à un corps en équilibre

mécanique statique implique la nullité de la résultante et du moment résul-
tant des efforts externes appliqués (Eq. 4.1a et 4.1b).

X
F =0 (4.1a)
X
M =0 (4.1b)

a) Si le nombre d’efforts externes est égal au nombre d’équations, le corps

est en état d’équilibre isostatique et les efforts sont entièrement détermi-
nés.

b) Si ce nombre d’efforts est supérieur au nombre d’équations, le corps est

en état d’équilibre hyperstatique et la détermination des efforts implique
la prise en compte de la déformation du milieu pour obtenir les équations
complémentaires.

4.1.2 Etat mécanique local

Pour caractériser l’état de sollicitation mécanique (efforts internes) en

tout point du matériau indépendamment de la géométrie du corps, la des-
cription en termes de force est insuffisante, il faut introduire la notion de
densité surfacique de force, appelée encore traction de surface.

27
 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

Contraintes

Elles caractérisent l’état de sollicitation mécanique (efforts internes) en

tout point du matériau, indépendamment de sa nature et de ses caractéris-
tiques de résistance.

Traction de surface

La figure (4.1) illustre un principe extrêmement important appelé principe

de la coupe 1 . Si la force de contact résultante f représentant l’action exercée
par la partie A sur la section de coupe de la partie B doit bien équilibrer la
force externe F , soit f + F = 0. Cette équation n’a pas de raison, sauf le cas
particulier où la distribution de l’effort est uniforme (en module et/ou en
direction) sur cette section. Chaque élément d’aire dA supportera une force


élémentaire df = t dA telle que A t dA = f. La densité surfacique de force
t, appelé vecteur contrainte ou traction de surface, est une fonction du point
courant P de la section de coupe. Plutôt que d’employer le terme section
de coupe, nous utiliserons désormais le terme de facette plus adapté au fait
que les coupes sont virtuelles. Étant un point P 0 donné, considérons deux

Figure 4.1 – Le vecteur contrainte t dépend du point P et de l’orientation

n de la facette

sections de coupe A et A0 caractérisé respectivement par leurs orientations n

et n0 . Pour chacune s’elles le vecteur contrainte t vérifie les deux équations
(4.2a) et (4.2b).

t dA = f sur A (4.2a)
A
t dA0 = f sur A0 (4.2b)
A
Au point P0 : t(P0 ; n) 6= t(P0 ; n0 )
Le vecteur contrainte t est également une fonction de l’orientation n de la
facette.
1. Le principe d’équilibre est énoncé comme suit : "Tout fragment d’un solide soumis à
un système de forces en équilibre est lui-même soumis à un système de force en équilibre,
différent du premier".

Mécanique des Solides 28 1ère année du cycle d’Ingénieurs

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

4.1.3 Repère local, contraintes normales et de tangentielles

Au point P , pour une facette de direction n, on définit le repère local

(trièdre direct) (P; n; r; w) de telle manière que les vecteurs t, n et r soient
coplanaires (Fig. 4.2).
Dans ce repère les composantes du vecteur contrainte t(P; n) sont ffnn ,ffrn
et ffwn , le premier indice faisant référence à l’axe de projection, le second à
l’orientation n du plan de la facette. Par construction on a toujours ffwn = 0
et les composantes non nulles de t sont sa composante normale ffnn selon
n appelée contrainte normale notée ff et sa composante tangentielle ffrn
située dans le plan de la facette perpendiculaire à n appelée contrainte de
cisaillement notée fi .

Figure 4.2 – Le repère local de la facette : P; n; r; w

Contrainte normale ff

C’est la composante de t qui tire lorsqu’ell est parallèle à n. Le signe de

ff est donc intrinsèque, ff > 0 pour la traction et ff < 0 pour la compression.

Contrainte de cisaillement fi

Appelée aussi contrainte tangentielle. C’est la composante de t qui agit

à la manière d’une cisaille, d’où,son nom.
La dimension des composantes de contrainte est [ml`1 t`2 ], rapport de la
force [mlt`2 ] à l’aire [l2 ] de la section sur laquelle elle s’exerce.. Son unité
est le Pascal (ou Nm`2 ) . En pratique, cette unité est trop petite et l’unité
courante du mécanicien est le Mégapascal (MP a ou Nmm2 ), voire le Gi-
gapascal (GP a ou GNm`2 ).
L’état mécanique local au point P , parmi toutes les directions de facettes
possibles, certaines supporteront un cisaillement maximal tandis que d’autres
supportent une traction maximale. Il est important de faire une distinction
entre les deux types de composantes normale et tangentielle car l’expérience

1ère année du cycle d’Ingénieurs 29 Mécanique des Solides

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

montre que tout matériau présente, selon sa nature, un seuil caractéristique

de contrainte ffC ou fiC selon son mode de ruine, seul au-delà duquel il cède.
- Les matériaux ductiles -tels les métaux- cèdent par déformation plas-
tique lorsque le cisaillement maximal fiM atteint le seuil fiC
- Les matériaux dits fragiles -tels le vers- cèdent lorsque la traction
maximale ffM atteint le seuil ffC .

4.1.4 Tenseur de contraintes

Équilibre des forces

L’équilibre d’un tétraèdre qui doit s’écrire en termes de force et non en

termes de contrainte (densité surfacique de force) implique la relation (4.3).

t(P; n)AABC + t1 (P; `e1 )AP BC + t2 (P; `e2 )AP AC + t3 (P; `e3 )AP AB = 0
(4.3)
Comme, en vertu du principe de l’action et de réaction, t(P; `ei ) = `t(P; ei ),
cette relation s’écrit en tenant compte des relations (4.3).

t(P; n) = t1 (P; e1 )n1 + t2 (P; e2 )n2 + t3 (P; e3 )n3 (4.4)

En introduisant les composantes ff1i ,ff2i ,ff3i des vecteurs contraintes ti (P; ei )

Figure 4.3 – Équilibre d’un tétraèdre

avec i 2 1; 2; 3 représentant les actions des faces P BC, P AC et MAB sur

Mécanique des Solides 30 1ère année du cycle d’Ingénieurs

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

le reste du corps, on trouve

t1 (P; e1 ) = ff11 e1 + ff21 e2 + ff31 e3

t2 (P; e2 ) = ff12 e1 + ff22 e2 + ff32 e3 (4.5)
t3 (P; e3 ) = ff13 e1 + ff23 e2 + ff33 e3

Finalement la relation entre t(P; n) et n s’écrit de manière générale comme

suit
t(P; n) = ff n (4.6)
A l’aide du tenseur ff(P ) dont les composantes dans le repère (P; e1 ; e2 ; e3 )
sont dans l’équation (4.7).
0 1
B ff11 ff12 ff13 C
B C
ff = B
B ff21 ff22 ff23 C
C (4.7)
@ A
ff31 ff32 ff33

La connaissance des vecteurs contraintes agissant sur trois facettes ortho-

gonales de normales e1 , e2 et e3 suffit à déterminer au point P l’état
de contrainte caractérisé par le tenseur des contraintes ff, donc le vecteur
contrainte t(P; n) sur toute facette contenant P .

Équilibre des moments

Cette fois, on considère une particule parallélépipédique élémentaire de

volume dx1 dx2 dx3 . L’équilibre de cette particule impose que moment résul-
tant soit nul. On va calculer le moment par rapport à l’axe x3 .

ff12 (dx1 dx3 )dx1 ` ff21 (dx1 dx3 )dx2 = 0 (4.8a)

ff12 = ff21 (4.8b)

De manière analogue on obtient avec les deux autres composantes du mo-

ment (Eq. 4.9a et 4.9b).

ff23 = ff32 (4.9a)

ff13 = ff31 (4.9b)

La symétrie du tenseur des contraintes (Eq. 4.10a et 4.10b).

ff = ff T (4.10a)
ffij = ffji (4.10b)

est une conséquence directe de la nullité du moment résultant.

Les composantes de ff forment deux groupes distincts :
- Les composantes de la diagonale principale ffij avec i = j compo-
santes normales (traction > 0 ou compression < 0) dirigés selon les
normales aux facettes.

1ère année du cycle d’Ingénieurs 31 Mécanique des Solides

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

- Les composantes hors diagonales ffij avec i 6= j composantes tan-

gentielles (contraintes de cisaillement) contenues dans le plan des
facettes.

4.1.5 Propriétés du tenseur contraintes

Le tenseur ff(P ) définissant l’état de contrainte en un point P est un

opérateur intrinsèque qui lie la normale n à une facette passant par ce
point, au vecteur contrainte qui s’applique à cette facette (Eq. 4.11).

t = t(P; n) (4.11)

Cette définition est indépendante du repère choisi.

Comme pour tout tenseur symétrique d’ordre 2, il existe toujours un repère
particulier, le repère principal 2 dans lequel le tenseur des contraintes ff(P ) a
pour représentation une matrice diagonale dont les composantes sont appe-
lées contraintes principales ff1 , ff2 et ff3 . Ce sont les valeurs propres toujours
réelles de l’opérateur. Notons que le tenseur ff est toujours défini par six
quantités indépendantes, ses six composantes dans un repère quelconque ou
ses trois contraintes principales auxquelles il faut associer les trois angles
d’Euler définissant la position du repère principal par rapport à celui du la-
boratoire.
Les contraintes principales ff1 , ff2 et ff3 , ainsi que le tenseur ff lui-même,
vérifient les équations (4.12a), (4.12b) et (4.12c).

det(ff ` –I) (4.12a)

–3 ` I1 –2 + I2 – ` I3 = 0 (4.12b)
ff 3 ` I1 ff 2 + I2 ff ` I3 I = 0 (4.12c)

Ces équations sont indépendantes du repère choisi de sorte que les trois co-
efficients I1 , I2 et I3 sont des combinaisons des composantes de ff invariantes
par changement de repère. Ces coefficients représentent le jeu d’invariants
du tenseur : sa trace (4.13a), la somme des déterminants mineurs des com-
posantes (4.13b) et son déterminant (4.13c).

I1 = trff = ffii = ff1 + ff2 + ff3 (4.13a)

2 2 2
I2 = ff11 ff22 ` ff12 + ff22 ff33 ` ff32 + ff33 ff11 ` ff31 = ff1 ff2 + ff2 ff3 + ff3 ff1
(4.13b)
I3 = det ff = ff1 ff2 ff3 (4.13c)

2. Soit A et B deux points arbitraires d’un solide. Si l’on constate, après une déformation,
que c’est uniquement la distance AB qui a changé sans aucun changement d’angle, alors les
points A et B sont situés dans une direction principale. Pratiquement, lorsqu’on veut mesurer
les contraintes, les jauges d’extensométrie doivent être collées dans la direction principale.

Mécanique des Solides 32 1ère année du cycle d’Ingénieurs

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

Il existe une infinité de jeux d’invariants combinaisons du jeu précédents. Les

combinaisons la plus utilisées sont définies dans les équations (4.14).

I2˜ = trff 2 ou I2˜ = I12 ` 2I2 (4.14a)

I3˜ = trff 3 ou I3˜ = I13 ` 3I1 I2 + 3I3 (4.14b)

4.2 Équations de mouvement

4.2.1 Parallélépipède élémentaire

Pour fixer les idées,on considère une portion d’un solide de forme parallé-
lépipédique rectangle de côtés dxi parallèles aux axes xi (Fig. 4.4). Sur ce
corps isolé agissent des tractions de surface sur les six faces et deux effets
de volume. La figure ne montre que les composantes parallèles à l’axe x1 ,
par simplicité. On remarque que, sur deux faces parallèles, les composantes
des contraintes diffèrent de l’accroissement correspondant. Prenons la com-
posante ff11 par exemple ; en x1 , x2 , x3 et indépendamment du temps, on
a ff11 (x1 ; x2 ; x3 ) ; un rien plus loin sur l’axe x1 , soit en x1 + dx1 ; x2 ; x3 ,
la contrainte vaut ff11 (x1 + dx1 ; x2 ; x3 ). La formule de Taylor limitée aux
premiers termes donne

@ff11 (x1 ; x2 ; x3 )
ff11 (x1 + dx1 ; x2 ; x3 ) = ff11 (x1 ; x2 ; x3 ) + + O(dx21 ) (4.15)
@x1

@ff11
Dans la figure, on écrit brièvement ff11 et ff11 + dx1 sur facettes x1
@x1
négatives et positives. On applique sur cette portion du solide les lois d’Euler
(Eq 1.7a et 1.7b).

4.2.2 Équations de mouvement

La loi de la quantité de mouvement selon x1 s’écrit sous forme de l’équa-

tion (4.16), qui donne après simplification l’équation (4.17a). On trouve deux
équations semblables selon x2 et x3 . L’ensemble de ces équations aux déri-

1ère année du cycle d’Ingénieurs 33 Mécanique des Solides

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

Figure 4.4 – Portion parallélépipédique élémentaire d’un solide

vées partielles s’écrit sous forme de l’équation (4.17b) 3 .

@ff11
`ff11 dx2 dx3 + (ff11 + dx1 )dx2 dx3
@x1
@ff21
`ff21 dx1 dx3 + (ff21 + dx2 )dx1 dx3
@x1 (4.16)
@ff31
`ff31 dx1 dx2 + (ff31 + dx3 )dx1 dx2
@x3
+b1 dx1 dx2 dx3 = ȷa1 dx1 dx2dx3

@ff11 @ff21 @ff31

+ + + b1 = ȷa1 (4.17a)
@x1 @x2 @x3
@ffji
+ bi = ȷai (4.17b)
@xj
div ff + b = ȷa (4.17c)

On note que les composantes du tenseur contrainte sont des fonctions sont
des fonctions des coordonnées et du temps (ffij (x1 ; x2 ; x3 ; t)), dérivables et
continues (de classe C 1 ).

3. Les équations de mouvement @ffji;j +bi = ȷai sont appelés les équations de mouvement
de Cauchy. On signale aussi que les la contrainte ff et la déformation " sont appelés les
variables de Cauchy. Par ailleurs, on verra en deuxième année que le moment M et la courbure
» (l’inverse du rayon de courbure) sont appelés, respectivement, contrainte généralisée et
déformation généralisées.

Mécanique des Solides 34 1ère année du cycle d’Ingénieurs

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

4.2.3 Symétrie du tenseur contrainte et principe de réciprocité

des ffij

Pour exprimer la loi du moment cinétique, il est avantageux de prendre les

axes par le centre du parallélépipède élémentaire. La figure (4.5) montre les
seules forces qui interviennent quand on considère l’axe x1 . En effet, dans


(1.7b), le dernier terme A ~ xˆ~ tdA correspond aux couples des forces de
surface. Les forces ff23 dx1 dx3 et ff32 dx1 dx2 donnent autour de x1 l’équation

(ff23 dx1 dx3 )dx2 ` (ff32 dx1 dx2 )dx3 = (ff23 ` ff32 )dx1 dx2 dx3 (4.18)

La loi du moment cinétique se réduit donc, selon l’axe x1 à

(ff23 ` ff32 )dx1 dx2 dx3 = 0 (4.19)

d’où
ff23 = ff32 (4.20)

On procède semblablement selon les deux autres axes et on obtient l’en-

semble des trois équations (Eq. 4.10b), parfois appelés secondes équations
de mouvement de Cauchy. Ceci prouve que le tenseur des contraintes est
symétrique. Par suite, l’état de contrainte en un point d’un solide ou d’un
milieu continu ne dépend plus que de six paramètres. D’où la dénomina-
tion "Principe de réciprocité des contraintes tangentielles" (Fig. 4.5). Sur
deux facettes orthogonales, les composantes tangentielles 4 perpendiculaires
à l’arête commune sont égales et simultanément dirigées vers cette arête ou
opposées à celle-ci (Fig. 4.5).

Figure 4.5 – Réciprocité des contraintes tangentielles

4. Dans l’objectif de différencier physiquement la nature de la contrainte normale et

la contrainte tangentielle, cette dernière peut être notée fi au lieu de ff, c’est la notation
de l’ingénieur. e.g. le glissement des terrains est effectué à un angle de 45‹ , i.e. c’est la
contrainte tangentielle qui prédomine, l’ingénieur utilise la notation fi pour enlever toute
confusion.

1ère année du cycle d’Ingénieurs 35 Mécanique des Solides

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

4.2.4 Équilibre statique

Dans le cas d’équilibre statique, les termes fonction du temps dispa-

raissent (accélération ai = 0). Les équations (4.17b) deviennent
@ffji
+ bi = 0 (4.21)
@xj
Les équations (4.21) s’appelles les équations d’équilibre en volume.

4.3 Tenseur déviateur et tenseur sphérique

On appelle contrainte normale moyenne la quantité définie par l’équation

(4.22).
1 1 1
ff0 = ffii = (ff1 + ff2 + ff3 ) = I1 (4.22)
3 3 3
C’est un invariant, il permet d’écrire le tenseur contrainte ff comme la somme
de deux tenseurs (Eq. 4.23 et 4.24).

ffij = ff0 ‹ij + sij (4.23)

ou encore, matriciellement
0 1 0 1
B ff0 0 0 C B ff11 ` ff0 ff12 ff13 C
B C B C
ff = B
B 0 ff0 0 C
C + B
B ff21 ff11 ` ff0 ff23 C
C (4.24)
@ A @ A
0 0 ff0 ff31 ff32 ff11 ` ff0

Dans le premier tenseur, dit tenseur contrainte sphérique, toutes les contraintes
normales sont égales et les contraintes tangentielles sont nulles. Ce ten-
seur caractérise par exemple un état de contrainte hydrostatique. Le second
tenseur, dit tenseur déviateur des contraintes, s’écrit

sij = ffij ` ff0 ‹ij (4.25)

Le tenseur déviateur sij a, pour direction principales, les mêmes directions

que celles du tenseur contrainte ffij , et pour contraintes principales

s1 = ff1 ` ff0 (4.26a)

s2 = ff2 ` ff0 (4.26b)
s3 = ff3 ` ff0 (4.26c)

Les invariants du déviateur sont définis dans les équations (4.27).

Is1 = 0 (4.27a)
1
Is2 = ` sij sij (4.27b)
2
1
Is3 = sij sjk ski = det s (4.27c)
3

Mécanique des Solides 36 1ère année du cycle d’Ingénieurs

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

Vu la simplification apportée par le premier invariant dans l’équation ca-

ractéristique, on peut trouver les contraintes principales d’un tenseur ffij en
cherchant d’abord celles de son déviateur sij , puis en utilisant les équations
(4.27).

4.4 Espace des contraintes principales

Le tenseur des contraintes (comme tout tenseur symétrique) peut être

caractérisé par les trois contraintes principales et l’orientation du repère
principal. Dans de nombreux cas, l’orientation du repère principal ne joue
pas un rôle essentiel, et on pourra caractériser le tenseur des contraintes
par les trois contraintes principales ff1 , ff2 et ff3 . On peut donc représen-
ter un tenseur des contraintes par un point d’un espace à trois dimension
(O; ff1 ; ff2 ; ff3 ). Au tenseur ff on associe le point ˚ (état de contrainte)
ayant comme coordonnées 5 les contraintes principales ff1 , ff2 et ff3 de ff.
Cette représentation, très utile, exige néanmoins certaines précautions : on

Figure 4.6 – Espace des contraintes principales

représente géométriquement l’espace des contraintes principales par une es-

pace vectoriel, alors que ce n’est pas le cas 6 .
Dans cet espace, les tenseurs sphériques sont représentés par les points de
l’axe hydrostatique (´) de cosinus directeurs c1 = c2 = c3 = p1 . Les ten-
3
seurs déviateurs sont représentés par les points du plan déviatorique (ı) (Eq.
4.28), perpendiculaire en O à l’axe hydrostatique (´).

ff1 + ff2 + ff3 = 0 (4.28)

La décomposition d’un tenseur en partie sphérique et déviateur correspond à

la projection orthogonale sur ´ et (ı). En particulier, la projection sur (´)
5. Le repère (O; ff1 ; ff2 ; ff3 ) est considéré orthonormé.
1 2
6. La somme de deux tenseurs ffij et ffij n’est pas une somme vectorielle.

1ère année du cycle d’Ingénieurs 37 Mécanique des Solides

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

Figure 4.7 – Projection des vecteurs ffi sur le plan déviatorique

est caractérisée par la trace de ff. Dans le plan déviatorique (ı) on trace la
projection des axes (Off1 , Off2 et Off3 ), qui font entre eux un angle de 2ı
3
et un tenseur ff sera représenté par le point S

S = ff1 h1 + ff2 h2 + ff3 h3 (4.29)

h1 , h2 et h3 étant les trois vecteurs unitaires portés par les axes Off1 , Off2
et Off3 -plutôt, par leurs projections- néanmoins, on les note Off1 , Off2 et
Off3 . En particulier, on vérifie que le point S ainsi construit caractérise le
déviateur, puisque, si l’on rajoute à ff un tenseur sphérique arbitraire, le point
S ne change pas, car h1 + h2 + h3 = 0.
On peut montrer que la position de point S est complètement caractérisée
par les deux invariants Is2 et Is3 (Eq. 4.27b et 4.27c).
Finalement, on constate que la position du point ˚ dans l’espace des contraintes
principales est complètement caractérisé par I1 , Is2 et Is3 :

— I1 fixe la projection sur la trisectrice (´) ;

— Is2 fixe la distance à (´) ;
— Is3 Fixe l’orientation de la projection de O˚ ~ sur le plan (ı)

Exercice 4.1 La matrice représentant le tenseur contrainte au point O d’un

solide vaut :
0 1
B 2 `2 0 C
B C
ff = B
B `2 3 1 C
C (4.30)
@ A
0 1 2
Trouver les composantes du vecteur contrainte (ou traction de surface T (x; n))
agissant en O sur un plan parallèle au plan ABC.

Mécanique des Solides 38 1ère année du cycle d’Ingénieurs

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

Figure 4.8 – Exercice 4.1

Exercice 4.2 Les contraintes dans un solide valent :

0 1
B 0 `2x1 3x23 C
m`2 ]
B C
ff(x) = B
B `2x1 0 0 C [N
C (4.31)
2
@ A
3x3 0 0

Trouver :
1. Les composantes de la force de volume au point A(2; 4; 3) [m] dont
l’accélération vaut a(3; 3; 3) [m s`2 ] sachant que ȷ = 2kg m`3 ;
2. Le vecteur contrainte en A sur le plan d’équation 12x1 +3x2 +4x3 `48 =
0.

1ère année du cycle d’Ingénieurs 39 Mécanique des Solides

 CHAPITRE 4. TENSEUR DES CONTRAINTES ET ÉQUATIONS D’ÉQUILIBRE

Mécanique des Solides 40 1ère année du cycle d’Ingénieurs

Chapitre 5

Tenseur des déformations et

équations de compatibilité

5.1 Introduction

Sous l’action d’un chargement, tous les corps subissent un changement

de forme, de position et d’orientation entre un état initial avant application
de la charge et un état final consécutif à son application. Les déforma-
tions caractérisent le changement de forme local en tout point du matériau,
indépendamment de sa nature et de ses caractéristiques de résistance.

5.2 État local de déformation

La transformation géométrique qui décrit le passage d’un petit élément

de matière de l’état initial à l’état final se décompose en Translation, Rota-
tion et Déformation, cette dernière étant seule responsable du changement
de forme du petit élément. Translation et Rotation sont des mouvements de
corps rigides qui traduisent les changements de position et d’orientation. Du
point de vue de l’état final, Rotation et Déformation ne sont pas commu-
tatives. Si l’on ne s’intéresse qu’au changement de forme du petit élément
et non plus à sa position absolue les deux configurations finales sont équi-
valentes pour un chargement donné. Seule la Déformation responsable des
variations de longueur et des distorsions angulaires caractérise le changement
de forme. Elle seule sera liée à la contrainte appliquée au travers de la loi
de comportement du matériau. C’est la variable pertinente du mécanicien.

5.2.1 Petites et grandes déformations

Cependant, lors de l’application de deux chargements successifs, il faut

tenir compte de l’ordre des chargements puisque le second agit sur un état

41
CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

Figure 5.1 – Transformation géométrique de passage entre l’état initial et

l’état final

intermédiaire – Position, orientation et forme – induit par le premier, état

qui devient le nouvel état de référence pour le calcul de la répartition des
contraintes induit par le second chargement. Le calcul de la répartition des
contraintes induit par le second chargement doit s’effectuer à partir de la
forme modifiée par le premier. L’état de contrainte est donc différent de
celui de la Déformation induite par le premier.
Ce n’est que dans le cas où les variations de position, orientation et forme
seront suffisamment faibles que l’on pourra les négliger et considérer qu’en
première approximation le second chargement agit également sur la géo-
métrie initiale, négligeant les variations de géométrie induites par le premier
chargement. Il faudra donc distinguer le cas des petites déformations comme
celles des matériaux cristallins travaillant en régime élastique de celui des
grandes déformations comme celles qui interviennent dans les opérations de
formage et d’emboutissage des métaux en régime plastique ou celles des
élastomères supportant de grandes déformations en régime élastique.

5.2.2 Principe de superposition

Une approche classique de résolution d’un problème mécanique – connue

sous le nom de principe de superposition – consiste à décomposer la sollici-
tation appliquée à une structure en une somme de sollicitations élémentaires
pour lesquelles il est aisé de déterminer, à partir de la géométrie initiale,
l’état de contrainte et la déformation induite correspondante, puis de déter-
miner l’état de déformation global par addition de ces états élémentaires
de déformation. Cette approche se réfère toujours à la géométrie initiale et
n’est donc applicable qu’en petites déformations.

Mécanique des Solides 42 1ère année du cycle d’Ingénieurs

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

5.2.3 Extension et glissement

A l’échelle microscopique, le changement de forme se caractérise par

un changement de longueur et un changement d’orientation des liaisons
atomiques. A l’échelle macroscopique, lorsqu’une structure se déforme, un
segment de longueur l voit varier sa longueur, mais aussi sa direction.

Figure 5.2 – Extension et glissement

l + ´l
- La variation de longueur se caractérise par l’extension e = ou
l
´l
l’allongement relatif " = .
l
- La variation de direction est pour sa part le résultat de la composition
des deux opérations de Rotation et de Déformation. Afin de s’affranchir
de la rotation et de ne retenu que la distorsion de forme, il suffit de
suivre l’évolution au cours de la déformation de l’angle ’ formé par
deux directions initialement orthogonales. Le glissement ‚ est défini
ı
comme l’angle complémentaire ’ = .
2
L’allongement relatif " étant défini par le rapport de deux longueurs et le
glissement ‚ étant défini par un angle, la déformation est une grandeur sans
dimension, donc sans unité de mesure.

5.3 Les grandes déformations

Nous nous limiterons à la présentation simplifiée dans un essai classique

de traction. On enregistre l’allongement ´l d’une éprouvette de section ini-
tiale S
0 et de longueur initiale l0 en fonction de la force appliquée F .
Pour caractériser le comportement du matériau on s’affranchit de la géomé-

1ère année du cycle d’Ingénieurs 43 Mécanique des Solides

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

trie de l’éprouvette en introduisant la contrainte nominale de traction (Eq.

5.1).
F
ffn = (5.1)
S0

L’extension est donnée par l’équation (5.2).

l + ´l ´l
e= =1+ = 1 + "n (5.2)
l l

´l
"n = étant la déformation nominale d’extension. La courbe d’essai F =
l
F (´l) tracée dans le plan est de caractère brut, elle ne caractérise que
l’éprouvette étudiée. Par contre ffn = ffn ("n ) est une courbe représentative
du comportement en traction du matériau étudié. Tant que l’allongement ´l

Figure 5.3 – Déformation nominale et vraie

reste petit devant la longueur initiale l0 , typiquement tant "n < 1%, on peut
négliger les variations de forme de l’éprouvette et les grandeurs nominales
ffn et "n sont représentatives de l’état de contrainte et de déformation vraies
du matériau. C’est le régime des petites déformations.
Cependant, lorsque l’éprouvette s’allonge de manière importante, comme
c’est le cas pour les métaux en régime plastique pour lesquels la déformation
nominale à rupture atteint quelques dizaines de % ou les élastomères en
régime élastique pour lesquels elle peut atteindre quelques centaines de
%, on constate expérimentalement que l’aire de section droite S décroît
au fur et à mesure de l’allongement de l’éprouvette. Il est donc naturel,
pour s’affranchir de la géométrie de l’éprouvette, de prendre à tout instant
du chargement les dimensions vraies S < S0 et l > l0 pour définir les
vraies contrainte et déformation. La définition de la déformation vraie devient
incrémentale au sens où, sous l’action de la contrainte vraie ff = FS à un
instant donné du chargement, le matériau subit un accroissement d’extension
de = l+dl
l
auquel correspond l’accroissement de déformation vraie d" = dll . La
déformation vraie " lorsque l’éprouvette s’est allongée de ´l s’obtient par
l’intégrale définie dans l’équation (5.3).

Mécanique des Solides 44 1ère année du cycle d’Ingénieurs

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

5.4 Les petites déformations

La plupart des structures étant constituées de matériaux cristallins peu

déformables travaillant en régime élastique, l’approximation des petites dé-
formations est largement suffisante. Elle offre de plus l’avantage du principe
de superposition. Nous nous limiterons à la présentation de la théorie des
petites déformations.
 l0 +´l
dl ´l
"= = ln(1 + ) = ln(1 + "n ) (5.3)
l0 l l

5.4.1 Le vecteur déplacement

Plutôt que caractériser la transformation affine tangente par les positions

finales x des points matériels de position a dans l’état initial, il est plus
commode d’étudier le champ de déplacement u induit par la transformation
et défini par l’équation (5.4).

Figure 5.4 – Vecteur déplacement et affinité tangente

u(a) = x ` a (5.4)

En appelant G = grad(u) le tenseur gradient du champ de déplacement au

point P0 , la transformation affine tangente dx = F (a) da s’écrit alors, en
termes de déplacement, conformément à l’équation (5.5a et 5.5b).

du = G(a) da (5.5a)
@ui
dui = daj (5.5b)
@Aj
G =F `I (5.5c)

Notons que G, comme F , contient lui aussi les informations locales de

changement de forme (déformation) et d’orientation (rotation) de la boule
initiale.

1ère année du cycle d’Ingénieurs 45 Mécanique des Solides

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

La déformation étant seule responsable des changements de longueur, le

segment da de longueur initiale L (Eq. 5.6) prend par déformation pure la
longueur finale (Eq. 5.6).
L2 = (da)T da (5.6)

l2 = (dx)T dx avec dx = da + du (5.7)

On trouve
l2 = L2 + (du)T da + (da)T du + (du)T du (5.8)

comme
du = G da

On trouve
l2 = L2 + (da)T (G T + G + G T G)da (5.9)
2 2
Le rapport l `L
L2
associé à la variation de longueur du segment da est fonction
du tenseur symétrique
GT + G + GT G

qui caractérise la déformation pure.

5.4.2 Gradient de déplacement et petites déformations

La description mathématique des déformations se simplifie considérable-

ment lorsque la forme initiale de la structure est peu modifiée sous l’action
du chargement. Dans ce cas, la variation de forme de la boule initiale de
centre P0 est telle que dx = da + du ı da soit F = I + G ı I. Les varia-
2 2
tions relatives de longueur sont très petites l +L
L2
ı l`L
L
et les composantes
Gij de G sont petites devant l’unité. On peut donc négliger le terme G T G
par rapport au terme G T + G et la déformation pure est représentée par la
partie symétrique du gradient de déplacement G. Sa partie complémentaire
antisymétrique représentera la rotation locale.
L’approximation des petites déformations s’applique à la plupart des struc-
tures travaillant en domaine élastique pour lesquelles Gij dépasse rarement
10`3 .

5.4.3 Translation, rotation et déformation

On décompose G en deux parties, symétrique et antisymétrique (Eq.

5.10). On verra par la suite que la partie symétrique présente da déformation
et la partie antisymétrique présenta la rotation(Eq. 5.10).
" ˙
z }| { z }| {
1 1
G = (G + G T ) + (G ` G T ) (5.10)
2 2

Mécanique des Solides 46 1ère année du cycle d’Ingénieurs

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

Autrement-dit, les tenseurs " et ˙ sont liés à G par les équations (5.11 et
5.12). 0 1
1 T 1 @ @ui @uj A
" = (G + G ) ou "ij = + (5.11)
2 2 @aj @ai
0 1
1 1 @ui @uj A
˙ = (G ` G T ) ou ˙ij = @ ` (5.12)
2 2 @aj @ai
C’est clair que le tenseur " est symétrique alors que ˙ est un tenseur anti-
symétrique (Eq. 5.13 et 5.14).

" = "T ou "ij = "ji (5.13)

˙ = `˙T ou ˙ij = `˙ji (5.14)

Partie antisymétrique du tenseur de déplacement

La partie antisymétrique ˙ qui ne porte que trois composantes indé-

pendantes induit dans la boule de centre P0 un champ de déplacement
du = ˙ da, qui correspond à la rotation locale de la boule de vecteur !
(Eq. 5.15).
du = ˙ da = ! ˆ da (5.15)

Le vecteur rotation locale ! définissant l’axe, le sens et l’amplitude de la

rotation locale a pour composantes !1 , !2 et !3 définies dans les équations
(5.16a, 5.16b et 5.16c).

!1 = `˙23 (5.16a)
!2 = `˙31 (5.16b)
!3 = `˙12 (5.16c)

Soit !i = `˙jk , avec i; j; k la permutation circulaire du groupe 1; 2; 3. C’est

l’analogue du vecteur tourbillon de l’hydrodynamique.

Figure 5.5 – Notion de déplacement : translation Q0 Q0 , rotation Q0 Q00 et

déformation Q00 Q

1ère année du cycle d’Ingénieurs 47 Mécanique des Solides

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

Partie symétrique du tenseur gradient de déplacement

La partie symétrique ", qui ne comporte que six composantes indépen-

dantes, induit dans la boule, de centre P0 , un champ de déplacement du =
" da correspondant à la déformation locale. cette partie du champ de dépla-
cement, associée à la déformation pure, est la seule responsable du change-
ment de forme de la boule.
Le déplacement u d’un point Q0 de la boule (Fig. 5.5 et Eq. 5.17) est la
somme :

- D’un déplacement de corps rigide composé de la translation u0 du

centre P0 de la boule et de la rotation locale ! ˆ a de la boule
autour du vecteur !(P0 ) associé à la partie antisymétrique ˙.
- D’un déplacement " da responsable du changement de forme par la
déformation associée à la partie symétrique "(P0 ).
T ranslation Rotation Dformation
z }| { z }| { z }| {
u(Q0 ) = u(P0 ) + !(P0 ) ˆ da + "(P0 )da (5.17)

5.5 Tenseur des déformations

5.5.1 Définition

Dans le cadre de la théorie des petites déformations, l’état de déformation

en un point P est donc caractérisé par le tenseur symétrique "(P ). La relation
donnant le déplacement d’extrémité u(P; l) d’un segment unitaire d’origine
P et de direction l est définie par l’équation (5.18).

u(P; l) = "(P ) l (5.18)

Comme dans le cas du tenseur des contraintes, la connaissance des vecteurs

déplacement d’extrémité u1 , u2 et u3 de trois segments orthogonaux de
longueur unitaire d’origine P respectivement orientés selon les axes Px1 , Px2
et Px3 suffit pour déterminer en ce point P l’état de déformation. La relation
(5.18) entre u(P; l) et l s’écrit de manière générale à l’aide du tenseur
"(P ) dont les composantes dans le repère (P; x1 ; x2 ; x3 ) sont définies par la
formule (5.19). 0 1
B "11 "12 "13 C
B C
"(P ) = B
B "21 "22 "23 C
C (5.19)
@ A
"31 "32 "33

5.5.2 Repère local, extension, distorsion et glissement

Au point P où l’état local de déformation est caractérisé par le tenseur

"(P ), le déplacement d’extrémité u(P; l) d’un segment unitaire d’origine P

Mécanique des Solides 48 1ère année du cycle d’Ingénieurs

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

orienté dans la direction l induit, d’une part, une variation de longueur du

segment et d’autre part un changement d’orientation. Pour caractériser ces
deux effets, introduisons le repère local trièdre direct (P; l; t) associé à la
direction l telles que les vecteurs l, r et u soient coplanaires (Fig. 5.6).
Dans ce repère, les composantes du vecteur déplacement u(P; l) sont "ll , "rl
et "tl . Le premier indice, faisant référence à l’axe de projection et le second
à l’orientation l du segment. Par construction, on a toujours "tl = 0 et les
composantes non nulles de u sont sa composante normale "ll = u(P; l) ´ l
selon l appelée extension et sa composante perpendiculaire "rl = u(P; l) ´ r
située dans le plan de glissement l; r appelée distorsion.

5.5.3 Propriétés des déformations

Le tenseur "(P ) définissant l’état de déformation en un point P est un

opérateur intrinsèque, symétrique, par construction, qui lie le vecteur dépla-
cement d’extrémité u(P; l) du segment unitaire l d’origine P à son orienta-
tion. Cette définition est indépendante du repère choisi, seule sa représenta-
tion, les valeurs de ses composantes, dépend du repère.
Comme le tenseur des contraintes, il possède les propriétés communes aux
tenseurs symétriques d’ordre 2. Toutes les propriétés du tenseur des contraintes
ff -en particulier : valeurs propres, directions principales, invariants, sphérique
et déviateur- s’appliquent au tenseur des déformations ". Il suffit de rem-
placer la normale n à la facette par le segment l et le vecteur contrainte
T (P; n) par le vecteur déplacement d’extrémité u(P; l).

5.6 Interprétation physique

Pour saisir le signification des quantités "ij , il suffit de raisonner dans un

plan (x; y) par exemple (Fig. 5.6). Soit AB et AC deux éléments de ligne de
la configuration initiale, perpendiculaires entre eux, qui se déplace en A0 B 0
et A0 C 0 dans la configuration déformée. Les déplacements sont petits, bien
que pour la clarté du dessin, il n’apparaissent pas tels. La figure montre les
deux composantes du déplacement des trois points A, B et C. On note ‚1
et ‚2 les petits angles dont tournent les éléments dx et dy, de sorte que
l’angle droit en A varie de la quantité totale ‚ = ‚1 `‚2 dans le changement
de configuration (‚ = ı2 ` ¸ est positif si l’angle droit diminue).
Calculons la dilatation des éléments de ligne,

@u
dx
A0 B 0 ` AB @u
= @x =
AB dx @x

1ère année du cycle d’Ingénieurs 49 Mécanique des Solides

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

Figure 5.6 – Interprétation géométrique des composantes de la déformation.

Les déplacements de leurs accroissements étant positifs, ‚2 est négatif

@v
dy
A0 C 0` AC @y @v
= =
AC dy @y
avec : cos ‚1 ı 1 et cos ‚2 ı 1.
On constate alors que les déformations normales représentent les dilatations
habituelles. Calculons maintenant la variation de l’angle droit,
@v
v+ dx `v @v
‚1 ı tan ‚1 = @x
ı avec dx ı dx0 (5.20a)
dx0 @x
@u
u+ @y
dy `v @u
‚2 ı tan ‚2 = ı avec dy ı dy 0 (5.20b)
dy 0 @y
@u @v
‚ = ‚1 ` ‚2 = + (5.20c)
@y @x

5.7 Glissement

La quantité ‚ (Eq. 5.20c), très employée par les ingénieurs, s’appelle

le glissement. La raison en est la suivante : Considérons un carré (ou rec-
tangle) élémentaire dans le plan (x; y) (Fig. 5.7) ; si l’on envisage que la
déformation tangentielle, à l’exclusion de toute dilatation des côtés, ce carré
(ou rectangle) se transforme en un losange (ou parallélogramme), comme
indiqué d’une manière exagérée dans la figure.
Déplaçons maintenant, en bloc, la configuration déformée A0 B 0 C 0 D 0 de ma-
nière à amener A0 B 0 (ou A0 C 0 ) en coïncidence avec AB ou AC (Fig. 5.7b ou
5.7c). On voit que l’angle CAC 0 (ou BAB 0 ) est égal à ‚ = ‚1 + j‚2 j ; par

Mécanique des Solides 50 1ère année du cycle d’Ingénieurs

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

(a) (b) (c)

Figure 5.7 – Glissement ‚

conséquent, on peut caractériser la déformation tangentielle par le glisse-

ment du côté CD (ou BD) parallèlement au côté AB (ou AC), glissement
défini par la quantité sans dimension
CC 0 BB 0
‚ = tan ‚ = =
AC AB
qui donne, sous forme d’un petit angle [rad], le déplacement d’un coté par
rapport à l’autre par unité de distance entre ces deux côtés (si AC = 1,
‚ = CC 0 ou si AB = 1, ‚ = BB 0 ).
En conclusion, quel que soit le plan (xi ; xj ) considéré (i 6= j), le glissement
représente le double de la déformation tangentielle.

‚xy = 2"xy (5.21a)

‚yz = 2"yz (5.21b)
‚xz = 2"xz (5.21c)

Au tenseur déformation peuvent donc s’associer les deux formes matricielles

0 1 0 1
1 1
B "11 "12 "13 C B "11 ‚
2 12
‚
2 13 C
B C B 1 1 C
"= B
B "21 "22 "23 C
C = B 2 ‚21
B "22 ‚
2 23
C
C (5.22)
1 1
@ A @ A
"31 "32 "33 ‚
2 31
‚
2 32
"33

Exercice 5.1 Un état de déformation vaut, dans le système d’axes (xi ),

0 1
B 2 0 `4 C
B C
" = kB
B 0 3 0 C
C (5.23)
@ A
`4 0 `4

où k est une constante (très petite).

1. Calculer le tenseur déformation dans les axes (x0 ) définis par la matrice
des cosinus directeurs.
0 1
1 0 1
1 B
B p C
C
[C] = p B
B 0 2 0 C
C (5.24)
2 @ A
`1 0 1

1ère année du cycle d’Ingénieurs 51 Mécanique des Solides

CHAPITRE 5. TENSEUR DES DÉFORMATIONS ET ÉQUATIONS DE COMPATIBILITÉ

2. Obtenir les invariants de l’état de déformation dans chacun des systèmes

d’axes.
3. Trouver les dilatations principales.

Exercice 5.2 On donne, dans le plan (x1 ; x2 ), trois champs de déplacement

où k est un petit paramètre pour chacun d’eux :
(a) u1 = `kx1 et u2 = kx1 (b) u1 = `kx2 et u2 = 0 (c) u1 = kx2 et
u2 = kx1
1. Dessiner la configuration déformée du carré unité situé dans le premier
quadrant ; visualiser le paramètre k.
2. Calculer la rotation matérielle.

Mécanique des Solides 52 1ère année du cycle d’Ingénieurs

Chapitre 6

Loi de comportement

6.1 Essais de la traction simple

6.1.1 Cas de l’acier doux

Dans l’objectif de mettre en évidence la loi de comportement, on préfère

la présenter d’abord dans le cas de chargement uniaxial, puis, on généralise.
En pratique, l’acier doux occupe une place importante dans le domaine de
la construction métallique, la construction mécanique, la fabrication des
pipelines, les tripodes de transport d’énergie électrique, les réservoirs, les
appareils à pression, les silos, les racks, etc... Vu son omniprésence dans
beaucoup de domaines, on prend ce matériau métallique comme un exemple
pour étudier la courbe ff = ff("). La figure (6.1) présente la variation de

Figure 6.1 – Courbe contrainte-déformation pour un acier doux

la contrainte nominale en fonction de la déformation (dilatation dans ce

cas), suite à un chargement axial monotone appliqué sur une éprouvette
normalisée d’un acier doux. Au début du chargement, la contrainte augmente
linéairement en fonction de la déformation en traçant la droite OA appelée
la droite de Hooke (Eq. 6.1).

ff = E" (6.1)

53
 CHAPITRE 6. LOI DE COMPORTEMENT

E est le module d’Young, appelé parfois, le module de rigidité longitudi-

nale.
Le matériau est dit élastique lorsqu’il revient à son état initial le temps
quand les efforts extérieurs cessent. Ainsi, la limite de proportionnalité ne
coïncide pas toujours avec la limite élastique. Dans ce cas, la contrainte
au point A est appelé limite de proportionnalité notée ffp . Par contre, la
contrainte au point B est appelée limite élastique ou limite d’étirage notée
ffe .
Entre les points B et D, on constate une dilatation importante de l’éprou-
vette de traction avec une contrainte nominale constante. C’est le palier
d’étirage qui caractérise l’acier doux (ou extra-doux). On dit que le maté-
riau est plastifié. En réalité, dans ce palier, l’acier ne présente plus un état
homogène. Il est composé des particules élastiques et d’autres qui sont to-
talement plastifiées.
A partir du point D, le métal rentre dans un nouvel état appelé l’écrouissage.
i.e. la contrainte nominale augmente avec l’augmentation de la déformation
totale dans la phase plastique. Au niveau du point E, la machine de traction
enregistre une valeur extrémale de la contrainte. Cette valeur est appelée
la contrainte de rupture notée ffR . Cette appellation ne reflète pas la réa-
lité puisque la rupture réelle s’effectue au point ultérieur F . La déformation
enregistrée à ce point est appelée l’allongement à la rupture.
Dans la zone de striction, la section de l’éprouvette diminue d’une façon
remarquable, i.e. la contrainte vraie est supérieure à la contrainte nominale
qui continue à se décroitre jusqu’à la rupture. Néanmoins, l’hypothèse des
petites déformations est toujours adoptée. En état de service de la structure
ou de la machine, l’élément structural ou la pièce en acier ne travaille pas
dans la zone de striction. Ainsi, on n’étudiera pas ici cette partie de la courbe
ff = ff(").
On suppose maintenant que l’effort de la traction cesse au niveau du point
G du diagramme ff = ff("). Dans ce cas, l’allure du déchargement 1 suit
le chemin d’une ligne droite GH parallèle à la droite de Hooke. Lorsque le
matériau est chargé à nouveau, on remarque que le métal conserve le même
module d’Young initial mais avec une nouvelle limite élastique qui coïncide
cette fois avec la nouvelle limite de proportionnalité. Cette nouvelle limite
élastique s’appelle la limite élastique actuelle ou le seuil plastique actuel.
Ceci prouve que le métal a une mémoire qui lui rappelle, à chaque nouvel
état de chargement, sa dernière limite élastique.

1. On note que le déchargement est toujours élastique bien que le matériau soit plastifié

Mécanique des Solides 54 1ère année du cycle d’Ingénieurs

 CHAPITRE 6. LOI DE COMPORTEMENT

6.1.2 Contraction latérale

Vu l’homogénéité du matériau objet de l’essai de traction, on va choisir

une des dimensions transversales de l’éprouvette objet de l’essai. Si cette
dernière est de forme cylindrique, le diamètre c’est la dimension de référence.
Pour étudier la contraction latérale, on note respectivement D et D ˜ les
diamètres de l’éprouvette avant et au moment du chargement.
On appelle "lat la déformation latérale la quantité définie dans l’équation
(6.2).
D˜ ` D
"lat = (6.2)
D
Cette grandeur est négative puisqu’avec l’effort de traction on observe un
rétrécissement de la section de l’éprouvette. La figure (6.2) illustre la va-
"
riation du rapport ” = lat en fonction de la déformation longitudinale ".
"
On remarque que la zone OA correspond à la phase élastique linéaire, i.e.
dans cette phase, le rapport ‌ est constant. Cette constante est appelée
coefficient de Poisson noté ‌. C’est l’une des grandeurs mécaniques qui ca-
ractérisent le matériau. A partir du point A, le rapport ‌ tend vers la valeur
0; 5. Pour les métaux, le coefficient de Poisson est approximativement égal
à 0; 33.

Figure 6.2 – Illustration du coefficient de Poisson

6.1.3 Cas de l’alliage d’aluminium

La figure (6.3) présente la courbe de traction ff = ff(") d’un alliage d’alu-

minium. Le parcours OA présente la phase élastique linéaire dans laquelle
la contrainte nominale est proportionnelle à la déformation. Le chemin BE
présente le phénomène d’écrouissage. A partir du point E la courbe ff = ff(")
décroit jusqu’à atteindre la rupture au point F . Comme dans le cas de l’acier

1ère année du cycle d’Ingénieurs 55 Mécanique des Solides

 CHAPITRE 6. LOI DE COMPORTEMENT

doux, la phase EF s’appelle la zone de striction.

Contrairement à la courbe de traction relative relative à l’acier doux, celle de
l’alliage d’aluminium n’admet pas un palier d’étirage. De ce fait, on définit
une limite élastique à 0; 2% notée ff0;2 . Cette grandeur généralise la notion
de la limite élastique ffe . Concernant la déformation latérale, l’aluminium
se comporte de la même façon que l’acier doux. Si à partir d’un point G

Figure 6.3 – Courbe contrainte-déformation pour un acier doux

appartenant à la zone d’écrouissage on effectue un déchargement dont le

chemin suit une ligne droite parallèle au chargement élastique initial. Lors-
qu’on recharge une deuxième fois à partir du point H, on atteint le point G
avant de poursuivre le chemin de chargement plastique, i.e. le déchargement
n’a pas d‘influence sur le chargement monotone.

6.2 Hypothèses de linéarité

Dans ce paragraphe, on va citer trois hypothèses de linéarité qui peuvent

être supposées dans un problème mécanique. Si les deux premières hypo-
thèses sont déclarées en même temps, on peut juger que le problème est
linéaire. Par contre, si l’une des deux n’a pas lieu, alors le problème est dit
non linéaire.

6.2.1 Linéarité géométrique

Prenons l’exemple d’une poutre en porte à fou avec un effort vertical

appliqué au niveau de son extrémité (Fig. 6.4). Lorsqu’on suppose que l’effort
appliqué reste perpendiculaire à la poutre avant et après la déformation
de la poutre, ceci implique l’hypothèse de la linéarité géométrique. C’est
l’exemple de la pression qui reste toujours perpendiculaire à la paroi interne

Mécanique des Solides 56 1ère année du cycle d’Ingénieurs

 CHAPITRE 6. LOI DE COMPORTEMENT

d’un réservoir ou d’un appareil à pression, on n’a pas besoin de réactualiser

les calculs grâce à la géométrie supposée invariable. Pour les poutres, on
suppose que la flexion ne cause pas la dilatation (on néglige l’extension de la
poutre). Autrement-dit, la linéarité géométrique c’est l’hypothèse classique
des petites déformations.

Figure 6.4 – Poutre en porte à fou soumise à un chargement de singularité

P

6.2.2 Linéarité matérielle

On dit qu’un problème mécanique est linéaire de point de vue maté-

riau lorsque le module d’Young et le coefficient de Poisson caractérisant
ce dernier restent constants pendant tout le processus de chargement de la
structure ou de la pièce étudiée. i.e. tout état d’élasticité non linéaire ou de
plastification du matériau est exclu.

6.2.3 Linéarité des conditions aux limites

Soit un solide ou un système mécanique soumis à un chargement extérieur.

Si l’inversion des conditions aux limites implique l’inversion des résultats,
alors, le problème mécanique est linéaire de point de vue conditions aux
limites. e.g. la poutre en porte à fou de la figure (6.4) est soumise à un
chargement dirigé vers le bas, par conséquent, le déplacement transversal
est dirigé aussi vers le bas. Par contre, si on remplace l’effort P est par –P ,
le déplacement transversal s’inverse et la poutre fléchie vers le haut.
On considère maintenant une poutre dont l’une des extrémités et encastrée
et l’autre est posée sur un appui simple, soumise à un effort uniformément
répartit dirigé vers le bas. Ceci cause une flexion de la poutre vers le bas.
Naturellement, lorsqu’on inverse l’effort, i.e. on remplace l’effort répartit
p par –p, le déplacement transversal w 0 est dirigé vers le haut mais il est
différent de –w (Fig. 6.5). Ceci implique que ce problème n’est pas linéaire
de point de vu conditions aux limites.

1ère année du cycle d’Ingénieurs 57 Mécanique des Solides

 CHAPITRE 6. LOI DE COMPORTEMENT

Figure 6.5 – Poutre simplement appuyée

6.3 Aspects phénoménologiques

Tous les solides, à une température pas trop proche de la température

de fusion et pour une vitesse de sollicitation pas trop élevée, exhibent un
comportement ff = ff(") réversible, linéaire dans un premier temps, éventuel-
lement non linéaire ensuite. Il s’agit de mouvements réversibles d’atomes,
de molécules ou de cellules, qui correspondent à des déformations n’excé-
dant pas 0; 2 à 0; 5 10`2 pour les matériaux métalliques, les composites, les
bétons, les verres, les céramiques, les bois, ... Ce qui justifie l’hypothèse des
petites déformations pour ces matériaux. La dissipation est nulle. Par contre,
pour les matériaux caoutchoutiques les déformations de nature entropique
peuvent atteindre 500% et plus avec de fortes non-linéarités.
On ne considère que les matériaux homogènes, i.e. les caractéristiques du
matériau ne dépendent pas du point étudié. On se limite ici au domaine élas-
tique. Si le domaine subit un cycle de charge, le travail des forces extérieures
est nul, i.e. le matériau restitue intégralement l’énergie qui lui a été fournie
pour qu’il se déforme. L’élasticité est synonyme de réversibilité. On rajoute
une hypothèse supplémentaire, c’est celle de l’élasticité linéaire, c’est une
autre manière de dire que les contraintes varient proportionnellement à la
déformation. Les déformations doivent être suffisamment petites pour rester
dans la partie linéaire de la loi constitutive et présenter le même compor-
tement en traction et en compression. D’autre part, on suppose que la loi
reliant les déplacements aux déformations est également linéaire, i.e. on se
place donc dans les hypothèses de la linéarité géométrique et matérielle.
Le tenseur des contraintes et celui des déformations sont déterminés par la
connaissance de neufs éléments qui les caractérisent.

6.4 Élasticité linéaire

6.4.1 Domaine de validité et d’emploi

La limite du domaine d’élasticité dépend de la température et des sol-

licitations mécaniques antérieures, dues aux procédés d’élaboration ou au
service des pièces mécaniques. La température, favorisant les mouvements
irréversibles, diminue la taille de ce domaine jusqu’à la rendre nulle à la

Mécanique des Solides 58 1ère année du cycle d’Ingénieurs

 CHAPITRE 6. LOI DE COMPORTEMENT

fusion.
L’hypothèse d’isotropie élastique est très bien vérifiée pour les matériaux
métalliques et les céramiques, elle l’est un peu moins pour le béton, alors
que pour les composites et les bois, il faut utiliser des théories à anisotropie
particulière.
L’élasticité linéaire s’applique aux déformations élastiques infinitésimales et
linéaires (relation linéaire entre contraintes et déformations).

6.4.2 Formulation

Potentiel d’état

Pour obtenir une théorie linéaire, il suffit de choisir comme potentiel

d’état convexe, une forme quadratique définie positive des composantes du
tenseur des déformations (Eq. 6.3 et 6.4).
1
ȷ = ":E :" (6.3)
2
ou encore
1
ȷ =
"ij Eijkl "kl (6.4)
2
ȷ étant la masse volumique, E le tenseur d’ordre 4 des rigidités élastiques
ou tenseur de Hooke, dont les composantes sont les coefficients d’élasticité
isotherme éventuellement fonctions de la température.
Le tenseur des contraintes ff dérive d’un potentiel pour donner la loi
définie dans les équations (6.5) et (6.6).

@
ff = ȷ =E :" (6.5)
@"

ou encore
@
ffij = ȷ = Eijkl "kl (6.6)
@"ij
Le tenseur E vérifie certaines propriétés. En effet, la symétrie des tenseurs
des contraintes ffij = ffji et des déformation "ij = "ji conduit à une première
symétrie mineure (Eq. 6.7a et 6.7b).

Eijkl = Ejikl (6.7a)

Eijkl = Ejilk (6.7b)

L’existence du potentiel implique les symétries majeures (Eq. 6.8).

@ffij @2 @2
Eijkl = =ȷ =ȷ = Eklij (6.8)
@"kl @"ij @"kl @"kl @"ij

ff = E : " (6.9)

1ère année du cycle d’Ingénieurs 59 Mécanique des Solides

 CHAPITRE 6. LOI DE COMPORTEMENT

ou encore
ffij = Eijkl "kl (6.10)
Les équations (6.9) et (6.10) présentent la loi de Hooke généralisée du milieu
élastique linéaire, qui dépend de 21 coefficients, constantes indépendantes,
du tenseur symétrique élastique.
Par inversion, ou par définition du potentiel dual (Eq. ??), on introduit de
même le tenseur de souplesse, du quatrième ordre, symétrique, S = E `1
(Eq. 6.12).
1
ȷ ˜= ff : S : ff (6.11)
2
@ ˜
"=ȷ = S : ff (6.12)
@ff
Le cas de l’élasticité non linéaire se traite par le choix d’autres formes du
potentiel (") ou ˜ ("), qui doivent assurer leur convexité, par rapport aux
déformations ou par rapport aux contraintes selon la formulation.

Élasticité linéaire isotrope : loi de comportement

L’isotropie et linéarité imposent que le potentiel soit un invariant qua-

dratique du tenseur des déformations, i.e. une combinaison linéaire du carré
1 1
du premier invariant "2I = (tr"2 ) et du second invariant "II = " : " = (tr"2 )
2 2
(Eq. 6.13, 6.14 et 6.15).
1 1
ȷ = (–"2I + 4—"II ) = (–(tr") + 2—" : ") (6.13)
2 2
@
ff = ȷ = –(tr")I + 2—" (6.14)
@"
ffij = –"kk ‹ij + 2—"ij (6.15)
– et — sont les coefficients de Lamé, I est le tenseur unité de composantes
‹ij . Les composantes du tenseur de Hooke sont détaillées dans l’équation
(6.16).
Eijkl = –‹ij ‹kl + —(‹ik ‹jl + ‹il ‹jk ) (6.16)
En travaillant avec la formule duale " = S ff on obtient

1+‌ ‌
"= ff ` (trff)I (6.17)
E E
ou
1+‌ ‌
ffij ` ffkk ‹ij
"ij = (6.18)
E E
E est le module d’Young et ‌ le coefficient de Poisson.
Les composantes du tenseur de souplesse S = E `1 s’écrivent suivant l’équa-
tion (6.19).
`1 1+‌ ‌
Sijkl = Eijkl = (‹ik ‹jl + ‹il ‹ik ) ` ‹ij ‹kl (6.19)
E E

Mécanique des Solides 60 1ère année du cycle d’Ingénieurs

 CHAPITRE 6. LOI DE COMPORTEMENT

Les relations classiques entre les coefficients de Lamé, le module d’Young,

le coefficient de Poisson, le module de cisaillement G et le module de
compressibilité hydrostatique K.
‌E
–= (6.20)
(1 + ‌)(1 ` 2‌)

E
—=G = (6.21)
2(1 + ‌)
E 3– + 2—
K = = (6.22)
3(1 ` 2‌) 3
3– + 2—
E =— (6.23)
–+—
–
‌ = (6.24)
2(– + —)
Pour que la loi d’élasticité dérive d’un potentiel quadratique défini positif, la
seule condition que doivent vérifier les coefficients est définie dans l’équation
(6.25). Le module E étant supposé positif.
1
`1 < ‌ < (6.25)
2
Dans un repère quelconque, la relation de compliance et la loi de Hooke
généralisée pour un matériau isotrope (Eq. 6.26 et 6.27).
0 1 0 10 1
"11 1 `‌ `‌ 0 0 0 ff11
B C B CB C
B "22 C B `‌ 1 `‌ 0 0 0 CB ff22 C
B C B CB C
"33 C = B `‌ `‌ 1 0 0 0 ff33
B C B CB C
B C B CB C
B CB C (6.26)
B
B 2"23 C
C B 0
B
0 0 2(1 + ‌) 0 0 CB
CB ff23 C
C
B C B CB C
B
@ 2"13 C
A
B 0
@ 0 0 0 2(1 + ‌) 0 CB
A@ ff13 C
A
2"12 0 0 0 0 0 2(1 + ‌) ff12
0 1
0 1 1`‌ ‌ ‌ 0 0 0 0 1
ff11 B C "11
B C
B
B ‌ 1`‌ ‌ 0 0 0 C
CB C
B ff22 C B CB "22 C
B C B ‌ ‌ 1`‌ 0 0 0 CB C
ff33 "33
B C B CB C
B C
C= E B 1 ` 2‌ CB C
B B 0 0 0 0 0 CB C
B
B ff23 C
C (1 + ‌)(1 ` 2‌) B
B 2
CB
CB 2"23 C
C
B C B 1 ` 2‌ CB C
B ff13 C B 0 0 0 0 0 CB 2"13 C
2
@ A B C@ A
B C
ff12 @ 1 ` 2‌ A 2"12
0 0 0 0 0
2
(6.27)
Dans le cas d’un essai de traction uniaxiale sur une éprouvette de faible
diamètre vis-à-vis de sa longueur, dans le domaine élastique linéaire du ma-
tériau, E est le coefficient de proportionnalité entre la contrainte axiale
nominale ff et la déformation axiale mesurée ", c’est une caractéristique de
l’éprouvette. C’est la pente à l’origine de la courbe ff = ff(") pour un essai
de traction uniaxiale. Le coefficient qui relie les contraintes tangentielles
aux déformations angulaires est le module d’élasticité en cisaillement ou
le module de glissement noté G. Pour un matériau isotrope, on démontre

1ère année du cycle d’Ingénieurs 61 Mécanique des Solides

 CHAPITRE 6. LOI DE COMPORTEMENT

que le module de glissement est une fonction du module d’élasticité et du

coefficient de Poisson (Eq. 6.28).
E
G = (6.28)
2(1 + ‌)
Le module d’Young et le module de glissement sont homogènes à des
contraintes et s’expriment en NM `1 dans le système international. Quand
on exerce une traction sur une éprouvette, sa section diminue. C’est l’effet de
Poisson. Le coefficient de Poisson ‌ij mesure la déformation obtenue dans la
direction j quand on impose une déformation dans la direction i (Eq. 6.29).
"jj
‌ij = ` (6.29)
"ii

Mécanique des Solides 62 1ère année du cycle d’Ingénieurs

Chapitre 7

Forme intégrale de l’équilibre

et de la cinématique

7.1 Forme intégrale de l’équilibre

7.1.1 L’intégration par partie est la clef de la forme intégrale

Soit !i (x; y; z) un champ vectoriel continu et dérivable (de classe C 1 ).

Pour tout !j , on a 
@ffij
( + bj )!j dV = 0 (7.1)
V @xi
L’intégration par partie donne
    
@ffij @!j
!j dV + bj !j dV = ffij !j ni dA ` ffij dV + bj !j dV = 0
V @xi V V V @xi V

d’où   
@!j
dV = ffijtj !j dA + bj !j dV (7.2)
V @xi A V
Exploitons le fait que la forme différentielle de l’équilibre est symétrique,
on aura
0 1
@!j 1 @!j 1 @!i 1 @!j @!i A
ffij = ffij + ffji = ffij @ +
@xi 2 @xi 2 @xj 2 @xi @xj
On obtient finalement
 0 1  
1 @!i @!j A
ffij @ + dV = bi !i dV + ti !i dA (7.3)
V 2 @xj @xi V A

C’est la forme intégrale de l’équilibre.

Remarque

La forme intégrale des équations de mouvement s’écrit

  0 1  
1 @!i @!j A
ȷai !i dV + ffij @ + dV = bi !i dV + ti !i dA (7.4)
V V 2 @xj @xi V A

63
 CHAPITRE 7. FORME INTÉGRALE DE L’ÉQUILIBRE ET DE LA CINÉMATIQUE

La démonstration n’est pas difficile si l’on se rappelle du principe de d’Alem-

bert et on remplace bi par bi ` ȷai .

7.2 Forme intégrale de la cinématique

La forme différentielle de la cinématique est transformée en forme inté-

grale par l’expression
 2 0 13
1 @ui @uj A5
4"ij ` @ + !ij dV = 0 (7.5)
V 2 @xj @xi

où !ij (x; y; z) est un champ de tensoriel symétrique quelconque (!ij = !ji ),

continu et dérivable (de classe C 1 ). L’intégration par partie donne
  0 1 
@`
@!ji A
"ij !ij dV = ui dV + (nj !ji ) ui dA (7.6)
V V @xj A

C’est la forme intégrale de la cinématique.

7.3 Principe des déplacements virtuels (ou des travaux

virtuels)

Le champ vectoriel !i doit être continu et dérivable quelconque. On

va l’interpréter maintenant comme un champ de déplacement virtuel noté
‹ui . On peut remarquer que dans le terme de gauche de l’équation (7.3),
apparaissent les déformations virtuelles ‹"ij
0 1
1 @‹ui @‹uj A
‹"ij = @ + (7.7)
2 @xj @xi

L’équilibre intégral peut s’écrire

  
ffij ‹"ij dV = bi ‹ui dV + ti ‹ui dA (7.8)
V V A

En effet, le terme de droite de l’équation (7.8) s’appelle le travail virtuel

des forces extérieures et noté ‹Wext , car il représente le travail des forces
bi dV et ti dA sur les déplacements fictifs ‹ui (Fig. ).
 
‹Wext = bi ‹ui dV + ti ‹ui dA (7.9)
V A

Le terme de gauche de (7.8)


‹W int = ffij ‹"ij dV (7.10)
V

est appelé travail virtuel intérieur, est susceptible d’une interprétation de

même nature (Fig. ). Sur la facette dx2 dx3 d’une brique élémentaire, agit

Mécanique des Solides 64 1ère année du cycle d’Ingénieurs

 CHAPITRE 7. FORME INTÉGRALE DE L’ÉQUILIBRE ET DE LA CINÉMATIQUE

la force normale ff11 dx2 dx3 , dont le déplacement associé est ‹"11 dx1 . Ces
deux termes fournissent le travail élémentaire

(ff11 dx2 dx3 )(‹"11 dx1 ) = ff11 ‹"11 dV (7.11)

En procédant de même avec les autres composantes et en intégrant sur le

volume entier, on trouve l’expression (7.11). Le travail de l’accroissement
@ff11
@x1
dx1 donne une contribution négligeable (infiniment petit d’ordre supé-
rieur).
En conclusion, l’égalité
‹Wint = ‹Wext (7.12)

s’exprime dans le principe suivant :

Principe des déplacements virtuels ou des travaux virtuels : un solide
déformable est en équilibre si le travail virtuel extérieur est égal au travail
virtuel intérieur, pour tout champ de déplacement virtuel compatible.

Remarques

L’égalité des travaux garantit l’équilibre du solide. Le principe exprime

l’équilibre et fournit des équations d’équilibre.
Le principe des travaux virtuels est indépendant de la nature du matériau. Il
s’applique même dans le cas d’une dissipation d’énergie, e.g. en plasticité.
Le champs virtuel doit être compatible avec les appuis (Fig. ), alors l’inté-
grale de surface de l’équation (7.8) devient
   
ti ‹ui dA = ti ‹ui dA + ti ‹ui dA = ti ‹ui dA (7.13)
A At Au At

1ère année du cycle d’Ingénieurs 65 Mécanique des Solides

 CHAPITRE 7. FORME INTÉGRALE DE L’ÉQUILIBRE ET DE LA CINÉMATIQUE

Mécanique des Solides 66 1ère année du cycle d’Ingénieurs

Chapitre 8

Introduction à la Plasticité

8.1 Introduction

Lorsqu’un matériau est sollicité jusqu’à la rupture, les essais montrent

que la contrainte de rupture ffR est une grandeur présentant de fortes fluc-
tuations pouvant même dépasser la décade pour certains matériaux et que
le mode de ruine dépend de la nature du matériau. Ainsi la rupture peut
intervenir brutalement sans déformation préalable pour les matériaux quali-
fiés aujourd’hui de fragiles, tandis qu’elle n’intervient qu’après une étape de
grande déformation permanente pour les matériaux qualifiés aujourd’hui de
ductiles. Nous savons maintenant que les matériaux fragiles rompent bru-
talement au-delà d’une certaine tension, tandis que les matériaux ductiles
s’écoulent plastiquement avant de rompre sous cisaillement. Si la rupture
est toujours l’étape ultime de la ruine des structures, elle est précédée d’une
étape de plastification pour les matériaux ductiles.

8.2 Comportement plastique en sollicitation uni-

axiale

Le comportement plastique est celui d’un corps solide qui prend des
déformations permanentes sans se fissurer. On admet généralement que ces
déformations permanentes se produisent au-delà d’un seuil de contrainte
appelé seuil de plasticité ou limite d’élasticité (§ 8.2.1).

8.2.1 Seuil de la plasticité

En dessous du seuil de plasticité les déformations ne restent pas toujours

parfaitement élastiques (frottements internes), mais après décharge, le corps
doit reprendre progressivement mais intégralement sa forme initiale. Le seuil
de plasticité doit être vu comme une schématisation commode : au-delà du

67
 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.1 – Seul de plasticité

seuil les déformations permanentes deviennent beaucoup plus importantes

que les déformations élastiques. La définition du seuil est conventionnelle
(Fig. 8.1) et dépend de la précision des moyens de mesure utilisés. La limite
d’élasticité vraie (point A) définie par la valeur de la contrainte ffV qui pro-
duit une déformation permanente dans le premier grain de la microstructure
du matériau est inaccessible à la mesure expérimentale. Lorsque le nombre
de grains plastifiés augmente la déformation permanente devient accessible
à la mesure. En général, la résistance à la déformation plastique dépend
plus ou moins, selon les matériaux et la température, de la vitesse de défor-
mation (visco-plasticité) et de l’histoire des déformations permanentes déjà
subies par le matériau soit au cours du processus d’élaboration, soit au cours
des sollicitations antérieures (écrouissage).
On constate expérimentalement que lorsqu’un solide est soumis à une
contrainte non isotrope, il possède vis-à-vis d’un nouveau chargement un
seuil égal à la plus grande valeur de la charge antérieure. Ainsi, sur la figure
(8.1) représentant la courbe de charge (force F , déplacement u), le seuil
de plasticité à la première charge OB correspond au point A. A la seconde
charge CB 0 , le seuil de plasticité correspond au point B 0 . L’absence de coïn-
cidence rigoureuse entre B et B 0 et la boucle décharge-recharge BCB 0 sont
la conséquence d’un effet de frottement interne. Pour des chargements infi-
niment lents, les frottements internes n’interviennent pas, les points B et B 0
et les trajets BC et CB’ sont confondus, la plasticité est alors indépendante
du temps, absence d’effet visqueux. Dans ce cas, le nouveau seuil de plasti-
cité au point B 0 confondu avec B est bien égal à la valeur de la plus grande
charge atteinte lors de la charge précédente, point B. C’est cette évidence
expérimentale qui justifie la notion de seuil d’écoulement.
Grace à leur ductilité, les matériaux plastifiables peuvent être moulées sous
l’action de contraintes suffisantes et gardent ensuite leur forme sous des
contraintes plus faibles. Tous les métaux non alliés et certains alliages pré-
sentent une phase de comportement plastique. La déformation plastique est
beaucoup utilisée dans certains procédés de déformation plastique, e.g. for-

Mécanique des Solides 68 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

geage, estampage, emboutissage, laminage, filage, pliage, etc

Les terres, l’argile et la plupart des pâtes sont plastifiables. Enfin, les corps
réputés fragiles, e.g. bétons, roches, verres, etc, sous pression hydrostatique
suffisamment élevée cessent d’être fragiles et acquièrent un comportement
plastique vis-à-vis des contraintes additionnelles.

8.2.2 Influence des déformations plastiques sur les propriétés

mécaniques

Déformation plastique modérées

Des déformations plastiques modérées n’entraînent pas de changement

des modules d’élasticité dans les corps compacts. Leurs trois effets essentiels
sont :
- Durcissement : Elévation du seuil de plasticité pour des contraintes de
même nature et de même signe que celles qui ont produit l’écoulement
plastique. La courbe de charge (sans décharge) ff = f(") est une
fonction non décroissante de " traduisant le durcissement du matériau
écrouissage ou work-hardening (Eq. 8.1)

@ff
>0 (8.1)
@"
- Réduction de la ductilité : C’est la diminution de l’allongement pos-
sible avant la rupture. La position du point de rupture R n’est quasi-
ment pas modifiée par les décharges AO2 , BO3 , ... Ces décharges ont
pour effet de déplacer l’origine des déformations de O en O2 puis O3 ,
... réduisant progressivement l’allongement restant à parcourir jus-
qu’à la rupture. On observe généralement une diminution progressive
du durcissement au fur et à mesure que l’on se rapproche du point de
rupture R. La courbe de charge ff = f(") est une fonction croissante
de " (Eq. 8.1). Mais cette dérivée décroit.
- L’effet Baushinger : On observe généralement une diminution (en va-
leur absolue) du seuil de plasticité pour des contraintes de même
nature mais de signes opposés à celles qui ont produit la déformation
plastique (Eq. 8.2).
jff20 j <jff10 j (8.2)

Dans tous les cas on observe une réduction de l’amplitude du domaine élas-
tique (Eq. 8.3).
jff20 j`jff10 j < ff2 ` ff1 (8.3)

Dans le cas des métaux pour lesquels ff10 = `ff1 , on peut écrire

ff20 > `ff2 (8.4)

1ère année du cycle d’Ingénieurs 69 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.2 – Durcissement, ductilité et effet Bauschinger

Déformations plastiques importantes

Des déformations plastiques importantes ont pour effet d’étirer et d’orien-

ter les grains cristallins ou les chaines de polymères. Sous l’effet du lami-
nage intense, le matériau initialement isotrope devient anisotrope entraînant
le changement de ses modules d’élasticité. Les métaux fortement plastifiés
prennent une texture, les roches une schistosité.
Les modifications produites par les déformations plastiques dans tous les cas
être rattachées à un changement de structure à une certaine échelle. Dans
certains cas ces changements sont définitifs, mais il arrive assez souvent que
le corps ait tendance à reprendre très lentement sa structure initiale. Dans
le cas des métaux ce phénomène est appelé restauration, i.e. après décharge
d’une éprouvette, plus on attend pour la recharger plus l’écrouissage s’es-
tompe car le durcissement disparaît par annihilation lente des dislocations.

8.2.3 Comportement plastique pur

Dans la pratique la vitesse de déformation est suffisamment lente pour

que l’on puisse négliger les effets de frottement interne, i.e. absence de
déformation différée ou de recouvrance retardée. La plasticité pure est indé-
pendante du temps.

L’élasto-plasticité parfaite

Elle correspond au cas idéal où le seuil d’écoulement est une constante

matérielle indépendante de la valeur de la déformation plastique atteinte.
L’écrouissage ou le durcissement est nul. Le comportement plastique parfait
est schématisé par le modèle d’un patin idéal à frottement sec (Fig. 8.3)
dont le seuil de glissement correspond au seuil d’écoulement ffe indépendant
de "P . Au-delà du seuil, quand ff – ffe le patin n’oppose aucune résistance
additionnelle au glissement. Ce modèle ne permet pas de définir la déforma-

Mécanique des Solides 70 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.3 – Plasticité parfaite : Modèle du patin à frottement sec

tion plastique, il est par contre très utile pour représenter la notion de seuil
d’écoulement plastique libre (non contenu) le Yield point. L’élasticité est
représentée par le ressort de raideur E monté en série. Tant que la contrainte
ff < ffe , le patin est bloqué, seul le ressort travaille et la relation ff = ff(")
s’écrit suivant l’équation (Eq. 8.5).

ff = E"e (8.5)

La déformation totale est d’origine purement élastique jusqu’au point A.

Lorsque l’on cherche à augmenter la contrainte ff appliquée au-delà de
ffe , le patin se met à glisser sans opposer de résistance, le matériau se
plastifie, la contrainte reste plafonnée à la valeur ffe et le matériau s’écoule
plastiquement. Sa déformation totale (Eq. 8.6).

" = "e + "p (8.6)

C’est la somme de la déformation élastique constante "e = Effe et de la

déformation plastique "p dont la loi d’évolution est indéterminée dans ce
modèle élémentaire. Si l’on procède à une décharge à partir de l’état de
déformation atteint au point B le patin se bloque instantanément en arrêtant
l’écoulement plastique de sorte que "p = "p (B) tandis que le ressort se
détend progressivement et la déformation élastique suit la loi "e = Eff (ff <
ffe ) de sorte qu’après décharge complète la déformation permanente est
"p = "p (B). Un nouveau chargement retendra le ressort jusqu’à revenir au
point B au-delà duquel l’écoulement plastique peut reprendre à ff = ffe
fixé. Le fait que le seuil d’écoulement soit indépendant de l’amplitude de la
déformation plastique ›p traduit bien l’absence d’écrouissage du matériau.
En élasto-plasticité parfaite il n’existe pas de relation biunivoque entre ff et
"p .

L’élasto-plasticité écrouissante

Le modèle le plus simple consiste à mettre en parallèle sur le patin un

ressort de raideur K appelée module d’écrouissage (ou module plastique)
qui correspond au mécanisme de création des contraintes internes nécessaires

1ère année du cycle d’Ingénieurs 71 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

pour contenir l’écoulement plastique. Tant que la contrainte ff < ffe , le

patin est bloqué, seul le ressort E travaille et la relation ff = ff(") s’écrit
ff = E"e . La déformation totale est d’origine purement élastique jusqu’au
point A.
Lorsque l’on cherche à augmenter la contrainte ff appliquée au-delà de ffe ,
le patin se met à glisser opposant toujours la résistance ffe , le matériau se
plastifie mais le ressort K vient limiter l’écoulement possible de sorte que
ff = ffe +K"p est la somme de la contrainte seuil du patin et de la contrainte
ffI = K"p = ff ` ffe du ressort K. Le ressort élastique se déforme également
selon la loi d’élasticité ff = E"e . La déformation totale " = "e + "p au point
"
B est alors la somme de la déformation élastique "e = et de la déformation
E
ff ` ffe
plastique "p = de sorte que la réponse ff = f(") s’écrit suivant l’équation
K
(8.7).
1 1 1
"=( + )ff ` ffe (8.7)
E K K
soit
1
"(B) ` "(A) = (ff(B) ` ff(A) (8.8)
E0
avec
1 1 1
0
= + (8.9)
E E K
E 0 : Pente du segment AB est appelée module tangent. Comme E >> K
alors E 0 ı K. Ce modèle simple met en évidence le caractère incrémental
des lois de la plasticité. De manière générale, il y’aura une relation biuni-
voque ff = ff(") dont la linéarisation locale au voisinage d’une point B de
la zone plastique conduit immédiatement à la définition du module tangent
(Eq. 8.10).
dff
E0 = (8.10)
d"
d’où son appellation et à l’écriture incrémentale (Eq. 8.11).

dff = K d"p (8.11)

de la loi de plasticité, le module tangent E 0 (ff) et le module plastique K(ff)

étant deux des fonctions décroissantes de la contrainte ff rendant compte de
la diminution progressive du durcissement lors de la plastification.

Mécanique des Solides 72 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.4 – Plasticité écrouissante : Modèle du patin à frottement sec

8.3 Domaine d’élasticité en sollicitation quelconque

8.3.1 Généralités

L’essai uni-axial ne fournit que le cas particulier des sollicitations réelles

traction-compression que peut subir un élément de volume d’un milieu continu.
Une sollicitation quelconque de cet élément de volume est caractérisée par
le tenseur de contraintes ff, élément de l’espace vectoriel des contraintes à
six dimensions noté R6 (ffij = ffji avec j 2 f1; 2; 3g). A l’aide des résultats
expérimentaux, on aboutit à l’existence d’un domaine d’élasticité initial pour
tout trajet de charge en un point strictement à l’intérieur du matériau.
Par rapport à la sollicitation uni-axiale, le domaine d’élasticité est l’inter-
valle [`ff0 ; ff0 ] tracé sur l’axe des contraintes ff. Autrement dit, cet intervalle
est l’intersection du domaine d’élasticité initial avec l’axe des sollicitations
traction-compression. De même, dans le cas d’une sollicitation plus com-
plexe telle que celle réalisée dans l’expérience de traction-compression et
de torsion du tube mince (Fig. 8.5), où les seules composantes non nulles
du tenseur des contraintes sont la contrainte normale ffzz et la contrainte
tangentielle ffz„ variant indépendamment l’une de l’autre (Fig. 8.5), on ob-
tient le domaine d’élasticité formé par l’intersection du domaine défini dans
l’espace R6 avec le sous-espace correspondant à la sollicitation appliquée,
i.e. le plan d’équation (8.12).

ffzr = ffrr = ffr„ = ff„„ = 0 (8.12)

Dans le cas d’un matériau écroui, la charge peut être poursuivie au-delà
du domaine d’élasticité initial, l’état de charge ffij ayant atteint la frontière
initiale d’élasticité peut la franchir, il y a alors apparition des déformations
permanentes. On définit, en chaque point du trajet de charge, le domaine
d’élasticité actuel. Lorsque le point de charge franchit la frontière initiale
d’élasticité, il entraîne avec lui la frontière de ce domaine d’élasticité (Fig.
8.6). Dans le cas des sollicitations tridimensionnelles, le non d’écrouissage

1ère année du cycle d’Ingénieurs 73 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.5 – Eprouvette pour l’expérience de traction-compression et torsion

d’un tube mince. Exemple de domaine initial d’élasticité déterminé expéri-
mentalement (H.D. Bui, 1970)

Figure 8.6 – Domaine d’élasticité initial et actuel dans l’espace à six dimen-
sions R6 pour le matériau écroui

est donné à l’entraînement et à la déformation du domaine d’élasticité. L’ef-

fet Bauschinger n’est donc qu’un cas très particulier du phénomène d’écrouis-
sage.

Matériau parfaitement plastique

Dans le cas du matériau parfaitement plastique, le domaine d’élasticité

n’est pas modifié par l’apparition des déformations plastiques. Le point de
charge ffij ne peut pas franchir ce domaine. Les déformations plastiques ne
se produisent que si l’état de contrainte ffij ne peut pas franchir ce domaine.
Les déformations plastiques ne se produisent que si l’état de contrainte ffij
est sur la frontière d’élasticité (Fig. 8.7).

8.3.2 Quelques résultats expérimentaux concernant les mé-

taux

A partir de la frontière initiale d’élasticité il y a apparition des défor-

mations irréversibles. Ainsi, pour déterminer la frontière actuelle on doit
effectuer un déchargement à partir de l’état de charge situé au point B

Mécanique des Solides 74 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.7 – Domaine d’élasticité initial et actuel dans R6 pour le matériau

parfaitement plastique

(Fig. 8.6), puis un chargement dans une direction différente et chercher l’ap-
parition des premières nouvelles déformations irréversibles. L’appareillage
conçu par H.D. Bui pour ses expériences sur des éprouvettes métalliques,
avait une résolution de 2 10`5 . Pour la détermination de la frontière initiale
d’élasticité, il est courant et commode d’explorer le domaine élastique en
parcourant des trajets de chargement radiaux à partir de l’état naturel en
cherchant à détecter la perte de la linéarité en fonction de la contrainte.
La figure (8.8) présente les frontières d’élasticité bidimensionnelle suivant
des sollicitations de traction-compression et torsion obtenues après écrouis-
sage en compression pure d’une éprouvette en aluminium. On constate que
l’écrouissage en compression relève le seuil de plasticité en traction et aussi
le seuil de cisaillement pur, c’est l’effet d’expansion. On remarque aussi
que l’effet Baushinger contrebalance l’effet d’expansion et devient prépon-
dérant. En effet, les seuils de traction et de cisaillement pur sont réduits
après écrouissage suffisant en compression.
La figure 6.9 montre l’évolution du domaine d’élasticité dans le cas d’un
métal pour deux types de trajets de charge. On remarque l’apparition d’une

Figure 8.8 – Frontière d’élasticité successives après écrouissage en compres-

sion pour Al99:5 d’après (Bui 1970)

1ère année du cycle d’Ingénieurs 75 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.9 – Ecrouissage d’un alliage de Fer (d’après Bui 1970)

bosse au voisinage du point de charge. On peut dire que l’incrément de

déformations plastiques est perpendiculaire au domaine de chargement au
point de charge. C’est la preuve expérimentale de la loi de normalité.

8.4 Théorie de l’élasto-plasticité

Dans ce paragraphe, on présente les hypothèses de la théorie générale de

la plasticité incrémentale, avec ou sans écrouissage, des solides homogènes
isotropes.
Hypothèse 1 : Les accroissements totaux de déformation d"ij sont dé-
composables en incréments de déformation élastique d"eij et de déformation
p
plastique d"ij (Eq. 8.13).
p
d"ij = d"eij + d"ij (8.13)

Les accroissements de déformation élastique d"eij obéissent à la loi de Hooke

généralisée (Eq. 8.14).
d"eij = Sijkl dffkl (8.14)
Hypothèse 2 : Existence des surfaces de charge indépendantes du temps.
Considérons un élément de matière soumis à une sollicitation quelconque
définie par le tenseur contrainte ff (ou ffij ). La théorie de la plasticité sup-
pose l’existence d’une fonction f = f(ffij ; ¸mn ) définissant le domaine du
chargement, tel que :
p
f = f(ffij ; ¸mn ) < 0 : Domaine élastique d"ij = 0 ;
p
f = f(ffij ; ¸mn ) = 0 : Apparition des déformations plastiques d"ij 6=
0;
f = f(ffij ; ¸mn ) > 0 : Domaine inaccessible par le matériau.
¸mn sont des paramètres internes qui quantifient l’ensemble de phénomènes
irréversibles. Les expressions f = f(ffij ; ¸mn ) sont appelées fonctions de

Mécanique des Solides 76 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

charge dans l’espace des états de contraintes dans le solide. Dans le domaine
plastique cette fonction devient la surface qui est définie par l’équation (Eq.
8.15).
f = f(ffij ; ¸mn ) = 0 (8.15)

f est appelée surface de charge.

En l’absence des paramètres internes, la surface de charge devient (Eq. 8.16).

f(ffij ) = 0 (8.16)

Elle présente le critère de plasticité ou le seuil de plasticité initial. Pour

chaque valeur dans le champ des paramètres internes ¸mn correspond à
une surface de charge. Cette surface est une fonction qui dépend à la fois
du champ de contraintes et des paramètres internes. On rappelle que la
différentielle totale de f s’écrit suivant l’équation (8.17).
@f @f
df = dffij + d¸mn (8.17)
@ffij @¸mn

avec i; j; m et n sont des indices muets.

Il existe trois types de chargements :
1. Déchargement, df < 0 : L’état de contrainte est purement élastique
p
(d"ij = d¸mn = 0). Il s’agit d’un processus de déchargement (Eq.
8.18).
@f
dffij < 0
@ffij (8.18)
et f(ffij ; ¸mn ) = 0
p
2. Chargement neutre, df = 0 mais d"ij = d¸mn = 0.

@f
dffij = 0
@ffij (8.19)
et f(ffij ; ¸mn ) = 0

p
3. Chargement plastique, df = 0 mais d"ij 6= 0 et d¸mn 6= 0.

@f
dffij > 0
@ffij (8.20)
et f(ffij ; ¸mn ) = 0

La représentation géométrique de f(ffij ; ¸mn ) = 0 dans l’espace des contraintes

principales est une hyper-surface délimitant l’ensemble des états de contraintes
dans le solide en élasticité et en plasticité. Le vecteur gradient de la fonction
f au point de chargement a le sens de la normale extérieure de la surface
@f
de charge. Il est donné par l’expression dffij .
@ffij
Hypothèse 3 : La déformation plastique est proportionnelle à la projection
de l’accroissement de contrainte sur la normale à la surface de charge au

1ère année du cycle d’Ingénieurs 77 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.10 – Domaine élastique convexe et loi de normalité

point de chargement.
L’incrément de déformation plastique est donné par la relation (8.21).

p @f
d"ij = ˛ij dffkl (8.21)
@ffkl
où — est un coefficient de proportionnalité. On obtient la forme générale de
la déformation plastique d’un matériau écroui qui est donnée par l’équation
(8.22).
p @h @f
d"ij = — dffkl (8.22)
@ffij @ffkl
Hypothèse 5 : Postulat de Drucker ou hypothèse de stabilité des matériaux
écrouisssables, il est annoncé comme suit "Le travail accompli lors d’un
cycle de charge complet quelconque (cycle de charge entièrement fermé
dans l’espace des contraintes) est non négatif" (§ 8.4.1).

8.4.1 Hypothèse de stabilité des matériaux, postulat de Dru-

cker

En partant des considérations générales sur le comportement des maté-

riaux inélastique, Drucker a présenté une définition générale d’un matériau
écrouissable stable, dit standard (qui exclut les matériaux instables ou non
standards). La figure (?? illustre les différents comportements des matériaux
inélastiques. Au sens de Drucker, le matériau (a) est stable car, soumis à un
´ff´"
incrément de contrainte ´ff, il fournit le travail positif. Les matériaux
2
(Fig. 8.11b et c) sont considérés comme instables car le travail est néga-
tif. En termes généraux, le postulat de stabilité de Drucker peut s’énoncer
comme suit :
Le travail accompli lors d’un cycle de charge complet quelconque (cycle de
charge entièrement fermé dans l’espace des contraintes) est non négatif.
L’hypothèse précédente entraîne des conséquences très importantes pour la
nature des déformations plastiques des matériaux écrouissables. Examinons

Mécanique des Solides 78 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

un cycle de charge particulier ABCA0 fermé dans l’espace des contraintes

0
(Fig. 8.12). L’état de contrainte initial ffij est représenté par un point A
quelconque situé à l’intérieur de la surface de charge initiale. Suite à l’aug-
mentation de la sollicitation, les contraintes atteignent la valeur ffij sur la
surface de charge (point B). En ce point, l’accroissement de contrainte dffij ,
représenté par le vecteur BC~ dirigé vers l’extérieur de la surface de charge
p
initiale, provoque des déformations plastiques d"ij et une légère modification
de la surface de charge.
Le trajet de retour CA est un déchargement élastique par rapport au nou-
veau seuil (f+df) qui admet toujours A comme un point qui se trouve à
l’intérieur. Pour ce cycle de charge, le postulat de Drucker s’exprime par
l’inégalité (8.23).
0 0
A= (ffij ` ffij )d"ij – 0 (8.23)

Où ffij sont les contraintes appartenant au trajet de charge. On peut expli-

citer cette inégalité de la manière suivante (Eq. 8.24).
  
0 0 p 0
A= (ffij `ffij )d"eij + 0
(ffij +dffij `ffij )(d"ij +d"eij )+ (ffij 0
`ffij )d"eij – 0
AB BC CA
(8.24)
En regroupant les termes semblables et en négligeant les termes du second
ordre, on trouve l’équation (8.25).
 
0 0 p
A= (ffij ` ffij )d"eij + 0
(ffij ` ffij )d"ij – 0 (8.25)
BC

Le premier terme donne l’équation (8.26).

0 0 0
(ffij ` ffij )d" = ffij d"eij ` ffij
0
d"eij (8.26)

est nul car, d’une part, le travail élastique total d’un cycle fermé est nul
et d’autre part, la déformation élastique est déterminée à chaque état de
contrainte. Comme le trajet BC peut être rendu aussi petit que l’on veut
(8.27).
0 p
(ffij ` ffij )d›ij – (8.27)

Figure 8.11 – Stabilité des matériaux au sens de Drucker

1ère année du cycle d’Ingénieurs 79 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.12 – Cycle de charge fermé dans l’espace des contraintes

Cette inégalité est encore appelée principe de Hill du travail local maximum.
Si A est très proche du point B, alors :

p
dffij d"ij – 0 (8.28)

Cette dernière inégalité est connue par la condition d’unicité car elle a pour
conséquence l’unicité au sens incrémental des solutions plastiques.

8.4.2 Propriétés

Loi de normalité

Le vecteur des accroissements de déformations plastiques est normal à la

surface de charge et dirigé vers l’extérieur (Eq. 8.29).

p @f
d"ij = d– (8.29)
@ffij

Avec d– – 0 est un facteur de proportionnalité.

Convexité

La surface de charge f = 0 est une fonction convexe.

D’après ce qui précède, on peut écrire :

@f
d– = — dffkl – 0 (8.30)
@ffkl

D’où
@f @h
= (8.31)
@ffij @ffij
Si la fonction de charge et le potentiel plastique sont identiques à une
constante près, le matériau est dit standard ou associé.

f =h (8.32)

Mécanique des Solides 80 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Pour un matériau standard, Prager a proposé la loi de comportement définie

dans l’équation (8.33).

p @f @f
d"ij = ˛— dffkl (8.33)
@ffij @ffkl

avec — = —(ffij ; ¸mn ) – 0.

@f
8
1 si f = 0 et dffkl – 0
>
>
>
>
>
>
>
< @ffkl
˛ = > 0 si f < 0 (8.34)
@f
>
>
>
1 si f = 0 et dffkl < 0
>
>
>
:
@ffkl
Cette loi de comportement est établie pour des surfaces de charge régulières.
Pour le cas de la présence des points anguleux en lesquels la fonction de
charge n’est pas dérivable, Koiter a proposé une loi constitutive associée aux
surfaces qui possèdent des points singuliers. Cette loi permet la construction
de la surface de charge par l’assemblage d’un nombre fini d’arrêts dont les
intersections sont des points de singularité (Eq. 8.35).

fr (ffij ; ¸mn ) (8.35)

La loi de normalité associée à f est donnée par l’équation (8.36).

p @fr @fr
d"ij = ˛r —r dffkl (8.36)
@ffij @ffkl

avec —r = —r (ffij ; ¸mn ) – 0.

@fr
8
1 si fr = 0 et dffkl – 0
>
>
>
>
>
>
>
< @ffkl
˛r = > 0 si fr < 0 (8.37)
@fr
>
>
>
1 si fr = 0 et dffkl < 0
>
>
>
:
@ffkl
Dans le cas des surfaces non régulières, la condition de stabilité de Drucker
p
est satisfaite lorsque l’incrément de déformation d"ij se trouve à l’intérieur
de l’hyper-cône formé par les vecteurs normaux des surfaces adjacentes.

8.4.3 Puissance de dissipation plastique

La dissipation plastique par unité de volume en un point dans un solide

déformable est définie par le scalaire définit par l’équation (8.38).

˜
D("_ij ) = max(ffij "_ij ) (8.38)
ff
_ ij

˜
où ffij est un champ de contrainte statiquement admissible.
Les propriétés de la puissance de la dissipation plastique D("_ij ) sont :

1ère année du cycle d’Ingénieurs 81 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

- D est fonction univoque de déformation, i.e. elle ne dépend que de

la façon dont le corps se déforme plastiquement et des grandeurs des
vitesses de déformation.
- est une fonction homogène d’ordre 1.

8— – 0 D(—"_ij ) = —D("_ij ) (8.39)

- D est une fonction convexe pour toute vitesse de déformation "ij

cinématiquement admissible :

D("_˜ij ) ` D("_ij ) – ffij ("_˜ij ` "_ij ) (8.40)

8.4.4 Règles de l’écrouissage

Paramètres d’écrouissage

Les règles de l’écrouissage sont décrites par les relations liant les seuils
plastiques actuels (ou limites élastiques actuelles) ffY du matériau et les
déformations plastiques accumulées au cours du chargement plastique. Pre-
mièrement, nous avons besoin de définir quantitativement les variables d’état
qui peuvent représenter la déformation plastique. Ces variables sont appelés
paramètres d’écrouissage ou variables plastiques internes souvent notées ».
Le seuil plasticité actuel ffY , le module tangent Et et le module plastique
K peuvent être exprimés comme des fonctions du paramètre d’écrouissage
».
Le paramètre d’écrouissage » peut être défini de plusieurs façons. Les plus
fréquentes sont définies comme suit :

p
» = "p = d"p d"p : Déformation plastique équivalente ;


» = Wp = ff d"p : Travail plastique (Energie de déformation plas-
tique) ;


» = "p = d"p : Déformation plastique.
On note que la déformation plastique "p est un mauvais paramètre d’écrouis-
sage car cette quantité est réversible et ne peut pas être accumulé.

Idéalisation des règles de l’écrouissage

La notion de l’écrouissage est expliquée par les relations qui lient le seuil
de plasticité actuel ffY , le module tangent Et ou le module plastique K et
le paramètre d’écrouissage ». La déformation plastique accumulée pendant
le chargement en traction (ou compression) a une influence, à la fois, sur
le seuil plastique actuel en traction et aussi en compression, c’est une autre
explication de l’effet Bauschinger. L’écrouissage peut être idéalisé suivant
les trois règles citées ci-après :
Écrouissage isotrope : Lorsque du seuil de plasticité actuel est le même

Mécanique des Solides 82 1ère année du cycle d’Ingénieurs

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

dans le cas de la traction ou de compression, on dit que l’écrouissage est

de nature isotrope (Fig. 8.13). Cette règle d’écrouissage peut être exprimée
par l’équation (8.41).
jffj =jff(»)j (8.41)

Écrouissage cinématique : Lorsque la différence algébrique entre le seuil

plastique actuel de la traction et celui de la compression reste constante,
alors la nature de l’écrouissage est de nature cinématique (Fig. 8.13). Soit
ffyt et ffyc sont les seuils plastiques actuels de traction et de compression
respectivement. Cette règle d’écrouissage peut être exprimée par l’équation
(8.42).
ffyt (») ` ffyc (») = 2ffy0 (8.42)

avec ffy0 la limite élastique initiale (premier seuil plastique) qui est suppo-
sée la même, à un signe mois près, en traction et en compression. Conven-
tionnellement, la règle de l’écrouissage cinématique peut être exprimée par
l’équation (8.43).
jff ` c(»)j = ffy0 (8.43)

avec c(0) = 0.
Écrouissage indépendant : Dans cette dernière règle d’écrouissage, les
seuils plastiques actuels de traction et de compression sont indépendants l’un
de l’autre (Fig. 8.13). Les déformations plastiques accumulés pendant l’essai
de traction excessive n’ont pas une influence sur le seuil plastique actuel de
compression et l’inverse est vrai. Dans ce cas, il faut définir deux paramètres
d’écrouissage »t et »c pour définir le comportement plastique en traction et
en compression respectivement. La règle de l’écrouissage indépendant est
exprimée mathématiquement par les équations (8.44)

ff = ffyt (»t ) pour ff > 0

(8.44)
ff = ffyc (»c ) pour ff < 0

Écrouissage mixte : Dans le cas d’un écrouissage mixte, on enregistre une

contribution à la fois de l’écrouissage isotrope et celui cinématique. i.e.
l’incrément de déformation plastique peut être décomposé en deux parties
(Eq. 8.45).
d"p = d"pi + d"pk
(8.45)
ou d" = d"e + d"pi + d"pk

avec d"pi et d"pk représentent respectivement l’incrément de déformation

plastique d’origine isotrope et celui d’origine cinématique. On suppose que
le ratio entre ces deux incréments est fixe, on peut écrire alors (Eq. 8.46).

d"pi = a d"pk
(8.46)
ou d"pk = (1 ` K)d"p

1ère année du cycle d’Ingénieurs 83 Mécanique des Solides

 CHAPITRE 8. INTRODUCTION À LA PLASTICITÉ

Figure 8.13 – Illustration des règles d’écrouissage en se basant des char-

gements cycliques uni-axiaux. Du gauche à droite : Écrouissage isotrope,
Écrouissage cinématique et Écrouissage indépendant

avec 0 » a » 1. On peut généraliser les deux premières règles d’écrouissage

en une seule équation (8.47).

jff ` c("pk )j = ffy (»i ) (8.47)

»i est un paramètre d’écrouissage qui a une relation avec "pi , c(0) = 0 et

ffy0 = ff0 .
Pour simplifier, on définit la fonction c comme suit :
 
pk pk
c(" )= K d" = (1 ` a) K d"p (8.48)

jff ` c("pk )j = ff0 + a(ffy (») ` ff0 ) (8.49)

Pour a = 1, cette équation traduit la règle de l’écrouissage isotrope (8.50).

jffj = ffy (») (8.50)

et pour a = 0, elle traduit la règle de l’écrouissage cinématique (8.51).

jff ` c("pk )j = ff0 (8.51)

Mécanique des Solides 84 1ère année du cycle d’Ingénieurs

Bibliographie

[1] Livre : François Frey, Analyse des structures et milieux continus - Méca-
nique des structures, Volume 2, Édition : Presses polytechniques et uni-
versitaires romandes, ISBN 2-88074-4-434-2, Deuxième édition, 2006,
Lausanne.
[2] Livre : François Frey, Analyse des structures et milieux continus - Méca-
nique des solides, Volume 3, Édition : Presses polytechniques et univer-
sitaires romandes, ISBN 2-88074-358-3, Deuxième édition, 2006, Lau-
sanne.
[3] Livre : Serge Étienne, David Laurent, Emile Gaudry, Philippe Lagrange,
Julian Ledieu, Jean Steinmetz, Les Matériaux, Édition DUNOD, ISBN
978-2-10-051524-0, 2008, Paris.
[4] Livre : Charles MASSONET et Serge DESCOTTO, Mécanique des Ma-
tériaux, Édition De Boeck , ISBN 978-2-8041-2021-4, 1994, Bruxelles.
[5] Livre : Catherine Potel, Principe et Applications de Mécanique Analy-
tique, Édition Cépaduès, ISBN 2-85428-742-8, 2006, Toulouse.
[6] Livre : Milan Jirasek and Zdenek P. Bazant, Inelastic Analysis of Struc-
tures, Édition JOHN WILEY & SONS LTD, ISBN 0-471-98716-6, 2002,
England.
[7] Livre : Jean Salençon , De l’Élasto-plasticité au calcul de la rupture,
Édition polytechnique, ISBN 978-2-7302-0915-1, 2002, Paris.
[8] Article : Jean Lemaitre, Petite Histoire de l’Expérimentation en Méca-
nique des Solides, Journal : Meccanica, Edition : Springer, ISSN Print :
0025-6455 (Print), ISSN Online 1572-9648, Volume 36, Issue 1, pp
13–35, January 2001.
[9] Livre : François Sidoro, Mécanique des milieux continus, Engineering
school. HAL Id : cel-00530377, École Centrale de Lyon, 1980, Lyon.

85