Rhéologie des fluides complexes

Elisabeth Guazzelli

To cite this version:

Elisabeth Guazzelli. Rhéologie des fluides complexes. École d’ingénieur. Rhéologie des fluides com-
plexes, France. 2001. �cel-01522165�

HAL Id: cel-01522165

https://cel.hal.science/cel-01522165
Submitted on 13 May 2017

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est

archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents
entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non,
lished or not. The documents may come from émanant des établissements d’enseignement et de
teaching and research institutions in France or recherche français ou étrangers, des laboratoires
abroad, or from public or private research centers. publics ou privés.
 Rhéologie des fluides complexes

Élisabeth Guazzelli
Octobre 2001

Ce cours présente les notions de base de la rhéologie des fluides non-newtoniens. J’ai voulu don-
ner une présentation des phénomènes nouveaux rencontrés ainsi que leur description quantitative.
Je me suis néanmoins limitée à des situations simples en choisissant en particulier le cisaillement
simple comme écoulement de base pour décrire les lois constitutives des matériaux, ce qui offre un
cadre mathématique simplifié et facilite l’interprétation physique. J’ai donné également quelques
principes de rhéométrie et quelques calculs d’écoulements simples. Par contre, je n’ai pas voulu
fournir une analyse structurelle microscopique qui permettrait de déduire les lois rhéologique
de certains de ces liquides, car cela me semblait bien au-delà des perspectives de ce cours qui
n’est qu’introductif. Enfin, bien qu’une présentation des contraintes et déformations soit faite en
début, ce cours nécessite quelques connaissances préalables en mécanique des fluides et/ou en
mécanique des milieux continus. Quelques bons livres de base et la référence d’un film sur la
rhéologie sont donnés à la fin du cours. Remarques, corrections et suggestions sont les bienvenues
à [email protected].

Table des matières

1 Introduction 3

2 Phénoménologie 3
2.1 Comportements non linéaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2 Comportement dépendant du temps . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.3 Résistance à l’étirement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.4 Effets élastiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.5 Contraintes normales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

3 Contrainte et déformation 9
3.1 Contrainte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3.2 Déformation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.3 Décomposition du champ de gradient de vitesse au voisinage d’un point . . . . . . 11

4 Lois constitutives 11
4.1 Fluides newtoniens et non-newtoniens . . . . . . . . . . . . . . . . . . . . . . . . . 11
4.2 Viscoélasticité linéaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
4.2.1 Fonctions fluage et relaxation . . . . . . . . . . . . . . . . . . . . . . . . . . 13
4.2.2 Solide élastique parfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
4.2.3 Liquide visqueux newtonien . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
4.2.4 Solide de Kelvin-Voigt . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
4.2.5 Liquide de Maxwell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
4.2.6 Comportement viscoélastique linéaire général . . . . . . . . . . . . . . . . . 18
4.3 Liquides nonlinéaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.3.1 Liquides plastiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.3.2 Fluide rhéofluidifiant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21

1
 4.3.3 Fluide rhéoépaississant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.4 Comportement dépendant du temps . . . . . . . . . . . . . . . . . . . . . . . . . . 23

5 Rhéomètrie 24
5.1 Rhéomètres simples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
5.1.1 Rhéomètre plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
5.1.2 Rhéomètre de Couette . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
5.1.3 Rhéomètre cône-plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
5.1.4 Autres types de rhéomètres . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.2 Cisaillement oscillant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.3 Contraintes normales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
5.4 Viscosité élongationnelle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29

A Coordonnées cartésiennes, cylindriques et sphériques 30

B Dimensions et unités 31

C Exercices et problèmes 32
C.1 Course entre un fluide newtonien et un fluide non-newtonien . . . . . . . . . . . . 32
C.2 Loi rhéologique à trouver . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
C.3 Écoulement dans un tube cylindrique d’un fluide en loi de puissance . . . . . . . . 32
C.4 Écoulement de boues de forage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
C.5 Effet Weissenberg . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35

2
1 Introduction
⊂
La rhéologie, mot inventé par Bingham en 1929 à partir du verbe grec ρ ǫ́ω qui veut dire couler
(voir l’historique sur la figure 1), est l’étude des écoulements et des déformations. Devant l’impuis-
sance de la théorie de l’élasticité et de la mécanique des fluides (théories élaborées au 19e siècle) à
décrire et à expliquer les propriétés de matériaux aux comportements mal définis et intermédiaires
entre celui du solide élastique parfait (où les contraintes sont proportionnelles aux déformations) et
celui du fluide newtonien (où les contraintes sont proportionnelles aux vitesses de déformation), il
est apparu nécessaire d’élaborer cette nouvelle discipline. Les études expérimentales s’attachent à
mesurer les propriétés de l’écoulement des matériaux tandis que les approches théoriques cherchent
les équations constitutives reliant contraintes et déformations.

2 Phénoménologie
2.1 Comportements non linéaires
Un des comportements pratiques les plus intéressants des fluides non-newtoniens est leur re-
lation non linéaire entre contrainte et vitesse de déformation. La structure interne du fluide est
complexe et peut être influencée par l’écoulement.
Certains fluides ne s’écoulent qu’à partir d’une certaine contrainte seuil (liquide plastique de
Bingham, comme sur la figure 2). Cette propriété est particulièrement utile pour le transport de
particules en empêchant la sédimentation (boues de forage) et se rencontre dans la vie pratique
dans les pâtes dentifrices, le ketchup, la graisse et les peintures non-coulantes.
La plupart des fluides sont rhéofluidifiants : suspensions diluées de particules solides, suspen-
sions de vésicules déformables (comme le sang), encres, peintures, solutions diluées de polymères,
polymères liquides (acétate de cellulose), pâte à papier (voir la figure 2). Leur viscosité effec-
tive diminue lorsqu’on augmente la contrainte. Cet effet est dû en général à une brisure de la
structure interne par l’écoulement. Quelques fluides sont rhéoépaississants comme les suspensions
concentrées ou encore le sable mouillé (voir encore la figure 2).

2.2 Comportement dépendant du temps

Les fluides thixotropes ont une viscosité effective qui diminue avec le temps quand une contrainte
constante est appliquée (par exemple, le ketchup, le yoghourt, certaines peintures). Cela s’explique
par des changements de structures intervenant dans le fluide avec des temps caractéristiques
comparables aux temps d’observation. Comme la réponse du fluide dépend de son histoire, cela
complique les mesures rhéologiques car il existe des effets d’hystérésis (effets de mémoire). Il existe
aussi, mais plus rarement, des fluides antithixotropes (voir un exemple sur la figure 3).

2.3 Résistance à l’étirement

La résistance bien plus forte à un étirement du fluide que celle à laquelle on s’attendrait
venant d’une viscosité de cisaillement est attribuée à une viscosité élongationnelle très grande
bien qu’elle soit peut-être due, en fait, à un effet élastique. Cet effet rend les bulles d’air pointues à
l’arrière lorsqu’elles montent dans un shampooing : la pointe convertit l’écoulement élongationnel
en écoulement de cisaillement (voir la figure 4). Ce phénomène réduit la pulvérisation de fines
gouttes hors des jets avec les applications qui en découlent : maintenir les jets plus concentrés dans
les lances à incendie et réduire la tendance à faire une brume explosive pour le fuel des avions
(voir la figure 5). Cela aussi explique la réduction de la traı̂née dans les écoulements turbulents
résultant de l’addition d’une très faible quantité d’un polymère de très grande masse moléculaire
au fluide en écoulement.

3
Fig. 1 – Rhéologie et nombre de Deborah, article de M. Reiner, Physics Today, janvier 1964,
page 62.

4
 

Contrainte

6

 (2)
 
 (4) 

  (1)
 
  


 (3)



 -
Taux de cisaillement

Fig. 2 – Fluide (1) newtonien, (2) de Bingham, (3) rhéofluidifiant, (4) rhéoépaississant.

Fig. 3 – Fluide antithixotrope : (a) Au repos ; (b) Après une forte agitation (photographies tirées
de [2]).

2.4 Effets élastiques

On trouve de nombreux fluides élastiques dans l’industrie alimentaire. Une démonstration des
effets élastiques dans les liquides est l’expérience du siphon ouvert dans laquelle on peut siphonner
un liquide élastique d’un bécher à un autre sans utiliser un tuyau (voir la figure 6). Cela est dû
à une forte viscosité élongationnelle qui induit une contrainte très forte au niveau du bec du
bécher qui tire le fluide. Un autre expérience est celle de la détente élastique (voir la figure 7).
Un matériau, que l’on verse d’un récipient, peut être coupé à l’aide de ciseaux et remonte dans le
récipient.

2.5 Contraintes normales

