1484

SEMESTRE 3
2006-2007
MODULE F312 :
ELASTICITE LINEAIRE
PROBLMES et CORRIGES
Bienvenue vous
au laboratoire de :
Dimensionnement Des Structures
du Dpartement :
Gnie Mcanique et Productique
Ce livre lectronique est destin complter le cours enseign durant la deuxime
anne du module F312, et relatif au Troisime Semestre.
Il reprend le plan suivi en amphithtre avec davantages de dtails, d'illustrations
ainsi que des corrigs des Problmes du fascicule Travaux Dirigs qui, nous lesprons,
vous permettront de mieux comprendre cette matire qui nest pas si terrible quelle
peut laisser paratre.
A prsent, choisissez sur votre gauche dans longlet signet un chapitre du programme
que vous dsirez voir ou revoir .
BON TRAVAIL-BON COURAGE
G.LHERMET-G.VESSIERE
TD-2A-S3-F312
F312
CONTRAINTES et CRITERES
DEFORMATIONS
LOIS DE COMPORTEMENT
PROBLEMES DE SYNTHESE
CONTACT DE HERTZ
IUTB-LYON1-GMP-DDS
Page 2
TD-2A-S3-F312
1. CONTRAINTES ET CRITERES
PROBLME N1
En un point P, situ sur la surface extrieure dun cylindre soumis un moment de
torsion C, l'tat de contraintes est dfini par la matrice plane suivante, o est la
contrainte de torsion :
z
P
G
[ ]xyz
0 0
= 0 0
0 0 0
1) Reprsenter l'tat de contraintes sur les faces d'un cube infiniment petit isol en P
dans les axes xyz.
2) Dterminer les lments principaux des contraintes.
3) Tracer le tri-cercle de MOHR. En dduire les cisaillements principaux en module.
RPONSES N1
C
y
P
z
z
IUTB-LYON1-GMP-DDS
Page 3
TD-2A-S3-F312
Composante suivant y du
vecteur contrainte x
yx=
Facette x
y
45
La facette z est principale, et la contrainte principale associe est nulle. (Bord libre)
Cest donc un tat de contraintes planes en un point situ la surface dune poutre.
Les lments principaux des contraintes nous sont donns par :
Y =
X
1
2 + 4 2
2
do : X = Y =
avec la contrainte normale =0
[ ]
P
xyz
xx xy xz
= yx yy yz
zx zy zz
[ ]
P
xyz
r r
X
avec = x , X =45
0 0
0
P
[ ]XYZ = 0
= 0 0
Les axes principaux sont positionns par tan =
0 0 0
0
0
0
th t
X = Y =
Z =
/2
X
mY(-)
IX
P IZ
/2
Y
mZ(0)
mX()
IY
hh
CY
CX
CZ
IUTB-LYON1-GMP-DDS
Page 4
TD-2A-S3-F312
PROBLME N2
Considrons l'tat plan de contraintes
au point P reprsent sur la figure cicontre. Les contraintes sont en MPa.
28
1) Ecrire la matrice des contraintes
36
h
en P dans le repre (xyz).
36
2)
Dterminer
les
lments
30
principaux des contraintes.
7
7
x
P
3) On se place dans le plan des
contraintes (plan xy).
36
a) Calculer les composantes du
36
vecteur contrainte, ainsi que les
contraintes normale et tangentielle
28
s'exerant sur un plan de coupe dont
la normale fait un angle de 30 par rapport laxe x.
b) Dterminer les normales aux plans de coupe sur lesquels sexercent les contraintes
de cisaillement maximum. En dduire la valeur de ce cisaillement et la valeur de la
contrainte normale correspondante.
4) Reprendre le problme avec le cercle de MOHR.
RPONSES N2
y
28
36
[ ]Pxyz
xx
= yx
zx
xy
yy
zy
xz
yz
zz
[ ]
P
xyz
7 36 0
= 36 28 0 MPa
7
0 0 0
36
P
36
36
28
Etat plan de contraintes : la facette z est donc principale et Z=0. Les 2 autres valeurs
propres sont donnes par : YX =
xx + yy
2
1
2
La direction propre X est calcule par tan =
IUTB-LYON1-GMP-DDS
yy ) + 4 xy2
2
xx
X=55MPa
r r
X xx
avec = x , X
xy
Y=-20MPa
)
Page 5
TD-2A-S3-F312
r r
55 7
tan =
do = 5313 = x , X
36
[ ]PXYZ
55
20
5313
0
0
55
= 0 20 0 MPa
0
0
0
P
20
55
En utilisant la relation de CAUCHY :
[ h ]/ xyz
[ h ]/ xyz
7 cos 30 + 36 sin 30
= 36 cos 30 + 28 sin 30
[ h ]/ xyz
[ h ] = [ ][h]
y
t
7 36 0 cos 30
= 36 28 0 sin 30
0 0 0 0
28
tt
ht
36
yh
xh 24.062
= yh = 45.177 MPa
zh 0
th
hh
xh
30
r
r
h = 24,062 x + 45,177 y
36
7
7 36 cos 30
Compte tenu des relations : hh = [h ][ ][h ] = [cos 30 sin 30 0] 36 28 0 sin 30
0 0 0 0
t
7 36 sin 30
th = [t ][ ][h] = [ sin 30 cos 30 0]36 28 0 cos 30
0 0 0 0
et
do : hh = 7 cos 2 30 + 28 sin 2 30 + 2.36 sin 30 cos 30
th = (7 28)sin 30 cos 30 + 36(cos 2 30 sin 2 30)
IUTB-LYON1-GMP-DDS
hh=43,42MPa
th=27,09MPa
Page 6
TD-2A-S3-F312
Compte tenu de ltat plan de contraintes, il est possible dutiliser les relations :
xx + yy
hh =
xx yy
+
cos 2 + xy sin 2
2
2
yy
th = xx
sin 2 + xy cos 2
2
t
27,09
8,42
qui mnent au mme rsultat.
27,09
43,42
30
x
[ ]
P
htz
43,42 27,09 0
= 27,09 8,42 0 MPa
0
0
0
43,42
27,09
27,09
8,42
tt=peut etre calcul avec la relation donnant hh avec langle +90, ou laide du
premier invariant linaire de la matrice : I1 = tr [ ]
tt=-8,42MPa
Les facettes sur lesquelles sexercent les cisaillements maximums sont les plans
bissecteurs des didres principaux. La valeur de ces cisaillements principaux est gale, en
module, la demie diffrence des contraintes principales dans chacun des plans principaux.
La valeur de la contrainte normale associe cette facette est la demie somme.
= 45
max =

[ ]z
P
55 + 20
= 37.5MPa
2
55 20
=
= 17.5MPa
2
17,5 37,5 0
= 37,5 17,5 0 MPa
0
0
0
IUTB-LYON1-GMP-DDS
y
37,5
17,5
37,5
17,5
17,5
x
37,5
17,5
Page 7
TD-2A-S3-F312
Rsum des diffrents rsultats obtenus.
t
8,42
Y
55
17,5
28
37,5
36
27,09
27,09
5313
45
20
30
P
43,42
37,5
17,5
36
7
Reprendre le problme avec le cercle de MOHR :

Cest une reprsentation plane du faisceau des contraintes. Sa construction ncessite la
connaissance de 2 vecteurs contraintes sexercant sur 2 facettes orthogonales.
Dans notre cas les donnes du problme sont:
[ ]
P
xyz
7 36 0
= 36 28 0 MPa
0 0 0
Les facettes x et y sont orthogonales, et nous devons tout retrouver laide du cercle
partir de la connaissance des trois valeurs :
xx= 7 MPa
IUTB-LYON1-GMP-DDS
yy= 28 MPa
xy= 36 MPa
Page 8
TD-2A-S3-F312
[ x ]/ xyz
1)
2)
3)
4)
5)
7
= 36 MPa y
0
[ ]
/ xyz
36
= 28 MPa
0
Tracer les axes hh et th et adopter une chelle 1cm=10MPa

Positionner les points mx (7,36) et my (28,-36).
Tracer le diamtre mxmy du cercle. Son intersection avec laxe des abscisses donne le
centre du cercle IZ.
Tracer le cercle CZ de centre IZ et de rayon IZmx.
Les intersections du cercle CZ avec laxe des abscisses donnent les points mX et mY qui
matrialisent les facettes principales du plan des contraintes. OmX = 55 = X
OmY = 20 = Y
6)
2 = 10626
Toutes les autres rponses sont obtenues en positionnant les points sur le cercle qui
matrialisent les facettes associes et en mesurant les valeurs correspondantes,
compte tenu de lchelle.
th
36
mx(7,36)
m (17.5,+37.5)
mh(43.42,27.09)
60
Z
10626
mY(-20)
O
mZ(0)
hh
28
7
IZ
mX(55)
90
y
mt(-8.42,27.09)
CZ
-36
IUTB-LYON1-GMP-DDS
my(28,-36)
m (17.5,-37.5)
1 cm
5 MPa
Page 9
TD-2A-S3-F312
PROBLME N3
P
G
Soit un arbre section circulaire de

diamtre 42mm. La section la plus sollicite
est soumise un moment de flexion M de 800
Nm et un moment de torsion C de 600 Nm.
1) Calculer la matrice des contraintes au
point le plus sollicit ( > 0) dans le repre
xyz.
2) Dterminer les lments principaux des
contraintes.
3) Calculer le cisaillement maximum en
module.
4) Reprendre le problme avec le cercle de
MOHR.
RPONSES N3
C=600Nm
M=800Nm
y
P
Isolons le point P de la section G soumise

de la flexion-torsion, dans des axes xyz
parallles aux axes G,xyz de la thorie
des poutres.
La facette z est perpendiculaire au bord libre xy.
D=42mm
IUTB-LYON1-GMP-DDS
Page 10
TD-2A-S3-F312
Compte tenu du module F213 :

Cours sur la flexion : xxP =
+ M Y P 800 .10 3 64 42
Z =
= 109,98 MPa
IY
2 42 4
Cours sur la torsion : P max i =
C
600.103 32 42
R=
= 41,24 MPa
IG
2 42 4
x
La matrice des contraintes au point P est donc :
109,98
41,24
[ ]
P
xyz
109,98 41,24 0
= 41,24
0
0 MPa
0
0
0
41,24
41,24
41,24
Elments Principaux des contraintes :
109,98
La direction z est principale et Z=0.

Etat plan de contraintes : Les 2 autres valeurs propres sont donnes par :
Y =
xx + yy
2
1
2
X=123,74MPa
yy ) + 4 xy2
2
xx
Y=-13,74MPa
La direction propre X est calcule par

123,74
(r )
r
= x , X
xx
tan = X
xy
1844
(r )
r
= 1844 = x , X
[ ]PXYZ
0
0
123,74
= 0
13,74 0 MPa
0
0
0
IUTB-LYON1-GMP-DDS
13,74
P
13,74
123,74
Page 11
TD-2A-S3-F312
Le cisaillement maximum en module : thmax i = Z =
X Y
2
123,74 + 13,74
= 68,74MPa
2
Cercle de MOHR :
La construction, et la mthodologie est identique au problme prcdent.
m(55,68.74)
th
mx(109.98,41.24)
41.24
Z
3688
O
109.98
mZ(0)
mY(-13.74)
IZ
hh
mX(123.74)
90
my(0,-41.24)
CZ
m(55,-68.74)
IUTB-LYON1-GMP-DDS
Page 12
TD-2A-S3-F312
PROBLME N4
y
50
100
Considrons l'tat plan de contraintes au

point P reprsent sur la figure ci-contre.
100
Les contraintes sont en MPa.
x
50
50
P
1) Ecrire la matrice des contraintes en P
dans le repre xyz.
100
2) Dterminer les lments principaux des
100
contraintes. Reprsenter les lments
50
principaux des contraintes dans le plan xy
3) Tracer le tri-cercle de MOHR. En dduire la valeur du cisaillement maximum. Dans
quel plan se trouvent les normales relatives aux plans de coupe soumis au cisaillement
maximum ?
RPONSES N4
50
Matrice des contraintes en P dans le repre xyz :
[ ]
P
xyz
100
100
50 100 0
= 100 50 0 MPa
0
0
0
50
50
100
100
50
Elments principaux des contraintes :
X=50MPa
150
(r )
r
[ ]
P
XYZ
150
Z=0
= 45 = x , X
50
45
Y=-150MPa
0
0
50
= 0 150 0 MPa
0
0
0
50
IUTB-LYON1-GMP-DDS
Page 13
TD-2A-S3-F312
Tricercle de MHOR :
Les trois cercles CX, CY, et CZ sont construits partir des 3 valeurs propres de la matrice des
X=50MPa
Y=-150MPa Z=0
Positionner les points mX(50), mY(-150), et mZ(0), et tracer les 3 cercles.
contraintes.
th
mx(-50,100)
Z
X
90
O
mY(-150)
mZ(0)
IZ
IX
IY
hh
mX(50)
90
CZ
my(-50,-100)
La valeur du cisaillement maximum est donne par le rayon du plus grand des 3 cercles :
max = Z =
50 + 150
= 100MPa
2
Les normales relatives aux plans de coupe soumis au cisaillement maximum sont les bissectrices du plan
principal XY. = 45 .Ce sont donc les facettes x et y.
IUTB-LYON1-GMP-DDS
Page 14
TD-2A-S3-F312
PROBLME N5
Considrons les trois tats plans de contrainte ci-dessous. Les axes sont principaux. Les
contraintes sont en MPa.
1) Ecrire les matrices relatives ces trois tats de contrainte.
2) Daprs le critre de VON MISES, quel est ltat de contrainte le plus dangereux ?
Y
700
Y
600
400
100
-400
-100
Cas (b)
Cas (a)
Cas (c)
RPONSES N5
Y
Y
700
600
400
100
Cas (a)
[ ]
P(a)
XYZ
100
Cas (b)
400
Cas (c)
100 0 0
100 0 0
400 0 0
P (b )
P(c)
600 0 MPa ; [ ]XYZ = 0

400 0 MPa
= 0 700 0 MPa ; [ ]XYZ = 0
0
0
0
0 0
0 0
0 0
La contrainte
VM
VM
2
2
qu
en 2D compte tenu que dans les 3cas Z=0 scrit : qu = X + Y X Y
VM
qu
= 100 2 + 700 2 100.700
VM ( a )
qu
= 655,74 MPa
VM
qu
= 100 2 + 600 2 + 100.600
VM ( b )
qu
= 655,74 MPa
VM
qu
= 400 2 + 400 2 + 400.400
VM ( c )
qu
= 692,82 MPa
Ltat de contrainte le plus dangereux, au sens de VON MISES, est donc le cas (c) .
IUTB-LYON1-GMP-DDS
Page 15
TD-2A-S3-F312
PROBLME N6
P
G
Un arbre de diamtre D est soumis un

moment de flexion M de 117.81 Nm et un
moment de torsion C de 157.08 Nm. Les limites
lastiques en traction et compression du
matriau sont identiques et valent 750 MPa.
1) Dterminer les expressions de la
contrainte maximum de flexion et de la
contrainte maximum de torsion en fonction de
D3.
2) En appliquant le critre de VON MISES,
calculer (en mm) la valeur du diamtre D pour
que lon ait un coefficient de scurit de 3.
RPONSES N6
C=157,08Nm
M=117,81Nm
y
P
G
y
:D
z
z
La formulation de dpart est identique celle du Problme 3 :
IUTB-LYON1-GMP-DDS
Page 16
TD-2A-S3-F312
Cours sur la flexion :
xxP =
Cours sur la torsion :
+ M Y P 117 ,81 .10 3 64 D 1, 2 6

Z =
= 3 10
IY
2 D4
D
max i
C
157 ,08 .10 3 32 D 0,8 6
=
R=
= 3 10
IG
2 D4
D
x
La matrice des contraintes au point P est donc :
[ ]
P
xyz
1,2 6
10
3
0 D
0,8
= 0 0 = 3 106
D
0 0 0 0
0,8 6
10
D3
0
0
La contrainte quivalente de VON MISES pour une poutre :

VM
qu
= 2 + 3 2
Et le coefficient de scurit :
3=
750
2
0,8 6
1,2 6
3 10 + 3 3 10
IUTB-LYON1-GMP-DDS
D=3
e
e
=
VM
2
qu
+ 3 2
3
750
(1,2.10 ) + 3(0,8.10 )
6 2
6 2
D= 20 mm
Page 17
TD-2A-S3-F312
PROBLME N7
Diamtre D
Figure 1
N
G
Figure 2
Un arbre de diamtre D = 20 mm est soumis suivant son axe un effort de traction de

