Filière Génie Civil

GC2 2020-2021

Mécanique des sols

Chapitre 3:
L’écoulement de l’eau dans les sols

Khamlichi Abdellatif

1
 Plan
• Introduction
• La charge hydraulique
• Le gradient hydraulique
• La loi de Darcy
• Equation de l’hydraulique souterraine
• Cas des écoulements à deux dimensions en milieu homogène et isotrope
• Conditions aux limites
• Méthode de Kozeny/Casagrande pour la localisation de la ligne de suintement
d’un barrage en terre
• Formation des sables boulants et des renards
• Force d’écoulement volumique
• La détermination du gradient hydraulique critique
• Le renard hydraulique
• Profondeur critique d’une excavation
• Les filtres granulaires
• Les géotextiles

2
1 Introduction

Nous étudierons les charges hydrauliques responsables de l’écoulement de l’eau

à travers le sol.

Nous étudierons l’influence des forces d’infiltration sur la stabilité interne des sols.

Nous verrons les méthodes de calcul des écoulements dans un sol.

Nous traiterons aussi des filtres protecteurs utilisés pour faire obstacle à l’érosion
interne des sols.

3
2 La charge hydraulique

D’après l’équation de Bernoulli, la charge hydraulique ou charge est donnée par:

2 v2 charge de vitesse
v u →
h= + +z 2g
2g γ w
u
→ charge de pression
γw
v: vitesse de l’eau charge d’élévation
z →
g: accélération gravitationnelle
u: pression de l’eau
γw: poids volumique de l’eau
z: hauteur du point considéré par rapport à une surface de référence (axe des z
dirigé vers le haut)

h exprime l’énergie par unité de poids d’eau.

La charge de vitesse étant négligeable, l’écoulement de l’eau dans le sol dépend
principalement de la charge de pression en un point donné:

u
h
+z
γw
4
2 Le gradient hydraulique

On définit le gradient hydraulique i en un point d’une ligne d’écoulement par

le rapport:
dh
i=−
d

où h exprime la charge hydraulique suivant la ligne d’écoulement (Bernoulli).

dh = −i d

i représente la perte de charge par unité de longueur suivant la ligne de courant.

Cette perte de charge traduit le frottement exercé par l’eau sur le squelette solide,
c'est-à-dire la viscosité.

La poussée d’écoulement qui en résulte est à l’origine de nombreux sinistres:

glissement de terrain, formation de renard hydraulique, …
5
2 Le gradient hydraulique

La relation définissant le gradient hydraulique dans un milieu unidirectionnelle

se généralise aisément dans un milieu à deux ou trois dimensions. On a alors

 
i = −grad(h)

 
Si i = 0 la charge hydraulique est la même en tout point du milieu; l’eau interstitielle
est dite en équilibre hydrostatique.

6
3 La loi de Darcy

La loi de Darcy est la loi fondamentale de l’hydraulique des sols. Publiée en 1856
elle exprime la proportionnalité entre la vitesse d’écoulement et le gradient hydraulique.
C’est une loi expérimentale qui s’écrit dans le cas d’un milieu isotrope sous la forme:

 
v = ki


En tout point du sol, le vecteur i et la ligne de courant sont tangents et sont
orientés dans le même sens.

7
3 La loi de Darcy

La loi de Darcy se vérifie en général très bien à condition de rester en régime

laminaire, c’est-à-dire tant que les vitesses restent faibles.
Le coefficient de proportionnalité k est appelé coefficient de perméabilité du sol.
Il s’exprime en cm/s.
Dans le cas d’un milieu anisotrope, la loi de Darcy se généralise en:

   
[ v] =  k   i 

 k xx k xy k xz 
 
k  = k k yy k yz 
   xy 
 k xz k yz k zz 

8
3 La loi de Darcy

On peut aussi écrire la loi de Darcy en terme de la perméabilité intrinsèque sous la

forme suivante:

 k 0  
v = − grad(u) − ρw g 
η

Dans le cas d’un milieu anisotrope, un tenseur de perméabilité intrinsèque

est introduit et la loi de Darcy prend la forme:

 1  
v = − k 0 grad(u) − ρ w g 
η

9
4 Equation de l’hydraulique souterraine

Considérons un sol isotrope soumis à un écoulement. En combinant la condition