Sous l’effet d’un cisaillement simple, l’élasticité des liquides se manifeste par l’apparition de
contraintes normales : il peut y avoir une tension le long des lignes de courant et une pression per-
pendiculaire à celles-ci. Une démonstration spectaculaire de ce phénomène est l’effet Weissenberg
(voir les figures 8 et 9). Le fluide élastique remonte le long de l’axe d’un batteur rotatif alors qu’il
y a creusement de la surface libre dans le cas d’un fluide newtonien à cause des forces centrifuges.
Cela est dû à une tension le long des lignes de courant qui augmente lorsqu’on se rapproche du
batteur et qui provoque la remontée du fluide le long de l’axe. Une autre expérience frappante
de ce phénomène est celle du gonflement à l’extrusion d’un tube (voir les figures 10 et 11). À la
sortie de l’orifice du tube, la tension le long des lignes de courant se relâche en se contractant
dans la direction longitudinale, ce qui résulte en une expansion latérale du liquide à cause de son
incompressibilité.

5
Fig. 4 – Bulles d’air ascendantes dans un tube contenant différents fluides : (a) Huile newtonienne ;
(b) Fluide de Boger ; (c) Solution de 6% de CMC ; (d) Solution de 3% de Polyox (photographies
tirées de [2]).

Fig. 5 – Jet à nombre de Reynolds élevé (photographies tirées de [2]) : Eau (photo du haut) ;
Solution aqueuse de 200 ppm d’oxyde de polyéthylène (photo du bas).

6
Fig. 6 – Siphon ouvert réalisé avec une solution aqueuse de 0,75% de Polyox WSR 301 (photogra-
phies tirées de [2]). La séquence montre le développement du siphon à partir du versement initial
hors du bécher.

Fig. 7 – Détente élastique (photographies tirées de [2]).

7
 Fig. 8 – Effet Weissenberg (photographie tirée de [2]).

Fig. 9 – Instabilité de la remontée du liquide (photographie tirée de [2]).

Fig. 10 – Gonflement à l’extrusion pour une solution aqueuse de 2% de polyacrylamide (photo-

graphie tirée de [2]).

8
Fig. 11 – Gonflement retardé à l’extrusion pour une solution aqueuse de 5% de polyacrylamide
lorsque le nombre de Reynolds augmente, montrant la compétition entre élasticité et inertie (pho-
tographies tirée de [2]).

3 Contrainte et déformation
3.1 Contrainte
Il est possible de classer les forces qui s’exercent sur un élément de fluide de volume dV en
deux catégories selon leur portées :
– Les forces en volume qui sont des forces à longues portées (gravité, forces électriques,
magnétiques). Par exemple, la gravité s’exerce sur des distances très grandes devant les
dimensions moléculaires et de fait la force de gravité sur un élément de volume est propor-
tionnelle à son volume.
– Les forces de surfaces à courtes portées dont l’origine est dans les interactions moléculaires
qui assurent la cohésion du liquide. Ces interactions à courtes portées ne vont concerner
qu’une mince couche externe sur un élément de volume donné. La force globale exercée par
ces interactions est proportionnelle à l’aire de la surface limitant l’élément de fluide et est
indépendante de son volume.
Sur un élément de surface dS, de normale unité n, cette force de surface va s’écrire :
dF = σ · ndS. (1)
Le tenseur des contraintes σ (ou σij ) est symétrique et donc σij = σji 1 . Les termes diagonaux σii
sont les contraintes normales tandis que les termes non diagonaux sont les contraintes tangentielles.
Dans un fluide au repos, le tenseur des contraintes est 2 :
σij = −pδij , (2)
3
où p est la pression. Le tenseur des contraintes est alors isotrope :
σ11 = σ22 = σ33 = −p. (3)
Dans un fluide en mouvement, le tenseur des contraintes s’écrit :
σij = −pδij + dij , (4)
où le premier terme qui est la partie isotrope continue à s’appeler la pression tandis que le second
terme, la partie anisotrope due à la viscosité du fluide, est appelé le déviateur.
1. Ici et par la suite, nous utilisons une notation indicielle ainsi que la convention de sommation d’Einstein (un
indice répété implique une sommation, par exemple ui ui = u21 + u22 + u23 ).
2. Le symbole de Kronecker δij est défini par

0 si i 6= j
δij =
1 si i = j.

3. Il y a un signe négatif car le fluide est en compression.

9
 U
- 6σyy

6 σxy
y -
6 -
-x -
σxx
h z

σzz

Fig. 12 – Cisaillement simple.

-
α
uM
uM
y
′
e-
6

uM0 u
M0′
z x

Fig. 13 – Déformation.

Considérons le cas où les contraintes sont appliquées par l’intermédiaire d’un plan animé d’une
vitesse U et d’un plan parallèle fixe distants de h comme indiqué sur la figure 12. Les contraintes
ne dépendent que de y par symétrie. De plus seules existent la contrainte de cisaillement σ =
σxy = σyx (σxz = σyz = 0) et les contraintes normales σxx , σyy et σzz . C’est cette configuration
de cisaillement simple que nous utiliserons par la suite.

3.2 Déformation
Sous l’action des contraintes, un matériau va se déformer. On considère l’écoulement de cisaille-
ment laminaire de la figure 12. Si un élément matériel est localisé au point M0 de coordonnées
(x,y,z), il sera localisé au point M0′ de coordonnées (x + ξ,y,z) à l’instant t puisque la vitesse
du fluide est parallèle à l’axe des x (voir le schéma de la figure 13). Le déplacement ξ dépend
seulement de y. L’élément matériel situé au point M de coordonnées (x,y + dy,z), sera localisé
au point M ′ de coordonnées [x + ξ + (dξ/dy)dy,y + dy,z] à l’instant t. On appelle déformation, la
variation du déplacement lorsqu’on passe d’une couche à une couche infiniment voisine :
dξ
γ= = tan α. (5)
dy
C’est cette variation qui caractérise un mouvement de cisaillement pour lequel les différentes
couches ont des déplacements relatifs les unes par rapport aux autres.
La vitesse de déformation est définie par :
dγ d dξ d dξ dvx
γ̇ = = = = , (6)
dt dt dy dy dt dy
car la composante suivant x de la vitesse est vx = dξ/dt. On voit que la vitesse de déformation est
égale au taux de cisaillement. Dans le cas particulier du cisaillement simple, le taux de cisaillement

10
n’a qu’une seule composante suivant x, dvx /dy = γ̇. Dans le cas général, le taux de cisaillement
(ou vitesse de déformation) associé(e) aux directions i et j est :
∂vi ∂vj
+ = 2eij , (7)
∂xj ∂xi
où le tenseur e est appelé le tenseur taux de déformation.

3.3 Décomposition du champ de gradient de vitesse au voisinage d’un

point
Considérons un élément matériel de fluide situé au point r avec une vitesse v(r,t) et un élément
situé au point r + dr avec une vitesse v + dv, un développement au premier ordre donne :
dv = grad v · dr, (8)
où grad v est le tenseur des gradients de vitesse.
On peut décomposer ce tenseur ∂vi /∂xj en une partie symétrique, le tenseur taux de déformation
eij , et une partie antisymétrique, le tenseur taux de rotation ωij :
∂vi 1 ∂vi ∂vj 1 ∂vi ∂vj
= ( + )+ ( − )
∂xj 2 ∂xj ∂xi 2 ∂xj ∂xi
= eij + ωij . (9)
L’équation (8) peut s’écrire alors :
dv = e · dr + ω · dr, (10)
Lorsque le fluide est incompressible, la trace du tenseur symétrique eij est nulle, soit ekk =
exx + eyy + ezz = 0. L’écoulement correspondant à e · dr est un écoulement de déformation pure
encore appelé écoulement élongationnel. Le tenseur antisymétrique ωij n’a que trois composantes
indépendantes. Il est possible de lui associer un vecteur vitesse de rotation Ω = rotv/2 où rotv est
la vorticité. L’écoulement correspondant à ω · dr = Ω ∧ dr est un écoulement de rotation en bloc.
De manière générale, pour un fluide incompressible, un écoulement peut être considéré comme la
superposition d’un écoulement de déformation pure et d’un écoulement de rotation pure 4 .
Dans le cas particulier du cisaillement simple de la figure 12, la décomposition du tenseur des
gradients de vitesse en une partie symétrique et une partie antisymétrique s’écrit :
     
0 γ̇ 0 0 γ̇/2 0 0 γ̇/2 0
grad v =  0 0 0  =  γ̇/2 0 0  +  −γ̇/2 0 0 . (11)
0 0 0 0 0 0 0 0 0
Le terme symétrique a une trace nulle, il n’y a donc pas de variation en volume des éléments de
fluides. Le terme symétrique a des axes propres à 45◦ et les valeurs propres sont γ̇/2 et −γ̇/2.
Le terme antisymétrique donne une rotation à la vitesse angulaire −γ̇/2. La déformation d’un
élément de fluide carré dans un écoulement de cisaillement simple est présentée sur la figure 14.
On retrouve que le taux de cisaillement ou vitesse de déformation n’a qu’une seule composante
γ̇ = 2exy = 2eyx.

4 Lois constitutives
4.1 Fluides newtoniens et non-newtoniens
Les équations qui gouvernent l’écoulement ou la déformation d’un matériau sont de deux types :
– Les équations de conservation (conservation de la masse et de la quantité de mouvement)
qui sont les mêmes pour tous les matériaux.
4. Lorsque le fluide est compressible, le terme symétrique isotrope ekk 6= 0 va donner un accroissement relatif
de volume d’un élément de fluide, alors que le terme symétrique anisotrope, eij avec i 6= j, correspond à une
déformation pure sans changement de volume.