62.83 kN et un moment de torsion de 157.08 mN (figure 1). Les limites lastiques en
traction et compression du matriau sont identiques et valent 400 MPa.
1) Calculer (en MPa) la contrainte de traction et la contrainte maximum de
torsion.
2) En appliquant le critre de VON MISES, calculer le coefficient de scurit.
3) L'effort de traction est excentr de x par rapport laxe du cylindre
(figure 2). Le cylindre est donc soumis de la torsion, de la traction et de la flexion.
Calculer, en appliquant toujours le critre de VON MISES la valeur limite de la
contrainte de flexion au-del de laquelle le matriau devient plastique.
En dduire :
a)
La valeur du moment de flexion (en Nm).
b)
La valeur de x (en mm) x (en mm).
RPONSES N7
La contrainte de traction due leffort normal :
max i
N
62,83.103 4
= x =
= 200 MPa
A
20 2
157,08 Nm
G
62,83 kN
La contrainte de torsion due au couple :
max i
C
157,08.103 32 20
R=
= 100 MPa
IG
2 20 4
IUTB-LYON1-GMP-DDS
Page 18
TD-2A-S3-F312
La contrainte quivalente de VON MISES pour une poutre :

VM
qu
= 2 + 3 2
Et le coefficient de scurit :
Do :
400
(200)2 + 3(100)2
e
e
=
VM
2
qu
+ 3 2
= 1,5
Dans cette deuxime partie leffort de traction est excentr de x par rapport laxe du
cylindre. La section G est donc soumise de la traction, de la flexion et de la torsion.
62,83 kN
La contrainte normale est donc la superposition

de la contrainte de traction et de celle de flexion
157,08 Nm
x
= flexion + traction
=
e
e
=
VM
qu
2 + 3 2
En ce plaant la limite le coefficient de scurit vaut 1.
1=
400
+ 200 ) + 3(100)
2
flexion
flexion = 160,55MPa
Nous pouvons en dduire le moment de flexion :
Pflexion =
M flexion
IY
ZP
M flexion =
M flexion = N x x
IUTB-LYON1-GMP-DDS
160,55 20 3
= 126100 Nmm
32
x=
126100
= 2mm
62830
Page 19
TD-2A-S3-F312
PROBLME N8
P
G
Une barre section circulaire de rayon r est

soumise simultanment un couple de torsion C
et un moment de flexion M.
1) Calculer les contraintes de flexion et de
torsion au point P en fonction de M, C et r.
2) Donner la relation qui lie M et C lorsque le
matriau travaille la limite en appliquant le
critre de VON MISES. Mettre la relation sous
M2
C2
la forme
+
= 1 en prcisant les
a2
b2
expressions de a et b en fonction de r et e.
3) Reprsenter la relation prcdente en
fonction de R dans des axes M,C.
4) Calculer les valeurs limites de M et C (en
Nm) pour les deux critres dans les conditions
suivantes : M = C, r = 10 mm, e = 500 MPa
RPONSES N8
C
x
y
P
G
Rayon : r
IUTB-LYON1-GMP-DDS
Page 20
TD-2A-S3-F312
Cours sur la flexion : xxP =
+ M Y P M 4 r 4M
= 3
Z =
r4
r
IY
Cours sur la torsion : P max i =
C
C 2 r 2C
R=
= 3
IG
r4
r
x
4M
r 3
=4M/ r3
=2C/ r3
2C
r 3
=2C/ r3
y
=2C/ r3
=2C/ r3
=4M/ r3
Si le matriau travaille la limite, le coefficient de scurit vaut 1. Appliquons le critre de

VON MISES dans le cas des poutres :
e
e
=
VM
2
qu
+ 3 2
1=
4M
2C
2
3 + 3 3 = e
r
r
( )
+ 3( )
M2
r 3 e
1=
C2
r 3 e
+
2 3
e
2C
4M
3 + 3 3
r
r
r 3 e
4
=1 ;
Nous en tirons les expressions de a et b en fonction de r et e :
a=
b=
M2
C2
+
=1
(a )2 + (b )2
r 3 e
2 3
Cest lquation dune ellipse de grand et petit axe a et b:
IUTB-LYON1-GMP-DDS
Page 21
TD-2A-S3-F312
CONTRAINTES
1. Vecteur contrainte en un point P :
1.1. Dfinition :
Considrons un solide en quilibre sous l'action d'un torseur extrieur et un point P
l'intrieur du solide.
Plan de coupe
Effectuons une coupe suivant un plan quelconque

passant par le point P en dehors des points
d'application du chargement et isolons l'une des
deux parties rsultant de la coupe. Soit h la

normale au plan de coupe (la normale tant dirige
vers l'extrieur de la partie isole).
f
h
La partie isole est en quilibre sous l'action d'une partie
du torseur extrieur et des forces intrieures f

s'exerant sur des lments de surface d'aires A
entourant chaque point P de la surface sectionne. Ces
forces lastiques inconnues sont les forces de cohsion
qui remplacent l'action de la partie retire.
Par dfinition, le vecteur contrainte au point P relatif au plan de coupe dont la
normale extrieure (au solide) est h , est dfini par :

z
df=hdA
P
r
r
f
h ( P ) = lim
A
A 0
r
df
=
dA
dA
x
A 0
Unit de mesure :
La norme du vecteur contrainte tant homogne une pression (ou une tension),
lunit SI de contrainte est le Pascal (1Pa = 1N / m). Le m tant une unit d'aire
beaucoup trop grande, le Pascal est une unit beaucoup trop petite. On utilise
gnralement le MgaPascal (1MPa = 106 Pa) ou le GigaPascal (1GPa = 109 Pa).
1 MPa = 106 Pa = 1N / mm
1 GPa = 109 Pa = 1kN / mm
IUTB-LYON1-GMP-DDS
Page 22
TD-2A-S3-F312
1.2. Composantes du vecteur contrainte :

On note :
ih
les composantes de
r
h
si la normale h est quelconque par rapport
aux axes du repre.

z
xh
r
h = yh
dA
x
ij
zh
les composantes de
r
h
r
1) Si la normale h est en concidence avec lun des axes du repre xyz.
La notation est alors la suivante pour les composantes du vecteur contrainte
axe de projection
i j
normale au plan de coupe
xy
r
y / xyz = yy
y
P
r
h .
zy
//
2) Si on utilise comme repre le repre mobile qui contient la normale h
au plan de coupe (le deuxime axe t du repre appartient au plan form par h
et
r
h ).
P
x
IUTB-LYON1-GMP-DDS
th
hh
Plan de Coupe
r
hh
h / ht =
th
Page 23
TD-2A-S3-F312
2. tat de contraintes en un point. Matrice des contraintes :

2.1. Dfinition :
Considrons un paralllpipde infiniment
petit de ct dx, dy, dz isol autour dun
point. Il est en quilibre lastique. Coupons
le paralllpipde par un plan passant par
ses sommets Ax, Ay, et Az. Le ttradre
infiniment petit ainsi obtenu est encore en
quilibre.
z
zdxdy
Az
dz
ydxdz
dy
Ay
dx
La dmonstration de la nature matricielle

des contraintes a t faite par CAUCHY 1
en mettant en quilibre le ttradre
infiniment petit.
Ax
x
xdydz
dAh
Az
-xx
hx
h = hy
hz
dz
-yy
Gx
-yx
-xy
Gy
dAy
-zy
-yz
dx
dAx
-zx
dy
-xz
Ay
Gz
Ax
dAz
IUTB-LYON1-GMP-DDS
-zz
Page 24
TD-2A-S3-F312
n nglige le poids ou les forces dinertie qui sexercent sur le ttradre. En effet
le ttradre tant infiniment petit elles sont dordre suprieur (dxdydz) par
rapport aux forces qui sexercent sur les faces (dxdy).
La rsultante des forces sexerant sur le ttradre est nulle.
r
R = 0
La relation scrit sous forme matricielle :
xx
xh
=
yx
yh
zx
zh
xy
yy
r
x
r
y
r
h
zy
xz hx

yz hy
zz hz
r
z
r
h
[ ] = [][h]
h
Le vecteur h est unitaire (sa norme vaut 1).
La matrice [ ] est la matrice des contraintes au point P. Le vecteur h est le

vecteur contrainte sexerant sur un plan de coupe passant par P et de normale
sortante h (dirige vers lextrieur de la matire).

Le moment rsultant (par rapport G) des forces sexerant sur le ttradre
est nul.
r
MG = 0
Lquation dquilibre du moment donne :
xy = yx
yz = zy
zx = xz
IUTB-LYON1-GMP-DDS
Page 25
TD-2A-S3-F312
Ces trois galits constituent le thorme de rciprocit des contraintes

tangentielles.
Consquence : la matrice des contraintes [ ] est symtrique.
xx xy xz hx
xh
=

=
yy
yz hy
yx
yh
=
hz
zh
zx
zy
zz

Symtrique
[] = []
zz
z
Facette z
yz
xz
zy
zx
1
xy
P
yx
Facette y
y
yy
y
Facette x
xx
x
1
IUTB-LYON1-GMP-DDS
Page 26
TD-2A-S3-F312
2.2. Elments Principaux des .

On recherche sil existe des plans de coupe particuliers passant par P tels que la
normale h au plan de coupe et le vecteur contrainte h soient colinaires.
cas gnral
r
r
h et h non colinaires
cas particuliers
r
r
r
r
h et h colinaires h = h
Le cas particulier o h et h sont colinaires revient chercher en algbre
linaire les vecteurs propres h de lespace, tant les valeurs propres de [ ] .
r
r
h = h
[ ] = [h] = [][h]
h
[[] [I]][h] = [0]

r
Les vecteurs propres h de lespace devant tre non nuls, le systme dquations
admet des solutions diffrentes de 0 que si son dterminant est nul.
xx
xy
xz
yx
yy
yz
= 0
zx
zy
zz
Le dveloppement du dterminant caractristique
me
lquation caractristique du 3
[ ] [I ] = 0
conduit
degr en :
3 I1 2 + I2 I3 = 0
Les termes I1 , I2 , I3 sont des invariants de la matrice des contraintes [ ] (la
recherche des vecteurs propres ne dpend pas de la base initiale xyz).
IUTB-LYON1-GMP-DDS
Page 27
TD-2A-S3-F312
I1 = tr[ ]
est linvariant linaire de
I2 = tr[co ]
[].
est linvariant quadratique de
I3 = []
[].
[] .
est linvariant cubique de
tr [ ] est la trace de la matrice [ ] . La trace dune matrice est gale la somme

des termes situs sur sa diagonale principale.
tr[ ] = xx + yy + zz
La matrice
[] tant symtrique, les 3 valeurs propres sont relles et les 3
vecteurs propres sont orthogonaux et constituent une base dans laquelle la

matrice
[] sera diagonale.
On note X, Y et Z les 3 valeurs propres de [ ] . Ce sont les contraintes

principales. On note X, Y et Z les 3 axes construits sur les 3 vecteurs propres de
lespace. Ce sont les axes principaux des contraintes.
Z
1
Facette Z
Facette Y
Facette X
X
1
[]
XYZ
IUTB-LYON1-GMP-DDS
= 0
0
0
Y
0
0
Z
Page 28
TD-2A-S3-F312
2.3. Equations dquilibre local :

Lorsquon se dplace de dx, dy, dz suivant les axes x, y et z partir dun point P
appartenant un solide, les contraintes varient. De plus le solide peut tre soumis
des forces dinertie variables. Dsignons par f la force par unit de volume
(poids ou forces dinertie).
Un paralllpipde isol autour du point P doit tre en quilibre lintrieur du
corps.
z
zz +d zz
- xx
dx
xz +d xz
- y
- yx
- xy
zx +d zx
G'y
d
y +d y
Gz
- z
Gx
G'x z +d z
-zx
Gy
y +d y
xy +d xy
yx +d yx
- y
G'z
xx +d xx
x
- xz
fdv
dy
-zz
Les quations dquilibre de la somme des forces sexerant sur le paralllpipde

conduisent :
xy
xz
xx
+
+
+ fx = 0
x
y
z
yx
yy
yz
+
+
+ fy = 0
x
y
z
zy
zz
zx
+ fz = 0
+
+
z
x
y
Ce sont les quations dquilibre local.
IUTB-LYON1-GMP-DDS
Page 29
TD-2A-S3-F312
2.4. Expression des composantes de h sur un repre ht li llment

daire :
y
t
h
h
th
hh
hh =
xx + yy
th =
xx yy
2
xx yy
cos 2 + xy sin 2
sin 2 + xy cos 2
Ce sont des formules identiques celles qui existent pour les inerties dune
surface plane.
2.5. Elments principaux des contraintes :

a) A la surface dun corps quelconque, dans une membrane ou dans
une plaque :
xx
[]xyz = yz
0
z Z
xy
yy
0
0
0
Contrainte principale
Etant donne la forme plane de la matrice on peut dire que 0 est une contrainte
principale et que laxe z est un ax principal ( z Z ).
Le calcul se limite donc la recherche des deux contraintes principales ( X et Y )
et des deux axes principaux (X et Y) dans le plan xy.
IUTB-LYON1-GMP-DDS
Page 30
TD-2A-S3-F312
X =
Y =
xx + yy
1
( xx yy ) 2 + 4 2xy
2
1
( xx yy ) 2 + 4 2xy
2
+
2
xx + yy
Z = 0
X
Y
tg =
= (x, X)
X xx
xy
Le troisime vecteur propre Y est perpendiculaire X (de +

contraintes.
b) Dans une poutre :
[ ]xyz
=
0
) dans le plan des

2
0 0
0 0
1
2 + 4 2
2
2
1
Y =
2 + 4 2
2
2
Z = 0
X
tg =
= (x, X )
X =
IUTB-LYON1-GMP-DDS
Page 31
TD-2A-S3-F312
Critres de Rsistance:
1 Critres de dpassement de la limite lastique pour des
matriaux ductiles :
1.1 Limites lastiques dun matriau homogne et isotrope :
Traction
Compression
Torsion
Seule la limite lastique en traction est connue et impose par la norme.
Les limites lastiques sont considres comme des nombres positifs.

On classe les matriaux en deux catgories : les ductiles et les fragiles.
Matriaux ductiles
e e/2
ee
e>e
e e
Matriaux fragiles
1.2 Critres lmentaires (pour un matriau ductile) :

a) Traction :
F
X
[] = 0
XYZ
0
0 0
0 0
0 0
Avec : X =
F
A
Le coefficient de scurit se calcule par : =
IUTB-LYON1-GMP-DDS
e
X
Page 32
TD-2A-S3-F312
b) Compression :
F
X
[] = 0
XYZ
0
0 0
0 0
0 0
Avec x =
F
A
Le coefficient de scurit se calcule par : =
e
X
c) Torsion :
0 0
[ ] = 0 0
xyz
0 0 0
La contrainte maximum de torsion se calcule par :
=
z
C
R
IG
On montre que les valeurs propres en torsion

sont opposes :
0
0
X = +
[] = 0
Y = 0
XYZ
0
0
0
Le coefficient de scurit se calcule par (puisquon a suppos que le matriau est

ductile et que par consquent e =
e
):
2
e = e 2
= ( X Y ) 2
IUTB-LYON1-GMP-DDS
e
( X Y )
Page 33
TD-2A-S3-F312
1.3 Critres de VON MISES 2 et de TRESCA 3 :

Ce sont deux critres de dpassement de ltat lastique du corps, pour des
matriaux ductiles. Ils sont bass tous les deux sur le comportement en traction
du matriau. Ils sappliquent des tats de contraintes quelconques (3D ou 2D). Ils
sont trs voisins lun de lautre (lcart maximum entre les deux critres est
denviron 10%). Le critre le plus rcent et le plus utilis est celui de VON MISES.
On recherche une contrainte dite quivalente au sens de la traction telle que le
coefficient de scurit puisse se mettre sous la forme :
e
equ
La contrainte equ est une fonction f dterminer des contraintes principales X, Y,

et Z.
Par exemple, pour les critres lmentaires :
En traction : f(X,Y,Z) = X
En compression : f(X,Y,Z) = |X|

En torsion : f(X,Y,Z) = X - Y
1.4 Critres de VON MISES ou critre de la contrainte tangentielle

octadrique:
1) Contrainte tangentielle octadrique :
Z
oct
h=1/3
1
1
Y
X
2
3
VON MISES Richard (1883-1953)

TRESCA Henri (1814-1885)
IUTB-LYON1-GMP-DDS
Page 34
1
1
1
TD-2A-S3-F312
Considrons loctadre rgulier (les huit faces sont des triangles quilatraux)
isol autour dun point P dans les axes principaux. Calculons le vecteur contrainte
r
h
sexerant sur la face ombre de loctadre.
[ ] = [ ]
oct
= [ ][h] = 0
0
0
Y
0
0
1
X
1
1

0
Y
1 =
3
3
1
Z
Z
Calculons la contrainte normale octadrique, puis la contrainte tangentielle

