Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Dimensionnement Des Structures

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 87

Dimensionnement des

structures
Antoine Legay
Matre de confrence
2015-2016
Cnam-Paris

Table des matires

I Poutre et torseur de cohsion

I.1 Introduction au dimensionnement des structures 1


I.1.1 Modle mcanique 2

I.2 Modle de poutre 3


I.3 Poutre dans son environnement 5
I.4 Torseur de cohsion 6
I.4.1 Dfinition 6
I.4.2 Dtermination 7
I.4.3 Classification des sollicitations 8

II Sollicitations simples sur les poutres

II.1 Traction 9
II.1.1
II.1.2
II.1.3
II.1.4
II.1.5

Torseur de cohsion 9
Contrainte normale 10
Allongement, dformation et dplacement 10
Relation contrainte-dformation 12
Relation entre effort normal et chargement 12

II.2 Torsion 13
II.2.1
II.2.2
II.2.3
II.2.4
II.2.5
II.2.6

Torseur de cohsion 13
Moment quadratique polaire de section 14
Contrainte tangentielle 14
Dformation et rotation des sections 15
Relation contrainte-dformation 16
Relation entre moment de torsion et chargement 16

II.3 Flexion 17
II.3.1 Torseur de cohsion 17
II.3.2 Moment quadratique de section 18

ii
II.3.3
II.3.4
II.3.5
II.3.6
II.3.7

Contrainte normale 18
Dformation 21
Dplacement 21
Relation contrainte-dformation 21
Relations moment de flexion - effort tranchant - chargement 22

III Calcul de treillis

25

35

III.1 Hypothses et critre de dimensionnement 25


III.1.1 Hypothses sur les liaisons 25
III.1.2 Rgles de construction dun treillis 26
III.1.3 Critre de dimensionnement 28

III.2 Mthode des nuds 28


III.3 Flambage des poutres droites 28
III.3.1 Introduction 28
III.3.2 Charge critique de flambage dune poutre droite 29
III.3.3 lancement et rayon de giration 31
III.3.4 Critre de dimensionnement 32
III.3.5 Autres conditions aux limites 33

IV Contraintes et dformations

IV.1 Introduction 35
IV.2 Caractrisation des contraintes et des dformations tridimensionnelles 36
IV.2.1 Oprateur des contraintes et des dformations 38
IV.2.2 Thorme de superposition 39

IV.3 Problme plan 39


IV.3.1
IV.3.2
IV.3.3
IV.3.4
IV.3.5
IV.3.6
IV.3.7

Hypothses 39
Etat de contraintes planes 40
Expressions des contraintes subies par un carr non align avec x et y 41
Expressions des dformations dun carr non align avec x et y 44
Relation entre les contraintes et les dformations dun carr non align avec x et y
Directions principales 45
Cercle de Mohr des contraintes 45

V Critres de dimensionnement
V.1
V.2
V.3
V.4
V.5

44

49

Objectifs 49
Matriaux ductiles : critre de Tresca 49
Matriaux ductiles : critre de Von Mises 51
Comparaison des critres de Tresca et de Von Mises 52
Fatigue des matriaux 52

VI Enveloppes minces

55

59

VI.1 Action dun fluide au repos sur un solide 55


VI.2 Application un rservoir cylindrique 56

VII Initiation au calcul lments finis

VII.1 tude de llment de barre 59


VII.1.1 quilibre de llment barre 59
VII.1.2 Exemple dapplication 60
VII.1.3 Remarques sur la mthode des lments finis 61

iii

VII.2 tude de deux barres 61


VII.2.1 Assemblage des matrices de rigidit lmentaires 61
VII.2.2 Mise en uvre pratique 63

VII.3 lment barre pour le calcul des treillis 64


VII.4 Elment de poutre pour le calcul des portiques 65

VIII Moyens exprimentaux

VIII.1 Jauges de dformation 67


VIII.1.1
VIII.1.2
VIII.1.3
VIII.1.4
VIII.1.5
VIII.1.6

Principe 67
Pont de Wheatstone 68
Utilisation du boitier 70
Diffrents montages 71
Capteurs jauges 72
Exploitation dune rosette de 3 jauges 45o

VIII.2 Photolasticit 74
VIII.2.1 Principes 74
VIII.2.2 Mise en quation 77
VIII.2.3 Rseaux de courbes caractristiques

79

73

67

I Poutre et torseur de cohsion

I.1 Introduction au dimensionnement des structures


Dimensionnement des structures
Une structure est un assemblage intelligent dlments et de matriaux afin dassurer une
fonction. La figure I.1 montre par exemple la structure en balsa dun avion daromodlisme
permettant dassurer la forme de la voilure portante, ainsi que la structure dun pylne
lectrique qui permet de maintenir les lignes lectriques une certaine hauteur.
Le but du dimensionnement est de dterminer les formes, dimensions, matriaux afin de
satisfaire la fonction demande dans toutes les conditions de vie de la structure. Par exemple
la structure en balsa de lavion daromodlisme doit rsister aux efforts arodynamiques

Figure I.1 Exemples de structures : structure en balsa dun avion daromodlisme, pylne
lectrique

Poutre et torseur de cohsion

Figure I.2 Problme rel : dimensionnement des pieds dune table.

en vol, la structure du pylne lectrique doit rsister des vents forts et des surcharges de
neige et de verglas.
Deux principales mthodes existent pour dimensionner une structure :
Mthode non prdictive "essai-erreur" : on construit un prototype rel (ou une maquette
chelle rduite), puis on le teste en condition relle ; cette mthode a lavantage de
ne faire appel aucune connaissance a priori de la mcanique mais est coteuse.
Mthode prdictive : on fait un modle mcanique "virtuel" bas sur des quations mathmatiques, puis on le teste ; cette mthode est moins couteuse, mais a linconvnient
de faire appel des connaissances de mcanique et de mathmatiques.
Cest cette deuxime mthode qui est dveloppe dans ce cours. On se limite au dimensionnement des structures en statique et en lasticit linaire.

Problme rel
Le problme rel fait intervenir (Fig. I.2) :
Une structure, comprenant des incertitudes sur sa gomtrie et son matriau ;
Des liaisons avec lextrieur, souvent assez mal maitrises ;
Des efforts appliqus, parfois assez complexes.
Lors de la phase de conception, la solution relle de ce problme nest pas accessible (dplacements, contraintes, ...). Une fois la structure fabrique et place dans son environnement,
la solution est partiellement accessible par des mesures (jauges de dformation, photolasticit,... ).

I.1.1 Modle mcanique


Afin de trouver une solution approche du problme rel, on utilise un modle mathmatique
du problme rel. Les modles gnralement utiliss en mcanique sont :
le modle de poutre,
le modle de plaque,

I.2 Modle de poutre


1111
0000
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111

1 variable

a) poutre

2 variables

b) coque

111
000
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111

3 variables

c) tridimensionnel

Figure I.3 Trois modles du pied de table.

le modle de coque,
le modle plan en contraintes planes,
le modle plan en dformations planes,
le modle axisymtrique,
le modle tri-dimensionnel.
Pour lexemple prcdent dun pied de table, on peut par exemple choisir :
Le modle de poutre (Fig. I.3 a) :
hypothse cinmatique de poutre
1 variable le long de laxe de la poutre dcrit le problme
encastrement de type poutre
torseurs defforts quivalents
Le modle de coque (Fig. I.3 b) :
hypothse cinmatique de coque
2 variables sur la surface moyenne de la coque dcrivent le problme
encastrement de type coque
torseurs defforts quivalents distribus
Le modle tri-dimensionnel (Fig. I.3 c) :
encastrement tri-dimensionnel
3 variables dans les 3 directions de lespace dcrivent le problme
forces surfaciques distribues
Pour les trois modles proposs, lencastrement est modlis de faon parfaite alors que la
liaison relle est ralise par une pice intermdiaire souple. Ces modles ne permettent pas
de dimensionner cette pice intermdiaire. Cest au concepteur de choisir le modle le plus
adapt par rapport aux critres de dimensionnement quil pense tre les plus judicieux.

I.2 Modle de poutre


Hypothses gomtriques
La figure I.4 montre un assemblage de poutres permettant de construire une charpente
mtallique. Une poutre est un solide dont une dimension est plus grande que les 2 autres.

Poutre et torseur de cohsion

Figure I.4 Charpente constitue dun assemblage de poutres

11
00
00
11
00
11
00
11

MO

s
O
Figure I.5 Modle de poutre et coupure fictive en deux parties.

Dun point de vue plus gomtrique, une poutre est un solide engendr par une surface
plane appele section droite constante ou lgrement variable dont le plan reste orthogonal
une courbe de grand rayon appele ligne moyenne dcrite par le centre de surface O
de la section droite (Fig. I.5). La plus grande dimension transversale est petite devant la
longueur de la fibre moyenne (rapport de 5 10 au moins).
Dans le cadre de ce cours, on ne sintresse quaux poutres droites, cest dire celles
dont la ligne moyenne est une droite. De plus, dans la majorit des cas, on ne prend que
des sections constantes.

Hypothses sur le matriau


On suppose que le matriau est :
homogne : les proprits sont les mmes en tout point,
isotrope : les proprits sont les mmes dans toutes les directions, ce qui nest pas le
cas dun matriau composite,
continu : le matriau ne contient pas dasprits, il y a continuit de la matire.

Hypothses sur les dformations et les dplacements


On suppose que les dformations sont petites et restent dans le domaine lastique. De plus
le dplacement est aussi considr petit devant la taille de la poutre, cela permet de faire

I.3 Poutre dans son environnement

F
2

F
2

Figure I.6 Principe de Saint-Venant.

Figure I.7 Hypothse de Bernoulli.

tous les calculs sur la configuration initiale et non dforme de la poutre.

Hypothses sur le chargement


Le chargement est appliqu progressivement, on nglige les effets dinertie et on ne sintresse
qu la configuration finale statique.

Principe de Saint-Venant
Les contraintes et les dformations dans une rgion loigne du point dapplication des
efforts ne dpendent que du torseur des efforts de cohsion au point considr (Fig. I.6).
Autrement dit, la faon dont on applique le chargement na pas dinfluence loin de
lapplication de la charge.

Hypothse de Bernoulli
Une section de la poutre initialement plane et perpendiculaire la ligne moyenne reste plane
et perpendiculaire la ligne moyenne aprs dformation (Fig. I.7).

I.3 Poutre dans son environnement


La poutre dans son environnement subit des actions mcaniques extrieures (Fig. I.8). Ces
actions sont partages en deux groupes : actions des liaisons et actions des efforts extrieurs.
La poutre est en liaison avec lextrieur aux points Ai :

y
A1

x
B1

B2

An Bm

Figure I.8 Poutre dans son environnement.

Poutre et torseur de cohsion


n

torseurs des actions mcaniques nots LAi


ces actions sont inconnues

Les efforts extrieurs sont appliqus aux points Bi


n

torseurs des actions mcaniques nots Bi


ces actions sont connues

Dans le plan (
x ,
y ), les liaisons classiquement rencontres
sont
:

Fx M x

n o
encastrement :

L
=
0
F
x
y

yA

articulation :

appui simple :

x nLo =

yA

x n o

L =

Fx

M z A

Mx

0 A

Mx

0 A

La poutre est en quilibre sous laction des efforts cits ci-dessus. Le principe fondamental
de la statique scrit
i=nb.X
liaisons n
i=1

LAi +

i=nb.
efforts n
X

n o

Bi = 0

i=1

ou plus simplement avec les notations de ce cours


n o

n o

n o

L + = 0 .

La poutre est en quilibre isostatique si toutes les inconnues de liaisons peuvent tre
dtermines en crivant lquilibre de la poutre. La poutre est en quilibre hyperstatique si
certaines inconnues de liaisons ne peuvent pas tre dtermines en crivant lquilibre de la
poutre : le nombre dinconnues non dtermines donne le degr dhyperstatisme. Enfin, si il
y a plus dquations que dinconnues de liaisons, la poutre est mobile : le nombre dquations
supplmentaires par rapport au nombre dinconnues de liaisons est le degr de mobilit.

I.4 Torseur de cohsion


I.4.1 Dfinition
Une poutre est coupe en deux parties par une section fictive de centre de section
O (Fig. I.5). Le torseur de cohsion est le torseur des actions mcaniques de 2 sur 1
travers la surface exprim en O .
Laction surfacique de 2 sur 1 travers la surface est une force surfacique appele

"vecteur contrainte" pour la direction normale


n ; cette action surfacique est note

T (M, n ) (Fig. I.9). Le torseur de cohsion est la somme sur la surface de T (M,
n ) o

n est le vecteur tangent laxe de la poutre en O .

I.4 Torseur de cohsion

T (M,
n)

O
Figure I.9 Vecteur contrainte.

La rsultante vaut :

ZZ

s =

T (M,
n )d.

Le moment au point O de vaut :

MO =

ZZ

O M T (M,
n )d.

Finalement le torseur de cohsion en O vaut :


n

KO =

MO O

I.4.2 Dtermination
La dtermination du torseur de cohsion se fait en crivant lquilibre de 1 ou de 2 . On
utilise les notations suivantes :


L : torseur des actions mcaniques inconnues sur la poutre provenant des liaisons,

: torseur des actions mcaniques connues sur la poutre ,

L 1 : torseur des actions mcaniques inconnues sur la partie 1 ,

L 2 : torseur des actions mcaniques inconnues sur la partie 2 ,

2 1 = KO : torseur des actions mcaniques de 2 sur 1 , cest dire le torseur




1 : torseur des actions mcaniques connues sur la partie 1 ,


2 : torseur des actions mcaniques connues sur la partie 2 ,


de cohsion.

Lquilibre de la poutre scrit

n o

n o

n o

L + = 0 .
n o

Ceci permet de dterminer le torseur inconnu L dans le cas ou le problme est isostatique.
Lquilibre de la partie 1 scrit
n

n o

L1 + 1 + 21 = 0