11
 -

Fig. 14 – Déformation d’un élément de fluide carré dans un écoulement de cisaillement simple :
déformation sans changement de volume avec des axes propres à 45◦ puis rotation.

– Les relations constitutives qui dépendent des propriétés et de la nature du matériau.

Le matériau est traité comme un milieu continu et est caractérisé par une masse volumique
ρ(r,t) et un champ de vitesse v(r,t), qui sont (en général) des fonctions continues de la position
r = (x,y,z) et du temps t. Dans la plupart des écoulements, on peut négliger la compressibilité du
matériau et la condition d’incompressibilité s’écrit :
div v = 0. (12)
On écrit la conservation de la quantité de mouvement sous la forme de Cauchy :
Dv ∂v
ρ = ρ[ + (v · grad)v] = div σ + ρg, (13)
Dt ∂t
où g est l’accélération de la gravité.
Pour l’équation constitutive, on a besoin d’une relation entre la contrainte σ et l’écoulement
v(r,t).
Pour un fluide newtonien, comme l’eau dans des conditions ordinaires de température et de
pression par exemple, la relation entre contrainte et taux de cisaillement est linéaire :
σij = −pδij + 2µeij , (14)
c’est-à-dire que le déviateur (partie non isotrope du tenseur des contraintes) ou encore tenseur des
contraintes de viscosité s’écrit :
dij = 2µeij , (15)
où la constante µ est la viscosité du fluide qui ne dépend pas des contraintes ou du taux de
cisaillement. Pour le cisaillement simple de la figure 12, cette relation s’écrit :
σ = µγ̇. (16)
Pour les fluides newtoniens, l’équation de conservation de la quantité de mouvement (13) devient
l’équation de Navier-Stokes :
Dv ∂v
ρ = ρ[ + (v · grad)v] = −grad pd + µ∆v, (17)
Dt ∂t
où la pression dynamique pd est en fait un champ de pression modifié de façon à intégrer les forces
extérieures conservatives (c’est à dire qui dérivent d’un potentiel comme la gravité).
De nombreux fluides ne sont pas newtoniens. Par exemple, pour les cristaux liquides, il y a
une relation tensorielle entre σij et eij . Pour d’autres fluides isotropes, le rapport des contraintes
aux vitesses de déformation dépend des contraintes, des vitesses de déformation et éventuellement
du temps. Le but de la rhéologie est de déterminer la relation entre contrainte et déformation ou
vitesse de déformation qui détermine toutes les propriétés rhéologiques du matériau. On cherche
une relation σij = f (eij ).
Dans la suite on note σ la contrainte et γ̇ la vitesse de déformation car nous considérons que
les caractérisations du matériau ont été faites en utilisant un cisaillement simple. Nous cherchons
la relation σ = f (γ̇).

12
4.2 Viscoélasticité linéaire
Un comportement non-newtonien important est la viscoélasticité. La réponse du fluide à une
déformation présente à la fois un aspect visqueux (où les contraintes sont proportionnelles aux
vitesses de déformation) et un aspect élastique (où les contraintes sont proportionnelles aux
déformations). Un exemple de fluide viscoélastique est la pâte de silicone connue sous le nom
de “silly-putty”. Une boule de “silly-putty” rebondit sur le sol comme une balle élastique (aux
temps courts) mais s’étale comme un fluide visqueux (aux temps longs) si on la pose sur une
surface horizontale.

4.2.1 Fonctions fluage et relaxation

Les fonctions fluage et relaxation sont les fonctions essentielles en viscoélasticité linéaire. Elles
sont définies de la façon suivante :
– La fonction fluage, f (t), est la déformation subie par le matériau lorsqu’on impose à ce dernier
une contrainte d’amplitude unité au temps t = 0, contrainte qui est maintenue constante
(voir la figure 15).
– La fonction relaxation, g(t), est la contrainte résultant de l’application d’une déformation
d’amplitude unité à l’instant initial t = 0, déformation qui est maintenue constante au cours
du temps (voir la figure 16).
σ f (t)
6 6

- -
t t

Fig. 15 – Fonction fluage.

La connaissance de la fonction fluage ou relaxation d’un matériau permet de déterminer

toutes les propriétés viscoélastiques du matériau. Lorsqu’on dispose d’un rhéomètre imposant
une contrainte constante, σ0 , on mesure une déformation γ(t) = σ0 f (t). Lorsqu’on dispose d’un
rhéomètre imposant une déformation constante, γ0 , on mesure une déformation σ(t) = γ0 g(t).

γ g(t)
6 6

- -
t t

Fig. 16 – Fonction relaxation.

13
4.2.2 Solide élastique parfait
L’équation rhéologique d’un solide élastique parfait (loi de Hooke) s’écrit :
γ(t) = Jσ(t), (18)
où J s’appelle la complaisance élastique. On utilise également le module de rigidité ou module de
cisaillement G = 1/J. La relation (18) s’applique au cisaillement simple de la figure 12. Pour un
solide élastique parfait, la déformation et la contrainte sont reliées par une relation linéaire mais
le coefficient de proportionnalité dépend du type de déformation imposée (c’est-à-dire, dépend du
module d’Young, du coefficient de Poisson ou du module de compression uniforme).

f (t)
6

-
t

Fig. 17 – Fonction fluage d’un solide élastique parfait.

6σ

σ
?

Fig. 18 – Ressort de complaisance élastique J.

De la loi de Hooke, on déduit les propriétés suivantes :

– Dès qu’une contrainte est appliquée, instantanément une déformation prend naissance, pro-
portionnelle à la contrainte. Inversement, si la contrainte est ramenée à zéro, immédiatement
la déformation s’annule. La déformation élastique est instantanée et récupérable.
– Le comportement est bien solide. Soumis à une contrainte constante, le matériau atteint
instantanément un état d’équilibre.
– La fonction fluage est f (t) = J pour t ≥ 0 (voir la figure 17).
Un solide élastique parfait est symbolisé par un ressort de complaisance élastique J (voir la
figure 18).

4.2.3 Liquide visqueux newtonien

L’équation rhéologique d’un liquide visqueux newtonien s’écrit :
dγ(t) 1
= σ(t). (19)
dt µ

14
En tenant compte des conditions initiales, à t < 0, σ(t) = 0 et γ(t) = 0, l’équation (19) s’intègre :

1 t
Z
γ(t) = σ(t′ )dt′ . (20)
µ 0

f (t)
6

1
pente µ

-
t

Fig. 19 – Fonction fluage d’un liquide visqueux newtonien.

6σ

σ
?

Fig. 20 – Amortisseur de viscosité µ.

On en déduit les propriétés suivantes :

– Le liquide visqueux newtonien se souvient de toutes les contraintes qui lui ont été im-
posées dans le passé. En effet, l’expression de γ(t) dépend de toutes les valeurs prises par la
contrainte de 0 à t.
– Si la contrainte est ramenée à 0 à un instant donné, la déformation demeure constante et
égale à la valeur à cet instant. La déformation est irrécupérable.
– Le comportement est bien celui d’un liquide. Soumise à une contrainte constante, σ0 , la
déformation, γ(t) = σ0 t/µ croı̂t linéairement avec le temps. Le matériau s’écoule indéfiniment.
– La fonction fluage est f (t) = t/µ pour t ≥ 0 (voir la figure 19).
Un liquide visqueux newtonien est symbolisé par un amortisseur de viscosité µ (voir la figure
20).
Nous avons présenté les deux comportements viscoélastiques linéaires élémentaires. Le compor-
tement viscoélastique linéaire général se définit en construisant un modèle constitué d’un assem-
blage analogique et symbolique de ressorts (solides élastiques parfaits) et d’amortisseurs (liquides
visqueux newtoniens) en série ou en parallèle.

4.2.4 Solide de Kelvin-Voigt

Le solide de Kelvin-Voigt est constitué par l’association en parallèle d’un ressort de complai-
sance élastique J et d’un amortisseur de viscosité µ (voir la figure 21). La contrainte imposée à

15
l’ensemble est la somme des contraintes de chaque branche :

σ = σe + σv , (21)

où σe = γe /J dans la branche du ressort et σv = µdγv /dt dans la branche de l’amortisseur. La

déformation subie est identique dans chaque branche et égale à la déformation totale subie par
l’ensemble :
γ = γe = γv . (22)
L’équation rhéologique du solide de Kelvin-Voigt est :
dγv 1
σ=µ + γe . (23)
dt J

6σ

µ
J

σ
?

Fig. 21 – Association en parallèle d’un ressort de complaisance élastique J et d’un amortisseur

de viscosité µ.

f (t)
6
J

-
t

Fig. 22 – Fonction fluage d’un solide de Kelvin-Voigt.

La fonction fluage de ce modèle se trouve par la résolution de l’équation différentielle suivante :

df 1
µ + f = 1. (24)
dt J
On obtient :
−t
f (t) = J(1 − exp ), (25)
θ
où θ = Jµ est appelé temps de retard. La fonction fluage du solide de Kelvin-Voigt est une exponen-
tielle croissante possédant une asymptote horizontale d’amplitude J. On a donc un comportement
élastique aux temps longs après un régime transitoire de durée θ (voir la figure 22). Il s’agit d’une
élasticité retardée par la présence de l’élément visqueux.