octadrique.
r
r
+ Y + Z
hh = oct = h h = X
3
2oct 2oct =
th = oct =
oct =
1
3
( X + Y + Z ) 2
2X + 2Y + 2Z
3
9
Y )2 + ( Y Z )2 + ( Z X )2
2) Enonc du critre :
Pour quun solide reste dans ltat lastique il faut que la contrainte tangentielle
octadrique soit infrieure la limite quelle atteint dans un essai de traction.
oct lim
oct
3) Calcul de
lim
oct
La matrice des contraintes en traction scrit lorsquon atteint la limite lastique :
e
[ ] = 0
XYZ
0
lim
oct =
1
3
Y )2 + ( Y Z )2 + ( Z X )2 =
lim
oct =
IUTB-LYON1-GMP-DDS
0 0
0 0
0 0
1
2 2e
3
2
e
3
Page 35
TD-2A-S3-F312
4) Formulation du critre en 3D. Contrainte quivalente de VON

MISES :
En introduisant un coefficient 1 (le coefficient de scurit) :
oct = lim
oct
1
3
2
e
3
Y )2 + ( Y Z )2 + ( Z X )2 =
( 1
2 )2 + ( 2 3 )2 + ( 3 1 )2
est le coefficient de scurit.
VM
equ =
1
2
( 1
2 )2 + ( 2 3 )2 + ( 3 1 )2
5) Formulation et image du critre en 2D :
Pour un corps quelconque :

La formule en 3D de VON MISES scrit en 2D compte tenu que Z = 0 :
X
[ ] = 0
XYZ
0
0
Y
0
0
2X + 2Y X Y
VM
equ =
2X + 2Y X Y
La courbe limite (lorsque le coefficient de scurit vaut 1) est une ellipse

dquation : 2X + 2Y X Y = 2e
Y
Z=0
Ellipse de VON MISES
M
TRACTION
TRACTION
45
COMPRESSION
COMPRESSION
Tr
T
T
TORSION
Hexagone de TRESCA
IUTB-LYON1-GMP-DDS
Page 36
TD-2A-S3-F312
Ltat du corps sera lastique si le point M reprsentatif dun tat de contraintes

planes se trouve lintrieur de lellipse. Ltat limite est reprsent par le point
M situ sur lellipse.
Le coefficient de scurit peut alors de mettre sous la forme : =
Pour une poutre : [ ]xyz
=
0
0 0
0 0
OM'
OM
On reporte les expressions des contraintes principales pour une poutre dans
2
X
+ 2Y X Y
. On rappelle que les contraintes principales pour une
poutre valent (en fonction de la contrainte de flexion et de la contrainte de

torsion ) :
YX =
+ 4 2
On obtient ainsi lexpression du coefficient de scurit de VON MISES pour une

poutre qui est la fois flchie et tordue :
2 + 3 2
VM
equ =
2 + 3 2
6) Vrification du critre la compression et la torsion :

Le critre de VON MISES tant bas sur la traction il faut quil vrifie la
compression et la torsion.
Compression
Daprs lellipse de VON MISES, on constate que les limites lastiques en traction
(point Tr) et en compression (point C) doivent tre identiques. Autrement dit le
matriau doit tre symtrique, ce qui est caractristique dun matriau ductile.
'e = e
Torsion
En torsion les valeurs propre sont opposes (point T sur lellipse). La relation de
scurit pour les poutres donne la limite puisque =1 et = e .
1=
0 + 3 2e
2
On obtient donc la relation entre les limites lastiques en traction et en torsion :
e =
e
3
0.577 e
Cette relation caractrise un matriau ductile.
IUTB-LYON1-GMP-DDS
Page 37
TD-2A-S3-F312
N.B. Le critre de VON MISES est donc un critre bien adapt aux
matriaux ductiles, mais on peut quand mme lappliquer aux matriaux
fragiles. On obtient alors une scurit un peu plus faible quelle nest en
ralit.
1.5 Critres de TRESCA ou critre de la contrainte tangentielle maximum :

1) Enonc du critre :
Pour quun solide reste dans ltat lastique il faut que la contrainte tangentielle
maximum soit infrieure la limite quelle atteint dans un essai de traction.
lim
Maxi
Maxi
lim
2) Calcul de Maxi :
th
lim
Maxi =
Cercle de MOHR
en traction
2
hh
X=e la limite de traction
lim
Maxi =
e
2
3) Formulation du critre en 3D. Contrainte quivalente de TRESCA :
Maxi
e
2
Ou encore, en introduisant un coefficient 1 (le coefficient de scurit) :
Maxi =
e
2
Le cisaillement maximum est le rayon du plus grand cercle de MOHR :
Maxi = Max
Maxi = Max
IUTB-LYON1-GMP-DDS
i j
2
i j
2
e
2
Page 38
TD-2A-S3-F312
e
e
=
2 Maxi
Max i j
est le coefficient de scurit.

i, j sont les contraintes principales
T
equ
= 2 Maxi
4) Formulation et image du critre en 2D :
Pour un corps quelconque :

[ ] = 0 Y 0
XYZ
On montre que la courbe limite est un
hexagone
0
0
0
inscrit dans lellipse de VON MISES.
Ltat du corps sera lastique si le point M reprsentatif dun tat de contraintes
planes se trouve lintrieur de lhexagone. Ltat limite est reprsent par le
point M situ sur lhexagone.
OM''
Le coefficient de scurit peut alors de mettre sous la forme : =
OM
0
Pour une poutre : [ ]xyz = 0 0

0 0 0
Y
Z
2
2
1
=
2
2
= 0
X =
+ 4 2 > 0
+ 4 2 < 0
Etant donn le signe des contraintes principales, le cisaillement maximum est
Z =
X Y
1
=
2 + 4 2
2
2
La relation de scurit scrit :
IUTB-LYON1-GMP-DDS
e
=
2 Maxi
2 + 4 2 T =
equ
e
+ 4 2
2
2 + 4 2
Page 39
TD-2A-S3-F312
5) Vrification du critre la compression et la torsion :

Compression
Daprs lhexagone de TRESCA, on constate que les limites lastiques en traction

(point Tr) et en compression (point C) doivent tre identiques. Autrement dit le
matriau doit tre symtrique, ce qui est caractristique dun matriau ductile.
'e = e
Torsion
En torsion les valeurs propre sont opposes (point T sur lhexagone de TRESCA).
La relation de scurit pour les poutres donne la limite puisque =1 et = e .
1=
0 + 4 2e
2
On obtient donc la relation entre les limites lastiques en traction et en torsion :
e =
e
2
Cette relation caractrise un matriau ductile.
IUTB-LYON1-GMP-DDS
Page 40
TD-2A-S3-F312
2. DEFORMATIONS
PROBLME N9
Considrons le carre infiniment petit, de ct unit, dans les
axes xy.
Les dformations dans le plan x,y valent :
xx = yy = 50d
xy = -100 d
1
x
1) Tracer la figure dforme du carr.

2) Dterminer les lments principaux de la dformation.
Tracer la figure dform dun carr isol dans les axes
principaux et de ct
RPONSES N9
y
Matrice de dformation :
[ ]
P
xyz
xx
= yx
zx
xy xz
yy yz
zy zz
Dans notre cas :
[ ]
P
xyz
100.10-6
50.10-6
xy
1
100 0
50
0 d
= 100 50
0
zz
0
l
l x
xx =
l = 1 xx
or l=1
50.10-6
xy
100.10-6
l y
De mme pour yy =
Le carr se transforme aprs dformation en un paralllogramme.

Langle initialement droit entre les axes x et y vaut aprs dformation
IUTB-LYON1-GMP-DDS
xy =
2 xy
Page 41
TD-2A-S3-F312
Elments principaux de la dformation :
[ ]
P
xyz
100 0
50
0 d
= 100 50
0
zz
0
Y
La direction z est principale et Z= zz.
2 .50.10-6
Les 2 autres valeurs propres sont donnes par :
=
Y
X
xx + yy
2
1
m
2
X=150d
yy ) + 4
2
xx
2
2
xy
P
2
Y=-50 d
45
La direction propre X est calcule par :
xx
tan = X
= 1
xy
r r
dou = 45 = x , X
2 .150.10-6
X
[ ]
P
xyz
0
0
150
= 0 50 0 d
0
0 Z
Le carr se transforme aprs dformation, dans les axes principaux, en un rectangle. Langle
droit reste droit aprs dformation, car la distorsion angulaire est nulle. Les dformations
principales X et Y reprsentent les extremums des dilatations linaires unitaires au point P.
IUTB-LYON1-GMP-DDS
Page 42
TD-2A-S3-F312
PROBLME N10
F
Considrons une barre en traction.

1) Ecrire la matrice des dformations dans le
repre principal en fonction de la dformation longitudinale
X et du coefficient de POISSON . du matriau.
2) Dterminer dans le repre principal XY la
position angulaire des axes h pour lesquelles hh=0.
Application numrique : Le matriau est un AU4G. Calculer
la position angulaire en des axes h par rapport laxe de
traction.
RPONSES FN10
Matrice de dformation dans le repre principal:
[ ]
P
XYZ
X
= 0
0
0
0
Z
Y
0
Laxe de traction X est principal et enregistre

le maximum des longations unitaires : X
Laxe Y est principal et enregistre
le maximum des raccourcissements unitaires : Y
Daps la loi de POISSON : Y = X et Z = X
[ ]
P
XYZ
X
= 0
0
IUTB-LYON1-GMP-DDS
0
X
0
0
0
X
Page 43
TD-2A-S3-F312
[ ]
P
XYZ
0
0

1 0
= X 0
0 0
Ecrivons lexpression donnant hh dans des axes quelconques :
hh =
xx + yy
xx yy
2
Y
P
cos 2 + xy sin 2
Adaptons cette relation dans les axes principaux :
hh =
X + Y
2
X Y
2
cos 2
car la distorsion angulaire est nulle .
Nous cherchons la position de la direction h qui enregistre une dilatation linaire unitaire nulle
. Si une jauge tait colle dans la direction h, elle mesurerait zro, quelque soit la valeur de
leffort de traction.
hh =
X + Y
2
cos 2 =
cos 2 =
X Y
X + Y
X Y
cos 2 = 0 ;
cos 2 =
1
X X
=
X + X
1+
or cos 2 = 1 tan
1
1 + tan 2
1+
tan =
do
tan 2 =
Le coefficient de POISSON dun AU4G est gal 0,33
tan =
IUTB-LYON1-GMP-DDS
1
0,33
60
Page 44
TD-2A-S3-F312
PROBLME N11
y
On considre un point P la surface d'un corps
en un endroit o ne s'applique aucune force
extrieure. Les rsultats enregistrs sur chaque
J2
jauge d'une rosette 45 colle dans le plan
J
45
tangent en P sont respectivement :
jauge J1: 950 d
J1
x
P
jauge J:-175 d
jauge J2:-475dd
1) Quels sont les lments principaux de la dformation.
2) Dans quelle direction ' enregistrerait-on une dilatation linaire nulle?
3) Dans quel systme d'axes enregistrerait-on une distorsion extremum ?
RPONSES N11
Le point P se situant la surface dun corps o ne sapplique aucune force extrieure, la forme
de la matrice des dformations est la suivante :
y
[ ]
P
xyz
xx
= yx
0
xy
yy
0
0
0
zz
Les jauges Jx et Jy mesurent respectivement

Les dilatations linaires unitaires xx et yy .
La demi distorsion angulaire xy est dterminer
indirectement par la troisime jauge J .
Jy
45
J
Jx
P
La dilatation linaire est calcul par les lois
de comportement et ncessite la connaissance du coefficient de Poisson
[ ]
P
xyz
IUTB-LYON1-GMP-DDS
0
950 xy
= yx 475 0 d
0
zz
0
Page 45
TD-2A-S3-F312
Pour calculer xy , utilisons la relation hh =
xx + yy
2
xx yy
2
cos 2 + xy sin 2 en faisant
coincider h avec .
=
xy =
xx + yy
2
xx + yy
xx yy
2
cos 90 + xy sin 90
xy = 175
[ ]
P
xyz
950 475
2
xx + yy
2
+ xy
xy=-412,5.10-6
412,5 0
950
0 d
= 412,5 475
0
zz
0
Elments principaux de la dformation :

La direction z perpendiculaire au plan o est coller la rosette est principale.
Nous recherchons donc les directions propres, et les valeurs propres dans le plan de la
rosette xy.
Les 2 autres valeurs propres sont donnes par :
Y =
xx + yy
1
2
X=1061 d
yy ) + 4 xy2
2
xx
Y =
X
950 475 1
m
2
2
(950 + 475)2 + 4 412,52
Y=-586 d
La direction propre X est calcule par :
tan =
X xx 1061 950
=
0,269
xy
412,5
r r
dou = 15 = x , X
[ ]
P
XYZ
0
0
1061
= 0
586 0 d
0
Z
0
IUTB-LYON1-GMP-DDS
Jy
45
J
Jx
x
15
X
Page 46
TD-2A-S3-F312
Cherchons la direction pour laquelle on enregistrerait une dilatation linaire nulle :
hh =
X + Y
2
X Y
2
cos 2 = 0 ; cos 2 =
X + Y
X Y
cos 2 =
1061 586
0,288
1061 + 586
534
Cherchons le systme daxes pour lequel on enregistrerait une distorsion extremum :
( th )max i = (2 th )max i
th =
xx yy
2
sin 2 + xy cos 2
Dans les axes principaux,

la relation devient :
th =
X Y
2
sin 2
Jy
=45
J
534
ht est maximum pour sin2 maxi

sin2=1
45
Jx
x
15
X
534
IUTB-LYON1-GMP-DDS
Page 47
TD-2A-S3-F312
DFORMATIONS:
1. Dplacements et Dformations
1.1. Vecteurs dplacement dun point P :
CHARGEMENT
Trajectoire de P
u(P
y
S
Soit un point P quelconque du solide dans son tat initial S. Aprs chargement le
solide passe dans son tat final S et le point P dcrit une certaine trajectoire qui
lamne dans la position P.
Soit u(P ) = P P' le vecteur dplacement du point P. Cest une fonction vectorielle
des trois coordonnes initiales x, y, et z du point P.
Pour simplifier au maximum on suppose que le dplacement u(P ) = P P':

Est petit (1% de la plus grande dimension du solide). Cest lhypothse
des petits dplacements. De ce fait on peut considrer que le dplacement
de P (qui seffectue thoriquement le long de la trajectoire) est tangent
cette trajectoire.
Quil possde des drives partielles par rapport au variables x, y et z

faibles (10-3). Les dformations dans le solide qui dpendent des drives
seront donc du mme ordre. Cest lhypothse des petites dformations.
IUTB-LYON1-GMP-DDS
Page 48
TD-2A-S3-F312
1.2. Dcomposition d'une matrice carre quelconque :
z
Trajectoire de P
P
P
u(P)
r
k
v
i
u
x

OP = y u( P ) = v
w
z
r
j
O
Une matrice carre quelconque [T] peut toujours tre dcompose en une somme
d'une matrice symtrique []et d'une matrice antisymtrique.[r]
[T ] = [ ] + [r ]
S
AS
[T ]
carre quelconque de terme Tij .
[ ]
carre symtrique de terme ij ( ij = ji ).
[r ]
carre
antisymtrique
de
terme
rij
( rij = rji
et
par
consquent rii = 0 )
Avec :
1) [ ] =
1
([T ]+ t [T ])
2
La matrice [ ] est donc symtrique.

2) [r ] =
1
([T ] t [T ])
2
La matrice [r ] est donc antisymtrique.
IUTB-LYON1-GMP-DDS
1
( T + Tji )
2 ij
1
= ( Tji + Tij ) = ij
2
ij =
ji
1
( T Tji )
2 ij
1
rji = ( Tji Tij ) = rij
2
rij =
Page 49
TD-2A-S3-F312
1.3. Petite rotation autour d'un axe quelconque :

Considrons un solide soumis une petite rotation autour dun axe pendant un
intervalle de temps dt.
dx
PQ = dy
d
P
dz
Q
du
du
du = dv
dw
r
r
V( Q ) = PQ
r
r
du = V( Q )dt = dt PQ = d PQ
du
rx
dx
ry dz rz dy
du = dv = ry dy = rz dx rx dz
dw
rz
dz
rx dy ry dx
Soit sous forme matricielle:
0
du

dv
=
rz

ry
dw
rz
0
rx
ry dx

rx dy
0 dz
[du] = [r ][da ]
antisymtrique
Une petite rotation est donc caractrise par une matrice antisymtrique [r ] .
IUTB-LYON1-GMP-DDS
Page 50
TD-2A-S3-F312
1.4. Etat de dformations en un point. Matrice des