o 2 1 est le torseur de cohsion KO . On en dduit que le torseur de cohsion vaut


n

KO = L 1 1 .

Poutre et torseur de cohsion

Mt

Traction

Torsion

Mfz

Ty x
Flexion

Figure I.10 Diffrentes sollicitations pour les poutres.


Traction

o
n
KO =

N
0

Torsion

0
0

0 O

Mt
0

0
O

0
0

Flexion pure

(plan (
x ,
y ))

0
0

Mfz O

0
0

Tableau I.1 Sollicitations simples.

Flexion simple

(plan (
x,
y ))

0
T
y
0

0
0

Mfz O

Lquilibre de la partie 2 scrit


n

n o

L2 + 2 + 12 = 0

o 1 2 = 2 1 = KO

dduit que le torseur de cohsion vaut

daprs le thorme des actions mutuelles. On en

KO = L 2 + 2 .

I.4.3 Classification des sollicitations

Dans le cas dune poutre droite, on utilise habituellement le repre orthonorm (


x ,
y ,
z)

o x est suivant la ligne moyenne et les vecteurs y et z sont dans les directions transversales. Dans ce repre, le torseur de cohsion vaut
n

KO =

T
z

Mt

M fy

M f z O

N est leffort normal,

Ty est leffort tranchant suivant


y,

T est leffort tranchant suivant


z,
z

Mt est le moment de torsion,

Mfy est le moment de flexion autour de


y (flexion dans le plan (
x ,
z )),

Mfz est le moment de flexion autour de


z (flexion dans le plan (
x ,
y )).
La figure I.10 prsente graphiquement ces diffrents cas. Les sollicitations simples tudies dans la suite du cours sont rpertories dans le tableau I.1.

II Sollicitations simples sur les poutres

II.1 Traction
II.1.1 Torseur de cohsion

On suppose une poutre daxe


x soumise un chargement de traction : F
x lextremit

droite en A et F x lextrmit gauche en O (Fig. II.1). Laire de la section est note S.


En isolant la partie (2) pour calculer le torseur de cohsion,
n

KO = L 2 + 2

F
x y

l
A

F
x

Figure II.1 Poutre en traction.

10

Sollicitations simples sur les poutres

avec
n

et

n o

L2 = 0

o
o n
n

2 = F
x|0
= F
x | O A F
x

{z

= 0 car colinaires

on a

= F
x|0

o n

x|0
KO = F

donc leffort normal N vaut ici F .

II.1.2 Contrainte normale


On suppose que le vecteur contrainte dans une section de la poutre (perpendiculaire la

ligne moyenne) est port par


x et quil est uniforme sur toute la surface . On note le

vecteur contrainte x ou est la contrainte normale de traction.


Le torseur de cohsion en O centre de la section vaut :
n

avec

S =

ZZ

TO =

)
S

MO O
ZZ

x d =

x = S
x

| {z }

=S aire section

et

MO =

ZZ

(y
y + z
z )
x d = 0

car y et z sont des fonctions impaires intgres sur des intervalles symtriques. Finalement :
n

TO =

N
x

o N est leffort normal qui vaut N = S ou encore


=

N
.
S

On parle de traction lorsque N et sont positifs, on parle de compression quand N et


sont ngatifs.

II.1.3 Allongement, dformation et dplacement


Allongement
Lorsque la poutre est soumise de la traction, le matriau tant dformable, elle sallonge.
Cet allongement, not l, a une unit (millimtre ou mtre). La valeur de l, pour un
mme matriau et une mme section S, dpend de la longueur l.

11

II.1 Traction

x + x

x
u(x)

u(x + x)

u(x + x) u(x) + x
Figure II.2 Dformation dun tronon de poutre en traction

Dformation
Afin de pouvoir comparer les rsultats entre eux, on prfre utiliser une grandeur adimensionne que lon note et qui vaut
=

l
.
l

On appelle lallongement axial unitaire.


La dformation axiale entraine pour la plus part des matriaux des dformations dans

les directions transversales


y et
z notes y et z . Elles valent
y = z = .
Le coefficient est appell coefficient de Poisson, il caractrise le rtrcissement transversal
dune poutre sollicite en traction. Ce coefficient a t caractris par Simon Denis Poisson
(franais, 1781-1840).
Dplacement

La section fictive situe labscisse x (OO = x


x ) se dplace dans la direction axiale
x

dun vecteur u(x) x . Autrement dit, le vecteur u(x) x est le vecteur dplacement du point
O .
On montre quil existe une relation entre la dformation (x) labscisse x et le dplacement u(x) donne par
du(x)
.
dx
La dformation est la drive du dplacement par rapport x.
(x) =

Preuve On isole un petit tronon de poutre de longueur x compris entre les abscisses x et x + x (Fig.
II.2). Lextrmit de gauche se dplace dune valeur u(x), la section de droite se dplace dune valeur
u(x + x). La longueur finale du tronon est u(x + x) u(x) + x, son allongement est u(x + x) u(x).
La dformation du tronon est alors
u(x + x) u(x)
=
x
En faisant tendre x vers 0, on obtient exactement la dfinition de la drive :
lim

x0

u(x + x) u(x)
du(x)
=
x
dx

On a bien
=

du(x)
dx

12

Sollicitations simples sur les poutres

F
L

Re

rupture
zone plastique
zone lastique

Figure II.3 Essai de traction

II.1.4 Relation contrainte-dformation


Les essais sur des prouvettes soumises de la traction (Fig. II.3) montrent que et sont
proportionnels et obissent la loi
= E
dans le domaine lastique du matriau (voir aussi la figure V.1 dans la partie sur le critre
de Tresca). Cette constante E ne dpend que du matriau, cest le module de Young ou
module dlasticit longitudinal. Ce module a t caractris par Thomas Young (anglais,
1773-1829).
Son unit est le Pascal.
La limite lastique Re est la contrainte maximale admissible du domaine lastique. Le
critre de dimensionnement dune poutre en traction est donc
<

Re
s

o s est le coefficient de scurit adopt.


En remplaant par son expression en fonction de u, puis de N , on a
=

N
=
,
E
ES

soit la relation entre le dplacement et leffort normal


du(x)
N (x)
=
.
dx
ES(x)

II.1.5 Relation entre effort normal et chargement

On tudie un tronon de poutre qui nest soumis qu un effort axial linique p


x . Par

exemple, pour une poutre verticale o laxe x est vers le bas, cette force linique vaut
p = gS
o est la densit du matriau composant la poutre, g est la gravit et S est laire de la
section de la poutre.
Dans ce tronon, on montre que lon a la relation suivante entre leffort normal et la
charge linique :

dN
+ p = 0.
dx

13

II.2 Torsion

x + x

N(x)

N(x + x)

p
x

Figure II.4 quilibre dun tronon de poutre en traction

111
000
000
111
000
111
000
111
000
111
000
111
000
111

er

M
r

C
x

Figure II.5 Poutre sollicite en torsion.

Preuve On isole un petit tronon de poutre de longueur x compris entre les abscisses x et x + x (Fig.
II.4). Lextrmit de gauche subit la force N (x) car la matire est droite et, par convention, leffort
normal est laction de la partie de droite sur la partie de gauche (en supposant que laxe positif soit vers la
droite). Lextrmit de droite subit la force N (x+x). Le tronon subit la force linique p, soit la rsultante
px. Lquilibre du tronon scrit
N (x) + N (x + x) + px = 0
soit en divisant par x

N (x + x) N (x)
+p = 0
x
En faisant tendre x vers 0, on obtient exactement la dfinition de la drive :
lim

x0

dN (x)
N (x + x) N (x)
=
x
dx

On a bien

dN (x)
+p=0
dx

II.2 Torsion
II.2.1 Torseur de cohsion

On suppose une poutre de section circulaire daxe


x , encastre son extrmit gauche O

et soumise un moment C x son extrmit droite A (figure II.5).


En isolant la partie (2) pour calculer le torseur de cohsion,
n

KO = L 2 + 2

avec
n

n o

L2 = 0

14

Sollicitations simples sur les poutres

er y

Figure II.6 Rpartition de la contrainte tangentielle dans la section dune poutre sollicite en
torsion.

et

2 =

on a

0 | C
x

0 | C
x + O A 0

KO =

donc le moment de torsion Mt vaut ici C.

0 | C
x

0 | C
x

II.2.2 Moment quadratique polaire de section


Le moment quadratique polaire de la section dabscisse x est not I0 (x). Il peut dpendre
de x si la section de la poutre varie, son expression est :
I0 (x) =

Z Z

r2 dS

S(x)

Pour une section circulaire de rayon R (diamtre D = 2R), le calcul est le suivant :
I0 =

r=R Z =2

r=0

r2 (r d dr) = 2

=0

D4
R4
=
4
32

Pour une section circulaire creuse de diamtre extrieur De et de diamtre intrieur Di ,


le moment quadratique polaire vaut
I0 =

II.2.3 Contrainte tangentielle

De4 Di4 .
32

On utilise le repre cylindrique (


er ,
e ,
x ) pour reprer un point M de la poutre, ainsi

OM = x x + r e r . On suppose que le vecteur contrainte au point M appartenant une

section de la poutre et pour la direction axiale


x est port par
e et quil est proportionel

r (Fig. II.6) :

T (M,
x ) =
e

On peut montrer que la contrainte tangentielle en un point M de la poutre vaut :


(x, r) =

Mt (x)
r
I0 (x)

o I0 (x) est le moment quadratique polaire en O de la section.


Preuve On suppose que la rpartition de contrainte tangentielle est proportionnelle r :
(x, r) = max

r
R

15

II.2 Torsion

(x + x)

(x)

(x)

(x)

x + x

(x + x) (x)

x
Figure II.7 Dformation et rotation des sections
o max est la valeur de la contrainte sur la peau de la poutre et R est la rayon de la poutre. Le moment

du torseur de cohsion MO vaut

MO =

ZZ

O M T (M,
n ) d

soit ici

MO =

ZZ

Donc on a

r
er max

r
max

e d =
R
R
Mt =

soit aussi

ZZ

max

I0
x
r 2 d
x =
R

I0
max
R

R
Mt
I0
En remplaant dans lexpression de , on a bien le rsultat attendu.
max =

II.2.4 Dformation et rotation des sections


Une section dabscisse x tourne dun angle (x) sans se dformer. La dformation est une
distorsion angulaire entre deux sections trs proches caractrise par (x). Elle est proportionnelle r et la variation de (x), soit
=r

d
.
dx

Preuve Deux sections distantes de x tournent respectivement des angles (x) et (x + x). La rotation
relative entre les deux sections est de (x + x) (x) (Fig. II.7). Cette rotation relative provoque la
distorsion angulaire (x) qui caractrise la dformation en torsion. La relation entre (x) et (x) est
tablie gomtriquement en dessinant deux triangles rectangles :
un triangle dans le plan tangent au cylindre, qui fait apparaitre langle (x) ( gauche sur la figure
II.7),
un triangle dans le plan de la section, qui fait apparaitre langle (x + x) (x) ( droite sur la
figure II.7).

Ces deux triangles partagent un cot de longueur a. On peut crire dans le triangle faisant intervenir (x)
que
a
tan (x) =
x

16

Sollicitations simples sur les poutres

et dans le triangle faisant intervenir (x) que


tan((x + x) (x)) =

a
r

Les angles tant petits, on a


tan (x) (x)
et
tan((x + x) (x)) (x + x) (x)
En crivant lgalit de a, on a
(x) x = r((x + x) (x))
En divisant par x et en faisant tendre x vers 0, on a la dfinition de la drive
(x) = r lim

x0

(x + x) (x)
x

On a bien

d
.
dx

(x) = r

II.2.5 Relation contrainte-dformation


La contrainte tangentielle est proportionnelle la distorsion angulaire et au module
dlasticit transversal G :
= G.
En remplaant par son expression en fonction de , on a
= rG

d
.
dx

En utilisant lexpression de en fonction du moment de torsion, on a donc


Mt (x)
d(x)
=
.
dx
GI0 (x)
Le module dlasticit transversal G est fonction du module dYoung E et du coefficient
de Poisson :
G=

E
2(1 + )

II.2.6 Relation entre moment de torsion et chargement


On tudie un tronon de poutre qui nest soumis qu une rpartition linique de moment

c
x.
Dans ce tronon, on peut montrer que lon a la relation suivante entre le moment de
torsion et la rpartition linique de moment :
dMt
+ c = 0.
dx

Preuve Soit un petit tronon de poutre de longueur x (Fig. II.8). Lquilibre en moment autour de
x
scrit :
Z
x+x

Mt (x) + Mt (x + x) +

c(x) dx = 0

En supposant que le couple linique c(x) est constant on a :

x+x

c(x) dx = c x

17

II.3 Flexion

c
x

Mt (x)
x

Mt (x + x)
x

x + x

x
Figure II.8 quilibre dun tronon de poutre en torsion

1111
0000
0000
1111
0000
1111
0000
1111
0000
1111

1111
0000
0000
1111
0000
1111
0000
1111
0000
1111

A
x

F
y

allongement des fibres, traction

x
v(x)

Dforme de la section
(coefficient de Poisson)

rtrcissement des fibres, compression


Figure II.9 Poutre sollicite en flexion.
donc lquilibre devient
Mt (x) + Mt (x + x) + c x = 0

Mt (x + x) Mt (x)
= c
x

En faisant tendre x vers 0, on a :


dMt
= c
dx

soit aussi

dMt
+ c = 0.
dx

II.3 Flexion
II.3.1 Torseur de cohsion

On suppose une poutre daxe


x , encastre son extrmit gauche O et soumise un

chargement transversal F y son extrmit droite A (Fig. II.9).


En isolant la partie (2) pour calculer le torseur de cohsion, on a :
n

KO = L 2 + 2

avec

n o

L2 = 0

18

Sollicitations simples sur les poutres

et
n
o
o n