16
4.2.5 Liquide de Maxwell
Le liquide de Maxwell est constitué de l’association en série d’un ressort de complaisance
élastique J et d’un amortisseur de viscosité µ (voir la figure 23). La contrainte imposée à l’ensemble
est supportée en totalité par chaque élément :

σ = σe = σv , (26)

où σe = γe /J dans l’élément composé du ressort et σv = µdγv /dt dans dans l’élément composé de
l’amortisseur. La déformation totale est la somme des déformations subies par chaque élément :

γ = γe + γv . (27)

L’équation rhéologique du liquide de Maxwell est :

dγ 1 dσ
= σ+J . (28)
dt µ dt
Sa résolution est immédiate en tenant compte des conditions aux limites, à t < 0, σ = 0 et γ = 0 :

1 t
Z
γ(t) = Jσ(t) + σ(t′ )dt′ . (29)
µ 0

6σ

σ
?

Fig. 23 – Association en série d’un ressort de complaisance élastique J et d’un amortisseur de

viscosité µ.

La fonction fluage de ce modèle est :

t
f (t) = J + . (30)
µ
Aux temps longs, le comportement est celui d’un fluide visqueux newtonien car la déformation
augmente linéairement avec le temps (voir la figure 24). La différence avec un fluide newtonien est
la présence d’une élasticité instantanée donnée par la complaisance élastique J.

17
 f (t)
6

1
pente µ
J

-
t

Fig. 24 – Fonction fluage d’un liquide de Maxwell.

4.2.6 Comportement viscoélastique linéaire général

Toutes les différentes associations de ressorts et d’amortisseurs donnent lieu à un comporte-
ment viscoélastique linéaire. A partir des lois d’association en série ou en parallèle, on peut déduire
l’équation rhéologique et les fonctions fluage et relaxation. On peut montrer que toutes les asso-
ciations, aussi complexes soient-elles, peuvent être représentées par des modèles de structure bien
définie : le modèle de Kelvin-Voigt généralisé et le modèle de Maxwell généralisé.
Le modèle de Kelvin-Voigt généralisé est constitué par l’association en série d’un liquide de
Maxwell de viscosité µ0 et complaisance J0 et de N solides de Kelvin-Voigt de viscosité µi et
complaisance Ji . Comme la fonction fluage d’une association en série est la somme des fonctions
fluages des différents éléments, on peut en déduire immédiatement la fonction fluage de ce modèle
général (représentée sur la figure 25) :
N
t X −t
f (t) = J0 + + Ji (1 − exp ), (31)
µ0 i=1 θi

où θi = Ji µi est le temps de retard du ie solide de Kelvin-Voigt. On a une élasticité instantanée

donnée par la complaisance J0 , puis une variation assez rapide correspondant à l’élasticité retardée
des différents solides de Kelvin-Voigt et enfin un comportement visqueux newtonien aux temps
longs caractérisé par une variation linéaire de pente 1/µ0 . Ce modèle permet d’obtenir aussi un
comportement solide aux temps longs en posant µ0 → ∞. La fonction fluage, donnée par :
N
X −t
f (t) = J0 + Ji (1 − exp ), (32)
i=1
θi

présente une asymptote horizontale (voir la figure 26).

Le comportement général peut être également décrit par d’autres modèles équivalents, de struc-
tures différentes mais conduisant à la même équation rhéologique. Une autre structure équivalente
est le modèle de Maxwell généralisé. Il est constitué par l’association en parallèle de N liquides
de Maxwell ayant un module de rigidité Gi et une viscosité µi . Ce modèle décrit un comporte-
ment liquide avec élasticité instantanée et retardée et est équivalent au modèle de Kelvin-Voigt
généralisé précédent.
Selon que l’une des branches se réduit à un amortisseur ou à un ressort libre, on peut décrire
d’autres types de comportement : liquide sans élasticité instantanée, solide avec élasticité instan-
tanée et retardée, solide sans élasticité instantanée, par exemple.
Avec des rhéomètres à fluage, on privilégiera le modèle de Kelvin-Voigt dont la fonction fluage
est simple, tandis qu’avec des rhéomètres à relaxation on utilisera le modèle de Maxwell.

18
 ← régime newtonien
f (t)
6

← élasticité retardée

J0 ← élasticité instantanée
-
t

Fig. 25 – Fonction fluage du modèle de Kelvin-Voigt généralisé.

f (t)
6

J0
-
t

Fig. 26 – Fonction fluage du modèle de Kelvin-Voigt généralisé avec µ0 → ∞.

4.3 Liquides nonlinéaires

Au delà de la viscoélasticité linéaire et plus loin dans le régime nonlinéaire, on considère des
modèles ad hoc qui sont des approximations plus ou moins bonnes du comportement de certains
matériaux. On utilise principalement des modèles pour lesquels la viscosité dépend du taux de
cisaillement :
σ = µ(γ̇)γ̇, (33)
sans effet dépendant du temps.

4.3.1 Liquides plastiques

Ce sont des liquides présentant une contrainte seuil d’écoulement que nous avons évoqués dans
la section 2.1.
La représentation la plus simple d’un fluide à seuil est le modèle de Bingham qui donne la
relation suivante entre contrainte et taux de cisaillement :

σ = σc + µp γ̇, (34)

où σc est la contrainte seuil et µp la viscosité plastique. Le rhéogramme d’un fluide de Bingham est
représenté sur la figure 27. Le matériau ne commence à s’écouler qu’au-delà du seuil et présente
ensuite un comportement newtonien.
Certains matériaux ne présentent pas un comportement newtonien au-delà du seuil. C’est le
cas du fluide de Casson pour lequel la relation utilisée est :
√ √ p
σ = σc + µp γ̇. (35)

19
 σ
6

pente µp
σc

-
γ̇

Fig. 27 – Rhéogramme d’un fluide de Bingham.

On rencontre d’autres relations comme celle de Herchel-Buckley :

σ = σc + k γ̇ n , (36)

où le coefficient k et l’exposant n sont à déterminer.

Des modèles analogiques peuvent aussi rendre compte des déformations d’écoulement plastique.
On est amené à introduire dans ces modèles un élément non linéaire : le patin représenté sur la
figure 28. C’est un limiteur de contrainte. Le patin représente le frottement solide : il glisse au delà
d’un seuil de contrainte σc , la contrainte reste alors constante et égale à la valeur seuil (voir la
figure 29).

6σ

σc

-
σ
?

Fig. 28 – Patin ayant un seuil de contrainte σc .

σ
6

σc

-
γ

Fig. 29 – Variation de la contrainte avec la déformation aux extrémités d’un patin.

20
 La représentation symbolique d’un fluide de Bingham est l’association en parallèle d’un patin
ayant un seuil de contrainte σc et d’un amortisseur de viscosité µp représentée sur la figure 30. La
contrainte imposée à l’ensemble est la somme des contraintes de chaque branche :

σ = σp + σv , (37)

où σp = σc dans la branche du patin lorsqu’il glisse et σv = µp γ˙v dans la branche de l’amortisseur.
La déformation subie est identique dans chaque branche et égale à la déformation totale subie par
l’ensemble :
γ = γp = γv . (38)
L’équation rhéologique correspondante est donnée par l’équation (34) et le rhéogramme est celui
de la figure 27.

6σ

µp
σc

σ
?

Fig. 30 – Association en parallèle d’un patin ayant un seuil de contrainte σc et d’un amortisseur
de viscosité µp .

4.3.2 Fluide rhéofluidifiant

Ce sont des liquides qui ont une viscosité qui décroı̂t avec le taux de cisaillement et que nous
avons présentés dans la section 2.1.
Une loi empirique très utilisée pour la variation de la viscosité avec le taux de cisaillement est
la loi de puissance (proposée pour la première fois par Ostwald en 1925) :

σ = k γ̇ n , (39)

où le coefficient k et l’exposant n < 1 sont à déterminer empiriquement. La viscosité s’écrit donc :

µ = k γ̇ n−1 . (40)

Le cas où n = 1 correspond au comportement newtonien. Il est judicieux de tracer le rhéogramme

en coordonnées log-log comme sur la figure 31 :

log σ = log k + n log γ̇, (41)

où la pente de la droite donne l’exposant n et l’ordonnée à l’origine la constante k.

Une extension du modèle en loi de puissance est la relation de Carreau qui fait intervenir 5
paramètres :
µ − µ∞ n−1
= [1 + (λγ̇)a ] a , (42)
µ0 − µ∞
où µ0 est la viscosité à cisaillement nul, µ∞ est la viscosité à cisaillement infini, λ une constante
de temps et a un paramètre qui décrit la transition entre le comportement à faible cisaillement et
la région en loi de puissance. Ces différents coefficients sont à déterminer empiriquement.
Une autre approche est celle de Shandraw qui admet la relation :

σ = σ0 + µ∞ γ̇ − σ0 exp(−bγ̇). (43)

21
 log σ
6



 

log k  pente n