dformations :
1.4.1. Relations dplacements-dformations :
Pour valuer les dformations dans le solide, considrons deux points voisins P et
r
r
Q dans le solide. Soient u ( P ) et u ( Q ) les dplacements (tangents aux
trajectoires) des deux points P et Q. Lcart entre u ( Q ) et u ( P ) est la
diffrentielle du du vecteur u ( P ) .
Trajectoire de Q
u(Q)
du
Q
z
dur
dx
TRANSLATION
ROTATION
DEFORMATION
//
dx
x
du
u(P)
Trajectoire de P
La petite transformation de dx en dx' est la somme de trois transformations

lmentaires et petites (lordre dans lequel seffectue les transformations
nintervient pas) :
Une translation caractrise par u ( P0 ) .

Une rotation autour de P caractrise par dur .
Une dformation caractrise par du .
IUTB-LYON1-GMP-DDS
Page 51
TD-2A-S3-F312
La rotation et la dformation sont contenues dans du .
du = dur + du
r
La diffrentielle du dun vecteur u ( P ) se calcule par :
r
u u r v r w r
=
k
j+
i +
x x
x
x
r
u u r v r w r
=
k
j+
i +
y
y
y
y
r
u u r v r w r
k
=
j+
i +
z
z
z
z
r
r
r
r
u
u
u
du =
dx +
dy +
dz Avec :
x
y
z
On reporte les drives partielles de u dans la diffrentielle du .
r
r
u
u
r
u
v
v
v
du = (
dx +
dy +
dz)i + (
dx +
dy +
dz) j +
x
y
z
x
y
y
r
w
w
w
(
dx +
dy +
dz)k
x
y
z
u
u
u
du =
dx +
dy +
dz
x
y
z
r
r
r
r
v
v
v
dx +
dy +
dz
Or : du = du i + dvj + dwk en identifiant : dv =
a1
a2
a 3
dw =
w
w
w
dz +
dy +
dz
x
y
z
Soit sous forme matricielle :
du
x

v
dv = x
dw
w
u
y
v
y
w
y
z dx
v
dy
z

w dz
[du] = [T][dx]
La matrice carre (non symtrique) [T] est la matrice de transformation.
x
[T ] = v
x
w
IUTB-LYON1-GMP-DDS
u
y
v
y
w
y
z
v
z
w
Page 52
TD-2A-S3-F312
La matrice [T] est toujours dcomposable en une somme d'une matrice

symtrique [] et d'une matrice antisymtrique [r].
La matrice symtrique [] caractrise la dformation. La matrice
antisymtrique [r] caractrise la rotation.
[du] = [T ][dx ]
[du] = [[] + [r ]][dx ] avec : [] = 1 ([T ]+ t [T ]) et [r ] = 1 ([T ] t [T ])

2
2
[du] = [][dx ] + [r ][dx ] do [du] = [du ] + [dur ] avec [du ] = [][dx ] et [dur ] = [r ][dx ]
Le vecteur du reprsente la diffrentielle pendant la dformation. Le vecteur du r

reprsente la diffrentielle pendant la rotation.
La matrice carre symtrique [ ] = 1 ([T ]+ t [T ]) caractrise donc la dformation

2
pure au voisinage du point P. Cest la matrice des petites dformations de

GREEN 4 au point P. Ses termes sont petits (de lordre de 10-3).
xx xy xz
=
=
[ ] = yx yy yz
xyz
zx zy = zz
avec :
[du ] = [][dx ]
De la relation [ ] =
dformations :
[]= []
Symtrique
1
([T ]+ t [T ]) on tire les relations dplacements2
xx =
xy = yx =
1 u v
(
+
)
2 y x
yy
yz = zy =
1 v w
(
+
)
2 z
y
zx = xz =
1 w u
(
+
)
2 x
z
zz
u
x
v
=
y
w
=
z
GREEN George (1793-1841)
IUTB-LYON1-GMP-DDS
Page 53
TD-2A-S3-F312
1.4.2. Signification gomtrique des termes de [] :

a) Dforme dun cube de ct unit :
Considrons un cube de ct unit isol en P dans le repre xyz avant dformation
et cherchons la dforme de ce cube.
Les vecteurs diffrentiels des cts PAx, PAy et PAz pendant la dformation se
[du ] = [][da]
calculent avec :
Ct PAx du cube :
Ct PAy du cube :
Ct PAz du cube :
[du ]x
[du ]y
[du ]z
xx
= yx
zx
xy
yy
xx
= yx
zx
xy
yy
xx
= yx
zx
xy
yy
zy
zy
zy
xz 1
xx

yz 0 = yx
zx
zz 0
xy
xz 0

yz 1 = yy
zy
zz 0

xz 0
xz

yz 0 = yz
zz
zz 1
z
r
( du )x = AxAx'
r
( du )y = AyAy'
r
( du )z = AzAz'
xz
yz A'z zz
Az
P
Ay
A'x
xx
x
IUTB-LYON1-GMP-DDS
yx
Ax
yy
A'y
zy xy
zx
Page 54
TD-2A-S3-F312
Dans un repre quelconque un cube infiniment petit de ct unit

(PAx=PAy=PAz=1) se transforme en un paralllpipde faces obliques.
b) Interprtation des termes xx , yy , zz , :
Calculons le rapport (
l
)x .
l
l
)x est la dilatation linaire unitaire (sans dimension) suivant laxe x du
l
repre.
Les termes sans dimension xx , yy , zz , reprsentent donc sensiblement les
variations de longueurs des cts du cube, quand les cts de celui-ci valent 1.
Si les cts du cube, ne valent pas 1 mais l, on a les relations :
xx (
l
)x
l
yy (
l
)y
l
zz (
l
)z
l
Remarque :
Le rapport
V V0
V0
est la dilatation cubique unitaire (ou dilatation volumique) en
un point P.
Volume initial du cube : V0=1
Volume final V du paralllpipde faces obliques :
V = ( PAx' PAy' ) . PAz'
1 + xx
xy
zx
PAx' PAy' =
yx 1 + yy
zy
zx
zy
1 + xx + yy
V = ( PAx' PA'y ) . PAz' 1 + xx + yy + zz 1 + tr[ ]
V V0
tr[ ]
V0
La dilatation cubique est une constante physique qui ne dpend pas du repre o
elle est calcule. La trace tr [] de la matrice des dformations est donc
invariante.
IUTB-LYON1-GMP-DDS
Page 55
TD-2A-S3-F312
c) Interprtation des termes xy , yz , zx :
Calculons langle Ax' P Ay' entre les vecteurs PAx' et PAy' en posant :
Ax' P A'y =
xy
2
xy est la distorsion de langle droit x,y du repre (en rd).
xy 2 xy
yz 2 yz
zx 2 zx
Les termes sans dimension xy , yz , zx , reprsentent donc sensiblement les demi-
distorsions des faces du cube.
z
(1+xx)
1
(1+zz)
(/2-yz)
(/2-zx)
(/2-xy)
y
(1+yy)
REMARQUE : On peut, puisque la matrice des dformations [] est

symtrique, et quelle ne comporte donc que six termes indpendants, utiliser
pour la matrice des dformations [] une notation, assez pratique, en colonne :
[]
xyz
IUTB-LYON1-GMP-DDS
xx

yy
zz
=
xy
yz

zx
Attention !
ij = 2 ij
Page 56
TD-2A-S3-F312
1.4.3. Elments principaux de [] :

On recherche sil existe des axes h particuliers passant par P tels que le vecteur
unitaire h de laxe h et la diffrentielle de la dformation du de h soient

colinaires.
Le schma suivant ne reprsente que la dformation pure dun vecteur qui est le
vecteur unitaire h dun axe h passant par P (la rotation et la translation ne sont
pas reprsentes).
(dform)
du
(dform)
du
P
(norme 1)
y
cas gnral
r
du et h non colinaires
(norme 1)
y
cas particuliers :
r
du = h
Comme pour les contraintes, le problme revient chercher les vecteurs propres
r
h de lespace, tant les valeurs propres de [ ] .
r
Les vecteurs propres h de lespace devant tre non nuls, le systme dquations
admet des solutions diffrentes de 0 que si son dterminant est nul.
xx
xy
xz
yx
yy
yz
= 0
zx
zy
zz
Le dveloppement du dterminant caractristique de [ ] [I ] conduit

lquation caractristique :
3 I1 2 + I2 I3 = 0
IUTB-LYON1-GMP-DDS
Page 57
TD-2A-S3-F312
I1 = tr[] est linvariant linaire (la dilatation cubique).

I2 = tr[co ] est linvariant quadratique.
I3 = [ ]
est linvariant cubique.
La matrice [] tant symtrique, les 3 valeurs propres sont relles et les 3

vecteurs propres sont orthogonaux et constituent une base dans laquelle la
matrice [] sera diagonale.
On note, X, Y et Z les 3 valeurs propres de [ ] . Ce sont les dformations
principales. On note X, Y et Z les 3 axes construits sur les 3 vecteurs propres de
lespace. Ce sont les axes principaux des dformations
Z
1
P
1
X
Dans le repre principal un cube infiniment petit de ct unit se transforme
en un paralllpipde faces rectangulaires.
X
[ ] = 0
XYZ
0
0
Y
0
0
Z
NOTATION EN COLONNE POUR [] :
IUTB-LYON1-GMP-DDS
[]
XYZ
X

= Y
Z
Page 58
TD-2A-S3-F312
1.5. Etat de dformations planes :

Tout ce qui a t crit propos de ltat de contraintes planes est
intgralement transposable ce paragraphe consacr ltat de dformations
planes. Aucune dmonstration ne sera donc faite dans ce paragraphe. Il suffit
de faire le parallle entre le vecteur contrainte
r
h et le vecteur diffrentiel
r
du quand la norme du vecteur qui se dforme dx vaut 1 ( dx = h ).
h
h
(dform)
du
hh
th
h
hh
th
(norme 1)
[du ] = [][h]
(norme 1)
[ ] = [][h]
[ ] =
hh

th
[du ] =
hh

th
1.5.1. Dfinition :
Un tat de dformations planes est tel que toutes les dformations sont contenues
dans un plan (xy par exemple).
La matrice des dformations scrit :
[]xyz
IUTB-LYON1-GMP-DDS
xx
= yx
0
xy
yy
0
0
0
[]
xyz
xx

= yy
xy

Page 59
TD-2A-S3-F312
La visualisation dun tat de dformations planes xy peut se faire laide dun carr
de ct unit isol dans le plan xy.
xy
yy
1
(/2-2xy)
(1+yy)
(1+xx)
yx
P
xx
NB. Si zz 0 ltat de dformation est dit pseudo plan mais cet tat est
rgit par les mmes quations que ltat plan.
1.5.2. Expression des composantes de du sur un repre ht li au point P :

y
t
(dform)
xx + yy
th =
IUTB-LYON1-GMP-DDS
hh
th
hh =
du
xx yy
2
(norme 1)
xx yy
2
cos 2 + xy sin 2
sin 2 + xy cos 2
Page 60
TD-2A-S3-F312
1.5.3. Elments principaux des dformations :
xx
= yz
0
[]xyz
z Z
xy
yy
0
0
0
Dformation principale
Etant donne la forme plane de la matrice on peut dire que 0 est une dformation
principale et que laxe z est un axe principal ( z Z ).
Le calcul se limite donc la recherche des deux dformations principales ( X et
Y ) et des deux axes principaux (X et Y) dans le plan xy.

Le calcul donne :
xx + yy
X =
Y =
xx
Z = 0
tg =
2
+ yy
2
X xx
xy
1
( xx yy ) 2 + 4 2xy
2
1
( xx yy ) 2 + 4 2xy
2
+
= (1, X )
X
Y
Le deuxime vecteur propre Y est perpendiculaire X (de +

dformations.
IUTB-LYON1-GMP-DDS
) dans le plan des

2
Page 61
TD-2A-S3-F312
EXTENSOMTRIE
1. But :
z
Dtermination exprimentale d'un

tat plan de contrainte en un point P.
P
x
2. Principe :
[]
Diagonalisation
/xyz
[]
Lois de Lam ( E et )
/XYZ
[]P/XYZ
Mesure
3. Mthode :
Grille
Support isolant
Protection
Pour mesurer
l'tat de
dformation au
point P, on
utilise des
"jauges de
dformation".
Ce sont des
rsistances
(R=120) trs
minces que l'on
colle sur la
structure au
point P.
Cosse de liaison
IUTB-LYON1-GMP-DDS
Page 62
TD-2A-S3-F312
3.1 Une jauge de dformation est constitue d'une grille trs fine (2m)
imprim selon une technique de photogravure sur un support isolant trs mince
(30m). Le fil qui constitue la grille est gnralement en constantan (Cu+Ni) qui est
un alliage dont la rsistance est stable aux alentours de 20C. Le support isolant
est un polyimide.
3.2 Facteur de jauge :

Il est dtermin en traction.
Les brins de la grille tant aligns suivant x, le fil subit une dformation de
traction (comme s'il tait tendu ses extrmits). Le fil subit donc une
dformation xx suivant x et une dformation -xx suivant y cause de l'effet de
POISSON.
La rsistance du fil vaut :
R =
L
S
R
L S
=
+
L
S
S
D
L
= 2
= 2
L
S
D
En outre, on montre que lors de la dformation la rsistivit du fil varie de faon

proportionnelle son volume (loi de BRIDGMAN). C tant la constante de
BRIDGMAN du matriau constituant le fil, on peut crire :
V
L S
L
= C
= C(
+
) = C (1 2 )
V
L
S
L
R
L
= ( C (1 2 ) + (1 + 2 ))
R
L
R
L
= K
R
L
R
= K xx
R
avec
K = C (1 2 ) + (1 + 2 )
La constante K est le facteur de jauge. Pour le constantan C1.13 et 0.285

K2.05. La valeur de K donne par le constructeur est dtermine
exprimentalement en flexion 4 points par prlvements statistiques. En flexion 4
points, la partie suprieure de la barre est soumise de la traction.
IUTB-LYON1-GMP-DDS
jauge
Page 63
TD-2A-S3-F312
3.3 Les instruments de mesure :

Pour mesurer la variation de rsistance R/R on utilise un pont de WHEATSTONE.
On montre que :
R1
R2
e = E(
e ( xx)
Le pont est dit en quilibre si e=0 R1R3=R2R4

Supposons que les quatre rsistances varient.
On obtient :
R3
R4
R4
R1
)
R1 + R2 R3 + R4
e = E
R1R2
(R1 + R2 )
R1
R1
R2
R2
R3
R3
R 4
R4
On remarque que deux rsistances adjacentes agissent en sens oppos, tandis que
deux rsistances opposes agissent dans le mme sens.
3.4 Les cblages :

Le cblage le plus simple est le "montage quart de pont - deux fils". On remplace
une rsistance du pont par une jauge J, de rsistance R. Aprs avoir ralis
l'quilibre du pont (e=0) et dform la structure sur laquelle est colle la jauge, la
rsistance de la jauge varie de R. La tension de sortie e est donne par :
J
x (xx)
R2
e ( xx)
E R
( )
4 R
E
(K xx )
4
4
xx =
e
KE
La valeur de xx se lit directement sur le
cadran du voltmtre en
microdformations (1d=10-6 m/m).
e =
R3
R4
e =
L2
J
x (xx)
R2
L1
e ( xx)
R3
R4
IUTB-LYON1-GMP-DDS
La dilatation des fils de liaison, si la

temprature varie pendent le mesure,
peut entraner une erreur non
ngligeable. On utilise alors le "montage
quart de pont - trois fils". Si L1=L2 les
variations de rsistance des fils dues
des effets thermiques se compensent
puisqu'ils sont situs dans deux branches
adjacentes du pont.
Page 64
TD-2A-S3-F312
3.5 Les effets parasites :

Le principal est la sensibilit latrale de la jauge. La rsistance d'une jauge ne se
modifie pas uniquement en raison de la dformation longitudinale xx suivant laxe
de ses brins mais aussi, pour une faible part, en raison de la dformation
transversale yy. (en effet la jauge ne travaille pas, lorsquelle est colle sur une
structure quelconque, dans les mmes conditions que lors de son talonnage en
traction). On peut effectuer une correction de la dformation xx en utilisant le
facteur de sensibilit transversale Kt de la jauge (ce facteur est fourni par le
constructeur).
4. Mesure de l'tat de dformation la surface d'un corps :

y
xy
yy
1
(1+xx)
(/2-2xy)
(1+yy)
yx
P
xx
A la surface d'un corps la matrice des dformations s'crit :
[ ]xyz
xx
= yx
0
xy
yy
0
0
zz
La dformation zz n'est pas indpendante. Elle vaut :
zz =
( xx + yy )
1
Au voisinage d'un point, et la surface d'un corps, un carr infiniment petit de

ct unit se transforme en un paralllogramme.
IUTB-LYON1-GMP-DDS
Page 65
TD-2A-S3-F312
Il faut donc mesurer trois quantits pour dfinir compltement la dformation :
xx qui est la dilatation linaire unitaire suivant l'axe x.
yy qui est la dilatation linaire unitaire suivant l'axe y.

xy qui est la demi-distorsion xy de l'angle droit xy.
xx et yy peuvent se mesurer avec deux jauges places suivant les axes x et y. La
demi-distorsion xy tant un angle, on ne peut pas la mesurer directement avec
une jauge. On effectue une troisime mesure de dilatation linaire unitaire dans le
plan xy suivant un axe h.
y
t
(dform)
du
hh
th
P
(norme 1)
On a montr au 3.5.3. que la dilatation linaire unitaire hh suivant laxe h se

calcule par :
hh =
xx + yy
2
xx yy
2
cos 2 + xy sin 2
La relation prcdente, dans laquelle xx , yy et hh sont mesures et fix,

permet de calculer xy .
IUTB-LYON1-GMP-DDS
Page 66
TD-2A-S3-F312
On utilise pour mesurer xx , yy et
hh des dispositifs appeles

"rosettes", constitues de trois
jauges colles sur le mme support
et espaces angulairement de 45 ou
120.
Rosette 45 (ou rosette rectangulaire) :
Jauges distinctes
Jauges superposes
Quantits mesures
xx
yy
= ( )
hh = 45

Jy
45
J
IUTB-LYON1-GMP-DDS
Jx
Page 67
TD-2A-S3-F312
De la relation : hh =
xx + yy
2
tenu que =45 :
xx yy
2
cos 2 + xy sin 2 on tire compte
yy
xx
+
+ xy
2
2
xy =
xx + yy
2
Rosette 120 (ou rosette delta) :

y
120
Quantits
mesures
J
P
Jx
xx
= ( hh ) = 120
= ( )
hh = 120

120
De la relation : hh =
tenu que :
xx + yy
2
xx yy
2
cos 2 + xy sin 2 on tire compte
xx
3
3
+
yy
xy
4
4
2
3
3
= xx +
yy +
xy
4
4
2
1) = 120 =
2) = 120
IUTB-LYON1-GMP-DDS
Page 68
TD-2A-S3-F312
xy =
yy =
et
1
(2 + 2 xx )
3
Pour les deux rosettes :