= F
y | O A(F
y)
2 = F
y | 0

2 =

on a

KO =

= F
y | (lx)
x (F
y)
o

F
y |(x l)F
z

F
y |(x l)F
z

donc leffort tranchant Ty et le moment de flexion Mfz valent ici :


Ty = F

Mfz = (x l)F.

II.3.2 Moment quadratique de section

Le moment quadratique autour de laxe (O ,


z ) de la section dabscisse x est not I(O ,
z ) (x).
Il peut dpendre de x si la section de la poutre varie, son expression est :

I(O ,
z ) (x)

Z Z

y 2 dS
S(x)

On dfinit le moment quadratique autour de laxe (O ,


y ) comme :

I(O ,
y ) (x) =

Z Z

z 2 dS.

S(x)

Pour une section rectangulaire dpaisseur h suivant


y et de largeur b suivant
z , le
calcul est le suivant :

I(O ,
z ) (x) =

h
2

h2

b
2

2b

y 2 dzdy = b

h
2

h2

y 2 dy = b

h y 3 i h2

h2

bh3
.
12

De mme, on a :

hb3
.
12
Le thorme de transport de Huygens permet de calculer le moment quadratique dune

I(O ,
y ) (x) =

section plus complexe :


I = I + Sd2
o I est le moment quadratique autour dun axe passant par le barycentre de la section, S
est laire de la section et d est la distance entre les deux axes et . Le calcul de I pour
la section en "i" scrit :

 bh3
b(a h)3
a 2 
+2
+ bh
12
12
2
Pour les sections les plus courantes, les moments quadratiques sont donns sur la figure

I=

II.10.

II.3.3 Contrainte normale


Une poutre en flexion subit des contraintes normales la section et des contraintes tangentielles . Dans le cas tudi, les fibres suprieures sont tires et les fibres infrieures
sont comprimes. La fibre moyenne ne subit pas de contraintes normales.

19

II.3 Flexion

D
b

I=

a4
12

I=

bh3
12

I=

y
h

d
D
(D4 d4 )
I=
64

b
section en I

Figure II.10 Moments quadratiques des sections courantes

D4
64

20

Sollicitations simples sur les poutres

max

>0

M(x, y)
y

x
<0
max

Figure II.11 Contrainte normale dans une section de poutre sollicite en flexion

On peut montrer que la contrainte normale en un point M (x, y) de la poutre vaut :


(x, y) =

Mfz (x)
y

I(O ,
z ) (x)

o I(x) est le moment quadratique de la section.


De plus, mme si la contrainte tangentielle nest pas uniforme dans une section, une
bonne approximation de sa valeur est
(x) =

T (x)
.
S(x)

La contrainte tangentielle est gnralement trs petite devant la contrainte normale, elle est
donc souvent nglige.
Preuve On suppose que le vecteur contrainte dans une section de la poutre scrit sous la forme(Fig.
II.11) :

T (M,
n ) = max
x = (x, y)
x
h/2
o (x, y) est la contrainte normale due la flexion, La rsultante du torseur de cohsion vaut :

ZZ

max

T (M,
n ) d =
h/2

s =

ZZ

y d
x = 0

car y est une fonction impaire intgree sur un intervalle symtrique (de h/2 h/2).
Le moment du torseur de cohsion vaut :

MO =

ZZ

O M T (M,
n ) d =

max
MO =
h/2
Or

ZZ

zy d = 0

donc :

et

ZZ

ZZ


y

(y
y + z
z ) max
x d
h/2

z + zy
y ) d
(y 2

ZZ

y 2 d = I(O ,
z ) (x),

max

z
MO = I(O ,
z ) (x)
h/2

Par ailleurs, le moment du torseur de cohsion vaut MO = Mfz (x)


z , donc :
max =

h/2
Mfz (x)

I(O ,
z ) (x)

En remplaant dans lexpression de la contrainte normale, on a :


(x, y) =

Mfz (x)
y

I(O ,
z ) (x)

21

II.3 Flexion

(x)

(x + x)

x
y
x

Figure II.12 Rotation (x) des sections pour une poutre sollicite en flexion

II.3.4 Dformation
Pour une poutre en flexion, la dformation est due la variation de langle de rotation (x)
des sections (Fig. II.12). Ltat de contrainte tant quivalent de la traction-compression
dans une fibre de la poutre, la dformation est un allongement unitaire dans la direction
axiale. Son expression en fonction de la rotation (x) des sections est :
(x, y) =

d(x)
y
dx

Preuve Soit un petit tronon de poutre de longueur x (Fig. II.12). La section dabscisse x tourne dun
angle (x), la section dabscisse x + x tourne dun angle (x + x). La fibre neutre de ce tronon a pour
longueur x, la fibre situe la distance y de la fibre neutre a pour longueur x aprs dformation. En
considrant que les angles de rotation des section sont petits, cette longueur x a pour expression :
x = x + y(x) y(x + x)
soit aussi
x x = y (x + x) (x)

Par ailleurs, la dformation dans la direction axiale de la fibre situe la distance y de la fibre neutre
vaut :
x x
=
x
Soit aussi
(x + x) (x)
= y
x
En faisant tendre x vers 0, on a
d(x)
y
=
dx

II.3.5 Dplacement
La poutre se dforme sous laction du chargement. Les points de la ligne moyenne se d

placent suivant
y de la valeur v(x). Autrement dit, le vecteur v(x)
y est le vecteur dplacement du point O .

La drive premire de v(x) reprsente la rotation de section dabscisse x autour de


z :
(x) =

dv(x)
dx

II.3.6 Relation contrainte-dformation


Ltat de contrainte tant quivalent de la traction-compression, la relation entre la

contrainte normale et lallongement unitaire suivant


x est
= E

22

Sollicitations simples sur les poutres

Mf (x)

x + x

x
p

T (x + x)
Mf (x + x)

T (x)

Figure II.13 quilibre dun tronon de poutre en flexion : effort tranchant

o E est le module dYoung du matriau.


Or on sait que
(x, y) =
donc :

En simplifiant par y, on a :

Mfz (x)
y

I(O ,
z ) (x)

et =

d(x)
y
dx

d(x)
Mfz (x)
y = E
y

I(O ,
dx
z ) (x)
Mf (x)
d(x)
=
.

dx
EI(O ,
z ) (x)

ou bien en fonction de dplacement de la ligne moyenne


Mf (x)
d2 v(x)
=
.

dx2
EI(O ,
z ) (x)

II.3.7 Relations moment de flexion - effort tranchant - chargement

On tudie un tronon de poutre qui nest soumis qu un effort transversal linique p


y.

Par exemple, pour une poutre horizontale o laxe y est vers le haut (gravit vers le bas),
cette force linique vaut
p = gS
o est la densit du matriau composant la poutre, g est la gravit et S est laire de la
section de la poutre.
Dans ce tronon, on peut montrer que lon a la relation suivante entre leffort tranchant
et la charge linique :

dTy
+ p = 0.
dx
On peut aussi montrer que lon a la relation suivante entre le moment de flexion et leffort

tranchant :

dMfz
+ Ty = 0.
dx

Preuve On isole un petit tronon de poutre de longueur x compris entre les abscisses x et x + x (Fig.
II.13).

Lextrmit de gauche subit la force T (x)


y car la matire est droite et, par convention, leffort
tranchant est laction de la partie de droite sur la partie de gauche (en supposant que laxe positif soit

vers la droite). Lextrmit de droite subit la force T (x + x)


y . Le tronon subit la force linique p, soit

la rsultante px y . Lquilibre du tronon scrit en rsultante sur


y :
T (x) + T (x + x) + px = 0
soit en divisant par x

T (x + x) T (x)
+p = 0
x

23

II.3 Flexion
En faisant tendre x vers 0, on obtient exactement la definition de la drive :
lim

x0

dT (x)
T (x + x) T (x)
=
x
dx

On a bien

dT (x)
+p=0
dx
Lextrmit de gauche subit par ailleurs le moment de flexion Mf (x + x). Lextrmit de droite subit

le moment de flexion Mf (x). Lquilibre en moment autour de


z au centre de section situe en x + x

(point O ) scrit :

Mf (x + x)
z Mf (x)
z + O
O (T (x)
y)= 0

avec O O = x
x . Le calcul du produit vectoriel donne :

O O (T (x)
y ) = T (x)x
z

On a alors en projetant sur


z et en divisant par x on a :
Mf (x + x) Mf (x)
+ T (x) = 0
x
En faisant tendre x vers 0, on trouve alors :
dMf
+ T = 0.
dx

24

Sollicitations simples sur les poutres

Traction
n
o
KO =

N
0

Torsion

0
0

0 O

0
0

Flexion

Mt
0

0
O

0
T
y
0

0
0

Mfz
O

dplacement

translation u
x

rotation
x

translation v
y

dformation

allongement

distorsion

allongement

relation

dplacement-dformation

relation
contrainte-dformation

relation

du
dx

= E

N
S

=r

d
dx

= G

Mt
r
I0

d2 v
y
dx2

= E

Mfz
y
I

contrainte-efforts

relation

du
N
=
dx
ES

d
Mt
=
dx
GI0

d2 v
Mf
=
2
dx
EI

dN
+p=0
dx

dMt
+c=0
dx

dTy
+p=0
dx

dplacement-efforts

quations reliant
les efforts
gnraliss

dMfz
+ Ty = 0
dx

Tableau II.1 Rsum des formules pour les sollicitations simples

III Calcul de treillis

III.1 Hypothses et critre de dimensionnement


III.1.1 Hypothses sur les liaisons
Un treillis est une structure compose de barres rotules entre elles (Fig. III.1). Les barres
sont connectes entre elles par des nuds, centres des liaisons rotules. Les articulations sont
supposes parfaites. Cette simplification permet de rsoudre relativement facilement le problme. Mme si les liaisons ne sont pas rellement des rotules mais des liaisons boulonnes,
on peut considrer dans une premire approche la structure comme un treillis de barres
rotules. Cela permet de trouver une bonne approximation des efforts normaux dans les
barres et donc de dimensionner.
La barre AB est connecte aux nuds A et B. Le torseur des actions mcaniques de A
sur la barre vaut en A
n

TA

o
= FA | 0

FA
A
barre AB

B
FB

Figure III.1 Structure de type treillis de barres

26

Calcul de treillis

De mme, le torseur des actions mcanique de B sur la barre vaut en B


n

TB

o
= FB | 0

Lquilibre de la barre scrit en A


n
n

o

o

o
= 0 | 0
FA | 0
+ FB | AB FB

Soient

FA = FB
et

AB FB = 0

Leffort extrieur FB est donc port par laxe de la barre. En notant i le vecteur unitaire
allant de A vers B, on a

i =

AB

k AB k

et

FB = FA = NAB i
o NAB est leffort normal dans la barre AB.
Une barre rotule ses deux extrmits ne subit que de la traction-compression. Chaque
barre reste un segment de droite aprs dformation, la figure III.1 montre la dforme globale
dun treillis (barres en tirets).
Si lon veut considrer les liaisons comme des encastrements, des moments sont alors
transmis et de la flexion apparat dans les barres ainsi que de la torsion dans le cas dune
structure tridimensionnelle : cest un portique. La modlisation dune structure comme un
treillis rotul a longtemps tait utilise dans les bureaux dtudes avant larrive des moyens
de calculs informatiques.

III.1.2 Rgles de construction dun treillis


Treillis plan isostatique
Le systme le plus simple est constitu par un triangle, soient 3 barres et 3 nuds (Fig.
III.2a). En notant n le nombre de nuds et b le nombre de barres, partir dun triangle
(n = 3, b = 4), chaque ajout de x nuds impose lajout de 2x barres, soit
n=3+x ;

b = 3 + 2x ;

2n = 3 + b

Le nombre de barres b est li au nombre de nuds n pour que le systme soit isostatique. Il
faut de plus que les encastrements du treillis nimposent pas dhyperstatisme la structure.

27

III.1 Hypothses et critre de dimensionnement

a) n=3 ; b=3
b) n=4 ; b=5
Treillis isostatiques

Treillis hyperstatiques

Figure III.2 Construction dun treillis

Figure III.3 Treillis avec mobilit interne

Treillis plan hyperstatique


Si 2n < 3 + b alors le treillis est hyperstatique, la mthode de rsolution des efforts normaux
prsente ensuite ne suffit pas dterminer elle seule les efforts dans les barres (Fig.
III.2b)). Il faut rsoudre le problme en crivant que les allongements des barres ne sont pas
indpendants pour que les barres restent articules entre elles.
Treillis plan avec mobilit interne
Si 2n > 3 + b alors il y a indtermination des efforts normaux et des mouvements sont
possibles sans efforts extrieurs (Fig. III.4).
Treillis tridimensionnel
Le treillis tridimensionnel le plus simple est compos de 4 nuds et de 6 barres. A chaque
nouveau nud, il faut ajouter 3 nouvelles barres pour garder lisostatisme. La rgle dun
treillis isostatique tridimentionnel est alors :
n = 4+x ;

n=4
b=6

b = 6 + 3x ;

b = 3n 6

n=4+1
b=6+3

Figure III.4 Treillis tridimensionnel

28

Calcul de treillis

2
1
(2)

(n)
n

u ()

(1)

x
Figure III.5 Mthode des nuds, quilibre dun nud.

III.1.3 Critre de dimensionnement


Les barres tant en tat de traction (NAB >0) sont dimensionnes la traction, les barres
tant en tat de compression (NAB <0) sont dimensionnes au flambage.

III.2 Mthode des nuds


La mthode permet de dterminer les efforts normaux dans les barres du treillis. On note
Ni leffort normal dans la barre i. La mthode de rsolution consiste crire successivement
lquilibre de chaque nud.

Soit un nud N connect n barres et soumis leffort extrieur F (Fig. III.5). On

note
u () le vecteur faisant un angle avec laxe
x (sens trigonomtrique positif). On