-
log γ̇

Fig. 31 – Rhéogramme d’un fluide rhéofluidifiant en loi de puissance.

σ
6 pente µ∞ 

 

σ0 
6
Q
?

-
γ̇

Fig. 32 – Rhéogramme du fluide de Shandraw.

6σ

σc

σ
?

Fig. 33 – Association en série d’un patin ayant un seuil de contrainte σc et d’un amortisseur de
viscosité µ.

Le rhéogramme correspondant présenté sur la figure 32 a une asymptote oblique, l’écart Q =

σ0 exp(−bγ̇) avec celle-ci étant exponentiellement décroissant.
On peut utiliser le patin pour représenter le comportement fluidifiant par un modèle sym-
bolique. Le schéma de la figure 33 comportant en série un patin et un amortisseur donne le
rhéogramme de la figure 34. Tant que la contrainte est inférieure au seuil du patin σc , la réponse

22
 σ
6
σc











-
γ̇

Fig. 34 – Rhéogramme correspondant à l’association en série d’un patin et d’un amortisseur.

est celle de l’amortisseur seul. Quand le patin se met à glisser (σ = σc = µγ̇), la contrainte demeure
constante. Cette représentation est sommaire mais on peut la raffiner en associant en parallèle un
ensemble de branches de ce type.

4.3.3 Fluide rhéoépaississant

Ce sont des liquides qui ont une viscosité qui croı̂t avec le taux de cisaillement et que nous
avons mentionnés dans la section 2.1. On utilise la loi de puissance pour les représenter :

σ = k γ̇ n , (44)

avec un exposant n > 1 qui est d’autant plus grand que le matériau s’écarte du comportement
newtonien.

4.4 Comportement dépendant du temps

Certains liquides peuvent avoir un écoulement dont les caractéristiques dépendent des traite-
ments antérieurs. Dans ce cas, la viscosité apparente n’est plus fixée pour une valeur donnée de
la contrainte ou du taux de déformation mais dépend également du temps.
C’est le cas des fluides thixotropes que nous avons mentionnés dans la section 2.2 et pour lesquels
la viscosité apparente a tendance à décroı̂tre dans le temps quand on lui applique une contrainte (ou
une vitesse de cisaillement) constante. Ce comportement complique considérablement les mesures
rhéologiques car on voit se manifester des effets d’hystérésis.

σ
6 A

?

B
 
9
 -
O γ̇

Fig. 35 – Boucle d’hystérésis de référence d’un gel fluidifiant.

23
 6

z
h
6
H- y
x jr
H
θ

Fig. 36 – Rhéomètre de Couette.

On procédera de préférence en réalisant des rhéogrammes dans des conditions expérimentales

fixées et bien reproductibles, comme la boucle de référence d’un gel fluidifiant présentée sur la
figure 35 tirée de la référence [3]. Cette boucle est effectuée de la manière suivante :
– Charge en contrainte (ou taux de cisaillement) pendant un temps tc court par rapport à
toutes les autres durées expérimentales.
– Attente sous cisaillement maximal pendant un temps ta .
– Décharge pendant un temps total td .
Au cours de la charge (partie OA), le gel est rhéofluidifiant. Pendant l’attente (partie AB), la
viscosité apparente diminue à cause du changement de structure (transformation de l’état gel à
l’état fluide). Pendant la décharge (partie BO), la viscosité est constante et le retour est linéaire
car la structure est détruite.
Dans quelques cas plus rares, on trouve aussi le comportement opposé, antithixotrope. À
contrainte (ou vitesse de cisaillement) donnée, la viscosité apparente augmente avec le temps.
Le terme d’antithixotropie est à préférer à celui de rhéopexie, moins bien défini, non totalement
équivalent et qui désigne la propriété de solidification progressive par agitation.

5 Rhéomètrie
5.1 Rhéomètres simples
5.1.1 Rhéomètre plan
La situation idéale est celle où contrainte et vitesse de déformation sont constantes dans tout
le volume du matériau. C’est ce que l’on rencontre dans un matériau homogène soumis à un
écoulement de cisaillement simple comme celui de la figure 12. Les contraintes tangentielles, σ =
σxy , sont constantes et égales aux contraintes sur les plans. La vitesse de déformation est la même
en tout point et vaut γ̇ = U/h. La viscosité de cisaillement est définie comme le rapport de la
contrainte tangentielle et du taux de cisaillement µ = σ/γ̇. Si le matériau est inhomogène, le profil
de vitesse n’est plus celui de l’écoulement de cisaillement simple car la viscosité varie de point en
point. On en est réduit à mesurer une viscosité effective que l’on peut exploiter si l’on dispose d’un
modèle théorique du milieu.
Les autres écoulements viscométriques sont issus de l’écoulement de cisaillement simple et sont
des généralisations pratiques de celui-ci.

24
5.1.2 Rhéomètre de Couette
Dans le rhéomètre de Couette, le fluide est placé dans l’entrefer de deux cylindres coaxiaux de
rayons R1 et R2 avec R1 < R2 et de hauteur h tournant avec des vitesses angulaires respectives
Ω1 et Ω2 (voir la figure 36). En pratique, un des cylindres est fixe et l’autre tourne. C’est souvent
le cylindre intérieur qui tourne, bien que cela puisse conduire à une instabilité de Taylor-Couette
au delà d’une vitesse de rotation critique. On peut soit imposer cette vitesse angulaire et mesurer
le moment M du couple à appliquer pour qu’elle soit maintenue, soit appliquer un couple donné
et mesurer la vitesse angulaire prise par le cylindre mobile.

Détermination de la contrainte Considérons le système de coordonnées cylindriques (r,θ,z)

représentées sur la figure 36. Compte tenu de la symétrie et de l’absence de gradient de pression
axial, vz = 0. Par ailleurs la vitesse vθ est indépendante de θ. La relation d’incompressibilité (12)
donne (voir l’appendice A) :
1 ∂
(rvr ) = 0. (45)
r ∂r
La vitesse vr ne peut donc être que de la forme vr = C/r où C est une constante. Or vr s’annule
sur les parois solides en R1 et R2 ; donc vr est nulle dans tout le fluide. Compte tenu de la symétrie
du champ de vitesse (0,vθ ,0), σrθ est la seule composante non nulle du tenseur des contraintes. La
composante angulaire de l’équation de la quantité de mouvement, équation (13), se réduit à (voir
l’appendice A) :
1 ∂ 2
(r σrθ ) = 0. (46)
r2 ∂r
Le moment est donc conservé :

M = M1 = M2 = M (r) = 2πhr2 σrθ , (47)

où M1 et M2 sont les moments des couples appliqués aux cylindres intérieur et extérieur. La
contrainte vaut :
M
σ = σrθ = . (48)
2πhr2

Détermination du taux de cisaillement Le taux de cisaillement est donné par (voir l’appen-
dice A) :
∂vθ vθ d vθ
γ̇ = 2erθ = − = r ( ), (49)
∂r r dr r
où vθ /r est la vitesse angulaire.

Entrefer étroit Lorsque l’entrefer utilisé est étroit, en pratique (R2 − R1 )/R1 < 0.15, on peut
considérer que la contrainte et le taux de cisaillement sont constants dans tout le volume de
l’entrefer. La valeur constante de la contrainte sera prise égale à :
M
σ ≈ σ(R1 ) = . (50)
2πhR12

Le taux de cisaillement est égal à une constante K :

d vθ
γ̇ = r ( ) = K. (51)
dr r
En intégrant cette équation entre le cylindre intérieur et le cylindre extérieur,
Ω2 R2
vθ dr
Z Z
d( ) = K , (52)
Ω1 r R1 r

25
on trouve l’expression de cette constante :

Ω2 − Ω1 (Ω2 − Ω1 )R1
γ̇ = K = R2
≈ . (53)
ln R1 R2 − R1

À partir de la mesure du moment appliqué et de la vitesse angulaire entre les deux cylindres, on
en déduit la viscosité :
σ M R2 − R1
µ= ≈ . (54)
γ̇ 2πh(Ω2 − Ω1 ) R13

Entrefer large Il faut se donner a priori un modèle rhéologique (matériau newtonien, plastique,
rhéofluidifiant, ou autre), en déduire le taux de cisaillement en tout point, intégrer dans l’entrefer
et vérifier a posteriori que le modèle est le bon. Prenons deux exemples :
– Dans le cas d’un modèle de fluide newtonien, l’équation rhéologique σ = µγ̇ s’écrit :
M d vθ
= µr ( ). (55)
2πhr2 dr r
En intégrant cette équation entre le cylindre intérieur et le cylindre extérieur,
Ω2 R2
vθ M dr
Z Z
d( )= , (56)
Ω1 r 2πhµ R1 r3