a) Les dformations principales se calculent par :
X =
Y =
xx + yy
xx
2
+ yy
2
1
( xx yy ) 2 + 4 2xy
2
1
( xx yy ) 2 + 4 2xy
2
+
a) La position du repre principal dans le plan se calcule par :
tg =
X xx
xy
= (1, X )
c) Les contraintes principales se calculent par les lois de Lam dans un tat
de contraintes planes :
E
( + Y )
1 - 2 X
E
( + X )
=
1 - 2 Y
X =
Y
IUTB-LYON1-GMP-DDS
Page 69
TD-2A-S3-F312
3. LOIS DE COMPORTEMENT
PROBLME N12
Z
Face Y
Y
X
Dans une plaque suppose indformable, on pratique un logement
paralllpipdique de 10 mm par 10,005 mm par 10 mm. Dans ce logement, on introduit un
cube en aluminium, qui lui est susceptible de se dformer.
Les cts du cube daluminium valent 10 mm.
Les constantes lastiques de laluminium valent : E =70 GPa, =1/3.
Le cube daluminium est comprim sur sa face suprieure avec une pression de 60 MPa.
Le frottement est nglig. Le repre XYZ est le repre principal des contraintes et des
dformations.
1) Le cube vient-il en contact avec la face Ydu logement ?
On supposera pour rsoudre cette question que le cube vient en contact avec la face Y
du logement. Donnez dans ces conditions les valeurs de X, Y et Z, et calculez la
contrainte Y.
Pour que le contact se ralise il faut que la contrainte Y soit ngative.
Dans le cas contraire prendre Y = 0 puisquil ny a pas contact et calculer X, Y et Z.
2) Quelle devrait tre la valeur minimum de Z pour que le contact se ralise ?
IUTB-LYON1-GMP-DDS
Page 70
TD-2A-S3-F312
RPONSES N12
Z
60 MPa
10*10*10
10
10,005
Face Y
10
Y
Ecrivons la loi de HOOKE [ ] =
1+
[ ] tr[ ][I ] sous forme dveloppe et dans les axes
E
E
principaux :
( X ( Y + Z ))
=
X
1
Y = ( Y ( Z + X ))
E
Z = 1 ( Z ( X + Y ))
(1)
(2)
(3)
Supposons que le cube vient en contact avec la face Y du logement . Dans ces conditions :
Y = 5.10 4
X = 0
(1)
0 = X ( Y 60 )
1
1
4
( 60 + X )
5.10 =
3 Y
70.10
3
IUTB-LYON1-GMP-DDS
(2)
Z = 60 MPa
0 = X 3 ( Y 60)
1
35 = Y ( 60 + X )
3
(1)
(2)
Page 71
TD-2A-S3-F312
La rsolution des quations (1) et (2) conduit :
Y = 9,375MPa
Cette contrainte tant positive, le contact ne se ralise pas.

Le contact ne se ralisant pas, prenons Y = 0 et toujours X = 0 et Z = 60 MPa
60
(1)
0 =X
1
1
0 ( 60 + X )
Y =
3
3
70.10
X
Z = 70.103 60 3
(2)
X = 20MPa
6
Y = 381.10
= 762.10 6
Z
La rsolution donne :
(3)
La valeur minimum de Z pour que le contact se ralise conduit aux conditions :
X = 0
(1)
0 =X Z
1
1
4
0 ( Z + X )
5.10 =
3
70.10
3
(3)
Z =
X
3 Z
70.10
3
La rsolution donne :
Y = 0
Y = 5.10 4
(2)
X = 26,25MPa
Z = 78,75MPa
= 10 3
Z
La valeur minimum de Z pour que le contact se ralise : -78,75 MPA
IUTB-LYON1-GMP-DDS
Page 72
TD-2A-S3-F312
PROBLME N13
On donne une matrice de dformation dans les axes principaux :
1000
0
0
200 0
ZZ
0
[] XYZ = 10 6 0
0
Le matriau a comme caractristiques :E= 200 GPa et =0.3.

Calculer la matrice des contraintes correspondante.
a) En contraintes planes.
b) En dformations planes.
RPONSES N13
Ecrivons la loi de LAME:
[ ] = 2G[ ] + tr [ ][I ]
G=
Sachant que :
E
2(1 + )
E
(1 + )(1 2 )
Dans les axes principaux :
X
0
Y
0
0
X
E
0
0
=
(1 + )
0
Z XYZ
P
0
E
( + + )[I ]
0 +
(1 + )(1 2 ) X Y Z
Z XYZ
Y
0
Ecrivons la loi de LAME en contraintes planes :
X
0
0
X
E
0
0
=
(1 + )
0
0 XYZ
P
Y
0
Y
0
0
1 0 0
E
( X + Y + Z )0 1 0
0 +
(1 + )(1 2 )
0 0 1
Z XYZ
Ecrivons la loi de LAME en dformations planes :
X
0
Y
0
IUTB-LYON1-GMP-DDS
0
X
E
0
0
=
(1 + )
0
Z XYZ
P
Y
0
0
1 0 0
E
( X + Y )0 1 0
0 +
(1 + )(1 2 )
0 0 1
0 XYZ
Page 73
TD-2A-S3-F312
Contraintes planes
X
0
0
0
0
1 0 0
1000
E
E
0
(1000 200 + Z )0 1 0106
0
=
200 0 10 6 +
(1 + )
(1 + )(1 2 )
0 0 1
0
Z
0 XYZ
0
Y
0
X
0
Y
0
0
0
0
1 0 0
1000
200.103
0,3 200.103
6
(800 + Z )0 1 0106
0
0
=
200 0 10 +
1,3
1,3 0,4
0 0 1
0
Z
0 XYZ
0
Do le systme :
200.103
0,3 200.103
6
6
X = 1,3 1000.10 + 1,3 0,4 (800 + Z ).10
200.103
0,3 200.103
6
(800 + Z ).106
200
.
10
=
Y
1,3
1,3 0,4
3
200
.
10
0
,
3
200.103
6
(800 + Z ).106
0
.
10
=
+
Z
1,3
1,3 0,4
De lquation (3) nous tirons :

Valeur reporte dans les quations (1) et (2) :
[ ]
P
XYZ
207 0 0
= 0 22 0 MPa
0
0 0
(1)
(2)
(3)
Z = -343.10-6
X = 207 MPa et Y = 22 MPa
[ ]
P
XYZ
0
0
1000
= 0
200
0 10 6
0
343
0
Autre approche possible :

Dans des axes principaux :
X = 1 2 ( X + Y )
E
Y =
( Y + X )
1 2
200.103
6
=
X 1 0,32 (1000 + 0,3( 200 )).10
3
Y = 200.10 ( 200 + 0,3.1000 ).10 6
1 0,32
IUTB-LYON1-GMP-DDS
et
Z =
( X + Y )
1
Z =
0,3
(1000 200).106
0,7
Page 74
TD-2A-S3-F312
Dformations planes
X
0
0
0
0
1000
1 0 0
E
E
6
(1000 200)0 1 0106
0
=
200 010 +
0
(1 + )
(1 + )(1 2 )
0
0 0 1
Z XYZ
0
0
X
0
0
200.103
=
0
1,3
Z XYZ
Y
0
0
0
1000
1 0 0
0,3 200.103
0
6
(800)0 1 0106
200 010 +
1,3 0,4
0
0 0 1
0
0
Do le systme :
200.103
0,3 200.103
6
(800).106
1000
.
10
+
X
1,3
1,3 0,4
3
200.10
0,3 200.103
6
(800).106
200.10 +
Y =
1,3
1,3 0,4
3
0,3 200.10
(800).106
(3)
Z =
1,3 0,4
(1)
(2)
De la rsolution du systme nous tirons : X = 246 MPa ; Y = 62 MPa ; Z = 92 MPa
[ ]
P
XYZ
246 0 0
= 0 62 0 MPa
0
0 92
[ ]
P
XYZ
0
0
1000
= 0
200 010 6
0
0
0
Vrification :
Z = ( X + Y )
IUTB-LYON1-GMP-DDS
Z = 0,3(246 + 62)
Page 75
TD-2A-S3-F312
LOIS DE COMPORTEMENT
1. Relations Contraintes-Dformations
On cherche, dans ce paragraphe, tablir les relations entre la matrice des
contraintes [ ] et la matrice des dformations [ ] . On se place dans les axes
principaux des contraintes et des dformations.
Le matriau est suppos homogne et isotrope.
E est le module de YOUNG 5 du matriau (en MPa).
est son coefficient de POISSON 6 (sans dimension)
G est son module de COULOMB 7 (en MPa).
[]
[]
Lois de HOOKE [] = f([])

Lois de LAM [] = f-1([])
1.1. Elasticit trois dimensions :

1.1.1. Lois de HOOKE 8 :
1) Loi de HOOKE en traction uniaxiale :
X =
X
E
E est le module de YOUNG du

matriau (en MPa).
Y = Z = X =
Facette Z
E X
est le coefficient de POISSON

du matriau (sans dimension).
Facette Y
Facette X
YOUNG Thomas 1773-1829)

POISSON Simon-Denis (1781-1840)
7
COULOMB Charles (1736-1806)
8
HOOKE Robert (1635-1703)
6
IUTB-LYON1-GMP-DDS
Page 76
TD-2A-S3-F312
2) Loi de HOOKE en traction triaxiale :
X/E
Y -Y/E
-X/E -X/E
Y/E
Z -Z/E -Z/E
-Y/E
Y
Z/E
Par superposition, on obtient :
X
Y
Z
1
=
( X ( Y + Z ))
E
1
=
( Y ( Z + X ))
E
1
=
( Z ( X + Y ))
E
X
Z
[ ]=
En introduisant la trace de la matrice [ ] ( tr

prcdentes s'crivent :
+ Y + Z ) les relations
1+
X tr[ ]
E
E
1+
=
Y tr[ ]
E
E
1+
=
Z tr[ ]
E
E
X =
Y
Z
Sous forme matricielle :
[] =
XYZ
1+
[ ] tr[ ][I ]
E XYZ E
1 0 0
[I] = 0 1 0 est la matrice unit.
0 0 1
IUTB-LYON1-GMP-DDS
Page 77
TD-2A-S3-F312
3) Invariance de la loi de HOOKE :

Nous allons montrer que la forme de la loi de HOOKE est indpendante du repre.
Prenons la trace des deux membres de la loi de HOOKE. Compte tenu que la trace
de [I ] vaut 3, on trouve la relation entre les traces de [ ] de [ ] :
1+
[ ] tr[ ][I ]
XYZ
E XYZ E
1+
tr[ ] =
tr[ ] 3 tr[ ]
E
E
[ ] =
tr [ ] =
1 2
tr [ ]
E
La trace de la matrice des dformations tr[] est invariante (cest la dilatation

volumique). E et sont des constantes physiques qui ne dpendent pas du repre o
elles sont calcules puisque le matriau est isotrope. La trace de la matrice [] est
donc invariante. La loi de HOOKE qui a t tabli dans le repre principal
conserve donc la mme forme quel que soit le repre.
[ ] =
xyz
1+
[ ] tr [ ][I ]
E xyz E
4) Module de COULOMB G du matriau :
La relation entre xy et xy tire de la loi de HOOKE permet d'crire :
xy =
2 xy = xy =
2(1 + )
xy
E
1+
xy
E
Or en torsion on a :
xy =
1
xy
G
Do :
G =
E
2(1 + )
G est le module de COULOMB (en MPa) du matriau.
IUTB-LYON1-GMP-DDS
Page 78
TD-2A-S3-F312
1.1.2. Lois de LAM 9:

C'est la prsentation sous forme inverse des lois de HOOKE. La relation de
HOOKE s'inverse facilement compte tenu de la relation entre les traces des
matrices des contraintes et des dformations.
La loi de HOOKE scrit :
[] =
1+
[ ] tr [ ][I]
E
E
Compte tenu de la relation entre les trace de [ ] et [ ] :
tr[] =
1 2
tr[]
E
Nous obtenons en reportant dans la loi de HOOKE :
[] =
[ ] =
1+
E
[]
tr [][I ]
E
E 1 2
E
E
[] +
tr[ ][I ]
1+
(1 + )(1 2 )
Ou encore sous une autre forme en introduisant le module de COULOMB G et le

module de LAM :
[] = 2G[] + tr[][I]
G est le module de COULOMB (en MPa)
est le module de LAM (en MPa)
E
2(1 + )
E
=
(1 + )(1 2 )
G =
LAM Gabriel (1795-1870)
IUTB-LYON1-GMP-DDS
Page 79
TD-2A-S3-F312
NOTATION EN COLONNE POUR [ ] et [ ] :
[]
xyz
xx

yy
zz
=
xy
yz
zx
La loi de LAM [ ] =
xx

yy
zz
=
xy
yz

zx
[]
xyz
E
E
[ ] +
tr[ ][I ] scrit sous forme
1+
(1 + )(1 2 )
dveloppe dans des axes quelconques x,y,z :
xx =
yy =
zz =
xy =
yz =
zx =
E
(1 + ) xx
(1 + )(1 - 2 )
E
(1 + ) yy
(1 + )(1 - 2 )
E
(1 + ) zz
(1 + )(1 - 2 )
E
xy
1+
E
yz
1+
E
zx
1+
( e yy + ezz )
( ezz + exx )
( exx + e yy )
On peut crire les relations prcdentes sous cette forme matricielle :
1

xx
1
0

1
0
yy
1
2
zz
E
0
0
0
=
2
(1 + )(1 2 )
xy
0
0
0
0
yz
zx
0
0
0
0
IUTB-LYON1-GMP-DDS
0
0
0
0
1 2
2
0
0 xx

0 yy
0 zz
xy
0 yz

1 2 zx
2
Page 80
TD-2A-S3-F312
[ ] = [C ][ ]
La matrice symtrique [C] est la matrice des constantes lastiques du
matriau.
0
1

1
0
1
0
1
2
E
0
0
0
[C ] =
2
(1 + )(1 2 )
0
0
0
0
0
0
0
0
La loi de HOOKE [ ] =
0
0
0
0
1 2
2
0
0
0
0
1 2
2
1+
[ ] tr [ ][I ] scrit dans des axes quelconques x,y,z
E
E
sous forme dveloppe :
xx =
yy =
zz =
xy =
yz =
zx =
1
xx yy + zz
E
1
yy ( zz + xx )
E
1
zz xx + yy
E
1+
xy
1+
yz
1+
zx
))
)
))
0
0
0
1
xx
xx
1