note i langle entre x et le vecteur port par la barre (i) sloignant du nud N .
En notant Ni leffort normal dans la barre (i), lquilibre du nud scrit :
n
X

Ni
u (i ) + F = 0

i=1

En crivant ainsi lquilibre de chaque nud successivement, on aboutit un systme


dquations vectorielles dont le nombre est le nombre de nuds du treillis. En deux dimensions le nombre dquations scalaires est donc 2 fois le nombre de nuds du treillis (3 fois en
tridimensionnel). Ces quations permettent de calculer les efforts normaux dans les barres
du treillis ainsi que les ractions aux appuis si le problme est isostatique.

III.3 Flambage des poutres droites


III.3.1 Introduction
Le flambage est un phnomne dinstabilit de lquilibre. Cela peut tre reprsent simplement par une bille en quilibre (fig. III.6) : la bille est en quilibre stable au fond dune

29

III.3 Flambage des poutres droites


quilibre instable

1111111111111111
0000000000000000
0000000000000000
1111111111111111
0000000000000000
1111111111111111
0000000000000000
1111111111111111
0000000000000000
1111111111111111
quilibre stable

Figure III.6 Equilibre stable, quilibre instable

F < Fcr
F = Fcr
N = F

N = F

Mf = 0

Mf 6= 0

avant flambage aprs flambage


Figure III.7 Poutre en compression

valle, en quilibre instable au sommet dune montagne. Par analogie, lquilibre est dit
stable si la structure revient sa position dquilibre aprs une petite perturbation de sa
position.
Une poutre droite en compression garde sa forme droite tant que leffort normal N est
infrieur la charge critique de flambage Fcr (Fig. III.7).

III.3.2 Charge critique de flambage dune poutre droite


Soit une poutre droite (Fig. III.8), rotule son extrmit droite, simplement appuye son

extrmit gauche, et soumise un chargement de compression F


x son extrmit gauche.
Les torseurs defforts des liaisons scrivent :
n

et

o n

x + YA
y|0 .
LA = XA

Lquilibre statique scrit


n

o n

y|0
LO = Y0

n
n
n

= 0|0 ,
+ F
x|0
+ XA
x + YA
y|0
Y0
y|0
A

En dplaant le torseur LA de A en O on a :
n

LA

n
o

= XA
x + YA
y | OA (XA
x + YA
y)

30

Calcul de treillis

L
O

F
x

(1)

v(x)
O
x

(2)

Figure III.8 tude sur la configuration dforme.


n

LA

= XA
x + YA
y | L
x (XA
x + YA
y)

LA

= XA
x + YA
y | LYA
z)

Lquilibre statique scrit alors au point O :


n

n
o
n
n

Y0
y|0
+ XA
x + YA
y | LYA
z ) + F
x|0
= 0|0
O

ce qui donne les rsultats suivants aprs rsolution :


XA = F

YA = 0 ;

Y0 = 0

On suppose que la poutre possde une configuration dforme de type flexion en quilibre
sous laction dune charge de compression. Ltude du torseur de cohsion se fait sur cette
configuration dforme (Fig. III.8). Un point O de la ligne moyenne avant dformation

devient le point O aprs dformation.


On note le dplacement de ce point v(x) :

O O = v(x)
y

Le torseur de cohsion en O est laction de (2) sur (1), en isolant (1), on a :

{K} = {Text(1) } = {F
x | 0 }0

{K} = {F
x |O O (F
x )}0

O O (F
x ) = (x
x v(x)
y ) (F
x ) = v(x)F
z
Finalement

{K} = {F
x | v(x)F
z }O

donc leffort normal vaut N = F (ngatif si F est positif) et le moment de flexion vaut
Mf = v(x)F . Or le moment de flexion vaut aussi Mf = EI (x) avec v = donc
Mf = EIv (x). Les deux expressions du moment de flexion doivent tre gales, donc :
EIv (x) + v(x)F = 0.
Comme EI est positif, si F est positif, alors cette quation diffrentielle a une solution de
la forme
v(x) = sin x + cos x

31

III.3 Flambage des poutres droites

n=1

n=2

n=3

n=4
Figure III.9 Flambage dune poutre droite articule pour diffrentes valeurs de n.

Les conditions aux points O et A imposent que ces points ne se dplacent pas verticalement :
v(0) = 0 et v(L) = 0.
Cela entraine que = 0 et que sin L = 0 donc =
lexpression de v(x) est
v(x) = sin

n
o n est un entier. Finalement,
L

n
n
n
n2 2
n
x, donc v (x) =
cos
x, et v (x) = 2 sin
x.
L
L
L
L
L

En remplaant dans lquation diffrentielle de dpart, il vient


EI

n2 2
n
n
x + F sin
x = 0.
sin
L2
L
L

x, on trouve lexpression de F qui assure une solution lquation


En simplifiant par sin n
L
diffrentielle :

n2 2
.
L2
Si F est tel quil existe un n entier qui satifasse cette dernire quation alors la configuration
F = EI

dquilibre peut tre de la forme v(x) = sin n


L x.
Dans la pratique, la charge augmente en commenant par 0. Ds que la charge est
2

gale EI L 2 , cest dire la premire valeur de n entier (n = 1), alors la poutre prend la
configuration "courbe". La valeur de ntant pas donne, cette valeur peut devenir trs
grande et conduire la ruine de la poutre.
La figure III.9 montre les dformes pour diffrentes valeurs de n.

III.3.3 lancement et rayon de giration


Le rayon de giration r caractrise la "taille" de la section de la poutre. Le rayon de giration
(exprim en m) vaut par dfinition :
r=

I
S

o I est le moment quadratique de la section et S est laire de la section.


Llancement dune poutre caractrise le rapport de sa longueur par le rayon de giration
de sa section. Llancement caractrise aussi la flxibilit de la poutre : plus est grand plus

32

Calcul de treillis
Poutre lance, grand
Poutre courte, petit
Figure III.10 lancement dune poutre droite.

la poutre est lance, plus est petit plus la poutre est courte (Fig. III.10). Llancement
(sans dimension) vaut par dfinition
L
r
o L est la longueur de la poutre et r est le rayon de giration de la section.
=

III.3.4 Critre de dimensionnement


La contrainte normale dans la poutre quand N = Fcr , appele contrainte critique de flambage vaut (n = 1)
cr =

Fcr
EI 2
= 2 .
S
L S

En faisant apparatre llancement, on a :


cr =

E 2
.
2

En notant Re la limite lastique du matriau, il y a risque de ruine par flambage si cr < Re ,


soit encore :
Si cr < Re , ruine par flambage : la charge critique de flambage est atteinte avant la
limite lastique, dimensionnement au flambage,
Si Re < cr , ruine par compression : la limite lastique est atteinte avant la charge
critique de flambage, dimensionnement en compression.
Llancement critique cr est llancement pour le lequel la charge critique de flambage
est atteinte en mme temps que la limite lastique en compression :
cr

E 2
= Re 2 = Re cr =
cr

E 2
Re

Cette valeur ne dpend que du matriau, par exemple pour un acier dusage gnral :
E = 200GPa ; Re = 240MPa cr = 90
Finalement, le critre de dimensionnement dune poutre en compression se rsume ainsi :
Si < cr ( 100 pour lacier) alors on dimensionne en compression :
le critre scrit
=|

Re
N
|<
S
s

o s est le coefficient de scurit


Si > cr ( 100 pour lacier) alors on dimensionne au flambage :

33

III.3 Flambage des poutres droites


extrmit 1

extrmit 2

Le

rotul

rotul

libre

encastr

2L

encastr

encastr

0, 5 L

encastr

rotul

0, 7 L

rotul

rotul

0, 5 L

n=2

Tableau III.1 Longueurs quivalentes suivant les conditions aux extrmits.

le critre scrit
|N | <

Fcr
s

o s est le coefficient de scurit


Si est proche de cr , des mthodes danalyse plus fines existent mais ne sont pas
dtailles ici.

III.3.5 Autres conditions aux limites


Pour dautres appuis aux extrmits, les formules restent valables en remplaant la longueur
L par une longueur quivalente Le . Le tableau III.1 donne Le suivant les cas, L dsigne la
longueur relle de la poutre. La charge critique de flambage vaut alors :
Fcr = EI

2
L2e

IV Contraintes et dformations

IV.1 Introduction
Le solide est en quilibre sous laction de forces extrieures (Fig. IV.1).
Pour connaitre ltat de contrainte lintrieur du solide, on isole un petit cube d. Les
objectifs sont de :
Caractriser les forces agissant sur le petit cube : ce sont des forces internes la matire
qui sont vues par le cube comme des forces surfaciques (Pa ou MPa) agissant sur les
6 faces ; ces forces surfaciques sont appelles contraintes,
Caractriser les dformations du petit cube : le cube sallonge (ou se rtrci) dans
chaque direction et les angles initialement de 90o entre les arrtes du cube changent,
Etablir les relations entre les contraintes et les dformations,
Proposer un critre de dimensionnement.

Figure IV.1 Solide en quilibre sous laction de forces extrieures.

36

Contraintes et dformations

y
x

d
x

Figure IV.2 Traction suivant x du petit cube.

IV.2 Caractrisation des contraintes et des dformations tridimensionnelles


On isole le petit cube d en tudiant plusieurs tats simples :

3 tats de traction dans les 3 directions


x,
y et
z,
3 tats de cisaillement dans les 3 plans.

Le repre (
x ,
y ,
z ) est align avec les arrtes du cube. La superposition de ces 6 tats
donne ltat de contrainte gnral dans lequel peut se trouver un petit lment de matire.

Etat de traction suivant


x

Ltat de traction suivant


x (Fig. IV.2) est caractris par des forces surfaciques x
x et

x appliques sur les 2 faces ayant pour normales x et x . Il est facile de vrifier que
x

le cube est bien en quilibre. On a les relations suivantes :


x = Ex ;

y = x ;

z = x

o E est le module dYoung, le coefficient de Poisson, x est la contrainte normale ap

plique suivant
x et , et sont les allongements unitaires suivant
x,
y et
z . Ces
x

relations sont celles issues de la traction dune poutre. Le coefficient de Poisson entraine
des dformations dans les directions transversales. La contrainte x a pour unit le Pa, les
dformations sont sans dimension.
On peut alors crire les dformations en fonction des contraintes :
x =

1
x ;
E

y =

x ;
E

z = x
E

Etat de traction suivant


y

De la mme faon, ltat de traction suivant


y donne les relations suivantes :
y = Ey ;

x = y ;

z = y

et on peut alors crire les dformations en fonction des contraintes :


x =

y ;
E

y =

1
y ;
E

z =

y
E

IV.2 Caractrisation des contraintes et des dformations tridimensionnelles

37

Figure IV.3 Cisaillement dans le plan (x, y) du petit cube.

Etat de traction suivant


z

Enfin, ltat de traction suivant


z donne les relations suivantes :
z = Ez ;

x = z ;

y = z

et on peut alors crire les dformations en fonction des contraintes :

x = z ;
E

y =

z ;
E

z =

1
z
E

Cisaillement dans le plan (


x ,
y)
On applique sur le cube les forces surfaciques suivantes (Fig. IV.3) :


y sur la face de normale
x


x sur la face de normale
y


y sur la face de normale
x

x sur la face de normale


y

Lquilibre en rsultante est facile vrifier. Lquilibre en moment est vrifi en crivant
la somme des moments au centre du cube.
La contrainte tangentielle engendre une distorsion angulaire du cube : langle de
avant dformation devient un angle de 2 . La relation entre et scrit :

= G
o G est le module dlasticit transversale exprim en Pa. La contrainte de cisaillement
(ou tangentielle) a pour unit le Pa, est sans dimension. Lexpression de G en fonction
de E et est :
G=

E
2(1 + )

Etant donn que cet essai est effectu dans le plan (


x ,
y ), on utilise alors les notations
suivantes :
xy = Gxy
ou encore
xy =

1
xy
G

38

Contraintes et dformations

Cisaillement dans le plan (


x ,
z)
En faisant simplement une permutation des indicea xy en xz, on a
xz = Gxz
ou encore

1
xz
G

xz =

Cisaillement dans le plan (


y ,
z)
En faisant simplement une permutation des indicea xy en yz, on a
yz = Gyz
ou encore

1
yz
G

yz =

Superposition des 6 tats


En supposant que lon applique en mme temps les 6 sollicitations simples au petit cube,
on a alors en ajoutant les contributions de chaque chargement aux dformations :
x =

x y z
E
E
E
1
y
E

y +
E
1
= xz ;
G

x +
E

z = x
E

y =

xy =

1
xy ;
G

xz

z
E
1
z
E
yz =

1
yz
G

On rappelle que
G=

E
2(1 + )

IV.2.1 Oprateur des contraintes et des dformations


Afin de simplifier les notations et de regrouper dans un mme objet les 6 contraintes dune
part et les 6 dformations dautres part, on pose les deux matrices suivantes :
On appelle oprateur des contraintes, la matrice S dfinie par

S = xy

xz

xy xz
y

yz

yz

(
x ,
y ,
z)

Les quantits sont appeles les contraintes normales, les quantits sont appeles
les contraintes de cisaillement.

39

IV.3 Problme plan

y
d

S : plan moyen

Figure IV.4 Plaque sollicite dans son plan.

On appelle oprateur des dformations, la matrice E dfinie par


xy
E=
2


xz

xz
2

xy
2

yz
2

y
yz
2

(
x ,
y ,
z)

Les quantits sont appeles les allongements unitaires, les quantits sont appeles
les distorsions angulaires.

IV.2.2 Thorme de superposition


1

1 donnant ltat de contrainte S . Soit le chargement
2 donnant ltat
Soit le chargement

2
de contrainte S
.

Ltat de contrainte S associ la somme des deux chargements et la somme des deux
oprateurs des contraintes :
S

"

1
1
x

1
1

y

S = S +S

(
x ,
y)

"

"

2
2
x

2
2

y

2
1
+ x
x

1
2

+