on en déduit :
M R22 − R12
Ω2 − Ω1 = . (57)
4πhµ R12 R22
Si la différence de vitesse angulaire des deux cylindres est une fonction linéaire du couple
appliqué, le choix d’un modèle newtonien est correct. Sinon, il faut supposer un autre modèle
pour le fluide.
– Dans le cas d’un modèle de fluide en loi de puissance, l’équation rhéologique σ = k γ̇ n s’écrit :
M d vθ
2
= k[r ( )]n . (58)
2πhr dr r
En intégrant cette équation entre le cylindre intérieur et le cylindre extérieur,
Ω2 R2
vθ M 1 dr
Z Z
d( )=( )n 2 , (59)
Ω1 r 2πhk R1 r1+ n
on en déduit :
n M 1 − n2 −2
Ω2 − Ω1 = ( ) n (R1 − R2 n ). (60)
2 2πhk
Si ln(Ω2 − Ω1 ) est une fonction linéaire de ln M , l’hypothèse d’une loi de puissance est
correcte, sinon il faut choisir un autre modèle.
La plupart des rhéomètres de Couette utilisés en laboratoire ont un fonctionnement automatisé
et ont un entrefer étroit. La précision des mesures est bonne et est améliorée en diminuant les effets
d’extrémités qui perturbent l’écoulement. Ce type de rhéomètre nécessite néanmoins de disposer
d’un volume d’échantillon relativement important.

5.1.3 Rhéomètre cône-plan

Dans le rhéomètre cône-plan, le fluide est placé entre un plan et un cône de rayon R comme
présenté sur la figure 37. Le cône fait un angle ψ avec le plan et peut tourner à une vitesse angulaire
Ω. Ici encore, on peut travailler soit à vitesse imposée soit à moment imposé. À partir de la mesure
du moment M et de la vitesse angulaire Ω, on en déduit la viscosité. En effet, si l’angle du cône est

26
 Ω
6

R
```  -
``` ψ

Fig. 37 – Rhéomètre cône-plan.

assez petit (ψ < 5◦ ), la contrainte et la vitesse de déformation sont constantes dans tout l’espace
occupé par l’échantillon (toujours homogène) et sont respectivement égales à :
3M
σ= , (61)
2πR3
Ω
γ̇ = . (62)
ψ
La viscosité est donnée par :
3M ψ
µ= . (63)
2πΩR3
Ce rhéomètre est très utile lorsqu’on ne dispose que d’une faible quantité de fluide car en
général une goutte suffit, le fluide tenant entre le cône et le plan grâce à la tension capillaire. Il est
aussi très intéressant car on peut atteindre des vitesses de cisaillement élevées. Les inconvénients
sont que la mesure est très sensible à la position du cône (qui est tronqué et remplacé par un
petit méplat) et que ce type d’appareil ne convient pas pour des matériaux complexes et fragiles
comme les suspensions de particules. Un grand avantage est que l’on peut mesurer la poussée sur
la plaque ce qui permet une mesure des contraintes normales.

5.1.4 Autres types de rhéomètres

D’autres types de rhéomètres comme le rhéomètre de type Poiseuille et le viscosimètre à chute
de bille sont également utilisés. Mais ils ne donnent des résultats satisfaisants que pour des fluides
newtoniens.

5.2 Cisaillement oscillant

On peut étudier de façon dynamique un matériau en lui imposant une contrainte (déformation)
qui varie sinusoı̈dalement au cours du temps. La linéarité des équations entraı̂ne que la déformation
(contrainte) est également sinusoı̈dale et de même fréquence. On pose :

σ(t) = σ0 cos(ωt + δ),

γ(t) = γ0 cos(ωt), (64)

où σ0 et γ0 sont les amplitudes de la contrainte et de la déformation, ω la pulsation et δ le

déphasage entre les deux. Pour faciliter les calculs, on prend une notation complexe :

σ̄(t) = σ0 exp[i(ωt + δ)],

γ̄(t) = γ0 exp(iωt). (65)

La contrainte et la déformation sont reliées par un module de rigidité complexe :

σ̄ = Ḡγ̄ = (G′ + iG′′ )γ̄, (66)

27
où les parties réelle
√ et imaginaire de Ḡ sont le module de conservation et le module de perte avec
tan δ = G′′ /G′ et G′2 + G′′2 = σ0 /γ0 . On peut aussi définir une viscosité complexe :
σ̄ = µ̄γ̇¯ = (µ′ − iµ′′ )γ̇,
¯ (67)
′ ′′ ′′ ′ ′
avec µ = G /ω et µ = G /ω. Le module de conservation G est associé à la réponse élastique,
en phase avec la déformation. Le module de perte G′′ est associé à la réponse de type visqueux,
en quadrature de phase avec la déformation. Ce type de mesure dynamique est particulièrement
adapté à l’étude rhéologique de la viscoélasticité qui a été présentée dans la section 4.2.
Pendant la déformation, la contrainte extérieure appliquée met en jeu une certaine puissance
mécanique par unité de volume :
P = σ γ̇, (68)
qui s’écrit :
P (t) = −σ0 γ0 ω cos(ωt + δ) sin(ωt),
= −σ0 γ0 ω cos(ωt) sin(ωt) cos δ + σ0 γ0 ω sin2 (ωt) cos δ,
γ 2ω
= − 0 G′ (ω) sin(2ωt) + γ02 ωG′′ (ω) sin2 (ωt). (69)
2
La puissance moyenne fournie au cours d’un quart de période peut se décomposer en deux parties :
π
Z 2ω
< P (t) > = P (t)dt,
0
= Pe + Pd . (70)
Le premier terme Pe = −G ′
(ω)γ02 ω/π
change de signe tous les quarts de cycle. C’est l’énergie
de déformation élastique stockée puis restituée. Le second terme Pd = G′′ (ω)γ02 ω/2 est toujours
positif. C’est l’énergie dissipée par les frottements visqueux.

Modules de conservation et de perte

1

0.9
Module de conservation
0.8

0.7

0.6

0.5

0.4

0.3 Module de perte

0.2

0.1

0
-2 -1 0 1 2
10 10 10 10 10
Fréquence réduite

Fig. 38 – Modules de conservation, G′ , et de perte, G′′ , (normalisés par G) pour un liquide de

Maxwell en fonction de la fréquence réduite, ωτ .

Le module de rigidité complexe se calcule facilement pour :

– Un solide élastique parfait, Ḡ = σ̄/γ̄ = 1/J = G avec G′ = G et G′′ = 0.
– Un liquide visqueux newtonien, Ḡ = σ̄/γ̄ = µγ̇/γ̄ ¯ = iωµ avec G′ = 0 et G′′ = ωµ.
– Un liquide de Maxwell, Ḡ = σ̄/γ̄ = iωµ/(1 + iωτ ) où τ = µ/G avec une dépendance en
fréquence analogue à celle d’un circuit électrique résistance-capacité, G′ = Gω 2 τ 2 /(1 + ω 2 τ 2 )
et G′′ = Gωτ /(1 + ω 2 τ 2 ), présentée sur la figure 38.

28
 '$
 '$
-

&% &%

'$ '$

&%
 &%
-

Fig. 39 – Appareil à quatre rouleaux pour produire un écoulement purement élongationnel.

5.3 Contraintes normales

Dans un fluide newtonien soumis à un écoulement de cisaillement simple comme celui de la
figure 12 pour lequel v = (γ̇y,0,0), seules les contraintes tangentielles, σxy = σyx , sont modifiées
par l’écoulement et les contraintes normales restent isotropes et égales à −p. Pour des fluides non-
newtoniens, l’écoulement de cisaillement peut induire des contraintes normales qui ne contribuent
pas à la dissipation et sont donc appelées contraintes élastiques. En éliminant le terme de pression
isotrope, on peut définir les différences entre contraintes normales :

N1 = σxx − σyy , (71)

et
N2 = σyy − σzz . (72)
Elles sont indépendantes de la direction de l’écoulement et, à faible taux de cisaillement, sont
fonctions quadratiques de γ̇, N = Ψγ̇ 2 . À plus fort taux de cisaillement, la dépendance est en
loi de puissance. En général, la première différence N1 est positive et peut être du même ordre,
et même plus grande, que la contrainte de cisaillement. La seconde différence N2 est d’habitude
négative et de 10% plus petite que N1 .
L’anisotropie des contraintes normales a deux effets spectaculaires que nous avons mentionnés
dans la section 2.5 :
– L’ascension le long de l’axe d’un batteur rotatif (effet Weissenberg). Dans l’écoulement en-
gendré par le cylindre tournant, la vitesse est essentiellement azimutale avec un gradient
radial. L’anisotropie des contraintes normales conduit à σθθ > σrr , dans un système de co-
ordonnées cylindriques (r,θ,z), z étant l’axe du cylindre. Il y a une tension le long des lignes
de courant qui tend à pousser le fluide vers l’axe et à le faire remonter.
– Le gonflement à l’extrusion d’un tube. Cet effet est du à la relaxation des contraintes
normales σxx (x étant la direction de l’axe du tube) accumulées pendant l’écoulement à
l’intérieur du tube.

5.4 Viscosité élongationnelle

Un écoulement élongationnel peut être réalisé par un appareil à quatre rouleaux rotatifs ou par
un système de jets opposés (voir les figures 39 et 40). Nous considérons l’écoulement élongationnel

29
 - 

Fig. 40 – Système de jets opposés pour produire un écoulement purement élongationnel.