0
0
0
yy
yy
zz
zz
0
0
0
1 1
0
0 2(1 + )
0
0
E 0
xy
xy
0
yz
yz
0
0
0
2(1 + )
0
0
0
0
0
2(1 + ) zx
zx
0
[] = [C ]1[]
IUTB-LYON1-GMP-DDS
SYM
Page 81
TD-2A-S3-F312
La matrice symtrique [C]-1 est la matrice inverse des constantes lastiques

[C] du matriau.
[C ] 1
0
0
0
1
0
0
0
1
0
0
0
1
=
0
0 2(1 + )
0
0
E 0
0
0
0
2(1 + )
0
0
0
0
0
2(1 + )
0
1.2. Elasticit deux dimensions :
SYM
1.2.1. Etat de contraintes planes :
0
Y
X
[] = 0
XYZ
0
Lapplication de la loi de HOOKE [ ] =
0
0
1+
[ ] tr [ ][I] montre que la matrice
E
E
des dformations aura la forme :
X
[] = 0
XYZ
0
0
Y
0
0
Z
La matrice des contraintes est plane, mais celle des dformations ne lest pas.
Le dveloppement la dernire ligne de la loi de HOOKE permet de calculer Z :
[] =
1+
[ ] tr [ ][I]
E
E
Z =
tr[ ]
E
Or compte tenu de la relation entre les traces des matrices des contraintes et des
dformations :
tr[] =
E
tr[]
1 2
Z =
( + Y )
1 X
La dformation Z suivant une perpendiculaire au plan des contraintes nest

donc pas indpendante.
IUTB-LYON1-GMP-DDS
Page 82
TD-2A-S3-F312
NOTATION EN COLONNE POUR [ ] et [ ] :
xx
[ ] = yy
xy
xx

[] = yy
xy

zz =
La loi de LAM [ ] =
( xx + yy )
1
E
E
[ ] +
tr[ ][I ] scrit sous forme
1+
(1 + )(1 2 )
dveloppe dans des axes quelconques x,y,z :
E
xx + yy
1 - 2
E
=
yy + xx
1 - 2
E
=
xy
1+
xx =
yy
xy
xx
E
yy = 1 2
xy
0
1
1
0
1
0 0
2

xx
yy

xy
[] = [C ][]
La matrice symtrique [C] est la matrice en contraintes planes des constantes
lastiques du matriau.
0
1
E
1
[C ] =
0
1 2
1
0 0
2
SYM
IUTB-LYON1-GMP-DDS
Page 83
TD-2A-S3-F312
La loi de HOOKE [ ] =
1+
[ ] tr [ ][I ] scrit dans des axes quelconques x,y,z

E
E
sous forme dveloppe :
1
xx yy
E
1
yy xx
=
E
=
( xx + yy )
1
1+
xy
=
xx =
yy
zz
xy
xx
0
1
xx
1
0
yy = E 1
yy
xy
0
0 2(1 + ) xy

[ ] = [C ]1[ ]
La matrice symtrique [C]-1 est la matrice inverse, en contraintes planes, des
constantes lastiques [C] du matriau.
[C ]
0
1
=
1
0
E
0
0 2(1 + )
SYM
1.2.2. Etat de dformations planes :
X
[] = 0
XYZ
0
IUTB-LYON1-GMP-DDS
0
Y
0
0
0
Page 84
TD-2A-S3-F312
Lapplication de la loi de LAM [ ] =

matrice des
X
[] = 0
XYZ
0
E
E
[ ] +
tr[ ][I ] montre que la
1+
(1 + )(1 2 )
0
Y
0
0
Z
contraintes aura la forme :
La matrice des dformations est plane, mais celle des contraintes ne lest pas.
Le dveloppement la dernire ligne de la loi de LAM permet de calculer Z :
[ ] =
E
E
[ ] +
tr[ ][I ]
1+
(1 + )(1 2 )
Z =
E
( + Y + Z )
(1 + )(1 2 ) X
Or compte tenu de la relation entre les traces des matrices des contraintes et des
dformations :
tr[] =
1 2
tr[]
E
On trouve aprs calculs :
Z = ( X + Y )
La contrainte Z suivant une perpendiculaire au plan des dformations nest donc
pas indpendante.
La dmonstration de la forme des matrices symtrique [C ] et [C ] , en
1
dformations planes, est identique celle qui a t faite en contraintes planes. On

montre que :
E

[C ] =
1
0
(1 + )(1 2 )
1 2
0
0
2
[C ]
IUTB-LYON1-GMP-DDS
1 0
1+
=
1 0
E
0
0
2
Page 85
TD-2A-S3-F312
1.2.3. Liaisons contraintes planes - dformations planes :
DFORMATION PLANE (w=0)
CONTRAINTE PLANE ( z=0 )
[u] =
v
xx

[] = yy
xy

zz = ( xx + yy )
( xx + yy )
1
xx
[ ] = yy
xy
1
0
E
1
[C ] =
0
1 2
1
0 0
2
0
1
E

[C ] =
1
0
(1 + )(1 2 )
1 2
0
0
2
[C ] 1
zz =
1 0
1+
=
1 0
E
0
0
2
[C ]
0
1
=
1
0
E
0
0 2(1 + )
On montre facilement que pour passer de l'tat de contraintes planes celui de dformations
planes il suffit de changer :
E
1 2
en
1
E en
IUTB-LYON1-GMP-DDS
Page 86
TD-2A-S3-F312
4. PROBLEMES DE SYNTHESE
PROBLME N14
15
J1
J2
Y
Un exprimentateur maladroit colle une rosette de dformations sur une prouvette en
vue dun essai de traction. La rosette est constitue de 2 jauges 90: J1 & J2. Laxe de
la jauge J2 nest pas comme prvu colle suivant laxe de traction, mais elle fait un angle
de 15 avec cet axe. Il dcide de ne pas reprendre le collage et de raliser lessais de
traction en vue de caractriser le matriau constituant lprouvette.
Pour une force de traction de 10 000 N, les mesures ralises indiquent :
J1 553.10-6
J2 -129.10-6
Le diamtre de lprouvette est de 10 mm.
1) Ecrire la forme littrale de la matrice des dformations dans le repre de la
rosette xyz et dans le repre de lprouvette XYZ.
2) Calculer les dformations longitudinales suivant laxe x, suivant laxe y, ainsi
que la distorsion de langle droit (x,y).
3) Ecrire la matrice des dformations dans le repre (x,y,z).
4) Calculer la contrainte de traction, en dduire le module de YOUNG ainsi que
le coefficient de POISSON du matriau.
IUTB-LYON1-GMP-DDS
Page 87
TD-2A-S3-F312
RPONSES N14
10000 N
15
Jx
15
Jx
Jy
Jy
D=10mm
10000 N
Matrice de dformation, dans le repre de la rosette xyz : [ ]
P
xyz
xx
= yx
0
Matrice de dformation, dans le repre de lprouvette XYZ : [ ]
P
XYZ
xy
yy
0
X
= 0
0
0
0
zz
0
Y
0
0
0
Z
Les dformations longitudinales suivant laxe x, etsuivant laxe y,sont mesures directement
par les jauges Jx et Jy :
xx= 553.10-6 et yy= -129.10-6
La distorsion de langle droit (x,y) peut etre calcule partir de la relation :
IUTB-LYON1-GMP-DDS
Page 88
TD-2A-S3-F312
th =
xx yy
2
sin 2 + xy cos 2
En effet faisons coincider les axes ht avec les axes de lprouvette XY.
XY =
xx yy
2
sin ( 30) + xy cos( 30)
Or la distorsion angulaire dans les axes principaux est nulle : XY=0
0=
xy
xx yy
sin ( 30 ) + xy cos( 30) ;
2
553 + 129
sin ( 30 )
2
=
;
cos( 30 )
0=
553 + 129
sin ( 30 ) + xy cos( 30)
2
xy= -196,87.10-6do la distorsion angulaire : xy=2xy= -393,74.10-6
Do la matrice de dformation dans le repre xyz de la rosette :
[ ]
P
xyz
553 197 0
= 197 129 0 d
0
0
zz
X =
Calcul de la contrainte de traction :
N X 10000
=
= 127,3MPa
52
A
Dtermination du module de YOUNG :

Utilisons la loi de HOOKE en traction pure : X = E X
Nous devons donc calculer la dformation longitudinale de lprouvette :
xx + yy
xx yy
cos( 30) + xy sin ( 30)

2
553 129 553 + 129
X =
+
cos( 30 ) 197 sin ( 30 ) = 606d
2
2
127,3
E= X =
E=210100 MPa
X 606.10 6
X =
Dtermination du coefficient de POISSON :

Utilisons la loi de POISSON : =
Y
X
Nous devons donc calculer la dformation transversale de lprouvette :

La trace de la matrice de dformation tant le premier invariant :
xx + yy + zz = X + Y + Z = 553 129 + zz = 606 + Y + Z

182
= Y =
X 606
Dtermination du module de COULOMB : G =
IUTB-LYON1-GMP-DDS
E
2(1 + )
Y = 182 d
=0,3
G= 80807 MPa
Page 89
TD-2A-S3-F312
PROBLME N15
On considre un point P la surface d'un
corps en un endroit o ne s'exerce aucune
force extrieure. On mesure dans le plan
tangent la surface, les dilatations linaires
unitaires suivant les trois axes x, et
choisis comme l'indique la figure ci-contre
(rosette 120). Les rsultats des mesures
sont :
jauge J1 : 200 .10-6
jauge J : 80 .10-6
jauge J: -40 .10-6
1) Calculer xy et yy.
2) Dterminer graphiquement avec le
cercle de MOHR les dilatations principales
J
120
J1
-120
J
et les axes principaux.

3) E = 200000 MPa , = 0.3 calculer la matrice des contraintes dans les axes principaux.
RPONSES N15
y
120
J
P
Jx
120
IUTB-LYON1-GMP-DDS
Page 90
TD-2A-S3-F312
Les jauges Jx, J et J mesurent directement les dilatations unitaires suivant les axes x
Les quantits mesures
xx
= ( hh ) = 120
= ( )
hh = 120

Nous devons dpouiller la rosette, afin de dterminer la dilatation linaire suivant laxe y,
ainsi que la demi distorsion angulaire xy.
Utilisons la relation : hh =
xx + yy
2
xx yy
2
cos 2 + xy sin 2 en faisant coincider h avec la
direction puis avec la direction .
xx + yy xx yy
cos(2.120 ) + xy sin (2.120)

+
= ( hh ) =120 =
2
2
xx + yy xx yy
= ( hh )
+
cos( 2.120 ) + xy ( 2.120)
=120 =
2
2

3
3
=
= xx + yy
xy
xy
4 4
2
3
la rsolution du systme :
1
xx 3
3
yy = (2 + 2 xx )
=
+ yy +
xy
3
4 4
2
xy =
1
40 80
= 69,28d yy = (2.80 2.40 200) = 40d
3
3
La matrice de dformation dans les axes xyz est donc : [ ]
P
xyz
69,28 0
200
= 69,28
40
0 d
0
0
zz
Dtermination graphique, laide du cercle de MOHR des dformations, des lments

principaux.
La technique est identique celle que nous avons dvelopp sur les contraintes. Les deux
approches contraintes et dformations sont en effet rgient par les mmes lois
quadratiques.
Les directions x et y sont orthogonales, et nous devons construire le cercle partir de la
connaissance des trois valeurs
xx= 200 10-6
IUTB-LYON1-GMP-DDS
yy= -40 10-6
xy= -69,28 10-6
Page 91
TD-2A-S3-F312
[ x ]/ xyz
1)
2)
3)
4)
5)
200
= 69,28 d
0
[ ]
y / xyz
69,28
= 40 d
0
Tracer les axes hh et th et adopter une chelle 1cm=10d

Positionner les points mx (200,-69.28) et my (-40,69.28).
centre du cercle IZ.
Tracer le cercle CZ de centre IZ et de rayon IZmx.
matrialisent les facettes principales du plan des dformations. OmX = 219 = X
OmY = 59 = Y
6)
2 = 30
matrialisent les directions associes et en mesurant les valeurs correspondantes,
compte tenu de lchelle.
th
my(-40,69,28)
69,28
Z/2
y
mY(-59)
-69,28
hh
200
IZ
-40
30
mX(219)
mx(200,-69,28)
CZ
20 d
IUTB-LYON1-GMP-DDS
Page 92
TD-2A-S3-F312
120
120
Jx
15
X
La matrice de dformation dans les axes XYZ est donc : [ ]
P
XYZ
0
0
219
= 0 59 0 d
0
0 Z
Matrice des contraintes dans les axes principaux :
X = 1 2 ( X + Y )
;
E
Y =
( Y + X )
1 2
2.105
6
=
X 1 0,32 (219 0,3.59)10

;
5
2
.
10
6
Y =
( 59 + 0,3.219)10
1 0,32
La matrice des contraintes dans les axes XYZ est donc : [ ]
P
XYZ
IUTB-LYON1-GMP-DDS
X = 44,24MPa
Y = 1,47 MPa
44,24 0 0
= 0
1,47 0 MPa
0 0
0
Page 93
TD-2A-S3-F312
PROBLME N16
Considrons un point P situ la surface d'un cylindre
de rayon R soumis une force de traction P et un
couple de torsion M. Le plan est tangent en P au
cylindre. Le matriau est isotrope de constantes
lastiques E et .
1) Calculer les longations unitaires (x,)= +45) et
(x,)= -45) suivant les directions et dues :
a) P ( et en fonction de P, R, E et )
b) M ( et en fonction de M, R, E et )
c) P et M
2) En dduire les expressions de P et M en fonction de
R , E , , et .
3) Application numrique :
R = 10 mm, E= 200 GPa, = 0,3, e = e' = 240 MPa,
= 100 10-6, = 50 10-6
Calculer P en N, M en Nm et le coefficient de scurit
en appliquant le critre de VON MISES.
C
P
P
C
RPONSES N16
La section G est soumise un effort normal de traction P,

et un couple de torsion Mt=C .
Le plan est tangeant en P au cylindre.
Il est donc orient par la normale extrieure z.
Les axes xy et sont donc situs dans le plan .
Pour calculer , utilisons la relation :
hh =
xx + yy
2
xx yy
cos 2 + xy sin 2
en faisant coincider h avec .
xx + yy
2
xx yy
2
cos 90 + xy sin 90
P
z
45
45
Do :
De mme :
xx + yy
2
xx + yy
2
+ xy
xy
P
C
IUTB-LYON1-GMP-DDS
Page 94
TD-2A-S3-F312
Calcul des longations unitaires suivant et due P :

Le cylindre est donc soumis de la traction pure. Les axes x, et y sont donc principaux, et la
demie distorsion angulaire xy est nulle. La loi de HOOKE, Y = E Y ainsi que la loi de
POISSON, =
yy =
yy
E
X
nous permettent dcrire :
Y
Do :
xx =
P
P
=
;
EA ER 2
xx + yy
2
+ xy =
xx + yy
2
yy
P
E R 2
xy = 0
P
P
+
2
ER
ER 2 = P (1 )
2
2 ER 2
La demie distorsion angulaire xy tant nulle, le rsultat est identique pour la direction .
P
(1 )
2 ER 2
P
(1 )
2 ER 2
Calcul des longations unitaires suivant et due Mt=C :

Le cylindre est donc soumis de la torsion pure. Les axes , et sont donc principaux, La loi
de HOOKE en torsion, = G et compte tenue que dans les axes xy le point P est en
cisaillement pur, nous pouvons crire :
= G = xy = G xy ;
xy
; =
M t R CR
2C
=
= 3 ;
4
R
IG
R
2
2C
3
C
=
=
= R =
;
2
2G
2G GR 3
xy
xy
Pour la direction :
xx + yy
2
IUTB-LYON1-GMP-DDS
2C (1 + )
ER 3
xx = yy = 0
+ xy = xy ;
xy = xy ;
La relation liant les 3 constantes lastiques G =
xx + yy
xy = 2 xy ;
C
GR 3
C
GR 3
E
entraine :
2(1 + )
2C (1 + )
ER 3
Page 95
TD-2A-S3-F312
Calcul des longations unitaires suivant et due P et Mt=C :

Il nous suffit de superposer les rsultats prcdents :
P
(1 ) + 2C (1 +3 )
2
2 E R
E R
P
2C (1 + )
(
)
1
2 E R 2
E R 3
Calcul de P et Mt=C en fonction de R, E, , , et :

