See discussions, stats, and author profiles for this publication at: [Link]

net/publication/354177796

Influence de l’érosion interne sur la stabilité des ouvrages (digues et talus) :

analyse des processus hydromécaniques

Thesis · July 2021

DOI: 10.13140/RG.2.2.19251.55846

CITATION READS

1 1,418

1 author:

Abderrezak Bouziane
Université de Béjaïa
7 PUBLICATIONS 6 CITATIONS

SEE PROFILE

All content following this page was uploaded by Abderrezak Bouziane on 28 August 2021.

The user has requested enhancement of the downloaded file.

 République Algérienne Démocratique et Populaire
Ministère de l’Enseignement Supérieur et de la Recherche Scientifique
Université A. MIRA-BEJAIA

Faculté de technologie
Département de génie civil
Laboratoire de génie de la construction et Architecture

THÈSE
EN VUE DE L’OBTENTION DU DIPLOME DE
DOCTORAT
Domaine : Sciences et Technologies Filière : Génie civil
Spécialité : Génie civil

Présentée par
BOUZIANE Abderrezak

Thème

Influence de l’érosion interne sur la stabilité des ouvrages

(digues et talus) : analyse des processus hydromécaniques

Soutenue le : 08/07/2021 Devant le Jury composé de :

Nom et Prénom Grade

Mr MAZA Mustapha Professeur Univ. de Bejaia Président

Mr BENAMAR Ahmed [Link] Univ. Le Havre Normandie Rapporteur
- France
Mr TAHAKOURT Abdelkader Professeur Univ. de Bejaia Co-rapporteur
Mr KADRI Mohammed Professeur Univ. de Boumerdès Examinateur
Mr CHELOUAH Nasser Professeur Univ. de Bejaia Examinateur
Mr BOUZIDI Mohamed Amin MCA Univ. de Bejaia Examinateur

Année Universitaire : 2020/2021

 Remerciements

En tout premier lieu, je remercie ALLAH, tout puissant, qui m'a donné la capacité, la
patience et le désir d'achever ce travail en dépit de tous les obstacles et contraintes sur la voie
de son achèvement.

Je remercie sincèrement et spécialement mes deux directeurs de thèse Mr BENAMAR

Ahmed et Mr TAHAKOURT Abdelkader pour leur aide constante, leur inspiration et leurs
précieux conseils. Je remercie Mr TAHAKOURT pour le suivi de la thèse, pour sa disponibilité
et ses encouragements. Mes sincères remerciements s’adresse à Mr BENAMAR de m’avoir
accueilli au sein du laboratoire LOMC du Havre où j’ai pu réaliser une grande partie de la thèse,
pour ses critiques pertinentes et ses remarques qui m’ont été indispensables.

Mes sincères remerciements vont aux membres du jury, à Mr MAZA Mustapha pour
m’avoir fait l’honneur de présider le jury. A Mr CHELOUAH Nasser et Mr BOUZIDI
Mohamed Amine pour avoir accepté d'examiner et de commenter cette thèse. A Mr KADRI
Mohammed, pour avoir accepté d'examiner cette thèse ainsi que pour son déplacement à Béjaia.
Je suis très reconnaissant à tous les membres du jury pour avoir accepté de donner de leur temps.

Enfin, je voudrais exprimer ma gratitude à mes parents, mes frères Rafik et Merouane,
ma sœur Asma et à toute ma famille et amis sans exception, pour leurs encouragements et leur
soutien continu, ainsi qu'à tous ceux qui ont contribué à la réalisation de ce projet.

Merci à tous

i
Dédicaces

Je dédie ce travail :

A mes chers parents ;

A mes chers frères et à ma chère sœur ;

A ma chère fiancée et toute sa famille ;

A toute la famille BOUZIANE et BELHADDAD ;

A tous mes amis

ii
Table des matières

Remerciements ................................................................................................................ i

Dédicaces ........................................................................................................................ ii

Table des matières ......................................................................................................... iii

Liste des Figures ............................................................................................................ ix

Liste des tableaux ........................................................................................................ xiii

Introduction générale ...................................................................................................... 1

Chapitre I : Stabilité des ouvrages en terre : Généralités

I. Introduction .............................................................................................................. 4

II. Les ouvrages en terre............................................................................................ 4

III. Mécanismes de rupture des ouvrages en terre ...................................................... 6

III.1 Rupture par glissement ................................................................................. 6

III.2 Rupture par érosion externe .......................................................................... 7

III.3 Rupture par érosion interne ........................................................................... 9

Mécanismes d’érosion interne .................................................................. 9

Problème d’érosion interne dans les ouvrages en terre ........................... 11

IV. Critères d’instabilité ........................................................................................... 12

IV.1 Critères géométriques ................................................................................. 12

Critère de Kezdi [23] .............................................................................. 12

Critère de Kenney et Lau [31] ................................................................ 13

Critère de Li et Fannin [28] .................................................................... 15

Critère de Istomina [24] .......................................................................... 16

Critère de Burenkova [41] ...................................................................... 16

Critère de Mao [26] ................................................................................ 17

Critère de Wan et Fell [43] ..................................................................... 17

iii
 IV.2 Résumé des critères les plus utilisés ........................................................... 19

IV.3 Critères hydrauliques .................................................................................. 20

Critère de Terzaghi ................................................................................. 20

Critère de Skempton et Brogan [47] ....................................................... 21

Critère de Wan et Fell [27] ..................................................................... 22

Résumé sur les critères d’instabilité ....................................................... 24

V. Stabilité des pentes ............................................................................................. 24

V.1 Equilibre plastique limite................................................................................ 25

V.2 Méthode d’analyse de la stabilité des pentes .................................................. 26

Méthodes basées sur l'équilibre du moment global ................................. 26

Méthodes basées sur l'équilibre global des moments et des forces ......... 26

Méthodes basées sur l'équilibre des forces de chaque tranche ................ 27

Méthodes basées sur l'équilibre des moments et des forces de chaque

tranche 27

V.3 Méthode des tranches ..................................................................................... 27

Détermination du facteur de sécurité ....................................................... 29

Méthode ordinaire de Fellenius ............................................................... 30

Méthode simplifiée de Bishop ................................................................. 30

Méthode de Spencer ................................................................................. 31

Méthode simplifiée de Janbu ................................................................... 31

Méthode de Morgenstern Price ................................................................ 32

Chapitre II : Revue bibliographique

I. Introduction ............................................................................................................ 34

II. Influence de l’érosion interne sur le comportement mécanique des sols ........... 34

III. Influence de l’érosion interne sur la stabilité ..................................................... 36

IV. Les approches de modélisation de l’érosion interne .......................................... 42

IV.1 Approches numériques de l’érosion interne ............................................... 42

iv
 Approche continue .................................................................................. 42

Approche discrète ................................................................................... 45

V. Instabilités (ruptures) de pentes liées à l’érosion interne ................................... 48

VI. Conclusion .......................................................................................................... 51

Chapitre III : Etude expérimentale

I. Introduction ............................................................................................................ 52

I. Matériaux et méthode ............................................................................................ 52

II. Description des essais......................................................................................... 53

III. Matériaux étudiés ............................................................................................... 53

III.1 Reconstitution des matériaux ...................................................................... 54

III.2 Préparation des échantillons ....................................................................... 57

III.3 Essai de cisaillement à la boite ................................................................... 57

Principe de l’essai de cisaillement .......................................................... 58

Description de la boite de Casagrande .................................................... 59

III.4 Essai érosion-cisaillement ........................................................................... 61

Boite de Casagrande modifiée pour l’érosion interne ............................ 62

IV. Evaluation de la susceptibilité à l’érosion .......................................................... 62

V. Résultats et discussion ........................................................................................ 63

V.1 Avant érosion .................................................................................................. 64

V.2 Effet de la teneur en fine sur la cohésion et l’angle de frottement interne ..... 65

V.3 Résultats post-érosion ..................................................................................... 66

Conclusion .................................................................................................................... 70

Chapitre IV : Etude numérique : Modèle 1D

I. Introduction ............................................................................................................ 71

II. Modélisation du processus d’érosion interne ..................................................... 71

III. Ecoulement dans les milieux poreux .................................................................. 72

IV. Théorie des milieux poreux ................................................................................ 72

v
 IV.1 Processus d’érosion interne basé sur le continuum ..................................... 73

IV.2 Modélisation 1D de l’érosion interne.......................................................... 74

IV.3 Formulation mathématique ......................................................................... 74

IV.4 Equations d’équilibre de masse................................................................... 76

IV.5 Dégradation de la cohésion ......................................................................... 80

IV.6 Paramètres d’entrée ..................................................................................... 80

IV.7 Conditions aux limites ................................................................................ 81

V. Résultats et discussions ...................................................................................... 81

V.1 Etude paramétrique ......................................................................................... 81

Effet du coefficient d’érosion .................................................................. 81

Effet de la charge hydraulique ................................................................. 84

Effet de longueur de l’échantillon............................................................ 86

Effet de la porosité critique ncr................................................................. 88

V.2 Confrontation des résultats expérimentaux au modèle numérique ................. 90

V.3 Géométrie 1D pour l’analyse de l’érosion interne .......................................... 90

V.4 Simulation numérique..................................................................................... 90

V.5 Paramètres d’entrée du modèle....................................................................... 91

V.6 Calcul de la porosité initiale et critique des différents échantillons ............... 91

Porosité initiale ........................................................................................ 91

Porosité critique ....................................................................................... 92

VI. Résultats numériques.......................................................................................... 93

VI.1 Résultats de simulation ............................................................................... 93

VII. Conclusion .................................................................................................... 100

Chapitre V : Etude numérique : Modèle 2D

I. Introduction .......................................................................................................... 102

II. Formulation du modèle 2D .............................................................................. 102

II.1 Dégradation des paramètres de résistance du sol ......................................... 103

vi
 II.2 Implantation des équations du modèle dans comsol multiphysics ............... 104

II.3 Type d’étude ................................................................................................. 104

II.4 Les modules utilisés...................................................................................... 104

L’interface de darcy (dl) dans une étude transitoire (darcy’s law (dl)) .. 105

L’interface (PDE) de la forme « coefficients » (coefficient form pde(c))

avec la variable « n » ...................................................................................................... 105

L’interface (PDE) de la forme « coefficients » (coefficient form pde(c))

avec la variable « c » ...................................................................................................... 106

III. Exemple d’application...................................................................................... 106

III.1 Choix des points de mesure ...................................................................... 108

III.2 Maillage du domaine................................................................................. 108

III.3 Conditions aux limites .............................................................................. 108

III.4 Résultats et discussions ............................................................................. 109

III.5 Variation temporelle des différents paramètres au pied de la digue ......... 115

IV. Etude paramétrique .......................................................................................... 118

IV.1 Effet du coefficient d’érosion ................................................................... 118

IV.2 Effet de la porosité critique ....................................................................... 121

Effet de la charge hydraulique .............................................................. 124

V. Influence de l’érosion interne sur la stabilité vis-à-vis au glissement.............. 127

V.1 Analyse de stabilité ....................................................................................... 127

V.2 Principes fondamentaux de l'équilibre limite ............................................... 128

Bases de la méthode ............................................................................... 128

V.3 Méthodes de calcul du facteur de sécurité .................................................... 129

V.4 Paramètres de résistance au cisaillement ...................................................... 129

Cohésion ................................................................................................ 129

Angle de frottement interne ................................................................... 129

V.5 Cohésion et d’angle de frottement interne .................................................... 130

V.6 Evaluation du facteur de sécurité .................................................................. 132

vii
VI. Conclusion ........................................................................................................ 136

Conclusion générale.................................................................................................. 138

viii
 Liste des Figures

Figure I. 1 Barrage en enrochement de Gouhou après la brèche [5] .......................................... 5

Figure I. 2 Vue du barrage de Teton après la rupture [6] ........................................................... 5
Figure I. 3 Vue aérienne du barrage de Tangjiashan après rupture ............................................ 6
Figure I. 4 Deux types de glissement : glissement peu profond et glissement profond [10] ..... 7
Figure I. 5 Rupture d’un barrage en terre par surverse [11] ....................................................... 8
Figure I. 6 Riprap : (a) Renforcé et (b) Mélange sol-ciment [12] .............................................. 8
Figure I. 7 Illustration des quatre modes d'initiation de l'érosion interne : érosion interne
régressive, érosion interne de trou ou de fissures, érosion interne de contact entre deux sols et
la suffusion. .............................................................................................................................. 10
Figure I. 8 Les facteurs affectant l’initiation de l’érosion interne [21] .................................... 11
Figure I. 9 Illustration du critère de Kézdi et celui de Kenney et Lau [29, 41] ....................... 13
Figure I. 10 Critère de Kenney et Lau [29] .............................................................................. 14
Figure I. 11 Critère de Li [48] .................................................................................................. 15
Figure I. 12 Limite de stabilité du modèle capillaire [48] ........................................................ 16
Figure I. 13 (a) Contours de la probabilité d'instabilité interne pour les sols limon-sable-gravier
et argile-limon-sable-gravier à teneur en argile et à plasticité limitées ; (b) Contours de la
probabilité d'instabilité interne pour les sols sable-gravier avec moins de 10% de fines non
plastiques passant 0,075 mm [42] ............................................................................................ 18
Figure I. 14 Méthode alternative pour l'évaluation de l'instabilité interne des sols silt-sable-
gravier à granulométrie étalée .................................................................................................. 19
Figure I. 15 Relation entre le gradient critique ic et l’indice de stabilité (H/F)min [46] ............ 22
Figure I. 16 Gradient hydraulique moyenne à laquelle débute l'érosion de particules fines dans
l'échantillon (istart) en fonction du rapport de stabilité minimum (H/F) min [27] .................... 23
Figure I. 17 Gradient hydraulique moyenne à laquelle débute l'érosion de particules fines dans
l'échantillon (istart) en fonction du coefficient d’uniformité Cu [27] ........................................ 23
Figure I. 18 Gradient hydraulique moyenne à laquelle débute l'érosion de particules fines dans
l'échantillon (istart = ic) en fonction de la porosité [27] .......................................................... 24
Figure I. 19 Forces agissantes pour la méthode des tranches appliquée à une surface de
glissement composite ............................................................................................................... 28
Figure I. 20 Variation fonctionnelle de la direction de la force latérale par rapport à la direction
x [67] ........................................................................................................................................ 33
ix
Figure I. 21 Désignation de la force latérale pour la méthode Morgenstern-Price .................. 33

Figure II. 1 schéma du dispositif expérimental utilisé par Chen et al. (2016) [71]. .... 35
Figure II. 2 Résistance du cône avant et après l'érosion interne [72] ........................... 36
Figure II. 3 Schéma de la digue sur fondation d'Orléans (France) [74]. ...................... 37
Figure II. 4 (a) Schéma du modèle physique de Horikoshi et Takahashi (b) Répartition
de la variation spatiale du pourcentage des particules fines [75] ............................................. 38
Figure II. 5 Installation expérimentale de van Beek et al (77) : (a) petite échelle ; (b)
moyenne échelle ; (c) grande échelle ....................................................................................... 39
Figure II. 6 Mise en place des deux canaux par Hu et al. [77] : (a) Infiltration d'eau par
le sommet de la pente ; (b) Augmentation du niveau d'eau à la frontière hydraulique supérieure
.................................................................................................................................................. 40
Figure II. 7 simulation de l’évolution du facteur de sécurité pour une pente de 60° [78]
.................................................................................................................................................. 41
Figure II. 8 Série chronologique typique d'effondrement dû à l'érosion interne [78] .. 41
Figure II. 9 Modèle éléments finis de l’exemple illustratif [87] .................................. 44
Figure II. 10 Modèle numérique de Zou et al. [99] ..................................................... 46
Figure II. 11 Migration des particules fines lors de la suffusion pour différents instants
: (a) 0s, (b) 2s, (c) 5s, (d) 30 s [96] ........................................................................................... 47
Figure II. 12 Prédiction du changement du facteur de sécurité par le modèle BMM [92].
.................................................................................................................................................. 47
Figure II. 13 Evolution du processus d’érosion selon le modèle de Gu et al. [93] ...... 48
Figure II. 14 Coupe transversale de la terrasse fluviale telle qu'elle a été obtenue à partir
des investigations sur le terrain [120] ...................................................................................... 50
Figure II. 15 Processus de rupture du barrage de Gouhou en Chine [124] .................. 51

Figure III. 1 Courbes granulométrique des matériaux de base ..................................... 55