axisymétrique pour lequel v = ǫ̇(x, − y/2, − z/2) où ǫ̇ est le taux de déformation élongationnel.
Dans cette configuration d’écoulement, on a :
σxy = σxz = σyz = 0, (73)
et, en éliminant le terme de pression isotrope, on définit la viscosité élongationnelle, µE :
σxx − σyy = σxx − σzz = µE ǫ̇. (74)
Pour un liquide newtonien le rapport de Trouton, Tr = µE /µ = 3, ne dépend pas du taux de
cisaillement ou du taux d’élongation. Ce n’est pas le cas pour des fluides non-newtoniens pour
lesquels µE (ǫ̇) et pour lesquels le rapport de Trouton peut être très grand Tr ≈ 103 − 104 .

A Coordonnées cartésiennes, cylindriques et sphériques

En coordonnées cartésiennes (x, y, z), l’équation de continuité s’écrit :
∂vi ∂vx ∂vy ∂vz
= + + = 0, (75)
∂xi ∂x ∂y ∂z
et les équations de conservation de la quantité de mouvement sont :
∂ ∂vj ∂σij
ρ( + )vi = ρFi + , (76)
∂t ∂xj ∂xj
où Fi est une force en volume comme la gravité. Sur chaque coordonnée, elles s’écrivent :
∂ ∂vx ∂vy ∂vz ∂σxx ∂σxy ∂σxz
ρ( + + + )vx = ρFx + + + , (77)
∂t ∂x ∂y ∂z ∂x ∂y ∂z
∂ ∂vx ∂vy ∂vz ∂σxy ∂σyy ∂σyz
ρ( + + + )vy = ρFy + + + ,
∂t ∂x ∂y ∂z ∂x ∂y ∂z
∂ ∂vx ∂vy ∂vz ∂σxz ∂σyz ∂σzz
ρ( + + + )vz = ρFz + + + .
∂t ∂x ∂y ∂z ∂x ∂y ∂z
Le tenseur taux de déformation est :
1 ∂vi ∂vj
eij = ( + ). (78)
2 ∂xj ∂xi

30
 En coordonnées cylindriques (r, θ, z), l’équation de continuité s’écrit :

1 ∂ 1 ∂vθ ∂vz
(rvr ) + + = 0, (79)
r ∂r r ∂θ ∂z
et les équations de conservation de la quantité de mouvement :

∂ ∂ vθ ∂ ∂ v2 ∂σrr 1 ∂σrθ ∂σrz σrr − σθθ

ρ[( + vr + + vz )vr − θ ] = ρFr + + + + , (80)
∂t ∂r r ∂θ ∂z r ∂r r ∂θ ∂z r
∂ ∂ vθ ∂ ∂ vr vθ ∂σθr 1 ∂σθθ ∂σθz σrθ
ρ[( + vr + + vz )vθ + ] = ρFθ + + + +2 ,
∂t ∂r r ∂θ ∂z r ∂r r ∂θ ∂z r
∂ ∂ vθ ∂ ∂ ∂σzr 1 ∂σzθ ∂σzz σzr
ρ( + vr + + vz )vz = ρFz + + + + .
∂t ∂r r ∂θ ∂z ∂r r ∂θ ∂z r
Le tenseur taux de déformation est :
∂vr 1 ∂vθ vr ∂vz
err = ,eθθ = + ,ezz = , (81)
r r ∂θ r ∂z
1 1 ∂vz ∂vθ 1 ∂vr ∂vz 1 ∂ vθ 1 ∂vr
eθz = ( + ),ezr = ( + ),erθ = [r ( ) + ].
2 r ∂θ ∂z 2 ∂z ∂r 2 ∂r r r ∂θ
En coordonnées sphériques (r, θ, ϕ), l’équation de continuité s’écrit :

1 ∂ 2 1 ∂(sin θvθ ) 1 ∂vϕ

2
(r vr ) + + = 0, (82)
r ∂r r sin θ ∂θ r sin θ ∂ϕ
et les équations de conservation de la quantité de mouvement :

∂ ∂ vθ ∂ vϕ ∂ v 2 + vϕ2 ∂σrr 1 ∂σrθ

ρ[(
+ vr + + )vr − θ ] = ρFr + + (83)
∂t ∂r r ∂θ r sin θ ∂ϕ r ∂r r ∂θ
1 ∂σrϕ 1
+ + (2σrr − σθθ − σϕϕ + σrθ cot θ),
r sin θ ∂ϕ r
∂ ∂ vθ ∂ vϕ ∂ vr vθ vϕ2 ∂σθr 1 ∂σθθ
ρ[( + vr + + )vθ + − cot θ ] = ρFθ + +
∂t ∂r r ∂θ r sin θ ∂ϕ r r ∂r r ∂θ
1 ∂σθϕ 1
+ + [(σθθ − σϕϕ ) cot θ + 3σrθ ],
r sin θ ∂ϕ r
∂ ∂ vθ ∂ vϕ ∂ vθ vϕ vr vϕ ∂σϕr 1 ∂σϕθ
ρ[( + vr + + )vϕ + cot θ + ] = ρFϕ + +
∂t ∂r r ∂θ r sin θ ∂ϕ r r ∂r r ∂θ
1 ∂σϕϕ 1
+ + (3σrϕ + 2σθϕ cot θ).
r sin θ ∂ϕ r
Le tenseur des taux de déformation est :
∂vr 1 ∂vθ vr 1 ∂vϕ vθ vr
err = ,eθθ = + ,eϕϕ = + cot θ + , (84)
r r ∂θ r r sin θ ∂ϕ r r
1 sin θ ∂ vϕ 1 ∂vθ
eθϕ = [ ( )+ ],
2 r ∂θ sin θ r sin θ ∂ϕ
1 1 ∂vr ∂ vϕ 1 ∂ vθ 1 ∂vr
eϕr = [ + r ( )],erθ = [r ( ) + ].
2 r sin θ ∂ϕ ∂r r 2 ∂r r r ∂θ

B Dimensions et unités
Les dimensions et unités du cisaillement, de la contrainte et des viscosités dynamique et
cinématique sont données dans le système d’unités CGS (centimètre, gramme et seconde) et dans

31
Fig. 41 – Viscosité d’une solution à 1% de polyéthylène oxyde dans un mélange 50%/50%
eau/glycérine pour diverses fractions massiques de billes de verre, à une température de 25◦ .

le système international d’unités SI. Notons que 1 Poise = 0,1 Pa.s et que 1 Stokes = 10−4 m2 .s −1 .
CGS SI
γ̇ s−1 s−1
σ g.cm−1 .s−2 Pascal = kg.m−1 .s−2
µ Poise = g.cm−1 .s−1 Pa.s = kg.m−1 .s−1
ν = µ/ρ Stokes = cm2 .s −1 m2 .s −1

C Exercices et problèmes
C.1 Course entre un fluide newtonien et un fluide non-newtonien
On considère deux tubes verticaux identiques ouverts sur le haut et comportant deux robinets
identiques en bas. Ces deux tubes sont remplis à la même hauteur d’un liquide newtonien pour
le tube de gauche et d’un liquide non-newtonien pour le tube de droite. Les robinets sont ouverts
simultanément. Au début le fluide non-newtonien s’écoule plus vite mais quand la hauteur de
liquide devient plus faible, il ralentit par rapport au fluide newtonien et à la fin c’est le fluide
newtonien qui gagne en s’écoulant plus vite. Quel est le type de fluide non-newtonien? Expliquer
l’expérience. Donner un exemple de fluide ayant un comportement rhéologique semblable.

C.2 Loi rhéologique à trouver

On considère les variations de viscosité avec le gradient de vitesse représentées sur la figure
41. Dans la partie décroissante des courbes, quelle est la loi rhéologique du fluide? Donner en la
caractéristique à partir des courbes de la figure 41.

C.3 Écoulement dans un tube cylindrique d’un fluide en loi de puissance

On considère un écoulement de Poiseuille dans un tube cylindrique de rayon R induit par une
différence de pression ∆p sur une longueur L (voir la figure 42).
1. En considérant que l’écoulement est axisymétrique et que σzz = −p, montrer que la compo-
sante axiale de l’équation de conservation de la quantité de mouvement s’écrit :
dp 1 ∂(rσzr )
0=− + . (85)
dz r ∂r

32
  -
L

6
R
6
r
-z

Fig. 42 – Coupe du tube cylindrique.

2. Montrer que la variation de la contrainte de cisaillement σzr avec la position radiale r dans
le tube est donnée par :
∆p
σzr = r. (86)
2L
3. Montrer que pour un fluide en loi de puissance (σzr = κγ̇ n avec un taux de cisaillement
γ̇ = dvz /dr) la variation de la vitesse du fluide est donnée par :
 n1 1 1
R1+ n − r1+ n

∆p
vz = − . (87)
2Lκ 1 + n1

4. Quel est le profil de vitesse pour un fluide newtonien (n = 1)? et pour un fluide rhéofluidifiant
(n < 1)?
5. Montrer que le débit volumique de fluide Q est donné par :
 n1
πR3

∆pR
Q=− . (88)
2Lκ 3 + n1

Montrer que l’on retrouve bien la loi de Poiseuille pour un fluide newtonien.

C.4 Écoulement de boues de forage