Nous obtenons donc un systme de 2 quations, quil suffit de rsoudre :
2C (1 + )
P
= 2 ER 2 (1 ) + ER 3
P
=
(1 ) 2C (1 +3 )
2
2 ER
ER
En ajoutant puis en retranchant les 2 quations, nous
tirons :
ER 2
( + )
P=
1
ER 3
( )
C=
4(1 + )
Application numrique :
P=
3
2
ER 2
( + ) = 200.10 10 (100 + 50).106 = 13464 N
1
1 0,3
P=13464 N
C=
200.103 103
ER 3
(
(100 50).10 6 = 6041,5 Nmm
) =
4(1 + )
4(1 + 0,3)
C=6,04 Nm
Pour calculer le coefficient de scurit, il nous faut calculer les contraintes :
yy =
13464
P
=
= 42,86MPa
2
R
10 2
2C 2.6041
=
= 3,85MPa
R 3
103
Appliquons le critre de VON MISES dans le cas des poutres : =
e
2 + 3 2
IUTB-LYON1-GMP-DDS
240
42,862 + 3.3,852
= 5,53
e
e
=
VM
2
qu
+ 3 2
= 5,53
Page 96
TD-2A-S3-F312
PROBLME N17
z
P
x
45
Soit un point P situ la surface d'un corps en un

endroit o ne s'exerce aucune force extrieure.
La matrice des contraintes dans le repre xyz (o
z est la normale extrieure au corps, et o x et y
sont dans le plan tangent au corps) vaut en MPa:
80 15 0
= 15 40 0
0
0
0
[ ]
xyz
1) Etude des contraintes au point P.

a) Reprsenter l'tat de contraintes sur les faces d'un cube infiniment petit isol en P.
b) Calculer la contrainte normale s'exerant sur une facette dont la normale est la
bissectrice de xy.
c) Dterminer les lments principaux des contraintes.
d) Calculer le cisaillement maximum en module au point P.
e) On veut un coefficient de scurit de 3 par rapport au critre de VON MISES.
Calculer la limite lastique du matriau.
2) Etude des dformations au point P.
On colle une rosette 45 x,,y au point P. Calculer la valeur des dilatations que l'on doit
mesurer suivant les axes x, et y sachant que E = 200 GPa et = 0.3.
RPONSES N17
P
x
45
[ ]
P
xyz
IUTB-LYON1-GMP-DDS
80 15 0
= 15 40 0 MPa
0
0 0
Page 97
TD-2A-S3-F312
Etude des contraintes au point P.

Reprsentation de ltat des contraintes en P sur les faces dun cube infiniment petit isol
dans les axes xyz.
z
P
x
15
15
40
80
La facette z est perpendiculaire au plan tangent xy de la surface du corps. En labscence de
force extrieure agissant en P, le vecteur contrainte z est nul, do la reprsentation dans le
plan xy.
Soit une facette dcale de 45 par rapport x.
Calculons la contrainte normale sexercant sur cette facette
Compte tenu de la relation : hh = [h][ ][h]
t
xx + yy
2
xx + yy
2
xx yy
2
+ xy =
15
15
En faisant coincider h avec :
40
cos(2.45) + xy sin (2.45)
80
80
45
15
80 + 40
15 = 45MPa
2
15
40
= 45 MPa
Elments principaux des contraintes :

Etat plan de contraintes : la facette z est donc principale et Z=0. Les 2 autres valeurs
xx + yy
1
( xx yy )2 + 4 xy2
2
2
80
40
1
+
(80 40)2 + 4.152
YX =
m
2
2
X=85MPa Y=35MPa
propres sont donnes par : YX =
IUTB-LYON1-GMP-DDS
Page 98
TD-2A-S3-F312
La direction propre X est calcule par tan =
r r
X xx
avec = x , X
xy
r r
85 80
do = 1843 = x , X
15
85 0 0
P
La matrice des contraintes dans les axes principaux : [ ]XYZ = 0 35 0 MPa
0 0 0
tan =
Le cisaillement maximum est gal la plus grande demie diffrence des contraintes
principales. Dans notre cas : max = Y =
X Z
2
85 0
2
max=42,5 MPa
VM
VM
qu
en 2D compte tenu que Z=0 scrit : qu =
X2 + Y2 X
e
Le coefficient de scurit : = VM
e = X2 + Y2 X Y
qu
La contrainte
e = 222 MPa
e = 3 852 + 352 85 .35
Etude des dformations au point P

Connaissant la matrice des contraintes, et les constantes lastiques du matriau nous utilisons
les lois de HOOKE pou calculer la matrice de dformation:
P
[ ]Pxyz = 1 + [ ]Pxyz tr [ ][I ]

E
Jx
Jy
y
J
45
xx
yx
zx
xy xz
80 15 0
1 0 0
1 + 0,3
0,3
(80 + 40)0 1 0
15 40 0
yy yz =
3
3
200.10
200.10
0
0 0 1
0 0 xyz
zy zz xyz
IUTB-LYON1-GMP-DDS
Page 99
TD-2A-S3-F312
Ecrivons les quations 1, 2, et 4 :
6
xx = 200.103 (80 0,3.40) = 340.10
1
(40 0,3.80) = 80.106
yy =
3
200.10
xy = 1 + 0,3 ( 15) = 97,5.10 6
200.103
hh =
xx + yy
2
xx yy
2
cos 2 + xy sin 2
340 + 80
97,5 = 112,5d
2
La mesure effectue par la jauge x : 340 d

La mesure effectue par la jauge y : 80 d
La mesure effectue par la jauge : 112,5 d
Remarque :
La relation de HOOKE [ ]xyz =
P
1+ P
[ ]xyz tr [ ][I ] est valable quelque soit le systme daxes
E
E
puisque matricielle. Il est donc possible de passer par les axes principaux.
X
0
Y
0
0
1 + 0,3
=
0
200.103
Z XYZ
85 0 0
1 0 0
0 35 0 0,3 (85 + 35)0 1 0
200.103
0 0 0 xyz
0 0 1
En crivant les quations 1, et 2 il est possible de calculer: X

valeurs mesures par les 3 jauges, en utilisant la relation :
hh =
IUTB-LYON1-GMP-DDS
X + Y
2
X Y
2
Y , puis den dduire les
cos 2
Page 100
TD-2A-S3-F312
PROBLME N18
50
Deux forces P et Q sont appliques simultanment

au centre de gravit de la section droite libre
d'une poutre encastre. Une rosette rectangulaire
colle au point M indique :
xx = -30 10-6
= 205 10-6
yy = 100 10-6
Le module de YOUNG du matriau vaut 200 GPa.
1) Dterminer la matrice des contraintes au
point M dans le repre xyz en fonction de P, Q et
A.
On rappelle que la contrainte tangentielle due Q
vaut 3Q/2A au point M.
150
=
=
y2
M
z3
45
1
x
2) Dterminer les expressions littrales de xx ,

yy et en fonction de P, Q, A, E et .
3) En dduire les valeurs numriques :
a) du coefficient de POISSON .
b) des forces P et Q (en kN).
RPONSES N18
50
300
M
150
IUTB-LYON1-GMP-DDS
Page 101
TD-2A-S3-F312
50
300
M
150
=
Q
y
z
Nous devons dterminer le Torseurde Section

sexercant dans la section G de la barre o est colle la rosette M
Considrons le systme daxes G/xyz de la thorie des poutres

Torseur de section en G :
Nx = +P
RG ' = FG
0
Tz = +Q
M G ' = G ' G FG M y = 300Q
Nx=P
MY=300Q
TZ=Q
y
z
La section G est donc soumise de la traction-flexion-cisaillement
Rappels du cours F213 :(Effort Normal-Flexion)
Nx MY P
Z
+
A
IY
N
P
= x =
A
A
La contrainte normale en un point P quelconque de la section G vaut : xxP =

Or le point M est situ sur laxe neutre de flexion : ZM=0 do xxM
Rappels du cours F213 : (Cisaillement)
TZ
b
max
3T
= Z
2A
G
z
G
z
max
g
C
Z
IUTB-LYON1-GMP-DDS
Page 102
TD-2A-S3-F312
Au point M, dans les axes xyz lis la rosette :

La contrainte normale due leffort de traction P :
G
yy =
P
A
La contrainte tangentielle due au cisaillement Q :
xy =
0
3Q
=
2A
0
Do la matrice des contraintes au point M : [ ]xyz

M
3Q
2A
3Q
0
2A
P
0 MPa
A
0 0
La mesure ffectue par les jauges donne :
xx = -30 10
-6
= 205 10-6
Jx
yy = 100 10-6
Utilisons la loi de HOOKE :
IUTB-LYON1-GMP-DDS
xx
yx
zx
xy xz
yy yz
zy zz xyz
45
[ ]Pxyz = 1 + [ ]Pxyz tr [ ][I ]

E
Jy
0
1 + 3Q
=
E 2A
0
3Q
0
2A
1 0 0
P
P
0
0 1 0
A
E A
0 0 1
0 0
xyz
Page 103
TD-2A-S3-F312
Ecrivons les quations 1, 2, et 4 :
1
P
xx = E A
1 P
yy =
E A
+
3Q
1
xy
E 2A
Nous en dduisons :
hh =
xx + yy
2
xx yy
cos 2 + xy sin 2
xx + yy
2
+ xy
1
P 1 P
+
E
A E A 1 + 3Q
=
+
2
E 2A
xx =
Do :
yy =
=
P
EA
P
EA
P
(1 ) + 3Q (1 + )
2 EA
2 EA
Application Numrique :
xx =
P
xx = yy
EA
=0,3
yy =
P
P = yy EA
EA
P= 150kN
P
(1 ) + 3Q (1 + )
2 EA
2 EA
Q= 130,77 kN
IUTB-LYON1-GMP-DDS
Page 104
TD-2A-S3-F312
PROBLME N19
Jy
Jx
Une rosette 45 (x,,y,) quatre jauges (avec une quatrime jauge

perpendiculaire la jauge centrale ) est colle en un point P situ la surface d'un
corps en un endroit o ne s'exerce aucune force extrieure (ltat de contraintes est
plan).
Les mesures indiques par les trois jauges Jx, J, et J sont :
xx = 10-3
= + 1.2 10-3
= - 0.4 10-3
La jauge JY ne donne aucune indication.
1) Montrez que : xx + yy = + . En dduire la valeur de yy.
2) Dterminez lexpression littrale de la dformation linaire unitaire en
fonction de xx , yy et xy . En dduire la valeur numrique de xy.
3) Le matriau a comme caractristiques : E= 200 GPa et =1/3.
Calculez la dformation linaire unitaire zz (on rappelle que ltat de contraintes
est plan), puis la matrice des contraintes [ ] dans le repre xyz.
4) Dterminez, avec le cercle de MOHR, les lments principaux des
contraintes dans le plan xy (prcisez langle entre x et X) et la valeur du
cisaillement maximum en module.
5) Calculez la contrainte quivalente de VON MISES.
IUTB-LYON1-GMP-DDS
Page 105
TD-2A-S3-F312
RPONSES N19
y
45
45
Jy
J
Jx
La forme de la matrice de dformation est : [ ]
P
xyz
xx
= yx
0
xy
yy
0
0

P
0 ou [ ]z =
0
zz
0
0
zz
Utilisons le premier invariant linaire de la matrice : I1 = tr [ ] .
xx + yy + zz = + + zz
Do :
xx + yy = + ;
yy = + xx ;
yy = (1,2 0,4 1)10 3
yy=-0,2.10-3
1
cos 45
1
hh = [h][ ][h] avec h coincidant avec : = [ ][ ][ ] , et [ ] = sin 45 =

1
2
0
0
t
Do :
xx
1
[1 1 0] yx
= [ ][ ][ ] =
2
0
t
xx + yy
2
+ xy
IUTB-LYON1-GMP-DDS
xy =
xx + yy
2
xy
yy
0
0 1
+ yy
1
= xx
+ xy
0 1
2
2
zz 0
1 0,2 3
= 1,2
10
2
xy=0,8.10-3
Page 106
TD-2A-S3-F312
Ecrivons la loi de LAME en contraintes planes :
[ ] = 2G[ ] + tr[ ][I ]

Sachant que :
G=
E
2(1 + )
[ ] =
xx xy xz
200.103
yy
yz =
yx
1
zx zy zz 1 +
(1 + )(1 2 )
E
E
[ ] +
tr [ ][I ]
(1 + )
(1 + )(1 2 )
1
0,8
0
1 0 0
1
200.103
0,8 0,2 0 10 3 + 3
(0,8 + zz )0 1 0103
1
1
1 + 1 2
0 0 1
0
zz
0
3
3
En labscence de force extrieure sappliquant au point, la contrainte zz est nulle.
xx xy
yx yy
0
0
0
0
1 0 0
100 80
0 = 1,5 80 20 0 + 1,5(80 + zz )0 1 0
0 0 1
0
zz
0
0
zz=-0,4.10-3
Ecrivons la troisime quation : zz = 0 = 1,5 zz + 1,5(80 + zz )
xx xy
yx yy
0
0
0
0
1 0 0 210 120 0
100 80
0 = 1,5 80 20
0 + 60 0 1 0 = 120 30 0 MPa
0 0 1 0
0
40
0
0
0 0
[ ]
P
xyz
[ ]
P
XYZ
La contrainte
VM
qu
210 120 0
= 120 30 0 MPa
0
0 0
0
0
270
= 0 30 0 MPa
0
0
0
VM
qu
en 2D compte tenu que Z=0 scrit :
= 270 + 30 + 270.30
2
IUTB-LYON1-GMP-DDS
VM
qu
eqVM=286,2
2
X
2
Y
MPa
Page 107
TD-2A-S3-F312
Elments principaux, en utilisant le cercle de MOHR
Les facettes x et y sont orthogonales :

Tracer les axes hh et th et adopter une chelle 1cm=20MPa
Positionner les points mx (210,120) et my (30,-120)
centre du cercle IZ
Tracer le cercle CZ de centre IZ et de rayon IZmx
matrialisent les facettes principales du plan des contraintes. Om X = 270 = X
2 = 532
OmY = 30 = Y
matrialisent les facettes associes et en mesurant les valeurs correspondantes, compte
tenu de lchelle
th
m(120,150)
120
mx(210,120)
Zmax =150
532
Y
mY(-30)
IZ
30
hh
X
210
mX(270)
-120
my(30,-120)
CZ
m(120,-150)
20 MPa
IUTB-LYON1-GMP-DDS
Page 108
TD-2A-S3-F312
5. CONTACT DE HERTZ
PROBLME N20
12
10
110
100
80
70
Le roulement une range de rouleaux est un roulement SNR NU 1014 M. Le

diamtre des rouleaux est de 10mm. La charge statique que peut supporter le
roulement est donne par le constructeur : P=67000 N. On admet que cette charge
sapplique sur trois rouleaux.
Le matriau est un acier 100C6 (E=205 GPa, =0.3). Cet acier peut supporter, daprs
la norme ISO, des pressions en surface de lordre de 4000 MPa sans matage. Calculez
la pression de HERTZ pmax entre le rouleau et les bagues intrieure et extrieure.
RPONSES N20
2
10
1
3
12
110
100
Nous sommes dans le cas du contact de deux cylindres
80
70
IUTB-LYON1-GMP-DDS
Page 109
TD-2A-S3-F312
Cas du contact entre le rouleau 1 et la bague intrieure 3 :

La surface de contact est un rectangle troit.
La pression pmax au centre du rectangle est donne par :
10
1
qE
r
1 1 1
= +
r r1 r3
pmax = 0,418
avec
12
1 1 1
9
= +
=
r 5 40 40
100
80
Si le contact seffectue suivant une gnratrice de longueur L :
q=
P 67000
=
L
12
Or le contact seffectue dans notre cas sur 3 rouleaux : q =
pmax = 0,418
67000 9 205.103
qE
= 0,418
r
3 12 40
P 67000
=
L 3 12
pmax=3872 MPa
Cas du contact entre le rouleau 1 et la bague extrieure 2 :
La surface de contact est un rectangle troit.
2
10
1
La pression pmax au centre du rectangle est donne par :
qE
pmax = 0,418
r
1 1 1
avec
=
r r1 r2
12
100
80
1 1 1
9
=
=
r 5 50 50
Si le contact seffectue suivant une gnratrice de longueur L :
q=
P 67000
=
L
12
Or le contact seffectue dans notre cas sur 3 rouleaux : q =
pmax = 0,418
IUTB-LYON1-GMP-DDS
67000 9 205.103
qE
= 0,418
r
3 12 50
P 67000
=
L 3 12
pmax=3464 MPa
Page 110
TD-2A-S3-F312
PROBLME N21
Dterminer pmax au contact de la roue et

du rail. Leffort de compression est
P=5000 N. Le matriau est un acier
(E=210 GPa, =0.3).
1
500
250
RPONSES N21
P
Le rayon de la roue porteuse r1 est gal :500mm.
Le rayon de courbure de la surface du rail r2 vaut 250mm.
Nous devons chercher les courbures principales 1/r et 1/r

de la roue 1 et du rail 2.
1
500
Les courbures principales tant les courbures maximum

et minimum qui se trouvent dans deux plans orthogonaux.
Dans notre cas langle entre les plans contenant