1
2

+

2
1
+ y
y

(
x ,
y)

(
x ,
y)

Il existe exactement le mme thorme pour loprateur des dformations.

IV.3 Problme plan


IV.3.1 Hypothses
On sintresse aux contraintes et dformations dune plaque mince charge dans son plan

(Fig. IV.5). On attache la plaque un repre (


x ,
y ,
z ) o
z est perdendiculaire la

plaque et
x et
y sont dans le plan Dans cette configuration, rien ne dpend de z.

40

Contraintes et dformations

y
y

-x
-

d
-

x
x

-y

Figure IV.5 Petit carr isol.

IV.3.2 Etat de contraintes planes

On isole un petit cube d de la plaque, dpaisseur h, align avec les axes


x et
y . Ce cube

ne subit aucune contraintes sur la faces de normales z et z , par consquent : = 0,


z

xz = 0, yz = 0.

Afin de simplifier les notations :


on note = xy ,

on reprsente le cube de dessus comme un carr dans le plan (


x ,
y ),
on reprsente les contraintes sur les cots du carr par une seule flche au centre de
larrte.
Le dessin de la figure IV.5 reprsente un carr isol de la plaque.
Le carr subit les contraintes suivantes :

contrainte normale x dans la direction


x,

contrainte normale dans la direction


y,
y

contrainte de cisaillement dans les directions


x et
y.

Loprateur des contraintes dans le plan (


x ,
y ) est not

S=

(
x ,
y)

La figure IV.6 montre le carr dform suite au chargement appliqu sur ses 4 cots.

Lallongement du carr dans la direction


x est caractris par lallongement unitaire .
x

Lallongement du carr dans la direction


y est caractris par lallongement unitaire y . La
distorsion angulaire du carr est caractris par .

A partir des relations tridimensionnelles entre les contraintes et les dformations, et en


utilisant le fait que z = 0, xz = 0, yz = 0, on obtient :
x =

x y
E
E

41

IV.3 Problme plan

x
x

Figure IV.6 Dformations du petit carr.

y =

x + y
E
E

G
Soit en inversant les quations (rsolution dun systme de 2 quations 2 inconnues) :
=

x =
y =


E
x + y
2
1


E
x + y
2
1

= G

On rappelle que
G=

E
2(1 + )

Ces relations sont valables dans le cas prsent des contraintes planes.

IV.3.3 Expressions des contraintes subies par un carr non align avec x et y

Tous les dveloppements prcdents sont faits sur un carr align avec les axes
x et
y.
Pourtant rien nempche disoler un carr non align avec ces axes. On isole par exemple

sur la figure IV.7 un carr d align avec les axes x et y inclins dun angle par rapport

x et
y.
Ce carr subit les contraintes suivantes :

contrainte normale x dans la direction x ,

contrainte normale y dans la direction y ,

contrainte de cisaillement dans les directions x et y .

On note S loprateur des contraintes dans la base ( x , y ) comme

S =

( x , y )

42

Contraintes et dformations

-x
-

x
x

d
-

-x

d
-

-y

-y

Figure IV.7 Contraintes sur un petit carr non align x et y.

ly
-x
-

M
-

-y
lx

Figure IV.8 Petit triangle en quilibre autour de M .

43

IV.3 Problme plan

Il existe des relations entre S et S . Afin de trouver ces relations, on isole autour du point
M un petit triangle dpaisseur h (Fig. IV.8). On suppose que ltat de contrainte est le
mme en tout point du triangle.
Le bilan des actions mcaniques agissante sur les 3 faces du triangle est

sur la face de normale


x : F = hl
x hl
y
1

sur la face de normale


y : F 2 = y hlx
y hlx
x

sur la face de normale x : F 3 = x hl x + hl y

Les relations entre les longueurs des cots sont


ly = cos l ,
lx = sin l .
Lquilibre en rsultante donne

F1+ F2+ F3 = 0,
soit aussi en simplifiant par h

x ly
x ly
y y l x
y lx
x + x l x + l y = 0 ,
soit en remplaant lx et ly par leurs expressions en fonction de l puis en simplifiant par l :

x cos
x cos
y y sin
y sin
x + x x + y = 0 ,

De plus, les directions x et y scrivent dans la base (


x ,
y ) sous la forme

x = cos
x + sin
y

y = sin
x + cos
y

En projetant lquilibre en rsultante sur x on a :


x cos2 cos sin y sin2 sin cos + x = 0
soit :
x = x cos2 + y sin2 + 2 sin cos

En projetant lquilibre en rsultante sur y on a :


x cos sin cos cos y sin cos + sin sin + = 0
soit :
= (y x ) cos sin + cos2 sin2

En utilisant les formules trigonomtriques suivantes


1 + cos 2
1 cos 2
cos2 =
; sin2 =
; 2 cos sin = sin 2
2
2
on a alors :
x y
x + y
x =
cos 2 + sin 2 +
2
2
x y

sin 2 + cos 2
=
2

En isolant un triangle tourn dun angle de 2 par rapport celui utilis, on en dduit que
y = x ( +

x y
x + y
)=
cos 2 sin 2 +
2
2
2

44

Contraintes et dformations

x
x

Figure IV.9 Dformations du carr non align x et y.

IV.3.4 Expressions des dformations dun carr non align avec x et y


Le carr non align avec x et y se dforme suite aux contraintes x , y et qui lui sont
appliques. La figure IV.9 montre le carr dform suite au chargement appliqu sur ses 4

cots. Lallongement du carr dans la direction x est caractris par lallongement unitaire

x . Lallongement du carr dans la direction y est caractris par lallongement unitaire y .


La distorsion angulaire du carr est caractris par .
Les relations entre les dformations (x , y , ) et (x , y , ) sont similaires celle obtenues
pour les contraintes :
x =

x y

x + y
cos 2 + sin 2 +
2
2
2

x + y
x y
cos 2 sin 2 +
2
2
2

x y

=
sin 2 + cos 2
2
2
2

y =

IV.3.5 Relation entre les contraintes et les dformations dun carr non align
avec x et y
Les relations entre les contraintes et les dformations dun carr inclin dun angle sont
les mmes que celles pour le carr non inclin :
x =

1
x y
E
E


1
x + y
E
E
1
=
G

y =

ou encore


E
x + y
2
1

E
x + y
y =
2
1

x =

= G

45

IV.3 Problme plan

x
-x


y =
n2
2

1
0

d
-


x =
n1

-2

-1

-y

Figure IV.10 Directions principales.

On rappelle que
G=

E
2(1 + )

IV.3.6 Directions principales

On dit que les directions x et y sont directions principales si le carr inclin ne subit pas

de cisaillement, autrement dit, si = 0 (Fig. IV.10). Les directions principales sont notes
et
, elles sont orthogonales.
dans la suite
n
n
1

Langle 0 tel que = 0 soit nul vrifie lquation

x y
sin 20 + cos 20 = 0
2

soit
0 =

 
1
arctan x y
2
2

et
sont alors
Les directions principales
n
n
1
2

= cos

n
1
0 x + sin 0 y

= sin

n
2
0 x + cos 0 y
Les contraintes principales 1 et 2 valent
1 = x (0 ) et 2 = y (0 )
soient

x y
x + y
cos 20 + sin 20 +
2
2
x + y
x y
cos 20 sin 20 +
2 =
2
2
1 =

IV.3.7 Cercle de Mohr des contraintes


Lorsque varie, cest dire lorsque lon fait vari linclinaison du petit carr isol, x et
voluent. Le but est de reprsenter les contraintes normale et tangentielle sur un graphique
quand varie en portant x en abscisse et en ordonne dans un repre orthonorm.

46

Contraintes et dformations

20

1 x

2
x

x +y
2

x
x y
2

Figure IV.11 Cercle de Mohr des contraintes

Pour donn, on a
x y
x + y
=
cos 2 + sin 2
2
2
x y
=
sin 2 + cos 2
2
soit aussi en ajoutant les carrs des deux expressions prcdentes
x

avec

R=

x + y  2
2
+ = R2
2
s
 2
x
y

+ 2

y
Ceci est lquation dun cercle dans le repre (x , ), de centre ( x +
2 , 0) et de rayon R.

Le cercle peut tre trac la rgle et au compas (Fig. IV.11) :


le centre C est le milieu du segment [x y ] sur laxe des abscisses
un point du cercle a pour coordonnes (x , )

Ce cercle reprsente lensemble des points (x , ) possibles dans le repre ( x , y ) quand


varie.
Les contraintes normales principales 1 et 2 sont les points dintersections du cercle avec
laxe des abscisses car la contrainte tangentielle est nulle en ces deux points. On appelle
et langle caractrisant la direction
langle caractrisant la direction principale
n
1

. On sait que langle tel que soit nul vrifie


principale
n
2
0

cot oppos
tan 20 = =
x
y
cot adjacent
2
Graphiquement, on mesure langle 20 comme langle reliant dans cet ordre les trois points
(x , 0), C puis (x , ). Le sens positif est le sens trigonomtrique.
Dans tous les cas, par convention, on prend 1 droite et 2 gauche.
Si y < x (Fig. IV.12) alors 1 = 0 et 2 = 0 + 2 ,
Si x < y (Fig. IV.13) alors 1 = 0 +

et 2 = 0 .

47

IV.3 Problme plan

20
2

1 x

2 y

C
0 > 0

0 < 0

Figure IV.12 Cas o y < x : 1 = 0 et 2 = 0 +

20

1 x

20
2

1 x

2 x

xx

yy

1 x

C
0 > 0

20

Figure IV.13 Cas o x < y : 1 = 0 +

et 2 = 0

0 < 0

V Critres de dimensionnement

V.1 Objectifs
Le but est de vrifier que les contraintes dans la structure restent acceptables pour ne pas
engendrer de rupture en fonctionnement. Pour les matriaux ductiles, les critres utiliss
couramment imposent que le matriau reste dans le domaine lastique en tout point de
la structure, cest le cas des critres de Tresca et de Von Mises. De plus, le phnomne
de fatigue est le critre de dimensionnement prendre pour les pices subissant un grand
nombre de cycles de chargement.

V.2 Matriaux ductiles : critre de Tresca


Ce critre est d Henri Edouard Tresca (1814-1885) qui fut professeur titulaire de la chaire
de mcanique du Cnam. Il a observ que le facis de rupture dune prouvette casse suite
un chargement de traction tait inclin 45o par rapport laxe de traction pour les
matriaux ductiles. Or le cisaillement tant maximal pour cet angle, il en a dduit que la
rupture se fait par glissement engendr par les contraintes de cisaillement. Le critre est
donc bas sur le cisaillement maximal.

Ltat de contrainte la limite lastique dune prouvette de traction daxe


x est

Re 0 0

S=
0
0

0 0

0 0 (

x ,
y ,
z)

Le cercle de Mohr associ cet tat de contrainte (Fig. V.2) donne le cisaillement maximal
Re
max =
. Le critre de Tresca est
2

50

Critres de dimensionnement

xx

zone utile de lprouvette

zone
plastique

zone de striction

xx

zone lastique

xx

zone
dcrouissage

45o

facis de rupture 45o

xx

Figure V.1 Essai de traction rupture.

max =

Re
2

Re
2

x
Re

Figure V.2 Cercle de Mohr la limite lastique pour un essai de traction.

max <

Re
.
2

Pour crire le critre correctement il faut raisonner en tridimensionnel. En tridimensionnel,


on peut trouver 3 directions principales associes 3 contraintes principales 1 , 2 et 3
,
,
), loprateur des contraintes scrit
(Fig. V.3). Dans la base principale (
n
n
n
1

S=
0

(
n
1 n2 ,n3 )

Pour les 3 couples de valeurs (1 , 2 ), (1 , 3 ) et (2 , 3 ) on peut tracer 3 cercles de Mohr


(Fig. V.4). Cest ce que lon appelle le tri-cercle de Mohr. La zone possible pour les couples

n
pla

p
- 3 lan
1 -
2
2

pla
n
1 -
3
Figure V.3 Trois plans des contraintes principales en tridimensionnel.

51

V.3 Matriaux ductiles : critre de Von Mises

Zone possible

max

dtat de contrainte
Cercle (1 , 3 )
Cercle (1 , 2 )

x
3

1
Cercle (1 , 3 )

Figure V.4 Tri cercle de Mohr.

(x , ) est la zone grise comprise entre les trois cercles. Le cisaillement maximal dans ce
cas est
max =

1 3
.
2

Suivant les valeurs et les signes des 1 , 2 et 3 , le critre scrit en tridimensionnel


 | | | | | | 
1
2
1
3
2
3

max = Max

<

Re
.
2

En notant T resca = 2max la contrainte quivalente de Tresca, le critre scrit




T resca = Max |1 2 |, |1 3 |, |2 3 | < Re .


Dans le cas des contraintes planes, tant donn que 3 = 0, le critre devient


T resca = Max |1 2 |, |1 |, |2 | < Re .

V.3 Matriaux ductiles : critre de Von Mises


Le critre de Von Mises (1883-1953), plus rcent, est bas sur lnergie de dformation que
le matriau peut stocker avant plastification. On note V.M. la contrainte quivalente de Von
Mises. Elle vaut dans le cas tridimensionnel
1 q
(1 2 )2 + (1 3 )2 + (2 3 )2 ,
V.M. =
2
soit dans le cas des contraintes planes avec 3 = 0
V.M.

1 q

=
(1 2 )2 + 12 + 22 .
2

Le critre de Von Mises scrit


V.M. < Re .

52

Critres de dimensionnement

Re

Re

Re
1

Re

Re

Re

a) zone de Tresca

Re

Re

1
Re

Re

b) zone de Von Mises

Re

Re

c) comparaison

Figure V.5 Comparaison des critres de Tresca et de Von Mises.

V.4 Comparaison des critres de Tresca et de Von Mises


On compare les critres de Tresca et de Von Mises dans le cas des contraintes planes. Dans
le plan (1 , 2 ), le lieu des points tels que T resca < Re est
|1 | < Re

