Rhologie

Georges Cailletaud
Centre des Matriaux
MINES ParisTech/CNRS

Plan

Plan

Essais mcaniques
Structures
Elments de volume

Modles rhologiques
Les briques de base
Plasticit
Viscolasticit
Elastoviscoplasticit

Bilan

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

2 / 52

Plan

Essais mcaniques
Structures
Elments de volume

Modles rhologiques
Les briques de base
Plasticit
Viscolasticit
Elastoviscoplasticit

Bilan

Essais mcaniques

Structures

Tests dun avion civil

[Link]

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

4 / 52

Essais mcaniques

Structures

Vibration dune aile

[Link]

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

5 / 52

Essais mcaniques

Structures

Structures biologiques

[Link]

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

6 / 52

Essais mcaniques

Structures

Structures biologiques

[Link]

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

7 / 52

Essais mcaniques

Structures

Alimentaire

[Link]

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

8 / 52

Plan

Essais mcaniques
Structures
Elments de volume

Modles rhologiques
Les briques de base
Plasticit
Viscolasticit
Elastoviscoplasticit

Bilan

Essais mcaniques

Elments de volume

Machines dessai

[Link]
Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

10 / 52

Essais mcaniques

Elments de volume

Compression du gypse

Doc. Mines Paris-Gosciences

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

Plus de dtails sur site [Link]

8 mars 2010

11 / 52

Essais mcaniques

Elments de volume

Rsultat de compression du gypse

Doc. Mines Paris-Gosciences

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

12 / 52

Essais mcaniques

Elments de volume

Traction sur fibre

[Link]

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

13 / 52

Essais mcaniques

Elments de volume

Traction sur alliage mtallique

[Link]

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

14 / 52

Essais mcaniques

Elments de volume

Essais mcaniques
Les essais de base
Matriaux dont le comportement est insensible la vitesse de sollicitation
Essai de traction, ou essai dcrouissage
Essai sous chargement cyclique, ou essai de fatigue

Matriaux dont le comportement est sensible la vitesse de sollicitation

Essai contrainte constante, ou essai de fluage
Essai dformation constante, ou de relaxation

Autres essais
Essais sous chargement multiaxial
Tractiontorsion
Pression interne ou externe

Essais en flexion
Essais de fissuration

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

15 / 52

Essais mcaniques

Elments de volume

Rsultat de traction sur un alliage daluminium

Domaine dlasticit vrai, limite dlasticit conventionnelle, 0.2 , qui donne
0.2% de dformation rsiduelle la dcharge
Contrainte ultime, u
600
500
(MPa)

400
300
200
Tension curve
Elastic slope
0.2% residual strain

100
0
0

0.01

0.02

0.03

0.04

(mm/mm)

E=78000 MPa, 0.2 =430 MPa, u =520 MPa

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

Doc. Mines Paris-CDM, Evry

8 mars 2010

16 / 52

Essais mcaniques

Elments de volume

Rsultat de traction sur un acier inoxydable

Matriau prsentant un crouissage important : possibilit de durcissement dans
le domaine plastique, augmentation de la limite dlasticit courante
600
500
(MPa)

400
300
200
Tension curve
Elastic slope
0.2% residual strain

100
0
0

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08

(mm/mm)

E=210000 MPa, 0.2 =180 MPa, u =660 MPa

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

Doc. ONERA-DMSE, Chtillon

8 mars 2010

17 / 52

Essais mcaniques

Elments de volume

Traction et compression sur un alliage daluminium

Essai dformation impose symtrique 0.3%
A contrainte nulle, il reste une dformation positive
A dformation nulle, la contrainte est devenue ngative
300
200
(MPa)

100
0
-100
-200
-300
-0.005

-0.003

-0.001

0.001

0.003

0.005

(mm/mm)
Doc. Mines Paris-CDM, Evry
Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

18 / 52

Essais mcaniques

Elments de volume

Modles schmatisant les rsultats prcdents

ET

a. lastiqueparfaitement plastique

b. lastiqueplastique linaire

Module lastoplastique, ET = d /d .
ET = 0 : matriau lastique-parfaitement plastique
ET constant : matriau lasto-plastique linaire
Et fonction de la dformation dans le cas gnral

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