de continuité et la loi de Darcy, on obtient le système suivant:
 ∂v x ∂v y ∂v z
 + + =0
 ∂ x ∂y ∂z
 ∂h
 v x = −k
∂x

 v = −k ∂h
 y ∂y

 v = −k ∂h
 z ∂z

D’où l’équation de Laplace qui gère l’écoulement en charge:

∂2h ∂ 2h ∂ 2h
2
+ 2
+ 2
=0 ∆h = 0 (1)
∂x ∂y ∂z

Cette équation concerne la charge hydraulique exclusivement; elle suffit à

caractériser tout écoulement souterrain dans un sol. 10
4 Equation de l’hydraulique souterraine

Dans le cas général où la perméabilité est anisotrope, l’équation de l’hydraulique

souterraine s’écrit:


( )
div k grad(h) = 0 (2)

ou en terme de la pression:

   u 
div  k grad  + z  = 0
  γw 

11
5 Cas des écoulements à deux dimensions en milieu homogène et
isotrope

La plupart des problèmes d’hydraulique des sols peuvent être ramenés à deux
dimensions.

Dans ce contexte l’équation de Laplace s’écrit:

∂ 2h ∂ 2h
2
+ 2
=0
∂x ∂z
12
5 Cas des écoulements à deux dimensions en milieu homogène et
isotrope

On appelle ligne de courant le trajet d’une particule de l’eau.

On appelle tube de courant l’espace compris entre deux lignes de courant.

On appelle lignes équipotentielles les courbes qui passent par les points de même
valeur du potentiel h. Elles ont pour équation générale: h=cste.

Les lignes de courant et les équipotentielles forment un réseau de courbes

orthogonales en tout point.

Equipotentielles
Lignes de courant

13
5 Cas des écoulements à deux dimensions en milieu homogène et
isotrope

Plusieurs méthodes permettent de résoudre l’équation de Laplace (1) ou (2):

- analogique (analogie électrique),
- graphique,
- analytique,
- numérique,...
Il faut savoir aussi comment bien écrire les conditions aux limites.

14
6 Conditions aux limites

Au contact de l’eau avec une interface perméable (comme c’est le cas au

contact de l’eau de la retenue d’un barrage avec le sol), les vitesses de filtration
sont en général très faibles. La pression dynamique (charge de vitesse) est
négligeable et on a affaire à un cas de l'hydrostatique. La condition aux limites
sur l’interface s’écrit alors:

u + ρgz = cste
soit

h = cte

Au contact de l’eau avec une surface imperméable, la vitesse de filtration

normale est quasi-nulle et la condition aux limites s’écrit:

∂h
=0
∂n

15
6 Conditions aux limites

Ligne de suintement (frontière inconnue)

u = patm = 0

h=z Frontière libre

∆h = 0

h = h1 h = h2
∂h
=0
∂n
16
7 Méthode de Kozeny/Casagrande pour la localisation de la ligne de
suintement d’un barrage en terre

f = 0.3 e
d = B − e + f = B − 0.7e
y
m = H2 + d2
x0 = m − d
− y 2 + x 02
x=
x 2x 0

 −1 − 1 + tan 2 (α)
x N = x0
 tan 2 (α)

 y N = − x N tan(α)
17
7 Méthode de Kozeny/Casagrande pour la localisation de la ligne de
suintement d’un barrage en terre

xN

α > 30° : a = (1 − n) x 1 + tan 2 (α)

∆a  N
n=  d d2 H2
a + ∆a α < 30° : a = − − 2
α 2
 cos( ) cos (α) sin (α)
18
7 Méthode de Kozeny/Casagrande pour la localisation de la ligne de
suintement d’un barrage en terre

Nf
q = kH (m3 / s / m 2 ) q = 38m3 / jour / m 2
Ne
Q = 172 m3 / jour / m 19
8 Formation des sables boulants et des renards

En s’écoulant dans le sol, l’eau exerce une pression sur les particules dans la
direction de l’écoulement: c’est la force d’infiltration. Sa valeur est directement
proportionnelle à la perte de charge, c’est-à-dire au gradient hydraulique.

La force d’infiltration modifie la contrainte qui s’exerce entre les particules

de sol (la contrainte effective). Quand l’écoulement est ascendant, la
force d’infiltration diminue la contrainte effective; au contraire, quand il est
descendant, la force d’infiltration augmente la contrainte effective.

20
8 Formation des sables boulants et des renards