Re < 1 < Re ,

|2 | < Re

Re < 2 < Re ,

|1 2 | < Re

1 2 < Re et 2 1 < Re .

Les deux premires lignes donnent un carr centr en (0, 0) de cot 2Re . Les deux conditions suivantes donnent une zone comprise entre 2 droites de pente 1 et ayant pour ordonnes
lorigine Re et Re (). La zone de Tresca est reprsente sur la figure V.5a).
Dans le plan (1 , 2 ), le lieu des points tels que V.M. < Re est
12 + 22 1 2 < Re2 .
Ces points sont lintrieur dune ellipse centre en (0, 0) dont quelques points sont donns
pour faciliter le trac :
(Re , 0) (0, Re ) (Re , Re ) (Re , Re )
La zone de Von Mises est reprsente sur la figure V.5b).
On remarque que la zone de Von Mises est plus grande que celle de Tresca, la diffrence
est reprsente sur la figure V.5c). Le choix du critre (Von Mises ou Tresca) se fait partir
de rsultats dessais.

V.5 Fatigue des matriaux


La fatigue est un phnomne de dtrioration dune pice conscutive un grand nombre de
cycles de chargements (> 1000) alors que les contraintes sont dans le domaine lastique (les
critres de Tresca et de Von Mises sont vrifis). Schmatiquement, chaque cycle de chargement, des micro fissures apparaissent et grandissent dans les zones fortes concentrations
de contraintes. Ces fissures se propagent et peuvent engendrer la ruine de la structure.

53

V.5 Fatigue des matriaux

0a 90%
MPa

F/S
0a

50%

temps

250
Rr 10%
200

0a

150
Re
100
50
Rendu.
nb. cycles

100

101

102

103

104

105

106

107

108

Figure V.6 Courbes de Whler : diffrentes probabilits de dfaillance.

F/S
moy + a
moy

temps

moy a
Figure V.7 Chargement cyclique moyenne non nulle

La courbe de Whler, obtenue de faon exprimentale sur une prouvette sollicite en


traction, donne la contrainte applique en fonction du nombre de cycle de chargements
rupture (Fig. V.6). Les courbes sont donnes pour un chargement cyclique de moyenne
nulle avec une amplitude de 0a . En ralit une tude statistique doit tre mene sur un
ensemble dprouvettes. On peut alors tracer plusieurs courbes pour diffrentes probabilits
de dfaillance.
On utilise le plus souvent une seule courbe, celle ayant 10% de probabilit de dfaillance
par exemple. On peut caractriser plusieurs valeurs sur cette courbe :
limite la rupture Rr : rupture aprs 1 cycle de chargement,
limite lastique Re : le nombre de cycles rupture se situe aux alentours de 103 104 ,
limite dendurance Rendu. : valeur pour laquelle le nombre de cycles devient infini,
elle est comprise entre 0, 3 Rr et 0, 6 Rr suivant les matriaux. Pour les al-

liages non ferreux, lasymptote horizontale nexiste pas, on prend alors par convention
Rendu. = 0, 5 Rr .

Les courbes de Whler dpendent cependant dautres facteurs plus ou moins difficiles
prendre en compte sans faire dessais supplmentaires : taille de la pice, tat de surface,
corrosion, traitement de surface, temprature.

54

Critres de dimensionnement

Dans le cas dun chargement cyclique contrainte moyenne moy non nulle (Fig. V.7),
on peut appliquer des formules empiriques permettant de trouver la contrainte quivalente
cyclique moyenne nulle. Par exemple, la formule de Soderberg donne une relation entre
lamplitude a de la contrainte appique de moyenne non nulle moy lamplitude de la
contrainte quivalente 0a du chargement moyenne nulle :
a
moy
+
=1
0a
Re
Dans le cas dun tat de contrainte plan ou tri-dimensionnel, on utilise la contrainte
quivalente de Von Mises ou bien la contrainte quivalente de Tresca suivant le matriau :
le critre adapt est dtermin par des essais.

VI Enveloppes minces

VI.1 Action dun fluide au repos sur un solide


On suppose une surface dun solide en contact avec un fluide au repos dans lequel la

pression est uniforme. En chaque point M de la surface , laction du fluide est p


n o

n est la normale unitaire sortante de en M (Fig. VI.1).

La rsultante de laction du fluide sur le solide travers vaut :

s f s =

ZZ

p
n d

Lorsque la pression est uniforme, cette rsultante devient :

s f s = p

ZZ

n d

Le calcul apparemment compliqu de lintgrale double est en fait simple. On peut mon

trer que la projection de


s
sur une direction arbitraire d vaut :
f s

s f s d = pS

o S est laire de la surface projete de sur un plan de normale unitaire d . Le signe +


ou est dfinir suivant les cas mais est souvent facile trouver par intuition. De faon

gnrale :

si d est un vecteur rentrant dans , on a


s f s . d = +pS ;

si d est un vecteur sortant de (cas du dessin), on a


s f s d = pS.

56

Enveloppes minces

Solide

Fluide, pression p

Surface projete plane S

Figure VI.1 Action dun fluide sur un solide - Surface projette.

VI.2 Application un rservoir cylindrique


On souhaite dimensionner un rservoir mince cylindrique sous pression (Fig. VI.2). Le cylindre a pour rayon moyen R, sa longueur est L, son paisseur est e, la pression intrieure
est note p. On note S la surface intrieure du rservoir en contact avec le fluide. Le repre

orthonorm (
e ,
e ,
z ) attach au point M est tel que
z soit port par laxe du rservoir.
r

On suppose que les vecteurs contrainte pour les directions


er ,
e et
z valent

T (M,
er ) = 0

T (M,
e )=
e

T (M,
z ) = zz
z

En coupant le rservoir par une section fictive AA perpendiculaire laxe (section A-A

de normale
z ), en isolant la partie et en crivant lequilibre de en rsultante, il vient
1

ZZ

AA

zz
z p

ZZ

S1

n dS = 0 .

En projetant suivant la direction


z , la surface projete de S1 est un disque de rayon R, son
aire est R2 . La surface de coupe fictive AA est un anneau de rayon moyen R et dpaisseur
e. Le rayon tant beaucoup plus grand que lpaisseur, une bonne approximation de son aire
est de remplacer lanneau par un rectangle de longueur gale au primtre et de largeur gale

lpaisseur du rservoir, soit 2Re. Finalement, lquation prcdente, projette sur


z,
devient
2Rezz R2 p = 0,
soit la contrainte axiale

pR
.
2e
Cette contrainte ne dpend pas de la forme du fond du rservoir et est nomme leffet de
zz =

fond.

57

VI.2 Application un rservoir cylindrique

En coupant le rservoir dans le sens de la longueur par une section fictive BB dans le

plan (
e ,
z ) (section B-B de normale
e ), en isolant la partie 1 et en crivant lquilibre de
r

1 en rsultante, il vient
ZZ

BB


e p

ZZ

S2

n dS = 0 .

Dans le calcul des aires qui suit, on nglige les aires des fonds sphriques par rapport
aux aires faisant intervenir la longueur du rservoir. Cela revient isoler un tronon du

reservoir autour du point M situ loin des fonds. En projetant suivant la direction
e , la

surface projete de S1 est un rectangle de longueur b et de largeur 2R, son aire est 2LR. La
surface de coupe fictive BB est compose de 2 rectangles de longueurs b et de largeurs e,

son aire est 2Le. Finalement, lquation prcdente, projette sur


e devient

2Le 2LRp = 0,
soit la contrainte circonfrentielle

pR
.
e
La contrainte circonfrentielle est deux fois plus importante que la contrainte axiale, cest
=

elle qui dimensionne le rservoir. Les formules prcdentes donne une bonne approximation
de ltat de contrainte dun rservoir mince loin des fonds si lpaisseur est petite devant les
dimensions du rservoir, soit peu prs e <

R
40

et R < 3L.

A-A

S1
AA

111111
000000
000000
111111
000000
111111

000000
111111
er
000000
111111
000000
111111
000000
111111
000000
111111
000000
111111
000000
111111
000000
111111
e

A
L
B-B
111111
000000
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111

z
S2
000000
111111
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111
BB
000000
111111
0000000000000000
1111111111111111
000000
111111
000000
111111
0000000000000000
1111111111111111
000000
111111

er
Figure VI.2 Rservoir mince cylindrique en pression interne.

VII Initiation au calcul lments finis

VII.1 tude de llment de barre


VII.1.1 quilibre de llment barre
On dsigne par barre une poutre travaillant seulement en traction-compression. Typiquement, les treillis de poutres sont souvent approxims dans un premier temps comme un
ensemble de barres rotules entre elles (Fig. VII.1). Les 2 extrmits sont les nuds de
llment.
Un lment de barre et reprsent par un segment de droite reliant les deux extrmits
de la barre (Fig. VII.2). La barre est caractrise par sa longueur l, par laire de sa section
S et par son module dYoung E. Les extrmits sont appels les nuds de la barre, ils ont

Figure VII.1 Exemple de calcul par lments finis dun pylne lectrique (lments barres).

60

Initiation au calcul lments finis

Avant dformation :

u1

u2
O

Aprs dformation :

x
F2

x
F1
Figure VII.2 Elment barre 2 noeuds.

pour abscisses x1 = 0 et x2 = l.

La barre est en quilibre sous les actions des forces F1


x et F2
x aux nuds 1 et 2. On

note les dplacements des nuds 1 et 2 respectivement u x , u x .


1

Le torseur de cohsion en O peut scrire de deux faons :


n

KO = ext. 2 N = F2

ou

KO = ext. 1 N = F1

On remarque que la barre est bien en quilibre car F1 + F2 = 0.


La relation entre N , u1 et u2 pour une poutre en traction est
N=

ES
ES
l N =
(u2 u1 )
l
l

En utilisant les rsultats prcdents on a alors


F1 =
F2 =

ES
ES
(u2 u1 ) F1 =
(u1 u2 )
l
l

ES
ES
(u2 u1 ) F1 =
(u1 + u2 )
l
l

Ces deux dernires relations scrivent sous la forme matricielle suivante :


"

F1
F2

ES
=
l

La matrice
ES
k=
l

"

"

1 1

1 1

#"

u1
u2

est appelle la matrice de rigidit de llment barre.

VII.1.2 Exemple dapplication


On prend une poutre encastre gauche (u1 = 0), et on applique un effort F lextrmit
droite (F2 = F ). Le systme matriciel rsoudre est le suivant :
ES
l

"

1 1

#"

u1
u2

"

F1
F2

61

VII.2 tude de deux barres

qui devient ici


ES
l

"

1 1

#"

0
u2

"

F1
F

o les inconnues sont le dplacement de lextrmit droite u2 et la force de raction de


lextrmit gauche F1 . Ce systme de 2 quations 2 inconnues donne la solution
u2 =

Fl
ES

et
F1 = F
Leffort normal est dtermin par
N=

ES
ES F l
(u2 u1 ) =
(
0) = F
l
l ES

VII.1.3 Remarques sur la mthode des lments finis


La mthode des lments finis est base sur lcriture de lquilibre des lments. La rsolution dun problme par lments finis permet de dterminer les inconnues defforts de liaisons
(rsolution du problme de statique) et les inconnues de dplacements et defforts normaux
(rsolution du problme de r.d.m.) ; que le problme soit isostatique ou hyperstatique.

VII.2 tude de deux barres


VII.2.1 Assemblage des matrices de rigidit lmentaires
On suppose deux barres de longueurs, de modules dYoung et de sections diffrentes colles
bout bout et soumises de la traction (Fig. VII.3). Les barres sont numrotes I et II,
elles sont relies trois nuds 1, 2 et 3. Ces trois nuds subissent les forces extrieures

F1
x , F2
x et F3
x . Lquilibre global scrit
F1 + F2 + F3 = 0.
A lquilibre, lensemble des deux barres sest dform, les nuds 1, 2 et 3 se sont dplacs

respectivement de u
x, u
x et u
x (u < u < u si les deux barres sont en traction).
1

Les quilibres des nuds scrivent :


quilibre du nud 1 (mthode des nuds)
F1 + N1 = 0 N1 = F1
quilibre du nud 2
F2 N 1 + N 2 = 0 N 1 N 2 = F2
quilibre du nud 3
F3 N 2 = 0 N 2 = F3
Pour chaque lment, on peut relier leffort normal lallongement de la barre :

62

Initiation au calcul lments finis

E1 , S1

E2 , S2

u1
x

u2
x

u3
x

tat dform

N1
x

II

F1
x 1 N1
x

N1
x

x
N1
x 2 N2

x
F2

N2 x
II

N2
x

N2
x 3 F3
x
Figure VII.3 Deux barres en traction.

lment I :
N1 =

E1 S1
(u2 u1 )
L1

N2 =

E2 S2
(u3 u2 )
L2

lment II :

En remplaant dans les expressions des quilibres des nuds, on a les trois quations
suivantes :

E1 S1
(u2 u1 ) = F1
L1

E2 S2
E1 S1
(u2 u1 )
(u3 u2 ) = F2
L1
L2
E2 S2
(u3 u2 ) = F3
L2
En notant k1 = E1 S1 /L1 et k2 = E2 S2 /L2 , on a alors le systme matriciel suivant :

F1
k1
k1
0
u1


k1 k2 + k1 k2 u2 = F2


F3
0
k2
k2
u3

{z
K

} | {z }
Q

| {z }
F

Cette opration est lopration dassemblage des matrices de rigidit lmentaires, la matrice
K est appele matrice de rigidit de la structure, le vecteur Q est le vecteur des inconnues

63

VII.3 lment barre pour le calcul des treillis

de dplacements et le vecteur F est le vecteur des forces extrieures :


KQ = F

VII.2.2 Mise en uvre pratique


La premire tape consiste crire les deux matrices de rigidit des deux lments

k1 =

"

u1
k1

u2
#
k1 u1