19 / 52

Essais mcaniques

Elments de volume

Fonctionnement dun modle de plasticit instantane

Rgime lastique
OA, OB
Ecoulement plastique
AB
Dformation rsiduelle
OO

Dcomposition de la dformation, = e + p ;
Domaine dlasticit, dfinir par une fonction de charge f
Ecrouissage, dfinir par des variables dcrouissage, AI .

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

20 / 52

Essais mcaniques

Elments de volume

Rsultat de traction sur un acier haute temprature

Effet de la viscosit : Comportement sensible la vitesse de dformation

725 C

(MPa)

80
60
40
= 2.4 104 s1
= 8.0 105 s1
= 1.6 105 s1

20
0
0

0.02

0.04

0.06

0.08

0.1

Doc. Ecole des Mines, Nancy

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

21 / 52

Essais mcaniques

Elments de volume

Essai de fluage sur fil tain-plomb

Voir lexercice
Mines Paris-CDM, Evry

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

22 / 52

Essais mcaniques

Elments de volume

Rsultat de fluage sur une fonte (1)

0.03
0.025

0.02
0.015
=25MPa
=20MPa
=16MPa
=12MPa

0.01
0.005
0
0

200

400

600

800

1000

t (s)
Doc. Mines Paris-CDM, Evry
Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

23 / 52

Essais mcaniques

Elments de volume

Reprsentation schmatique dune courbe de fluage

Fluage primaire , priode de durcissement du matriau (crouissage)
Fluage secondaire , ou stabilis : p est une fonction puissance de la
contrainte applique
Fluage tertiaire , ou perte de rsistance conduisant la rupture, dcrite
par des variables dendommagement

III
II
I

t
Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

24 / 52

Essais mcaniques

Elments de volume

Rsultat de fluage sur une fonte (2)

T=500 C
T=600 C
T=700 C
T=800 C

0.001

p (s1)

0.0001
1e-05
1e-06
1e-07
1e-08
1

10

100

(MPa)
Doc. Mines Paris-CDM, Evry
Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

25 / 52

Essais mcaniques

Elments de volume

Essai de relaxation

On impose une dformation constante en fonction du temps

Pendant lessai :

E
= 0 = p + /
d p = d /E
Pour une dformation positive, la dformation viscoplastique augmente pendant
que la contrainte diminue
La valeur asymptotique de la contrainte est nulle (relaxation totale) ou non
(relaxation partielle)
Relaxation partielle sil existe une contrainte interne ou contrainte seuil dans le
matriau,

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

26 / 52

Essais mcaniques

Elments de volume

Reprsentation schmatique dune courbe de relaxation

Le point reprsentatif est obtenu comme la somme de la contrainte seuil s et de
la contrainte visqueuse v
La contrainte seuil reprsente un comportement plastique qui peut tre atteint
asymptotiquement lorsque la vitesse tend vers zro

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

p
Rhologie

8 mars 2010

27 / 52

Plan

Essais mcaniques
Structures
Elments de volume

Modles rhologiques
Les briques de base
Plasticit
Viscolasticit
Elastoviscoplasticit

Bilan

Modles rhologiques

Les briques de base

Les briques de base pour les modles de matriau

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

29 / 52

Modles rhologiques

Les briques de base

Diffrents types de rhologies

Plasticit indpendante du temps

= e + p

d p = f (...)d

= e + p

d p = f (...)dt

Elasto-viscoplasticit
Viscolasticit