Lorsque l’écoulement de l’eau dans le sol est ascendant et que le gradient

hydraulique est suffisamment élevé, la force d’infiltration peut devenir assez
grande pour annuler la contrainte effective. On dit alors que le sol est dans un
état de boulance: les particules de sol flottent pour ainsi dire dans l’eau et ne
peuvent supporter aucune charge.

L’état de boulance survient surtout dans les sables, les sables silteux, les sols
pulvérulents dépourvus de cohésion et dont les particules sont assez petites
pour être emportées par l’écoulement de l’eau. Par contre, l’état de boulance
peut difficilement survenir dans les sables graveleux (particules grosses) et argileux
(faible perméabilité et cohésion liant leurs particules).

21
9 Force d’écoulement volumique

L’écoulement de l’eau entraine une poussée: −Sdu = γ w iSd

La force d’écoulement par unité de volume (force volumique) est:

 
j = γw i
  
j = −γ w grad(h) −ρw gVs
z

jVs


ρsat gVs

Ecoulement

22
10 La détermination du gradient hydraulique critique

Pour définir le gradient hydraulique critique, c’est-à-dire celui qui va provoquer

l’état de boulance, on exprime l’équilibre de l’échantillon de sol montré sur la figure.

Fb = σb A = γ w ( h + h w + L ) A
Fh = σh A = ( γ w h w + γ sat L ) A

γ s + eγ w
γ sat = (1 − n) γ s + nγ w =
1+ e
h
i=
 L
Fh
γ 1  γs 
 ic = sat − 1 =  − 1 
Fb γw 1+ e  γw 

1  ρs 
ic =  − 1 
1 + e  ρw 
23
10 La détermination du gradient hydraulique critique

Ecoulement ascendant:


−ρw gVs

 
jVs 
jVs
( ρsat − ρw ) gVs

ρsat gVs

jc = γ w i c
jc = γ sat − γ w
γ sat
ic = −1
γw

24
11 Le renard hydraulique

Quand le gradient hydraulique critique est atteint, le sol se trouve dans un état
d’équilibre précaire. Dès que cette valeur est dépassée, le sol entre dans un état
de boulance et ne peut plus porter de charge.

Le gradient hydraulique critique dépend de la densité relative du sol et de son indice

des vides: plus l’indice des vides est élevé, plus le gradient hydraulique critique est
faible et plus ce sol est susceptible d’atteindre rapidement l’état de boulance.

1  ρs 
ic =  − 1
1 + e  ρw 

Ce même phénomène peut se produire dans la nature si le gradient hydraulique

augmente pour une quelconque raison et atteint sa valeur critique, le courant d’eau
soulève les grains et le milieu perd toute consistance au chargement comme si
il était liquéfié. On appelle ce phénomène renard hydraulique (nom qui souligne
son caractère sournois).

25
 11 Le renard hydraulique

Palplanche
Perte de charge
totale
Réservoir
Réservoir

Profondeur
d’ancrage Zone de renard
hydraulique

Sable

Couche imperméable

Pour prévenir le phénomène de renard hydraulique, Terzaghi recommande

un facteur de sécurité de 3:
ic
i≤
3

On verra les prescriptions de l'Eurocode 7 dans ce domaine en Géotechnique II. 26

11 Le renard hydraulique

Palplanche

Filtre
Réservoir
Zone de renard
hydraulique

Filtre

Palplanche
Ligne d’écoulement

27
11 Le renard hydraulique

Ligne d’écoulement

Filtre drain Filtre drain

Ligne d’écoulement

Tablier en béton
Tablier en béton

Ligne d’écoulement
28
12 Profondeur critique d’une excavation

Des conditions de boulance peuvent se rencontrer dans un dépôt de sol argileux

agissant comme un écran étanche pour emprisonner une nappe captive.
Si on creuse cette argile à une profondeur critique, le fond de l’excavation
commencera à se soulever sous l’effet des forces d’infiltration découlant de
la perte de charge. Pour déterminer cette profondeur critique, on exprime
la condition d’équilibre critique:

ρt gL = ρw gh w
∆h c
E hw
h w = L + ∆h c H
L Argile
∆h c ρsat
= − 1 = ic
L ρw
Sable Nappe captive

L = H−E ∆h c h
E =H− =H− w
ic 1 + ic
29
13 Les filtres granulaires

L’eau s’écoule à travers le sol sous l’effet d’une charge hydraulique. Si la