k1

k1

et

k2 =

u2

"

u2
k2

k2

u3
#
k2 u2
k2

u3

en reprant les lignes et les colonnes de chaque matrice par les inconnues de dplacements
associes. On range ensuite dans la matrice de rigidit K de la structure chaque terme des
deux matrices la ligne et la colonne correspondante :

u1

u2

u3

k1

k1

K =
k1

k1 + k2

k2

u1

k2 u2

k2

u3

Le systme rsoudre est alors KQ = F .


La deuxime tape consiste faire le bilan des dplacements et des forces connus et
inconnus. En prenant un encastrement lextrmit gauche et en appliquant une force F
lextrmit droite, on a :

u1

= 0, connu

Q=
u2 inconnu

u3

u1

inconnu, raction lencastrement

F =
u2 = 0, connu

inconnu

u3

= F, connu

Si le dplacement est connu en un nud alors la force est inconnue, si la force est connue
alors le dplacement est inconnu.
La troisime tape est la rsolution du systme dquations complet afin de dterminer
toutes les inconnues

k1

k1

k1

F1

k2 + k1 k2
u2 = 0
F
u3
0
k2
k2

Une toutes les inconnus trouves, on peut calculer les efforts normaux dans chaque barre :
N1 =

E1 S1
(u2 u1 )
|{z}
l1
=0

et
N2 =

E2 S2
(u3 u2 ).
l2

64

Initiation au calcul lments finis

u2Y

u2X
u1Y

1
u1X

Figure VII.4 Elment barre dans une base globale.

VII.3 lment barre pour le calcul des treillis


Les barres composant un treillis sont positionnes arbitrairement dans lespace et font des

angles diffrents avec le repre global de la structure ( X , Y ) (Fig. VII.4).

On note langle entre laxe X du repre global et laxe


x du repre local la barre.
Le vecteur dplacement dun point de la barre scrit dans le repre local

u = u
x.
Il scrit dans le repre global

u = uX X + uY Y .

En projetant les deux quations prcdentes sur


x il vient
u = uX cos + uY sin .

Y
En notant uX
1 et u1 les dplacements suivant X et Y du nud 1 de la barre dans le repre
global, et en appliquant la formule prcdente au nud 1, on a
Y
u1 = uX
1 cos + u1 sin .

En utilisant les mmes notations pour le nud 2, on a


Y
u2 = uX
2 cos + u2 sin .

Ceci peut scrire sous la forme matricielle suivante


"

u1
u2

"
|

cos sin
0

cos
{z

=T

soit

uX
1

#
uY
1
X
sin
u2
}
Y

u2

| {z }
=Q

q = TQ
o Q est le vecteur des inconnus de dplacements aux nuds de llment dans le repre
global et T est la matrice de transformation passant du repre global au repre local.

65

VII.4 Elment de poutre pour le calcul des portiques

Il est possible dcrire les mmes relations pour les forces extrieures agissant aux nuds
de llment :

F 1 = F1
x = F1X X + F1Y Y et F 2 = F2
x = F2X X + F2Y Y
donc

F1X = F1
x X = F1 cos et F1Y = F1
x Y = F1 sin

de mme,
F2X = F2 cos et F2Y = F2 sin
ce qui scrit sous forme matricielle

F1X

cos


sin

=
0

FY
1
X
F
2

F2Y

"
#
F
1

cos
F2

sin