Figure III. 2 Courbes granulométriques des matériaux à base de Kaolinite et de Sable
NE34 (Sable#1) ........................................................................................................................ 55
Figure III. 3 Courbes granulométriques des matériaux à base de Limon et de Sable NE34
(Sable#1) .................................................................................................................................. 56
Figure III. 4 Courbes granulométriques des matériaux à base de Sable fin (Sable#2) et
de Sable NE34 (Sable#1) ......................................................................................................... 56

x
 Figure III. 5 Les trois principaux types de sollicitations .............................................. 58
Figure III. 6 Principe de l’essai de cisaillement directe ............................................... 59
Figure III. 7 Appareil de cisaillement........................................................................... 60
Figure III. 8 Schéma illustratif du dispositif expérimental de l’essai d'érosion-
cisaillement............................................................................................................................... 62
Figure III. 9 Résultat de l’essai de cisaillement avant érosion pour le mélange de
Kaolinite et sable NE 34 ........................................................................................................... 64
Figure III. 10 Résultat de l’essai de cisaillement avant érosion pour le mélange de Limon
et sable NE 34 .......................................................................................................................... 65
Figure III. 11 Résultat de l’essai de cisaillement avant érosion pour le mélange de sable
fin (Sable#2) et sable NE 34 .................................................................................................... 65
Figure III. 12 Effet de la teneur en fines sur la cohésion.............................................. 66
Figure III. 13 Effet de la teneur en fines sur l’angle de frottement interne .................. 66
Figure III. 14 Résultat de l’essai de cisaillement après érosion pour le mélange de
Kaolinite et sable NE 34 ........................................................................................................... 67
Figure III. 15 Résultat de l’essai de cisaillement après érosion pour le mélange de Limon
et sable NE 34 .......................................................................................................................... 67
Figure III. 16 Évolution de la réduction de la cohésion en fonction de la teneur en
amendes .................................................................................................................................... 69

Figure IV. 1 Schéma représentatif de l’état initial et après érosion ............................. 74

Figure IV. 2 Effet du coefficient d’érosion sur les différents paramètres : (a)
Concentration, (b) Porosité, (c) Cohésion. ............................................................................... 83
Figure IV. 3 Effet de la charge hydraulique sur les différents paramètres : (a)
Concentration, (b) Porosité, (c) Cohésion. ............................................................................... 85
Figure IV. 4 Effet de la longueur sur les différents paramètres : (a) Concentration, (b)
Porosité, (c) Cohésion. ............................................................................................................. 87
Figure IV. 5 Effet de la porosité critique sur les différents paramètres : (a) Concentration,
(b) Porosité, (c) Cohésion......................................................................................................... 89
Figure IV. 6 Schéma explicatif du modèle 1D ............................................................. 90
Figure IV. 7 Profil de la pression à la fin de l’érosion pour les différents matériaux .. 94
Figure IV. 8 Variation temporelle de la porosité au milieu des différents échantillons
.................................................................................................................................................. 96

xi
 Figure IV. 9 Profil linéaire de la porosité à la fin de l’érosion pour les différents
matériaux .................................................................................................................................. 97
Figure IV. 10 Variation temporelle de la concentration au milieu de l’échantillon pour
les différents matériaux ............................................................................................................ 97
Figure IV. 11 Profil linéaire de la perméabilité à la fin de l’érosion pour les différents
matériaux .................................................................................................................................. 98
Figure IV. 12 Variation temporelle de la vitesse de Darcy au milieu de l’échantillon pour
les différents matériaux ............................................................................................................ 98
Figure IV. 13 Variation temporelle de la cohésion au milieu de l’échantillon pour les
différents matériaux .................................................................................................................. 99
Figure IV. 14 Variation temporelle de la masse volumique du mélange au milieu de
l’échantillon pour les différents matériaux ............................................................................... 99
Figure IV. 15Variation temporelle de la masse viscosité de l’eau au milieu de
l’échantillon pour les différents matériaux ............................................................................. 100

Figure V. 1 coupe transversale et caractéristiques géométrique de l’exemple étudié 107

Figure V. 2 Maillage du domaine ............................................................................... 108
Figure V. 3 Les conditions aux limites ....................................................................... 109
Figure V. 4 Profil spatial de la pression interstitielle à la fin du processus d’érosion 109
Figure V. 5 Distribution spatiale de la porosité à différents instants ......................... 110
Figure V. 6 Position de la ligne de saturation avant et après érosion ......................... 111
Figure V. 7 Distribution spatiale de la teneur en fines à différents instants ............... 112
Figure V. 8 Distribution spatiale de la perméabilité à différents instants .................. 113
Figure V. 9 Distribution spatiale de la concentration à différents instants ................ 114
Figure V. 10 Distribution spatiale de la cohésion à différents instants ...................... 114
Figure V. 11 Distribution spatiale du champ de vitesse de Darcy à différents instants
................................................................................................................................................ 115
Figure V. 12 Evolution temporelle de la porosité au pied de l’ouvrage..................... 116
Figure V. 13 Evolution temporelle de la perméabilité au pied de l’ouvrage ............. 116
Figure V. 14 Evolution temporelle de la teneur en fines au pied de l’ouvrage .......... 117
Figure V. 15 Evolution temporelle de la concentration au pied de l’ouvrage ............ 117
Figure V. 16 Evolution temporelle de la vitesse de Darcy au pied de l’ouvrage ....... 118
Figure V. 18 Etude paramétrique du coefficient d’érosion ........................................ 121

xii
 Figure V. 19 Etude paramétrique de la porosité critique : (a)(b)(c)(d)(e)(f) .............. 124
Figure V. 20 Etude paramétrique de la charge hydraulique ....................................... 126
Figure V. 21 Discrétisation en tranche et les forces agissant sur une tranche ............ 129
Figure V. 22 Profil de la cohésion et de l’angle de frottement utilisés dans l’analyse de
stabilité pour différentes valeurs de λ et à différents instants : (a) λ=1, (b) λ=1.5, (c) λ=2... 132
Figure V. 23 Surface de glissement critique pour plusieurs instants : (a) à l’état initial
(b) à 500 jours (c) à 1500 jours (d) à 2500 jours (e) à 3600 jours .......................................... 134
Figure V. 24 Comparaison de la réduction du facteur de sécurité par les différentes
méthodes................................................................................................................................. 134
Figure V. 25 Variation temporelle du facteur de sécurité pour différentes valeurs de λ
................................................................................................................................................ 136

Liste des tableaux

Tableau III. 1 Caractéristiques physicochimiques et composition chimique du sable de

fontainebleau (NE34) ............................................................................................................... 56
Tableau III. 2 Caractéristiques physicochimiques et composition chimique des fines
utilisés....................................................................................................................................... 57
Tableau III. 3 Évaluation de la vulnérabilité des spécimens à l'érosion interne selon les
différents critères ...................................................................................................................... 63
Tableau III. 4 Récapitulatif des résultats de l'évolution de la cohésion et de l'angle de
frottement interne dû à l'érosion interne ................................................................................... 68

Tableau IV. 1 Paramètres d’entrée du modèle ............................................................. 80

Tableau IV. 2 Paramètres d’entrée du modèle d’érosion ............................................. 91
Tableau IV. 3 valeurs de λ et ncr pour les différentes simulations ............................... 92

Tableau V. 1 Paramètres d’entrée du modèle 2D ....................................................... 107

Tableau V. 2 Paramètres d’entrée pour évaluer la stabilité après érosion .................. 130

xiii
Tableau V. 3 Variation temporelle du facteur de sécurité .......................................... 135

xiv
Introduction générale
Introduction générale

Dans l'antiquité, les barrages en terre étaient construits dans le seul but de fournir de
l'eau ou d'irriguer. Avec le développement des civilisations, les besoins en matière
d'approvisionnement en eau, d'irrigation, de contrôle des inondations, de navigation, de qualité
de l'eau, de contrôle des sédiments et d'énergie se sont accrus. C'est pourquoi les barrages sont
construits pour répondre aux besoins de l'approvisionnement en eau, la lutte contre les
inondations, l'irrigation, la navigation, le contrôle de la sédimentation et l'énergie
hydroélectrique. Un barrage est la pierre angulaire du développement et de la gestion des
ressources en eau d'un bassin fluvial. Le barrage polyvalent est un projet très important pour
les pays en développement, car la population tire des avantages domestiques et économiques
d'un seul investissement.

Les barrages et les digues sont souvent construits dans les régions montagneuses mais
aussi dans les vallées près des villes. Si une érosion interne se produit dans le remblai, elle peut
constituer une menace énorme pour la sécurité de la population en aval de l’ouvrage. C'est
pourquoi la recherche sur l'érosion interne des ouvrages en terre est stratégique pour la
protection des vies humaines et des biens.

Ces ouvrages, qui sont utilisées pour prévenir les inondations, fournir de l'énergie
hydroélectrique ou stocker l'eau pour la consommation ou l'irrigation, sont construits dans des
vallées à la confluence des cours d’eau. Les êtres humains ont une longue histoire de
construction de barrages, et les plus anciennes structures remontent à 3000 avant J.-C. en
Mésopotamie, au Moyen-Orient et dans l'Égypte ancienne. Avec le développement de la
science et de la technologie, la construction dans ce domaine a également fait des progrès
considérables. Grâce à la construction de barrages, les ressources en eau peuvent être
pleinement utilisées pour résoudre les problèmes de pollution engendrés par l'énergie fossile.
Ainsi, au cours du siècle dernier, des barrages ont été construits à grande échelle dans le monde
entier.

A la suite de ruptures de digues, un nombre croissant de chercheurs ont mis en évidence

les risques d’instabilité qui peuvent être associés aux barrages affectés par l'érosion interne, et
se sont ensuite concentrés sur l’analyse du processus et l’évaluation des risques. L'érosion
interne est un phénomène complexe qui demeure l'une des raisons les plus fréquentes de la

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 1

défaillance des digues fluviales et des barrages en terre. Selon le comité Français des grands
barrages (CFGB) [1], l’érosion interne est définie par l’arrachement et le transport des particules
fines d’un sol sous l’effet d’un écoulement d’eau le traversant. Elle est considérée comme l’une
des causes principales de la rupture des digues et des barrages en terre. Dans une étude de Foster
et al. [2], parmi les 11192 barrages étudiés, 136 (1.2 %) ont révélé des dysfonctionnements,
environ 46 % ont présenté une érosion interne, 48 % une surverse et 5,5 % un glissement. Ces
statistiques confirment que l'érosion interne est une cause importante de défaillance.
Aujourd'hui, de nombreuses questions subsistent sur l'érosion interne, en particulier son
influence sur les propriétés mécaniques des sols, notamment les propriétés de résistance au
cisaillement.

La présente étude tente d’apporter une contribution à la conception des ouvrages

d’ingénierie géotechnique et hydraulique en mettant l’accent sur les digues en terre et les talus
naturels. Dans ce travail, on s’intéresse à l’effet de l’érosion interne sur les paramètres de
résistance au cisaillement des sols, et par conséquent, sur la stabilité des ouvrages en terre. Il
comporte deux parties principales : (i) une étude expérimentale au laboratoire de l’érosion
interne et ses conséquences sur la résistance au cisaillement des sols, cette partie est réalisée à
l’aide de la boite de Casagrande modifiée ; (ii) Un modèle numérique basé sur la théorie des
milieux poreux, en considérant un milieu composé de trois phases, a été ensuite formulé. Ce
modèle a été implémenté dans un logiciel de simulation dans deux conditions d’écoulement :
(i) modèle unidimensionnel où l’effet de la dispersion a été négligé, (ii) modèle bidimensionnel
incluant l’effet de la dispersion, qui a servi ensuite à la simulation du processus de l’érosion
interne à l’échelle de l’ouvrage. En tenant compte de l’érosion interne dans un ouvrage
hydraulique (digue), le facteur de stabilité de cette dernière a été mesuré à différentes étapes.
Les calculs numériques effectués tentent de fournir une idée sur l’instant où l’érosion interne
peut conduire à la rupture de cet ouvrage.

Le manuscrit de thèse est scindé en cinq chapitres, en plus de l’introduction générale et

de la conclusion générale. Ces chapitres se succèdent comme suit :

Introduction générale : Elle consiste en une introduction à la problématique de la présente

étude. Elle aborde globalement, l'érosion interne et de ses conséquences sur les structures en
terre, les insuffisances des recherches actuelles, et les objectifs de la présente étude sont
présentés.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 2

Chapitre 1 (Stabilité des ouvrages en terre : Généralités) : Ce chapitre présente une
description brève et non exhaustive des différents mécanismes de l'érosion interne. Ensuite, les
différents critères d’évaluation de la stabilité interne des sols ont été passés en revue. Ces
critères ont été utilisé pour évaluer le potentiel d’érosion des sols étudiés dans le cadre de cette
thèse. Les différentes méthodes d’évaluation de la stabilité des pentes, utilisées dans cette étude,
ont été ensuite présentées.

Chapitre 2 (Revue bibliographique) : Une revue de littérature en termes d'investigations

expérimentales et numériques sur l'érosion interne est fournie dans ce chapitre.

Chapitre 3 (Etude expérimentale) : Ce chapitre présente l’étude expérimentale effectuée au

laboratoire sur l’effet de l’érosion interne sur les paramètres de résistance au cisaillement des
sols. Des essais séquentiels d’érosion-cisaillement ont été effectués sur un ensemble
d’échantillon confectionnés à partir des mélanges de sable de fontainebleau (France) avec
différentes fractions d’argile ou de limon.

Chapitre 4 (Etude numérique : Modèle 1D) : Il décrit un modèle mathématique

unidimensionnel pour la simulation de l’érosion interne et sa validation. Le modèle d’érosion
est formulé à l'échelle du volume élémentaire représentatif (VER) pour simuler le transport de
particules dans les milieux poreux. Une confrontation aux résultats expérimentaux est présentée
afin de déterminer le coefficient d’érosion λ.

Chapitre 5 (Etude numérique : Modèle 2D) : Dans ce chapitre, nous présentons une
simulation numérique de l’érosion interne dans un écoulement bidimensionnel et ses
conséquences sur la stabilité d’un ouvrage en terre. Le modèle d’érosion est formulé à l'échelle
du volume élémentaire représentatif (VER) pour simuler le transport de particules au sein d’une
digue. Les paramètres du modèle ont été évalués dans le temps et dans l’espace et l’influence
de l’érosion interne sur la stabilité de l’ouvrage a été évaluée en reliant l’endommagement et le
transfert de masse par un paramètre de couplage qui est la porosité du milieu. Ce dernier
chapitre présente ainsi une étude sur la prédiction de la dégradation de la stabilité d’une pente
naturelle due à l’érosion interne. La stabilité globale de l’ouvrage est évaluée et le facteur de
sécurité est calculé en considérant la variation temporelle de la distribution spatiale de la
cohésion et de l’angle de frottement interne.

Conclusion générale : Enfin, une conclusion générale résume les principaux résultats de ce
travail de thèse et fournit des perspectives pour les futurs travaux de recherche sur le sujet.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 3

Chapitre I : Stabilité des
ouvrages en terre
(Généralités)
 Chapitre I Stabilité des ouvrages en terre : Généralités

I. Introduction

Ce chapitre passe en revue les différents mécanismes de l'érosion interne, à savoir

l'érosion interne de trou, l’érosion régressive, l’érosion de contact et la suffusion. À partir de
plusieurs études expérimentales, les chercheurs ont identifié les principaux paramètres de
l’initiation et de l'évolution de l'érosion interne, ainsi, ils ont proposé des critères pour
l’évaluation de la susceptibilité à l’érosion d’un sol donné, les critères les plus répandu dans la
littérature ont été présenté dans ce chapitre. Ces critères seront ensuite utilisés au chapitre 3
pour l’évaluation de la susceptibilité à l’érosion interne des matériaux étudiés dans le cadre de
cette étude. Les différentes méthodes d’évaluation de la stabilité des pentes, utilisées dans cette
étude, ont été présentées dans ce chapitre.

II. Les ouvrages en terre

Un barrage en remblai est une barrière permettant de retenir l'eau artificiellement, tandis
qu'un barrage en glissement de terrain désigne le blocage d'un bassin de drainage par des
glissements, des effondrements, des écoulements et des avalanches, naturellement [3, 4]. En
raison de problèmes tels que des matériaux de construction de mauvaise qualité, un contrôle de
qualité insuffisant pendant la construction, un mauvais entretien, un afflux extrême, des
glissements de terrain et des conditions hydrologiques et environnementales changeantes, un
barrage en remblai peut se rompre par débordement ou par érosion interne, ou présenter des
difficultés liées à l'érosion interne. Par exemple, le barrage en enrochement de Gouhou (Chine)
(Figure I. 1) a rompu en 1993 lors du remplissage initial du réservoir, en raison de l'érosion
interne. La rupture a commencé par une fuite dans le barrage par le haut de la paroi en béton.
L'eau s'est écoulée à une altitude élevée sur la pente aval. L'érosion interne s'est ensuite
développée de manière régressive de l'aval vers l'amont, entraînant la rupture du barrage [5].
La Figure I. 2 montre le processus de rupture du barrage de Teton (Etats unis) par « Piping » [6].
En raison de la présence de fractures hydrauliques dans le limon perméable utilisé dans le
noyau, une érosion interne a été initiée dans les fissures par l'écoulement d’infiltration. Après
qu'une grande quantité de matériaux du remblai ait été enlevée par l'écoulement, un "Pipe" s'est
formé. Ensuite, l'érosion progresse assez rapidement jusqu'à ce que la partie des matériaux
située au-dessus du "Pipe" devienne instable et s'effondre [7]. Charles [8] a souligné que

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 4

Chapitre I Stabilité des ouvrages en terre : Généralités

certains barrages en remblai dans les pays européens ont connu une érosion interne, ce qui est
d'une importance cruciale pour l'évaluation de la sécurité à long terme. Les risques potentiels
posés par les ruptures de barrages deviennent de plus en plus graves en raison de la croissance
rapide des villes en aval des ouvrages qui pourraient être touchées par les ruptures de barrages.

Figure I. 1 Barrage en enrochement de Gouhou après la brèche [5]

Figure I. 2 Vue du barrage de Teton après la rupture [7]

Un barrage de retenue en remblai est constitué d'une masse hétérogène de terre et de

roches non consolidées ou mal consolidées. Il ne comporte pas de barrières de protection contre

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 5

 Chapitre I Stabilité des ouvrages en terre : Généralités

les infiltrations, ni de déversoir ou d'exutoire protégé contre les débordements. Le risque de

défaillance est donc relativement élevé et les conséquences qui en découlent sont dévastatrices.
Par exemple, le barrage en terre de Tangjiashan (Chine) (Figure I. 3) déclenché par le
tremblement de terre de Wenchuan en 2008 a menacé plus de 1,3 million de personnes dans les
zones en aval et les villes de Xuanping et de Yancheng ont été inondées [9].

Figure I. 3 Vue aérienne du barrage de Tangjiashan après rupture

III. Mécanismes de rupture des ouvrages en terre

III.1 Rupture par glissement

Un glissement dans un remblai ou dans un sol naturel ou une roche est un mouvement
de masse de matériaux. Les glissements peuvent se produire en raison d'un mauvais compactage
du sol, d'une pente trop raide pour le matériau du talus, d’une infiltration, d'un rabattement
brusque du niveau du réservoir, d'un affaissement du pied du talus ou de la saturation et de
l'affaiblissement du talus ou de la fondation [10].

Les phénomènes de glissement peuvent être classés en deux types : le glissement peu
profond et le glissement profond (Figure I. 4). Le glissement peu profond se situe dans les 2 à 3
premiers mètres sous la surface du remblai. Le glissement peu profond ne menace généralement
pas la sécurité immédiate d’un barrage ce qui n’est pas le cas du glissement profond qui peut

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 6

Chapitre I Stabilité des ouvrages en terre : Généralités

impacter les fondations. Un glissement massif de matériaux peut avoir de lourdes conséquences
pour la stabilité d’un barrage.

Figure I. 4 Deux types de glissement : glissement peu profond et glissement profond [11]

III.2 Rupture par érosion externe

L’érosion externe ou surverse se matérialise par la déstructuration de la surface de
l’ouvrage par arrachement et déplacement des particules de sol dû à l’écoulement de l’eau sur
la surface de l’ouvrage. Il peut se produire lors d’une période de crue lorsque la capacité du
déversoir (ou évacuateur de crue) est insuffisante.

En général, les barrages en remblai ne peuvent pas supporter un débordement important,

en raison de la résistance à l’érosion limitée des matériaux de construction. Les barrages en
béton sont généralement plus résistants au débordement, en raison de la durabilité des matériaux
constitutifs ainsi que de forte résistance à l’érosion de la fondation rocheuse. Cependant, si les
roches de la fondation de barrage en béton sont faibles et fracturées, le support de fondation
peut être très altéré par un débordement et induire une déstabilisation de l’ouvrage.

La rupture d’un barrage en remblai due à un débordement dépend de certains détails

géométriques du barrage. L’érosion commencera aux endroits où il y a un changement de
géométrie ou une discontinuité. Ça pourra être au pied aval du barrage ou au bord aval de la
crête du barrage. La Figure I. 5 montre la progression de l’érosion d’un barrage qui s’initie au
pied du barrage et qui régresse vers l’amont jusqu’à la crête de l’ouvrage.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 7

Chapitre I Stabilité des ouvrages en terre : Généralités

Figure I. 5 Rupture d’un barrage en terre par surverse [12]

Le choix du type de protection dépend des conditions dans lesquelles l’installation doit
être faite, la disponibilité du matériel de protection, le coût et la maintenance de la protection.
Quelques solutions sont proposées pour améliorer la capacité du barrage à lutter contre l’érosion
de surverse :

✓ Agrandissement de la capacité de déversoir ;

✓ Protection moins coûteuse par couverture végétale du versant. L’efficacité du dispositif
dépend des propriétés des plantes, de la surface du barrage et du type de matériaux
constitutifs du barrage. Cependant, les racines des plantes pénètrent dans le barrage et
peuvent provoquer la formation de chemins d’infiltration dans le barrage après leur
dépérissement ;
✓ Protection par installation d’un riprap qui peut être renforcé, en enrochement, constitué
d’un mélange sol-ciment et ou en béton (Figure I. 6).

(a) (b)

Figure I. 6 Riprap : (a) Renforcé et (b) Mélange sol-ciment [13]

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 8

Chapitre I Stabilité des ouvrages en terre : Généralités

III.3 Rupture par érosion interne

Selon Foster et al. (2000)[2], Wan et Fell (2004) [14] et Zhang et al. (2009) [15], 46%
des dégâts observés sur les ouvrages en terre proviennent de l'érosion interne. En outre, Wan
[16] a rapporté que plusieurs barrages en France souffrent de l’érosion interne.

L’érosion interne est définie par Richards et Reddy [17] comme la migration de
particules initiées par les forces érosives de l’eau le long de fissures dans les matériaux cohésifs
ou des vides dans une structure de sol pulvérulent qui provoque une instabilité interne. Kenney
et Lau [18], Fell et al. [19] décrivent l’instabilité interne comme une incapacité pour des
particules plus grosses d’empêcher la migration de particules plus fines emportées par un
écoulement d’infiltration. Selon le comité Français des grands barrages, l’érosion interne se
développe lorsqu’il y a combinaison de deux conditions : l’arrachement des particules et leur
transport. Le phénomène d’érosion interne peut mettre en jeu différents phénomènes qui se
combinent entre eux. Cette dernière peut se produire au sein d’un ouvrage en terre ou dans sa
fondation. Quatre types de mécanisme d’initiation de l’érosion interne dans l’ouvrage
hydraulique sont différenciés par leurs conditions limites.

Mécanismes d’érosion interne

Le processus entier de l’érosion interne peut survenir selon quatre phases : initiation,
continuation, progression et développement. Quant à l’initiation de l’érosion interne, la
communauté scientifique, lors de la rencontre d’Aussois, France en 2005 a retenus quatre
processus d’initiation [20]: (a) Erosion interne régressive, (b) Erosion interne de trou ou de
fissures, (c) Erosion interne de contact entre deux sols, (d) La suffusion. La Figure I. 7 montre
ces quatre mécanismes d’érosion interne.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 9

Chapitre I Stabilité des ouvrages en terre : Généralités

Figure I. 7 Illustration des quatre modes d'initiation de l'érosion interne : érosion interne
régressive, érosion interne de trou ou de fissures, érosion interne de contact entre deux sols et la
suffusion.

Erosion interne régressive : Il s'agit de l'érosion des particules de sol à l'extrémité de

sortie d'un chemin d'infiltration en raison d'une vitesse de sortie élevée. La sortie du chemin
d'infiltration peut se situer sur une surface libre, comme la surface du sol en aval d'une fondation
en terre ou la face aval d'un remblai homogène, ou à l'interface entre le noyau d'argile et le filtre
aval d'un barrage en terre ou en enrochement à noyau central [16]. Le processus se déplace
progressivement vers le côté amont de la digue ou de sa fondation jusqu'à ce qu'un conduit
continu soit formé.

Erosion interne de trou ou de fissure : Le terme « érosion de trou » est utilisé la première
fois par Sherard (1979) [21]. Ce mode d’érosion interne se produit par un écoulement concentré
dans une fissure au sein d’un ouvrage en terre. Lorsqu’elle apparaît une fissure par laquelle
s’écoule de l’eau, cette dernière érode les parois de cette fissure et conduit à la formation d’un
trou. L’apparition des fissures dans les ouvrages en terre peut être due à plusieurs causes telles
que les tassements différentiels, compactage insuffisant des matériaux et la dessiccation.

Erosion interne de contact : Ce type d’érosion interne se manifeste à l’interface d’un

matériau fin en contact avec un matériau grossier. L’écoulement étant parallèle au plan contact,
il érode les particules fines dans et à travers le sol grossier. Ce processus affecte généralement
la surface de contact fondation-digue.

La suffusion : La suffusion est la forme d’érosion interne dans laquelle la matrice

grossière d’un sol ne satisfait pas la condition d’autofiltration [21]. Par conséquent, les
particules fines sont éliminées par les vides entre les particules grossières. Le sol est donc
considéré intérieurement instable.

Garner et Fannin (2010) [22] ont développé un diagramme de Venn, comme le montre
la Figure I. 8, pour illustrer que l'érosion commence lorsque les trois conditions défavorables
suivantes sont remplies: (1) la sensibilité des matériaux (ex : forme de la courbe
granulométrique, plasticité et dispersivité) ; (2) les conditions de contrainte (ex : La géométrie
des fondations et les pratiques de construction sont propices au développement de zones à faible
contrainte dans le remblai) ; et (3) la charge hydraulique (ex : le niveau d'eau remonte jusqu'à
la fissure et la vitesse d'écoulement dans la fissure est suffisante pour déclencher une érosion
de trou) se produisent.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 10

Chapitre I Stabilité des ouvrages en terre : Généralités

Susceptibilité matérielle

•Érosion de
trou
•Érosion de
contact
•Érosion
régressive
Condition
Condition
de
hydraulique
contrainte
critique
critique

Figure I. 8 Les facteurs affectant l’initiation de l’érosion interne [22]

Problème d’érosion interne dans les ouvrages en terre

L'érosion interne provoquée par l’écoulement d’infiltration est un processus progressif
et peut induire les incidents suivants [16] : (a) rupture de la pente de la face aval ; (b) tassement
excessif du barrage avec perte de la revanche ; et (c) formation de trous d'évacuation sur la crête
et en amont du barrage. Dans le premier cas, comme la perte des particules fines se produit à la
sortie de l’écoulement de la face aval, la pression locale de l'eau interstitielle augmente en plus
de la diminution de la résistance au cisaillement du sol ; la stabilité de la pente est alors affectée.
Une rupture locale de la pente peut se produire et le processus se poursuit alors en progressant
vers l'amont. Dans le second cas, en raison de la perte de particules fines dans le noyau ou les
filtres, la matrice grossière peut se réarranger et le volume total diminue et le tassement se
produit. Dans le troisième cas, en raison de la perte de particules fines dans les chemins
d'écoulement très perméables, les particules fines sont emportées le long de ces chemins. Par
conséquent, un changement de volume de la matrice grossière se produit le long de ces trajets
et un trou local se forme.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 11

 Chapitre I Stabilité des ouvrages en terre : Généralités

IV. Critères d’instabilité

Sachant que la phase de l’initiation de l’érosion interne est l’évènement le plus

importante et le plus difficile à évaluer (c'est-à-dire qu'il a tendance à avoir l'impact le plus
potentiel sur la probabilité de défaillance d’un ouvrage en terre), plusieurs chercheurs ont
proposé des critères d’évaluation de la stabilité interne des sols. Il existe principalement deux
types de critères (critères géométriques, critères hydrauliques) [23–29]. Les critères
géométriques sont basés sur la forme de la courbe granulométrique [23–28], tandis que les
critères hydrauliques sont basés sur le gradient hydraulique qui provoque l’initiation du
détachement des particules fines du sol . Cette phase est influencée par les trois principales
conditions citées ci-dessus: conditions géométriques du sol, conditions hydrauliques et
conditions mécaniques [30, 31, 40, 32–39]. Dans ce qui suit nous présenterons un ensemble de
critères les plus utilisé. Cependant, cette liste n’est pas exhaustive et représente essentiellement
les critères utilisés dans cette étude.

IV.1 Critères géométriques

La courbe granulométrique d’un sol donne des informations importantes sur sa structure
et représente sa carte d’identité. Une bonne analyse de cette courbe a permis aux chercheurs de
proposer plusieurs critères d’évaluation de la stabilité interne d’un sol. Bien que la stabilité
interne dépend de plusieurs paramètres tels que la densité, le gradient hydraulique et l’état de
contrainte ; les critères basés sur la courbe granulométrique ont montré leur efficacité dans
plusieurs études [21, 24, 28, 31, 41–44].

Critère de Kezdi [23]

La méthode de Kézdi pour l’évaluation de la susceptibilité d’un matériau à l’érosion
interne consiste à :

1. Diviser la courbe granulométrique du matériau en deux parties, à partir d’un

point (d) dont la position est variable le long de la courbe granulométrique et
arbitrairement choisie sur cette courbe. Le point (𝑑) permet d’obtenir une
fraction grossière (𝐹𝑐 ) et une fraction fine (𝐹𝑓 ). La fraction constituée des
particules grossières forme le squelette du matériau jouant le rôle de filtre. Elle
présente un réseau de pores à l’intérieur duquel se déplace la fraction constituée
de particules fines (Figure I. 9).

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 12

Chapitre I Stabilité des ouvrages en terre : Généralités

2. Déterminer les deux paramètres 𝑑85𝑓 et 𝑑15𝑐 qui sont respectivement, le

diamètre correspondant à 85% de passant sur la courbe granulométrique de la
partie fine, et le diamètre correspondant à 15% de passant sur la courbe
granulométrique de la partie grossière.
𝑑 𝑑
3. La stabilité du matériau est évaluée par le rapport (𝑑15𝑐 ). Si (𝑑15𝑐 ≤ 4), le
85𝑓 85𝑓

matériau est considéré stable.

Sachant que la limite entre la fraction grossière et fine n’est pas exactement connue,
Kézdi propose de balayer sur l’ensemble de la courbe granulométrique.

Figure I. 9 Illustration du critère de Kézdi et celui de Kenney et Lau [31, 42]

Critère de Kenney et Lau [31]

Kenney et Lau [31], ont effectué des essais expérimentaux et ont proposé un critère
permettant d’évaluer la stabilité interne. Seize échantillons de grave sableuse compactés entre
80 et 100% de leur densité maximale sont soumis à un écoulement descendant dans un
perméamètre à charge constante. Les essais ont été effectués à des gradients hydrauliques
inconnus. Le diamètre des échantillons est soit 245 mm ou 580 mm et leur longueur est soit 450
mm ou 860 mm. Une surcharge de 10 kPa et des vibrations manuelles ont été appliquées. La
comparaison des courbes granulométriques avant et après l’essai a permis de constater si
l’instabilité interne a été provoquée dans les matériaux par les différentes sollicitations
appliquées.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 13

Chapitre I Stabilité des ouvrages en terre : Généralités

A partir des résultats obtenus, ils ont un élaboré un critère d’évaluation de la stabilité
interne des sols granulaires dont le diamètre est supérieur à 80 µm. Cette approche consiste à
créer une courbe de forme H= f(F) à partir de la courbe granulométrique. F est le pourcentage
des passants à un diamètre d et H celui des particules dont le diamètre se situe entre d et 4d.
Cette courbe de forme ne dépend pas de la taille des particules mais l’allure de la courbe
granulométrique toute entière. La comparaison de H et F permet de déclarer un sol
potentiellement stable ou instable. Les valeurs de H entre d et 4d donnent le nombre de
particules capables de prévenir le passage des particules de dimensions d. Une insuffisance en
particules dont le diamètre est compris entre d et 4d rend possible la migration des particules
de dimensions d (Figure I. 10).

Figure I. 10 Critère de Kenney et Lau [31]

Si sur une partie de la fraction fine correspondant à F < 0.20 d’une granulométrique
étalée (Cu > 3) ou F < 0.30 pour une granulométrie serrée (Cu < 3), H>F alors il y a plus de
particules dont le diamètre est compris entre d et 4d, susceptibles de bloquer la mobilité des
particules de diamètre d. Dans ce cas, le matériau est classé comme étant stable sous un
écoulement selon Kenney et Lau. Si H<F, le matériau est classé dans ce cas comme instable.

En se basant sur les travaux de Lobochkov, Kenney et Lau ont défini initialement une
limite de stabilité de H/F=1.3. Après discussion avec Milligan [45] et Sherard et Dunnigan
[46], il ressort que cette limite a été trop conservative. Une nouvelle limite de H/F= 1 a donc
été retenue. Par rapport à ces deux limites de stabilité, Skempton et Brogan [47], pensent qu’un
rapport de H/F= 1 correspondrait mieux à des matériaux assez denses tandis qu'un rapport de
H/F=1.3 serait plus adapté aux matériaux plus lâches. Kenney et Lau font remarquer que les
effets combinés des différentes sollicitations auxquelles étaient soumis les échantillons lors des

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 14

Chapitre I Stabilité des ouvrages en terre : Généralités

essais constituent des conditions très sévères qu’on rencontre rarement dans la pratique. Cela
signifie tout simplement qu’il existe un potentiel d’instabilité dans ces sols mais n’implique pas
que l’instabilité s’y produira inéluctablement sur le terrain.

Critère de Li et Fannin [28]

Pour pallier aux imperfections des critères géométriques (Kezdi [23] et Kenney et Lau
[31, 48]) qui sont des critères empiriques utilisant une mesure indirecte de la taille des
constrictions, Li [49] s’est inspiré du critère de Kovacs [50] pour proposer un critère semi-
empirique d’évaluation de stabilité interne. Le critère de Kovacs compare la taille (d) de la
particule à celle des constrictions Dc. Si d<Dc alors la particule est susceptible d’être transportée
à travers le réseau des pores formés par les grosses particules. Le critère de Kovacs présente un
grand intérêt sur les autres par le fait qu’il tient à la fois compte de la courbe granulométrique ;
de la forme des grains et de la porosité du sol.

Figure I. 11 Critère de Li [49]

Une limite entre sols stable et instable a été proposée à l’aide du ratio de constriction
′ ⁄
𝑑85 𝑂50 . L’organigramme ci-dessus illustre cette méthode (Figure I. 11). Si le ratio>0.42, le
′
matériau est stable sinon il est instable (Figure I. 12). Dans ce ration 𝑑85 est le diamètre
correspondant à 85% en masse de la partie fine et 𝑂50, le diamètre moyen d'un tube capillaire.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 15

Chapitre I Stabilité des ouvrages en terre : Généralités

Figure I. 12 Limite de stabilité du modèle capillaire [49]

Critère de Istomina [24]

Estomina [24] a développé un critère simple pour l’évaluation de la stabilité interne basé
𝑑
sur le coefficient d’uniformité 𝐶𝑢 = 𝑑60 , d60 étant le tamis correspondant à 60% du passant
10

cumulé et d10 est le tamis correspondant à 10% du passant cumulé. Il évalue le potentiel de
transport des particules fines à travers les vides intergranulaires de telle sorte que le sol est
considéré stable si 𝐶𝑢 est inférieur à 10, tandis qu’il est instable si 𝐶𝑢 est supérieur à 20. Les
sols avec un coefficient d’uniformité entre 10 et 20 sont considérés en zone de transition.

Critère de Burenkova [41]

La méthode de Burenkova [41] est basé sur 22 essais expérimentaux réalisés sur des
matériaux granulaires. Deux paramètres sont essentiels pour évaluer la stabilité interne par ce
𝑑 𝑑
critère : ℎ′ = 𝑑90 et ℎ′′ = 𝑑90 , avec d15, d60 et d90 sont respectivement les tamis correspondant à
60 15

15%, 60% et 90% du passant cumulé. Le rapport d90/d60 représente la pente de la partie grossière
de la courbe granulométrique. Les valeurs élevées représentent des particules grossières de
taille presque unique qui auront de grands espaces de constriction par rapport à un sol bien
gradué. Le rapport d90/d15 peut être considéré comme une mesure de l'action filtrante entre la
fraction grossière et la fraction fine.

Selon ce critère le sol est considéré stable si ℎ′ et ℎ′′ satisfaisaient l’équation 1.

1 + 0.76 log(ℎ′′ ) ≤ ℎ′ ≤ 1 + 1.86log⁡(ℎ′′ ) (1)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 16

Chapitre I Stabilité des ouvrages en terre : Généralités

Critère de Mao [26]

Selon le critère de Mao [26], le sol est considéré instable vis-à-vis de l’érosion interne
si la relation 2 est satisfaite :

𝑝𝑓 ≥ 100⁄4(1 − 𝑛) (2)
Où pf est la fraction fine (Pour les sols à granulométrie discontinue, pf = le passant au
point de discontinuité (%); Pour les sols bien gradués, pf = le passant (%) correspondant au
diamètre de division df = 1.3√𝑑15 𝑑85 ) ; n représente la porosité du milieu. Li (2008) [51] a
rapporté que cette méthode est très recommandée en chine.

Critère de Wan et Fell [43]

La méthode de Wan et Fell est basée sur les travaux de Burenkova (1993) [41]. Leur
étude est réalisée sur des mélanges de limon-sable-gravier, argile-limon-sable-gravier et sable-
gravier. Deux approches ont été proposées par Wan et Fell, la première concerne à la fois les
sols à granulométrie étalée et discontinue et l’autre ne concerne que les sols à granulométrie
étalée.

a) Méthode pour les sols à granulométrie étalée et discontinue : Les auteurs ont défini
des contours de probabilité d'instabilité interne égale, comme le montre la Figure I. 13. La Figure
I. 13a doit être appliquée aux mélanges limon-sable-gravier et argile-limon-sable-gravier ayant
un indice de plasticité inférieur à 13 % et moins que 10 % de la fraction de taille de l'argile
passant à 0,002 mm, et la Figure I. 13b aux sols sable-gravier avec moins de 10 % de fines de
limon non plastique passant 0,075 mm. Les courbes de probabilité de la Figure I. 13 sont définies
par les équations suivantes :

𝑃 = exp(𝑍) /[1 − exp⁡(𝑍)] (3)

Avec

𝑍 = 2.378 log(ℎ′′ ) − 3.648ℎ′ + 3.701 (4)

Pour la Figure I. 13a

𝑍 = 3.875 log(ℎ′′ ) − 3.591ℎ′ + 2.436 (5)

Pour la Figure I. 13b, où P représente la probabilité d’instabilité interne.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 17

Chapitre I Stabilité des ouvrages en terre : Généralités

(a)

(b)

Figure I. 13 (a) Contours de la probabilité d'instabilité interne pour les sols limon-sable-
gravier et argile-limon-sable-gravier à teneur en argile et à plasticité limitées ; (b) Contours de la
probabilité d'instabilité interne pour les sols sable-gravier avec moins de 10% de fines non plastiques
passant 0,075 mm [43]

b) Méthode alternative pour les sols à granulométrie étalée : cette méthode est
représentée par les rapports d90/d60 et d20/d5. La Figure I. 14 montre les limites des zones
stables et instables. La zone stable peut être définie par les équations suivantes :

30/log⁡(𝑑90 ⁄𝑑60) < 80 (6)

Et

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 18

Chapitre I Stabilité des ouvrages en terre : Généralités

30/log⁡(𝑑90 ⁄𝑑60) < 80 ⁡𝑒𝑡⁡15/log⁡(𝑑20 ⁄𝑑5 ) > 22 (7)

Figure I. 14 Méthode alternative pour l'évaluation de l'instabilité interne des sols silt-sable-
gravier à granulométrie étalée

IV.2 Résumé des critères les plus utilisés

Tableau I. 1 résumé des critères d’évaluation de l'instabilité interne des sols

Références Matériaux Critère géométrique Remarque

Istomina (1957) [24] Sols 𝐶𝑢 ≤ 10 : Stable
sableux 10 < 𝐶𝑢 ≤ 20 : Transitoire
𝐶𝑢 > 20 : Instable

Kézdi (1979) [42] Tous les (𝑑15𝑐 ⁄𝑑85𝑓 )𝑚𝑎𝑥 ≤ 4 : Instable

sols
Kenney et Lau (1985) Sols (𝐻 ⁄𝐹 )𝑚𝑖𝑛 ≥ 1 : Instable Sols
[31] granulaires moyennement
denses à denses
Burenkova (1993) [41] Sol
0.76 lg(ℎ′′ ) + 1 < ℎ′ < 1.68 lg(ℎ′′ ) + 1
granulaires
Instable
gradués
Mao (2005) [26] Sol
𝑝𝑓 ≥ 100⁄4(1 − 𝑛) : Instable
granulaires

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 19

Chapitre I Stabilité des ouvrages en terre : Généralités

Wan et Fell (2008) [43] Sols bien 30/log (d90/d60) < 80, ou
gradués (1 30/log (d90/d60) < 80 et 15/log (d20/d5) >
≤ Cc ≤ 3) 22
Stable
Li and Fannin (2008) Sols H et F sont
Pour F < 15, (H/F)min ≥ 1.0; Stable
[28] granulaires obtenus à
Pour F > 15, H ≥ 15; Stable
(H/F)min
h′: d90/d60; h″: d90/d15; n: porosité; pf : la fraction fine (Pour les sols à granulométrie discontinue, pf =
le passant au point de discontinuité (%); Pour les sols bien gradué, pf = le passant (%) correspondant
au diamètre de division df = 1.3√𝑑15 𝑑85 ). Cc=(D30)2/(D10xD60)

IV.3 Critères hydrauliques

Les critères géométriques précédemment mentionnés peuvent vérifier le potentiel
d’érosion interne d’un sol, étant donné que les critères hydrauliques permettent d’identifier
l’initiation de l’érosion interne. En d’autres termes, les critères géométriques reflètent les
propriétés intrinsèques des sols, tandis que les critères hydrauliques expriment les facteurs
externes pour le déclenchement de l’érosion interne [52]. La suffusion est déclenchée si les
forces hydrodynamiques induites par le flux d’infiltration sur les grains du sol dépassent un
seuil critique, c.-à-d. le gradient hydraulique critique 𝑖𝑐𝑟 .

Critère de Terzaghi
Terzaghi [29] a proposé un gradient hydraulique ascendant critique (𝑖𝑐𝑟 = 𝛾 ′ ⁄𝛾𝑤 ) égale
à 1 pour les sols sans cohésion soumis à une contrainte effective nulle. Cependant, ce critère a
été initialement introduit pour décrire le soulèvement hydraulique et ne s’applique pas à
l’érosion interne des sols. Les expériences menées par Skempton et Brogan [47] ont indiqué
que dans les graviers sableux instables, des fuites concentrées peuvent se produire à des
gradients critiques suggérés par Terzaghi pour un matériau granulaire homogène de même
porosité. Un facteur de réduction des contraintes 𝛼 a été introduit dans la définition de 𝑖𝑐𝑟 pour
tenir compte de l’influence des particules fines présentes dans la matrice des contraintes du sol :
𝑖𝑐𝑟 = (1 − 𝛼)(𝛾 ′ ⁄𝛾𝑤 ). Cependant, la définition de 𝛼 reste floue. Li [51] a suggéré que 𝛼 est lié
′
à 𝑑85 correspondant à la fraction fine et à l’ouverture 𝑂50 par l’expression suivante : 𝛼 =
′ ⁄
3.85(𝑑85 𝑂50 ) − 0.616.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 20

Chapitre I Stabilité des ouvrages en terre : Généralités

Moffat et Fannin [39] ont observé une tendance à l'augmentation de 𝑖𝑐𝑟 avec la
′
contrainte effective verticale moyenne 𝜎𝑣𝑚 . Li et Fannin [28] ont proposé une enveloppe hydro-
mécanique théorique pour caractériser le gradient hydraulique critique d'un sol soumis à une
contrainte de surcharge :

𝛼 0.5𝛾 ′
𝑖𝑐𝑟 = ̅′
(𝜎 𝑣𝑚 + ) (8)
1−𝛼 𝛾𝑤

Où 𝛼 est le coefficient de transformation des contraintes, 𝜎̅′ 𝑣𝑚 est la contrainte verticale

effective dans le sol.

Critère de Skempton et Brogan [47]

Skempton et Brogan [47] ont effectué des essais d’écoulement au laboratoire sur des
granulométries uniformes et discontinues constituées d’une grave sableuse. Le diamètre et la
hauteur des échantillons sont respectivement de 139 mm et 155 mm. Le but de leur étude est de
déterminer le gradient hydraulique à partir duquel commence la migration des particules et donc
l’instabilité interne dans les matériaux. Dans un perméamètre, un écoulement ascendant à
intensité croissante a été appliqué aux matériaux jusqu’à la rupture. A l’issue des essais, le
gradient hydraulique (ic) est représenté en fonction de l’indice de stabilité (H/F) défini par
Kenney et Lau [31, 48] (Figure I. 15). Ils remarquent une augmentation rapide de (ic) à partir
d'une faible valeur à une valeur élevée à partir de la valeur (H/F)min = 1. Ils notent aussi que
pour les matériaux instables, le gradient hydraulique critique varie entre 1/3 et 1/5 du gradient
théorique seuil égal à 1 proposé par Terzaghi. Cependant, ils font aussi remarquer que la relation
entre le gradient hydraulique (ic) et l’indice de stabilité (H/F) est très influencée par la densité
de l’échantillon et par la proportion de sable dans du gravier. Cette relation nécessite donc une
étude plus poussée.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 21

Chapitre I Stabilité des ouvrages en terre : Généralités

Figure I. 15 Relation entre le gradient critique ic et l’indice de stabilité (H/F)min [47]

Critère de Wan et Fell [27]

Wan et Fell [27] ont étudié l’impact des écoulements sur l’instabilité interne des sols.
Les mêmes échantillons compactés et constitués à partir des mélanges limon-sable-gravier et
argile-limon-sable-gravier sont soumis à un écoulement ascendant. Au cours de ces essais, le
gradient appliqué est augmenté par palier jusqu’ à l’apparition de l’instabilité interne ou que le
gradient le plus élevé possible qui est de 20 soit atteint. Selon ces auteurs, l’instabilité qui
correspond à l’érosion des fines est initiée lorsque l’écoulement aval commence à devenir
tourbillon. Ils concluent que dans les sols instables le gradient hydraulique critique est inférieur
à 0.8 voire à 0.5, et inférieur à 1 dans les sols stables. Dans ce dernier, la suffusion se manifeste
par une érosion régressive ou renard. Contrairement aux auteurs comme Adel et al. [53] ;
Skempton et Brogan [47], ne définissent aucune relation entre le gradient hydraulique critique
et le rapport de stabilité minimum (H/F)min; le coefficient d’uniformité Cu ou le taux de fines
contenu dans l’échantillon. Des tentatives pour définir des telles relations ont montré une grande
dispersion des données (Figure I. 16 et Figure I. 17).

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 22

Chapitre I Stabilité des ouvrages en terre : Généralités

Figure I. 16 Gradient hydraulique moyenne à laquelle débute l'érosion de particules fines dans
l'échantillon (istart) en fonction du rapport de stabilité minimum (H/F) min [27]

Figure I. 17 Gradient hydraulique moyenne à laquelle débute l'érosion de particules fines dans
l'échantillon (istart) en fonction du coefficient d’uniformité Cu [27]

Wan et Fell [27] ont montré aussi que la présence de fines n’influence pas l’instabilité
et pourrait même renforcer la résistance à la suffusion. Cette résistance est d’autant plus
importante que ces fines sont de type argileux (kaolinite). Dans les granulométries discontinues
l’érosion débute généralement à des faibles gradients (ic < 0.2) [27]. Cependant le risque que
l’instabilité se produise à des faibles valeurs du gradient hydraulique (ic < 0.2) lorsque la

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 23

 Chapitre I Stabilité des ouvrages en terre : Généralités

porosité augmente a été relevé (Figure I. 18). Par exemple, une augmentation de la densité du
sol de 90% à 95% ferait passer le gradient hydraulique critique d'environ 0,2 à 3,2.

Figure I. 18 Gradient hydraulique moyenne à laquelle débute l'érosion de particules fines dans
l'échantillon (istart = ic) en fonction de la porosité [27]

Résumé sur les critères d’instabilité

Nous avons présenté dans la section précédente une liste de critères utilisés pour
l’évaluation de la susceptibilité à l’érosion interne des matériaux. Ces critères sont divisés en
deux types. Les critères géométriques basés sur la forme de la courbe granulométrique du
matériau, et les critères hydrauliques basés sur le gradient hydraulique pour l’initiation de
l’érosion. Cependant, la liste présentée ici n’est pas exhaustive et il existe d’autres critères peu
utilisés. Les critères présentés sont essentiellement ceux qui ont été utilisés dans ce travail afin
d’évaluer la susceptibilité des matériaux étudiés.

V. Stabilité des pentes

Dans la plupart des cas, l'objectif principal de l'analyse de la stabilité des pentes est de
contribuer à la conception sûre et économique des excavations, des remblais, des barrages en
terre, des digues, des décharges et des terrils. Les évaluations de la stabilité des pentes visent à
identifier les paramètres géologiques, matériels, environnementaux et économiques critiques
qui affecteront le projet, ainsi qu'à comprendre la nature, l'ampleur et la fréquence des
problèmes potentiels liés aux pentes. Pour les pentes en général et l'analyse de la stabilité des

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 24

Chapitre I Stabilité des ouvrages en terre : Généralités

pentes en particulier, l'expérience géologique et géotechnique acquise dans une région donnée
est précieuse.

Les objectifs principaux des analyses de stabilité des pentes sont les suivants

(l) Comprendre le développement et la forme des pentes naturelles et les processus

responsables des différentes évolutions naturelles.

(2) Évaluer la stabilité des pentes à court terme (souvent pendant la construction) et des
conditions de long terme.

(3) Évaluer la possibilité de glissements de terrain impliquant des pentes.

(4) Analyser les glissements de terrain et comprendre les mécanismes de défaillance et

l'influence de facteurs environnementaux.

(5) Permettre le réaménagement des pentes endommagées ainsi que la planification et

la conception de mesures préventives et correctives, le cas échéant.

(6) Étudier l'effet des charges sismiques sur les pentes et les remblais.

L'analyse des pentes tient compte de divers facteurs liés à la topographie, à la géologie
et aux propriétés des matériaux, en fonction, souvent, du fait que la pente a été façonnée
naturellement ou artificiellement.

V.1 Equilibre plastique limite

L'objectif principal de la plupart des analyses de stabilité est de déterminer le facteur de
sécurité de la pente en se basant sur le concept d'équilibre plastique limite. Tout d'abord, on
suppose une surface de rupture. On dit qu'il existe un état d'équilibre limite lorsque la contrainte
de cisaillement le long de la surface de rupture est exprimée comme suit

𝑠
𝜏=
𝐹

Où 𝜏 = contrainte⁡de⁡cisaillement⁡, 𝑠 = 𝑟é𝑠𝑖𝑠𝑡𝑎𝑛𝑐𝑒⁡𝑎𝑢⁡𝑐𝑖𝑠𝑎𝑖𝑙𝑙𝑒𝑚𝑒𝑛𝑡 et 𝐹=
𝑓𝑎𝑐𝑡𝑒𝑢𝑟⁡𝑑𝑒⁡𝑠é𝑐𝑢𝑟𝑖𝑡é. Selon le critère de Mohr-Coulomb, la résistance au cisaillement peut être
exprimée par

𝑠 = 𝐶 + 𝜎𝑛 tan 𝜑

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 25

Chapitre I Stabilité des ouvrages en terre : Généralités

Avec 𝐶= Cohésion, 𝜎𝑛 = Contrainte normale et 𝜑= Angle de frottement interne. 𝐶 et 𝜑

sont connus par les paramètres de résistance au cisaillement des sols.

V.2 Méthode d’analyse de la stabilité des pentes

La méthode utilisée dans cette étude est la méthode d’équilibre limite. Une autre
méthode d'analyse de la stabilité est la méthode des éléments finis, qui se base sur la mécanique
des solides en considérant non seulement les équations d'équilibre mais aussi celles de
compatibilité. L'avantage de la méthode des éléments finis est sa capacité à déterminer non
seulement le facteur de sécurité mais aussi les déplacements. Cependant, les déplacements
dépendent fortement du module d'élasticité, qui est difficile à déterminer. La méthode des
éléments finis ne sera pas abordée dans cette étude.

En raison du grand nombre de méthodes d'équilibre limite, il n'est pas possible de passer
en revue chacune d'entre elles. Par conséquent, seules les méthodes utilisées dans cette étude
seront abordées. Hopkins et al. [54] et Duncan [55] ont présenté une revue de plusieurs
méthodes. Les méthodes sont basées sur la méthode des tranches et peuvent être divisées
globalement en quatre catégories, en fonction du nombre d'équations d'équilibre à satisfaire.

Méthodes basées sur l'équilibre du moment global

Ces méthodes ne sont applicables qu'aux surfaces de rupture circulaires et ne prennent
en compte que l'équilibre global des moments. L'équilibre global des forces n'est ni considéré
du tout ni satisfait dans les deux directions. On trouve dans cette catégorie la méthode de
Fellenius [56], la méthode normale [57] et la méthode simplifiée de Bishop [58].

Méthodes basées sur l'équilibre global des moments et des forces

La méthode originale de Spencer [59] fait partie de cette catégorie qui suppose que
toutes les forces entre tranches sont parallèles et s'inclinent à un angle δ avec l'horizontale, où
δ est une inconnue à déterminer. Elle considère l'équilibre global du moment, l'équilibre global
des forces dans la direction δ, et l'équilibre des forces de chaque tranche dans la direction
perpendiculaire à δ. Comme l'équilibre global des forces dans deux directions perpendiculaires
est entièrement satisfait, la méthode peut également être utilisée pour les surfaces de rupture
non circulaires, où un centre de moment doit être choisi arbitrairement. La méthode originale
de Spencer a été affinée par la suite et est appelée ci-après la méthode de Spencer [60], dans
laquelle l'équilibre des moments est également satisfait pour chaque tranche. La méthode

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 26

Chapitre I Stabilité des ouvrages en terre : Généralités

originale de Spencer a l'avantage que le facteur de sécurité converge toujours, alors que la
méthode de Spencer peut parfois rencontrer des problèmes de convergence.

Méthodes basées sur l'équilibre des forces de chaque tranche

Ces méthodes ne prennent en compte que l'équilibre des forces dans chaque tranche.
Une fois l'équilibre des forces satisfait dans chaque tranche, l'équilibre global des forces sera
automatiquement satisfait. Bien que l'équilibre des moments ne soit pas pris en compte
explicitement, ces méthodes peuvent donner des résultats précis si l'inclinaison des forces entre
les tranches est supposée de manière à ce que l'équilibre des moments soit satisfait
implicitement. Des hypothèses arbitraires sur l'inclinaison des forces entre les tranches peuvent
avoir une grande influence sur le facteur de sécurité. Selon l'inclinaison des forces entre les
tranches, une gamme de facteurs de sécurité peut être obtenue dans de nombreux problèmes.
Les méthodes d'équilibre des forces doivent être utilisées avec précaution, et l'utilisateur doit
être bien conscient des hypothèses particulières de forces entre les tranches utilisées. Sont
incluses dans cette catégorie les méthodes suggérées par Janbu [61, 62], Lowe et Karafiath [63],
et l'U.S. Army Corps of Engineers [64].

Méthodes basées sur l'équilibre des moments et des forces de chaque

tranche
On trouve dans cette catégorie la méthode de Spencer [59], la méthode de Janbu [61,
62], et la méthode de Morgenstern et Price [65]. Toutes ces méthodes considèrent l'équilibre
des moments et des forces dans chaque tranche. Si l'équilibre des moments et des forces est
satisfait dans chaque tranche, l'équilibre global sera automatiquement satisfait. Le concept de
base de ces méthodes est le même ; la différence réside dans l'hypothèse des forces entre les
tranches. Si l'équilibre des moments et des forces est respecté, l'hypothèse des forces entre les
tranches ne devrait avoir qu'un faible effet sur le facteur de sécurité obtenu. Toutes ces méthodes
peuvent être appliquées à des surfaces de rupture circulaires et non circulaires.

V.3 Méthode des tranches

La méthode des tranches est un outil très puissant qui peut être utilisé pour analyser des
pentes de n'importe quelle configuration composée de différents sols et conditions des eaux
souterraines.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 27

Chapitre I Stabilité des ouvrages en terre : Généralités

Figure I. 19 Forces agissantes pour la méthode des tranches appliquée à une surface de
glissement composite

La figure I.19 montre les forces qui doivent être définies pour un problème général de
stabilité de pente. Les variables associées à chaque tranche sont définies comme suit :

W : poids total d’une tranche de largeur b et de hauteur h

P : force normale totale sur la base de la tranche sur une longueur l

Sm : force de cisaillement mobilisée sur la base de la tranche. Il s'agit d'un pourcentage

de la résistance au cisaillement telle que définie par l'équation de Mohr-Coulomb. Autrement
dit, 𝑆𝑚 = 𝑙{𝑐 ′ + [𝑃⁄𝑙 − 𝑢] tan 𝜑 ′ }/𝐹 où 𝑐 ′ est la cohésion effective, 𝜑 ′ angle de frottement
interne effectif, F est le facteur de sécurité et u est la pression interstitielle.

R : le rayon ou le bras de levier associé à l'effort de cisaillement mobilisé Sm

f : décalage perpendiculaire de la force normale par rapport au centre de rotation

x : distance horizontale de la tranche au centre de rotation

α : angle entre la tangente au centre de la base de chaque tranche et l'horizontale.

E : forces horizontales inter tranches

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 28

Chapitre I Stabilité des ouvrages en terre : Généralités

L : indice désignant le côté gauche

R : indice désignant le côté droit

X : forces verticales inter tranches

k : coefficient sismique pour tenir compte d'une force horizontale dynamique

e : distance verticale du centre de gravité de chaque tranche au centre de rotation

Une charge uniforme sur la surface peut être prise en compte comme une couche de sol
de poids unitaire et de densité appropriés. Les variables suivantes sont nécessaires pour définir
une charge linéaire :

L : charge linéaire (force par unité de largeur)

ω : angle de la charge linéaire par rapport à l'horizontale

d : distance perpendiculaire de la ligne de charge au centre de rotation

L'effet de la submersion partielle du talus ou des fissures de tension dans l'eau nécessite
la définition de variables supplémentaires :

A : forces résultantes de l'eau

a : distance perpendiculaire de la résultante de la force de l'eau au centre de rotation

Détermination du facteur de sécurité

Les éléments de la statique qui peuvent être utilisés pour calculer le facteur de sécurité
sont la somme des forces dans deux directions et la somme des moments. Ceux-ci, ainsi que les
critères de rupture, sont insuffisants pour rendre le problème déterminant. Il faut en savoir plus
sur la distribution des forces normales ou la distribution des forces inter tranches. Des éléments
supplémentaires de la physique ou une hypothèse doivent être invoqués pour rendre le problème
déterminant. Toutes les méthodes considérées dans cette étude utilisent cette dernière
procédure, chaque hypothèse donnant lieu à une méthode d'analyse différente. À des fins de
comparaison, chaque équation est obtenue en utilisant de manière cohérente les équations de la
statique.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 29

Chapitre I Stabilité des ouvrages en terre : Généralités

Méthode ordinaire de Fellenius

La méthode ordinaire est considérée comme la plus simple des méthodes de tranches
car c'est la seule procédure qui aboutit à une équation linéaire du facteur de sécurité. On affirme
généralement que les forces entre les tranches peuvent être négligées parce qu'elles sont
parallèles à la base de chaque tranche [56]. Cependant, le principe de Newton "action égale
réaction" n'est pas satisfait entre les tranches. Le changement de direction sans discernement de
la force inter-tranches résultante d'une tranche à l'autre entraîne des erreurs de facteur de
sécurité qui peuvent atteindre 60 % [66].

La force normale sur la base de chaque tranche est dérivée soit de la somme des forces
perpendiculaires à la base, soit de la somme des forces dans les directions verticale et
horizontale.

∑ 𝐹𝑉 = 0 ⇒ 𝑊 − 𝑃 cos 𝛼 − 𝑆𝑚 sin 𝛼 = 0 (9)

∑ 𝐹𝐻 = 0 ⇒ 𝑆𝑚 cos 𝛼 − 𝑃 sin 𝛼 − 𝑘𝑊 = 0 (10)

En substituant (10) dans (9) et en déterminant la force normale, on obtient

𝑃 = 𝑊 cos 𝛼 − 𝑘𝑊 sin 𝛼 (11)

Le facteur de sécurité est obtenu par la somme des moments autour d'un point commun
(c'est-à-dire un centre de rotation fictif ou réel pour l'ensemble de la masse).

∑ 𝑀𝑂 = 0 ⇒⁡ ∑ 𝑊𝑥 − ∑ 𝑆𝑚 𝑅 − ∑ 𝑃𝑓 + ∑ 𝑘𝑊𝑒 ± 𝐴𝑎 + 𝐿𝑑 = 0 (12)

En introduisant les critères de rupture et la force normale de (11) et en résolvant pour le facteur
de sécurité, on obtient

∑{𝑐 ′ 𝑙𝑅 + (𝑃 − 𝑢𝑙)𝑅 tan 𝜑 ′ }

𝐹= (13)
∑ 𝑊𝑥 − ∑ 𝑃𝑓 + ∑ 𝑘𝑊𝑒 ± 𝐴𝑎 + 𝐿𝑑

Méthode simplifiée de Bishop

La méthode simplifiée de Bishop néglige les forces de cisaillement entre les tranches et
suppose donc qu'une force normale ou horizontale définit adéquatement les forces entre les
tranches [144]. La force normale à la base de chaque tranche est obtenue en additionnant les
forces dans une direction verticale. En substituant les critères de rupture et en résolvant la force
normale, on obtient :

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 30

Chapitre I Stabilité des ouvrages en terre : Généralités

𝑐 ′ 𝑙 sin 𝛼 𝑢𝑙 tan 𝜑 ′ sin 𝛼

𝑃 = [𝑊 − + ]⁄𝑚𝛼 (14)
𝐹 𝐹
𝑚𝜶 = cos 𝛼 + (sin 𝛼 tan 𝜑 ′ )⁄𝐹

Le facteur de sécurité est obtenu de la somme des moments autour d'un point commun.
Cette équation est la même que (12) puisque les forces inter tranches s'annulent. Par conséquent,
l'équation du facteur de sécurité est la même que pour la méthode ordinaire (équation 13).
Cependant, la définition de la force normale est différente.

Méthode de Spencer
La méthode de Spencer suppose qu'il existe une relation constante entre les efforts
normaux et de cisaillement entre les tranches [67].

𝑋𝐿 𝑋𝑅
tan 𝜃 = = (15)
𝐸𝐿 𝐸𝑅
Où θ = angle de la force résultante entre les tranches par rapport à l'horizontale.

La force normale peut être déduite de l’équilibre des forces horizontales et verticales

𝑐 ′ 𝑙 sin 𝛼 𝑢𝑙 tan 𝜑 ′ sin 𝛼

𝑃 = [𝑊 − (𝐸𝑅 − 𝐸𝐿 ) tan 𝜃 − + ]⁄𝑚𝛼 (16)
𝐹 𝐹
Les forces entre les tranches (𝐸𝑅 − 𝐸𝐿 ) doivent s'annuler et l'équation du facteur de
sécurité par rapport à l'équilibre des forces se réduit à

∑{𝑐 ′ 𝑙 cos 𝛼 + (𝑃 − 𝑢𝑙) tan 𝜑 ′ cos 𝛼}

𝐹= (17)
∑ 𝑃 + sin 𝛼 + ∑ 𝑘𝑊 ± 𝐴 − 𝐿 cos 𝜔

Méthode simplifiée de Janbu

La méthode simplifiée de Janbu utilise un facteur de correction 𝑓0 pour tenir compte de
l'effet des forces de cisaillement inter tranches. Le facteur de correction est lié à la cohésion, à
l'angle de frottement interne et à la forme de la surface de rupture [61]. La force normale est
obtenue de la somme des forces verticales, les forces de cisaillement inter tranches étant
ignorées.

En utilisant l’équilibre des forces horizontales, on obtient l’expression du facteur de

sécurité tel que

∑{𝑐 ′ 𝑙 cos 𝛼 + (𝑃 − 𝑢𝑙) tan 𝜑 ′ cos 𝛼}

𝐹0 =
∑ 𝑃 + sin 𝛼 + ∑ 𝑘𝑊 ± 𝐴 − 𝐿 cos 𝜔

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 31

Chapitre I Stabilité des ouvrages en terre : Généralités

𝐹0 est utilisé pour désigner le facteur de sécurité non corrigé.

Le facteur de sécurité corrigé est donné par

𝐹 = 𝑓0 𝐹0

𝑓0 est un coefficient de correction pour prendre en compte l’effet des forces de

cisaillement inter tranches.

Méthode de Morgenstern Price

La méthode de Morgenstern Price suppose une fonction mathématique arbitraire pour
décrire la direction des forces inter tranches tel que

𝜆𝑓(𝑥) = 𝑋⁄𝐸

Où 𝜆 = une constante à évaluer dans la résolution du facteur de sécurité, et 𝑓(𝑥) est une
fonction de variation par rapport à x. La figure I.20 montre des fonctions typiques (c'est-à-dire
f(x)). Pour une fonction constante, la méthode de Morgenstern-Price est identique à la méthode
de Spencer. La figure I.21 montre comment la fonction demi-sinus et x sont utilisés pour
désigner la direction des forces inter-tranches.

Morgenstern et Price [65] ont basé leur solution sur la somme des forces tangentielles
et normales à chaque tranche. Les équations d'équilibre des forces ont été combinées, puis la
technique numérique de Newton-Raphson a été utilisée pour résoudre les équations de moment
et de force pour le facteur de sécurité et λ.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 32

Chapitre I Stabilité des ouvrages en terre : Généralités

Figure I. 20 Variation fonctionnelle de la direction de la force latérale par rapport à la direction

x [68]

Figure I. 21 Désignation de la force latérale pour la méthode Morgenstern-Price

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 33

 Chapitre II

Revue bibliographique
 Chapitre II Revue bibliographique

I. Introduction

Ce chapitre est consacré à la présentation d’un état de l’art des approches expérimentales
et numérique existantes dans la littérature. L’accent est mis principalement sur l’influence de
l’érosion interne sur le comportement mécanique des sols ainsi que sur la stabilité globale d’un
ouvrage en terre. Des exemples de cas de ruptures des ouvrages en terre où l’érosion interne a
été identifiée comme le facteur déclenchant principale ont été présentés dans ce chapitre.

II. Influence de l’érosion interne sur le comportement mécanique

des sols

Sterpi [69] a effectué plusieurs tests de suffusion sur des échantillons de sable limoneux
afin d'étudier la quantité de particules érodées sous des gradients hydrauliques contrôlés. Une
série de tests de compression triaxiale a également été réalisée pour évaluer l'effet de l'érosion
sur le comportement mécanique des sols testés. Ces essais triaxiaux ont été réalisés sur des
échantillons reconstitués avec différentes teneurs en particules fines afin de reproduire les
différentes étapes du développement de l'érosion. Les essais de compression ont montré que la
rigidité et la résistance au cisaillement des échantillons testés augmentent avec la diminution
de la teneur en particules fines. En d'autres termes, il semblerait que la rigidité et la résistance
au cisaillement ne soient pas dégradées par l'érosion des particules fines. Par ailleurs, Chang et
Zhang [6] ont effectué des tests de compression sur des échantillons ayant été soumis à l'érosion
par suffusion et ont montré que la résistance au cisaillement et la rigidité étaient, dans leur cas,
fortement réduites par le processus d'érosion. Par conséquent, l'influence de la suffusion sur les
propriétés mécaniques des sols est encore aujourd'hui une question ouverte et peut dépendre
fortement de la granulométrie initiale du sol [70].

Chen et al. (2016) ont effectué des essais de laboratoire visant à évaluer le comportement
mécanique des sols à granulométrie discontinue soumis à l’érosion interne. Les tests ont été
réalisés sur un appareil triaxial piloté par ordinateur pour contrôler l'état de contrainte du sol,
introduire des infiltrations dans l’échantillon et effectuer des tests de compression triaxiale. La
Figure II. 1 montre le dispositif expérimental utilisé et une description plus détaillée est donnée
par Chang et Zhang (2011) [71]. La particularité de cette étude est que les auteurs contrôlent le

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 34

Chapitre II Revue bibliographique

taux d’érosion en ajoutant une quantité prédéfinie de sel dans l’échantillon ce qui permets ainsi
d’accélérer le processus. Les déformations radiales et axiales au cours de l'érosion ont été
mesurées à l'aide d'une méthode photographique. Par la suite, des essais de compression
triaxiale drainée ont été réalisés pour étudier le comportement contrainte-déformation des sols
qui avaient perdu différentes quantités de particules fines. Les résultats des essais ont montré
que l'angle de frottement maximal et l'angle de frottement critique du sol ont diminué avec la
perte de particules fines. Après une perte significative de particules fines, le comportement
contrainte-déformation est passé d'un comportement initialement dilatateur à un comportement
plus contraignant.

Figure II. 1 schéma du dispositif expérimental utilisé par Chen et al. (2016) [71].

Ke et Takahashi (2012) [72] ont mis l'accent dans leur étude sur les changements de la
résistance des sols dus à l'érosion interne. Une série de tests d'infiltration unidimensionnels
ascendante à charge constante est effectuée pour provoquer l'érosion interne d'un échantillon
de sol en contrôlant trois paramètres, à savoir (a) la teneur en fines, (b) la densité relative du
sol, et (c) le gradient hydraulique maximal imposé à l'échantillon. Les conséquences
mécaniques de l'érosion interne sont examinées par des tests CPT. Les résultats des essais ont
montré que la perte des particules fines causée par l'érosion interne provoque une augmentation
significative du taux de vide et de la conductivité hydraulique, ce qui entraîne une diminution
de la résistance du sol par rapport à sa valeur initiale. L'influence de l'érosion interne sur les
essais CPT est montrée par les profils de résistance de pointe du cône obtenus avant et après

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 35

 Chapitre II Revue bibliographique

l'érosion interne. La Figure II. 2 montre les résultats CPT avant et après érosion d’un échantillon
étudié.

Figure II. 2 Résistance du cône avant et après l'érosion interne [72]

Benamar et al. (2019) [73] ont effectué une étude expérimentale de l’érosion interne sur
des échantillons prélevés dans les digues du Rhin en France. Ces digues sont composées de sols
à grains grossiers qui, d'après les critères empiriques disponibles, ont un fort potentiel de
suffusion. Le but de cette étude était d’évaluer le potentiel de suffusion de tels sols et à
caractériser leur stabilité et la variation des paramètres physiques du sol au cours du processus
de suffusion. Les tests ont été réalisés sur un dispositif expérimental permettant l’application
d’un flux ascendant à charge hydraulique contrôlée. Les auteurs ont constaté que les particules
érodées ont une taille inférieure à 500 µm. Les auteurs ont aussi déduit que les méthodes
habituelles d’évaluation de la stabilité interne basées sur des critères géométriques surestiment
la susceptibilité à la suffusion des sols du Rhin et, par conséquent, ils considèrent que, pour de
tels sols, il est préférable d'effectuer des tests de laboratoire pour évaluer le processus de
suffusion.

III. Influence de l’érosion interne sur la stabilité

Jie Yang [74] a rapporté dans sa thèse des investigations expérimentales et numériques
multi échelles. A l’échelle de laboratoire, il a formulé un modèle numérique pour les tests
d’érosion sous chargement mécaniques triaxiales afin d’analyser le comportement couplé
hydromécanique des sols soumis à l’érosion interne. Le modèle est formulé sur la base de la
théorie des milieux poreux en considérant un milieu composé de quatre phases. La variation

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 36

Chapitre II Revue bibliographique

spatiale de la porosité et de la teneur en fines a été pris en compte dans l’évaluation de

l’autofiltration lors de la suffusion unidimensionnelle (1D) et bidimensionnelle (2D). A
l’échelle de la géo-structure un modèle hydromécanique a été développé et implémenté dans le
code de calcul par éléments finis (ABAQUS). Des simulations numériques ont été réalisée sur
une digue à l’échelle réelle afin de quantifier l’effet de l’érosion interne sur sa stabilité. La digue
étudiée est celle de Val d’Orléans (France), elle a été construite sur une formation d'alluvions
sablo-limoneuses sur des calcaires karstiques (Figure II. 3). Un aquifère confiné sous pression
est situé sous la couche de calcaire. Elle peut s'élever à l'intérieur des grottes karstiques et même
atteindre la couche d'alluvions. L'érosion interne due au flux d'infiltration peut avoir lieu lorsque
le niveau d'eau est différent dans l'alluvion et dans l'aquifère confiné.

Figure II. 3 Schéma de la digue sur fondation d'Orléans (France) [74].

Horikoshi et Takahashi [75] ont réalisé des essais de suffusion sur des modèles
physiques réduits de remblais sous un gradient hydraulique constant (Figure II. 4). Les auteurs
ont tenté de comprendre l’évolution du phénomène de suffusion en visualisant le changement
temporel de la distribution spatiale des particules fine dans le remblai. Les résultats des tests
ont montré une forte diminution des particules fines dans le pied du remblai (zone B) et dans la
zone près de la ligne phréatique (Zone C) (Figure II. 4b). Après l’érosion du pied de la pente, les
auteurs ont remarqué qu’une érosion régressive se développe le long de la ligne phréatique. De
plus sous la ligne phréatique, les fines érodées se déplacent non seulement latéralement par
l’écoulement, mais aussi verticalement en raison de la force gravitationnelle, et se déposent
dans la fondation. Ce dépôt entraîne l'expansion de la région riche en fines dans la fondation et
provoque une diminution de la perméabilité de l'ensemble du remblai.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 37

Chapitre II Revue bibliographique

(a) (b)

Figure II. 4 (a) Schéma du modèle physique de Horikoshi et Takahashi (b) Répartition de la
variation spatiale du pourcentage des particules fines [75]

Dans le but d'étudier l'érosion régressive, van Beek et al [76] ont réalisé une série
d'expériences à petite, moyenne et grande échelle (Figure II. 5). Les phases observées dans les
expériences étaient très similaires : infiltration, érosion régressive, élargissement du canal,
rupture et formation de brèche dans la digue. Au départ, un transport de sable très limité et
presque négligeable a été observé. Une augmentation de la charge hydraulique a entraîné une
érosion régressive, marquée par l'apparition de dépôts de sable. Les auteurs ont suggéré que la
transition entre l'érosion régressive et l'élargissement du conduit ne pouvait pas être observée
par la quantité de sable transporté ou par le débit. La défaillance s'est produite par la suite dans
un court laps de temps.

(a)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 38

Chapitre II Revue bibliographique

(b)

(c)

Figure II. 5 Installation expérimentale de van Beek et al (77) : (a) petite échelle ; (b) moyenne
échelle ; (c) grande échelle

Hu et al [77] ont effectué des tests dans un canal pour étudier le déclenchement d'un
glissement de terrain fluidisé. Le matériau testé a été prélevé sur un dépôt de glissement de
terrain cosismique dans la zone du tremblement de terre de Wenchuan en 2008. La pression
interstitielle, le déplacement interne et les vibrations induites dans le sol ont été surveillés
pendant les expériences. Le sol a été progressivement humidifié soit par un débit uniforme au
sommet de la pente pour simuler l'infiltration des eaux de ruissellement, soit par une élévation
progressive du niveau de l'eau à la limite supérieure de la pente pour simuler une élévation de
la nappe phréatique (Figure II. 6). La microsismicité a été surveillée. Il semble que l'érosion
interne de la fraction fine du sol, provoquée par les forces d'infiltration, ait joué un rôle
important dans l'effondrement de la pente. Le réarrangement possible des particules ont été
enregistrés par les accéléromètres sous la forme d'un signal de vibration haute fréquence
croissant avant que l'on puisse observer une augmentation significative de la pression
interstitielle et un déplacement important de la pente.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 39

Chapitre II Revue bibliographique

(a)

(b)

Figure II. 6 Mise en place des deux canaux par Hu et al. [77] : (a) Infiltration d'eau par le
sommet de la pente ; (b) Augmentation du niveau d'eau à la frontière hydraulique supérieure

Les travaux de Fox et al. [78] ont porté sur l’étude au laboratoire du processus de rupture
des berges causé par l’érosion interne. Des expériences de laboratoire ont été réalisées sur des
sols non cohésifs à l'aide d'un lysimètre bidimensionnel. Une série de tests a été réalisée en
variant la pente des berges (90, 60, 45, 36 et 26°). Les résultats de cette étude ont mis en
évidence trois phases du processus de rupture causé par l’érosion : (i) l'érosion par infiltration ;
(ii) la formation de fissures de tension ; et (iii) l'effondrement de la berge (Figure II. 8). Afin de
prédire le taux d’érosion, les auteurs ont apporté une légère modification sur la loi d’érosion
déjà existante en négligeant la contrainte de cisaillement critique (𝑞𝑠∗ = 𝑘𝑑𝑠 (𝜏 ∗ − 𝜏𝑐𝑠
∗ 𝑎
) , avec

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 40

Chapitre II Revue bibliographique

𝑞𝑠∗ est le flux de sédiment, 𝑘𝑑𝑠 est le coefficient d’érosion, 𝜏 ∗ est la contrainte de cisaillement
∗
appliquée, 𝜏𝑐𝑠 est la contrainte de cisaillement critique et a est un paramètre déterminé
expérimentalement). En tenant compte de l’érosion du pied de la pente, les auteurs ont évalué
l’évolution du facteur de sécurité en fonction du temps à l’aide du logiciel SLOPE/W (Figure
II. 7).

Figure II. 7 simulation de l’évolution du facteur de sécurité pour une pente de 60° [78]

Figure II. 8 Série chronologique typique d'effondrement dû à l'érosion interne [78]

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 41

 Chapitre II Revue bibliographique

IV. Les approches de modélisation de l’érosion interne

IV.1 Approches numériques de l’érosion interne

L'une des difficultés de la modélisation du processus d'érosion est qu'il se produit très
lentement sur le terrain, sur une période de plusieurs décennies. Les études expérimentales en
laboratoire peuvent utiliser des gradients hydrauliques artificiels élevés afin d'accélérer
l'apparition du phénomène, mais les microstructures qui se développent ainsi dans les sols
seront probablement différentes de celles qui se développent naturellement à un rythme plus
lent [79]. Par conséquent, l’étude numérique s’avère une alternative intéressante pour la
compréhension du comportement hydromécanique des sols érodables. Des modèles
mathématiques et numériques ont été développés pour décrire l'initiation et la progression du
processus de l'érosion interne, soit dans le cadre de la mécanique du continuum en utilisant les
méthodes des éléments finis ou des différences finies [80, 81, 90, 82–89], ou dans une approche
associant la méthode des éléments discrets pour décrire la dynamique de la phase solide
granulaire avec la dynamique des fluides pour résoudre l'écoulement des fluides interstitiels
[91–99].

Approche continue
Les modèles éléments finis sont basés principalement sur la théorie des milieux poreux
où l’érosion du sol est représentée par un échange de masse entre la masse solide et le fluide.

Vardoulakis et al. [85] ont développé un modèle mathématique pour le problème de

production de sable. Le milieu poreux saturé est supposé être composé du solide, du fluide et
des particules solides fluidisées. L’érosion de trou et l'érosion de surface ont été étudiées sur la
base des équations du bilan de masse et du transport des particules, ainsi que sur l'écoulement
du fluide dans un milieu poreux continu contrôlé par la loi de Darcy. Les équations
gouvernantes dérivées ont été résolues numériquement en utilisant la méthode des différences
finies en condition 1D pour calculer l'évolution de la porosité, la pression interstitielle et la
concentration des particules fines fluidisées pendant l'érosion.

Stavropoulou et al. [84] ont étendu l’étude de Vardoulakis et al. [85] sur le problème de
production de sable au couplage de l’érosion et de l’endommagement mécanique d’une roche
autour d’un puits de forage. Cette étude suppose que l’élasticité (Module de Young) et la
résistance (Cohésion) de la roche dépendent de la porosité de telle sorte que le matériau

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 42

Chapitre II Revue bibliographique

s'affaiblit à mesure que la porosité augmente. Le modèle mathématique, composé d'équations

d'érosion, d'équations d'écoulement du mélange et d'équations d'équilibre des contraintes, est
résolu numériquement par la méthode des éléments finis de Galerkin. Les résultats numériques
montrent que l'érosion, qui entraîne la production de sable, est élevée à proximité de la surface
libre.

Des modèles mathématiques similaires ont été développés par Papamichos et al. [82],
Cividini et Gioda [80], Steeb et al (94), Uzuoka et al [100], Zhang et al [87]. Cependant, aucun
modèle ne peut prendre en compte avec précision la modification de la résistance au
cisaillement du sol après une érosion interne. La partie mécanique consistant en un modèle
élastique avec endommagement est trop simple pour cela.

Zhang et al. [87] se sont intéressé à l’étude numérique de l’effet du transport de

particules fines provoqué par l'eau de pluie sur la stabilité d’une pente naturelle. Les auteurs
ont proposés un modèle couplé d’écoulement-érosion dans les sols insaturés. La loi d’érosion
proposée ici est une version modifiée de celle proposée par Cividini et al. [101]. Les équations
couplées ont été résolu par la méthode des éléments finis en utilisant le programme Comsol
Multiphysics [102]. La Figure II. 9 montre le modèle éléments finis illustratif de la pente étudiée
avec un angle de pente β de 35 °. La pente est composée de sols résiduels. L'épaisseur de la
couche de sol est de 6 m. Le modèle des éléments finis est composé de 1100 éléments
quadrilatéraux. Les résultats obtenus montrent que la pente atteint un état instable plus
rapidement si l'on considère l'érosion interne. Après 68 min d'infiltration, la surface de
glissement critique est d'environ 3,0 m, et le facteur de sécurité de la pente est inférieur à
1. Cependant, le facteur de sécurité de la pente sans prise en compte de l’érosion interne est
toujours supérieur à 1. Ce résultat implique que l'érosion interne peut entraîner un glissement
de terrain plus brutal.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 43

Chapitre II Revue bibliographique

Figure II. 9 Modèle éléments finis de l’exemple illustratif [87]

Cividini et al. (2009) [101] ont effectué une étude expérimentale et numérique sur l’effet
de l’érosion interne sur les tassements des sols pulvérulents. L’étude expérimentale avait pour
objectif d’évaluer la variation de la densité du matériau en fonction de la masse érodée. Afin
d’estimer les tassements, une analyse numérique des contraintes dont laquelle les déformations
sont liées aux résultats expérimentaux de la densité du matériau érodé est effectuée par la
méthode des éléments finis. En conclusion, les auteurs ont donné un intervalle de déformations
volumiques liées à l’érosion interne (équation 18).

∆𝜌𝑓 . 𝛾0 /𝛾𝑔 ≤ 𝜀𝑣𝑜𝑙 ≤ ∆𝜌𝑓 (18)

Avec ∆𝜌𝑓 est la variation de la densité des particules fine, 𝛾0 est le poids volumique sec
initial, 𝛾𝑔 est le poids volumique des grains et 𝜀𝑣𝑜𝑙 est la déformation volumique totale.

Cividini et Gioda [80] ont présenté une approche par éléments finis pour l'analyse de
l'érosion et du transport de particules fines dans un sol granulaire soumis à un flux d'infiltration.
L'équation de continuité de la masse des particules transportées a d'abord été dérivée en
considérant un schéma conceptuellement similaire à celui applicable à l'analyse des problèmes
d'écoulement advectif. Ensuite, sa formulation par éléments finis et une procédure d'intégration
temporelle en deux étapes ont été décrites. Une loi incrémentale non linéaire, exprimant le taux
d'érosion en fonction de la vitesse et du temps d'infiltration, a ensuite été proposée sur la base
des résultats d'une série de tests d'érosion en laboratoire récemment présentés dans la littérature.
Cette loi a finalement été adoptée dans la simulation par éléments finis des tests expérimentaux
mentionnés pour un problème bidimensionnel. Sur la base des résultats numériques obtenus, et
de leur accord satisfaisant avec les données expérimentales, ils ont conclu que le modèle

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 44

Chapitre II Revue bibliographique

éléments finis développé pour les besoins de cette étude est potentiellement adapté à l'analyse
des problèmes réels d'érosion. Il convient toutefois de souligner que les données expérimentales
n'étaient pas suffisantes pour déduire les lois régissant le taux de dépôt et la variation de la
conductivité hydraulique. Par conséquent, seule la loi de l'érosion a été introduite dans les
calculs. Le coefficient de conductivité hydraulique a été supposé constant dans le temps et
l'apparition éventuelle de dépôts à certains endroits de l'échantillon n'a pas été prise en compte.

Approche discrète
Les propriétés physiques et mécaniques des matériaux granulaires tels que le sable
peuvent être très différentes de celles des matériaux continus. Cundall et Strack [103] ont
proposé la méthode des éléments discrets (DEM) pour décrire le comportement des matériaux
granulaires. La méthode des éléments discrets (DEM) permet d’étudier le mouvement de
chaque particule individuellement. Dans cette méthode, le mouvement de translation et de
rotation de chaque particule est régi par la seconde loi de Newton, où les forces agissantes sur
chaque particule sont : les forces externes, les forces gravitationnelles, les forces inter-
particulaires et les forces dues au fluide (en cas de présence d’un fluide). Le mouvement des
particules en présence d’un fluide est régi par le système d’équation suivant :

𝜕𝑢
⃑ 𝑓𝑒𝑥𝑡 + 𝑓𝑓𝑙𝑢𝑖𝑑𝑒
= +𝑔
𝜕𝑡 𝑚
𝜕𝜔
⃑ ⃑⃑
𝑀
=
𝜕𝑡 𝐼
Avec ⃑𝑢 est la vitesse de la particule, ⃑𝑓𝑒𝑥𝑡 est la somme des forces extérieures appliquées sur la

particule, ⃑𝑓𝑓𝑙𝑢𝑖𝑑𝑒 est la force appliquée par le fluide sur la particule, 𝑚 est la masse de la particule, 𝑔
⃑⃑

⃑⃑ est la vitesse angulaire de la particule, ⃑⃑⃑

est l’accélération de la pesanteur, ⃑𝜔 𝑀 est le moment appliqué
sur la particule et 𝐼 est le moment d’inertie de la particule.

Couplée à différentes approches pour simuler l'écoulement des fluides à l'échelle des
pores, l'approche par éléments discrets peut représenter assez bien la microstructure et décrire
de plus en plus précisément les mécanismes physiques au sein des matériaux granulaires.
Cependant, son coût de calcul limite les études à des volumes de petites échelles ; la plupart des
études ont été réalisées à l'échelle du volume élémentaire représentatif (VER).

Concernant son application sur le problème d’érosion interne, Zou et al. [99] ont étudié
le mécanisme de transport des particules dans un système filtre-sol (Figure II. 10). Ils ont analysé

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 45

Chapitre II Revue bibliographique

les distributions des particules érodées du sol dans les différents filtres. Les résultats de la
simulation ont montré que la masse érodée et la profondeur de pénétration des particules du sol
dans le filtre sont liées au rapport de taille des particules du sol et du filtre, au gradient
hydraulique et au temps d'érosion. La comparaison des résultats obtenus au critère empirique
de conception du filtre montre que le modèle discret fournit une approche efficace pour étudier
la micro-propriété de filtration et le mécanisme d'érosion dans un système de filtre-sol, ce qui
est utile pour la conception du filtre.

Figure II. 10 Modèle numérique de Zou et al. [99]

Zheng et al. [96] ont effectué une étude micro-macro de la suffusion sur des sols
idéalisés bien gradués et discontinus en utilisant la méthode couplée CFD-DEM (Figure II. 11).
L'interaction entre les particules du sol et l'écoulement d'infiltration est modélisée par l'échange
de quantité de mouvement entre deux phases. Afin d'identifier les changements des états
macroscopiques du squelette du sol, les auteurs ont surveillé l’évolution de la teneur en fines,
la contraction volumétrique et l'indice des vides. Plusieurs phénomènes de colmatage et de
décolmatage, responsables des changements soudains de la teneur en fines et de la réponse du
squelette, sont observés pendant la suffusion. L’étude paramétrique effectuée indique que la
teneur initiale en fines et le gradient hydraulique affectent de manière significative la cinétique
de la suffusion.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 46

Chapitre II Revue bibliographique

Figure II. 11 Migration des particules fines lors de la suffusion pour différents instants : (a) 0s,
(b) 2s, (c) 5s, (d) 30 s [96]

Vandamme et Zou [92] ont étendu leur modèle de particules appelé BMM (Bluff
Morphology Model) pour étudier l'effet de l'infiltration bidimensionnelle sur la stabilité et
l'effondrement des pentes et des levées. Cette étude a mis l’accent sur l'interaction entre
l'érosion, l'infiltration et l'instabilité des pentes en tenant compte de plusieurs phénomènes
complexes tels que : la contribution des sols partiellement saturés et de la succion. L'évolution
temporelle de la stabilité des pentes et de l'effondrement induit par l'infiltration ont été prédits
par le modèle (Figure II. 12). Ces résultats étaient en accord raisonnable avec les résultats
expérimentaux obtenus pour le cas d’un sable homogène non cohésif et les sols cohésifs à
couches multiples. Le rabattement rapide sur un sol sableux est également étudié, et le lieu et
le moment de l'effondrement de la digue sont bien déterminés. Ce modèle permet de bien
localiser les zones érodées ainsi que la quantité de masse érodée due à l'écoulement
d'infiltration.

Figure II. 12 Prédiction du changement du facteur de sécurité par le modèle BMM [92].

Gu et al. [93] ont développé un modèle couplé CFD-DEM pour simuler l’évolution
temporelle du comportement d’un sable argileux pendant l’érosion interne. Dans ce modèle, le

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 47

 Chapitre II Revue bibliographique

couplage particules-fluide est résolu par un schéma de ''grille grossière fixe'' dans le code de
calcul (PFC3D). La suffusion de la matrice argileuse à l'initiation de l'érosion est traduite par
un processus de dégradation de la force de liaison entre les particules selon une loi de
dégradation. Cette loi est dérivée de la loi bien connue du seuil de contrainte de cisaillement
traitant l'érosion interne du sol, qui est basée sur deux paramètres d'érosion : la contrainte de
cisaillement critique et le coefficient d'érosion tel que

𝑚 = 𝑘𝑒𝑟 (𝜏 − 𝜏𝑐𝑟 ) (19)

Avec 𝑚 est la masse érodée, 𝑘𝑒𝑟 est le coefficient d’érosion, 𝜏 est la contrainte de
cisaillement appliquée et 𝜏𝑐𝑟 est la contrainte de cisaillement critique. La Figure II. 13 montre
l’évolution de l’érosion prédite par ce modèle.

Figure II. 13 Evolution du processus d’érosion selon le modèle de Gu et al. [93]

V. Instabilités (ruptures) de pentes liées à l’érosion interne

La rupture d'un barrage en terre est essentiellement contrôlée par l'interaction de deux
processus hydrauliques: le flux d'infiltration à travers le corps du barrage et l'érosion de la face
aval due à l'écoulement de surface [104]. Par ailleurs, l’évaluation de la stabilité d’un ouvrage
en terre nécessite l’introduction de divers paramètres (hydraulique, géométriques et

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 48

Chapitre II Revue bibliographique

sédimentologiques). Gregoretti et al. [104] se sont intéressé dans leur étude à l’amélioration de
la compréhension des conditions hydrauliques et sédimentologiques conduisant à la rupture des
barrages en terre. A cet effet, ils ont réalisé une série logique d'expériences qui étendent les
résultats de Liao et Chou [105] à un certain nombre de configurations de barrages, ainsi qu'une
gamme plus large de tailles de sédiments et de lit de pentes. Ces tests ont permis de distinguer
trois principales typologies de rupture des barrages en terre : La surverse, la rupture par
glissement, rupture par érosion interne. Le mode de rupture dépends des conditions hydraulique
(ex : niveau de retenue d’eau), géométriques (ex : pente amont et aval du barrage) et des
matériaux (ex : plasticité).

Plusieurs études ont spécifiquement porté sur les mécanismes physiques qui déterminent
la rupture des barrages en terre. Les barrages de glissement sont habituellement constitués d'une
masse hétérogène de matériaux non consolidés ou mal consolidés et peuvent donc être soumis
à des écoulements d'infiltration intenses. Plusieurs études expérimentales [105–107] et
observations in situ [108–110] ont montrés que la rupture des barrages de glissement peut se
produire non seulement par la surverse mais aussi par érosion régressive qui se produit au pied
de la face aval du barrage et migre progressivement en amont ou par la combinaison des deux
mécanismes. L’érosion interne par suffusion et le renard sont aussi des processus d’érosion par
infiltration qui peuvent conduire à la rupture des barrages [110–112].

De multiples études définissent l'effondrement des pentes ou le glissement de terrain dû

à l'infiltration d'eau par différents termes tels que le « Piping », Sapping ou Spring sapping,
érosion interne, érosion d’infiltration ou souterraine et tunnel scouring [29, 47, 107, 113–119].
Certains chercheurs ont toutefois établi des distinctions entre les différents mécanismes
impliqués dans l'instabilité causée par les infiltrations [117, 120].

Crosta et di Prisco [120] ont présenté dans leur étude les phénomènes antécédents de
l’effondrement d’une pente qui s’est produit le long d’une ancienne terrasse fluviale dans la
Plaine de Po (Nord de l’Italie). La rupture est caractérisée par une chaîne complexe de
phénomènes (érosion, glissement, liquéfaction) se déroulant au cours d'une certaine période de
temps. En effet, l'évolution de l'infiltration avec le temps a joué le rôle de déclencheur de tous
ces phénomènes. Crosta et di Prisco relatent principalement cet effondrement dû à l’érosion qui
a débutée en saison de culture à cause des techniques agricoles. La présence d'incisions
superficielles causées par les roues des tracteurs en forme de disques tranchants qui coupent la
couche de sol superficielle a induit une infiltration d'eau localisée (Figure II. 14). Dans ce cas,

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 49

Chapitre II Revue bibliographique

ils favorisaient l'infiltration dans des couches plus profondes, non saturées et relativement plus
perméables. Afin de clarifier l'évolution de l'infiltration et les phénomènes d'érosion qui
précédaient et caractérisaient l'effondrement, des analyses numériques et des essais ont été
réalisés sur des modèles à petite échelle.

Figure II. 14 Coupe transversale de la terrasse fluviale telle qu'elle a été obtenue à partir des
investigations sur le terrain [120]

La rupture catastrophique en 1993 du barrage de Gouhou de 71 m de hauteur en province

de Qinghai en Chine a donné lieu à de nombreuses études visant à identifier les causes de la
rupture [5, 121, 122]. Selon le groupe d’investigation de la rupture du barrage de Gouhou [123],
le barrage soufrerait de plusieurs problèmes structuraux. Ces défauts structuraux ont favorisé
l’apparition des canaux de fuite. L’érosion interne concentrée dans un conduit qui en a résulté
a été considérée comme la cause directe qui a déclenché la rupture du barrage. Zhang et Chen
[124] ont apporté de nouvelles explications sur les modes et les processus de la rupture. Sur la
base des conditions géométriques des matériaux impliqués dans l’évènement et les conditions
hydrauliques dans le barrage, la susceptibilité à l’érosion interne a été évaluée par les auteurs.
Selon le critère de Kenney et Lau [31] et de [41], les trois matériaux du barrage étaient
susceptibles à l’érosion interne. L’analyse des conditions hydrauliques dans le barrage a été
réalisée en utilisant un programme en éléments finis SVFlux avec le code de résolution des
équations aux dérivées partielles FlexPDE. Après l’analyse des résultats, le processus possible
de la rupture du barrage a été présenté (Figure II. 15).

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 50

 Chapitre II Revue bibliographique

(1) Montée du niveau de l'eau (4) développement progressif d’un conduit

(2) Infiltration de l’eau dans le barrage (5) L'affouillement de la pente et la chute
de cailloux
(3) Initiation de l’érosion interne (6) Formation de la brèche

Figure II. 15 Processus de rupture du barrage de Gouhou en Chine [124]

VI. Conclusion

Les principaux résultats des recherches expérimentales et numériques sur l'érosion

interne ont été passés en revue dans ce chapitre. Cependant, peu d'études numériques à l'échelle
d'une structure entière sont disponibles dans la littérature [80, 84, 85, 88, 100]. De plus, les
réponses mécaniques d'un sol érodé ont souvent été négligées ou simplifiées en utilisant un
modèle élastique avec des paramètres de dommages induits.

Par conséquent, dans le but d'améliorer la prévision de la suffusion des sols et de faciliter
la détection fiable et l'alerte précoce du dysfonctionnement éventuel des structures
hydrauliques, des infrastructures urbaines, des ruptures de pente et des glissements de terrain,
il est nécessaire de développer une approche qui puisse mieux prendre en compte les
phénomènes hydromécaniques couplée avec la suffusion dans les problèmes de conception et
d'ingénierie à l'échelle réelle.

Le chapitre suivant concerne la présentation de la partie expérimentale de cette étude,

qui a pour objectif de déterminer l’effet de l’érosion interne sur les paramètres de résistance au
cisaillement des sols.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 51

 Chapitre III

Etude expérimentale
 Chapitre III Etude expérimentale

I. Introduction

Pour étudier le phénomène d’érosion interne au laboratoire, plusieurs dispositifs

expérimentaux ont été développé [16, 125–128]. La plupart de ces études portent sur
l’évaluation de la susceptibilité à l’érosion des matériaux étudiés, et sur l’étude des conditions
géométriques et hydrauliques conduisant à l’apparition de l’érosion. Ces dernières années,
plusieurs chercheurs se sont investis dans l’étude du comportement hydromécanique des sols
érodables [38, 71, 129–133]. Ils ont ainsi proposé différentes démarches expérimentales
pouvant évaluer la réponse hydromécanique permettant d’estimer l’effet de l’érosion interne
sur la stabilité des ouvrages en terre.

Dans ce chapitre, nous allons présenter la démarche expérimentale adoptée et les

résultats obtenus pour évaluer l’influence de l’érosion interne sur la résistance au cisaillement
des sols. Cette étude a été effectuée au niveau du laboratoire ondes et milieux complexes de
l’université Le Havre Normandie en France. Une série de tests a été effectuée sur des
échantillons reconstitués au laboratoire et composés principalement de sable mélangé avec des
sols fins (Kaolinite, limon et sable fin). Au total neuf échantillons ont été étudiés. Une
comparaison des caractéristiques avant et après érosion interne des sols sera présentée plus loin.

I. Matériaux et méthode

L’approche expérimentale retenue dans cette étude consiste à étudier le comportement

au cisaillement d’un sol soumis à l’érosion interne. L’essai est composé de deux parties
essentielles :

1. Le cisaillement d’un échantillon de sol et détermination de ces paramètres de

résistance au cisaillement.
2. L’application d’un écoulement interne afin de provoquer une érosion interne,
ensuite, le cisaillement de l’échantillon. L’objectif étant de comparer les
paramètres de résistance au cisaillement d’un sol avant et après érosion.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 52

 Chapitre III Etude expérimentale

II. Description des essais

Une procédure d'essai en plusieurs étapes a été menée pour évaluer les conditions
conduisant à la modification de la résistance du sol soumis à l'érosion interne, évaluée par la
boîte de cisaillement directe. La procédure d'essai détaillée est effectuée comme suit pour
chacun des neuf mélanges de sol :

(1) Confection de l’échantillon

(2) Après la confection de l'échantillon, celui-ci a été placé dans la boîte de cisaillement,
la boite elle-même a été placée sur l'appareil et un essai de cisaillement a été effectué afin de
déterminer les paramètres de résistance à l’état initial.

(3) Un échantillon identique a été fabriqué et placé dans la boîte de cisaillement. Avant
le test, l'échantillon était saturé (écoulement descendant à faible gradient) et la hauteur d'eau
dans la cuve était ajustée à une valeur constante. Ensuite, l'infiltration a commencé.

(4) L’écoulement a été maintenu pendant une heure pour assurer l'achèvement de
l'érosion interne, c'est-à-dire que le débit est stable et que les enregistrements de turbidité
deviennent constants.

(5) L'échantillon est ensuite soumis à un test de cisaillement direct comme décrit dans
l'étape 1 afin de déterminer les paramètres de résistance après érosion.

(6) Les étapes de (1) à (5) sont répétées pour chaque charge normale (verticale) (trois
contraintes normales ont été appliquées : 50, 100 et 200 kPa).

III. Matériaux étudiés

Le choix des matériaux est une étape clé dans l’étude expérimental de l’érosion interne.
On trouve dans la littérature trois types de matériaux essentiels à savoir : les matériaux naturels,
les matériaux artificiels et les matériaux reconstitués. Le choix du matériau dépend
principalement de l’objectif de l’étude. Les matériaux naturels sont des matériaux qui se
retrouve dans la nature à leur état réel. Ils sont utilisés dans les études des cas de rupture des
ouvrages spécifiques. Le but de leur utilisation est d’identifier les scénarios de ruptures et de

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 53

Chapitre III Etude expérimentale

les caractériser. Les matériaux artificiels utilisés dans la littérature pour étudier l’érosion interne
sont des mélanges de billes de verre de différentes tailles et de proportions. Les proportions
sont choisies de manière à représenter une granulométrie donnée. Bien que ces matériaux soient
peu représentatifs de la réalité, mais leur utilisation est largement répandue dans la littérature.
Le troisième type de matériau correspond aux sols reconstitués. Il s’agit de mélanger des
matériaux naturels de référence pour former un nouveau matériau. En jouant sur les proportions,
on peut reconstituer des matériaux de n’importe quelle granulométrie souhaitée. Ce type de
matériaux est le plus utilisé dans les études de l’érosion interne parce qu’il permet de fabriquer
une large gamme de matériaux qui peuvent représenter des matériaux réels.

III.1 Reconstitution des matériaux

La reconstitution des matériaux permet de fabriquer des matériaux de mêmes
caractéristiques qu’un matériau naturel susceptible à l’érosion interne.

Les sols testés dans la présente étude ont été obtenus à partir du mélange du Sable#1
(sable de Fontainebleau NE34) avec différentes fractions de sols fins. Afin d’élargir notre étude
et de mieux comprendre l'effet de l'érosion interne sur les caractéristiques de résistance au
cisaillement, trois types de particules fines (kaolinite, limon et sable# 2) ont été mélangés avec
le Sable#1 à différentes teneurs (10 %, 15 %, 20 %). Au total, neuf (9) matériaux ont été utilisés.
Le Sable#2 a été obtenue à partir de sable NE34 passant dans le tamis de 200 µm. Le Sable#1
représente la matrice granulaire tandis que la kaolinite, le limon et le sable#2 fonctionnent
comme les fines érodables. Tous les échantillons ont été compactés dans la boîte de cisaillement
à une densité sèche de 1.6 et une teneur en eau de 16 %. La Figure III. 1 montre les courbes
granulométriques des matériaux de référence et les figures (Figure III. 2, Figure III. 3, Figure III.
4) montrent les courbes granulométriques des matériaux reconstitués. Le Tableau III. 1 et le
Tableau III. 2 montrent les caractéristiques physico-chimiques des différents matériaux utilisés.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 54

Chapitre III Etude expérimentale

100

90 Kaolinite
80 Sable#1
Limon
70
Passant (%)

Sable#2
60

50

40

30

20

10

0
0.01 0.1 1 10 100
Diamètre (µm)

Figure III. 1 Courbes granulométrique des matériaux de base

100

90 90% Sable#1+10% Kaolinite

80
85% Sable#1+15%Kaolinite
70
Passant (%)

80% Sable#1+20% Kaolinite

60

50

40

30

20

10

0
0.1 1 10 100 1000
Diamètre (µm)

Figure III. 2 Courbes granulométriques des matériaux à base de Kaolinite et de Sable NE34
(Sable#1)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 55

Chapitre III Etude expérimentale

100
90 90% Sable#1+10% Limon
80
85% Sable#1=15% Limon
70
Passant (%)

80% Sable#1=20% Limon

60
50
40
30
20
10
0
0.1 1 10
Diamètre (µm) 100 1000

Figure III. 3 Courbes granulométriques des matériaux à base de Limon et de Sable NE34
(Sable#1)
100

90 90% Sable#1+10% Sable#2

80
85%Sable#1+15% Sable#2
70
Passant (%)

80% Sable#1+20% Sable#2

60

50

40

30

20

10

0
0.1 1 10 100 1000
Diamètre (µm)

Figure III. 4 Courbes granulométriques des matériaux à base de Sable fin (Sable#2) et de Sable
NE34 (Sable#1)

Tableau III. 1 Caractéristiques physicochimiques et composition chimique du sable de

fontainebleau (NE34)

Minéralogie SiO2 Al2O3 Fe2O3 TiO2 CaO K2O

(%) >99.700 <0.130 <0.030 <0.033 <0.03 <0.035
Poids Dry: 16 kN/m3, Solid: 26.5 kN/m3
volumique

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 56

Chapitre III Etude expérimentale

Autres pH=7; angularity coefficient = 1.1; Loss on ignition: <0.05%

Tableau III. 2 Caractéristiques physicochimiques et composition chimique des fines utilisés

Minéralogie (%)
SiO2 Al2O3 Fe2O3 TiO2 CaO MgO Na2O K2O
Silt 55.00 18.00 9.00 0.74 2.00 2 – 10.00
Kaolinite 56.40 28.70 1.40 1.20 0.30 0.30 0.10 0.80
Limite d’Atterberg MBV Sa
WL (%) WP (%) IP (g/100g) (m2/g)
Silt 32.0 13.0 19.0 1.5 _
Kaolinite 57.0 32.0 25.0 3.5 70.0

III.2 Préparation des échantillons

Les échantillons ont été compactés dans la boîte de cisaillement en trois couches avec
la méthode de compactage humide. Cette méthode permet d'obtenir une densité de sol uniforme
en compactant les différentes couches à une densité inférieure à la valeur souhaitée (1.6 g/cm3)
car le compactage d'une couche supérieure densifie également les couches inférieures ; c'est le
concept de "sous-compactage" qui a été prouvé fiable par d'autres chercheurs [134, 135]. La
teneur en eau optée pour la préparation des échantillons est de 16%. Cette valeur est choisie
afin de permettre une bonne mise en place des différents matériaux. En effet, on a pu constater
lors des essais expérimentaux qu’une teneur en eau faible nécessite une énergie de compactage
importante d’où la difficulté de la mise en place de l’échantillon. Tandis qu’une teneur en eau
importante ramène l’échantillon à son état de liquidité où la densité fixée ne peut être respectée
(c-à-d. le volume de l’échantillon se réduit au-dessous du volume de la boite).

III.3 Essai de cisaillement à la boite

Les principaux types de sollicitation qu’on peut rencontrer dans le sol sont : le
cisaillement, la compression et la traction.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 57

Chapitre III Etude expérimentale

Figure III. 5 Les trois principaux types de sollicitations

Etant donné que la cause principale de sollicitation dans les sols est la gravité terrestre,
créant des efforts verticaux dirigés vers le bas, on comprend intuitivement que les contraintes
de traction sont en générales absentes. De plus, on sait que la résistance à la traction des sols
est très faible voire nulle.

Par ailleurs, la résistance d'un sol à la compression est supérieure à sa résistance au

cisaillement. Ainsi pour dimensionner des ouvrages de génie-civil et d’hydraulique
(Fondations, Barrage en terre, Digues...), il est indispensable de connaître la résistance au
cisaillement du sol, car s'il y a rupture ce sera suivant ce type de sollicitation. Cependant, les
calculs de l’ingénierie géotechnique négligent les effets de l’érosion interne sur ces paramètres.
C’est pour cette raison que nous avons choisi d’étudier l’effet de l’érosion interne sur ces
paramètres.

Il existe différents dispositifs expérimentaux pour mesurer au laboratoire les paramètres

de résistance au cisaillement des sols qui sont : l’essai de cisaillement à la boite ou essai de
Casagrande, essai triaxial et l’essai de scissomètre.

Afin d’évaluer les paramètres de résistance au cisaillement avant et après érosion, nous
avons modifié la boite de Casagrande de manière à appliquer un écoulement à travers un
échantillon de sol sous une charge hydraulique constante. Le schéma explicatif du dispositif est
montré par la Figure III. 8Figure III. 8.

Principe de l’essai de cisaillement

Cet essai a pour objectif de déterminer les caractéristiques de résistance au cisaillement
d’un sol donné.

L’échantillon est mis dans une boite, la boite elle-même est composée de deux demi
boites, l’échantillon est soumis à une contrainte verticale (𝜎) avant d’être cisailler à une vitesse
constante. La contrainte de cisaillement (𝜏) appliquée s’accroît jusqu’à une valeur maximale.
On peut donc déduire l’état de contrainte de l’échantillon à la rupture (𝜎𝑟 , 𝜏𝑟 ). La Figure III. 6
montre un schéma explicatif du principe de cet essai.

Pour pouvoir déterminer le critère de rupture de Mohr-Coulomb et définir ses valeurs

de cohésion et d’angle de frottement, il faut réaliser plusieurs essais sur un même matériau avec

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 58

Chapitre III Etude expérimentale

différentes contraintes verticales. Dans notre cas, nous avons réalisé trois essais pour un
matériau et on a pris la moyenne des valeurs obtenues.

Description de la boite de Casagrande

L'échantillon de sol est placé à l'intérieur de deux demi-boites qui peuvent se déplacer
horizontalement l'une par rapport à l'autre. Un piston permet d'exercer sur le sol un effort normal
N constant pendant toute la durée de l'essai.

Figure III. 6 Principe de l’essai de cisaillement directe

Une demi-boîte est entraînée horizontalement à vitesse constante. A tout instant, on

mesure la force de cisaillement T. Un second comparateur vertical permet de mesurer la
variation de hauteur ∆ℎ de l'échantillon.

On exerce sur le plan de séparation des deux demi-boites une contrainte dont les
composantes normale et tangentielle ont pour valeur moyenne :

𝑁 𝑇
𝜎=𝐴, 𝜏=𝐴
𝑐 𝑐

Où 𝐴𝑐 est la section corrigée de l'eprouvette (cette section variant en fonction du

déplacement relatif ∆𝑙).

Un schéma global du bâti de chargement sur lequel est fixée la boîte de Casagrande est
présenté à la Figure III. 7.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 59

Chapitre III Etude expérimentale

Figure III. 7 Appareil de cisaillement

1 : bâti de la machine

2 : Bloc moteur

3 : Boite de cisaillement

4 : Anneau dynamométrique

5 : Comparateur de l’anneau

6 : Etrier de chargement

7 : comparateur de mesure de ∆ℎ

8 : Plateau supportant les masses

9 : Les masses

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 60

Chapitre III Etude expérimentale

10 : Volant d’approche de l’anneau

11 : Bague de contrôle

12 : Volant de réglage

III.4 Essai érosion-cisaillement

Plusieurs dispositifs expérimentaux ont été proposés par plusieurs chercheurs pour
l’étude de l’érosion interne des sols [16, 17, 126, 136–138]. La majorité de ces essais avait pour
objectifs d’étudier la susceptibilité des matériaux granulaires à l’érosion interne. Cependant, il
existe peu d’études permettant d’évaluer le comportement post-érosion des sols. Parmi les
chercheurs qui ont largement investi dans cette thématique, on trouve Takahashi qui visent dans
ces différents travaux expérimentaux à mettre l’accent sur l’effet de l’érosion interne sur la
résistance des sols granulaires [75, 129]. Dans ses études, Takahashi utilise un appareil triaxial
modifié pour le suivi du comportement des matériaux étudiés.

Cette étude constitue une alternative à l’étude de l’effet l’érosion interne sur les
paramètres de résistance au cisaillement des sols sableux. Les essais ont été réalisés sur un
appareil de cisaillement à la boite de Casagrande présenté ci-dessus. Pour atteindre notre
objectif, nous avons apporté une modification à la boite de cisaillement de manière à appliquer
une charge hydraulique avant de procéder au cisaillement de l’échantillon. Nous avons donc
placé un réservoir d’alimentation en eau qui permet l’application de la charge hydraulique.

Comme décrit ci-dessus, deux types de tests ont été effectués : un simple test de
cisaillement direct seul, et un test séquentiel de cisaillement et d'érosion où l'échantillon a été
cisaillé après le processus d'érosion pour connaître son effet de cette dernière sur la résistance
du sol au cisaillement. Le test est effectué sur un échantillon de sol mis dans la boîte de
Casagrande modifiée de section carrée (6 cm de large). Le sol est cisaillé horizontalement le
long d'un plan prédéterminé en déplaçant la moitié supérieure de la boîte par rapport à la partie
inférieure (Figure III. 6). L'échantillon est soumis à une contrainte normale contrôlée et la partie
supérieure de l'échantillon est tirée latéralement à une vitesse de déformation contrôlée (1,2
mm/min selon la norme française NF P94-071-1) et la charge latérale appliquée a été enregistrée
en continu.

Dans le test de cisaillement-érosion, afin de protéger l'échantillon contre l'érosion de

surface, un géotextile est placé sur la surface supérieure. Une charge hydraulique constante de

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 61

 Chapitre III Etude expérimentale

0,4 m a été appliquée sur l'échantillon pendant une heure pour provoquer l'érosion interne, puis
l'échantillon est maintenu au repos pendant une heure avant qu'un test de cisaillement direct ne
soit effectué comme décrit précédemment. Pour chaque échantillon, une courbe de Mohr-
coulomb est tracée afin d'en déduire les valeurs de cohésion et d'angle de frottement interne.

Boite de Casagrande modifiée pour l’érosion interne

Des tests à charge hydraulique constante avec un écoulement descendant ont été
effectués pour provoquer l’érosion interne des échantillons étudiés. La Figure III. 8 montre le
schéma explicatif du dispositif expérimental modifié pour permettre d’applique la charge
hydraulique. L’échantillon est mis dans la boite de 6 cm de côté et de 4 cm de profondeur, et la
charge hydraulique est fixée à 400 mm. Un géotextile est placé sur la face supérieure de
l’échantillon pour éviter son érosion à la surface, elle permet ainsi d’uniformiser l’écoulement
et éviter ainsi l’apparition d’un trou en surface.

Figure III. 8 Schéma illustratif du dispositif expérimental de l’essai d'érosion-cisaillement

IV. Evaluation de la susceptibilité à l’érosion

Afin de s'assurer que l'érosion interne se produira pendant les tests d'infiltration, la
susceptibilité à la suffusion des neuf mélange a été évaluée par cinq critères suggérés par Kezdi
[42], Istomina [24], Mao [26], Burenkova [25] et Wan&Fell [43]. Les résultats sont résumés
dans le Tableau III. 3. Selon la plupart de ces critères, les différents sols testés sont instables et
vulnérables à l'érosion interne.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 62

 Chapitre III Etude expérimentale

Tableau III. 3 Évaluation de la vulnérabilité des spécimens à l'érosion interne selon les
différents critères

Spécimen Critère d’évaluation de la stabilité interne

Kezdi Istomina Mao [26] Burenkova Wan & Fell
[42] [24] [25] [43]
10% Kaolinite U S U S U
15%Kaolinite U S U S U
20%Kaolinite U U U U U
10%Limon U S U S U
15%Limon U S U U U
20%Limon U S U U U
10%Sable#2 S S U S U
15%Sable#2 S S U S U
20%Sable#2 S S U S U
U : Instable, S : Stable

V. Résultats et discussion

Les échantillons ont été testés sous différentes contraintes verticales (55, 110 et 222
kPa). La valeur de la contrainte de cisaillement maximale à la rupture du sol est tracée en
fonction de la contrainte normale pour chaque test. Selon le critère de rupture de Mohr-
Coulomb, la ligne qui représente l'enveloppe de rupture peut être écrite comme suit

𝜏 = 𝐶 + 𝜎 tan 𝜑 (20)
Avec :

C : est la cohésion du sol

𝜑 : est l’angle de frottement interne du sol

τ : est la contrainte de cisaillement

σ : est la contrainte verticale

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 63

Chapitre III Etude expérimentale

V.1 Avant érosion

Les résultats des tests de cisaillement direct avant érosion sont présentés sur les figures
(Figure III. 9, Figure III. 10, Figure III. 11) pour les différents matériaux testés impliquant
différents types et teneurs en fines. Selon ces résultats, on remarque que plus on augmente la
teneur en fine plus le sol est cohérent, alors que l'angle de frottement tend à diminuer lorsque
la teneur en fines augmente. Cela implique que les particules fines ajoutées au sol lui offrent
une meilleure cohésion.

La comparaison des différents types de fines montre qu’à 10% de taux de fines, le
mélange à base de limon rend le sol beaucoup plus cohérent que la kaolinite et le sable fin,
tandis que, à partir de 15% de fines, le mélange à base de Kaolinite devient le plus cohérent.
Concernant l’angle de frottement interne, pour les taux de fines de 10 et 15%, le mélange de
kaolinite avait les valeurs les plus grandes, bien qu’à 20% de fines, le mélange à base de limon
avait l’angle de frottement le plus important.

200

150
τ (kPa)

100

y = 0,7371x + 17,917 10% Kaolinite

50 y = 0,6834x + 19,708 15% Kaolinite
y = 0,6066x + 20,903 20% Kaolinite
0
0 50 100 150 200 250
σ (kPa)

Figure III. 9 Résultat de l’essai de cisaillement avant érosion pour le mélange de Kaolinite et
sable NE 34

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 64

Chapitre III Etude expérimentale

200

τ (kPa) 150

100

y = 0,6496x + 18,514 10% Limon

50 y = 0,6266x + 19,111 15% Limon
y = 0,6081x + 20,306 20% Limon
0
0 50 100 150 200 250
σ (kPa)

Figure III. 10 Résultat de l’essai de cisaillement avant érosion pour le mélange de Limon et
sable NE 34

200

150
τ (kPa)

100

y = 0,625x + 16,125 10% Sable#2

50
y = 0,582x + 17,319 15% Sable#2
y = 0,5713x + 17,917 20% Sable#2
0
0 50 100 150 200 250
σ (kPa)

Figure III. 11 Résultat de l’essai de cisaillement avant érosion pour le mélange de sable fin
(Sable#2) et sable NE 34

V.2 Effet de la teneur en fine sur la cohésion et l’angle de frottement

interne
L’effet de la teneur en fine sur la cohésion et l’angle de frottement interne des
échantillons étudiés est représenté sur les figures (Figure III. 12, Figure III. 13). On remarque que
les deux paramètres présentent une relation polynomiale avec la teneur en fine. On remarque
sur ces deux figures que pour une teneur en fines nulle, l’échantillon de sable seul à une teneur

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 65

Chapitre III Etude expérimentale

en eau de 16% présente une cohésion de 9.6 kPa et un angle de frottement interne de 38.2°. Ces
mêmes valeurs ont été obtenues expérimentalement pour le sable seul.

25

20
Cohésion (kPa)

15

y = -0,0253x2 + 1,0671x + 9,6 Kaolinite

10 R² = 0,999
y = -0,0316x2 + 1,1535x + 9,6 Limon
5 R² = 0,9886
y = -0,0226x2 + 0,8628x + 9,6 Sable#2
R² = 0,9989
0
0 5 10 15 20
Teneur en fine (%)

Figure III. 12 Effet de la teneur en fines sur la cohésion

40
38
36
34
32
φ (°)

30
28 y = -0,0172x2 - 0,0042x + 38,2 Kaolinite
R² = 0,9999
26 y = 0,0162x2 - 0,6651x + 38,2 Limon
24 R² = 0,997
22 y = 0,02x2 - 0,8275x + 38,2 Sable#2
R² = 0,9997
20
0 5 10 15 20
Teneur en fines (%)

Figure III. 13 Effet de la teneur en fines sur l’angle de frottement interne

V.3 Résultats post-érosion

Les résultats des tests de cisaillement direct après avoir subir les échantillons à une
érosion interne sont présentés sur les figures (Figure III. 14, Figure III. 15) pour les différents
matériaux testés impliquant différents types et teneurs en fines.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 66

Chapitre III Etude expérimentale

Les résultats montrent que les échantillons qui avaient plus de teneur en fine reste
toujours plus cohérent après érosion. On remarque également que la cohésion du mélange de
limon se dégrade beaucoup plus que les autres mélanges. On peut également constater que le
taux de diminution de l'angle de frottement est plus évident pour la kaolinite que pour les autres
fines. On a constaté pour le mélange de (Sable#1+Sable#2) que l’angle de frottement interne
tend à augmenter après érosion. Le Tableau III. 4 résume les valeurs de la cohésion et de l'angle
de frottement interne avant et après érosion pour les différents matériaux.

200

150
τ (kPa)

100

y = 0,5582x + 17,083 10% Kaolinite

50 y = 0,6834x + 19,708 15% Kaolinite
y = 0,6066x + 20,903 20% Kaolinite
0
0 50 100 150 200 250
σ (kPa)

Figure III. 14 Résultat de l’essai de cisaillement après érosion pour le mélange de Kaolinite et
sable NE 34

160
140
120
100
τ (kPa)

80
60 y = 0,5739x + 16,528 10% Limon
40 y = 0,5507x + 16,944 15% Limon
20 y = 0,5221x + 17,222 20% Limon

0
0 50 100 σ (kPa)150 200 250

Figure III. 15 Résultat de l’essai de cisaillement après érosion pour le mélange de Limon et
sable NE 34

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 67

Chapitre III Etude expérimentale

Afin de pouvoir représenter graphiquement la comparaison entre les paramètres de

résistance au cisaillement avant et après érosion interne, le taux de réduction de la résistance
est défini par le rapport suivant :

tan 𝜑̅
∆𝑅𝜑 = [1 − ] ∗ 100 (21)
tan 𝜑0
tan 𝐶̅ (22)
∆𝑅𝐶 = [1 − ] ∗ 100
tan 𝐶0
Où ∆𝑅𝜑 est la réduction relative de l'angle de frottement interne, ∆𝑅𝐶 est la réduction
relative de la cohésion, 𝜑
̅ et 𝜑0 sont l'angle de frottement interne après et avant l'érosion
̅ et 𝐶0 sont la cohésion après et avant l'érosion respectivement.
respectivement, 𝐶

Tableau III. 4 Récapitulatif des résultats de l'évolution de la cohésion et de l'angle de frottement

interne dû à l'érosion interne

Spécimen Cohésion (kPa) Angle de frottement interne (°)

Avant Après ΔRC (%) Avant Après ΔR𝜑 (%)
10 % (K) 17.9 17.1 4.46 36.4 29.2 24
15 % (K) 19.7 17.9 9.13 34.3 28 22
20% (K) 20.9 18.5 11.48 31.2 26.8 16
10 % (SI) 18.5 16.5 10.81 33 29.8 11.89
15 % (SI) 19.1 16.9 11.51 32.1 28.8 12.36
20 % (SI) 20.3 17.2 15.27 31.3 27.6 14
10 % (SA) 16.1 15.4 4.04 32 33.2 -4.72
15 % (SA) 17.3 15.8 8.67 30.2 32.1 -7.78
20 % (SA) 17.9 16.1 10.05 29.7 30.9 -4.92
K : Kaolinite, SI : Limon, SA : Silica#2

Les résultats de la variation de la cohésion et de l'angle de frottement interne due à

l'érosion interne sont montrés sur le Tableau III. 4. On peut noter qu'après l’érosion interne, les
échantillons ayant une teneur en fines initiale plus élevée avaient la cohésion la plus élevée et
l'angle de frottement interne le plus faible. On peut également constater que pour tous les
échantillons, plus la teneur en fines initiales est élevée, plus la cohésion diminue après l’érosion.
Des résultats similaires ont été rapportés par Ke et Takahashi [139]. Cette diminution de la
cohésion est plus importante pour le mélange de limon (Tableau III. 4) et devient plus importante

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 68

Chapitre III Etude expérimentale

lorsque la teneur en fines augmente. En ce qui concerne l'angle de frottement interne, on peut
noter qu'il diminue après l'érosion pour tous les échantillons testés, sauf pour le sable#2.

L'étude indique une faible dégradation de l'angle de frottement pour le mélange

contenant plus de kaolinite, alors que la tendance inverse a été démontrée pour le limon. Par
conséquent, l'érosion interne influe plus fortement sur l'angle de frottement des mélanges
contenant plus de limon, alors que la dégradation de l’angle de frottement interne est légèrement
affectée par la teneur en kaolinite.

Pour rendre les résultats du Tableau III. 4 plus clair, la réduction de la cohésion en
fonction de la variation de la teneur en fines est illustrée sur la Figure III. 16 où une relation
linéaire entre la réduction de la cohésion et la teneur en fines est mise en évidence.

Figure III. 16 Évolution de la réduction de la cohésion en fonction de la teneur en amendes

La régression linéaire présentée à la Figure III. 16 montre que le taux de réduction de la

cohésion a une relation linéaire en fonction de la teneur en fines pour tous les mélanges.
Cependant, la variation de l'angle de frottement interne ne fournit pas une tendance bien définie.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 69

Chapitre III Etude expérimentale

Conclusion

Dans ce chapitre, nous avons présenté une étude expérimentale de l’influence de

l’érosion interne sur les paramètres de résistance au cisaillement (cohésion et angle de
frottement interne). L’étude a été réalisée sur l’appareil de cisaillement directe en modifiant la
boite de Casagrande de manière à pouvoir appliquer une charge hydraulique. Les matériaux
testés dans cette étude sont des matériaux reconstitués au laboratoire à partir de mélange de
différents matériaux de référence. L’étude de l’effet de la teneur en fines sur les paramètres de
résistance au cisaillement a été investie. Les résultats ont révélé que, l’augmentation de la teneur
en fine offre une meilleure cohésion au sol, tandis que l'angle de frottement tend à diminuer.
En comparant les résultats avant et après le test d’érosion, on remarque que la cohésion et
l’angle de frottement se dégradent à la fin du processus. L’analyse des résultats post-érosion-
interne a permis de constater que les échantillons ayant une teneur en fines initiale plus élevée
avaient la cohésion la plus élevée et l'angle de frottement interne le plus faible. On peut
également constater que pour tous les échantillons, la cohésion des sols qui avaient la teneur en
fines initiales élevée se dégrade beaucoup plus après l’érosion.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 70

 Chapitre IV

Etude numérique : Modèle

1D
 Chapitre IV Etude numérique : Modèle 1D

I. Introduction

Ce chapitre sera consacré à la modélisation du processus d’érosion interne en

s’appuyant sur la théorie des milieux poreux. Le modèle mathématique se base sur les mêmes
idées que les modèles de Vardoulakis et al. [85], de Stavropoulou et al. [84] et de Papamichos
et al. [83] qui traitent le problème de production de sable.

Le modèle présenté ici apporte une amélioration sur la loi d’érosion par la définition
d’un seuil d’érosion à travers un paramètre ncr appelé « porosité critique ». Ce paramètre influe
sur le processus d’érosion en considérant qu’au-delàs d’une certaine porosité, le processus
devient stable.

Une étude paramétrique montrant l’effet d’un nombre de paramètres sur la cinétique de
l’érosion est ensuite présentée, et l’effet de chaque paramètre est discuté. Nous allons, par la
suite, confronter les résultats expérimentaux du chapitre 2 au modèle mathématique en
déterminant pour chaque matériau, son coefficient d’érosion ‘λ’ et sa porosité critique ‘ncr’. Les
résultats seront présentés et discutés.

II. Modélisation du processus d’érosion interne

L’érosion du sol est modélisée par un transfert de masse entre la matrice solide et le
fluide interstitiel en se basant sur la théorie des milieux poreux. Les équations différentielles
utilisées pour modéliser ce transfert de masse en milieu poreux sont celles proposées par
Vardoulakis et al. [85] avec une modification de la loi gouvernant la variation de la porosité du
milieu. En effet, le modèle traite le problème d’érosion interne par suffusion où la porosité subit
une légère augmentation contrairement au problème d’érosion de trou traité dans le modèle de
Vardoulakis et al. [85] où la porosité atteint sa valeur limite de 100% (toute la masse solide est
érodée). La variation dans le temps et dans l'espace de la porosité, résultant de l'érosion interne,
est choisie comme paramètre de couplage. La cohésion du sol est reliée à la porosité du sol par
une loi linéaire de manière à ce que le matériau se dégrade à mesure que la porosité augmente.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 71

 Chapitre IV Etude numérique : Modèle 1D

III. Ecoulement dans les milieux poreux

L'écoulement à travers un milieu poreux est le principal mécanisme de nombreux

processus naturels dans les applications d'ingénierie, comme l'écoulement des eaux
souterraines, l'écoulement à travers les structures en terre (c'est-à-dire les barrages en remblai,
les digues et les levées) et tout processus qui implique des infiltrations [140].

La dynamique des fluides dans les milieux poreux est dominée par la structure de la
matrice poreuse qui occupe tout le volume et qui permet une dissipation efficace de l'énergie
cinétique du fluide. Le mouvement des fluides dans les milieux poreux peut être décrit par la
loi de Darcy (Equation. 23111), qui est l'équation de quantité de mouvement la plus
couramment utilisée. L'écoulement à travers un milieu poreux est défini comme une relation
linéaire entre la perte de charge de l'eau et la vitesse moyenne du fluide, proportionnée par la
conductivité hydraulique comme propriété du milieu :

𝑘
𝑣= (∇𝑝 + 𝜌𝑔)
𝜇 (23)

Où k est la perméabilité intrinsèque, v est la vitesse de Darcy, ∇𝑝 est le gradient de

pression, µ est la viscosité cinématique du fluide, ρ est la masse volumique du fluide et g est
l’accélération de la pesanteur.

La loi de Darcy est utile pour modéliser les d'écoulements souterrains saturés en régimes
laminaires.

IV. Théorie des milieux poreux

La théorie des milieux poreux (TMP) est utilisée dans les branches de l'ingénierie qui
traitent de l'écoulement des fluides dans le sol, les roches fracturées ou d'autres types de milieux
poreux où les problèmes ne peuvent pas être classés avec précision en mécanique des solides
ou en mécanique des fluides. La TMP a été introduite pour la première fois par Truesdell [141]
et étendue aux fractions volumiques comme décrit par Bowen [142, 143], de Boer [144], et
Ehlers et Bluhm [145]. Elle représente une extension de la mécanique classique du continuum
monophasique pour une approche de modélisation multiphasique, dans laquelle une large

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 72

Chapitre IV Etude numérique : Modèle 1D

sélection de modèles pour les matériaux granulaires a été développée, principalement dans le
contexte du matériau simple de Noll [146].

La théorie des milieux poreux, basée sur le continuum, offre la possibilité d'obtenir un
cadre de modélisation consistant, où toutes les équations de conservation classiques peuvent
être écrites et résolues mathématiquement. Dans ce cadre, on prend une caractéristique
moyenne des attributs du matériau du VER (volume élémentaire représentatif), où les
caractéristiques individuelles des éléments sont remplacées par cette valeur moyenne
représentative au lieu de considérer le comportement individuel des éléments dans ce volume.
Un exemple de cette méthode consiste à utiliser la vitesse moyenne, la porosité moyenne et la
résistance moyenne au lieu du comportement particulier de chaque particule [145].

IV.1 Processus d’érosion interne basé sur le continuum

Dans un cadre de modélisation basé sur la mécanique du continuum, les phénomènes
d'érosion peuvent être expliqués par la modélisation de l'échange de masse entre la phase liquide
et du squelette solide. Pour chaque continuum, des équations d'équilibre peuvent être écrites et
l'interaction de l'échange de masse entre les phases, par exemple la force motrice de l'érosion,
peut être incluse. Par conséquent, à mesure que l'érosion interne se poursuit, la partie continue
de la phase solide se réduit en masse et en volume au fil du temps et le continuum de la phase
liquide augmente en volume en remplissant l'espace créé (c'est-à-dire que le matériau devient
de plus en plus poreux). Cette approche continue peut-être décrite comme une modélisation
"macro-échelle". Pour les modèles ayant une perspective macro d'un système, l'accent est
principalement mis sur la production et le comportement global de ce système plutôt que sur le
comportement des éléments qui le composent. Une analogie avec ce concept est l'image sur un
écran de télévision ou un moniteur, qui est obtenue non pas en se concentrant sur chaque pixel
et en analysant son état et sa variation de couleur à chaque instant, mais plutôt par le résultat
combiné assemblé.

Vardoulakis et al. [85] ont été les premiers à proposer un cadre théorique basé sur un
continuum pour l'aspect hydromécanique des phénomènes d'érosion pour la production de
sable, qui est un phénomène d'érosion des parois des puits de forage dans les grès des réservoirs
d'hydrocarbures pendant l'extraction des fluides. Ils ont introduit une théorie de mélange à trois
composants pour un VER constitué d'un squelette solide, de particules fluidisées et du fluide.
Nous allons présenter en détail dans la section suivante, la formulation mathématique de ce

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 73

Chapitre IV Etude numérique : Modèle 1D

modèle avec une modification sur la loi d’érosion que nous avons proposé, une validation par
rapport aux résultats expérimentaux sera abordée.

IV.2 Modélisation 1D de l’érosion interne

Comme on avait déjà discuté au chapitre 2, de nombreux modèles numériques pour
étudier l’érosion interne ont été proposés [100, 147–151]. Plusieurs de ces modèles sont conçus
dans le cadre de la théorie des milieux poreux continus basés sur la méthode des éléments finis
ou la méthode des différences finis. La plupart de ces modèles décrivent efficacement le
détachement et le transport des particules fines à travers la matrice solide. Cette partie est
consacrée à la présentation d'une approche de modélisation numérique pour décrire le transport
des particules fines lors de la suffusion. Cette approche traite le matériau comme un milieu
continu. Le modèle est construit en supposant trois constituants. L'objectif d'un tel "problème
simple" est de voir l’évolution des principales caractéristiques physiques et mécaniques du sol
pendant le processus de suffusion.

IV.3 Formulation mathématique

On considère un élément de volume 𝑉 d’un milieu granulaire poreux saturé. A l’état
initial, le volume est composé de deux phases : une phase solide et une phase fluide avec des
volumes (𝑉𝑠 ), (𝑉𝑓 ) et masse (𝑀𝑠 ),(𝑀𝑓 ) respectivement. Lorsque le processus d’érosion
commence, une troisième phase apparait, cette phase représente les particules du sol érodables
avec un volume (𝑉𝑠𝑒 ) et une masse (𝑀𝑠𝑒 ) (Figure IV. 1). Le volume des vides 𝑑𝑉𝑣 est le volume
des pores interconnecté "qui est entièrement occupé par le mélange de particules fluides et fines-
érodées. Il est défini comme étant la somme du volume du fluide et le volume du solide érodé

𝑉𝑣 = 𝑉𝑓 + 𝑉𝑠𝑒 (24)

dVf dVf dVv

Après érosion dVse

dV
dVs
dVs

Figure IV. 1 Schéma représentatif de l’état initial et après érosion

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 74

Chapitre IV Etude numérique : Modèle 1D

La partie solide érodée représente des particules en suspension qui se déplacent avec le
fluide. De plus, il est supposé que le fluide et les particules érodées partagent à tout instant la
même vitesse. Le squelette solide est supposé être rigide. En conséquence, les vitesses des
différents composants sont

𝑓
𝑣𝑖𝑠𝑒 = 𝑣𝑖 = 𝑣̅𝑖 , 𝑣𝑖𝑠 = 0 (25)
La porosité du milieu est définie par n avec

𝑉𝑣 𝑉𝑓 + 𝑉𝑠𝑒
𝑛= = (26)
𝑉 𝑉
La concentration des particules solides érodées est définie par c.

𝑉𝑠𝑒
𝑐= (27)
𝑉𝑣
On rappelle que dans le problème d’érosion, les valeurs de 𝑐 et 𝑛 dépendent de la
position 𝑥𝑖 et du temps 𝑡.

Les densités de la phase fluide et solide-érodée sont les mêmes que les densités des
matériaux constituants initiaux c.-à-d. si le fluide est une eau et le solide est du quartz, 𝜌𝑓 =
1000⁡𝑘𝑔/𝑚3 et 𝜌𝑠𝑒 = 2650⁡𝑘𝑔/𝑚3 sont respectivement leurs densités.

𝑀𝑓 𝑀𝑠𝑒 𝑀𝑠
𝜌𝑓 = , 𝜌𝑠𝑒 = = = 𝜌𝑠 (28)
𝑉𝑓 𝑉𝑠𝑒 𝑉𝑠

La densité partielle du mélange (fluide-solide érodé) est définie par 𝜌̅ avec

𝑀𝑓 + 𝑀𝑠𝑒
𝜌̅ = (29)
𝑉𝑓 + 𝑉𝑠𝑒
En remplaçant les équations (27) et (28) dans l’équation (29), la formule de la densité
partielle du mélange devient comme suit

𝜌̅ = (1 − 𝑐)𝜌𝑓 + 𝑐𝜌𝑠 (30)

Les densités partielles des trois phases sont définies comme

𝑀𝑠𝑒 𝑉𝑠𝑒 𝑐𝑛𝑉

𝜌̅𝑠𝑒 = = 𝜌𝑠 = 𝜌𝑠 = 𝑐𝑛𝜌𝑠
𝑉 𝑉 𝑉

(31)
𝑀𝑠 𝑉𝑠 𝑉 − 𝑉𝑓 − 𝑉𝑠𝑒
𝜌̅𝑠 = = 𝜌𝑠 = 𝜌𝑠 = (1 − 𝑛)𝜌𝑠
𝑉 𝑉 𝑉

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 75

Chapitre IV Etude numérique : Modèle 1D

𝑀𝑓 𝑉𝑓 𝑛𝑉 − 𝑉𝑠𝑒
𝜌̅𝑓 = = 𝜌𝑓 = 𝜌𝑓 = (1 − 𝑐)𝑛𝜌𝑓
𝑉 𝑉 𝑉
La densité totale du mélange est donnée par la somme des densité partielles de chaque
phase

𝑀 𝑀𝑠 + 𝑀𝑓 + 𝑀𝑠𝑒
𝜌= = = 𝜌̅𝑠 + 𝜌̅𝑓 + 𝜌̅𝑠𝑒
𝑉 𝑉 (32)
= (1 − 𝑛)𝜌𝑠 + 𝑛𝜌𝑓 + 𝑐𝑛(𝜌𝑠 − 𝜌𝑓 )
Le débit volumique 𝑞̅𝑖 (le volume d’écoulement à travers la section transversale unitaire
à l’instant dt) et la vitesse 𝑣̅𝑖 du mélange sont relié par la relation suivante [85]

𝑉̅ 𝑉̅ 𝑉̅ 𝑞̅𝑖
𝑞̅𝑖 = , 𝑣̅𝑖 = = = (33)
𝑆𝑖 𝑑𝑡 𝑆𝑖̅ 𝑑𝑡 𝑛𝑆𝑖 𝑑𝑡 𝑛
Où 𝑉̅ est le volume de l’écoulement à travers la section transversale 𝑆𝑖 à l’instant dt. 𝑆𝑖̅
est une partie de 𝑆𝑖 , qui correspond à la surface des pores.

IV.4 Equations d’équilibre de masse

Considérant maintenant le milieu bi-phasique qui consiste en solide (indice 1) et
mélange fluide-solide érodé (indice 2). La densité partielle de ces deux phases sont

𝑀𝑠 𝑀𝑓 +𝑀𝑠𝑒
𝜌1 = = (1 − 𝑛)𝜌𝑠 , 𝜌2 = = 𝑛𝜌̅ (34)
𝑉 𝑉

Les vitesses des deux phases sont 𝑣𝑖1 = 𝑣𝑖𝑠 = 0 et 𝑣𝑖2 = 𝑣̅𝑖 . L’équilibre de masse de la
phase solide-érodé est donné par l’expression suivante

𝜕𝜌̅𝑠𝑒 𝜕 𝜕𝜌̅𝑠𝑒
+ 𝑑𝑖𝑣(𝜌̅𝑠𝑒 𝑣𝑖𝑠𝑒 ) − (𝐷𝑖𝑗 ) = 𝑚̇ (35)
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑖
Où 𝐷𝑖𝑗 (𝑚/𝑠 2 ) est le coefficient de dispersion mécanique. Ce paramètre sera négligé dans le
cas du processus d’érosion unidimensionnelle (1D) traité dans cette partie. Cependant, nous allons le
prendre en compte dans l’étude de l’érosion dans les digues (processus 2D). Le paramètre 𝑚̇ (𝑘𝑔/𝑚3 . 𝑠)
est le terme de génération de masse, qui correspond à la masse érodée et transportée à chaque instant et
à chaque point. 𝑚̇ consiste en général en deux termes

𝑚̇ = 𝑚̇𝑒𝑟 − 𝑚̇𝑑𝑒𝑝 (36)

Avec 𝑚̇ 𝑒𝑟 est la masse érodée et 𝑚̇ 𝑑𝑒𝑝 est la masse déposée

L’équation du bilan de masse de la phase solide exprimée en termes de densité partielle

et de la vitesse de cette phase est donnée par

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 76

Chapitre IV Etude numérique : Modèle 1D

𝜕𝜌̅𝑠
+ 𝑑𝑖𝑣(𝜌̅𝑠 𝑣𝑖𝑠 ) = 𝑚̇𝑠 (37)
𝜕𝑡
En utilisant l’équation (31), on obtient

𝜕𝑛 𝑚̇𝑠
=− (38)
𝜕𝑡 𝜌𝑠
L’équation de l’équilibre de masse de la phase solide-érodé s’écrit comme suit

𝜕𝜌̅𝑠𝑒 (39)
+ 𝑑𝑖𝑣(𝜌̅ 𝑒𝑠 𝑣𝑖𝑠𝑒 ) = 𝑚̇𝑠𝑒
𝜕𝑡
Compte tenu de la définition (31) de la densité partielle 𝜌̅𝑠𝑒 , l’équation de l’équilibre de masse
de la phase solide-érodé devient

𝜕(𝑐𝑛) 𝑚̇𝑠𝑒
+ 𝑑𝑖𝑣(𝑐𝑞̅𝑖 ) = (40)
𝜕𝑡 𝜌𝑠
Sachant que la masse de toutes les particules érodées provient de la phase solide, c-à-d.
𝑚̇𝑠𝑒 = −𝑚̇𝑠 et en éliminant le terme de génération de masse en utilisant l’équation (38), on
obtient l’équation reliant le changement de la porosité au flux des particules érodées

𝜕𝑛 𝜕(𝑐𝑛)
= + 𝑑𝑖𝑣(𝑐𝑞̅𝑖 ) (41)
𝜕𝑡 𝜕𝑡

D’autre part, l’équation d’équilibre de masse de la phase fluide est donnée par

𝜕𝜌̅𝑓 𝑓 (42)
+ 𝑑𝑖𝑣(𝜌̅𝑓 𝑣𝑖 ) = 0
𝜕𝑡
Et à travers les équations (22) et (16), elle devient

𝜕(1 − 𝑐)𝑛 (43)

+ 𝑑𝑖𝑣((1 − 𝑐)𝑞̅𝑖 ) = 0
𝜕𝑡
Finalement, en combinant (41) et (43), on obtient l’équation de continuité de débit

𝑑𝑖𝑣(𝑞̅𝑖 ) = 0 (44)
Les équations (38), (41) et (44) constituent l'ensemble des équations d'équilibre de
masse pour le problème considéré. Ces équations ne suffisent pas pour résoudre le problème de
l'érosion car il faut une relation constitutive pour le terme de génération de masse 𝑚̇⁄𝜌𝑠 .

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 77

Chapitre IV Etude numérique : Modèle 1D

Une loi constitutive de la quantité de masse érodée, suggérant que l'érosion est
provoquée par le flux des particules érodées "⁡𝑐𝑞̅𝑖 ", est donnée par la relation suivante [85, 152,
153]

𝑚̇ ′
𝑐2
= 𝜆 (𝑐 − ) |𝑞̅𝑖 | (45)
𝜌𝑠 𝑐𝑐𝑟

Vardoulakis et al. [85] ont supposé que 𝜆′ prend la forme suivante (équation 46), seulement
pour le problème de production de sable.

𝜆′ = 𝜆(1 − 𝑛) (46)
Cette relation entraine l’augmentation de la porosité du milieu jusqu'à 1 (érosion totale
du milieu), elle est utilisée par Vardoulakis et al. [85] pour résoudre le problème de production
de sable. Pour adapter cette relation au problème de suffusion, nous avons introduit un
paramètre 𝑛𝑐𝑟 , qui limite la diffusion de la porosité dans le milieu. Cette limitation est justifiée
par le fait que la suffusion entraine seulement les particules fines du sol et conduit à une
variation limitée de la porosité du milieu. En effet, l’équation (46) est remplacée par l’équation
suivante (équation 47)

𝑛 𝑛𝑐𝑟
𝜆′ = 𝜆 (1 − )=𝜆 (47)
𝑛𝑐𝑟 𝑛𝑐𝑟 − 𝑛
En introduisant l’équation (47) dans (45) donne

𝑚̇ 𝑛 𝑐2
= 𝜆 (1 − ) (𝑐 − ) |𝑞̅𝑖 | (48)
𝜌𝑠 𝑛𝑐𝑟 𝑐𝑐𝑟
Où 𝑐𝑐𝑟 est la valeur critique de 𝑐. 𝑛𝑐𝑟 est la valeur critique de la porosité (n). Le coefficient 𝜆 a
une dimension inverse de la longueur, [𝜆]=L-1, il est lié probablement à la fréquence spatiale
des points potentiels pour l’initiation d’érosion dans le squelette solide du milieu poreux et doit
être déterminé expérimentalement.

En combinant l’équation (38) avec les équations (41) et (48) et en tenant compte de
l’équation de continuité (44), on obtient la formulation de la cinétique d’érosion suivante

𝜕(𝑐𝑛) 𝜕𝑐 𝜕𝑛
= 𝑞̅ + (49)
𝜕𝑡 𝜕𝑥 𝜕𝑡
𝜕𝑛 𝑛 𝑐2
= 𝜆 (1 − ) (𝑐 − ) |𝑞̅ | (50)
𝜕𝑡 𝑛𝑐𝑟 𝑐𝑐𝑟

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 78

Chapitre IV Etude numérique : Modèle 1D

Par ailleurs, l’écoulement dans le milieu poreux est gouverné par la loi de Darcy qui
relie la vitesse de décharge au gradient de la pression, à la viscosité et à la structure du milieu
poreux :

𝜅
𝑢𝑑𝑙 = − 𝜇 ∇𝑝 (En négligeant l’effet de la gravité)
𝑘
(51)
𝑢𝑑𝑙 = − 𝜇 (∇𝑝 + 𝜌̅ 𝑔) (Avec prise en compte de l’effet de la gravité)

Avec, 𝜇 est viscosité dynamique du mélange de fluide et particules érodées donnée par
[154]

𝜇 = 𝜇0 (1 + 2.5𝑐) (52)
µ0 étant la viscosité initiale qui correspond à la viscosité dynamique du fluide seul.

𝜅 : est la perméabilité intrinsèque du milieu poreux. Comme la variation de la porosité

affecte considérablement la perméabilité, selon l’équation de Carman-Kozeny, elle est donnée
par :

𝑛3
𝜅 = 𝜅0 (53)
(1 − 𝑛)2
𝜅𝜌𝑔
𝐾= (54)
𝜇
Avec 𝜅0 est la perméabilité intrinsèque du matériau.

𝜌̅ est donné par l’équation (30)

𝜌̅ = (1 − 𝑐)𝜌𝑓 + 𝑐𝜌𝑠

Pour plus de commodité, Gravanis et al. [155] suggèrent que 𝑞̅ peut être remplacé dans
les équations (49) et (50) par la vitesse de Darcy (𝑢𝑑𝑙 ). On obtient donc

𝜕(𝑐𝑛) 𝜕𝑐 𝜕𝑛
= 𝑢𝑑𝑙 + (55)
𝜕𝑡 𝜕𝑥 𝜕𝑡
𝜕𝑛 𝑛 𝑐2 (56)
= 𝜆 (1 − ) (𝑐 − ) |𝑢𝑑𝑙 |
𝜕𝑡 𝑛𝑐𝑟 𝑐𝑐𝑟

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 79

Chapitre IV Etude numérique : Modèle 1D

IV.5 Dégradation de la cohésion

Afin d’évaluer l’effet de l’érosion interne sur la cohésion des sols, nous avons considéré
dans cette étude une loi linéaire qui relie la porosité du milieu à sa cohésion. Cette loi est utilisée
ailleurs par Gravanis et al. [155]. Cette relation considère que la cohésion du milieu se dégrade
au fur et à mesure que la porosité augmente. Elle est donnée par la relation suivante :

1 − 𝑛(𝑡)
𝐶̅ (𝑡) = 𝐶0 (57)
1 − 𝑛0
Avec

𝐶̅ (𝑡) est la cohésion du sol après érosion ;6363

𝐶0 est la cohésion initiale du sol ;

n(t) et n0 sont la porosité actuelle et initiale respectivement.

IV.6 Paramètres d’entrée

Les paramètres d’entrée du modèle sont montrés au Tableau IV. 1. La porosité initiale
est calculée en considérant un sol ayant une densité sèche de 1.6 et une densité absolue de 2.65.
la viscosité initiale est celle de l’eau, et la concentration initiale est prise si petite juste pour
déclencher le processus d’érosion. Les valeurs du coefficient d’érosion, la longueur de
l’échantillon, la charge hydraulique ainsi que la porosité critique seront discutées dans la section
de l’étude paramétrique.

Tableau IV. 1 Paramètres d’entrée du modèle

Paramètres Valeur Description

n0 0.396 Porosité initiale
c0 10-3 Concentration initiale
ccr 0.3 Concentration critique
k0 2.7436e-13m2 Perméabilité intrinsèque initiale
µ0 10-3 Pa.s Viscosité initiale
ρf0 1000 kg/m3 Masse volumique initiale de l’eau
ρs 2650 kg/m3 Masse volumique solide initiale

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 80

 Chapitre IV Etude numérique : Modèle 1D

IV.7 Conditions aux limites

Les équations décrivant la cinétique d’érosion ont été résolu en considérant la condition
à l’entrée de l’échantillon suivante :

✓ Condition de type Dirichlet, [c (0, t) = c0 = cste]. Cette condition, selon

Vardoulakis conduit à l’effet de l’érosion de surface [85]. Dans notre cas, en
introduisant le coefficient ncr, les simulations montrent une érosion interne
régressive.
✓ Une charge hydraulique (H) en amont et une charge nulle en aval.

V. Résultats et discussions

Nous allons dans ce qui suit présenter les résultats des simulations numérique issus du
modèle mathématique présenté précédemment. Pour résoudre le système d’équations non-
linéaire de la cinétique d’érosion, nous avons eu recours au logiciel Comsol multiphysics. Le
logiciel Comsol est un code de calcul par éléments finis qui offre la possibilité de coupler
plusieurs physiques au même temps et permet d’effectuer des études paramétriques sur les
différents paramètres.

V.1 Etude paramétrique

Afin de montrer la capacité du modèle à simuler le comportement érosif, nous avons
effectué plusieurs simulations en variant les différents paramètres du modèle, à savoir, le
coefficient d’érosion, la charge hydraulique, la longueur de l’échantillon et la porosité critique
qui est définie dans notre étude comme le seuil de l’érosion interne.

Effet du coefficient d’érosion

Le coefficient d’érosion (λ) est probablement lié à l’érodabilité interne du sol qui dépend
essentiellement des éléments structuraux du milieu poreux (granulométrie du sol, densité du
sol, plasticité du sol... etc.) Ce paramètre est utilisé largement dans les modèles d’érosion
proposé par plusieurs auteurs [83–85, 155, 156].Pour montrer l’effet du coefficient d’érosion,
nous avons effectué des simulations avec trois valeurs (500, 1000 et 1500 1/m). Les valeurs de
la charge hydraulique, la porosité critique et la longueur de l’échantillon sont prises égale à
H=1m, ncr=0.5 et L=5cm respectivement. Les valeurs des autres paramètres sont montrées au
Tableau IV. 1.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 81

Chapitre IV Etude numérique : Modèle 1D

La Figure IV. 2 montre la variation temporelle de la concentration, la porosité et de la

cohésion au point de sortie de l’écoulement pour les différentes valeurs de λ. La Figure IV. 2a
montre que le coefficient d’érosion amplifie le processus d’érosion et conduit à un entrainement
important de la concentration, d’où la porosité du milieu augmente aussi jusqu’à atteindre sa
valeur maximale qui est définie comme seuil d’érosion (Figure IV. 2). On peut constater aussi
que, pour λ=1000 et λ=1500[1/m], la valeur de la porosité critique est atteinte dans les deux
cas, cela veut dire que les deux valeurs ramènent le processus d’érosion jusqu’à sa fin, et la
seule différence est que la valeur de λ=1500 entraine plus rapidement les particules fines et
accélère le phénomène par rapport à la valeur de λ=1000. La valeur de λ=500 étant faible pour
entrainer une quantité importante de fines, en effet, la porosité critique ne peut pas être atteinte.
La Figure IV. 2c illustre l’effet du coefficient d’érosion sur la variation temporelle de la cohésion
au milieu de l’échantillon. Etant donné que la cohésion est supposée être inversement
proportionnelle à la porosité, elle présente un comportement similaire à celui de la porosité.

(a)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 82

Chapitre IV Etude numérique : Modèle 1D

(b)

(c)

Figure IV. 2 Effet du coefficient d’érosion sur les différents paramètres : (a) Concentration, (b)
Porosité, (c) Cohésion.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 83

Chapitre IV Etude numérique : Modèle 1D

Effet de la charge hydraulique

La Figure IV. 3 montre l’effet de la charge hydraulique sur la porosité, la concentration
et la cohésion. Dans cette simulation, nous avons pris λ=1500 [1/m], ncr=0.5 et L=5cm. Afin
de montrer l’effet de la charge hydraulique, trois valeurs ont été choisi : 0.5m, 1m et 2m. On
remarque sur la Figure IV. 3a que la valeur de la concentration maximale est atteinte dans les
trois cas de la charge hydraulique mais à des instants différents. Plus la charge est grande, plus
la concentration maximale est atteinte plus rapidement. Concernant la porosité, on remarque
aussi qu’elle a atteint la valeur critique choisie comme seuil d’érosion pour les différentes
valeurs de (H), cependant, la valeur critique est atteinte plus rapidement en augmentant la
charge. Ce constat nous amène à dire que la charge hydraulique tend à accélérer le processus.

(a)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 84

Chapitre IV Etude numérique : Modèle 1D

(b)

(c)

Figure IV. 3 Effet de la charge hydraulique sur les différents paramètres : (a) Concentration,
(b) Porosité, (c) Cohésion.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 85

Chapitre IV Etude numérique : Modèle 1D

Effet de longueur de l’échantillon

L’effet de la longueur de l’échantillon est étudié en choisissant trois longueur : 5, 10 et
15cm. Pour cette étude, la charge hydraulique est prise égale à H=2m, le coefficient d’érosion
λ=1000 [1/m] et une porosité critique égale à ncr=0.5. Les résultats sont montrés dans la Figure
IV. 4. On remarque que plus la dimension est grande, la concentration est plus importante. On
remarque aussi que le processus se fait plus rapidement lorsque la longueur est plus petite. On
voit dans la Figure IV. 4a que l’échantillon avec 5cm atteint la fin du processus alors que pour
l’échantillon de 15cm, le processus vient de commencer. Cela s’explique aussi en regardant la
courbe de la variation temporelle de la porosité (Figure IV. 4b).

(a)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 86

Chapitre IV Etude numérique : Modèle 1D

(b)

(c)

Figure IV. 4 Effet de la longueur sur les différents paramètres : (a) Concentration, (b) Porosité,
(c) Cohésion.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 87

Chapitre IV Etude numérique : Modèle 1D

Effet de la porosité critique ncr

L’effet de la porosité critique (ncr) est étudié en prenant trois différentes valeurs :
0.45,0.5 et 0.55.: les valeurs de la charge hydraulique, de la longueur et du coefficient d’érosion
sont de : 2m, 5cm et 1500 [1/m] respectivement. Les résultats obtenus pour la concentration et
la porosité ainsi que de la cohésion sont montrés dans la Figure IV. 5. On remarque que, plus la
valeur de ncr est grande, on a plus de concentration de fines tel que montré dans la Figure IV. 5a.
Le paramètre ncr étant choisi comme le seuil d’érosion, peut refléter en réalité l’érodabilité d’un
sol donné. Concernant la porosité, elle augmente dans tous les cas jusqu’à atteindre sa valeur
critique. En effet, la porosité critique doit être déterminée en fonction de l’érodabilité du
matériau, c-à-d. plus le sol est érodable, plus la porosité critique devient importante. Dans la
section suivante, nous allons déterminer pour chaque matériau étudié expérimentalement sa
porosité critique en tenant compte de la variation de la cohésion.

(a)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 88

Chapitre IV Etude numérique : Modèle 1D

(b)

(c)

Figure IV. 5 Effet de la porosité critique sur les différents paramètres : (a) Concentration, (b)
Porosité, (c) Cohésion.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 89

Chapitre IV Etude numérique : Modèle 1D

V.2 Confrontation des résultats expérimentaux au modèle numérique

Les essais d’érosion présentés dans le chapitre précédent ont été simulés numériquement
afin de pouvoir examiner la performance du modèle. La série de matériaux utilisée dans les
expérimentations a été implanté dans le modèle en choisissant les mêmes paramètres, à savoir :
la longueur de l’échantillon, la porosité initiale, la viscosité de l’eau et la pression appliquée.
Le coefficient d’érosion λ a été déterminé de manière à pouvoir reproduire les valeurs de la
cohésion après érosion obtenues expérimentalement. De même pour le paramètre ncr, qui
représente la valeur critique que peut atteindre le sol après érosion.

V.3 Géométrie 1D pour l’analyse de l’érosion interne

Vu que les essais expérimentaux ont été effectués sur des échantillons de petite longueur
(4cm), nous avons considéré que l’écoulement est unidirectionnel, la géométrie considérée
consiste donc en un segment de droite où l’écoulement se fait de manière descendante. La Figure
IV. 6 montre le schéma explicatif du modèle ainsi que les conditions de pression aux limites de
l’échantillon.

Figure IV. 6 Schéma explicatif du modèle 1D

V.4 Simulation numérique

Le système d'équations pour la cinétique de l'érosion interne a été résolu par la méthode
des éléments finis en utilisant le logiciel Comsol Multiphysics [102]. Comsol est un outil de
simulation multiphasique permettant de résoudre, de manière couplée, des équations

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 90

Chapitre IV Etude numérique : Modèle 1D

différentielles partielles par la méthode des éléments finis. Il fournit une base de données
d'équations permettant de modéliser différents phénomènes physiques. Les modules utilisés
dans Comsol pour notre étude sont les suivants : La loi de Darcy (dl) pour l'écoulement de l'eau
dans le milieu poreux, l'équation aux dérivées partielles (PDE) utilisant la porosité comme
variable dépendante et une seconde (PDE) pour la concentration en solides fluidisés.

V.5 Paramètres d’entrée du modèle

Pour résoudre le processus d'érosion interne dans la boîte de cisaillement, une
modélisation Comsol 1D a été réalisée. L'échantillon est modélisé par un segment de ligne
droite de longueur "L" qui représente la longueur de l'échantillon. Le Tableau IV. 2 montre les
paramètres d'entrée du modèle.

Tableau IV. 2 Paramètres d’entrée du modèle d’érosion

Paramètres Valeur Description

L 0.04 m Longueur de l’échantillon
n0 0.396 Porosité initiale
ncr Tableau IV. 3 Porosité critique
c0 10-3 Concentration initiale
ccr 0.3 Concentration critique
H 0.4 m Charge hydraulique
k0 2.7436e-13m2 Perméabilité initiale
λ Tableau IV. 3 Coefficient d’érosion
µ0 10-3 Pa.s Viscosité initiale
ρf0 1000 kg/m3 Masse volumique initiale de l’eau
ρs 2650 kg/m3 Masse volumique solide initiale

V.6 Calcul de la porosité initiale et critique des différents échantillons

Porosité initiale
Sachant que la masse volumique sèche de tous les échantillons est prise égale à 1.6
g/cm3, et que la masse volumique solide du sable et des autres fines utilisées (Kaolinite, Limon)
est presque de même grandeur et d’environ 2.65 g/cm3. Nous avons considéré une même
porosité initiale qui peut être calculée par la relation suivante :

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 91

Chapitre IV Etude numérique : Modèle 1D

𝛾𝑑
𝑛=1−
𝛾𝑠

1.6
On obtient donc 𝑛0 = 1 − 2.65 = 0.396.

Avec 𝑛0 est la porosité initiale considérée dans les simulations numériques.

Porosité critique
La porosité critique est calculée en considérant la relation entre la dégradation de la
cohésion et la porosité de chaque échantillon. Nous rappelons que nous avons considéré une loi
qui stipule que la cohésion du milieu se dégrade au fur et à mesure que la porosité augmente
(eq.42). Afin de trouver la porosité critique pour chaque échantillon, nous allons considérer sa
cohésion avant et après érosion obtenue expérimentalement.

On a donc

𝐶̅ (58)
𝑛𝑐𝑟 = 1 − (1 − 𝑛0 )
𝐶0
Le Tableau IV. 3 montre les valeurs de λ et ncr à considérer dans l’analyse numérique. La
porosité initiale est considérée constante tandis que la porosité critique est calculée par la
formule (48).

Tableau IV. 3 valeurs de λ et ncr pour les différentes simulations

Cohésion avant Cohésion après érosion (kPa)

Spécimen n0 λ ncr
érosion (kPa) Expérimentale Numérique
10 % (K) 0.396 4500 0.423 17.9 17.1 17.1
15 % (K) 0.396 3000 0.451 19.7 17.9 17.9
20% (K) 0.396 2500 0.465 20.9 18.5 18.5
10 % (SI) 0.396 2700 0.461 18.5 16.5 16.5
15 % (SI) 0.396 2500 0.465 19.1 16.9 16.9
20 % (SI) 0.396 2200 0.488 20.3 17.2 17.2
10 % (SA) 0.396 4500 0.422 16.1 15.4 15.4
15 % (SA) 0.396 3000 0.448 17.3 15.8 15.8
20 % (SA) 0.396 2700 0.456 17.9 16.1 16.1

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 92

 Chapitre IV Etude numérique : Modèle 1D

VI. Résultats numériques

Nous allons présenter dans cette section les différents résultats numériques obtenus
après simulation de chaque matériau. Dans un premier temps, nous allons comparer l’effet de
la teneur en fines sur les résultats obtenus pour les matériaux de même nature. Ensuite, nous
allons comparer entre les matériaux de nature différente.

VI.1 Résultats de simulation

Comme indiqué dans le paragraphe précédent, nous allons présenter ici les résultats
obtenus pour les simulations des différents matériaux étudiés. Les figures (Figure IV. 7, Figure
IV. 8, Figure IV. 9, Figure IV. 10, Figure IV. 11, Figure IV. 12, Figure IV. 13, Figure IV. 14, Figure
IV. 15) montrent les résultats pour les différents paramètres du modèle numérique. L’effet de la
teneur et du type des fines sur ces paramètres est représenté dans les mêmes figures. Les
résultats présentés concernent essentiellement la variation temporelle des différents paramètres
du modèle au milieu de l’échantillon (x=L/2) (cette zone nous intéresse parce que c’est l’endroit
où le cisaillement s’effectue dans l’essai de cisaillement directe) ainsi que leurs distributions
spatiales à la fin du processus d’érosion interne.

La Figure IV. 7 représente le profil de pression à la fin du processus d’érosion pour les
différentes simulations, on remarque que la courbe change d’allure après l’érosion interne, la
diminution de la pression interstitielle est due probablement à l’augmentation de la perméabilité
du milieu. Il faut noter que le phénomène de dépôt est négligé dans cette étude, les surpressions
interstitielles sont donc absente. En effet, la pression dans l’échantillon tend à diminuer
beaucoup plus pour l’échantillon qui contient plus de fines. Ces échantillons présentaient un
potentiel d’érosion important.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 93

Chapitre IV Etude numérique : Modèle 1D

Figure IV. 7 Profil de la pression à la fin de l’érosion pour les différents matériaux

La Figure IV. 8 montre la variation temporelle de la porosité au milieu des échantillons.

On remarque que la porosité augmente jusqu’à atteindre la valeur critique définie comme le
seuil d’érosion. Cette valeur est choisie de manière à obtenir les mêmes résultats expérimentaux
concernant la dégradation de la cohésion. Il faut signaler que plus on s’approche de la sortie de
l’échantillon, plus le processus d’érosion se fait d’une manière plus rapide. Cela se traduit par
une érosion interne régressive qui est définie par la propagation du processus d’érosion de l’aval
vers l’amont de l’ouvrage dans les cas réels. On remarque aussi que l’échantillon présentant
plus de teneur en Limon est le plus susceptible à l’érosion. Ce constat est déjà observé
expérimentalement tel qu’il est présenté dans le Tableau III. 4 où les valeurs de ∆𝑅𝐶 pour le
mélange de limon sont plus importantes.

La Figure IV. 9 montre le profil linéaire de la porosité à la fin du processus de l’érosion

pour les matériaux étudiés. On voit bien que la porosité est importante en s’approchant de l’aval
de l’échantillon, ceci peut se traduire par une érosion plus importante à cet endroit. La porosité
diminue en progressant vers l’amont de l’échantillon où nous avons considéré que la porosité
reste inchangée vu qu’on a protégé la surface avec une couche de géotextile. Ce phénomène est
déjà défini dans le chapitre 1 comme l’érosion régressive.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 94

Chapitre IV Etude numérique : Modèle 1D

La Figure IV. 10 donne l’évolution temporelle de la concentration de particules dans le

fluide. Les courbes de la concentration présentent la même forme dans tous les cas. Elles sont
caractérisées par trois phases principales :

1. La phase 1 [0-1400 s] : Elle représente la zone d’érosion. La variation brusque

de la concentration des particules fines dans le fluide jusqu’au pic provient du
détachement de fines particules de sol sous l’effet des contraintes hydrauliques
résultantes.
2. La phase 2 [1400-2500s] : Elle représente une zone de chute de la concentration
après avoir atteint son pic. Le taux d’érosion tend à diminuer rapidement jusqu’à
des valeurs très faibles.
3. La phase 3 [2500-3600s] : Elle représente la zone dans laquelle la concentration
solide se stabilise au palier atteint à la fin de la deuxième phase. Une quantité
minimale de particules fines continue à être érodée avant de tendre vers une
valeur insignifiante de la concentration.

Ces courbes caractéristiques de l’érosion interne décrivant la cinétique de l’érosion

interne sont en bon accord avec d’autres études numériques et expérimentales [85, 157]. Ce
résultat montre la fiabilité du modèle à simuler le problème de l’érosion interne.

Le pic d’érosion atteint pour chaque matériau mets en évidence la susceptibilité à

l’érosion de chacun d’eux. L’échantillon de 20% de limon montre un pic plus important qui est
d’environ 0.12. Les échantillons de kaolinite et de sable fin (sable#2) montrent un
comportement très similaire.

D’autres chercheurs [85, 131, 157] indique que le pic d’érosion est d’autant plus
important si la charge hydraulique est importante. Dans notre cas les simulations ont été
effectuées avec une charge hydraulique constante. En effet, la différence entre les pics est
provoquée par le coefficient d’érosion « λ ». Ce coefficient lui-même dépend du gradient
appliqué et de la plasticité du matériau [155].

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 95

Chapitre IV Etude numérique : Modèle 1D

Figure IV. 8 Variation temporelle de la porosité au milieu des différents échantillons

La conductivité hydraulique peut être utilisée comme un indicateur dans l'étude du

transport des particules à travers le milieu poreux et de la sensibilité à l'érosion interne. Le
comportement suffusif d’un sol dépend de la distribution granulaire et de la teneur en fines. La
Figure IV. 11 montre la distribution spatiale de la perméabilité le long de l’échantillon. Cette
distribution de la perméabilité, calculée avec l'équation de Carman-Kozeny (53), reflète
l'évolution significative de l’hétérogénéité de la perméabilité suite à l'érosion interne. Cette
inhomogénéité de la perméabilité se traduit par une distribution de la pression interstitielle avec
des gradients de pression significativement plus faibles dans la région érodée.

La Figure IV. 12 montre l’évolution temporelle de la vitesse de Darcy au milieu de

l’échantillon. On remarque que la vitesse augmente au fur et à mesure que les particules fines
s’érodent ce qui laisse plus de vides et permet au fluide de s’écouler plus vite.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 96

Chapitre IV Etude numérique : Modèle 1D

Figure IV. 9 Profil linéaire de la porosité à la fin de l’érosion pour les différents matériaux

Figure IV. 10 Variation temporelle de la concentration au milieu de l’échantillon pour les

différents matériaux

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 97

Chapitre IV Etude numérique : Modèle 1D

Figure IV. 11 Profil linéaire de la perméabilité à la fin de l’érosion pour les différents matériaux

Figure IV. 12 Variation temporelle de la vitesse de Darcy au milieu de l’échantillon pour les
différents matériaux

Le résultat de simulation de la dégradation de la cohésion au milieu de l’échantillon est

montré sur la Figure IV. 13. On peut observer que les valeurs finales prédites correspondent

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 98

Chapitre IV Etude numérique : Modèle 1D

parfaitement aux valeurs obtenues expérimentalement par l’essai d’érosion-cisaillement. Ces

valeurs sont obtenues en fixant un coefficient d’érosion (λ) pour chaque test.

Figure IV. 13 Variation temporelle de la cohésion au milieu de l’échantillon pour les différents
matériaux

Figure IV. 14 Variation temporelle de la masse volumique du mélange au milieu de

l’échantillon pour les différents matériaux

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 99

 Chapitre IV Etude numérique : Modèle 1D

Figure IV. 15Variation temporelle de la masse viscosité de l’eau au milieu de l’échantillon pour
les différents matériaux

Vu que la masse volumique du mélange (eau-solide) et la viscosité du fluide sont

directement liées à la concentration des particules érodées, on voit bien dans les figures (Figure
IV. 14, Figure IV. 15) que leurs courbes présentent la même allure que celle de la concentration.
La masse volumique du mélange augmente jusqu’à une valeur de 1150 kg/m^3 puis revient à
sa valeur initiale (1000 kg/m3) à la fin du processus d’érosion.

VII. Conclusion

Dans ce chapitre, nous avons présenté un modèle hydromécanique de l’érosion interne.

Le modèle est basé sur l’équilibre de masse de la masse solide érodée et l’écoulement de fluide
interstitiel, et sur l’écoulement de Darcy en milieux poreux. L’étude, présentée ici, est limitée
aux problèmes d’écoulement unidirectionnel (1-D) où la dispersion hydrodynamique a été
négligée. Pour pouvoir représenter l’érosion interne par suffusion qui entraine une
augmentation limitée de la porosité du milieu, une loi d’érosion a été proposée avec un
paramètre (ncr) qui tient compte de cet effet. L’étude paramétrique effectuée sur les différents
paramètres a permis de montrer que le modèle peut reproduire le mécanisme d’érosion à travers
l’évaluation temporelle non-linéaire des différents paramètres tels que : la porosité (n), la

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 100

Chapitre IV Etude numérique : Modèle 1D

perméabilité (k), la concentration (c) et d’autres paramètres qui ont été présentés. L’influence
de chaque paramètre a été présentée et discutée. Le modèle a permis d’identifier les zones les
plus érodables qui se concentrent principalement en aval (zone de fortes vitesses). La
confrontation des résultats expérimentaux présentés au chapitre 3 au modèle numérique a été
présentée et a permis de déterminer le coefficient d’érosion pour chaque matériau.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 101

 Chapitre V

Etude numérique :
Modèle 2D
 Chapitre V Etude numérique : Modèle 2D

I. Introduction

Dans le chapitre précédent, le modèle à trois composantes de l'érosion interne a été

amélioré en ajoutant le terme 𝑛𝑐𝑟 qui permet de définir un seuil pour le processus d’érosion
interne par suffusion. Une plate-forme de simulation a ensuite été mise en œuvre sur le logiciel
de simulation par éléments finis Comsol Multiphysics pour simuler les essais d'érosion
descendante sous charge hydraulique constante en mettant l'accent sur l'évolution de la
perméabilité, la porosité et la concentration des particules érodées pendant et après l'érosion.
Dans ce chapitre, un modèle 2D plus général qui introduit l’effet de la dispersion a été proposé
afin de réaliser des analyses à l’échelle d'un ouvrage entier. Le modèle a été aussi implémenté
dans le code de calcul par éléments finis Comsol Multiphysics afin d'évaluer comment l'érosion
interne peut avoir un impact sur la stabilité des structures en terre. Après avoir simulé quelques
cas numériques, nous avons ensuite considéré que la cohésion et l’angle de frottement interne
du sol constituant l’ouvrage se dégradent au fur et à mesure que la porosité du milieu augmente.
Les distributions spatiales de ces deux paramètres vont être implémentées ensuite dans le code
de calcul Geoslope afin d’évaluer l’évolution de la stabilité de l’ouvrage dans le temps.

II. Formulation du modèle 2D

Les équations qui gouvernent le processus d’érosion en 2D sont données comme suit tel
que décrit dans le chapitre 4. Ici, nous avons ajouté le terme de dispersion qui est associé, à
l’échelle microscopique, à la fois au mouvement du fluide et à la structure porale complexe du
milieu poreux, ce qui génère l’hétérogénéité des vitesses de l’écoulement. Il existe trois causes
de ce phénomène [158] :

✓ L’écoulement traverse les grands pores plus vite qu’à travers les petits ;
✓ L’écoulement peut parcourir des chemins significativement différents
correspondant à différents chenaux à cause de la tortuosité du milieu ;
✓ Lorsque l’écoulement traverse les pores, il va se déplacer plus rapidement au
centre du pore que le long des parois à cause de l’effet des forces de frottement.

Les équations de gouvernance du processus d’érosion sont données donc par :

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 102

Chapitre V Etude numérique : Modèle 2D
𝜕𝑐 𝜕𝑛
𝑛 + 𝑢𝑑𝑙 . ∇𝑐 − ∇. 𝐷. ∇(𝑐. 𝑛) = (1 − 𝑐) (59)
𝜕𝑡 𝜕𝑡
𝜕𝑛 𝑛 𝑐2
= 𝜆 (1 − ) (𝑐 − ) 𝑢𝑑𝑙 (60)
𝜕𝑡 𝑛𝑐𝑟 𝑐𝑐𝑟
Où 𝑛 est la porosité, 𝑐 est la concentration, 𝐷 est le coefficient de dispersion, 𝜆 est le
coefficient d’érosion, 𝑛𝑐𝑟 est la porosité critique et 𝑐𝑐𝑟 est la concentration critique. 𝑢𝑑𝑙 est la
vitesse de Darcy donnée comme suit

𝑘 (61)
𝑢𝑑𝑙 = − ∇𝑝
𝜇
Dans cette équation, 𝑢𝑑𝑙 est le champ de vitesse de Darcy ou le vecteur de décharge
spécifique exprimé en (m/s), k est la perméabilité intrinsèque du milieux poreux exprimé en
(m2), 𝜇 est la viscosité dynamique du fluide exprimé en (Pa.s), 𝑝 est la pression interstitielle
exprimé en (Pa).

Pour les applications à grande échelle, il peut être nécessaire de prendre en compte les
effets de la gravité. La loi de darcy s’applique alors lorsque le gradient du potentiel hydraulique
entraîne le mouvement du fluide dans le milieu poreux. Elle visualise le champ du potentiel
hydraulique en considérant la différence de potentiel de pression et de hauteur entre les points
de départ et d’arrivée de la ligne d’écoulement. Selon la loi de Darcy, le flux net à travers une
surface poreuse est donné par :

𝑘
𝑢𝑑𝑙 = − (∇𝑝 + 𝜌𝒈) (62)
𝜇
Dans cette équation, 𝒈 est le vecteur de l’accélération gravitationnelle. En outre, dans
le module de la loi de Darcy dans Comsol Multiphysics, 𝒈 peut être définie par (𝒈 = −𝑔∇𝐷).
Où g est la magnitude de l’accélération gravitationnelle en (m/s2), D est l’élévation en (m), ∇𝐷
est le vecteur unitaire dans la direction dans laquelle la gravité agit. 𝜌 est la masse volumique
du fluide exprimée en (kg/m3).

II.1 Dégradation des paramètres de résistance du sol

La structure interne d’une digue est influencée par les infiltrations et l'érosion interne,
ce qui pourrait affecter sa stabilité. Dans ce chapitre, nous avons supposé que la cohésion « C »
et l'angle de frottement « φ » ne sont pas affectés par l'humidité, c'est-à-dire que la cohésion
« C » et l'angle de frottement « φ » changent avec l'érosion interne tout au long du processus.
La cohésion « C » et l'angle de frottement « φ » sont influencés par la porosité et peuvent être
obtenus à l'aide des relations suivantes [159].

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 103

Chapitre V Etude numérique : Modèle 2D
1−𝑛
𝐶 = 𝐶0 (63)
1 − 𝑛0
1−𝑛 (64)
𝜑 = 𝜑0
1 − 𝑛0

II.2 Implémentation des équations du modèle dans comsol multiphysics

Le modèle est intégré sous l’environnement Comsol-multiphysics version 5.5. En effet,
Comsol est un outil de simulation multiphysique permettant de résoudre, de façon couplée, des
équations aux dérivées partielles par la méthode des éléments finis. Il dispose d’une base de
données d’équations permettant de modéliser différents phénomènes physiques comme le
transport des espèces chimiques en milieu poreux, l’écoulement des fluides, le transport de
chaleur en milieu poreux…etc. De manière générale, les étapes de la modélisation sont les
suivantes :

✓ Choix des dimensions de l’espace (2D) ;

✓ Détermination des interfaces utilisées ;
✓ Définition des propriétés physiques des différents paramètres ;
✓ Définition des équations des différentes variables ;
✓ Définition de la géométrie ;
✓ Définition des conditions aux limites et de l’état initial ;
✓ Maillage du domaine ;
✓ Résolution du problème par le solveur et post-traitements.

II.3 Type d’étude

L’étude Comsol permet d’accéder à l’évolution des différentes variables du modèle
telles que la porosité, la concentration des particules érodée en suspension, la perméabilité et la
vitesse de l’écoulement en fonction des différents éléments tels que le temps et l’espace. Dans
notre cas, nous avons effectué des simulations temporelles (en anglais : Time Depend Study)
ensuite une étude paramétrique pour mettre en évidence l’effet de chacun des paramètres.

II.4 Les modules utilisés

Les modules utilisés dans Comsol pour notre étude sont :

✓ L’interface de la loi de Darcy (dl) avec une variable dépendante

de la pression des pores (p), pour l’écoulement en milieu poreux
qui est décrit par l’équation (62).

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 104

Chapitre V Etude numérique : Modèle 2D
✓ L’équation aux dérivées partielles (PDE), avec une variable
dépendante de la porosité (n) pour l’équation de la loi d’érosion
(équation 59) ;
✓ L’équation aux dérivées partielles (PDE), avec une variable
dépendante de la concentration (c) qui décrit l’évolution de la
concentration de particules érodées (équation 60) ;

L’interface de darcy (dl) dans une étude transitoire (darcy’s law (dl))
Le mode de Darcy (dl) utilise la pression (p) comme une variable dépendante. Ce
module se compose d’une équation de la loi de Darcy (équation 62) et d’une équation de
conservation de la masse fluide (équation 65).

𝜕(𝜌̅ . 𝑛) (65)
+ ∇(𝜌̅ . 𝑢𝑑𝑙 ) = 𝑄𝑚
𝜕𝑡
Pour une étude temporelle, la perméabilité et la porosité du système sont évaluées en
fonction du temps. Dans ce module, on peut choisir entre la perméabilité intrinsèque (k) (m2)
ou la conductivité hydraulique (K) (m/s) afin de déterminer la perméabilité du milieu. Dans
notre cas, nous avons choisi la perméabilité intrinsèque qui est exprimée par la formule de
Kozeny-Carman (équation 53). La densité et la viscosité du fluide sont considérées comme la
densité et la viscosité du mélange fluide-particules fluidisées (équations 66 et 67). A noter que
l’effet de la gravité est pris en compte dans l’équation de Darcy (62).

𝜌̅ = (1 − 𝑐)𝜌𝑓 + 𝑐𝜌𝑠 (66)

𝜇 = 𝜇0 (1 + 2.5𝑐) (67)

𝑛3 (68)
𝜅 = 𝜅0
(1 − 𝑛)2

La résolution de cette physique permet de déterminer la vitesse selon la loi de Darcy

(udl) et la pression de l’écoulement (p) comme des fonctions dépendantes du temps et des
cordonnées spatiales.

L’interface (PDE) de la forme « coefficients » (coefficient form pde(c)) avec

la variable « n »
Dans un domaine 2D et pour une variable dépendante (u), l’équation différentielle
partielle d’une forme, qui contient différents coefficients, est présentée comme suit :

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 105

 Chapitre V Etude numérique : Modèle 2D
𝜕 2𝑢 𝜕𝑢
𝑒𝑎 2 + 𝑑𝑎 + ∇. (−𝑐. ∇𝑢 − 𝛼𝑢 + 𝛾) + 𝛽. ∇𝑢 + 𝑎𝑢 = 𝑓 (69)
𝜕𝑡 𝜕𝑡

Nous avons utilisé la forme de l’équation (69) afin de présenter l’équation de la loi
d’érosion (équation 60) en choisissant la variable dépendante (u) comme la porosité n et en
déterminant les coefficients ea, α, da, β, c, γ, a et f. Par identification, on obtient les coefficients
: ea=α=β= γ=a=0, da= 1, c=D (coefficient de diffusion) et avec un terme source (f) qui représente
la partie droite de l’équation (69). D’ailleurs, en terme d’unité, Comsol donne au terme source
(f) une unité par défaut qui est (1/m2), dans notre cas, cette unité doit être adaptée pour
𝜕𝑛
correspondre à celle du terme de la dérivée temporelle de la porosité ( 𝜕𝑡 ), c’est-à-dire (1/s).

L’interface (PDE) de la forme « coefficients » (coefficient form pde(c)) avec

la variable « c »
Le développement de l’équation (59) donne :

𝜕𝑐 𝜕𝑛
𝑛 − 𝑞̅ ∇𝑐 + ∇(D. n. ∇c + D. c. ∇n) = (1 − 𝑐) (70)
𝜕𝑡 𝜕𝑡
En identifiant l’équation (59) et (69), on obtient les coefficients ea = γ= a =0, c =D, β = udl, α =
D. ∇n, da= n et

𝜕𝑛
𝑓= (1 − 𝑐)
𝜕𝑡

La variable dépendante u dans ce cas est remplacée par la concentration (c).

III. Exemple d’application

Pour illustrer la capacité du modèle numérique à traiter des problèmes d'ingénierie à

l'échelle réelle, on a choisi d'analyser un exemple du problème de digue, où l'érosion peut avoir
lieu en raison de gradients hydrauliques élevés. La coupe transversale de la digue est montrée
sur la Figure V. 1. Les valeurs des différentes caractéristiques géométriques sont montrées sur
le Tableau V. 1Tableau V. 1. Cette étude vise à montrer comment l'érosion interne peut affaiblir
la digue, et conduire à l'effondrement de la digue. Dans cette simulation, on suppose que seules
les particules fines sont transportées dans le flux d'eau ; la formation éventuelle d'une érosion
de conduit ne peut donc pas être reproduite. Il convient de souligner que dans cette exemple,

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 106

Chapitre V Etude numérique : Modèle 2D
l'érosion est due à des gradients hydrauliques élevés. La dimension et le maillage par éléments
finis du domaine physique sont indiqués dans la Figure V. 2. Le niveau de l'eau est situé à 2
mètres au-dessous du sommet de la digue, c'est-à-dire à 10 m, on considère la configuration
d'une condition normale.

Figure V. 1 coupe transversale et caractéristiques géométrique de l’exemple étudié

Tableau V. 1 Paramètres d’entrée du modèle 2D

Paramètres Valeur Description

L1 24 [m] Caractéristique géométrique
L2 5 [m] Caractéristique géométrique
L3 24 [m] Caractéristique géométrique
H 12 [m] Caractéristique géométrique
Hw 10 [m] Charge Hydraulique en amont
Hd 5 [m] Charge Hydraulique en aval
λ 1 [1/m] Coefficient d’érosion
n0 0.396 Porosité initiale
k0 2.7436e-13[m2] Perméabilité initiale
c0 10-3 Concentration initiale
µ0 10-3 [Pa.s] Viscosité initiale
ρs 2650 [kg/m3] Masse volumique solide
ρf 1000 [kg/m3] Masse volumique de l’eau
ncr 0.500 Porosité critique
ccr 0.3 Concentration critique
D 3.56*10-6[m2/s] Coefficient de dispersion

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 107

Chapitre V Etude numérique : Modèle 2D
III.1 Choix des points de mesure
Nous avons choisi de suivre la variation locale des différentes variables du modèle
pendant le processus d’érosion interne au niveau du point de sortie de l’écoulement car c’est la
zone des vitesses relativement importante où l’érosion est plus importante que d’autres endroits.
La variation linéaire des différents paramètres est évaluée le long de la base de la digue pour
montrer la distribution linéaire.

III.2 Maillage du domaine

Le maillage adopté est de type triangulaire extra fin avec la méthode de triangulation de
Delaunay. Ce type de maillage a été appliqué après des tests numériques préliminaires de
précision du calcul. Le maillage final est constitué de 1289 éléments de domaines et 111
éléments de frontières triangulaires. La Figure V. 2 montre le maillage choisi pour le domaine
entier.

Figure V. 2 Maillage du domaine

III.3 Conditions aux limites

Dans le logiciel Comsol, les conditions aux limites peuvent être définies pour chaque
interface. Pour l’interface de Darcy, les conditions aux limites sont appliquées sur les frontières
indiquées dans la Figure V. 3 avec :

✓ Une charge hydraulique constante (Hw) en amont ;

✓ Une charge hydraulique constante (Hd) en aval ;
✓ Pas d’écoulement à travers les autres frontières.

Pour l’interface EDP, nous avons appliqué une condition de Dirichlet (c=c0) en amont
de la géométrie et un flux nul sur les autres frontières.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 108

Chapitre V Etude numérique : Modèle 2D

Figure V. 3 Les conditions aux limites

III.4 Résultats et discussions

La pression

Le profil des pressions obtenu est montré sur la Figure V. 4 avec les différentes lignes
d’isovaleurs. Le niveau de la ligne de saturation est aussi illustré dans la même figure.

Figure V. 4 Profil spatial de la pression interstitielle à la fin du processus d’érosion

La porosité

La Figure V. 5 montre la distribution spatiale de la porosité du milieu à différents

instants : 500, 1500, 2500 et 3600 jours. Les résultats montrent que la porosité augmente au
pied de la digue puis elle tend à se diffuser en amont, ce qui se traduit par une érosion régressive.
Les vitesses d’écoulement importantes au pied de la digue (Figure V. 11) impliquent des
contraintes de cisaillement importantes qui sont le facteur important pour le déclenchement de
l’érosion à cet endroit. Cependant, la définition du paramètre ncr a permet de limiter

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 109

Chapitre V Etude numérique : Modèle 2D
l’amplification de la porosité dans la digue. Cela permet de considérer que le type d’érosion
dans ce cas est une érosion interne par suffusion.

Figure V. 5 Distribution spatiale de la porosité à différents instants

Ligne phréatique

La position de la ligne de saturation avant et à la fin du processus d’érosion au sein de

la digue est montrée sur la Figure V. 6. On remarque une baisse du niveau de pression ce qui est
relié au départ de particules ce qui laisse plus de vide dans le milieu, ce vide est ensuite remplacé
par l’eau traversant la digue.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 110

Chapitre V Etude numérique : Modèle 2D

Figure V. 6 Position de la ligne de saturation avant et après érosion

Teneur en fines

Dans son étude, Yang [89] a développé un modèle basé sur la théorie des milieux poreux
en considérant un milieu constitué de quatre (4) phases : squelette solide, solide érodable,
particules solides fluidisées et un fluide pur. Il a donné la formule suivante pour le calcul
explicite de la teneur en fines (𝑓𝑐 ) au sein de l’ouvrage

(1 + 𝜀𝑣 )(1 − 𝑛0 )(1 − 𝑓𝑐0 )

𝑓𝑐 = 1 − (71)
1−𝑛
Où 𝑛 est la porosité du milieu, 𝑛0 et 𝑓𝑐0 sont la porosité et la teneur en fines initiales
respectivement et 𝜀𝑣 est la déformation volumique.

Dans la présente étude, nous avons négligé les déformations de volume, donc la relation (56)
devient

(1 − 𝑛0 )(1 − 𝑓𝑐0 )
𝑓𝑐 = 1 − (72)
1−𝑛
La Figure V. 7 montre la distribution spatiale de la teneur en fines calculée par la relation
(57). Les résultats montrent que le coté aval de l’ouvrage est plus touché par l’érosion interne
où on remarque une diminution importante. La teneur en fine à cet endroit a diminuée jusqu’à
environ 15%, ce qui implique une diminution de 50% de sa valeur initiale qui était de 30%.
Quant au coté amont de l’ouvrage, il n’a pas subi une érosion importante du fait que les vitesses

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 111

Chapitre V Etude numérique : Modèle 2D
d’écoulement à cet endroit sont moins importantes, ce qui implique que les contraintes de
cisaillement ne sont pas très importantes.

En conclusion, la perte continue de fines dans le côté aval de la digue pourrait entraîner
la rupture de la digue. Dans ce cas, la diminution de la teneur en fines près du pied de la digue
pourrait provoquer le glissement global de la digue. Ce problème sera d’ailleurs l’objet de la
dernière partie de ce chapitre.

Figure V. 7 Distribution spatiale de la teneur en fines à différents instants

Perméabilité intrinsèque

La variation de la perméabilité du milieu évaluée par l’équation de Kozeny-

Carman (équation 68) est montrée sur la Figure V. 8 aux instants : 500, 1500, 2500 et
3600 jours. Les résultats montrent une augmentation quatre (4) fois plus que sa valeur
initiale au niveau du pied de la digue. Cette augmentation diminue en allant vers l’amont
de l’ouvrage où la perméabilité reste prèsque inchangée.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 112

Chapitre V Etude numérique : Modèle 2D

Figure V. 8 Distribution spatiale de la perméabilité à différents instants

La concentration

La Figure V. 9 montre la répartition spatiale de la concentration des particules fines à

différents instants. On remarque qu’elle évolue dans le temps jusqu’à atteindre une
valeur maximale où le mélange fluide-solide devient plus visqueux et sa densité
augmente. Ensuite la concentration tend à diminuer jusqu’à ce qu’elle devient nulle suite
au lavage de toutes les particules érodées. La Figure V. 15 représente le suivi de
l’évolution temporelle de la concentration au point de sortie de l’écoulement situé au
pied de la digue, cette courbe montre plus clairement comment se fait le processus
d’érosion. Cette courbe contient deux partie principales, une première où la
concentration augmente jusqu’à une valeur maximale, c’est le rpocessus de détachement
des particule, et une deuxième où la concentration diminue, ce qui implique l’érosion
de toutes les particules.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 113

Chapitre V Etude numérique : Modèle 2D
Figure V. 9 Distribution spatiale de la concentration à différents instants

La cohésion

1−𝑛
La variation de la cohésion du milieu a été évaluée par l’équation 𝐶0 1−𝑛 . Cette équation
0

implique une variation proportionnelle de la cohésion en fonction de la porosité du

milieu. L’augmentation de la porosité est supposée un facteur qui dégrade la résistance
du sol. Cette hypothèse s’appui sur le fait que la cohésion d’un sol est assurée par ses
particules fines, ce qui implique que leurs départ conduit à la dégradation de celle-ci. En
effet, leurs départ donne un milieu plus poreux. Le profil spatial de la cohésion à
différents instants est monté sur la Figure V. 10. La cohésion en aval de l’ouvrage a atteint
une valeur de 16.90 kPa tandis qu’elle est restée constante au coté amont avec une valeur
de 20 kPa. Ces valeurs qui dépendent de la poroité dépendent aussi du coefficient
d’érosion λ et de la porosité critique ncr qui représente les facteurs conducteurs dans le
modèl. La dégradation de la cohésion au coté amont peut être un signe de début de
l’effondrment d’un ouvrage hydraulique tel que montré dans l’étude de cas effectué par
Zhang LM, Chen [124], où ils ont identifié les différents mécanismes qui ont causés la
rupture du barrage de Gouhou en Chine.

Figure V. 10 Distribution spatiale de la cohésion à différents instants

Vitesse de Darcy

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 114

Chapitre V Etude numérique : Modèle 2D
Le champ de vitesse de darcy est montré sur la Figure V. 11 pour différents instants.
L’’érosion qui a commencée au pied de la digue observée par l’augmentation importante de la
porosité a été le résultat des vitesses d’écoulement importantes à cet endroit comme il est montré
sur cette figure. Les vitesses sont d’autant plus faibles en s’approchant de l’amont de l’ouvrage
où on n’a pas observé d’érosion importante.

Figure V. 11 Distribution spatiale du champ de vitesse de Darcy à différents instants

III.5 Variation temporelle des différents paramètres au pied de la digue

La variation temporelle des différents paramètres a été évaluée au niveau du point de
sortie de l’écoulement situé au pied de la digue. Nous avons choisi cet endroit car c’est à ce
niveau que se produit plus d’érosion et permet de donner une idée générale sur la manière par
laquelle se produit le phénomène. Les figures (Figure V. 12, Figure V. 13, Figure V. 14, Figure V.
15, Figure V. 16) montrent les profils de la porosité, la perméabilité, teneur en fine, la
concentration et la vitesse de darcy respectivement. Ces résultats sont similaires aux résultats
obtenus par le modèle d’érosion 1D présenté au chapitre 4. Quant à la Erreur ! Source du r
envoi introuvable., elle montre l’interaction entre la concentration des particules fines érodées
et la porosité. Les deux paramètres augmentent simultanément jusqu’à ce que la porosité atteint
sa valeur critique qui est supposée comme le seuil d’érosion. Une fois la porosité critique est
atteinte, la concentration devient nulle pour la même valeur de la porosité où l’espace poreux
crée par le départ des particules est remplacé par le fluide (l‘eau dans notre cas).

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 115

Chapitre V Etude numérique : Modèle 2D

Figure V. 12 Evolution temporelle de la porosité au pied de l’ouvrage

Figure V. 13 Evolution temporelle de la perméabilité au pied de l’ouvrage

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 116

Chapitre V Etude numérique : Modèle 2D

Figure V. 14 Evolution temporelle de la teneur en fines au pied de l’ouvrage

Figure V. 15 Evolution temporelle de la concentration au pied de l’ouvrage

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 117

 Chapitre V Etude numérique : Modèle 2D

Figure V. 16 Evolution temporelle de la vitesse de Darcy au pied de l’ouvrage

IV. Etude paramétrique

Nous allons dans cette section présenter une étude paramétrique sur l’influences des
différents paramètres du modèle sur le processus d’érosion. Nous avons donc choisi d’étudier
l’effet du coefficient d’érosion (λ), de la porosité critique (ncr) ainsi que l’effet de la charge
hydraulique (Hw). Ces trois paramètres ont le plus d’impact sur le modèle, la charge hydraulique
est d’ailleurs le déclencheur du processus, la porosité critique définie le seuil de l’érosion au
bout duquel le système rentre en régime permanent et le coefficient d’érosion traduit le potentiel
d’érosion dans le milieu.

IV.1 Effet du coefficient d’érosion

Le coefficient d’érosion (λ) est probablement lié à l’érodibilité interne du sol qui dépend
essentiellement des éléments structuraux du milieu poreux (granulométrie du sol, densité du
sol... etc.) Ce paramètre est utilisé largement dans les modèles d’érosion proposé par plusieurs
auteurs. D’après Nguyen [160], Bendahmane [161] et Vardoulakis [85], ce coefficient peut

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 118

Chapitre V Etude numérique : Modèle 2D
plutôt avoir une influence sur l’évolution de la porosité du milieu (n) au cours du temps,
l’évolution de la concentration des particules érodées en suspension (c).

La Figure V. 17 montre les résultats de l’étude paramétrique du coefficient d’érosion λ.

L’étude paramétrique est effectuée avec quatre (4) valeurs différentes de λ : 1, 1.5, 2 et 3. Les
résultats obtenus montrent que le coefficient d’érosion tend à accélérer le processus d’érosion
et conduit à atteindre le seuil d’érosion plus rapidement. Comme nous avons supposé que les
particules soient déjà détachées, il suffit que le coefficient d’érosion ait une valeur qui puisse
entrainer les particules érodées, et le processus se poursuit jusqu’à atteindre le seuil de l’érosion.

(a)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 119

Chapitre V Etude numérique : Modèle 2D

(b)

(c)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 120

Chapitre V Etude numérique : Modèle 2D

Figure V. 17 Etude paramétrique du coefficient d’érosion : (a) Porosité, (b) Perméabilité

intrinsèque, (c) Concentration, (d) Vitesse de Darcy

IV.2 Effet de la porosité critique

La porosité critique ncr est un coefficient proposé pour définir le seuil de l’érosion. Nous
allons montrer dans cette partie l’influence de ce paramètre en considérant trois valeurs égales
à : 0.45, 0.5 et 0.55. Le paramètre ncr tient compte du fait que seules les particules fines d’un
sol peuvent s’érodées. En effet, la porosité du milieu atteint une valeur critique au-delà de
laquelle il n’y aura plus d’érosion. Cette valeur critique de porosité dépend du taux de fines
érodables, Yang [162] a présenté dans son étude un modèle d’érosion basé sur la théorie des
milieux poreux où il a considéré un milieu à quatre phases : phase solide non érodable, phase
liquide, phase solide fluidisée et une phase solide érodable.

La Figure V. 18 montre les résultats de l’étude paramétrique de la porosité critique. On

peut constater sur le profil de la teneur en fine à la fin du processus d’érosion que, pour un sol
avec une teneur initiale en fine égale à 30% atteint des valeurs de 25%, 16% et 6% pour des
valeurs de ncr égale à 0.45, 0.5 et 0.55 respectivement. Les valeurs de ncr se diffèrent d’un
matériau à un autre en fonction de son érodabilité.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 121

Chapitre V Etude numérique : Modèle 2D

(a)

(b)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 122

Chapitre V Etude numérique : Modèle 2D

(c)

(d)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 123

Chapitre V Etude numérique : Modèle 2D

(e)

Figure V. 18 Etude paramétrique de la porosité critique : (a) Porosité, (b) Perméabilité

intrinsèque, (c) Concentration, (d) Vitesse de Darcy, (e) Teneur en Fine

Effet de la charge hydraulique

Le troisième paramètre étudié dans cette partie est l’influence de la charge hydraulique
(Hw) sur l’érosion. Quatre (4) charges hydrauliques ont été testées : Hw=9m, Hw=10m, Hw=11m
et Hw=12m. L’évolution temporelle des différents paramètres évaluée au pied de l’ouvrage est
montrée sur la Figure V. 19. On remarque d’après ces résultats que la charge hydraulique influe
sur l’érosion par l’accélération du processus. On remarque que pour un instant donné, plus la
charge hydraulique est importante, plus les valeurs de la concentration, la porosité, la vitesse
de Darcy sont importantes.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 124

Chapitre V Etude numérique : Modèle 2D

(a)

(b)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 125

Chapitre V Etude numérique : Modèle 2D

(c)

(d)

Figure V. 19 Etude paramétrique de la charge hydraulique (a) Porosité, (b) Perméabilité intrinsèque,
(c) Concentration, (d) Vitesse de Darcy

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 126

 Chapitre V Etude numérique : Modèle 2D
Nous allons nous intéresser maintenant sur l’effet de l’érosion interne sur la stabilité au
glissement de la digue déjà présentée. La stabilité au glissement étant en relation avec les
paramètres de résistance au cisaillement (C et phi), nous avons déterminé la variation
temporelle de ces deux paramètres au fur et à mesure que l’érosion évolue au sein de la digue.

V. Influence de l’érosion interne sur la stabilité vis-à-vis au

glissement

Dans la partie précédente, nous avons présenté une simulation du phénomène de

l’érosion interne au sein d’un ouvrage hydraulique en terre (digue). Nous avons par la suite
tracé les variations spatiales et temporelles des différents paramètres qui interviennent dans le
modèle. Afin d’évaluer la dégradation de la stabilité globale d’un tel ouvrage, nous avons choisi
de le faire en considérant que, la cohésion et l’angle de frottement interne du sol constituant
l’ouvrage se dégradent au fur et à mesure que la porosité du milieu augmente. Pour ce faire,
nous avons suivi la variation dans le temps et dans l’espace de ces deux paramètres en
implantant une loi de dégradation linéaire. Vu que le phénomène se diffuse beaucoup plus en
longueur qu’en élévation, nous avons négligé les variations des paramètres en élévation, nous
avons donc évalué la cohésion et l’angle de frottement le long de la digue. Une fois les courbes
sont obtenues, on les introduit comme paramètres d’entrée de l’étude de stabilité sur le code de
calcul SLOPE/W. L’analyse de stabilité est effectuée ensuite par les diverses méthodes
d'équilibre limite pour déterminer le facteur de sécurité et sa variation en fonction du temps.

V.1 Analyse de stabilité

Les problèmes liés aux ruptures des pentes naturelles et artificielles posent souvent de
sérieux défis en matière de géotechnique. En général, une surface de sol inclinée et exposée qui
n’est pas retenue peut-être sujette à des mouvements de masse dus aux forces gravitationnelles.
Les contraintes de cisaillement résultantes, induites le long d'une surface de rupture potentielle
ou connue, pourraient dépasser la résistance au cisaillement du sol et provoquer la rupture de
la pente. Le rapport entre la résistance au cisaillement existante et les contraintes de cisaillement
induites sur une surface de rupture potentielle est le plus souvent appelé le facteur de sécurité.
L'objectif de toute analyse de stabilité à l'équilibre limite est de déterminer ce facteur, par lequel
la résistance du sol est divisée ou réduite, pour amener la pente au seuil d'instabilité.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 127

Chapitre V Etude numérique : Modèle 2D
L’analyse de la stabilité des ouvrages en terre est le plus ancien type d’analyse numérique dans
le domaine de l’ingénierie géotechnique. L'idée de discrétiser une masse glissante en tranches
a été introduite au début du XXe siècle. En 1916, Petterson [163] a présenté l'analyse de stabilité
du quai Stigberg à Göteborg, Suède où la surface de glissement a été considérée comme
circulaire et la masse de glissement a été divisée en tranches. Au cours des décennies suivantes,
Fellenius [164] a introduit la méthode des tranches ordinaires ou suédoises. Au milieu des
années 1950, Janbu [165] et Bishop [166] ont développé des avancées dans la méthode.
L’apparition des ordinateurs dans les années 1960 ont permis de traiter plus facilement les
procédures itératives inhérentes à la méthode qui a conduit à des formulations
mathématiquement plus rigoureuses telles que celles développées par Morgenstern et Price [65]
et par Spencer [67]. L'une des raisons pour lesquelles la méthode de l'équilibre limite a été
adoptée si facilement, c'est que les solutions pouvaient être obtenues par des calculs manuels.
Simplifier l'hypothèse a dû être adoptée pour obtenir des solutions, mais le concept de division
numérique d'un grand corps en petits pièces destinées à l'analyse étaient plutôt nouvelles à
l'époque.

V.2 Principes fondamentaux de l'équilibre limite

Bases de la méthode
De nombreuses techniques de solution différentes pour la méthode des tranches ont été
développées au fil des temps. En gros, elles sont toutes très similaires. Les différences entre les
méthodes dépendent des éléments suivants : quelles équations de la statique sont incluses et
satisfaites, et quelles forces entre les tranches sont incluses et quelle est la relation supposée
entre le cisaillement entre les tranches et les forces normales ? Figure V. 20 illustre une masse
de glissement typique discrétisée en tranches et les forces possibles sur la tranche. Les forces
normales et les forces de cisaillement agissent sur la base de la tranche et sur les côtés de la
tranche.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 128

Chapitre V Etude numérique : Modèle 2D

Figure V. 20 Discrétisation en tranche et les forces agissant sur une tranche

V.3 Méthodes de calcul du facteur de sécurité

Au fil du temps, de nombreuses méthodes différentes ont été développées pour calculer
les facteurs de sécurité. Toutes les méthodes sont basées sur des formulations d'équilibre limite,
à l'exception de la méthode des éléments finis, qui utilise des contraintes calculées par éléments
finis.

Les méthodes utilisées dans cette étude pour évaluer le facteur de sécurité sont les
suivantes :

✓ Méthode simplifiée de Bishop ;

✓ Méthode de Morgenstern-Price ;
✓ Méthode de Spencer ;
✓ Méthode de Janbu.

La description détaillée de ces méthodes peut être consulté sur le manuel d’utilisation
du code de calcul SLOPE/W utilisé dans cette étude [167].

V.4 Paramètres de résistance au cisaillement

Cohésion
Nous introduisons la relation (63) qui exprime la dégradation de la cohésion du matériau
qui est proportionnelle à la porosité [159].

Angle de frottement interne

Nous introduisons la relation (64) qui exprime la dégradation de l’angle de frottement
interne du matériau qui est proportionnelle à la porosité [159].

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 129

Chapitre V Etude numérique : Modèle 2D
V.5 Cohésion et d’angle de frottement interne
Le profil de cohésion et d’angle de frottement interne introduit dans l’étude de stabilité
de la digue sont ceux obtenus avec la simulation de la dégradation de la cohésion due à l’érosion
interne avec les paramètres d’entrés montrés dans le Tableau V. 1. La géométrie de la digue
étudiée dans cette partie est la même que celle étudiée dans la partie précédente. Les relations
qui expriment la dégradation des paramètres (C et φ) sont données par les relations (73, 74).
Les distributions spatiales de la cohésion et de l’angle de frottement interne sont introduites en
considérant qu’elles varient seulement dans la direction x. Ce profil est tracé pour chaque
instant, la Figure V. 21 montre les profils spatiaux de la cohésion « C » et d’angle de frottement
« φ » interne utilisé.

Tableau V. 2 Paramètres d’entrée pour évaluer la stabilité après érosion

Paramètres Valeur Description

Hw 10 [m] Charge Hydraulique en amont
Hd 5 [m] Charge Hydraulique en aval
λ 1 [1/m] Coefficient d’érosion
n0 0.396 Porosité initiale
k0 2.7436e-13[m2] Perméabilité initiale
c0 10-3 Concentration initiale
µ0 10-3 [Pa.s] Viscosité initiale
ρs 2650 [kg/m3] Masse volumique solide
ρf 1000 [kg/m3] Masse volumique de l’eau
ncr 0.500 Porosité critique
ccr 0.3 Concentration critique
D 3.56*10-6[m2/s] Coefficient de dispersion
C0 10 [kPa] Cohésion initiale
φ 25 [°] Angle de frottement interne initial

Une étude paramétrique sur l’effet du coefficient d’érosion λ a été effectuée pour
montrer son influence sur la réduction de la stabilité de la digue. Les valeurs de λ étudiées sont :
1, 1.5 et 2. La Figure V. 21 montre les profils de C et φ à différents instants et pour les différentes
valeurs de λ.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 130

Chapitre V Etude numérique : Modèle 2D

(a)

(b)

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 131

Chapitre V Etude numérique : Modèle 2D

(c)

Figure V. 21 Profil de la cohésion et de l’angle de frottement utilisés dans l’analyse de stabilité

pour différentes valeurs de λ et à différents instants : (a) λ=1, (b) λ=1.5, (c) λ=2.

V.6 Evaluation du facteur de sécurité

Pour évaluer l’influence de l’érosion interne sur la stabilité globale de la digue, nous
avons tout d’abord simuler le processus d’érosion dans l’ouvrage étudié. Ensuite, nous avons
tracé les profils de C et φ avec les relations (63,64) le long de la digue. La Figure V. 21 montre
le profil spatial de C et φ à différents instants.

Le code de calcul GEO-SLOPE offre une possibilité d’introduire des fonctions de

cohésion et d’angle de frottement interne. Nous avons donc introduit les courbes de C et φ
obtenues par la simulation de l’érosion interne dans l’ouvrage étudié.

La Figure V. 22 montre la surface de glissement de la pente aval de la digue à différents

instants ainsi que le facteur de sécurité correspondant obtenu en utilisant la méthode de
Morgensterne-Price. La Figure V. 22a montre aussi le profil des pressions interstitielles obtenue
par une simulation d’un écoulement transitoire par le module SEEP/W du code de calcul
GEOSTUDIO. Dans cette analyse, nous avons pris la perméabilité et la charge hydraulique
comme étant constantes. La ligne phréatique est aussi obtenue par la même analyse, on peut
l’observer sur toutes les figues. Les Figure V. 22b, V.22c, V.22d, V.22e montrent la distribution
spatiale de la cohésion du milieu après introduction des courbe obtenues après érosion interne
(Figure V. 21).

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 132

Chapitre V Etude numérique : Modèle 2D

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 133

Chapitre V Etude numérique : Modèle 2D

Figure V. 22 Surface de glissement critique pour plusieurs instants : (a) à l’état initial (b) à
500 jours (c) à 1500 jours (d) à 2500 jours (e) à 3600 jours

La Figure V. 23 illustre la variation temporelle du facteur de sécurité calculé par les

différentes méthodes. Le facteur de sécurité initial calculé par la méthode de Morgenstern-Price,
Bishp, Spencer et Janbu est égale à 1.118, 1.109, 1.113 et 1.045 respectivement.

Figure V. 23 Comparaison de la réduction du facteur de sécurité par les différentes méthodes

Les valeurs du facteur de sécurité calculé pour les différents instants sont montrées dans
le Tableau V. 3. En comparant les valeurs initiales et finales du coefficient de sécurité calculées
par les différentes méthodes, on constate une réduction de la stabilité de 16%. La valeur critique
du facteur de sécurité qui est de 1 calculée par les méthodes (Morgenstern-Price, Bishop et
Spencer) est atteinte presqu’à 1500j, ce qui correspond à environ quatre ans, tandis qu’elle est

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 134

Chapitre V Etude numérique : Modèle 2D
atteinte bien avant paour la méthode de Janbu. Ces résultats montrent clairement que l’érosion
interne peut conduire à la déstabilisation de la pente par son impact sur les paramètres de
résistance au cisaillement (Cohésion et angle de frottement interne).

Tableau V. 3 Variation temporelle du facteur de sécurité

Facteur de sécurité
Méthode de calcul
t = 0j t = 500j t = 1500j t = 2500j t = 3600j
Morgenstern-Price 1.118 1.105 1 0.947 0.932
Bishop 1.109 1.096 0.996 0.942 0.926
Spencer 1.113 1.1 0.996 0.943 0.928
Janbu 1.045 1.035 0.944 0.893 0.877

La Figure V. 24 montre l’effet du coefficient d’érosion λ sur la variation temporelle du facteur

de sécurité. La dégradation du facteur de sécurité s’amplifie en augmentant le coefficient λ.
Seule la valeur de λ=0.5 qui ne conduit pas à une valeur du facteur de sécurité inférieure à 1.
Les valeurs de λ supérieure à 1.5 conduisent tous à une même valeur finale du facteur de sécurité
qui est égale à 0.92. La comparaison des valeurs du facteur de sécurité obtenues à un instant
donné pour les différentes valeurs de λ montre que ce coefficient tend à accélérer le processus
de déstabilisation.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 135

 Chapitre V Etude numérique : Modèle 2D
Figure V. 24 Variation temporelle du facteur de sécurité pour différentes valeurs de λ

VI. Conclusion

Ce chapitre a été une tentative d’apporter une contribution à l’approche numérique pour
quantifier l’influence de l’érosion interne sur les infrastructures géotechniques. Le travail a
consisté à modéliser l'érosion interne du sol et le transport de particules fines par le fluide à
travers l'échange de masse entre les phases solide et fluide. Les équations différentielles
déterminantes ont été formulées sur la base du bilan massique de trois constituants supposés :
le squelette solide, les particules fines érodées et le fluide pur, qui ont ensuite été implémentés
dans le code de calcul par éléments finis Comsol Multiphysics.

L'exemple d'une digue homogène a démontré l'applicabilité du modèle pour simuler

l'évolution de l'érosion interne en fonction du temps pour un problème à l’échelle réelle. Les
effets de l'érosion interne ont entraîné une augmentation de la porosité du sol et une diminution
de la teneur en fines, ce qui a entraîné une diminution de la résistance au cisaillement de la zone
érodée, évaluée par une loi qui relie la cohésion et l’angle de frottement interne du sol à sa
porosité.

Une étude paramétrique sur l’effet de plusieurs paramètres a été présentée et discutée
dans ce chapitre, ce qui a permis de visualiser l’effet de chaque paramètre. Cette étude a montré
que le coefficient d’érosion a plus d’impact sur le processus d’érosion où sa détermination est
primordiale.

Nous avons tenté dans ce chapitre à travers l’étude d’un cas de digue de prédire la
dégradation de sa stabilité au glissement due à l’érosion interne. Pour ce faire, nous avons
effectué des simulations d’érosion interne au sein de cet ouvrage puis évaluer la dégradation
des caractéristiques de résistance au cisaillement (cohésion et angle de frottement) du sol la
constituant. Ensuite, nous avons procédé au calcul de la stabilité au glissement de la digue en
prenant les valeurs de la cohésion et angle de frottement obtenues après érosion comme
paramètres d’entrées de cette étude. Ensuite, une comparaison des différentes méthodes de
calcul de la stabilité des pentes a été effectuée.

Les résultats ont montré que la stabilité globale d’un tel ouvrage peut être altérée par le
phénomène d’érosion et conduire à la rupture de l’ouvrage. Une étude paramétrique sur l’effet

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 136

Chapitre V Etude numérique : Modèle 2D
du coefficient d’érosion λ sur le facteur de sécurité est effectuée. Cette étude a montré que ce
coefficient tend à accélérer le processus de déstabilisation.

La comparaison des différentes méthodes de calcul de stabilité des pentes a montré que
la méthode de Janbu sous-estime le facteur de sécurité de la pente par rapport aux autres
méthodes (Morgenstern-Price, Bishp et Spencer). La méthode de Janbu a montré un temps de
rupture plus antérieure que les autres méthodes.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 137

Conclusion générale
Conclusion générale
Les géo-structures telles que les barrages, les levées ou les digues sont soumises à des
flux d'infiltration variant dans le temps et dans l'espace. L'eau qui circule dans ces milieux
poreux peut entraîner le détachement et le transport d'une partie des particules de sol à l'intérieur
de ces structures ou de leurs fondations. Ce problème est généralement appelé érosion interne.
Le terme « suffusion », désignant l’érosion interne de volume, fait référence au détachement et
au transport de particules plus fines à travers une matrice de sol poreuse plus grossière en raison
de l'écoulement des eaux d'infiltration. La suffusion peut modifier avec le temps les propriétés
hydrauliques et mécaniques des sols et peut déclencher des dommages importants sur ces
structures et conduire à leur effondrement.

Les études expérimentales et numériques sur l'érosion interne ont été soigneusement
passées en revue au début de ce travail. Les tests de laboratoire et les analyses numériques à
l'échelle de l'échantillon jouent, sans aucun doute, un rôle important dans la pratique de
l'ingénierie en déterminant la susceptibilité du matériau à l'érosion interne et dans la
connaissance du mécanisme physique à l'origine des phénomènes d'érosion interne. Toutefois,
ils ne permettent pas de déterminer où et quand l'action du flux d'infiltration peut déclencher
l'instabilité interne du sol granulaire et les réponses mécaniques ultérieures à l'échelle d'une
structure hydraulique. Peu d'essais sur modèle à grande échelle ou d'études numériques à
l'échelle d'une structure entière sont disponibles dans la littérature.

Cette recherche a donc tenté de contribuer à la conception et à la durabilité des

structures d'ingénierie géotechniques et hydrauliques, avec un accent particulier sur les barrages
en remblai, les levées et les digues. Elle consistait principalement en une partie expérimentale
qui vise à étudier l’effet de la suffusion sur les paramètres de résistance au cisaillement des sols,
et une partie numérique qui vise à développer un modèle numérique de la suffusion en
introduisant, le couplage des phénomènes hydrauliques et mécaniques puis évaluer la stabilité
de la digue vis-à-vis du glissement.

L’effet de l’érosion interne sur les paramètres de résistance au cisaillement (cohésion et

angle de frottement interne) a été étudié en mettant en œuvre une campagne expérimentale, où
des essais de cisaillement directe avant et après érosion ont été réalisés, en utilisant l’appareil
de cisaillement directe. Afin de mesurer les paramètres de résistance au cisaillement après
érosion interne, la boite de Casagrande a été modifiée de manière à pouvoir appliquer un
écoulement à travers l’échantillon. Au total, neuf échantillons de sol ont été testés. Ils sont

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 138

composés principalement, de sable qui représente le squelette solide, mélangé avec différents
types de fines à différents pourcentages. Les résultats ont montré que l’érosion interne affecte
considérablement la cohésion et influe positivement sur l’angle de frottement interne.

La partie numérique est réalisée grâce au développement d'un modèle d'érosion en

considérant un milieu triphasique construit dans le cadre de la théorie des milieux poreux
continus. Les équations différentielles, régissant le problème physique, ont été implémentées
dans un code de calcul par éléments finis. En conjonction avec une mise en œuvre et des
simulations efficaces par éléments finis, il était important d'améliorer la capacité de prédiction
de la suffusion des sols et de faciliter la détection fiable et l'alerte précoce du dysfonctionnement
éventuel des structures hydrauliques, des infrastructures urbaines, des ruptures de pente et des
glissements de terrain.

A travers un exemple de digue, les possibilités du modèle pour simuler l’évolution de

l’érosion interne avec le temps ont pu être illustrées. En particulier, cet exemple a montré
comment les effets de l’érosion interne pouvaient conduire à une augmentation de la porosité
du sol et à la dégradation des propriétés mécaniques des digues dans le temps et dans l’espace.
La stabilité de la pente aval de la digue a été évaluée par la méthode d’équilibre limite. Les
résultats montrent une réduction du facteur de sécurité sous l’effet de l’érosion interne,
augmentant le risque d’instabilité mécanique de la digue.

1- L’érosion interne qui entraîne une évolution de la porosité du matériau conduit, de

ce fait, au changement de sa perméabilité. En conséquence, la pression interstitielle
et les contraintes effectives sont redistribuées au sein de différents matériaux ;
2- La perte de particules fines due à l’écoulement d’infiltration, phénomène qui affecte
le comportement mécanique du sol par une modification de la porosité (l’indice des
vides).

Sur la base de la théorie du milieu poreux, le matériau saturé a été considéré comme un
système composé de 3 constituants : le tissu stable du squelette solide, les particules fluidisées
et la phase fluide. Les fines se comportent comme un matériau fluide (décrit comme des
particules fluidisées). L'approche consistait à modéliser l'érosion du squelette du sol, le
transport par le flux d'eau et la filtration des particules fines par l'échange de masse entre les
phases solide et fluide.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 139

 Les effets de l'érosion interne ont conduit à une augmentation du taux de vide (porosité)
du sol et à une diminution de la teneur en fines, ce qui, à son tour, a entraîné une diminution de
la résistance au cisaillement de la zone érodée, calculée par une relation linéaire entre la porosité
et la cohésion.

La perte continue de fines dans le côté aval de la digue pourrait entraîner la rupture de
la digue. Dans ce cas, la diminution de la teneur en fines près du pied de la digue pourrait
provoquer le glissement global de la digue. L’étude de ce problème a montré que le facteur de
stabilité de la digue peut être réduit de 16% de sa valeur initiale.

Perspectives

Cette recherche a tenté de développer un modèle de suffusion à l'échelle d’une structure

entière pour faciliter la détermination de comment, où et quand l'action d'un flux d'infiltration
pourrait déclencher l'instabilité interne du sol granulaire et la réponse mécanique des structures
en terre. Il pourrait être appliqué pour analyser l'influence de la suffusion dans les structures
hydrauliques (barrages, digues, etc.), les infrastructures urbaines (pipelines, tunnels, etc.), les
ruptures de pente et les glissements de terrain.

✓ Les analyses de la suffusion à l'intérieur d'une digue ont été choisies comme
exemple dans cette recherche. Cependant, il convient de noter que ces analyses
ont été effectuées avec une hauteur d'eau constante. Les matériaux de la digue et
de la fondation ont été supposés être initialement homogènes. En réalité, la
hauteur d'eau peut varier, lentement ou rapidement, en fonction des
précipitations ou des activités humaines. Le chemin d'infiltration et la position
de la ligne de saturation peuvent être plus complexes et évoluer avec le temps
dans un domaine hétérogène, ce qui augmentera éventuellement l'incertitude
quant à l'emplacement de l'effondrement. Ces aspects devraient être pris en
compte dans les travaux futurs pour des études de cas réels.

Dans ce qui suit, des recommandations pour les futurs travaux de recherche sont
présentées concernant à la fois les aspects expérimentaux et numériques.

✓ Des essais au laboratoire plus poussés sont nécessaires pour identifier les
facteurs influençant la suffusion et pour saisir l'impact de l'érosion interne sur le
comportement mécanique des sols. Idéalement, la conception de l'appareil

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 140

 d'essai d'érosion devrait permettre de mesurer la déformation de l'échantillon, la
perte de masse de l'échantillon et l'évolution de la conductivité hydraulique le
long de l'échantillon afin d'améliorer la modélisation.
✓ La physique, à l'origine des phénomènes de suffusion, reste floue car les résultats
expérimentaux dans la littérature présentent des contradictions. D'une part, la
préparation de l'échantillon pourrait introduire des perturbations significatives
dans le spécimen ; d'autre part, la nécessité de disposer d'observations
microscopiques actualisées en temps réel du tissu du sol érodé pour mieux
comprendre le comportement hydro-mécanique des sols érodés.
✓ En tirant profit de l'approche discrète (DEM) qui permet une meilleure
considération et compréhension du comportement du matériau à l'échelle du
grain, une alternative aux études expérimentales est d'améliorer les approches
numériques hydro-mécaniques couplées existantes basées sur une formulation
discrète (DEM-CFD, DEM-LBM, DEM-PFV, etc.). Cela permettrait de
visualiser la migration des particules à l'échelle du grain et de réaliser des études
paramétriques à un coût relativement faible.
✓ Développer une approche de modélisation multi-échelle et multi-physique dans
le cadre de la théorie des milieux continus dans laquelle la multi-échelle désigne
un modèle basé sur la micromécanique.

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 141

 Références

[1] (CFGB) CF des GB (1997) Internal erosion: typology, detection, repair

[2] Foster M, Fell R, Spannagle M (2000) The statistics of embankment dam failures and
accidents. Can Geotech J 37:1000–1024. [Link]

[3] Costa JE, Schuster RL (1988) FORMATION AND FAILURE OF NATURAL DAMS.
Bull Geol Soc Am 100:1054–1068. [Link]
7606(1988)100<1054:TFAFON>[Link];2

[4] Casagli N, Ermini L (1999) Geomorphic analysis of landslide dams in the Northern
Apennine. undefined

[5] Chen Q, Zhang L (2006) Three-dimensional analysis of water infiltration into the
Gouhou rockfill dam using saturated-unsaturated seepage theory. Can Geotech J 43:449–
461. [Link]

[6] Bolton Seed H, Duncan JM (1987) The failure of Teton Dam. Eng Geol 24:173–205.
[Link]

[7] Bolton Seed H, Duncan JM (1987) The failure of Teton Dam. Eng Geol 24:173–205.
[Link]

[8] Charles JA, Skinner HD (2001) Compressibility of foundation fills. Proc Inst Civ Eng -
Geotech Eng 149:145–157. [Link]

[9] Fan X, Tang CX, Van Westen CJ, Alkema D (2012) Simulating dam-breach flood
scenarios of the Tangjiashan landslide dam induced by the Wenchuan Earthquake. Nat
Hazards Earth Syst Sci 12:3031–3044. [Link]

[10] Association of State Dam Safety. [Link]

instabilities#Slides. Accessed 15 May 2020

[11] Ohio Department of NATURAL RESOURCES (2011): Dam safety embankment

instabilities

[12] Sjah J (2013) Couplage SPH-DEM pour l’étude de l’érosion dans les ouvrages
hydrauliques. Ecole Centrale de Lyon

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 142

[13] MAREK MA (2011) Hydraulic Design Manual

[14] Wan CF, Fell R (2004) Investigation of Rate of Erosion of Soils in Embankment Dams.
J Geotech GEOENVIRONMENTAL Eng 130:373–380

[15] Zhang J, Yun LI, Guoxiang X, et al (2009) Overtopping breaching of cohesive

homogeneous earth dam with different cohesive strength. 52:.
[Link]

[16] Wan CF (2006) Experimental Investigations of Piping Erosion and Suffusion of Soils in
Embankment Dams and their Foundations

[17] Reddy KR, Richards KS (2012) Experimental investigation of initiation of backward

erosion piping in soils. Géotechnique 62:933–942. [Link]

[18] Kenney TC, Lau D (1986) Internal stability of granular filters: Reply. Can Geotech J
23:420–423

[19] Fell R, Wan CF, Foster M (2004) Methods for estimating the probability of failure of
embankment dams by internal erosion and piping : piping through the embankment.
University of New South Wales, School of Civil and Environmental Engineering,
Sydney

[20] Fell R, Fry JJ (Jean J. (2007) The state of the art of assessing the likelihood of internal
erosion of embankment dams, water retaining structures and their foundations. Taylor &
Francis, London, p 24

[21] Sherard JL (1979) Sinkholes in Dams of Coarse, Broadly Graded Soils. Proc 13th
ICOLD Congr India Vol. 2:25–35

[22] Garner S., Fannin R. (2010) Understanding internal erosion: a decade of research
following a sinkhole event. Int J HYDROPOWER DAMS 17:93–100

[23] Kézdi A (1979) Soil physics : selected topics. Elsevier Scientific Pub. Co., Amsterdam

[24] Istomina VS (1957) Filtration stability of soils. Gostroizdat, Moscow, Leningrad

[25] Burenkova V. (1992) Assessment of Suffusion in Non-Cohesive and Graded Soils. In:
Proceedings of the 1st International Conference “Geo-Filters.” Karlsruhe, pp 357–360

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 143

[26] Mao C. (2005) Study on piping and filters: Part I of piping. Rock Soil Mech 26:209–215

[27] Wan CF, Fell R (2004) Laboratory tests on the rate of piping erosion of soils in
embankment dams. Geotech Test J 27:295–303. [Link]

[28] Li M, Fannin RJ (2008) Comparison of two criteria for internal stability of granular soil.
Can Geotech J 45:1303–1309. [Link]

[29] Terzaghi K, B. Peck R, Mesri G (1996) Soil Mechanics in Engineering Practice, 3rd
Edition | Wiley. [Link]
fr/Soil+Mechanics+in+Engineering+Practice,+3rd+Edition-p-9780471086581.
Accessed 31 May 2020

[30] Sherard JL, Dunnigan LP, Talbot JR (1984) Basic Properties of Sand and Gravel Filters.
J Geotech Eng 110:684–700. [Link]
9410(1984)110:6(684)

[31] Kenney TC, Lau D (1985) Internal stability of granular filters. Can Geotech J 22:215–
225

[32] Schuler U (1995) How to deal with the problem of suffosion. Res Dev F dams, Swiss
Natl Comm Larges Dams 145–159

[33] Tomlinson SS, Vaid YP (2000) Seepage forces and confining pressure effects on piping
erosion. Can Geotech J 37:1–13. [Link]

[34] Wan CF, Fell R (2004) Investigation of Rate of Erosion of Soils in Embankment Dams.
J Geotech Geoenvironmental Eng 130:373–380. [Link]
0241(2004)130:4(373)

[35] Lim SS (2006) Experimental investigation of erosion in variably saturated clay soils. The
University of New South Wales

[36] Kakuturu S, Reddi LN (2006) Evaluation of the Parameters Influencing Self-Healing in

Earth Dams. J Geotech Geoenvironmental Eng 132:879–889.
[Link]

[37] Richards KS, Reddy KR (2007) Critical appraisal of piping phenomena in earth dams.
Bull Eng Geol Environ 66:381–402. [Link]

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 144

[38] Xiao M, Gomez J (2010) Effect of Piping on Shear Strength of Levees. In:
Geoenvironmental Engineering and Geotechnics. American Society of Civil Engineers,
Reston, VA, pp 51–56

[39] Moffat R, Fannin JR (2011) A hydromechanical relation governing internal stability of

cohesionless soil. Can Geotech J 48:413–424. [Link]

[40] Indraratna B, Nguyen VT, Rujikiatkamjorn C (2011) Assessing the Potential of Internal
Erosion and Suffusion of Granular Soils. J Geotech Geoenvironmental Eng 137:550–
554. [Link]

[41] Burenkova VV (1993) Assessment of Suffusion in Non-Cohesive and Graded Soils. In:
Brauns, J., Schuler, U. and Heibum, M., Eds., Proceedings of the 1st International
Conference “Geo-Filters”, Filters in Geotechnical Engineering, Balkema. pp 357–360

[42] Kézdi A (1979) Soil physics. Amsterdam, Netherland

[43] Wan C, Fell R (2008) Assessing the Potential of Internal Instability and Suffusion in
Embankment Dams and Their Foundations. J Geotech Geoenvironmental Eng 134:401–
407

[44] Sherman W (1953) Filter Experiments and Design Criteria. . ARMY ENGINEER
WATERWAYS EXPERIMENT STATION VICKSBURG MISS

[45] Milligan V (1986) Internal stability of granular filters: Discussion. Can Geotech J
23:414–418

[46] Sherard JL, Dunnigan LP (1986) Internal stability of granular filters: Discussion. Can
Geotech J 23:418–420. [Link]

[47] Skempton AW, Brogan JM (1994) Experiments on piping in sandy gravels.

Géotechnique 44:449–460. [Link]

[48] Kenney TC, Lau D (1986) Internal stability of granular filters. Can Geotech J 23:420–
423. [Link]

[49] Li M (2008) Seepage Induced Instability in Widely Graded Soils. 319.

[Link]

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 145

[50] Kovács G (1981) Seepage hydraulics. Elsevier Scientific Pub. Co.

[51] Li M (2008) Seepage induced instability in widely graded soils. University of British
Columbia, Vancouver, Canada

[52] Chang DS (2012) Internal erosion and overtopping erosion of earth dams and landslide
dams

[53] Den Adel H, Bakker KJ, Klein M, Breteler B (1988) Internal stability of minestone. In:
Proceedings International Symposium on Modelling Soil–Water–Structure Interaction,
International Association for Hydraulic Research (IAHR), Balkema Publishers. pp 225–
231

[54] Hopkins TC, Allen DL, Deen RC (1975) Kentucky Bureau of Highways Prepared in
cooperation with the US Department of Transp0rtation, Federal Highway Administration
Study Title: Investigation of Landslides on Highways

[55] Duncan JM (1996) State of the Art: Limit Equilibrium and Finite-Element Analysis of
Slopes. J Geotech Eng 122:577–596. [Link]
9410(1996)122:7(577)

[56] Fellenius W (1936) Calculation of the stability of earth dams. In: Trans., 2nd Congress
on Large Dams. U.S. Government Printing Office, Washington, DC, 4, pp 445–462

[57] Bailey WA, Christian JT (1969) ICES LEASE - I : a problem-oriented language for slope
stability analysis : users̓ manual - New Jersey Institute of Technology. MIT Soil
Mechanics Publ. No. 235, Massachusetts Inst. of Technology, Cambridge, MA

[58] Bishop AW (1955) The use of the slip circle in the stability analysis of slopes.
Geotechnique 5:7–17. [Link]

[59] Spencer E (1967) Embankments Assuming Parallel Inter-Slice Forces. Géotechnique

17:11–26

[60] Spencer E (1973) Thrust line criterion in embankment stability analysis. Geotechnique
23:85–100. [Link]

[61] Janbu N, Bjerrum L, Kjaernsli B (1956) Soil Mechanics Applied to some Engineering
Problems. Norwegian Geotechnical Institute, 1956

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 146

[62] Janbu N (1973) SLOPE STABILITY COMPUTATIONS. Publ Wiley Sons, Inc

[63] Lowe J, Karafiath L (1959) Stability of earth dams upon drawdown. In: 1st Pan American
Conf. on Soil Mechanics and Foundation Engineering. International Society of Soil
Mechanics and Geotechnical Engineering (ISSMGE), London, pp 537–552

[64] (USACE) USAC of E (1970) Engineering and Design Stability of Earth and Rock-fill
Dams. Engineer Manual EM 1110-2- 1902, Washington, DC

[65] Morgenstern NR, Price VE (1965) The Analysis of the Stability of General Slip Surfaces.
Géotechnique 15:79–93. [Link]

[66] Whitman R V., Bailey WA (1967) Use of Computers for Slope Stability Analysis. J Soil
Mech Found Div 93:475–498. [Link]

[67] Spencer E (1967) A Method of analysis of the Stability of Embankments Assuming

Parallel Inter-Slice Forces. Géotechnique 17:11–26.
[Link]

[68] Fredlund DG, Krahn J (1977) Comparison of Slope Stability Methods of Analysis. Can
Geotech J 14:429–439. [Link]

[69] Sterpi D (2003) Effects of the Erosion and Transport of Fine Particles due to Seepage
Flow. Int J Geomech 4:191–198. [Link]
3641(2004)4:3(191)

[70] Sibille L, Marot D, Sail Y (2015) A description of internal erosion by suffusion and
induced settlements on cohesionless granular matter. Acta Geotech 10:735–748.
[Link]

[71] Chang DS, Zhang LM (2011) A Stress-controlled Erosion Apparatus for Studying
Internal Erosion in Soils. Geotech Test J 34:579–589.
[Link]

[72] Ke L, Takahashi A (2012) Strength reduction of cohesionless soil due to internal erosion
induced by one-dimensional upward seepage flow. Soils Found.
[Link]

[73] Benamar A (2019) Suffusion evaluation of coarse-graded soils from Rhine dikes. Acta

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 147

 Geotech 1:. [Link]

[74] Yang J Analyses numériques de la problématique multi-physique des fontis au voisinage

d’une digue ou d’un ouvrage linéaire

[75] Horikoshi K, Takahashi A (2015) Suffusion-induced change in spatial distribution of

fine fractions in embankment subjected to seepage flow. Soils Found 55:1293–1304.
[Link]

[76] van Beek VM, Knoeff H, Sellmeijer H (2011) Observations on the process of backward
erosion piping in small-, medium- and full-scale experiments. Eur J Environ Civ Eng
15:1115–1137. [Link]

[77] Hu W, Hicher P-Y, Scaringi G, et al (2018) Seismic precursor to instability induced by

internal erosion in loose granular slopes. Géotechnique 68:989–1001.
[Link]

[78] Garey A Fox, Maria Chu-Agor, Glenn V Wilson (2007) Erosion of Noncohesive
Sediment by Groundwater Seepage: Lysimeter Experiments and Modeling. Soil Sci Soc
Am J 71:1822–1830. [Link]

[79] Muir Wood D, Maeda K, Nukudani E (2010) Modelling mechanical consequences of

erosion. Geotechnique 60:447–457. [Link]

[80] Cividini A, Gioda G (2004) Finite-element approach to the erosion and transport of fine
particles in granular soils. Int J Geomech 4:191–198.
[Link]

[81] Papamichos E, Stavropoulou M (1998) An erosion-mechanical model for sand

production rate prediction. Int J Rock Mech Min Sci 35:531–532.
[Link]

[82] Papamichos E, Vardoulakis I, Tronvoll J, Skjvrstein A (2001) Volumetric sand

production model and experiment. Int J Numer Anal Methods Geomech 25:789–808.
[Link]

[83] Papamichos E, Vardoulakis I (2005) Sand erosion with a porosity diffusion law. Comput
Geotech 32:47–58. [Link]

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 148

[84] Stavropoulou M, Papanastasiou P, Vardoulakis I (1998) Coupled wellbore erosion and
stability analysis. Int J Numer Anal Methods Geomech 22:749–769.
[Link]
K

[85] Vardoulakis I, Stavropoulou M, Papanastasiou P (1996) Hydro-mechanical aspects of

the sand production problem. Transp Porous Media 22:225–244

[86] Vardoulakis I, Papanastasiou P, Stavropoulou M (2001) Sand erosion in axial flow

conditions. Transp Porous Media 45:267–281.
[Link]

[87] Zhang L, Wu F, Zhang H, et al (2019) Influences of internal erosion on infiltration and

slope stability. Bull Eng Geol Environ 78:1815–1827. [Link]
017-1185-2

[88] Zhang XS, Wong H, Leo CJ, et al (2013) A Thermodynamics-Based Model on the
Internal Erosion of Earth Structures. Geotech Geol Eng 31:479–492.
[Link]

[89] Yang J, Yin ZY (2019) Internal erosion in dike‐on-foundation modeled by a coupled

hydromechanical approach. Int J Numer Anal Methods Geomech 43:663–683.
[Link]

[90] Chu-Agor ML, Wilson G V, Fox GA (2008) Numerical modeling of bank instability by
seepage erosion undercutting of layered streambanks. J Hydrol Eng 13:1133–1145.
[Link]

[91] Sjah J, Sph C, Pour DEM, et al (2013) COUPLAGE SPH - DEM POUR L’ÉTUDE DE
L’ÉROSION DANS LES OUVRAGES HYDRAULIQUES

[92] Vandamme JR Novel Particle Model for the Prediction of Stability and Episodic
Collapse of Coastal Cliffs and Levees

[93] Gu DM, Huang D (2019) A DEM-based approach for modeling the evolution process of
seepage-induced erosion in clayey sand. Acta Geotech 0:
[Link]

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 149

[94] Idelsohn SR (2006) Modeling bed erosion in free surface flows by the particle finite
element method. 237–252. [Link]

[95] Hosn RA, Nguyen CD, Sibille L, et al (2017) Microscale Analysis of the Effect of
Suffusion on Soil Mechanical Properties. [Link]

[96] Hu Z, Zhang Y, Yang Z (2019) Suffusion-induced deformation and microstructural

change of granular soils: a coupled CFD–DEM study. Acta Geotech 14:795–814.
[Link]

[97] Zhou H, Wang G, Jia C, Li C (2019) A novel, coupled CFD-DEM model for the flow
characteristics of particles inside a pipe. Water (Switzerland) 11:.
[Link]

[98] Cundall PA, Strack ODL (1979) A discrete numerical model for granular assemblies.
Géotechnique 29:47–65. [Link]

[99] Zou Y, Chen Q, Chen X, Cui P (2013) Discrete numerical modeling of particle transport
in granular filters. Comput Geotech 47:48–56.
[Link]

[100] Uzuoka R, Ichiyama T, Mori T, Kazama M (2012) Hydro-mechanical analysis of internal

erosion with mass exchange between solid and water. ICSE6 Paris 655–662

[101] Cividini A, Bonomi S, Vignati GC, Gioda G (2009) Seepage-Induced Erosion in

Granular Soil and Consequent Settlements. Int J Geomech 187–194.
[Link]

[102] COMSOL (2017) COMSOL Multiphysics User’s Guide, Version 5.3

[103] Cundall PA, Strack ODL (1979) A discrete numerical model for granular assemblies.
Geotechnique 29:47–65. [Link]

[104] Gregoretti C, Maltauro A, Lanzoni S (2010) Laboratory Experiments on the Failure of

Coarse Homogeneous Sediment Natural Dams on a Sloping Bed. J Hydraul Eng
136:868–879. [Link]

[105] Liao WM, Chou HT (2003) Debris flows generated by seepage failure of landslide dams.
77017

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 150

[106] Franca MJ, Almeida AB (2002) Experimental Tests on Rockfill Dam Breaching Process.
IAHR - Int Symp Hydraul Hydrol Asp Reliab Saf Assess Hydraul Struct

[107] Wörman A (1993) Seepage-Induced Mass Wasting in COARSE SOIL SLOPES.

119:1155–1168

[108] Cruden DM, Keegan TR, Thomson S (1993) The landslide dam on the Saddle River near
Rycroft, Alberta. Can Geotech J 30:1003–1015. [Link]

[109] Leps T (1973) Flow Through Rockfill. New York

[110] Meyer W, Schuster RL, Fellow, et al (1994) Potential for seepage erosion of landslide
dam. J Geotech Eng 120:1211–1229

[111] Ojha CSP, [Link], Singh VP, et al (2003) Determination of critical head in soil piping.
J Hydraul Eng 129:511–518. [Link]

[112] Takahashi T (2007) A Review of Japanese Debris Flow Research. Int J Eros Control Eng
2:1–14. [Link]

[113] Hutchinson JN (1968) Field meeting on the coastal landslides of Kent: 1–3 July 1966.
Proc. Geol. Assoc. 79:227–237

[114] Higgins CG, Schoner C (1997) Sinkholes formed by piping into buried channels.
Geomorphology 20:307–312. [Link]

[115] Iverson RM, Major JJ (1986) Groundwater Seepage Vectors and the Potential for
Hillslope Failure and Debris Flow Mobilization. Water Resour Res 22:1543–1548.
[Link]

[116] Jones JAA (1990) Piping effects in humid lands. Spec Pap Geol Soc Am 252:111–138.
[Link]

[117] Dunne T (1990) Hydrologyy mechanics, and geomorphic implications of erosion by

subsurface flow. Spec Pap Geol Soc Am 252:1–28. [Link]

[118] Hagerty DJ (1991) Piping/Sapping Erosion. I: Basic Considerations. J Hydraul Eng

117:991–1008. [Link]

[119] Sellmeijer JB, Koenders MA (1991) A mathematical model for piping. Top Catal

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 151

 15:646–651. [Link]

[120] Crosta G, Prisco C Di (1999) On slope instability induced by seepage erosion. Can
Geotech J 36:1056–1073. [Link]

[121] Lei L, Jin bao S (2000) Engineering behavior of gravel materials of Gouhou dam. J
NANJING Hydraul Res Inst

[122] Liu J, Dind L, Miao L, Yang K (1998) Model test for dam breach of Gouhou concrete
face sandy gravel dam. J Hydraul Eng 11:

[123] Gouhou Dam Failure Investigation Team (1996) Technical details of the Gouhou dam,
Gouhou concrete-Faced Rockfill Dam-design, Construction, Operation, and failure (ed.
by China National Flood and Drought Prevention Office). Beijing

[124] Zhang LM, Chen Q (2006) Seepage failure mechanism of the Gouhou rockfill dam
during reservoir water infiltration. Soils Found 46:557–568.
[Link]

[125] Douglas KJ, Fell R, Peirson WL, Studholme H (2019) Experimental investigation of
global backward erosion and suffusion of soils in embankment dams. 807:789–807

[126] Marot D, Rochim A, Nguyen HH, et al (2004) Assessing the susceptibility of gap-graded
soils to internal erosion: proposition of a new experimental methodology. Nat Hazards
83:365–388. [Link]

[127] Marot D, Rochim A, Nguyen HH, et al (2016) Assessing the susceptibility of gap-graded
soils to internal erosion: proposition of a new experimental methodology. Nat Hazards
83:365–388. [Link]

[128] Sharif Y a, Elkholy M, Chaudhry MH, Imran J (2015) Experimental Study on the Piping
Erosion Process in Earthen Embankments. J Hydraul Eng 141:1–9.
[Link]

[129] Ke L, Takahashi A (2012) Strength reduction of cohesionless soil due to internal erosion
induced by one-dimensional upward seepage flow. Soils Found 52:698–711.
[Link]

[130] Ke L, Takahashi A (2015) Drained Monotonic Responses of Suffusional Cohesionless

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 152

 Soils. J Geotech Geoenvironmental Eng 141:04015033.
[Link]

[131] Ouyang M, Takahashi A (2015) Influence of initial fines content on fabric of soils
subjected to internal erosion. Can Geotech J 53:299–313. [Link]
2014-0344

[132] Ke L, Takahashi A (2014) Experimental investigations on suffusion characteristics and

its mechanical consequences on saturated cohesionless soil. Soils Found 54:713–730.
[Link]

[133] Chang DS, Zhang L (2011) A stress-controlled erosion apparatus for studying internal
erosion in soils. Geotech Test J 34:103889. [Link]

[134] Bradshaw A, Baxter C (2007) Sample Preparation of Silts for Liquefaction Testing.
Geotech Test J 30:324–332. [Link]

[135] Yang ZX, Li XS, Yang J (2008) Quantifying and modelling fabric anisotropy of granular
soils. Géotechnique 58:237–248. [Link]

[136] Monnet J, Plé O, Broucke M, Nguyen D (2011) Internal Erosion Test: Design and
Preliminary Results. Ijcebm 1:

[137] Benamar A, Beaudoin A, Bennabi A, Wang H (2010) Experimental Study of Internal

Erosion of Fine Grained Soils. In: International Conference on Scour and Erosion. pp
368–377

[138] Hunter RP, Bowman ET (2018) Visualisation of seepage-induced suffusion and

suffosion within internally erodible granular media. Geotechnique 68:918–930.
[Link]

[139] Ke L, Takahashi A (2012) Strength reduction of cohesionless soil due to internal erosion
induced by one-dimensional upward seepage flow. Soils Found 52:698–711.
[Link]

[140] Ferdos F (2016) Internal erosion phenomena in embanlkment dams: Throughflow and
internal erosion mechanisms

[141] Truesdell C, Toupin R (1960) The Classical Field Theories. In: HDP. pp 226–858

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 153

[142] Bowen RM (1980) Incompressible porous media models by use of the theory of
mixtures. Int J Eng Sci 18:1129–1148. [Link]

[143] Bowen RM (1982) Compressible porous media models by use of the theory of mixtures.
Int J Eng Sci 20:697–735. [Link]

[144] De Boer R, Ehlers W (1986) On the problem of fluid- and gas-filled elasto-plastic solids.
Int J Solids Struct 22:1231–1242. [Link]

[145] Ehlers W, Bluhm J (2002) Porous Media. Springer-Verlag Berlin Heidelberg

[146] Noll W, Toupin RA, Wang C-C, Noll W (1968) A Mathematical Theory of the
Mechanical Behavior of Continuous Media. In: Continuum Theory of Inhomogeneities
in Simple Bodies. Springer Berlin Heidelberg, pp 25–54

[147] Hicher P (2012) Modeling the Behavior of Soil Subjected to Internal Erosion. ICSE6
Paris 647–654

[148] Govindaraju RS (1995) Non-dimensional analysis of a physically based rainfall-runoff-

erosion model over steep slopes. J Hydrol 173:327. [Link]
1694(95)02738-B

[149] Yang J, Yin ZY, Laouafa F, Hicher PY (2019) Modeling coupled erosion and filtration
of fine particles in granular media. Acta Geotech 14:1615–1627.
[Link]

[150] Chetti A, Benamar A, Hazzab A (2016) Modeling of Particle Migration in Porous Media:
Application to Soil Suffusion. Transp Porous Media 113:591–606.
[Link]

[151] Zhang DM, Gao CP, Yin ZY (2019) CFD-DEM modeling of seepage erosion around
shield tunnels. Tunn Undergr Sp Technol 83:60–72.
[Link]

[152] Sakthivadivel R (1967) Theory and mechanism of filtration of non-colloidal fines

through a porous medium

[153] Sakthivadivel R, Irmay S (1966) A review of filtration theories. Berkeley, Calif.:

Hydraulic Engineering Laboratory, College of Engineering, University of California,

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 154

 Berkeley.

[154] Einstein HA (1936) Der Geschiebetrieb als Wahrscheinlichkeitsproblem. ETH Zürich

[155] Gravanis E, Sarris E, Papanastasiou P (2015) Hydro-mechanical erosion models for sand
production. Int J Numer Anal Methods Geomech 39:2017–2036.
[Link]

[156] Vardoulakis I, Papamichos E (2003) A continuum theory for erosion in granular media.
Ital. Geotech. J. 3:83–90

[157] Seghir A, Benamar A, Wang H (2014) Effects of Fine Particles on the Suffusion of
Cohesionless Soils. Experiments and Modeling. Transp Porous Media 103:233–247.
[Link]

[158] Fetter CW (1993) Contaminant Hydrogeology. Macmillan Publishing Company, New

York

[159] Wan RG, Liu Y, Wang J (2006) Prediction of volumetric sand production using a
coupled geomechanics-hydrodynamic erosion model. J Can Pet Technol 45:34–40.
[Link]

[160] Nguyen DM, Plé O, Monnet J (2013) Experimental study of suffusion in granular soils.
In: Advanced Materials Research. Trans Tech Publications Ltd, pp 101–105

[161] BENDAHMANE F (2005) Influence des interactions mécaniques eau-sol sur l’érosion
interne. Université de Nantes

[162] Yang J, Yin ZY, Laouafa F, Hicher PY (2019) Internal erosion in dike-on-foundation
modeled by a coupled hydromechanical approach. Int J Numer Anal Methods Geomech
43:663–683. [Link]

[163] Petterson KE (1955) The Early History of Circular Sliding Surfaces. Géotechnique
5:275–296. [Link]

[164] Fellenius W (1936) Calculation of the Stability of Earth Dams. In: Proceedings of the
Second Congress of Large Dams. pp 445–463

[165] Janbu, N., Bjerrum, L., & Kjaernsli B (1956) Stabilitetsberegning for fyllinger

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 155

 skjaeringer og naturlige skraninger. Nor Geotech Publ 16:

[166] Bishop AW (1955) The use of the Slip Circle in the Stability Analysis of Slopes.
Géotechnique 5:7–17. [Link]

[167] GEO-SLOPE International Ltd (2012) Stability Modeling with SLOPE/W

THESE DE DOCTORAT : BOUZIANE ABDERREZAK 156

 Résumé
L’érosion interne est définie comme la migration de particules initiées par les forces d’écoulement au sein d’un ouvrage hydraulique
ou sa fondation. Ce processus peut affecter la stabilité des digues, de barrages en remblai ou les talus naturels. Les études expérimentales
présentes dans la littérature se concentrent principalement sur l’identification des gradients hydrauliques critiques pour l’initiation de l’érosion,
et sur l’analyse de la courbe granulométrique des sols érodables pour définir des critères permettant de quantifier le potentiel d’érosion d’un
sol donné. Cependant, il existe peu d’études sur l’effet de l’érosion sur le comportement mécanique des sols. Dans cette étude, on s’intéresse
dans un premier temps à l’effet de l’érosion interne sur les paramètres de résistance au cisaillement des sols sableux. Des essais de cisaillement
direct ont été effectués à l’aide d’un dispositif (boîte de Casagrande) modifié sur des sols avant et après érosion interne. Les résultats de cette
étude ont montré que l’érosion affecte considérablement la résistance au cisaillement des sols étudiés. On s’intéresse dans un deuxième temps
à la modélisation de l’érosion interne, un modèle mathématique avec une loi d’érosion modifiée, basée sur la diffusion de la porosité du milieu,
est proposé et résolu numériquement par la méthode des éléments finis. L’effet de l’érosion interne sur la stabilité au glissement a été étudié à
travers une co-simulation éléments finis-équilibre limite, en calculant la variation du facteur de sécurité en fonction du temps. Les résultats ont
montré que l’érosion interne entraine une réduction de la valeur du facteur de sécurité de 16% de sa valeur initiale, et peut conduire à la rupture
de l’ouvrage.

Mots clés : Erosion interne, Résistance au cisaillement, Stabilité, digues, Barrages en remblai, Talus naturels, Modélisation.

Abstract
Internal erosion is defined as the migration of particles initiated by flow forces within a hydraulic structure or its foundation. This
process can affect the stability of dikes, embankment dams or natural slopes. The experimental studies present in the literature focus mainly
on the identification of critical hydraulic gradients for erosion initiation, and on the analysis of the particle size distribution of erodible soils to
define criteria to quantify the erosion potential of a given soil. However, there are few studies on the effect of erosion on the mechanical
behavior of soils. In this study, the effect of internal erosion on the shear strength parameters of sandy soils is first investigated. Direct shear
tests were performed using a modified Casagrande box on soils before and after internal erosion. The results of this study showed that erosion
significantly affects the shear strength of the studied soils. In a second time, we are interested in the modelling of the internal erosion, a
mathematical model with a modified erosion law, based on the diffusion of the porosity of the medium, is proposed and solved numerically by
the finite element method. The effect of internal erosion on sliding stability was studied through a finite element-limit equilibrium co-
simulation, by calculating the variation of the safety factor over time. The results showed that internal erosion leads to a reduction of the safety
factor value by 16% of its initial value, and can lead to the failure of the structure

Key words: Internal erosion, Shear strength, Stability, Dikes, Embankment dam, Natural slopes, Modeling.

‫ملخص‬

‫ تركز الدراسات‬.‫ يمكن أن تؤثر هذه العملية على استقرار السدود أو المنحدرات الطبيعية‬.‫ُعرف التآكل الداخلي بأنه تنقل الجسيمات الناتجة عن قوى التدفق داخل هيكل هيدروليكي أو أساسه‬
َّ ‫ي‬

‫ وعلى تحليل منحنى حجم الجسيمات للتربة القابلة للتآكل لتحديد المعايير التي‬، ‫التجريبية المقترحة من قبل الباحثين بشكل أساسي على تحديد التدرجات الهيدروليكي ة الحرجة لبدء التآكل‬

‫ في هذه الدراسة تم دراسة تأثير التآكل الداخلي‬.‫ هناك القليل من الدراسات حول تأثير التآكل الداخلي على السلوك الميكانيكي للتربة‬، ‫ ومع ذلك‬.‫تسمح بتحديد إمكانية التآكل لتربة معينة‬

‫ أظهرت نتائج هذه الدراسة أن‬.‫ تم إجراء اختبارات القص المباشر باستخدام علبة كازاغروند المعدلة على التربة قبل وبعد التآكل الداخلي‬.‫على تغيرات مقاومة القص للتربة الرملية‬

‫ بنا ًء على‬، ‫ صببنا اهتمامنا الى نمذجة التآكل الداخلي و قد تم اقتراح نموذج رياضي مع قانون تآكل معدل‬، ‫ من جهة أخرى‬.‫االنجراف يؤثر بشكل كبير على مقاومة القص للتربة المدروسة‬

، ‫ تمت دراسة تأثير التآكل الداخلي على ثبات االمنحدرات من خالل محاكاة مشتركة بين حد التوازن و العناصر المحدود‬.‫ وحلها عدديًا بطريقة العناصر المحدودة‬، ‫انتشار مسامية الوسط‬

‫ ويمكن أن يؤدي إلى انهيار‬، ‫ من قيمته األولية‬٪16 ‫ أظهرت النتائج أن التآكل الداخلي يؤدي إلى انخفاض قيمة عامل األمان بنسبة‬.‫عن طريق حساب التباين في عامل األمان بمرور الوقت‬

.‫الهيكل‬

.‫ النمذجة‬، ‫ المنحدرات الطبيعية‬، ‫ السدود الترابية‬، ‫ الثبات‬، ‫ مقاومة القص‬، ‫ التآكل الداخلي‬:‫الكلمات المفتاحية‬

View publication stats