charge s’accroît suffisamment, la vitesse d’écoulement peut causer la
migration des particules fines de sol et ainsi mener à une érosion interne.

La stabilité du sol risque alors d’être sérieusement compromise, ce qui conduira

éventuellement à la dégradation des ouvrages supportés par ce sol.

L’érosion interne se produit surtout dans les sables fins et les silts.

Pour freiner l’érosion interne des sols, on a recours à deux types de filtres:
- les filtres granulaires;
- les géotextiles.

Ces filtres bloquent les particules fines tout en laissant passer l’eau. Le plus
souvent, on les place directement entre le sol qu’on doit protéger et le sol
migrateur, dont on doit maîtriser l’érosion interne.

30
13 Les filtres granulaires

Les filtres granulaires sont constitués d’une ou plusieurs couches de sols
pulvérulents, habituellement du sable et parfois du gravier. Chaque couche
doit mesurer au moins 150mm d’épaisseur. Les vides entre les particules
doivent être assez petits pour retenir les particules du sol migrateur, mais
suffisamment grands pour faciliter l’écoulement de l’eau.

Pour concevoir un filtre granulaire, on doit nécessairement connaître la
granulométrie des deux sols concernés par l’érosion interne. On détermine
d’abord les diamètres effectifs limites D15 et D50 du filtre à l’aide des critères
suivants:
D15,filtre
≤5 - empêcher le passage des
Rétention des particules: D85, sol migrateur
particules du sol migrateur;
D50,filtre - éviter le colmatage du filtre
≤ 25
D50, sol migrateur par les particules du sol migrateur.

Perméabilité: D15,filtre - assurer la plus grande

≥5 perméabilité possible du filtre
D15, sol migrateur - diminuer les pressions
hydrostatiques dans le filtre;
- assurer la stabilité du filtre.
31
13 Les filtres granulaires

En se référant aux mêmes critères, on vérifie ensuite la stabilité du filtre. Pour ce

faire, on effectue de nouveaux calculs, mais en considérant que le filtre agit
comme sol migrateur, tandis que le sol à protéger joue le rôle de filtre.

D15, sol à protéger Filtre Sol migrateur

≤5
D85, filtre
D50, sol à protéger
≤ 25
D50, filtre

Ecoulement
D15, sol à protéger
≥5
D15, filtre

Sol à protéger
d'une contamination
32
13 Les filtres granulaires
Pour terminer, on établit la courbe granulométrique du filtre, laquelle doit se situer
dans la zone délimitée par les deux opérations précédentes et être la plus
parallèle possible à la courbe du sol migrateur.

33
13 Les filtres granulaires

Si les dimensions des particules des deux sols concernés sont trop différentes, il
faut concevoir un filtre constitué de plusieurs couches de sol. La figure de la
diapositive 32 présente un filtre à deux couches. Dans ce cas, on devra s’assurer
de la stabilité de chaque couche de sol par rapport à la couche située en aval.

Pendant longtemps, on n’a disposé que des filtres granulaires pour faire échec à
l’érosion interne des dépôts de sol. Mais depuis qu’on utilise les géotextiles, plus
faciles à poser et beaucoup moins coûteux, ce type de filtre tend à disparaître.

34
13 Les géotextiles

Les géotextiles sont des membranes imputrescibles fabriquées à partir de

fibre synthétique comme le polypropylène, le polyester et le nylon; ils
peuvent être tissés ou non tissés. Les géotextiles non tissés sont confectionnés
par aiguilletage, liaison thermique ou liaison chimique.

Les géotextiles agissent comme des barrières perméables: il retiennent les

particules de sol tout en laissant l’eau s’écouler librement. Leur degré de
perméabilité est assez élevé et se compare à celui des sables fins à moyens
à granulométrie serrée. Contrairement aux filtres granulaires, ils sont extensibles
et résistent aux efforts de traction et aux déchirures.

Placés directement entre les sols migrateurs et les sols à protéger, ils conservent
longtemps leurs qualités filtrantes et continueront d’assurer une excellente
protection contre l’érosion interne, même si des déformations importantes
surviennent dans les dépôts de sol.

35
13 Les géotextiles

Filtre géotextile
36
13 Les géotextiles

Géomembrane d’étanchéité bitumineuse

37
13 Les géotextiles

Géotextile anti-érosion (érosion externe)

38