{z

=TT

La matrice qui apparait pour les forces est la transpose de celle prsente dans les relations
des dplacements. Finalement, lquilibre de la barre crit en fonction des dplacements et
des forces dans le repre local la barre
"

F1
F2

ES
=
l

"

1 1

#"

u1
u2

devient en fonction des dplacements et des forces dans le repre global

F1X

FY
1
X
F
2

F2Y

cos


sin

=
0

"
ES
1

l
1
cos

#"

cos sin
0

uY2

sin

Tous calculs faits, lquilibre de la barre en deux dimensions scrit

ES
l
|

cos2

cos sin

cos sin

sin2

cos2

cos sin

cos sin

sin2

{z

=kg

cos2

cos sin

cos2

cos sin

cos sin

sin2

cos sin

uX
1

#
uY
1
X
sin
u2

cos

sin2

uX
1

uY
1
X
u
2

uY2

F1X


FY
1
= X
F
2

F2Y

o la matrice kg est la matrice de rigidit de llment barre en deux dimensions.


Leffort normal est dtermin par
N=


ES
ES  X
Y
Y
(u2 u1 ) =
(u2 uX
)
cos

+
(u

u
)
sin

1
2
1
l
l

VII.4 Elment de poutre pour le calcul des portiques


Un portique est constitu dun ensemble de poutre assembles entre elles par des liaisons
encastrements (Fig. VII.5).

66

Initiation au calcul lments finis

V
II

IV

III

1111111
0000000
0000000
1111111

Figure VII.5 Exemple de portique discrtis par des lments poutres

(x)

1
w1

w2
x

w(x)
x

Figure VII.6 Elment de poutre en flexion

Soit un lment de poutre de longueur l, de section S, de module dYoung E et de


moment dinertie de section I (Fig. VII.6). Les deux extremits de la poutre sont les nuds
1 et 2 dabscisses respectives 0 et l.
On utilise les notations suivantes :

w(x) : dplacement suivant


y de la section dabscisse x,
(x) : rotation de la section dabscisse x,

w1 : dplacement suivant
y de la section au nud 1,

w2 : dplacement suivant y de la section au nud 2,


1 : rotation de la section au nud 1,
2 : rotation de la section au nud 2,

F1 : force extrieure suivant


y au nud 1,

F : force extrieure suivant


y au nud 2,
2

M1 : moment extrieur suivant


z au nud 1,

M2 : moment extrieur suivant


z au nud 2.
En suivant la mme dmarche que pour llment barre, on peut montrer que lquilibre
dun lment scrit

2EI
l3

3l 2l2

6 3l

3l

3l

{z
k

3l
6

3l

3l

} | {z }

w1

F1


l2
1 M 1


3l
w2 F 2

2l2

M2

2
q

{z
q

VIII Moyens exprimentaux

VIII.1 Jauges de dformation


VIII.1.1 Principe
Le principe utilis par les jauges de dformation est que la rsistance lectrique de certains
fils varie lorsquils sont tirs (fig. VIII.1). La jauge est colle la surface de la pice, la
dformation de la pice est alors relie directement la variation de rsistance lectrique de
la jauge.
La rsistance lectrique R dune jauge est relie la rsistivit du mtal composant la
jauge , la longueur du fil L et la section du fil S par :
L
.
S
Des variations , L et S des trois cractristiques de la jauge provoquent une variation
R=

R de la resistance. Cette variation R se calcule par


R =

R
R
R
+
L +
S

L
S

Aprs calcul des drives et division par R, on trouve


L S
R
=
+

L
S
En supposant que la section du fil est un carr de cot h, un allongement relatif du fil de
L
L

provoque une dformation transversale de h


o est le coefficient de Poisson du
h

matriau de la jauge. La section du fil S vaut S = h2 . La variation h de la largeur du fil


et la variation de section S sont relies par
S =

S
h = 2hh
h

68

Moyens exprimentaux
avant dformation

aprs dformation

allongement du fil
L

Figure VIII.1 Jauge de dformation.

Va
i1

Vd

Rd

Rb Vb

i2

i1
i2

Ra

Rc

Vc
Figure VIII.2 Pont de Wheatstone.

soit aussi

h
S
=2
S
h
En remplaant par la relation avec L, on trouve
L
S
= 2
S
L
Finalement, en ngligeant la variation de rsistivit du matriau qui est trs petite en
pratique, on a

R
L
=
(1 + 2)
R
L
On note k = 1 + 2 le facteur de jauge. Ce facteur est donn par le fabricant de la jauge, il
est de lordre de 0, 5 3. On a alors

L
R
=k
= k
R
L
o est lallongement unitaire dans la direction de la jauge.

VIII.1.2 Pont de Wheatstone


Les variations de rsistance lectrique des jauges sont trs petites, pour les mesurer, on
utilise un montage en pont de Wheatstone. Le pont de Wheatstone est un montage de
quatre rsistances pour lequel on impose la tension dentre V et on mesure la tension de
sortie S (Fig. VIII.2) :
Va = Ra i1 ; Vb = Rb i2 ; Vc = Rc i2 ; Vd = Rd i1 ; V = Req i

69

VIII.1 Jauges de dformation

avec

1
1
1
=
+
.
Req
Ra + Rd Rb + Rc

i = i1 + i2

S =V V

a
b

S =V V
c
d

V
Req

= i1 + i2

(1)

S = Ra i1 Rb i2 (2)
S = Rc i2 Rd i1 (3)
(1) i1 =

V
i2
Req

En remplaant dans (2) :


S = Ra

 V

Req

i2 =
En remplaant dans (3) :

R

Req

Ra
V (Ra + Rb )i2
Req

V S

1
Ra + Rb

S=



V
Rd  Ra
Rc  Ra
V S Rd
+
V S
Ra + Rb Req
Req Ra + Rb Req

Rc
Rd 
Ra Rd 
V  Ra Rc
+
Rd +
=
Ra + Rb Ra + Rb
Req Ra + Rb
Ra + Rb

S 1+
S

i2 Rb i2 =

Ra + Rb + Rc + Rd
Ra + Rb + Rc + Rd  Ra Rc + Ra Rd Ra Rd Rb Rd 
=V
Ra + Rb
(Ra + Rd )(Rb + Rc )
Ra + Rb

Finalement

S=V

Ra Rc Rb Rd
(Ra + Rd )(Rb + Rc )

Le pont est dit quilibr si S = 0, cest dire si Ra Rc = Rb Rd ou encore si


Ra
Rb
=
.
Rd
Rc
On suppose le pont quilibr et on cherche quelle est la variation de tension de sortie
S quand les rsistances varient respectivement de Ra , Rb , Rc et Rd . On sait que
au premier ordre on a
S =

S
S
S
S
Ra +
Rb +
Rc +
Rd .
Ra
Rb
Rc
Rd

Le calcul des drivs donne




S
Rb
Rb Rc
=V
+ (Ra Rc Rb Rd )
Rd
(Ra + Rd )(Rb + Rc )
(Ra + Rd )2 (Rb + Rc )2

mais comme le pont est quilibr, on a Ra Rc Rb Rd = 0 ce qui donne




Rb
S
=V
.
Rd
(Ra + Rd )(Rb + Rc )

De mme



Rc
S
=V
,
Ra
(Ra + Rd )(Rb + Rc )

70

Moyens exprimentaux

quilibrage du pont
rglage du facteur de jauge

Figure VIII.3 Boitier contenant le montage en pont de Wheatstone.




Rd
S
=V
,
Rb
(Ra + Rd )(Rb + Rc )



Ra
S
=V
.
Rc
(Ra + Rd )(Rb + Rc )

Finalement
S =



V
Rb Rd + Rc Ra Rd Rb + Ra Rc .
(Ra + Rd )(Rb + Rc )

Pour simplifier les calculs dans la suite, on prend les mmes rsistances nominales pour
les quatres rsistances :
R = Ra = Rb = Rc = Rd .
Ceci se justifie car les jauges branches sur le pont de wheatstone sont gnralement les
mmes. La valeur de S se simplifie
S =

V  Rd Ra Rb Rc 
+

.
4
R
R
R
R

En remarquant que
Ra
Rb
Rc
Rd
= ka ;
= kb ;
= kc ;
= kd ;
R
R
R
R
o k est le facteur de jauge, il vient
S =
ou encore


V 
k d + a b + c ,
4

4S
= a b + c d .
kV

VIII.1.3 Utilisation du boitier


Le boitier utilis permet de faire le montage du pont de Wheatstone et dafficher sur lcran
la grandeur caractristique de la variation de tension (Fig. VIII.3). Dans le cas o moins
de 4 jauges sont utilises, il est possible de les remplacer par des rsistances prsentes dans
le boitier et ayant des valeurs de rsistance nominales quivalentes aux jauges branches.

71

VIII.1 Jauges de dformation

Figure VIII.4 Exemple de montage en pont complet : essai de traction.

Aprs branchement des jauges, il faut indiquer le facteur de jauge k et quilibrer le pont.
Le boitier impose une tension V de 2 Volts.
La mesure affiche sur lcran la valeur A
A=

4S 6
.10
kV

soit aussi
A = (a b + c d ).106 .
Le facteur daffiche 106 est due au fait que les dformations sont souvent comprises entre
106 et 103 .

VIII.1.4 Diffrents montages


Exemple dun pont complet
Les 4 jauges sont actives. En quipant judicieusement une prouvette de traction de 4 jauges
montes en pont complet (Fig. VIII.4), on peut corriger les effets suivants :
dilatation due la temprature,
flexion parasite.
On note t la dformation due la traction, f la dformation due la flexion sur la face
suprieure et f sur la face infrieure, enfin d est la dformation due la dilatation.
La disposition des jauges est la suivante :

a : face suprieure, direction de traction a = t + f + d


b : face suprieure, direction transversale b = t + d
c : face infrieure, direction de traction c = t f + d
d : face infrieure, direction transversale d = t + d
laffichage indique alors
A=

4S 6
.10 = a b + c d = 2(1 + )t
kV

si k est le facteur de jauge. Si lon souhaite que laffichage donne t , il suffit de multiplier k
par k de faon avoir :
A=

4S
2(1 + )
.106 =
t = t
kk V
k

cest dire k = 2(1 + ). Le rglage du facteur de jauge sur le boitier est alors 2k(1 + ).

72

Moyens exprimentaux

a
b
Figure VIII.5 Exemple de montage en demi-pont : essai de flexion.

Exemple dun demi-pont


Seulement 2 jauges sont actives. En quipant judicieusement une prouvette de flexion de
2 jauges montes en demi-pont (fig. VIII.5), on peut corriger les effets suivants :
dilatation due la temprature,
traction parasite.
On note t la dformation due la traction, f la dformation due la flexion sur la face
suprieure et f sur la face infrieure, enfin d est la dformation due la dilatation.
La disposition des jauges est la suivante :

a : face suprieure, direction axiale a = t + f + d


b : face infrieure, direction axiale b = t f + d
laffichage indique alors
A=

4S 6
.10 = a b = 2f
kV

si k est le facteur de jauge. Si lon souhaite que laffichage donne f , il suffit de multiplier k
par k de faon avoir :

4S 6
2
.10 = f = f

kk V
k
cest dire k = 2. Le rglage du facteur de jauge sur le boitier est alors 2k.
A=

Quart de pont
Une seule jauge est active. Les 3 autres rsistance doivent avoir les mmes valeurs nominales
de rsistance que la jauge. Si seule la jauge a est active, laffichage indique
A=

4S 6
.10 = a .106
kV

Laffichage donne directement lallongement unitaire de la jauge a en m/m si k est le facteur


de jauge. Lutilisation dune seule jauge ne permet pas de corriger les effets parasites comme
linfluence des variations de temprature.

VIII.1.5 Capteurs jauges


En vue de mesurer des grandeurs physiques (forces, acclrations, pressions ...), il est possible
de les faire agir sur un corps dpreuve dont les dformations sont proportionnelles la
grandeur mesure. La figure VIII.6 montre un capteur de force utilis dans une balance de
cuisine.

73

VIII.1 Jauges de dformation

Figure VIII.6 Capteur de force dune balance de cuisine.

y
90o

45o 90 = y

45

0o

0 = x

3 jauges 0o , 45o et 90o


Figure VIII.7 Rosette de 3 jauges 45o .

Ce type de capteur est prcis, fidle et permet une certaine souplesse demploi. On
trouve ces capteurs dans des applications comme la pese, la mesure de pression, la mesure
de force, la mesure de couple. Les domaines dapplications sont vastes : automobile, mdical,
instruments de mesures...
Suivant les applications, les montages en quart de pont, demi-pont ou pont complet
peuvent tre utiliss.

VIII.1.6 Exploitation dune rosette de 3 jauges 45o


Une rosette 45o est constitue de 3 jauges colles en un mme point et mesurant les
allongements unitaires dans 3 directions distinctes 45o lune de lautre (Fig. VIII.7).

On place le repre (
x ,
y ) de telle manire aligner
x et
y avec les jauges perpendiculaires. On note 0 , 45 et 90 les allongements unitaires mesurs par les 3 jauges. Donc on
a les relations suivantes :
0 = x
90 = y
ainsi que
45 = ( = 45o )
x y
x + y
cos (2 45o ) + xy sin (2 45o ) +
2
2
x + y
45 = xy +
2

On peut finalement en dduire les dformations au point de mesure dans la base (


x ,
y ) en
45 =

fonction des allongements mesurs :


x = 0

74

Moyens exprimentaux

Figure VIII.8 Dispositif de photolasticit, tude des isochromatiques dun crochet soumis
de la traction.

y = 90
xy = 45

0 + 90
2

On peut ensuite en dduire loprateur des contraintes au point de mesure.

VIII.2 Photolasticit
VIII.2.1 Principes
La photolasticit permet dtudier les contraintes dans des pices planes en polymre transparent par un systme optique (Fig. VIII.8).
Dfinition lectromagntique de la lumire
Les phnomnes lumineux peuvent, selon la thorie lectromagntique, tre considrs comme
lis la propagation simultane dun champ lectrique E et dun champ magntique H,
constamment perpendiculaires entre eux ainsi qu la direction de propagation, et dont les
valeurs sont des fonctions sinusodales du temps.
Lumire polarise
Un filtre polarisant possde la proprit de ne laisser passer quune composante du champ
parallle une direction fixe dite axe de polarisation. Deux filtres polarisants successifs
axes parallles laissent passer la lumire ; sils ont croiss, cest dire axes perpendiculaires,
ils ne laissent pas passer la lumire, le faisceau polaris par le premier ayant une composante
nulle suivant laxe du second (Fig. VIII.9).

75

VIII.2 Photolasticit

polariseur 1

lumire
polarise

direction d 1

lumire quelconque

d1

polariseur 1
lumire quelconque

d1

polariseur 2

d2

polariseur 2
lumire
polarise

direction d 1
d

Sortie :
lumire

polarise

car d 1 k d 2

Sortie :
extinction

car d 1 d 2

Figure VIII.9 Filtres polarisants.

Birfringence accidentelle
Une lumire plane se prsentant suivant une direction de polarisation quelconque par rapport

aux axes dun corps birfringent se dcompose en deux composantes parallles b 1 et b 2


ces axes (Fig. VIII.10), chacune dentre elles se comportant comme une onde plane autonome
progressant une vitesse propre sa direction.
La plupart des corps transparents isotropes deviennent birfringents lorsquils sont soumis des contraintes ; cette birfringence accidentelle est telle que les axes de birfringence
coincident avec les directions principales des contraintes.
En plaant le milieu birfringent entre deux polariseurs croiss (Fig. VIII.10), on observe
alors une extinction de lumire lorsque les axes de birfringence sont parallles aux axes des
polariseurs : ce sont les isoclines.
De plus, chaque onde se propageant dans le milieu birfringent suivant chacune des

directions de birfrince b 1 et b 2 se propage avec une vitesse diffrente (Fig. VIII.11).

Londe suivant b 1 se propage la vitesse c1 et celle suivant b 2 se propage la vitesse c2 .


Les longueurs donde 1 et 2 des deux ondes sont diffrentes dans le milieu birfringent,
mais sont identiques dans lair aprs la traverse du milieu. Ce retard induit une extinction
de la lumire telle que (loi de Maxwell) :
1 2 = N

ce

o 1 et 2 sont les contraintes principales, N est lordre de frange, est la longueur donde,
c est la vitesse de la lumire et e est lpaisseur du milieu birfringent.
Le lieu des points o la lumire est teinte due au retard sont les isochromatiques. Lordre

76

Moyens exprimentaux

polariseur 1
lumire quelconque

milieu birfringent

b2

b1

d1

polariseur 1

polariseur 2

d2

milieu birfringent

lumire quelconque

polariseur 2

b2

b1

d1

Sortie :
lumire polarise

d2

Sortie :
extinction

b 1 est
direction principale

Figure VIII.10 Phnomne de birfringence accidentelle.

effet de la traverse de la lame dans la direction b 1


lumire dans le vide
Milieu
birfringent

retard

effet de la traverse de la lame dans la direction b 2


Figure VIII.11 Diffrence de phase entre les deux ondes qui sortent du milieu birfringent.

77

VIII.2 Photolasticit

de frange N est la composante de la lumire teinte en lumire blanche. Il y a extinction


totale en lumire monochromatique. Lobservation des couleurs permet de dterminer N .
Ordre N

couleur visible

passage rouge-bleu

passage rouge-vert

passage rouge-vert

Dans la pratique, leffet des isoclines perturbe lobservation des isochromatiques. On


ajoute alors dans le montage deux lames quart-donde pour faire disparaitre les isoclines
(lexplication du fonctionnement des lames quart-donde nest pas donne ici).

VIII.2.2 Mise en quation

En sortie du polariseur 1 le vecteur lumineux est parallle d 1 =


x :

2 
z

L =a
x sin
t .
T
c

Le vecteur
x se dcompose suivant b1 et b2 comme :

x = cos b1 sin b2 .
Donc le vecteur lumire se prsentant lentre du milieu birfringent vaut :

z
2 
z
2 
t
a sin b2 sin
t .
L = a cos b1 sin
T
c
T
c

Pendant la traverse du milieu birfringent dpaisseur e, la lumire se propage la vitesse

c1 suivant b1 et c2 suivant b2 . Les temps de traverse suivant chaque direction valent :


t1 =

e
e
et t2 = .
c1
c2

Le retard par rapport au temps quaurait mis chaque onde pour traverser le mileu est :
1 = c(t1 tvide ) = c
et

e

c1

2 = c(t2 tvide ) = e


c
e
1
=e
c
c1
c

1 .

c2
On appelle n1 et n2 les indices du milieu birfringent tels que :
n1 =

c
c
et n2 = .
c1
c2

Les retards 1 et 2 deviennent alors :


1 = e (n1 1) et 2 = e (n2 1).
A la sortie du milieu birfringent, le vecteur lumire a pour expression :

soit :

2 
z 1 
2 
z 2 
L = a cos b1 sin
t
a sin b2 sin
t
,
T
c
c
T
c
c

z + (n1 1)e 
2 
z + (n2 1)e 
2 
t
a sin b2 sin
t
.
L = a cos b1 sin
T
c
T
c

78

Moyens exprimentaux

En effectuant le changement de variable :


z = z + (n1 1)e
on peut crire :
z + (n2 1)e = z + (n1 1)e + (n2 1)e (n1 1)e = z + (n2 n1 )e = z
avec = e(n1 n2 ) qui est le retard entre les deux composantes du vecteur lumire
la sortie du milieu birfringent. En remplacant dans lexpression du vecteur lumire, son
expression devient :

2 
z 
2 
z 
L = a cos b1 sin
t
a sin b2 sin
t
T
c
T
c

Lanalyseur (ou polariseur 2) a pour direction de polarisation y . Le vecteur lumire la


sortie de lanalyseur vaut :
z 
2 
z 

2 

t
a sin cos
y sin
t
L = a cos sin
y sin
T
c
T
c
car

b1 .
y = sin et b2 .
y = cos

En simplifiant, on a :
h

2 
z 
2 
z i

L = a cos sin sin


y
t
sin
t
T
c
T
c

avec :
sin p sin q = 2 cos

pq
p+q
sin
2
2

on a :

2 

z
2  

L = a sin 2 cos
y
t +
sin
T
c
2c
T 2c
Lamplitude de sortie du vecteur lumire vaut :
A = a sin 2 sin

T c

En introduisant la longueur donde = c T , on a :


A = a sin 2 sin

Cette amplitude vaut zro dans deux cas diffrents :


sin 2 = 0 ou

sin

= 0.

Le premier cas sin 2 = 0 correspond aux isoclines, en effet ceci est quivalent = 0
ou = 2 .
Le deuxime cas correspond aux isochromatiques :
sin

=0
= N = N

79

VIII.2 Photolasticit

or la loi de Maxwell tant :


n1 n2 = c(1 2 )
et = (n1 n2 )e, on a :

N
ce
Dans le cas dune lumire monchromatique, on observe une extinction de lumire (bande
1 2 =

noire). Dans le cas dune lumire blanche, seule la radiation correspondante


c e (1 2 )
=N

est teinte.

VIII.2.3 Rseaux de courbes caractristiques


Isoclines
Lieu des points du plan ou les contraintes principales ont une direction constante, chaque
isocline est accompagne dune cote donnant cette orientation.
Exemple dapplication :
Sur lisocline dessine pour langle 30o , les directions principales des contraintes pour
tous les points de cette isocline sont 30o et 30 + 90 = 120o .
Isochromatiques
Lieu des points du plan pour lesquels la diffrence des contraintes principales est constante,
proportionnelle N : 1 2 = kN .
Exemple dapplication :
On tudie une pice rectangulaire troue en son milieu, sur laquelle on applique un effort
de traction (fig. VIII.12). On peut considrer que loin du trou, ltat de contrainte est celui
dune pice en traction. Cela entraine que 2 = 0 et donc par consquent 1 = kN . Si de
plus, on place le premier passage rouge-bleu (N = 1) dans cette partie de la pice, alors
1 = k, ou encore k = o est la contrainte de traction applique loin du trou.

Le long du trou, au point A, les directions principales des contraintes sont


x et
y . La
contrainte principale yy est nulle car aucun effort extrieur nest appliqu en ce point. On
a donc en ce point :
xx = kN = N
En connaissant le numro de lisochromatique passant par ce point, on connait N , et par
consquent la valeur de xx . Si le point A est situ entre les isochromatiques 4 et 5, on en
dduit que
4 < xx < 5
ou encore

xx
< 5.

Le coefficient de concentration de contraintes cause de la prsence du trou est compris


4<

entre 4 et 5.

80

Moyens exprimentaux

Isochromatique N=1

loin du trou, la pice est en tat de traction


Figure VIII.12 Plaque troue en "traction".

Isostatiques
Les isostatiques sont les courbes donnant en chaque point du plan, par leur tangente et leur
normale en ce point, la direction des contraintes principales. Ce rseau peut se construire graphiquement partir du rseau des isoclines. La mthodologie est la suivante (Fig. VIII.13) :
fixer le trac des isoclines sous une feuille de papier calque,
tracer le contour de la pice tudie,
placer une feuille de papier avec un quadrillage sous lensemble prcdent,
orienter le quadrillage pour quil soit parallle la direction 0o (tous les points de
lisocline 0o ont donc pour directions principales les axes du quadrillage),
tracer le long de lisocline 0o des petites croix parallles et perpendiculaires au quadrillage (par exemple avec un intervalle de 1 cm entre deux croix),
orienter le quadrillage pour quil soit parallle la direction 15o ,
tracer le long de lisocline 15o des petites croix parallles et perpendiculaires au quadrillage,
recommencer ces oprations pour les directions 30o , 45o , 60o et 75o ,
enlever le trac des isoclines de faon avoir seulement les croix et le contour de la
pice,
placer une nouvelle feuille de papier calque sur les croix pour pouvoir effacer les isostatiques sans effacer les croix,
tracer les isostatiques en suivant les trois rgles :
les isostatiques sont tangentes aux croix,
les isostatiques sont parallles et perpendiculaires entre elles,
le bord de la pice est une isostatique,
garder seulement le trac des isostatiques et vrifier que les trois rgles sont respectes.

VIII.2 Photolasticit

isoclines

croix

trac des isostatiques

isostatiques

Figure VIII.13 Etapes de construction des isostatiques.

81

Vous aimerez peut-être aussi