Les boues de forage utilisées dans l’exploitation pétrolière sont des suspensions concentrées
de particules solides. Leur comportement rhéologique est non newtonien et leur rhéogramme (ca-
ractéristique contrainte - gradient de vitesse) met en évidence l’existence d’une contrainte seuil, σc ,
en dessous de laquelle le fluide ne s’écoule pas. Au delà de la contrainte seuil, on peut représenter
le comportement du fluide soit par une relation linéaire (fluide de Bingham), soit par des relations
plus réalistes (comme celle du fluide de Casson).
Le problème de l’écoulement des boues de forage dans un milieu poreux (comme c’est le cas
dans un puit de forage) peut se ramener schématiquement à l’écoulement dans une série de tubes de
section circulaire. Pour simplifier encore, nous considérons l’écoulement stationnaire de Poiseuille
dans un tube cylindrique de rayon R induit par une différence de pression ∆p sur une longueur L
(voir la figure 42) d’un fluide de Bingham. Ce problème est tiré de la référence [5]

1. En considérant que l’écoulement est axisymétrique et que σzz = −p, montrer que la compo-
sante axiale de l’équation de conservation de la quantité de mouvement s’écrit :

dp 1 ∂(rσzr )
0=− + . (89)
dz r ∂r

33
Fig. 43 – Rhéogramme d’une boue de forage. La courbe en trait plein est l’ajustement d’un modèle
de Casson sur les points expérimentaux avec σc = 7,19 Pa et µp = 0,048 Pa.s. La courbe en trait
pointillé est un modèle de Bingham avec σc = 8,67 Pa et µp = 0,245 Pa.s.

2. En déduire que la contrainte de cisaillement varie linéairement avec la position radiale r :

σzr (r) = rσzr (R)/R, (90)

où σzr (R) = ∆pR/2L.

3. Montrer que le débit est donné par :
Z vz (R)
Q = −π r2 duz . (91)
vz (0)

4. Avec un taux de cisaillement γ̇ = dvz /dr, en déduire que :

σzr (R)
πR3
Z
2
Q=− 3 σzr γ̇(σzr )dσzr . (92)
σzr (R) 0

5. Montrer que pour un fluide de Bingham qui présente une contrainte seuil, σc , en dessous de
laquelle le fluide ne s’écoule pas, et une viscosité plastique, µp , le débit s’écrit :
σzr (R)
πR3
Z
2
Q=− 3 (R)
σzr (σzr − σc )dσzr . (93)
µp σzr σc

6. En rapportant le débit à son expression pour le fluide newtonien Qp (loi de Poiseuille),

montrer qu’on obtient :
4σc σ4
Q/Qp = 1 − + 4c . (94)
3σzr (R) 3σzr (R)
Représenter Q/Qp en fonction de σzr (R)/σc .
7. Montrer que pour un gradient de pression donné, le fluide ne s’écoule que dans des tubes de
rayon supérieur à un rayon critique que l’on déterminera.
8. On fait couler de la boue de forage, dont les caractéristiques sont données dans le rhéogramme
de la figure 43, dans un milieu poreux de 10 cm, avec une différence de pression entre entrée
et sortie de 7 bars (1 bar = 105 Pascal). Quel est le diamètre minimal des pores dans lesquels
s’écoule la boue? Quelles informations sur la distribution de taille des pores peut-on obtenir
à partir de la relation expérimentale débit-différence de pression?

34
C.5 Effet Weissenberg
Certaines solutions de polymère de grande masse moléculaire présentent un effet connu sous le
nom d’effet Weissenberg et qui est dû à l’anisotropie de contraintes normales sous cisaillement. Cet
effet est illustré sur la figure 44. Un barreau cylindrique de rayon R est plongé dans le liquide non-
newtonien de masse volumique ρ. Lorsqu’on fait tourner ce barreau à vitesse angulaire constante
ω, la surface libre du liquide s’élève le long du barreau. Une expérience similaire avec un liquide
newtonien montre une dépression de la surface libre autour du barreau tournant.

Fig. 44 – Ascension d’une solution de polyisobutylène le long d’un barreau tournant. Le rayon
du barreau est 0,32 cm. La vitesse de rotation est : (a) 0, (b) 2 tour/min, (d) 3 tour/min et
(g) 7 tour/min. En l’absence de rotation, l’ascension du liquide est due à des effets capillaires.
Photographies tirées de G. S. Beavers and D. D. Joseph, J. Fluid Mechanics, 69, 475 (1975)

Nous allons établir une théorie simplifiée de cet effet tirée de la référence [5]. Nous supposons
que la déformation de l’interface reste de petite amplitude et nous négligeons les effets de la tension
de surface du liquide. Nous supposons que l’interface est initialement plane (l’angle de contact entre
la surface du liquide et le barreau tournant est 90◦ ). Nous considérons des coordonnées cylindriques
où la direction suivant z est contraire à la gravité.
1. En considérant que le mouvement est purement azimutal (vz = vr = 0) et que l’écoulement
est axisymétrique, montrer que les équations de conservation de la quantité de mouvement
se réduisent à :
– Composante z
∂σzz
= ρg, (95)
∂z
– Composante r
v2 ∂σrr σrr − σθθ
−ρ θ = + , (96)
r ∂r r
qui peut s’écrire, en faisant apparaı̂tre les différences de contraintes normales σθθ − σrr
et σrr − σzz ,
v2 ∂(σrr − σzz ) ∂σzz σrr − σθθ
−ρ θ = + + , (97)
r ∂r ∂r r
– Composante θ
1 ∂r2 σrθ
= 0. (98)
r2 ∂r
2. En considérant que la contrainte de cisaillement est reliée au taux de cisaillement par l’ex-
pression newtonienne :
σrθ = µγ̇ = 2µerθ , (99)

35
 où µ est la viscosité et où le tenseur taux de déformation est donné en coordonnées cylin-
driques par :
1 ∂vθ vθ 1 ∂vr
erθ = ( − + ), (100)
2 ∂r r r ∂θ
et en cherchant des solutions de la forme vθ = An rn de l’équation (98), montrer que le champ
de vitesse est donné par :
ωR2
vθ = , (101)
r
et que le taux de cisaillement est :

ωR2
γ̇ = 2erθ = −2 . (102)
r2
3. En intégrant l’équation (95) à partir de la surface libre z = h(r) où σzz = 0, montrer
que la contrainte normale suivant l’axe de rotation est donnée par le champ de pression
hydrostatique :
σzz = −ρg[h(r) − z]. (103)
4. En intégrant l’équation (97), montrer que la forme de la surface libre est donnée par :

ω2 R4 ρR4
h(r) = (Ψ1 4 − 2 ), (104)
ρg r 2r
pour un fluide présentant une anisotropie de contraintes normales de la forme :

σθθ − σrr = Ψ1 γ̇ 2 , (105)

la seconde différence des contraintes normales, σrr − σzz , étant beaucoup plus petite et donc
négligée. Dans quelles conditions y a-t-il ascension le long du barreau tournant? Calculer la
hauteur d’ascension. Que se passerait-il si le fluide était newtonien?
5. Le fluide utilisé est 15000 fois plus visqueux que l’eau (pour un taux de cisaillement compris
entre 1 et 10 s−1 ). Sa masse volumique est 0,89 g cm−3 . La hauteur d’ascension en fonction
du carré de la vitesse de rotation, pour des barreaux de différents diamètres est présentée
sur la figure 45. Comparer le modèle aux résultats expérimentaux et estimer Ψ1 . Quel rayon
de barreau faut-il avoir pour observer une ascension du liquide ? Peut-on négliger les effets
d’inertie dans le modèle?

Fig. 45 – Hauteur d’ascension d’une solution de polyisobutylène en fonction du carré de la vitesse

de rotation ω, pour trois diamètres de barreau tournant. Données tirées de G. S. Beavers and
D. D. Joseph, J. Fluid Mechanics, 69, 475 (1975)

36
Références
[1] R. Byron Bird, R. C. Armstrong and Ole Hassager, Dynamics of Polymeric Liquids, Volume
1, Fluid Mechanics, 2nd Edition, Wiley, New York (1987).
[2] D. V. Boger and K. Walters, Rheological phenomena in focus, Elsevier, Amsterdam (1993).
[3] G. Couarraze et J. L. Grossiord, Initiation à la rhéologie, Technique et Documentation, La-
voisier, Paris (1983).
[4] P. Coussot et C. Ancey, Rhéophysique des pâtes et suspensions, EDP Sciences (1999).
[5] M. Fermigier, Hydrodynamique Physique - Problèmes résolus avec rappels de cours, Dunod,
Paris (1999).
[6] E. Guyon, J.-P. Hulin et L. Petit, Hydrodynamique Physique, Inter Editions / Éditions du
CNRS (1991), nouvelle édition revue et augmentée, EDP Sciences / CNRS Éditions (2001).
[7] H. Markovitz, Rheological behavior of fluids, film réalisé sous la direction du National Commit-
tee for Fluid Mechanics Films, et distribué par Encyclopaedia Britannica Educational Corp.
(425 N. Michigan Ave., Chicago, Il 60611, USA).

37