, ) = 0
F (, ,

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

30 / 52

Plan

Essais mcaniques
Structures
Elments de volume

Modles rhologiques
Les briques de base
Plasticit
Viscolasticit
Elastoviscoplasticit

Bilan

Modles rhologiques

Plasticit

Plasticit indpendante du temps

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

32 / 52

Modles rhologiques

Plasticit

Modle lastiqueparfaitement plastique

Le rgime de fonctionnement est dfini par la
fonction de charge f (de lespace des contraintes dans R)
f () = || y

Domaine dlasticit
si f < 0

E
= e = /

Dcharge lastique
si f = 0 et f < 0

E
= e = /

Ecoulement plastique
si f = 0 et f = 0

= p

La condition f = 0 est la condition de cohrence

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

33 / 52

Modles rhologiques

Plasticit

Modle de Prager
Fonction de charge deux variables, et X
f (, X ) = | X | y

avec X = H p

Ecoulement plastique si on vrifie la fois f = 0 et f = 0.

f
f
+
X =0

X
signe( X ) signe( X ) X = 0

soit : = X

Vitesse de dformation plastique fonction de la vitesse de contrainte

H
p = /
Vitesse de dformation plastique fonction de la vitesse de dformation

p =

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

E
E +H

Rhologie

8 mars 2010

34 / 52

Modles rhologiques

Plasticit

Ecriture des quations de llastoplasticit uniaxiale

Domaine dlasticit
si f (, Ai ) < 0

E
= /

Dcharge lastique
si f (, Ai ) = 0 et f (, Ai ) < 0

E
= /

Ecoulement plastique
si f (, Ai ) = 0 et f (, Ai ) = 0

E + p
= /

La condition de cohrence scrit :

f (, Ai ) = 0

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

35 / 52

Modles rhologiques

Plasticit

Illustration des deux types dcrouissage

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

36 / 52

Modles rhologiques

Plasticit

Modle dcrouissage isotrope

Fonction de charge deux variables, et R
f (, R ) = || R y
R dpend de p, dformation plastique cumule : p = | p |
dR /dp = H

soit R = H p

Ecoulement plastique ssi f = 0 et f = 0

f
f
+
R=0

R
sign() R = 0

soit sign() H p

Vitesse de dformation plastique fonction de la vitesse de contrainte

H
p = sign() /

H
soit p = /

Modles classiques
Ramberg-Osgood : = y + Kpm
Loi exponentielle : = u + (y u ) exp(bp)
Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

37 / 52

Plan

Essais mcaniques
Structures
Elments de volume

Modles rhologiques
Les briques de base
Plasticit
Viscolasticit
Elastoviscoplasticit

Bilan

Modles rhologiques

Viscolasticit

Viscolasticit

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

39 / 52

Modles rhologiques

Viscolasticit

Rponses lmentaires en viscolasticit

E0 + /
Elments en srie, modle de Maxwell : = /
Fluage sous une contrainte 0 :

= 0 /E0 + 0 t /

Relaxation la dformation 0 :

= E0 0 exp[t /]

Elments en parallles, modle de Voigt : = H + ou = ( H )/

Fluage sous une contrainte 0 : = (0 / H )(1 exp[t /0 ])
Les constantes = /E0 et 0 = /H sont homognes un temps,
dsignant le temps de relaxation du modle de Maxwell

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

40 / 52

Modles rhologiques

Viscolasticit

Modles composs
(E1 )

(H)
(E0 )
()

(E2 )

a. KelvinVoigt

()

b. Zener

Rponses en fluage et en relaxation


(1 exp[t /f ]) 0
E0 H


H
E0
(t ) = E (t ) 0 =
+
exp[t /r ] E0 0
H + E0 H + E0
(t ) = C (t ) 0 =

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

41 / 52

Plan

Essais mcaniques
Structures
Elments de volume

Modles rhologiques
Les briques de base
Plasticit
Viscolasticit
Elastoviscoplasticit

Bilan

Modles rhologiques

Elastoviscoplasticit

Elasto-viscoplasticit

Schma du modle
X = H vp

Rponse en traction

v = vp

|p | 6 y

= X + v + p
Domaine dlasticit, dont la frontire est |p | = y
Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

43 / 52

Modles rhologiques

Elastoviscoplasticit

Equations du modle
Trois rgimes de fonctionnement

(a) vp = 0
(b) vp > 0
(c ) vp < 0

|p | = | H vp |

6 y

p = H vp vp = y
p = H vp vp = y

(a) intrieur ou frontire du domaine dlasticit (|p | < y )

(b), (c ) coulement (|p | = y et | p | = 0 )
On peut rsumer les trois quations (en posant < x >= max (x , 0)) par

vp = h| X | y i signe( X )
ou :

vp =

<f >
signe( X ) ,

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

avec f (, X ) = | X | y

8 mars 2010

44 / 52

Modles rhologiques

Elastoviscoplasticit

Fluage avec un modle de Bingham

vp

o - y

t
Dformation viscoplastique
en fonction du temps

vp =

o y

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

vp

Evolution dans le plan

contraintedformation
viscoplastique

1 exp

Rhologie


avec : f = /H

8 mars 2010

45 / 52

Modles rhologiques

Elastoviscoplasticit

Relaxation avec un modle de Bingham

Transitoire : OA = BC

A
H

-E

Relaxation : AB

H
vp

vp

Effacement

Relaxation

= y

E
E +H

1 exp


+

avec : r =

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Effacement
incomplet : CD

Rhologie

E o
E +H

H + E exp



E +H

8 mars 2010

46 / 52

Modles rhologiques

Elastoviscoplasticit

Ingrdients des modles classiques de viscoplasticit

Modle de Bingham

vp =

<f >
signe( X )

Plus gnralement

vp = (f )
(0) = 0

et monotone croissante

vp est nulle si le point courant se trouve dans le domaine dlasticit ou sur le

bord de celui-ci

vp est non nulle si le point courant se trouve lextrieur du domaine dlasticit

On distingue des modles avec/sans seuil et avec/sans crouissage

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

47 / 52

Modles rhologiques

Elastoviscoplasticit

Modles viscoplastiques sans crouissage

Modles sans seuil : le domaine dlasticit peut se rduire lorigine ( = 0)

Modle de Norton

vp =

||

n
signe()

Modle de Sellars-Tegart

vp = A sh

||
K


signe()

Modles seuil
Modle de Perzyna

vp =

|| y

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

n
signe()

Rhologie

vp = 0

||
1
y

n
signe()

8 mars 2010

48 / 52

Modles rhologiques

Elastoviscoplasticit

Modles viscoplastiques avec crouissage

Notion dcrouissage additif : le durcissement provient des variables qui
expriment le seuil (X et R)

vp =

| X | R y

n

signe( X )

X dsigne la contrainte interne, internal stress (crouissage cinmatique)

R + y dsigne la contrainte de friction, friction stress (crouissage isotrope)
v est la contrainte visqueuse, drag stress

Notion dcrouissage multiplicatif : on fait varier la contrainte visqueuse, par

exemple :





vp =

||
K (p )

signe() =

||
K0 |p |m

signe()

(crouissage par la dformation, ou strain hardening)

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

49 / 52

Bilan

En plasticit et en viscoplasticit...
Domaine dlasticit dfini par une fonction de charge f < 0
Variables dcrouissage isotrope et cinmatique
En plasticit :
Ecoulement dfini par la condition de cohrence si f = 0, f = 0
Ecoulement plastique instantan :
d p = g (, . . . )d
En viscoplasticit :
Ecoulement dfini par la fonction de viscosit si f > 0
Possibilit dcrouissage sur la contrainte visqueuse
Ecoulement viscoplastique retard
d vp = g (, . . . )dt

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

50 / 52

Bilan

Identification des essais sur le fil de brasure

0.1
exp
sim

14

0.08

12

stress (MPa)

creep strain

10
0.06

0.04

4
1534 g
1320 g
1150 g
997 g
720 g

0.02

0
0

1000

2000

3000

4000

5000

time (s)

5000

10000

15000

20000

25000

time (s)

Essais de fluage

Relaxation =20%

Courbes obtenues avec un modle de Norton

p =

 2.3
800

Jessaie tout(e) seul(e) sur le site [Link] O

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

51 / 52

Bilan

Identification du fluage du sel gemme

0.008
0.007
0.006

strain

0.005
0.004
0.003
0.002
0.001
exp
sim
0
0

Eprouvette

0.5

1.5

2
time (Ms)

2.5

3.5

Essai 3 niveaux (3, 6, 9 MPa)

Courbes obtenues avec un modle de Lemaitre (strain hardening)

p =

 n
K

(p + v0 )m

Jessaie tout(e) seul(e) sur le site [Link] O

Georges Cailletaud (Centre des Matriaux/UMR 7633 )

Rhologie

8 mars 2010

52 / 52