les courbures principales est gal : /2.
IUTB-LYON1-GMP-DDS
250
Page 111
TD-2A-S3-F312
1
r1=
Z
Z
r1=500
r2=250
r2=
La surface de contact est une ellipse dont les demi axes valent :
Y
a = 3
Pm
n
b = 3
Pm
n
2a
2b
Avec : m =
Do : m =
4
1 1 1 1
+ + +
r1 r '1 r2 r '2
4
1
1
+
500 250
et
2000
3
n=
n=
4E
3 1 2
4 210.103
= 307692,3
3 1 0.32
Les constantes et sont dtermines par le tableau partir de la connaissance de :
cos =
IUTB-LYON1-GMP-DDS
B
A
Page 112
TD-2A-S3-F312
En sachant que
A=
1
B=
2
1
B=
2
2
3
=2
= 3.10 3
m
2000
2
1 1 1 1
1 1 1 1
+ + 2 cos 2
r1 r '1 r2 r '2
r1 r '1 r2 r '2
1 1
2
1 1 1
1 1 1 1 1
1
3
+ +
= =
cos =
= 10
r
r
r
r
r
r
r
r
2
2
2
250
500
1 2
1 2
1 2
2 1
= 1,284
= 0,802
B
103
1
= = 0,333
cos = =
3
A 3.10
3
Do :
= 70
On en dduit donc le grand axe et le petit axe de lellipse !
a = 3
5000 2000
Pm
= 2,8
= 1,2843
307692,3 3
n
b = 3
et
5000 2000
Pm
= 1,8
= 0,8023
307692,3 3
n
Ltat de contrainte est un tat triaxial de compression.
Cest au centre de lellipse que les contraintes sont maximum.

z
La contrainte de compression maximum

est suivant laxe Z :
pmax = 1,5
pmax=Z=1,5P/ab
P
5000
= 1,5
475MPa
ab
2,84 1,8
X
x
IUTB-LYON1-GMP-DDS
Y
y
Page 113
TD-2A-S3-F312
ANNEXES MATHEMATIQUES
1. DRIVES.DIFFRENTIELLES
1.1. Drives, diffrentielle dune fonction :
1) Driv, diffrentielle dune fonction y dune variable x :
y
Considrons la fonction y = x 2 .
La drive y de la fonction y = x 2
y=f(x)=x2
vaut y' = 2x = tg .
Tangente en Mo la
courbe y
Pour x=x0=1 la drive de la fonction

y = x 2 vaut :
y' = tg = 2
= 63.4 0
M0
d
La drive peut scrire sous la forme :
x0+dx
dy
y' = tg =
dx
x0=1
dy = y'dx est la diffrentielle de la fonction y. Cest laccroissement de la

fonction y prs quand la variable x augmente de dx.
2) Drivs, diffrentielle dune fonction z de deux variables x et y :

Soit z = f( x, y ) une fonction de deux variables x et y.
Les deux drives partielles par rapport x et y de la fonction z se notent :
z
x
z
z'y =
y
z'x =
La diffrentielle dz de la fonction z se calcule par :
dz =
IUTB-LYON1-GMP-DDS
z
z
dx +
dy
x
y
Page 114
TD-2A-S3-F312
Cest laccroissement de la fonction z

et quand la variable y augmente de dy.
prs quand la variable x augmente de dx
3) Drivs, diffrentielle dune fonction w de trois variables x, y et z :

Soit w = f( x, y, z) une fonction de trois variables x, y et z.
Les trois drives partielles par rapport x, y et z de la fonction w se notent :
w
x
w
w 'y =
y
w
w 'z =
z
w 'x =
La diffrentielle dw de la fonction w se calcule par :
dw=
w
w
w
dx+ dy+ dz
x
y
z
Cest laccroissement de la fonction w prs quand la variable x augmente de dx,

quand la variable y augmente de dy et quand la variable z augmente de dz.
1.2. Drives, diffrentielle dun vecteur :

Drives d'un vecteur:
Plaons nous dans le plan xy (c'est la
mme chose dans l'espace avec une
composante w en plus pour le vecteur et
une variable supplmentaire z).
y
u(P)
v
y
r
j
O
Les drives
r
i
r
i
x
u
x
et
s
j
x
r
r
r
r
u ( P ) = u ( x, y ) = u( x, y )i + v ( x, y ) j
r
r
r
r
i
u
j
v r
u
j + v
i + u
+
=
x x
x
x
x
sont nulles car la base xy est fixe. Do :
r
u r
v r
u
=
i +
j
x
x
x
IUTB-LYON1-GMP-DDS
Page 115
TD-2A-S3-F312
De mme, pour la driv par rapport y :
r
u r
v r
u
=
j
i +
y
y
y
Diffrentielle d'un vecteur:
Plaons-nous dans le plan (c'est la mme chose dans l'espace avec une composante
w en plus et une variable supplmentaire z).
r
r
r
r
r
r
r
u ( P ) = u ( x , y ) = u ( x, y )i + v ( x, y ) j
du = du i + dvj
r
r
u
u
r
v
v
du = (
dx +
dy )i + (
dx +
dy ) j
x
y
x
y
dv
du
u r v v
u r
r
v r
du = (
i +
j )dx + (
i +
j )dy
x
x
y
y
Or (voir drives d'un
vecteur) :
r
u r
u
v r
=
i +
j
x
x
x
dv
u(Q)
Q
y+dy
r
u r
v r
u
=
i +
j
y
y
y
du
P
y
r
r
r
u
u
du =
dx +
dy
x
y
u(P)
r
j
O
r
i
du
x
x+ x
La diffrentielle d'un vecteur dans le plan se calcule comme la diffrentielle

totale dz d'une fonction de deux variables z=f(x,y).
dz =
IUTB-LYON1-GMP-DDS
z
z
dx +
dy
x
y
Page 116
TD-2A-S3-F312
2. CALCUL MATRICIEL
2.1. Dfinition
La notion de matrice t introduite par le mathmaticien Anglais SYLVESTER 10.
Une matrice peut tre considre comme un tableau contenant des nombres ou des
symboles.
m12=2
[M] =
0 2 4 6 8
1 3 5 7 9
m25=9
Colonne 2
La matrice [M] est

rectangulaire dordre(2,5)
(2 lignes et 5 colonnes)
Ligne 1
2.2. Exemples de matrices :
Matrice carre symtrique dordre 3 :
[M ]
1= 2 3
=
= 2 4 5
=
3 5 6
Diagonale
principale
Symtrique
10
SYLVESTER James Joseph (1814-1897)
IUTB-LYON1-GMP-DDS
Page 117
TD-2A-S3-F312
Matrice carre diagonale dordre 3 :
[M ]
= 0
0
0
4
0
0
6
Matrice unit :
1 0 0
[I] = 0 1 0
0 0 1
Matrice vecteur :
1
[V] = 2
3

1
[i] = 0
0

vecteur
vecteur unitaire de laxe des x
Matrice rotation d'un vecteur dans le plan xy :
y
Vy
Uy
Vx
[V] = [R][U] avec :
IUTB-LYON1-GMP-DDS
x
Ux
cos sin 0
[R] = sin cos 0
0
0
1
Page 118
TD-2A-S3-F312
2.3. Oprations sur les matrices :
2.3.1. Multiplication :
1 2
[A] = 3 4
5 6
[B] =
1 2
2 3
La matrice [A] est dordre (3,2), la matrice [B ] est dordre (2,2).

Le matrice produit [P ] de [A] par [B ] dans cet ordre scrit :
[P] = [A][B]
Un terme quelconque pij de la matrice produit [P ] sobtient en faisant la somme des
produits des termes de la ligne i de [A] par les termes de la colonne j de [B ] .
p11 = 1 1 + 2 2
1 2 + 2 3
5 4
3 2 + 4 3 = 11 6
17 8
5 2 + 6 3
1 1 + 2 2
1 2
1 2
[P] = 3 4
= 31+ 42
2 3
5 1 + 6 2
5 6
ordre (3,2) ordre (2,2)
ordre (3,2)
Lordre le la matrice [P ] sera donc (3,2).
Pour que le produit des deux matrices soit possible il faut que le nombre de
colonnes de [A] soit identique au nombre de lignes de [B]. Avec lexemple
choisi le produit [B][A] est impossible.
2.3.2. Transposition :
Il suffit de changer les lignes en colonnes.
1 2
[M] = 3 4
5 6
IUTB-LYON1-GMP-DDS
[M] =
1 3 5
2 4 6
Page 119
TD-2A-S3-F312
Si la matrice est symtrique (donc carre), elle est gale sa transpose.

Attention la transposition d'un produit de matrices !!
[P ]
= [M ] [N ]
( m,n )
( m,p ) ( p,n )
[P ]= t [N]t [M ]
( n,m )
t
( n,p ) ( p,m )
[P ]= t [M] t [N]
( n,m )
( p,m ) ( n,p )
Explication : Pour que le produit de deux matrices soit possible il faut que le
nombre de colonnes de la premire matrice et le nombre de lignes de la
seconde soit identique.
2.3.3. Inversion :
La matrice inverse [M]
dune matrice [M] est telle que :
[M][M]1
= [M ]
[M] = [I]
1 0 0
[I ] = 0 1 0 est la matrice unit.
0 0 1
On montre en algbre linaire que linverse de [M ] se calcule par :
[M]
[coM]
M
|M| est le dterminant de la matrice.
[coM] est la transpose de la comatrice de [M] (ou adjointe de [M]).
[M] =
mij
M
ij
Le cofacteur de rang i,j se calcule par : comij = ( 1) i + j Mij o |Mij| est le mineur.
IUTB-LYON1-GMP-DDS
Page 120
TD-2A-S3-F312
Remarque : la matrice inverse d'une matrice rotation est gale sa transpose
( [R] = [R] ). En effet, la rotation inverse d'une rotation d'angle est une rotation
1
d'angle -.
[R] 1
cos sin 0
[R] = sin cos 0
0
0
1
cos sin 0
cos( ) sin( ) 0
= sin( ) cos( ) 0 = sin cos 0 = [R]

0
0
0
1
0
1
2.3.4. Diagonalisation d'une matrice symtrique d'ordre 3 (ou

recherche des vecteurs propres d'un espace vectoriel) :
Soit une matrice symtrique [M] qui fait correspondre deux vecteurs de l'espace
x,y,z
[V] = [M][U]
On montre en algbre linaire que si la matrice est symtrique d'ordre 3 :
1) Il existe toujours 3 vecteurs U orthogonaux tel que
[V] = [M][U] = [U] (autrement dit il existe 3 couples U et
r
V qui sont
colinaires 2 par 2).
2) Les 3 valeurs de (les valeurs propres de la matrices, autrement dit les
coefficients de proportionnalit entre les 3 couples U et V ) sont relles.
3) La matrices [M] est diagonale dans le repre construit sur les vecteurs
propres (les termes de la diagonale principale tant les valeurs propres).
Exemple : Diagonalisez la matrice symtrique [M] ci-dessous et donnez les 3
vecteurs propres (ordonnez les 3 valeurs propres par ordre dcroissant). Vrifiez
que les 3 vecteurs propres sont orthogonaux.
[M]xyz
Si [U] est un vecteur propre on a :
1 2 0
= 2 1
0
0
0 2
[V] = [M][U] = [U]
[[M] [I]][U] = [0]

IUTB-LYON1-GMP-DDS
Page 121
TD-2A-S3-F312
Pour que ce systme homogne (les seconds membres sont nuls) admette des
solutions diffrentes de 0 ( [U] [0 ] ) il faut que son dterminant soit nul.
[M] [I]
1
2
0
= 0
1
0
= 0
0
2
En dveloppant le dterminant par rapport sa troisime ligne on obtient :
(2 )( 3 )(1 + ) = 0
La rsolution de lquation donne les trois valeurs propres :
1 = 3
0
3 0
= 0 1 0
0 0 2
[M]XYZ
2 = 1
3 = 2
Recherche du premier vecteur propre U1 relatif la valeur propre 1 = 3 :
[[M] [I]][U ] = [0]

1
0
2 2 0 x

2 2 0 y = 0
0
0
0 5 z
x + y = 0
[U ]
1
z = 0
1
1
=
1
2
0
Recherche du deuxime vecteur propre U2 relatif la valeur propre
2 = 1 :
[[M] [I]][U ] = [0]

2
0
2 2 0 x

0 y = 0
2 2
0
0
0 1 z
x y = 0
[U ]
2
z = 0
r
1
1
=
1
2
0
Recherche du troisime vecteur propre U3 relatif la valeur propre 3 = 2 :
[[M] [I]][U ] = [0]

3
IUTB-LYON1-GMP-DDS
Page 122
TD-2A-S3-F312
0
3 2 0 x

2 3 0 y = 0
0
0
0 0 z
3x 2 y = 0
[U ]
3
0z = 0
0

= 0
1
On peut facilement vrifier que les 3 vecteurs propres sont orthogonaux.
2.4. Transformation dun vecteur et d'une matrice lors de la rotation d'un

repre dans le plan :
1) Transformation d'un vecteur :
[U']= t [R][U]
y
cos sin 0
[R] = sin cos 0
0
0
1
Uy
U
Ux
[R] est la matrice de rotation dun
Uy
vecteur dans le plan (voir rotation

d'un vecteur dans le plan).
Ux
2) Transformation d'une matrice :
On montre que : [M']= [R ][M ][R ]

t
2.5. Forme matricielle du produit scalaire :

1
[V] = 2
( 3,1)
3

4
[W] = 5
( 3,1)
6

r
r
V W = 1 4 + 2 5 + 3 6 = 32 est un scalaire, autrement dit une
"matrice" d'ordre 1.
r
r t
V W = [V] [W ]
( 1,1)
(1,3 ) (3,1)
ou
[W][V]
4
r
r

V W = [1 2 3] 5 = 1 4 + 2 5 + 3 6 = 3
6
IUTB-LYON1-GMP-DDS
Page 123

1484

Transféré par

Informations du documentcliquez pour développer les informations du document

Droits d'auteur :

Formats disponibles

1484

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

1484

Transféré par

Droits d'auteur :

Formats disponibles

SEMESTRE 3

avec la contrainte normale =0

Les axes principaux sont positionns par tan =

La direction propre X est calcule par tan =

En utilisant la relation de CAUCHY :

Compte tenu des relations : hh = [h ][ ][h ] = [cos 30 sin 30 0] 36 28 0 sin 30

do : hh = 7 cos 2 30 + 28 sin 2 30 + 2.36 sin 30 cos 30

th = (7 28)sin 30 cos 30 + 36(cos 2 30 sin 2 30)

qui mnent au mme rsultat.

Rsum des diffrents rsultats obtenus.

Reprendre le problme avec le cercle de MOHR :

Tracer les axes hh et th et adopter une chelle 1cm=10MPa

Soit un arbre section circulaire de

Isolons le point P de la section G soumise

La facette z est perpendiculaire au bord libre xy.

Compte tenu du module F213 :

Cours sur la torsion : P max i =

La matrice des contraintes au point P est donc :

Elments Principaux des contraintes :

La direction z est principale et Z=0.

La direction propre X est calcule par

Le cisaillement maximum en module : thmax i = Z =

Considrons l'tat plan de contraintes au

Matrice des contraintes en P dans le repre xyz :

Elments principaux des contraintes :

600 0 MPa ; [ ]XYZ = 0

Un arbre de diamtre D est soumis un

La formulation de dpart est identique celle du Problme 3 :

Compte tenu du module F213 :

Cours sur la flexion :

Cours sur la torsion :

+ M Y P 117 ,81 .10 3 64 D 1, 2 6

La matrice des contraintes au point P est donc :

La contrainte quivalente de VON MISES pour une poutre :

Un arbre de diamtre D = 20 mm est soumis suivant son axe un effort de traction de

La contrainte de traction due leffort normal :

La contrainte de torsion due au couple :

La contrainte quivalente de VON MISES pour une poutre :

La contrainte normale est donc la superposition

En ce plaant la limite le coefficient de scurit vaut 1.

Nous pouvons en dduire le moment de flexion :

Une barre section circulaire de rayon r est

Compte tenu du module F213 :

Cours sur la flexion : xxP =

Cours sur la torsion : P max i =

Si le matriau travaille la limite, le coefficient de scurit vaut 1. Appliquons le critre de

Nous en tirons les expressions de a et b en fonction de r et e :

Cest lquation dune ellipse de grand et petit axe a et b:

Effectuons une coupe suivant un plan quelconque

deux parties rsultant de la coupe. Soit h la

La partie isole est en quilibre sous l'action d'une partie

du torseur extrieur et des forces intrieures f

Par dfinition, le vecteur contrainte au point P relatif au plan de coupe dont la

normale extrieure (au solide) est h , est dfini par :

1.2. Composantes du vecteur contrainte :

si la normale h est quelconque par rapport

aux axes du repre.