Memoire 10

RÉPUBLIQUE DU BÉNIN
**********
MINISTÈRE DE L’ENSEIGNEMENT SUPERIEUR ET DE LA
RECHERCHE SCIENTIFIQUE
**********
UNIVERSITÉ NATIONALE DES SCIENCES, TECHNOLOGIES,
INGENIERIES ET MATHEMATIQUES (UNSTIM) D’ABOMEY
**********
ECOLE NATIONALE SUPERIEURE DES TRAVAUX PUBLICS (ENSTP-ABOMEY)
**********
Département : Génie Civil
MEMOIRE DE FIN DE FORMATION
POUR L’OBTENTION DU DIPLOME D’INGENIEUR DE CONCEPTION
Option : Ingénierie Routière
THEME
Influence des cycles de retrait-gonflement des

sols expansifs sur le comportement des
chaussées : cas du tronçon Pobè-Kétou
Présenté par :
Mr AGASSOUSSI Metognissè Ludovic
Encadré par :
Dr. Ing. AHLINHAN Marx Ferdinand
Maître de conférence des Universités Soutenu publiquement le 16 décembre

CAMES ENSTP/UNSTIM 2022 devant un jury composé de :
Dr. Ing ADEOTI Oyéniran Guy
Maître-Assistant des Universités

CAMES ENSTP/UNSTIM Président : Dr GBAGUIDI Hubert Frédéric
HOUEMAVO Donald Membres

Ingénieur de conception en Génie Civil 1. Dr. Ing ADEOTI Oyéniran Guy
2. Dr Ing CHANHOUN Marius
Sous la supervision de :
3. Ing RABENOMANANA ZO
Dr. GBAGUIDI Hubert Frédéric
Maître de conférence des Universités

CAMES ENSTP/UNSTIM
Année Académique : 2021-2022
2ème Promotion
Sommaire
Sommaire............................................................................................................................................................................................. i
Dédicaces ........................................................................................................................................................................................... ii
Remerciements ................................................................................................................................................................................iii
Liste des sigles et abréviations .................................................................................................................................................. v
Liste des figures ............................................................................................................................................................................ vii
Liste des tableaux ........................................................................................................................................................................viii
Liste des photos ............................................................................................................................................................................. ix
Résumé ................................................................................................................................................................................................x
Abstract ........................................................................................................................................................................................... xi
Introduction ...................................................................................................................................................................................... 1
Chapitre I : Généralités sur les argiles et le phénomène de retrait-gonflement. ......................................................7
I.1 Les argiles ..............................................................................................................................................................................8
I.2 Généralités sur les chaussées revêtues .................................................................................................................... 20
I.3 Comportement et pathologie des structures de chaussées sur sols gonflants ..............................................21
Chapitre II : Matériels, matériaux et méthodologie ........................................................................................................ 24
II.1 Milieu d’étude ................................................................................................................................................................... 25
II.2 Matériels d’étude ........................................................................................................................................................... 27
II.3 Méthodologie de recherche ........................................................................................................................................ 29
II.4 Description du mode opératoire de l’humidification et de séchage ............................................................... 38
Chapitre III : Caractérisation au laboratoire des sols naturels ................................................................................... 42
III.1 Essais d’identification ................................................................................................................................................. 43
III.2 Essai mécanique ............................................................................................................................................................ 48
III.3 Estimation du gonflement par la méthode indirecte [21] .................................................................................51
III.4 Présentation des résultats des essais sur les échantillons ayant subi le cycle d’humidification et de
séchage ....................................................................................................................................................................................... 53
Chapitre IV : Modélisation et prédiction du comportement des sols expansifs sous léffet des sollicitations
cycliques : humidification et séchage. .................................................................................................................................... 60
IV.1 Historique de la structure ............................................................................................................................................61
IV.2 Présentation du logiciel Plaxis ................................................................................................................................... 62
IV.3 La modélisation numérique de la chaussée ............................................................................................................. 76
IV.4 Les résultats de calcul ................................................................................................................................................. 86
CONCLUSION ET PERSPECTIVES ........................................................................................................................................ 95
RÉFÉRENCES BIBLIOGRAPHIQUES .................................................................................................................................... 98
Annexes .......................................................................................................................................................................................... 101
Table des Matières ..................................................................................................................................................................... 111
i
Influence des cycles de retrait-gonflement des sols expansifs sur le comportement des
Dédicaces
A:
 Mon Père AGASSOUSSI Ferdinand et ma Mère ANAGONOU
Clémentine à qui je dois l’éducation et le désir de toujours faire
mieux. Puisse Dieu établir son règne dans vos vies et vous combler.
 Ma sœur jumelle AGASSOUSSI Lurindic, ce travail est le tien
aussi.
 Ma famille
 Ma chère amie AïGO Blandine
LUDOVIC METOGNISSE AGASSOUSSI II

Remerciements
Merci à toi mon Dieu pour ses bienfaits dans ma vie. Gloire au tout puissant qui a veillé
sur moi jusqu’à l’aboutissement de ce travail. Adoration et louange à toi qui relève les
pauvres.
La réalisation de ce mémoire a été possible grâce au concours de plusieurs personnes à

qui je voudrais témoigner toute ma gratitude.
Une note spéciale de gratitude revient à mes encadreurs, Dr. Ing. Marx Ferdinand
AHLINHAN, Maitre de conférences des Universités du CAMES et Dr. Ing. Guy
Oyéniran ADEOTI, Maître Assistant des Universités du CAMES pour s’être rendus aussi
disponibles qu’accessibles en m’orientant dans l’accomplissement de ce travail. Puisse le
Très Haut vous accorder une longévité, vous honorer et vous combler de grâces dans
toutes vos entreprises. Merci également à mon superviseur Dr GBAGUIDI Hubert
Frédéric, Maitre de conférences des Universités du CAMES pour avoir accepté
superviser ce mémoire.
J’adresse également toute ma gratitude à l’endroit de Mr HOUEMAVO Donald, Ingénieur

de conception en Génie Civil pour son soutien et sa disponibilité. Que le seigneur bénisse
vos projets.
Je tiens à remercier :
Le Directeur de l’ENSTP Dr ALAMOU Adéchina Éric, le Directeur Adjoint de

l’ENSTP Dr AROUNA Ousséni et tous les membres de l’administration de l’ENSTP ;
Le Chef Service des Stages et Insertion Professionnelle, Ing Chéissou KOTO
TAMOU et Dr Peace HOUNKPE, chef département de Génie Civil à l’ENSTP pour
les démarches menées auprès de l’entreprise ADEOTI SA dans le but de nous
trouver des structures d’accueil pour nos stages académiques de fin de formation ;
Le Directeur technique du Centre National d'Etudes et de Recherches en Travaux
Publics (CNERTP), Monsieur Orphée DANSOU de nous avoir permis de faire
quelques mois de stage au sein du Centre National d’Essais et de Recherche des
Travaux Publics (CNERTP) afin de réaliser nos essais ;
Tout le personnel du CNERTP et particulièrement à Monsieur Patrick
KINNINSSODE, Chef Section essais mécaniques ;
Tout le personnel de ADEOTI SA et particulièrement au directeur des travaux
adjoint sur le projet de Construction de la route Guessou sud-Sinendé;
Toute ma gratitude à l’endroit de tout le corps professoral de l’ENSTP qui n’a ménagé
aucun effort pour nous donner le savoir auquel nous nous accrochons aujourd’hui pour faire
valoir notre formation. Il s’agit de :
Pr Edmond ADJOVI, Professeur Titulaire des Universités du CAMES ;

Pr Mohamed GIBIGAYE, Professeur Titulaire des Universités du CAMES ;
LUDOVIC METOGNISSE AGASSOUSSI III

Pr Aristide HOUNGAN, Professeur Titulaire des Universités du CAMES ;

Dr Eric ALAMOU, Professeur Titulaire des Universités du CAMES ;
Dr Adolphe TCHEHOUALI, Maître de Conférences des Universités du CAMES ;
Dr Ousséni AROUNA
Dr Yvette KIKI TANKPINOU, Maître de Conférences des Universités du CAMES;
Dr Peace HOUNKPE, Maître de Conférences des Universités du CAMES;
Dr Marx AHLINHAN, Maître de Conférences des Universités du CAMES ;
Dr Hubert Frédéric GBAGUIDI, Maître de Conférences des Universités du
CAMES ;
Dr Guy ADEOTI, Maître Assistant des Universités du CAMES ;
Dr Prudence Crespin YABI, Maître Assistant des Universités du CAMES ;
Dr Daniel AGOSSOU, Docteur en Génie Civil.
Je pense enfin à :
Mon grand-père Joseph ZOHA pour ses conseils et son aide. Que le Seigneur vous
accorde ces grâces ;
Mes tantes et mes oncles ;
Aux familles AGASSOUSSI et ZOHA, mes frères, sœurs et amis ;
Mes amis et tous ceux qui ont été d’une aide capitale dans mon cursus ;
A tous les camarades de la deuxième promotion de l’INSPEI et de l’ENSTP.
Vous qui m’avez toujours encouragé et soutenu tout au long de ma formation. Plein
d’amour et d’attention à mon égard, recevez ici ma sincère gratitude pour vos soutiens
matériels et moraux. Puisse Dieu établir sa seigneurie dans vos vies et vous comble de
ses bienfaits.
LUDOVIC METOGNISSE AGASSOUSSI IV

Liste des sigles et abréviations

PM : Puits manuel
HRB : Highway Research Board
GTR : Guide de terrassement routier
WL : Limite de Liquidité
IP : Indice de Plasticité
MO : Matière Organique γs : Poids Spécifique
γh : Poids humide
γsat : Poids volumique saturé
γd : Poids volumique sec CU : Cohésion non drainé
φU : Angle de flottement non drainé
eo : Indice des vides du sol en place
ei : Indice des vides initial
Cc : Indice de compression
Cs : Indice de gonflement
Cv : Coefficient de consolidation verticale
ơ'vo : Contrainte effective verticale
ơ'p : Contrainte de préconsolidation
σv : Contrainte verticale
Ψ : Angle de dilatance
ĸ * : Indice de gonflement
λ * : Indice de compression
µ * : Indice de fluage
𝜈𝑢𝑟 : Coefficient de poisson en charge-décharge (par défaut 0,15)
Ko = σ΄xx /σ΄
M : pente de la ligne de l'état critique
H : Hauteur
D : Diamètre
ΔH : Tassement
LUDOVIC METOGNISSE AGASSOUSSI V

W : Teneur en eau
% : Pourcentage
kN : kiloNewton
S : Surface
kPa : kiloPascal
NF P : Norme Française
EN : Norme Européenne
LUDOVIC METOGNISSE AGASSOUSSI VI

Liste des figures

Figure I-1: Couches tétraédriques de slice source :[9] .......................................................................... 10
Figure I-2: Couches octaédriques centré sur l’aluminium source : [9] ................................................ 10
Figure I-3: : Représentation schématique de quelques groupes de minéraux argileux TO (1:1) [1] 11
Figure I-4: Différents niveaux structuraux de sols argileux [9] ........................................................... 11
Figure I-5: Vue isométrique de la structure atomique de la kaolinite [13] ......................................... 12
Figure I-6: Structure élémentaire du réseau cristallin de la kaolinite [28] ....................................... 12
Figure I-7: Vue isométrique de la structure atomique de la montmorillonite [13] ............................ 13
Figure I-8: Structure élémentaire du réseau cristallin de la montmorillonite [28] .......................... 13
Figure I-9: Vue isométrique de la structure atomique et schématique de l’illite [13] ...................... 14
Figure I-10: Structure élémentaire du réseau cristallin de l’illite ........................................................ 14
Figure I-11: Carte de la dépression médiane [29] .................................................................................... 19
Figure I-12:Structure type d’une chaussée (SETRA & LCPC, 1994)[1] ................................................20
Figure I-13: Profil théorique initial [16] ......................................................................................................22
Figure I-14:Période de sécheresse : tassement et fissuration des zones asséchées. Le sol
durcit.[16] ...........................................................................................................................................................22
Figure I-15: Période de pluies : gonflement, disparition des fissures, le sol ramollit, la circulation
crée des ornières et des fissures dans la chaussée[16] ..........................................................................22
Figure I-16: Période de sécheresse : tassement et fissuration des zones asséchées. L’état de la
chaussée empire.[16] ........................................................................................................................................22
Figure II-1: Carte Pédologique de Pobè .......................................................................................................27
Figure II-2: Diagramme de plasticité .......................................................................................................... 31
Figure II-3: Diagramme de consistance ......................................................................................................32
Figure II-4: Classification GTR .....................................................................................................................35
Figure IV-1 : Modèle monodimensionnel du comportement élastoplastique. .......................................66
Figure IV-2: Représentation du comportement élastique-parfaitement plastique. ..........................67
Figure IV-3: Interface de saisie des paramètres du sol dans le modèle élastique linéaire ............68
Figure IV-4:Courbe intrinsèque du modèle de Mohr-Coulomb................................................................69
Figure IV-5: Interface de saisie des paramètres du sol dans le modèle Mohr-Coulomb .................69
Figure IV-6:Interface de saisie des paramètres du sol dans le modèle Soft Soil Model ...............70
Figure IV-7: Interface de saisie des paramètres du sol dans le modèle Hardening Soil ................ 71
Figure IV-8: Définition du module initiale Eo et du module à 50% de la rupture E 50. ......................72
Figure IV-9:Interface de saisie des paramètres avancé du sol ............................................................75
Figure IV-10:Modèle géométrique de la structure ...................................................................................76
Figure IV-11:Maillage du modèle ................................................................................................................... 81
Figure IV-12:Maillage du modèle ................................................................................................................... 81
Figure IV-13 : Génération des pressions interstitielles initiales ...........................................................82
Figure IV-14: Conditions hydraulique initiale .............................................................................................82
Figure IV-15 Contrainte effective ...............................................................................................................83
Figure IV-16:Désactivation de la charge .....................................................................................................84
Figure IV-17: Activation des couches de chaussée...................................................................................85
Figure IV-18: Activation de la charge ..........................................................................................................85
Figure IV-19: Déplacement total ...................................................................................................................87
Figure IV-20: Déformation de maillage .......................................................................................................87
Figure IV-21:Déplacement horizontal total ................................................................................................87
LUDOVIC METOGNISSE AGASSOUSSI VII

Figure IV-22: Contrainte totale ...................................................................................................................87

Figure IV-23: Déformation de maillage .......................................................................................................89
Figure IV-24 : Déplacement total .................................................................................................................89
Figure IV-25: Déplacement horizontal total ..............................................................................................89
Figure IV-26:Contrainte totale .....................................................................................................................89
Figure IV-27 : Déplacement total ................................................................................................................. 91
Figure IV-28 : Déformation de maillage ...................................................................................................... 91
Figure IV-29: Contrainte totale .................................................................................................................... 91
Figure IV-30 : Déplacement horizontal total ............................................................................................. 91
Liste des tableaux

Tableau I-1: Caractéristiques microscopiques des principales argiles [14] ........................................ 15
Tableau I-2: Différents types de chaussées selon la classification LCPC-SETRA (SETRA & LCPC,
1994) ....................................................................................................................................................................20
Tableau II-1:Localisation des sites de prélèvements ...............................................................................27
Tableau II-2:Potentiel de gonflement [21] ................................................................................................. 31
Tableau II-3:Désignation géotechnique des sols en fonction de la teneur en matières
organiques[9] .....................................................................................................................................................34
Tableau II-4: Classification H.R.B ................................................................................................................35
Tableau II-5: Nature du sol compressible en fonction de l’indice de compression Cc .....................37
Tableau II-6: Nature des sols en fonction de l’indice de compression Cc ..........................................37
Tableau II-7:Nature des sols gonflants en fonction de l’indice de gonflement Cs ...........................37
Tableau II-8: Tableau récapitulatif du nombre de carotte pour l'humidification ............................38
Tableau III-1: Les résultats des essais d’identification .........................................................................43
Tableau III-2:Récapitulatif des résultats qui émanent de l’essai de cisaillement ............................50
Tableau III-3:Récapitulatif des valeurs d’essai œdométrique de l’argile d’Issaba .......................... 51
Tableau III-4: Valeurs de l’amplitude et de la pression de gonflement ..............................................52
Tableau III-5: Tableau récapitulatif des essais de cisaillement ..........................................................56
Tableau IV-1: Valeurs des paramètres de compressibilité et de gonflement λ et k ........................75
Tableau IV-2:Paramètres caractérisant l’argile gonflant utilisée pour le modèle Mohr-coulomb
non drainé ...........................................................................................................................................................78
Tableau IV-3: Caractéristiques de graves latéritiques traitées au ciment pour Plaxis ...................80
Tableau IV-4:Caractéristiques de graves non traitées pour Plaxis ......................................................80
Tableau IV-5: Caractéristiques du béton bitumineux utilisé en couche de base ...............................80
Tableau IV-6: Valeurs de Ko pour quelques sols[27] ................................................................................82
Tableau IV-7: Récapitulatif des sollicitations de la simulation de la structure de chaussée .........94
LUDOVIC METOGNISSE AGASSOUSSI VIII

Liste des photos

Photo I-1: Exemple N°2 de fissure de retrait[16] .................................................................... 23
Photo I-2: Exemple N°1 d'une fissure de retrait [16] ............................................................. 23
Photo I-3:Affaissement longitudinal de la route Kétou-Pobè dans la zone de la Lama,
Bénin[9] ............................................................................................................................................... 23
Photo I-4:Fissuration au KP 39+000 de la route Kétou-Pobè, dépression de la Lama,
Bénin ..................................................................................................................................................... 23
Photo 5: Etuve .................................................................................................................................... 28
Photo 6: Bocal à humidification ...................................................................................................... 28
Photo 7: Bâtis œdométrique ........................................................................................................... 28
Photo 8: Appareil de cisaillement .................................................................................................. 28
Photo 9: Bague pour l'essai œdométrique ................................................................................... 28
Photo 10: bague rectangulaire pour l'essai de cisaillement ..................................................... 28
Photo II-6:Bague de carottage pour l'essai œdométrique ..................................................... 38
Photo II-7:Bague de carottage pour l'essai de cisaillement .................................................. 38
Photo II-8: Bac pour l'humidification .......................................................................................... 39
Photo II-9 :Support des bagues .................................................................................................... 39
Photo III-2: échantillons du cycle 0 ............................................................................................ 53
Photo III-3: échantillons du cycle 1 a) et b) .............................................................................. 54
Photo III-4 : Les échantillons du cycle 2 a) et b) ..................................................................... 54
LUDOVIC METOGNISSE AGASSOUSSI IX

Résumé
La zone sud du Bénin est traversée par la dépression de la Lama qui est fortement
constituée de sols hydromorphes et plus spécifiquement gonflant. Ces sols expansifs sont
traversés par plusieurs axes routiers dont celui du tronçon Pobè-Kétou. Malgré les
techniques envisagées comme solution pour permettre la construction des chaussées sur
ces types de sols, ces derniers continuent de causer des dégradations précoces sur les
chaussées construites. Le présent travail a étudié l’influence du cycle de retrait-
gonflement des sols expansifs sur le comportement des chaussées: cas du tronçon Pobè-
Kétou. Le but de ce travail est d’étudier la variation des paramètres mécaniques du sol au
fur et à mesure qu’il subit le cycle de retrait-gonflement dû à l'alternance des saisons
humides et sèches afin d’évaluer les conséquences de cette variation sur le comportement
des chaussées. Après une étude bibliographique sur le phénomène de retrait-gonflement
et sur les pathologies observées sur les bâtiments et les routes érigés sur sols gonflants;
nous avons effectué des descentes dans la zone d’étude qui est la commune de Pobè située
dans la dépression de la Lama afin d’appréhender l’état actuel des constructions.
L’abondance des pathologies recensées, nous ont permis de choisir au mieux les zones de
prélèvement de nos échantillons. Les sites retenus se trouvent ainsi dans l’arrondissement
d’Issaba. Après avoir confirmé que ces sols prélevés sont de l’argile plus précisément de
l’argile très plastique, surconsolidée, très compressible avec un pouvoir de gonflement
élevé, nous avons réalisé les essais de cisaillement et œdométrique sur des échantillons
de sol argileux qui ont subi des cycles de séchage-humidification. Ces essais nous ont
permis d'analyser l'influence des cycles de séchage-mouillage sur le comportement des
paramètres du sol. Cette étude révèle que les cycles de séchage-mouillage ont une
certaine influence sur les sols argileux. La structure de chaussée du tronçon Pobè-Kétou
précisément du PK 2+480 au PK 18+860 a été modélisée par le logiciel Plaxis pour une
détermination de déplacement. Les résultats des déplacements obtenus après calcul nous
ont montré que les déformations dans le corps de chaussée augmentent lorsque le nombre
de cycle augmente.
Mots clés : sol argileux, cycle, retrait, gonflement, déformation.
LUDOVIC METOGNISSE AGASSOUSSI X

Abstract
The southern zone of Benin is crossed by the Lama depression which is strongly made up
of hydromorphic and more specifically expansive soils. These expansive soils are crossed
by several roads including that of the Pobè-Kétou section. Despite the techniques
envisaged as a solution to allow the construction of pavements on these types of soil, the
latter continue to cause early deterioration of the pavements built. This work has studied
the influence of the phenomenon of shrinkage-swelling of expansive soils on the behavior
of pavements: case of the Pobè-Kétou section. The purpose of this work is to study the
variation of the mechanical parameters of the soil as it undergoes the phenomenon of
shrinkage and swelling due to the alternation of wet and dry seasons in order to assess
the consequences of this variation on the behavior of the Pobè-Kétou roadway. After a
bibliographic study on the shrinkage-swelling phenomenon and on the pathologies
observed on buildings and roads erected on swelling soils; we made descents in the study
area which is the municipality of Pobè located in the depression of the Lama in order to
understand the current state of the constructions. The abundance of identified
pathologies allowed us to choose the best areas for sampling our samples. The selected
sites are thus located in the districts of Issaba. After confirming that these sampled
soils are clay, more precisely very plastic, over consolidated, very compressible clay with
a high swelling power, we carried out the shear and oedometric tests on these clay soil
samples having undergone drying-wetting cycles. These tests allowed us to analyze the
influence of drying-wetting cycles on the behavior of soil parameters. This study reveals
that drying-wetting cycles have some influence on clay soils. The pavement structure of
the Pobè-Kétou section precisely on PK 2+480 to PK 18+860 was modeled by the Plaxis
software for a displacement determination. The results of the deformations obtained
after calculation show us that the deformations in the pavement body increase when the
number of cycles increases.
Keywords: clay soil, cycle, shrinkage, swelling, deformation.
LUDOVIC METOGNISSE AGASSOUSSI XI

Introduction
1
Introduction
L’innovation technologique apporte des bienfaits considérables à un nombre incalculable
de gens, et promet plus encore à l’avenir. Aujourd’hui, grâce à un téléphone mobile, un
agriculteur peut connaitre la demande et le coût de ses cultures. Cependant, en absence
d’une route pour transporter sa production jusqu’au marché, ces informations ne sont
guère utiles. Les routes ouvrent la voie aux transformations que permet la technologie.
Pour qu’une économie se développe, que de la richesse soit créée et que la prospérité soit
partagée, il faut que les personnes et les biens puissent circuler. Elle constitue à ce titre
un facteur important pour l’économie des pays surtout ceux en voie de développement
comme le Bénin.
Par ailleurs, si les routes exercent un impact bénéfique sur le plan socio-économique, il
faudra assurer leur stabilité, leur durabilité et leur bonne tenue au cours de leur
construction. Cette bonne tenue de la route dépend de la bonne tenue de la structure de
chaussée. La chaussée est sur le plan structural, une succession de couches empilées les
unes sur les autres pour résister au trafic. Ainsi il est important de maîtriser la
composition, le dimensionnement de cette structure afin de permettre à la route de jouer
son rôle prépondérant dans le développement.
Les principales sollicitations auxquelles doivent résister le sol sont donc celles du trafic
et du poids propre de la structure routière. Les dégradations précoces des chaussées
sont généralement dues aux surcharges du trafic ou à la faible portance du sol [1].
Les routes sont construites sur différents types de sols entre autres sur des sols argileux
qui subissent des variations de volume importantes quand leur teneur en eau varie : quand
on humidifie un tel sol à partir d’un état initial sec et dense, son volume augmente et l’on
parle de gonflement. Quand on assèche le même sol à partir d’un état initial plus lâche et
saturé, son volume diminue et l’on parle de retrait. Ces deux mécanismes (gonflement-
retrait) caractérisent les sols dits gonflants.
Les sols gonflants posent des problèmes aux constructions sur tous les continents et ce
sont pour l’essentiel des sols argileux. Pour ces sols, les variations de volume sont liées à
la présence de particules argileuses [1].
Certaines argiles contenues dans les sols naturels, les marnes ou les roches argileuses se
gonflent en présence d’eau et se rétractent de manière importante lorsqu’elles perdent
de l’eau. Il s’en suit des mouvements de terrains non uniformes pouvant aller jusqu’à
provoquer la fissuration de certaines habitations ou certaines routes. C’est ce que l’on
appelle le phénomène de retrait-gonflement des argiles [1].
LUDOVIC METOGNISSE AGASSOUSSI 2

Ces phénomènes cycliques de retrait-gonflement provoquent des contraintes parasites

dans les structures provenant des fondations. Ces contraintes engendrent des dommages
dans les bâtiments, les autoroutes, les pavages, les pistes d’aéroport et dans les
structures légèrement chargées [2].
Les méfaits des sols à supports instables, notamment ceux qui sont susceptibles de
gonflement, sont indiscutables. Les infrastructures routières et autres ouvrages
subissent des désordres, plus ou moins aggravés par les aléas climatiques de leurs sites
d’implantation. Ces dégâts handicapent le pays en général et particulièrement la région
concernée, tant sur le plan social qu’économique, dans la mesure où les infrastructures
constituent des facteurs importants de développement [1]. En république du Benin,
nombreuses sont les régions qui connaissent particulièrement les dommages engendrés
par le phénomène de retrait-gonflement des argiles gonflantes.
La route Pobè-Kétou (42,8 km) qui fait partie de la RN3 qui relie la côtière (Frontière
Togo-Cotonou-Frontière Nigeria) à la transversale Frontière Togo-Abomey-Kétou-Illara-
Frontière Nigeria, traverse cette zone de sols gonflants. Cette route à impact sous
régional, contribue au désenclavement des localités et à l’accroissement des échanges
entre le Bénin et le Nigeria. Elle améliore l’accessibilité de la zone du projet en favorisant
la libre circulation des personnes et des biens dans de meilleures conditions de sécurité.
Elle accroit également la production agricole, améliore la compétitivité du secteur routier
et stimule l’utilisation des ressources locales, humaines et matérielles.
Par ailleurs des dégradations précoces et très poussées sont observées sur cette route
une année après sa construction.
Ces dégradations peuvent êtres dues au trafic mais surtout à la portance et aux variations
de volume observées au cours des phénomènes de gonflements suite à la saturation du sol
et de retraits enregistrés après assèchement du sol.
Ainsi, il est nécessaire de trouver une solution aux problèmes que pose cette route afin
d’éviter les dégradations précoces. Vu qu’il est impossible de contourner ces sols, il
importe de les étudier afin de trouver une solution durable aux problèmes que posent
ceux-ci.
Face à cette situation de retrait-gonflement, nombreux sont les chercheurs qui ont initié
des recherches sur ce phénomène.
TANKPINOU KIKI YVETTE (2016) dans sa thèse en 2016 a fait la caractérisation à

l’échelle macroscopique et microscopique des sols argileux de quatre régions à savoir :
Adjaïgbonou, Lokossa, Ouèdèmè et Kpinnou. Son étude révèle l’insuffisance des méthodes
de caractérisation géotechnique, en particulier dans la recherche de leur potentiel

gonflant par les méthodes indirectes de caractérisation basées sur les paramètres
d’identification géotechnique à savoir la granulométrie, les limites d’Atterberg et la valeur
au bleu de méthylène [3].
Ghislain HOUNLELOU (2018) pour son mémoire de fin de formation en 2019 a travaillé
sur les sols compressibles dans la dépression de la Lama. Il a proposé une méthode de
mise en œuvre des remblais sur sol compressible [4].
Isidore DHOSSOUVI (2017) pour sa soutenance de fin de formation a travaillé sur les
sols gonflants de la commune de HOUEYOGBE dans la dépression de la LAMA. Des
résultats des essais d’identification effectués sur les échantillons intacts prélevés sur
09 sites ont montré que les propriétés des sols varient d’un sol à un autre et qu’on note la
présence des sols très plastique qui sont des argiles ou limon argilo-sableux. Cette étude
lui a permis de faire une base de donnée des différents paramètres physique et mécanique
des sols [5].
AGBELELE et al (2020) ont prédit la contrainte de rupture des sols argileux de TCHI
dans la dépression de la LAMA au Benin. Cette étude a permis l’identification de modèles
prédictifs de contrainte de rupture en fonction de paramètres physiques. Ils ont montré
que le modèle unifié basé sur le modèle de Larson Rolf, (1980) peut être utilisé pour
prédire la contrainte de rupture des sols de la dépression TCHI [6].
Zodier Virtus TOHOUNGBA et al (2020) ont travaillé sur l’influence du gonflement sur
le tassement des sols argileux de la commune de HOUEYOGBE de la dépression de la
LAMA. Des essais de gonflement libre à l’œdomètre ont été effectués sur les échantillons
pour déterminer le tassement et le gonflement libre. Il ressort que les sols étudiés de la
commune de HOUEYOGBE sont des argiles plastiques. Ils ont retenu que le gonflement
libre des sols a très peu d’influence sur la valeur du tassement final à l’œdomètre [8].
Donald HOUEMAVO (2020) a fait une étude exploratoire du comportement des

fondations superficielles sur sols gonflants dans la zone d’Issaba. Pour parvenir il a
d’abord fait la caractérisation des argiles d’Issaba et d’Ahoyèyè et a constaté que les
sols d’Issaba et d’Ahoyèyè sont des sols argileux plus précisément de l’illite. Ses deux
études réalisées avec les logiciels Plaxis et Géofond montrent que pour une même charge
de 200 kN centrée et non inclinée à l’état limite ultime et pour également une même
fondation de section 0,8x1,2 m², la capacité portante calculée avec Géofond n’est pas
vérifiée à 2m alors qu’avec le logiciel Plaxis, est vérifiée [9].
GBAFFONOU et al (2021) dans son article présentent les résultats d’une étude
expérimentale en labo réalisée à l’œdomètre et à la presse Richards sur un échantillon de
sol argileux de classe A3 ayant subi des cycles de séchage à 105°C et mouillage. Ces essais
les ont permis d'analyser l'influence des cycles de séchage-mouillage sur la

compressibilité et le comportement en succion du sol. Cette étude révèle que les cycles
de séchage-mouillage ont une certaine influence sur le comportement des sols argileux.
En effet, l'indice de compressibilité du sol augmente tandis que l'indice de gonflement
diminue selon les cycles ; la pression de préconsolidation et le module œdométrique
présentent une variation plus ou moins constante ; la perméabilité du sol augmente
nettement lorsque le nombre de cycles augmente. La succion du sol diminue à mesure que
le nombre de cycles de séchage-mouillage augmente dans des conditions extrêmes [10].
Maryam E. HAZOUME (2021) dans son travail a analysé les performances des
techniques de construction de route les plus utilisées sur sols gonflants au Bénin afin de
dégager les techniques ayant limité les dégradations précoces sur les chaussées et
pouvant assurer la durabilité des différentes structures de chaussée construites sur ces
sols. De son travail il ressort que la géogrille mise en place sur la section allant du PK
2+480 au PK 10+000 du tronçon Pobè-Onigbolo a joué un rôle de renforcement dans la
chaussée [11].
Les études antérieures réalisées sur les sols gonflants n’ont toujours pas permis de mieux
appréhender le comportement des chaussées sur ces sols. Ainsi dans le but de trouver les
causes des dégradations observées sur les chaussées construites sur sols gonflants, nous
avons décidé d’étudier le comportement des chaussées construites sur les sols gonflants
du Bénin soumis aux cycles de retrait-gonflement.
D’où l’importance de cette étude intitulée << Influence des cycles de retrait-
gonflement des sols expansifs sur le comportement des chaussées : cas du tronçon
Pobè-Kétou.>>

Objectifs
Objectif général
Lóbjectif général de cette étude est de contribuer à la compréhension de línfluence

des cycles de retrait-gonflement sur le comportement des sols expansifs dans la
dépression de la Lama, de quantifier cette influence afin de mieux cerner l´état limite
ultime et l´état limite de service voire la pérennisation des chaussées dans cette
dépression face aux cycles de saison pluvieuse (humidification, gonflement) et de saison
sèche (dessèchement et retrait).
Objectifs spécifiques
 Caractériser les échantillons pour confirmer ou infirmer síl ságit dárgile

gonflante
 Déterminer les caractéristiques des matériaux soumis aux cycles d’humidification
et séchage
 Faire la modélisation numérique d’une chaussée reposant sur une argile gonflante
par la méthode des éléments fins (Plaxis)
 Prédire le comportement des sols expansifs sous léffet des sollicitations
cycliques dues au gonflement (l´humidification pendant la saison pluvieuse et au
retrait (dessèchement pendant la saison sèche).

Chapitre I : Généralités sur les argiles et le

phénomène de retrait-gonflement.

Chapitre I : Généralités sur les argiles et le phénomène de retrait-

gonflement
I.1 Les argiles
Introduction
La variation du volume de certains sols notamment ceux de la dépression de la LAMA
constitue un phénomène très important en géotechnique car elle est à l’origine de
nombreux dommages pour les ouvrages. Ce phénomène dépend des caractéristiques
minéralogiques des minéraux argileux et se déclenche lorsque la teneur en eau du matériau
varie.
En saison pluvieuse lorsque l’eau pénètre dans ces sols, ils augmentent de volume : c’est le
phénomène de gonflement. Mais en saison sèche, ces sols perdent de l’eau et leurs volumes
diminuent, il se produit le phénomène de retrait. Le gonflement et le retrait sont
généralement accompagnés d’un changement des caractéristiques mécaniques du sol qui
peuvent influencer la durée de vie de l’ouvrage.
De nombreux travaux ont été réalisés afin d’étudier le phénomène de retrait-gonflement

qui ne garantit pas la pérennité des ouvrages construits sur les sols gonflants.
I.1.1 Définition
Il n’existe pas une définition unique du terme argile. Le mot argile englobe deux
connotations, l’une liée à la taille des grains et l’autre à la minéralogie. La définition dépend
de la discipline concernée. Il désigne selon les cas, un ensemble d’espèces minérales, une
famille de roches, une catégorie de sols ou encore une classe granulométrique. On utilise
le mot argile pour désigner les matériaux dont les composantes majeures sont des
minéraux argileux. L’argile se forme le plus souvent par altération chimique des roches qui
affleurent à la surface de l’écorce terrestre sous l’action des eaux [3]. Le domaine
privilégié de l’argile se trouve donc dans les zones superficielles, où elle est associée à
des dépôts sédimentaires, provenant des roches magmatiques, métamorphiques ou
sédimentaires rencontrées dans le substratum régional. Les argiles ont trois origines
possibles (Millot 1964) [12]:
 L’héritage, c’est l’apport de minéraux formés ailleurs sur le continent et

transportés soit par l’eau, soit par le vent. Au cours du transport, les argiles déjà
formées ne subissent que de faibles modifications de structure et de composition.
En général, elles arrivent ensuite dans la mer ou dans les bassins lacustres ;
 La transformation, les minéraux argileux conservent certains caractères
structuraux des matériaux originels, mais ils subissent au cours du transport des
transformations qui peuvent être importantes et qui les font passer d’un type
minéralogique à un autre par élimination de certains cations ou par fixation d’ions;

 La néogènese, les minéraux argileux sont dans ce cas formés par précipitation
d’ions en solution (Si4+, Al3+ et K+) susceptibles de réagir entre eux quand ils sont
mis en contact par diffusion ou mélange. Les composés ainsi formés sont stables
dans les conditions physicochimiques du milieu où ils ont pris naissance.
I.1.2 Aspect microstructure des argiles

I.1.2.1 Minéraux argileux
Les argiles sont constituées de fines particules de dimensions inférieures à 2 m chargées
négativement, leur présence rend l’argile fortement sensible aux molécules d’eau. Cette
sensibilité se traduit par le gonflement, le retrait et le changement des paramètres
mécaniques. C’est la variation de la teneur en eau qui est le facteur déclenchant des
variations de volume. Une augmentation de la teneur en eau entraîne le phénomène de
gonflement tandis que sa diminution induit le phénomène de retrait [13].
I.1.2.2 Les caractéristiques minéralogiques des argiles

Les minéraux argileux font partie du groupe des phyllosilicates qui contiennent une grande
proportion d'eau emprisonnée entre des feuillets de silice. Ils sont formés par la
dégradations lente de minéraux préexistants. La grande diversité des argiles découle des
conditions dans lesquelles ont eu lieu la dégradation des minéraux argileux et leur état
d’avancement. La cellule de base des minéraux argileux est appelée cristallite. Elle est
constituée d’un feuillet et d’un interfeuillet appelé aussi espace interfoliaire ou interstice.
Une particule d’argile résulte de l’empilement face-à-face de quelques cristallites
élémentaires. La phase solide des argiles est structurée suivant une organisation spatiale
particulière : le feuillet, la particule et l’agrégat.
 Le feuillet :
Les feuillets sont organisés en plusieurs couches superposées. On distingue deux
types de couches à partir desquelles sont construites les argiles :
 La couche Tétraédrique (T) : de formule générale SiO4, association de
quatre oxygènes formant les sommets d’un tétraèdre (Figure I-1). Dans la
couche tétraédrique, le cation dominant est Si4+ mais de fréquentes
substitutions ont lieu avec Al3+ et, plus rarement, avec Fe3+.
 La couche Octaédrique (O) : de formule générale Al2(OH)6 ou Mg3(OH)6, est
formée par deux plans d’anions oxygène O2- ou d’anions OH- aux sommets
de chaque octaèdre, enserrent aux centre un cation tel que : Al3+, Fe2+, Mg2+
ou Fe3+ (Figure I-2).

Figure I-1: Couches tétraédriques de slice source :[9]
Figure I-2: Couches octaédriques centré sur l’aluminium source : [9]
 La particule
La particule encore appelée cristallite ou tactoïde désigne une superposition de
feuillets argileux. La disposition des feuillets et le nombre de feuillets par
particule, sont variable suivant le type d’argile et son état hydrique. Lorsque l’on
observe une particule argileuse on peut différencier deux types de porosités :
 La porosité interfoliaire définie entre deux surfaces internes de feuillets
d’une même particule.
 La porosité interparticulaire définie entre deux surfaces externes de deux
particules.

Figure I-3: Représentation schématique de quelques groupes de minéraux

argileux TO (1:1) [1]
 L’agrégat
L’agrégat est un assemblage désordonné entre les particules argileuses, dont la
forme et les dimensions sont variables. A l’échelle de l’agrégat, le nombre de
niveaux structuraux est limité à deux:
 La microstructure qui correspond aux agrégats constitués par l’assemblage
des particules argileuses
 La macrostructure désigne l’assemblage macroscopique des agrégats.[9]
Figure I-4: Différents niveaux structuraux de sols argileux [9]
I.1.3 Les différentes familles de minéraux argileux

En prenant en compte la combinaison des couches d’octaèdres et de tétraèdres dans un
feuillet élémentaire, la charge de la couche et la nature du matériel dans l’espace
interfoliaire, on peut distinguer différentes familles de minéraux argileux disposant de
caractéristiques structurales, morphologiques et microscopiques spécifique, mais trois

grands groupes d’argiles sont souvent rencontrés en géotechnique : les kaolinites, les
illites et les smectites.
I.1.3.1 Les Kaolinites

Chaque feuillet de cette famille résulte de l’assemblage d’un feuillet tétraédrique SiO4
avec un feuillet octaédrique AlO2(OH)4 Figure I-5. Les argiles de cette famille sont
dites argiles T-O. Les couches de base successives sont maintenues par les liens
d’hydrogène entre les hydroxydes du feuillet en octaèdre et les atomes d’oxygène du
feuillet en tétraèdre. Ces liaisons considérées comme très fortes, empêchent
l’hydratation entre les feuillets élémentaires permettant un empilement relativement
important de couches et la formation d’un cristal d’assez grande dimension. L’espace inter-
foliaire de cette famille est de 7.2 Å. La kaolinite est une argile stable, non-gonflante et
sa structure élémentaire n’est pas affectée par la présence d’eau et ne peut être détruite
que par chauffage à 450°C. La formule générale d’un feuillet de kaolinite est donnée par
la formule suivante :
n[Al2Si2O5(OH)4] [13]
Figure I-5: Vue isométrique de la structure atomique de la

kaolinite [13]
Figure I-6: Structure élémentaire du réseau cristallin de la

kaolinite [28]

I.1.3.2 Les Montmorillonites

Ces minéraux sont également formés de feuillets à trois couches et sont des argiles T-
O-T. L’épaisseur d’un feuillet est de 10A°, il est constitué par l’association de deux
couches tétraédriques de type TO4(OH)2. T pouvant être AL3+, Fe3+, Mg2+, K+. (Figure I-7)
La formule générale de la montmorillonite est la suivante :
n[(Al1.67Mg0.33)Si4O10(OH)2] [13]
Figure I-7: Vue isométrique de la structure atomique de la montmorillonite [13]
Figure I-8: Structure élémentaire du réseau cristallin de la

montmorillonite [28]

Les liaisons d’un feuillet à l’autre sont faibles et l’eau pénètre facilement entre les
feuillets. Pour cette raison, les ions échangeables peuvent pénétrer dans la structure et
séparer les couches. La faiblesse de ces liaisons leur confère de grandes surfaces
spécifiques et de forts potentiels de gonflement. L’argile contenant un fort pourcentage
de montmorillonite est considérée comme une argile de nature gonflante.
I.1.3.3 L’illite
Ces argiles sont donc aussi appelées argiles T-O-T. La structure de l’illite est analogue à
celle de la montmorillonite, mais la présence d’ions potassium K+ dans la structure,
rééquilibre et renforce la liaison entre les feuilles, rendant ainsi l’argile moins sensible
aux phénomènes de retrait-gonflement. L’illite a un comportement intermédiaire (les
latérites font partie de cette famille d’argile).
Figure I-9: Vue isométrique de la structure atomique et

schématique de l’illite [13]
Figure I-10: Structure élémentaire du réseau cristallin de l’illite

Les différentes caractéristiques de chaque argile citée ci-dessus sont résumées dans le
tableau suivant :
Tableau I-1: Caractéristiques microscopiques des principales argiles [14]
Argile Kaolinite Illite montmorillonite

Type de feuillet 1 :1 2 :1 2 :1
Épaisseur de 7 10 9.6
feuillet (A°)
Nombre de feuillet 10-150 5-20 1-10 (montmorillonite
par particule sodique)
10-40 (montmorillonite
calcique)
CEC (méq/100 g) 5-15 10-40 80-150
Surface spécifique 10-45 80-100 700-900
(m2 /g)
Matériau - K+ H20, Na+, Ca2+ …..
interfoliaire
Comportement dans Non-gonflant Peu-gonflant Gonflant
l’eau
I.1.4 Les phénomènes de gonflement et de retrait des argiles

Le processus de gonflement des sols argileux au contact de l’eau met en jeu séparément
ou de façon combinée des phénomènes physico-chimiques et mécaniques variés.
I.1.4.1 Gonflement au sens physico-chimique

Les particules argileuses sont de type colloïdal, c’est-à-dire qu’elles ne sont pas contrôlées
par des forces gravitationnelles, mais par des forces de surface. Celles-ci provoquent
donc le phénomène de gonflement dont le mécanisme peut se produire à deux niveaux : au
niveau interfoliaire et interparticulaire[13].
 Le gonflement interfoliaire ou intraparticulaire: Il s'agit de l'infiltration de

molécules d'eau à l’intérieur des particules d’argiles, entre les feuillets
élémentaires. Lors d'une hydratation inter particulaire, les changements de volume
qui s'ensuivent sont généralement faibles. Les caractéristiques de ce gonflement
dépendent essentiellement de la structure du minéral gonflant. Le gonflement
interfoliaire est dû à l’hydratation et à la répulsion entre les particules, soit
l’augmentation de la distance interfeuillet [15].
Lors de l’hydratation d’une éprouvette d’argile sèche, l’eau pénètre à l’intérieur
des particules (dans l’espace interfoliaire) et s’organise en couches moléculaire,

puis interagit avec les feuillets et les cations déjà présents. Cette hydratation
favorise la création d’une pression, appelée pression de disjonction, qui a pour
conséquences une augmentation de la distance interfoliaire et donc un
gonflement de l’ensemble des deux feuillets : c’est le gonflement cristallin. Le
phénomène de gonflement dépend de la nature des cations interfoliaires. Par
exemple la montmorillonite sodique (Na+) gonfle plus que la montmorillonite calcique
(Ca2+) à cause de la facilité qu’à l’eau de se répartir sur les feuillets. Ce type de
gonflement ne représente environs qu’une dizaine de pourcentage du gonflement
macroscopique total. Il est donc nécessaire d’envisager un autre mécanisme de
gonflement que celui relevant des seules variations interfoliaires.
 Gonflement interparticulaire : Le gonflement interparticulaire est dû à l’attraction

entre les particules d’argile chargées négativement en surface et la charge
positive du dipôle des molécules d’eau. Sur la surface des particules d’argile, se
fixe ainsi une première couche d’eau dont les extrémités négatives attirent à leur
tour une seconde couche et ainsi de suite. Mais l’orientation préférentielle des
premières couches d’eau autour de la particule d’argile tend à diminuer pour les
couches d’eau plus lointaines, la transition entre l’eau liée et l’eau libre
interstitielle étant graduelle. En effet l’absorption des cations par les particules
argileuses et la formation de la double couche diffuse sont responsables de la
répulsion de deux particules. En effet l’absorption des cations par les particules
argileuses et la formation de la double couche diffuse sont responsables de la
répulsion de deux particules. Dans une argile sèche, les cations sont très
fortement absorbés à la surface des feuillets argileux. Les cations en excès de
ceux nécessaires à assurer l’électroneutralité des particules argileuses et les
anions associés sont présents sous forme de sels. Si l’argile sèche est placée au
contact de l’eau, les sels se dissolvent dans la solution, mais les différents ions se
répartissent de manière non homogène : la concentration en cation devient très
forte au voisinage des surfaces argileuses (surfaces externes), tandis que les
anions ont plutôt tendance à s’éloigner des surfaces argileuses. La zone perturbée,
ainsi créée par cette répartition non-homogène de cation et d’anion est la double
couche diffuse.
En résumé le gonflement au sens physico-chimique relève de deux cause :
 L’hydratation des espaces interfoliaires entrainant un accroissement de leurs

dimensions (augmentation de la porosité interfoliaire) : gonflement cristallin ;
 La répulsion entre les particules (augmentation de la porosité interfoliaire et
interparticulaire) : gonflement osmotique.

I.1.4.2 Gonflement au sens mécanique

Le phénomène de gonflement de même que de tassement peut provenir d’une modification
de l’état de contraintes dans le sol en présence d’eau. En effet, considérons un élément
de sol saturé à l’équilibre.
Si une contrainte extérieure σe est appliquée à cet élément de sol, l’eau étant moins
compressible que le squelette solide, la contrainte σe est immédiatement reprise par la
phase liquide et au temps initial t=0, on a :
𝑢𝑤 = σ𝑒 , σ’ = 0
Avec uw désignant la pression interstitielle à savoir la pression de l’eau des pores et σ’

contrainte effective, c'est-à-dire la contrainte qui s’exerce réellement sur le squelette
solide [15].
Si l’on permet alors au sol de se drainer, un phénomène de consolidation va se développer,

correspondant à l’expulsion de l’eau et au transfert de contrainte de l’eau sur le squelette
solide. Un nouvel équilibre va alors s’établir pour l’état de contraintes :
𝑢𝑤 = 0, σ’ = σ𝑒
Si la contrainte σe est alors supprimée, le même phénomène se produit en sens inverse et

immédiatement :
𝑢𝑤 = −σ𝑒 , σ’ = σ𝑒
La pression interstitielle devient négative et opposée à la variation de contrainte totale.

Dans des conditions de libre circulation de l’eau, un phénomène de gonflement va alors se
développer exprimant l’adsorption de l’eau et le transfert de contrainte négative de l’eau
sur le squelette solide jusqu’à l’état final :
𝑢𝑤 = 0, σ’ = 0
Il est donc possible d’affirmer que si la consolidation exprime une diminution de la

pression interstitielle jusqu’à son annulation, le gonflement exprime quant à lui la
diminution de la succion (𝑢𝑤 ) jusqu’à son annulation [15].

I.1.5 Le retrait de l’argile

Le phénomène de retrait des sols est un phénomène bien décrit. Les sols fins,
communément appelé argiles, marnes, limons en géotechnique routière contiennent une
proportion variable de minéraux argileux (au vrai sens du terme, correspondant à des
phyllosilicates). Ces minéraux, de par leur structure interfoliaire ont la faculté de pouvoir
procéder à des échanges d’ions ou de molécules. Selon le type de structure et donc le type
d’argile, la capacité d’échange s’accompagne d’une variation volumique : très faible voire
nulle dans le cas des kaolinites et des illites, ces variations deviennent très importantes
pour le groupe des smectites (parmi lesquelles on trouve les montmorillonites, les
saponites) qui contiennent des molécules d’eau entre les feuillets. La perte d’eau se
traduit par une diminution volumique du sol et par un retrait se manifestant par des
fissures dites de retrait[16].
I.1.6 Localisation des sols gonflants au Bénin

Les sols gonflants se rencontrent dans les régions semi-arides situées dans les zones
tropicales ainsi que dans les régions tempérées où les précipitations sont plus faibles
pendant toute ou une grande partie de l'année par rapport à l'évaporation.
Au Bénin, ces sols se rencontrent dans la dépression de la Lama où affleurent des argiles
marneuses ou calcaires traversant d'est en ouest le bassin sédimentaire (sud du pays). La
dépression de la Lama (Lama signifie boue en portugais) encore appelée dépression
médiane est formée par une bande orientée W.SW-ENE dont la largeur est de 25 km. Elle
couvre une aire qui s’étend de l’Est (Pobè, Adja-Ouèrè, Issaba, Onigbolo) à l’Ouest (Allada,
Toffo, Ouèdèmè) du Bénin et est estimée à plus de 3000 km2.
On retrouve également les sols gonflants sur certaines alluvions argileuses des fleuves
Mono, Ouémé et Niger et dans le centre du pays où ils se localisent sur des secteurs bien
circonscrits : le secteur des gabbros de Mbétékourou, à l'est de Dassa-Zoumè, sur des
diorites, à l'est de Setto et au sud de Logozohè sur des embréchites et basiques divers.
Aussi, on dénombre des poches d’argiles et de tourbe autour du lac Nokoué, de la rivière
Djonou à Godomey, et de la lagune de Porto-Novo.
Dans ces régions, les dégâts sont multiples pendant et après la pluie : soulèvement des
voies, boursouflement des chaussées, rupture des canalisations, fissuration des murs,
basculement des édifices, glissement des terrains et des fondations. Les bâtiments
conçus sans précautions sur des sols gonflants sont quasi-détruites après quelques
saisons, par l’action conjuguée d’une part du retrait-fissuration qui favorise la pénétration
de l’eau et par conséquent le gonflement en profondeur, puis d’autre part, des lézards qui
suivent la saison des pluies lorsque l’argile gorgée d’eau ne présente plus aucune résistance
à la charge verticale[1].

Figure I-11: Carte de la dépression médiane [29]

I.2 Généralités sur les chaussées revêtues

Les chaussées se présentent comme des structures multicouches mises en œuvre sur un
ensemble appelé plate-forme support de chaussée constituée du sol terrassé (dit sol
support) surmonté généralement d'une couche de forme. D’une manière générale, la
structure type d’une chaussée revêtue se présente comme ci-après :
Figure I-12:Structure type d’une chaussée (SETRA & LCPC, 1994)[1]
Globalement on distingue deux grandes familles de types de chaussées : les chaussées

flexibles et les chaussées rigides. Selon le guide français de dimensionnement des
chaussées (SETRA & LCPC, 1994), ces familles se subdivisent à leur tour en six types de
chaussées comme suit :
Tableau I-2: Différents types de chaussées selon la classification LCPC-SETRA

(SETRA & LCPC, 1994)
Nature des couches constitutives

Famille de Types de
Couche de Couche de
chaussées structures Couche de Base
Surface Fondation
Souples
Bitumineuse Matériaux granulaires
Chaussées
Flexibles Bitumineuses
Matériaux bitumineux
épaisses

Semi-rigides Matériaux Traités aux Liants

Bitumineux
Hydrauliques (MTLH)
Mixtes Bitumineux MTLH
Graves
Inverses Bitumineux Non MTLH
Traitées
Béton de ciment (avec

Chaussées
Rigides ou sans dispositif de MTLH ou Béton maigre
rigides
liaison, les goujons)
I.3 Comportement et pathologie des structures de chaussées sur sols

gonflants
I.3.1 Une route sur sol gonflant
La sécheresse a un impact sur les structures de chaussées par le biais du retrait des sols
supports. Ce phénomène se manifeste assez classiquement par une fissuration
longitudinale, très peu ramifiée et proche de la bande de roulement en rive de chaussée.
Les figures suivantes représentent les stades successifs d’évolution d’une route
construite sur des sols argileux soumis à des phénomènes de retrait et gonflement. La
sécheresse produit un retrait du sol à partir des accotements, vers le centre de la
chaussée. La surface du terrain se fissure et le retrait du sol sous la chaussée peut aussi
créer des fissures dans le revêtement (Figure I-14). À la saison des pluies suivante, le
sol se réhumidifie et regonfle, mais il devient plus mou et le trafic peut déformer et
fissurer encore plus la chaussée (Figure I-15). La sécheresse suivante réactive le
tassement et les fissures, ce qui dégrade progressivement la route (Figure I-16). Ce
processus se produit dans les régions argileuses de beaucoup de pays, notamment en
Afrique.

Figure I-13: Profil théorique initial [16]
Figure I-14:Période de sécheresse : tassement et fissuration des zones

asséchées. Le sol durcit.[16]
Figure I-15: Période de pluies : gonflement, disparition des fissures, le sol

ramollit, la circulation crée des ornières et des fissures dans la chaussée[16]
Figure I-16: Période de sécheresse : tassement et fissuration des zones

asséchées. L’état de la chaussée empire.[16]

I.3.2 Pathologie des structures de chaussées sur sols gonflants

Les figures suivantes montrent des exemples de routes fissurées par l’effet du retrait-
gonflement du sol de fondation argileux. La position des fissures dépend de l’étendue de
la zone atteinte par les variations de teneur en eau lors des cycles de sécheresse et
d’humidification.
Photo I-2: Exemple N°1 d'une

Photo I-1: Exemple N°2 de fissure de
fissure de retrait [16]
retrait[16]
Photo I-3:Affaissement longitudinal de Photo I-4:Fissuration au KP 39+000 de

la route Kétou-Pobè dans la zone de la la route Kétou-Pobè, dépression de la
Lama, Bénin[9] Lama, Bénin

Chapitre II : Matériels, matériaux et

méthodologie

Chapitre II : Matériels, matériaux et méthodologie

Notre étude s’intéresse à la commune de Pobè située dans la dépression de Lama dans la
république du Bénin. En se penchant sur l’étude générale des sols du Bénin, les recherches
et réflexions issues de la quatorzième réunion du Sous-comité ouest et centre africain
de corrélation des sols qui s’est tenu à Abomey du 09 au 13 Octobre 2000 stipulent que
ces sols se sont formés à partir des roches sédimentaires, du socle cristallin et des dépôts
alluviaux et ont été corrélés avec la BRM [17]. Ce chapitre présente le milieu d’étude à
travers sa situation géographique, sa population, le type du relief et le type du sol
rencontré dans cette zone. Les dispositifs utilisés lors de l’étude expérimentale ainsi que
la méthodologie mise en place pour la caractérisation des sols étudiés y sont également
présentés.
II.1 Milieu d’étude

II.1.1 Situation géographique
La commune de Pobè située au Sud-est du Bénin et au centre – est du département du
Plateau à la limite frontalière avec le Nigeria, est limitée au Nord par la Commune de
Kétou, au Sud et à l’Ouest par la Commune d’Adja-ouèrè et à l’Est par le Nigeria. Elle a
une superficie de 400 km2, représentant 11% de la superficie du département et 0,46%
de celle du Bénin[9].
II.1.2 Population et structure administrative

La population de la commune de Pobè est de 54 181 habitants en 1992 et 82 910 habitants
en 2002 soit un taux d’accroissement naturel de 4,35% sur cette période. Lors du
recensement de 2013 (RGPH-4), la commune comptait 123 677 habitants[9]. La densité
de population est de 309 habitants/km2. La Commune de Pobè est composée de quarante-
quatre (44) villages/quartiers de ville qui sont répartis dans cinq (05) arrondissements
dont un urbain (Pobè centre) et quatre ruraux (Towe, Igana, Ahoyèyè, Issaba). La
prolifération de la démographie observée au cours de ces dernières années dans cette
commune requiert donc la mise en place d’infrastructures pérennes dans l’optique de
favoriser l’épanouissement de cette population.
II.1.3 Climat
La Commune de Pobè jouit d’un climat de type subéquatorial avec deux (02) saisons
pluvieuses et deux (02) saisons sèches qui s’alternent. Les précipitations atteignent 1100
à 1200 mm par an. A l’instar des variations climatiques dans la sous-région, la commune de
Pobè ne bénéficie plus de ce type de climat. La zone de dépression est une zone de forte
rosée [18].
Les statistiques météo de ces 20 dernières années dans la commune de Pobè nous donnent
une humidité relative moyenne de 78,1%, une température moyenne sous abri de 27,6°C.

La durée annuelle moyenne de la journée H est de 12,5 ; la vitesse moyenne annuelle du

vent KMH (basculer vers MPH) est de 9,5 km/h et les précipitations moyennes annuelles
sont de 112 cm [9].
II.1.4 Relief [9]

Le milieu physique est situé globalement dans une zone de plateaux de 50 à 200 mètres
au-dessus du niveau de la mer et d’une altitude moyenne de 100 mètres. La Commune de
Pobè comporte une dépression qui n’est rien d’autre que la continuité de la dépression
médiane qui traverse tout le Bénin d’Ouest à Est. Elle découpe la Commune en deux zones
orographiques : une zone de dépression (altitude <50m) et une zone de plateau. La zone
de dépression (altitude <50m) comporte une grande partie de l’Arrondissement d’Issaba
jusqu’à la latitude d’Onigbolo et les arrondissements d’Ahoyèyè, d’Igana et de Towé. La
zone de plateaux regroupe tout le reste de la commune. Sur celle-ci on rencontre plusieurs
plateaux dont le plateau de Pobè d’altitude moyenne 100 mètres.
II.1.5 Sols [9]

Du point de vue physique, dans la commune de Pobè, on distingue deux types d’unité
pédologique : un sol hydromorphe et un sol ferralitique. Le sol hydromorphe très fertile
où se pratiquent les principales cultures agricoles est situé dans la dépression d’Issaba.
C’est la zone de prédilection des habitants Holli (ou hollidjé : un groupe socio culturel
proche des Yorubas qui en grande partie vit de l’agriculture). L’habitat rencontré dans
cette zone est de type dispersé à travers des cases disséminées dans les brousses. Il
occupe les (3/4) de la superficie de la commune. Cette zone est composée de roche
calcaire, ce qui a permis l’implantation d’une usine cimentière gérée par la SCB Lafarge à
Onigbolo. Le sol ferralitique est situé sur le plateau Pobè-Sakété. Il est composé de sable,
de grès et d’argile. C’est un sol rouge qui occupe le quart (1/4) de la superficie de la
commune.

Figure II-1: Carte Pédologique de Pobè
II.1.6 Présentation des sites de prélèvement

Les coordonnées géographiques en projection sont prises par rapport à la référence
mondiale UTM WGS 84 et en OXY. Ainsi le tableau suivant nous renseigne sur les
coordonnées géographiques de chaque site de prélèvements. Ces données sont consignées
dans le Tableau II-1.
Tableau II-1:Localisation des sites de prélèvements
Coordonnées
Localités
X Y Latitude Longitude
ISSABA 459660 783707 N 07°05'23,8" E 02°38'04,9"
AHOYEYE 465361 774802 N 07°00'33,9" E 02°38'04,9"
II.2 Matériels d’étude

Le matériel utilisé comprend l’appareils de cisaillement, la cellule œdométrique, l’étuve,
les bagues et le bocal en verre pour l'humidification.

Photo II-2: Bocal à humidification Photo II-1: Etuve
Photo II-4: Appareil de cisaillement

Photo II-3: Bâtis œdométrique
Photo II-5: Bague pour l'essai Photo II-6: bague rectangulaire pour
œdométrique l'essai de cisaillement

II.3 Méthodologie de recherche

Pour la détermination des caractéristiques géotechniques des sols de notre milieu d’étude,
plusieurs essais ont été effectués au Centre National d’Essai et de Recherches des
Travaux Publics (CNERTP). Il s’agit des essais d’identification physique et mécanique des
sols. Nous présenterons dans ce chapitre les différents matériels et méthodes de ces
essais qui se basent sur les normes françaises et les méthodes expérimentales appliquées
dans cette structure d’accueil.
II.3.1 Recherche et localisation de site de prélèvement

La recherche des sites a été effectuée en nous servant au préalable de la carte
pédologique de la commune ensuite nous avons effectué une descente sur les lieux afin de
définir des sites de prélèvement des échantillons et en choisissant les zones où le
phénomène se constate le plus dans la commune. Les arrondissements Ahoyèyè et Issaba
sont les zones les plus critiques et présentent de ce fait plus de pathologies. Compte tenu
des essais à réaliser (essais d’identification et mécaniques) nous avons procédé aux
prélèvements d’échantillons remaniés et non-remaniés.
II.3.2 Prélèvements
Nous avons procédé aux prélèvements dans les différentes régions de la commune de Pobè
le jeudi 15 Septembre au Vendredi 16 Septembre 2022. Compte tenu des essais à réaliser,
nous avons effectué des prélèvements d’échantillons intacts et remaniés. Pour ce qui
concerne les prélèvements intacts, nous avons utilisé des caisses cubiques métalliques de
section 20*20cm2, que nous avons battues avec une dame. Dans chaque zone d’étude, nous
avons fait 4 prélèvements intacts à raison d’un prélèvement par couche de 0.5m de
hauteur, donc de 0 à 0.50m, de 0.5 à 1.0m, de 1.0 à 1.50m et de 1.50 à 2.0m. A chaque
couche, nous enfonçons une caisse cubique métallique de 0,2 m de hauteur par battage
dans le sol en procédant de manière que l’enfoncement soit le plus uniforme possible et
sans endommager la caisse. La terre située autour de la caisse est soigneusement enlevée.
Pour disposer des échantillons dont les caractéristiques hydriques sont les plus proches
possibles des conditions in situ au moment du prélèvement, de la paraffine est appliquée
sur les deux faces libres du bloc. Arrivé au laboratoire ces échantillons sont découpés et
prélevés afin de réaliser les essais. Les échantillons remaniés quant à eux sont prélevés
dans des sacs de jute étiquetés selon la profondeur.
II.3.3 Description des essais

II.3.3.1 Essais d’identification
II.3.3.1.1 Essai de détermination de la teneur en eau pondérale (NFP 94 -050) [19]
Cet essai a pour but de déterminer la quantité d’eau (mw) contenue dans une masse humide
(mh) d’échantillon prélevé et étuvé. Cette quantité d’eau est exprimée en pourcentage
d’eau par rapport à sa masse sèche (ms) ; c’est la teneur en eau.

𝑚ℎ − 𝑚𝑠
𝑤= ∗ 100
𝑚𝑠
II.3.3.1.2 Limites d’Atterberg (NF P 94-051)[20]

Les limites d’Atterberg sont des paramètres d’identification des sols fins nous
permettant de connaître la limite de liquidité (ωL), de plasticité (ωP) et l’indice de
plasticité (IP). Cela permet d’obtenir les teneurs en eau qui délimitent les différents états
de consistance d’un sol. Ces limites qui s’expriment en pourcentages (%), sont les
suivantes:
 La limite de retrait (ωs)

C’est la teneur en eau maximale que le sol peut avoir sans changer de volume. Cette
limite sépare l’état solide sans retrait de l’état solide avec retrait.
 La limite de plasticité (ωP)
Elle est définie comme la teneur en eau d’un sol qui a perdu sa plasticité et se
fissure en se déformant lorsqu’il est soumis à de faibles charges. Cette limite
sépare l’état plastique de l’état semi-solide. Elle varie de 0 % à 100 %, mais elle
demeure généralement inférieure à 40 %.
 La limite de liquidité (ωL)
C’est la teneur en eau qui sépare l’état liquide de l’état plastique. Dans la plupart
des cas elle ne dépasse pas 100 %. Les limites de plasticité et de liquidité servent
à identifier et à classifier les sols à grains fins. Quant à la limite de retrait, elle
sert à l’étude de certains sols dont le volume varie fortement à cause des
changements de la teneur en eau. C’est surtout dans les régions arides, où se
trouvent des dépôts de la montmorillonite (problème de retrait-gonflement).
Ces différentes limites d’Atterberg permettent de déterminer les indices de plasticité

et de liquidité.
 L’indice de plasticité (IP)

Il s’exprime en pourcentage, correspond à la différence entre la limite de liquidité
et la limite de plasticité :
𝐼𝑝 = ω𝐿 − ω𝑝
Il permet de déterminer le degré de plasticité d’un sol.

Figure II-2: Diagramme de plasticité
Le Tableau II-2 nous donne une interprétation du potentiel du gonflement d’après

Snethen en fonction de la valeur de l’indice de plasticité IP.
Tableau II-2:Potentiel de gonflement [21]
Indice de plasticité IP (%) Potentiel de Gonflement

>35 Très élevé
22 – 48 Elevé
22 – 32 Moyen
<18 Faible
 L’indice de consistance (Ic)

Cet indice permet de savoir l’état hydrique d’un sol. Il correspond au rapport de la
différence entre la limite de liquidité et la teneur en eau naturel in-situ du sol (w)
à l’indice de plasticité.
ω𝐿 − ω
𝐼𝑐 =
𝐼𝑝

Figure II-3: Diagramme de consistance
 L’indice de liquidité (IL)

Cet indice permet de savoir si un sol est à l’état liquide, plastique, semi-solide ou
solide. Pour établir cet indice, on calcule le rapport de la différence de la teneur
en eau naturelle in situ (w) d’un sol et ses limites de plasticité à l’indice plasticité.
ω − ω𝑝
𝐼𝑙 =
𝐼𝑝
II.3.3.1.3 Détermination du Poids spécifique des grains solides (NF P 94-054)[22]

L’essai consiste en la détermination de la masse volumique des grains solides d’un
échantillon de sol à l’aide de pycnomètre. Une petite quantité d’échantillon sec est
introduite dans le pycnomètre. On y ajoute de l’eau jusqu’au trait de jauge. Après 72
heures de repos, temps requis pour que l’eau soit complètement absorbée, les masses
suivantes sont prises au cours de l’essai :
M1, Masse du pycnomètre vide ;
M2, Masse du pycnomètre plus échantillon sec ;
M3, Masse du pycnomètre rempli d’eau plus échantillon ;
M4, Masse du pycnomètre rempli d’eau.
La masse volumique des grains solides est ensuite calculée à l’aide de la formule suivante:
𝑀2 − 𝑀1
𝜌𝑠 = 𝜌𝑤 ∗
𝑀4 + 𝑀2 − 𝑀1 − 𝑀3
Par convention, la masse volumique de l’eau 𝜌𝑤 est égale à 1000 kg/m3. L’essai est réalisé
à deux reprises. Le résultat constitue la moyenne des deux valeurs.
II.3.3.1.4 Essai de détermination de la teneur en matière organique e (XP P 94-055)[23]

a. Définition et but de l’essai

Cet essai a pour objectif la détermination de la proportion de matière organique contenu

dans 80g de matériau sec pour le tamisât du tamis 400 µm d’un matériau humide au moyen
d’eau oxygéné. La pollution d’un sol par les débris végétaux entraine souvent des
conséquences dommageables si celui-ci doit être utilisé pour la construction
d’infrastructures importantes. Il existe plusieurs méthodes de détermination de la teneur
pondérale en matières organiques parmi lesquelles :
 La méthode chimique suivant la norme NF P 94-055 (Décembre 1993) ;

 La méthode par calcination selon la norme XP P 94-047(Décembre 1998) ;
 La méthode utilisant l’eau oxygénée. Elle est non normalisée mais se retrouve dans
la littérature.
b. Principe de l’essai
La méthode d’eau oxygénée est couramment utilisée pour la détermination du dosage en

matières organiques. L’eau oxygénée a une action limitée sur les restes de végétaux non
décomposés tels que les racines et les fibres. Elle est souvent utilisée dans le
prétraitement des sols avant l’analyse granulométrique par sédimentation en vue
d’éliminer les matières organiques colloïdales. On détermine ainsi le pourcentage de
matière organique :
𝑃1 − 𝑃2
𝑀𝑂 = ∗ 100
𝑃1
Avec :
- 𝑃1 , le poids initial
- 𝑃2 , le poids après essai de l’échantillon
- 𝑃1 − 𝑃2 , la perte de poids
c. Interprétation des résultats
Les seuils retenus sont présentés par le Tableau II-3 :

Tableau II-3:Désignation géotechnique des sols en fonction de la teneur en

matières organiques[9]
Teneur en matière organique Désignation géotechnique

MO ≤ 3 Non organique
Faiblement
3<MO ≤10 Vase
organique
Moyennement
0 <MO ≤ 30 Sol tourbeux
organique
MO > 30 Très organique Tourbe
II.3.3.1.5 Analyse granulométrique par tamisage (NF P 94-056) [24]

L’analyse granulométrique par tamisage est un essai d’identification qui permet de
connaitre la granulométrie d’un matériau tout en le faisant passer par une série de tamis
normalisés dont les ouvertures vont en décroissant du haut vers le bas. Les particules les
plus grosses restent donc emprisonnées sur les tamis les plus hauts (refus ou retenu),
tandis que les particules plus fines se dirigent vers les tamis inferieurs (tamisât ou
passant). Les résultats sont portés sur un graphique semi-logarithmique qui permet
d’obtenir une courbe granulométrique.
II.3.3.1.6 Analyse granulométrique par sédimentométrie (NF P 94-057)[24]

L’essai d’analyse granulométriques par sédimentation est réalisé afin de compléter la
courbe granulométrique par tamisage sur les particules de tailles inférieures à 80 µm.
L’essai est réalisé sur les passants issus du tamisage par voie humide lors de l’exécution
de l’analyse granulométrique par tamisage.
II.3.3.2 Classifications utilisées pour la caractérisation des matériaux

Les classifications utilisées afin de préciser la nature de ces sols sont la classification
américaine des sols intitulée Highway Research Board (HRB) et la classification GTR (NFP
11-300).
a. La classification américaine "Highway Research Board" (HRB)
La classification américaine des sols Highway Research Board (HRB) prend en compte le
pourcentage des passants aux tamis 0,2 mm, 0,4 mm et 80μm ; l’indice de plasticité Ip
ainsi que la limite de liquidité WL.

Tableau II-4: Classification H.R.B
b. La classification GTR
La classification GTR est basée sur les caractéristiques granulométriques : Dmax, passant
à 80µm et sur la valeur au bleu du sol et/ou sur l’indice de plasticité IP. La classification
établie pour les sols de Dmax ≤ 50 mm à partir de ces différents paramètres est résumée
par la Figure II-4.
Figure II-4: Classification GTR

II.3.3.3 Essais mécaniques

II.3.3.3.1 Essai de cisaillement direct (NF P 94-071-1) [25]
Il s’agit de déterminer expérimentalement la courbe intrinsèque d’un sol puis d’en déduire
les paramètres de cisaillement (angle de flottement interne et la cohésion) qui permettent
d’estimer par exemple la contrainte de rupture sous une fondation superficielle. Dans
l’essai de cisaillement à la boîte de Casagrande, on cherche à provoquer la rupture de
l’échantillon suivant un plan imposé. L’échantillon est placé dans une boîte constituée de
deux parties qui peuvent glisser horizontalement l’une sur l’autre. On applique à
l’échantillon un effort normal de compression N, verticalement, par l’intermédiaire d’un
piston et un effort tranchant T, horizontalement, en déplaçant la demi-boîte inférieure.
Un comparateur mesure la variation de hauteur de l’échantillon.
II.3.3.3.2 Essai œdométrique (XP P 94-090-1)[26]

a. Essai œdométrique (XP P 94-090-1)
L’œdomètre est un appareil de laboratoire utilisé en géotechnique. Cet appareil permet

de mesurer le tassement d’un échantillon de sol (intact de préférence) sous un effort
donné, et son gonflement au cours du déchargement des efforts soumis auparavant.
L'opération de chargement et déchargement s'effectue par paliers et elle est
chronométrée, afin d’en déterminer les propriétés mécaniques lors de calculs de
mécanique des sols.
La manipulation de cet essai a pour but de déterminer les caractéristiques de

compressibilité d’un sol qui permettent d’estimer le tassement d’un massif de sol, par
exemple sous une fondation superficielle. Sous l’effet des charges appliquées, le sol va se
déformer et va subir un tassement ; pour en évaluer l’ampleur, on reproduit un phénomène
semblable au laboratoire. On placera l’échantillon dans une enceinte cylindrique empêchant
toutes déformations latérales lors d’un chargement par pallier. On suppose donc que dans
le sol, les déformations horizontales sont nulles.
b. Principe de l’essai
Le sol est placé dans une enveloppe rigide Figure 3-7, on exerce sur sa partie supérieure
une pression variable à l’aide d’un piston et on mesure les affaissements observés après
stabilisation. On détermine ainsi la relation entre les contraintes effectives et les
déformations verticales.
2𝑣 2
𝐸 = 𝐸𝑜𝑒𝑑 ∗ [1 − ( )] 𝑎𝑣𝑒𝑐 𝑣 = 0.33 𝐶𝑜𝑒𝑓𝑓𝑖𝑐𝑖𝑒𝑛𝑡 𝑑𝑒 𝑝𝑜𝑖𝑠𝑠𝑜𝑛 𝑝𝑜𝑢𝑟 𝑙𝑒𝑠 𝑠𝑜𝑙𝑠
1−𝑣

Par analogie avec le module de Young E (théorie de l’élasticité linéaire), on va définir le

module œdométrique : 𝐸𝑜𝑒𝑑 exprimé en MPa ou kPa.
NB : Cet essai consiste en la détermination pour un sol, entre autres indice de compression
Cc, indice de gonflement Cs, la contrainte verticale effective Ơ’vo, la contrainte de
préconsolidation Ơ’p, l’indice des vides du sol en place e0, le coefficient de consolidation
vertical Cv.
c. Interprétation des résultats
Tableau II-5: Nature du sol compressible en fonction de l’indice de compression Cc
Indice de compression Cc Nature

Cc < 0,02 Sol incompressible
0,02 < Cc < 0,05 Sol très peu compressible
0,05 < Cc < 0,1 Sol peu compressible
0,1 < Cc < 0,2 Sol moyennement compressible
0,2 < Cc < 0,3 Sol assez fortement compressible
0,3 < Cc < 0,5 Sol très compressible
Cc > 0,5 Sol extrêmement compressible
Tableau II-6: Nature des sols en fonction de l’indice de compression Cc
Indice de compression Cc Nature

0,01 < Cc < 0,1 Sable
0,1 < Cc < 0,25 Kaolinite
0,25 < Cc < 0,8 Illite
0,8 < Cc < 2,5 Montmorillonite
Tableau II-7:Nature des sols gonflants en fonction de l’indice de gonflement Cs
Indice de gonflement Cs Gonflement potentiel

0,04 < Cs < 0,25 Gonflement élevé
Dès que Cs dépasse 0,07 On peut considérer que le soulèvement des
fondations risque de poser de graves
problèmes

II.4 Description du mode opératoire de l’humidification et de séchage

II.4.1 Mode opératoire de l’humidification
a. On fait le carottage des échantillons intacts dans les bagues de cisaillements et
les bagues œdométriques. Ces bagues sont pesées à vide au préalable et après le
carottage. Pour chaque type d’échantillon on fait une carotte qui va servir de
contrôle pour le calcul de la teneur en eau d’humidification. Dans notre cas pour
les quatre (04) échantillons en considérant qu’on va s’arrêter à 3 cycle le nombre
total de carottes et la répartition se présentent dans le Tableau II-8 comme suit:
Tableau II-8: Tableau récapitulatif du nombre de carotte pour l'humidification
échantillons 0-0.5m 0.5-1m 1-1.5m 1.5-2m

Nombre de
carottes pour
12 12 12 12
essai de
cisaillement
Nombre de
carottes pour
3 3 3 3
essai
œdométrique
Nombre de
carottes pour
la vérification 6 6 6 6
de la teneur en
eau
Nombre total
21 21 21 21
de carottes
Photo II-8:Bague de carottage pour Photo II-7:Bague de carottage pour

l'essai de cisaillement l'essai œdométrique

b. On introduit les carottes dans les bocaux en mettant sous chaque carotte une
étiquette sur laquelle sont inscrits les informations du matériau, le poids de la
bague et le poids de l’ensemble bague et matériau. Les bagues pour le contrôle de
la teneur en eau et la saturation sont mises dans un autre bocal et celles dédiées
aux essais dans un autre bocal. Le bocal est muni d’un outil servant de support pour
les carottes. Ce support est perforé de plusieurs trous pour permettre à ce que
l’eau à l’intérieur du bocal circule librement et humidifie les carottes par le bas
également.
Photo II-10 :Support des bagues Photo II-9: Bac pour l'humidification
c. On verse l’eau dans le bocal jusqu’à ce que l’eau submerge les carottes et on ferme
ce dernier avec son couvercle.
d. Les carottes qui servent de contrôle sont pesées chaque une heure 1h pendant 48h.
Afin d’éviter au maximum la variation de la masse de matériau due au mouvement
de l’eau sur les carottes on sort l’ensemble support et matériau avec soin et
douceur. Pour le suivi on trace au fur et à mesure qu’on pèse les carottes, la courbe
de l’évolution de la masse en fonction du temps pour chaque carotte. Les carottes
dédiées à l’essai sont intouchables. Dès qu’on constate que trois masses
consécutives ne varient ou les variations ne dépassent (0.2), on considère que le
degré de saturation est à 100%. Pour chaque échantillon, une carotte servant de
contrôle de ce dernier ira à l’étuve pour la détermination de la teneur en eau du
mouillage et celles dédiées à l’essai vont passer au soleil pour subir la phase de
séchage avec une carotte servant de contrôle pour la teneur en eau de séchage.
Par exemple pour l’échantillon d’Issaba (0-0.5m) pour le premier cycle on aura au
total sept (07) carottes dont quatre (04) carottes de cisaillement, une (01) pour
l’essai œdométrique et deux (02) qui serviront de contrôle pour la vérification de
la teneur en eau. Après l’humidification les quatre (04) carottes de cisaillement et
celle de l’essai œdométrique avec une (01) carotte de contrôle seront soumises à
la phase de séchage. La seule qui reste passe à l’étuve pour le calcul de la teneur
en eau d’humidification.

II.4.1.1 Cinétique de l’humidification

Le Graphe II-1 nous présente la diminution de la variation de la masse en fonction du
temps.
Variation de la masse au cours de

l'humidification
370
350
330
Masse
310
290
270
250
0 5 10 15 20 25 30
Temps (Heure)
Graphe II-1: Cinétique de l'humidification

II.4.2 Mode opératoire du séchage

Après l’humidification, on lance la phase du séchage qui nous permet de simuler le
phénomène qui se produit en saison sèche au laboratoire. Le séchage est fait au soleil à
une température variant entre ±26°C et ±36°C. La température minimale le matin tourne
autour de 26°C et celle maximale s’obtient dans l’après-midi et tourne autour de 36°C
d’après les valeurs de température obtenues après des mesures de température faites
pendant 3 jours. Après l’humidification, les carottes sont exposées au soleil pour subir la
phase du séchage pendant une durée de 9h de temps par jours. La teneur en eau de
séchage est fixée à 14%. Cette valeur de teneur en eau est choisie après une revue de
littérature sur les teneurs en eau des matériaux expansifs dans certaines zones en saison
sèche.
II.4.2.1 Cinétique du séchage

Le Graphe II-1 nous présente la cinétique du séchage
Cinétique du séchage
70
60
50
Teneur en eau %
40
30
20
10
0
0 10 20 30 40 50 60
Temps (h)
Graphe II-2: Cinétique du séchage

Chapitre III : Caractérisation au laboratoire des sols

naturels

Chapitre III : Caractérisation au laboratoire des sols naturels

Introduction
Le but de ce chapitre est d’élaborer une synthèse des caractéristiques géotechniques du
sol support prélevé et des échantillons ayant subi le cycle d’humidification et de séchage
afin de disposer d’une base complète pour l’analyse et la modélisation de leurs
comportements vis-à-vis des chaussées.
L’ensemble des essais conduits sur des échantillons de sols argileux est réparti en deux
catégories pour fournir les deux types de caractéristiques :
 Physiques
 Mécaniques.
III.1 Essais d’identification

III.1.1 Résultats des paramètres d’identification
L'ensemble des résultats d'identification simple sont résumés dans le Tableau III-1 :
Tableau III-1: Les résultats des essais d’identification
Issaba0.0-0.50 Issaba 0.50-1.0 Issaba 1.0-1.50 Issaba 1.50-2.0

Passant au
92 93 95 96
tamis 80µm (%)
Teneur en eau
51.53 56.58 46.1 42.41
naturelle w (%)
Limite de
liquidité (WL)% 92 93 93 95
Limite de
plasticité (WP) 51 52 52 57
Indice de
plasticité (IP) 41 41 41 38
Indice de
consistance 0.95 0.89 1.84 1.38
(Ic)
Matière
0.8 0.7 0.6 0.8
organique MO
(%)
Densité réelle
2.23 2.29 2.4 2.29
γs

Classe selon
A4 A4 A4 A3
GTR
Classe selon
A7-5 A7-5 A7-5 A7-5
HRB
III.1.2 Analyse et interprétation des résultats des essais d’identification

III.1.2.1 Teneur en eau in situ
Le Graphe III-1 présente les valeurs de teneur en eau in situ du sol prélevé à
différentes profondeurs dans la zone d’Issaba. La teneur en eau de cette localité varie
de 46.1% à 56.58%.
60 56.58
51.53
50 46.1
42.41
TENEUR EN EAU IN SITU
40
30
20
10
0
1
0-0.5 0.5-1 1-1.5 1.5-2
Graphe III-1 : les valeurs de teneur en eau in situ du sol prélevé à

différentes profondeurs dans la zone d’Issaba
III.1.2.2 Analyse granulométrique par tamisage et par sédimentométrie

Les courbes du Graphe III-2 représentent la répartition dimensionnelle (obtenue par
tamisage puis par sédimentométrie) des particules des sols étudiés.
Les courbes granulométriques (Graphe III-3) du sol d’Issaba à différentes profondeurs

nous montrent que ce sol a une granularité étalée et sa granularité varie en fonction de la
profondeur. L’analyse du Graphe III-3 nous permet de dire que le sol d’Issaba est un sol
fin car le pourcentage des passants au tamis 80 μm à différentes profondeurs est
supérieur à 35 %.

Courbes granulométriques NF EN ISO 17892-4 (2018)
Limon
GRAVIER Sable Argile
Grossier Moyen Fin Moyen Fin Grossie Moyen Fin
r
100
90
80
% Passants
70
ISSABA 0,00/0,50
60 ISSABA 0,50/1,00
50 ISSABA 1,00/1,50
ISSABA 1,50/2,00
40
30
20
10
0
20 2 0.2 0.02 0.002 0.0002
Ouverture des tamis (mm)
Graphe III-2: Récapitulatif des courbes granulométriques de l’argile d’Issaba
97
PASSANT AU TAMIS 80 MICRON
96
96
95
95
94
93
(%)
93
92
92
91
90
1
0-0.5 0.5-1 1-1.5 1.5-2
Graphe III-3:Valeurs du passant au tamis 80 micron (%)

III.1.2.3 Matière organique

Les valeurs des teneurs en matière organique obtenues sur le site d’Issaba à différentes
profondeurs et présentées dans le Tableau III-1 varient entre 0.6g/100g et 0.8g/100g.
Nous pouvons, en se basant sur le Tableau II-3, conclure que le sol d’Issaba est non
organique.
III.1.2.4 Limite d’Atterberg

Les résultats de l’essai sont présentés dans le Tableau III-1 et nous ont permis de tracer
les courbes de liquidité suivantes :
Graphe III-4:Récapitulatif des courbes de liquidité de l’argile d’Issaba de 0 à 2m

96 58 57
95
LIMITE DE PLASTICITE
LIMITE DE LIQUIDITE
95 56
Photo
94 III-1 54
93 93
93 52 52
92 52 51
92
50
91
90 48
1 1
0-0.5 0.5-1 1-1.5 1.5-2 0-0.5 0.5-1 1-1.5 1.5-2
Graphe III-6 :Valeurs de l’indice de liquidité Graphe III-5 Valeurs de la limite de plasticité
42
41 41 41
INDICE DE PLASTICITÉ (IP) %
41
40
39
38
38
37
36
1
0-0.5 0.5-1 1-1.5 1.5-2
Graphe III-7: Valeurs de l’indice de plasticité

Le récapitulatif des courbes de liquidité et de plasticité des Graphe III-5 et Graphe

III-6, nous montre que la limite de liquidité et de plasticité augmente en fonction de la
profondeur mais elle est restée constante de 0.5 à 1.5. La valeur de l’indice de plasticité
de chaque localité qu’exhibe le Graphe III-7, nous a permis de caractériser d’une part le
comportement plastique par le biais du diagramme de plasticité (Figure II-2), puis d’autre
part le potentiel de gonflement au moyen du Tableau II-2. Les valeurs de l’indice de
plasticité étant supérieures à 50% et la limite de liquidité étant supérieure à 60%, donc
l’argile d’Issaba à la profondeur 0 à 2m est très plastique et a un potentiel de gonflement
très élevé d’après les Tableau II-2.
Les différents paramètres mesurés au cours des essais d’identification ont permis
d’étayer de façon plus précise les caractéristiques physiques du sol prélevé à différentes
profondeurs. En effet, ces paramètres ont favorisé la connaissance de la nature et
d’autres propriétés à savoir la plasticité et le potentiel de gonflement.
III.2 Essai mécanique

Dans cette partie nous présenterons les résultats des essais de cisaillement et
œdométriques. Ces différents essais nous permettront de connaître les caractéristiques
mécaniques du sol argileux de la commune faisant l’objet de notre étude.
III.2.1 Essai de cisaillement

Les résultats provenant de l’essai de cisaillement précisent en occurrence la cohésion et
l’angle de frottement interne de ces sols argileux. Les résultats de l’essai sont présentés
dans les Tableau III-2 et nous ont permis de tracer les droites de Coulomb ci-dessous :

40.00 y = 0.0122x + 21.944

Série1 Linéaire (Série1) y = 0.0282x + 16.377 R² = 0.855
60.00 R² = 0.973 35.00
CONTRAINTE DE CISAILLEMENT (KPA)

CONTRAINTE DE CISAILLEMENT
50.00 30.00
25.00 Série1
40.00
20.00
(KPA)
30.00 Linéaire
15.00
(Série1)
20.00 10.00
5.00
10.00
0.00
0.00 0.00 100.00 200.00 300.00 400.00 500.00
0.00 100.00 200.00 300.00 400.00 500.00 CONTRAINTE NORMALE (KPA)
CONTRAINTE NORMALE (KPA)
Graphe III-8:Courbe intrinsèque de l’argile d’Issaba de 0 à 0 ,5 m Graphe III-9:Courbe intrinsèque de l’argile d’Issaba de 0,5 à 1,0 m
Série1 Linéaire (Série1) y = 0.0541x + 22.017 200.00 y = 0.0959x + 25.978

R² = 0.7645

180.00 R² = 0.8671
60.00
160.00
50.00 140.00 Série1

40.00 120.00
30.00 100.00
Linéaire
80.00 (Série1)
20.00
60.00
10.00 40.00
0.00 20.00
0.00 200.00 400.00 600.00 0.00
CONTRAINTE NORMALE (KPA) 0.00 100.00 200.00 300.00 400.00 500.00
CONTRAINTE NORMALE (KPA)
Graphe III-10:Courbe intrinsèque de l’argile d’Issaba de 1.0 à 1.5 m

Graphe III-11:Courbe intrinsèque de l’argile d’Issaba de 1.5 à 2 m

Tableau III-2:Récapitulatif des résultats qui émanent de l’essai de cisaillement
Profondeur W% W% 𝜸h 𝜸d Vitesse CU 𝝋U
(m) initiale finale (kN/m3) (kN/m3) (mm/min) (kPa) (dégrés)
0 – 0,5 51.53 50.07 17.08 11.27 1.5 16.4 1.6
0,5 – 1,0 49.59 50.42 16.70 11.17 1.5 21.9 0.7
1,0 – 1,5 46.10 47.97 16.94 11.60 1.5 22 3.1
1,5 – 2 42.73 42.73 17.25 12.11 1.5 26 5.5
Le Tableau III-2 nous révèle que le sol d’Issaba à la profondeur de 1.5 à 2.0m est plus
cohésif.
III.2.2 Essai œdométrique

Contraintes normales 1.0 Contraintes 100.0
10.0 normales 1000.0
1.0 10.0 100.0 1000.0
1.10 1.20
1.00 1.10
0.90 1.00
Indices des vides
Indices des vides
0.80 0.90
0.70 0.80
0.60 0.70
0.50 0.60
0.40 0.50
Courbe oedométrique Courbe oedométrique
Graphe III-13:Courbe de compressibilité de l’argile Graphe III-12:Courbe de compressibilité de l’argile
d’Issaba de 0,0 à 0,5 m d’Issaba de 0,5 à 1 m
Contraintes normales Contraintes normales

1.0 10.0 100.0 1000.0 10.0 100.0 1000.0 10000.0
1.25 1.04
1.15 1.02
Indices des vides
1.05 1.00
Indices des vides
0.98
0.95
0.96
0.85
0.94
0.75 0.92
0.65 0.90
Courbe oedométrique Courbe oedométrique
Graphe III-14:Courbe de compressibilité de l’argile Graphe III-15: Courbe de compressibilité de l'argile
d’Issaba de 1 à 1.5 m d'Issaba de 1.50 à 2.0m

Les valeurs issues de cet essai sont consignées dans le tableau ci-dessous.
Tableau III-3:Récapitulatif des valeurs d’essai œdométrique de l’argile d’Issaba
Profondeur
ei eo 𝝈′𝒗𝒐 𝝈′𝒑 cc cs
(m)
0 – 0,5 1.016 1.032 7.2 115 0.334 0.103
0,5 – 1,0 1.130 1.160 16.2 83 0.377 0.118
1,0 – 1,5 1.247 1.250 26.1 130 0.423 0.139
1,5 – 2 1.024 1.02 34.2 69.0 0.069 0.024
Du Tableau III-3, nous constatons que les contraintes de préconsolidation ơ’p sont
supérieures aux contraintes effectives verticales ơ’vo (ơ’p > ơ’vo) donc l’argile d’Issaba à
la profondeur 0 à 2 m est surconsolidée. De plus l’indice de compression (Cc) compris entre
0,3 < Cc < 0,5 nous permet de dire que l’argile d’Issaba est très compressible à la
profondeur 0 à 1.5 m d’après le Tableau II-5 et est peu compressible de 1.5 à 2m et
l’illite d’après le Tableau II-5. Enfin l’indice de gonflement (Cs) compris entre 0,04 < Cs
< 0,25 l’argile d’Issaba a un gonflement élevé d’après le Tableau II-7.
III.3 Estimation du gonflement par la méthode indirecte [21]

 D’après Seed et al ajusté par (Bekkouche et al 2001) nous pouvons estimer le
potentiel de gonflement εg (%) par :
𝜀𝑔 = 0.1823723 ∗ 𝐼𝑝 + 0.4559623
 D’après Komornik et David ajusté par (Bekkouche et al 2001) nous pouvons estimer
la pression de gonflement ơg (kPa) par :
log 𝜎𝑔 = 0,0079182 𝑤𝐿 + 0,00634245γd – 0,00161536w + 1,1328183

Tableau III-4: Valeurs de l’amplitude et de la pression de gonflement
Issaba
Profondeur (m) Potentiel de gonflement Pression de gonflement
εg (%) ơg (kPa)
0-0.5 7.93 70.72
0.5-1 7.93 70.57
1-1.5 7.93 76.59
1.5-2 7.38 85.87

III.4 Présentation des résultats des essais sur les échantillons ayant subi
le cycle d’humidification et de séchage
III.4.1 Déroulement de l’essai
Afin d'évaluer l'influence des différents cycles de séchage-mouillage sur le sol étudié,
nous avons procédé à la détermination des caractéristiques œdométriques et de
cisaillement des échantillons. Chaque échantillon a subi des cycles :
C0 : échantillon intact n'ayant pas subi de cycle
C1 : échantillon intact ayant subi un (01) cycle de séchage-mouillage
C2 : échantillon intact ayant subi deux (02) cycles de séchage-mouillage consécutifs
C3 : échantillon intact ayant subi trois (03) cycles de séchage-mouillage consécutifs
Des échantillons de chacun des cycles C1, C2, C3, C4, C5 ont été soumis à des essais
œdométrique et de cisaillement.
Cette partie nous présente les résultats des essais de cisaillement et œdométriques
effectués sur les éprouvettes ayant subi le cycle d’humidification et de séchage. Ici nous
allons faire une comparaison entre les résultats trouvés au cycle 1, 2 et 3 vu que les
échantillons du cycle 0 sont à l’état humide avant les essais. L’idéal serait de faire passer
ces échantillons par le séchage et qu’ils atteignent une teneur en eau de 14% comme les
autres échantillons du cycle 1, 2 et 3.
III.4.2 Texture et photos des échantillons en fonction du nombre de cycles

Cette partie nous présente la texture des échantillons à chaque cycle
 Cycle 0
a) b)
Photo III-2: échantillons du cycle 0

 Cycle 1
b) Après séchage
a) Après Humidification
Photo III-3: échantillons du cycle 1 a) et b)
 Cycle 2
a) Après humidification b) Après séchage
III.4.3 Photo III-4 : Les échantillons du cycle 2 a) et b)

 Cycle 3
a) Après Humidification
b) Après séchage
Photo III-5 : Les échantillons du cycle 3 a) et b)
III.4.4 Description des textures des échantillons au cours des cycles

Les échantillons du cycle 0 sont à l’état humide et ne présentent aucun gonflement ni
fissure. Au cours du cycle 1, après l’humification, nous notons un léger gonflement des
échantillons mais après séchage, nous constatons que les échantillons sont fissurés et
sont rétractés. Ce même constat est fait au deuxième et troisième cycle sauf que les
échantillons sont plus gonflants après réhumidification et présentent plus de fissures
après séchage par rapport au cycle 1. Ce qui montre que les échantillons se détériorent
au fur et à mesure qu’ils sont soumis aux cycles de retrait-gonflement.
III.4.5 Les résultats des essais de cisaillement

Le Tableau III-5 nous présente les résultats des essais de cisaillement. L’analyse des
Graphe III-17 et le Graphe III-16 nous révèle que, les valeurs de cohésion au niveau de
chaque profondeur ont augmenté du cycle 1 au cycle 2 et ont diminué du cycle 2 au cycle
3 alors que les valeurs de l’angle de frottement ont diminuées du cycle 1 au cycle 3.

Tableau III-5: Tableau récapitulatif des essais de cisaillement
Cycle
Cycle 0 Cycle 1 Cycle 2 Cycle 3
0- 0.5- 1- 1.5- 0- 0.5- 1- 1.5- 0- 0.5- 1- 1.5- 0- 0.5- 1- 1.5-

Profondeur
0.5m 1m 1.5m 2m 0.5m 1m 1.5m 2m 0.5m 1m 1.5m 2 0.5m 1m 1.5m 2m
CU (kPa) 16.4 21.9 22 26 24.5 15.7 9 20.3 26.6 27.7 18.9 25.4 17.5 9.5 16.0 7.8
𝝋U (dégrés) 1.6 0.7 3.1 5.5 21.5 26.1 24.9 24.4 10.1 14.5 12.2 15.3 5.34 9.7 12.2 10.8
GRAPHE ANGLE DE FROTTEMENT

GRAPHE COHESION 30
0.0 0.5 m
Angle de frottement
30 Cohésion 25
Cohésion kpa
20 0.5 1.0 m 20
0-0.5m
Cohésion 15
10 0.5-1m
1.0 1.5 m 10
Cohésion 1-1.5m
0 5
0 1 2 3 4 1.5 2.0 m 1.5-2m
Cohésion 0
Cycle C
0 1 2 3 4
Cycle C
Graphe III-17: Graphe de l'évolution de la cohésion Graphe III-16: Graphe sur l'évolution de l'angle de frottement au fil des
au fil des cycles cycles
III.4.6 Les résultats des essais œdométriques

Cycle
Eprouvette
0- 0.5- 1- 1.5- 0- 0.5- 1- 1.5- 0- 0.5- 1- 1.5- 0- 0.5- 1- 1.5-
0.5m 1m 1.5m 2m 0.5m 1m 1.5m 2m 0.5m 1m 1.5m 2 0.5m 1m 1.5m 2m
Paramètre
𝝈′𝒑 (Kpa) 115.0 83 130 69 1550 700 210 700 50 210 170 209 40 49 49 50
Eoed (MPa) 8.2 7 6.9 10 10 11.1 10.8 10.3 6.6 6.1 6.2 6.9 6.5 7.4 6.9 6.6
cc 0.334 0.377 0.423 0.366 0.271 0.290 0.315 0.323 0.48 0.42 0.41 0.37 0.45 0.40 0.39 0.34
cs 0.103 0.118 0.139 0.123 0.106 0.118 0.107 0.134 0.14 0.14 0.15 0.13 0.1 0.1 0.1 0.09
GRAPHE Indice de compressibilité GRAPHE Indice de gonflement

0.6 0.16
0.5 0.14
0.12
0.4
0-0.5m 0.1
0-0.5m
Cc
0.3
Cs
0.08
0.5-1m 0.5-1m
0.2 0.06
1-1.50m 0.04 1-1.5m
0.1
1.50-2m 0.02 1.50-2m
0 0
0 1 2 3 4 0 1 2 3 4
Cycle C Cycle C
Graphe III-19: Graphe de l'évolution de l’indice de compressibilité au fil Graphe III-18 : Graphe de l'évolution de l’indice de gonflement au fil des
des cycles cycles

Le Graphe III-19 montre une croissance de l’indice de compressibilité pour le sol lors du
passage du cycle 1 au cycle 2 et une décroissance du cycle 2 au cycle 3. En ce qui concerne
l’indice de gonflement, Le Graphe III-18 montre une augmentation de l'indice de
gonflement du cycle 1 au cycle 2 et une diminution de cet indice du cycle 2 au cycle 3. Il
est constaté généralement lors de la réhumidification du cycle 2 que le volume des
carottes augmente.
Le module oedométrique
12
10
8
0.0 0.5 m
Eoed
6
0.5 1.0 m
4 1.0 1.5 m
2 1.5-2.0 m
0
0 1 2 3 4
Cycle C
Graphe III-20:Courbe d'évolution du module œdométrique
Le Graphe III-20 traduisant l'évolution du module œdométrique montre une

décroissance du cycle 1 au cycle 3 à l’exception du sol à la profondeur 0.5 à 1.50m qui
décroit du cycle 1 à 2 et croit légèrement du cycle 2 à 3.

Conclusion
Un échantillon de sol prélevé à différentes profondeurs au Sud-est du Bénin précisément
dans la commune de Pobè à Issaba a été soumis à différents types d'essais afin de
caractériser les sols argileux de cette région, sensibles au phénomène de retrait-
gonflement. Les résultats d’essais d’identification révèlent que les matériaux prélevés à
Issaba sont très fins, plastiques et non organiques et ce en fonction de la profondeur. En
utilisant comme guide de classification, les classifications GTR et HBR, le sol d’Issaba est
de l’argile de classe A7-5 selon HRB et A4 selon le GTR. Les essais mécaniques nous
montrent que le sol du milieu d’étude est un sol argileux très compressible et surconsolidé
avec un pouvoir de gonflement élevé. Ceci lui confère sa propriété de sol gonflant en
présence de l’eau. Le cycle de séchage et mouillage nous a permis de mettre en évidence
l’effet du phénomène de retrait-gonflement sur les paramètres du sol.

Chapitre IV : Modélisation et prédiction du

comportement des sols expansifs sous léffet des
sollicitations cycliques : humidification et séchage.

Chapitre IV : Modélisation et prédiction du comportement des sols

expansifs sous léffet des sollicitations cycliques : humidification et
séchage.
Introduction
Ce chapitre est consacré à l’étude du comportement d’un système de chaussée reposant
sur un sol support gonflant, avec une élaboration d’un modèle numérique, conçu pour mieux
comprendre le mécanisme des déformations induites par le sol support sous l’influence de
l’effet des saisons pluvieuses et saisons sèches.
Pour bien s’approcher le plus près de l’état réel de sol, le modèle numérique est inspiré
d’une structure existante qui comprend un corps de chaussée de type semi-rigide et qui
repose sur un sol support gonflant, la chaussée a subi des dégradations majeures après
une année de sa conception. La section d’étude de la chaussée dégradée se trouve au
niveau de la route nationale RN3. Il s’agit de la route Pobè-Kétou. La modélisation a été
effectuée par le programme d’élément finis Plaxis, version 2012.
IV.1 Historique de la structure

Les travaux de construction de la route Pobè-Kétou ont été répartis en deux phases à
savoir :
-La phase n°1 : de la sortie de Pobè (PK 0+000) jusqu’à Onigbolo (PK 25+000).
-La phase n°2 : d’Onigbolo (PK 25+000) à Kétou (PK 46+000).
La route Pobè-Kétou a été construite en route bitumée standard avec une chaussée
bitumée de 7,0 m de largeur et des accotements non bitumés de 2×1,5 m, avec des tracés
horizontal et vertical conformes aux exigences d’une vitesse de référence de 100 km/h ;
les ouvrages d’assainissement étaient des dalots et ponceaux en béton armé. Pour
empêcher la remontée capillaire de l’eau, il a été posé sur la plateforme, un feuillard de
type « polyfelt TS 400 » fait en matière plastique spéciale avec un coefficient de
perméabilité K=3,9.10-2 et une charge surfacique de 0,1 kgf/cm2. La pose du feuillard a
été immédiatement suivie d’une couche de fondation en latérite de 20 cm d’épaisseur
compactée à 95% de l’Optimum Proctor Modifié (OPM). La couche de fondation a été
suivie d’une couche de base améliorée au ciment à concurrence de 4% sur les sections
allant du PK 0+150 au PK 3+500 puis du PK 18+500 au PK 25+00, le tout sur une épaisseur
de 25 cm. En ce qui concerne le tronçon du PK 3+500 au PK 18+500, les 25 cm de couche
de base ont été confectionnés en concassé de 0/60 à partir des roches calcaires se
trouvant aux abords immédiats de l’usine d’Onigbolo ; cette couche a été suivie par une
couche d’imprégnation de 0,8 kg d’émulsion à 65% de bitume. La compacité de la couche
de base a été fixée à 97% de l’OPM. La couche de revêtement a été faite en enrobés fins.

Malgré toutes ces précautions prises lors de la construction de la chaussée, de

nombreuses dégradations sont apparues un an après sa mise en service, montrant ainsi
l’inadéquation de la technique utilisée par rapport au type de sol rencontré dans la zone.
L’état des lieux avait donc permis de répartir les travaux en fonction des caractéristiques
du sol le long du projet et de l’état du revêtement en place. Pour les travaux, la route a
donc été scindée en quatre (4) sections :
 Section 1 : Pobè PK 0+000 au Carrefour Obélé PK 2+480 soit 2,480 km

 Section 2 : Carrefour Obélé PK 2+480 à Onigbolo PK 18+860 soit 16,380 km
 Section 3 : Onigbolo PK 18+860 à Odomèta PK 24+350 soit 5,490 km
 Section 4 : Odomèta PK 24+350 à Kétou PK 43+135 soit 18,785 km.
Notre étude s’intéresse à la deuxième section qui présente un état dramatique dont en
lieu et place d’une route bitumée on a une route en terre avec un revêtement presque
inexistant.
Par conséquent, durant les travaux de réhabilitation, les solutions proposées sur cette
section sont :
 La reconstruction sur la section de sol argileux expansif avec mise en place de

dispositions pour réduire les éventuels effets négatifs du gonflement des
matériaux argileux et plastiques du sol support du tronçon
 Aménagement d’accotements revêtus
A cet effet la structure de la nouvelle chaussée de la section 2 est constituée d’une

couche de roulement de 0.07 m en béton bitumineux, une couche de base en grave
latéritique de 0.15 m, une couche de fondation en grave concassée de 0.20 m et une couche
de forme en grave latéritique de 0.20 m [11].
IV.2 Présentation du logiciel Plaxis

IV.2.1 Généralité sur Plaxis
Plaxis est un logiciel d’éléments finis d’analyses de contraintes et de déformations dans
les massifs de sols sous chargement. Il permet d’effectuer une analyse des éléments finis
bidimensionnels de déformation et de stabilité dans le génie géotechnique et la mécanique
rocheuse.
Les entreprises et institutions du génie civil option géotechnique s’appuient sur Plaxis pour
différents projets, notamment des fouilles, des digues et des fondations, sans oublier les
tunnels, le pétrole et le gaz, les exploitations minières et la géomécanique des réservoirs.

L’ensemble des options simplifiées (initialisation des contraintes, pressions

interstitielles) permettent d’aller au but (prévoir le comportement d’un ouvrage), quitte
à réaliser ultérieurement, avec le même code et les mêmes données, un calcul affiné.
IV.2.2 Fonctionnalités principales du logiciel sont :

 Étalonnez avec précision les modèles de matériaux
 Interopérabilité avancée avec l’écosystème Bentley
 Automatisez les tâches pour gagner en efficacité grâce aux scripts Python
 Importez des fichiers CAO pour une modélisation simplifiée, afin de gagner du
temps
 Renforcez la fiabilité grâce à une bibliothèque de modèles constitutifs
prééminente
 Accédez à davantage de fonctionnalités grâce à l’analyse de sensibilité et à la
variation des paramètres.
IV.2.3 Les points forts de Plaxis 2D

Le logiciel PLAXIS 2D effectue une analyse des éléments finis bidimensionnels de
déformation et de stabilité dans le génie géotechnique et la mécanique rocheuse. Ces
points forts sont les suivants :
 Convivialité de l’interface (Windows©), pour la saisie des données, la définition des

phases de projets et l’analyse des résultats ;
 Génération automatique des maillages ;
 Des jeux prédéfinis de lois de comportement de sol (linéaire élastique, Mohr-
Coulomb, Hardening Soil, Soft Soil, …) et la possibilité de définir ses propres lois
de comportement ;
 Couplage avec les calculs d’écoulement et de consolidation ;
 Prise en compte des éléments de structure et de l’interaction sol-structure ;
 Calcul de coefficients de sécurité (procédure c/phi réduction).
IV.2.4 Les options par défaut et les solutions approchés

Le système d’options par défaut et de solutions approchées spécifiques, qui sont un des
faits de l’outil pour propulser le projet de la géotechnique, est destiné à faire gagner du
temps à l’opérateur, à lui éviter de devoir faire des choix tracassant, et enfin à améliorer
la convivialité du logiciel. Ce système est inséparable du traitement à partir d’un menu
arborescent. Chaque branche du menu réalise une tâche précise, bien définie et est
évidemment figée. Mais la diversité des branches rend l’outil extrêmement souple.
IV.2.4.1 Le maillage
Les options par défaut commencent dès le maillage : si seules les grandes lignes de celui-
ci importent, le détail des éléments, agencé de manière optimale du point de vue

numérique, sera entièrement généré par le logiciel à partir d’un petit nombre de nœuds-
clé, avec contrôle permanent à l’écran (le maillage généré automatiquement : ce qui est un
point fort de Plaxis), mais l’opérateur peut aussi bien entendu spécifier un maillage très
détaillé : paramétrer la finesse du maillage entre différentes options (très grossier,
grossier, moyen, fin, très fin), l'opérateur peut également décider de mailler plus
finement une certaine région ou/et le voisinage d'un élément grâce aux option « refine »
dans le « mesh » menu.
IV.2.4.2 Les conditions aux limites

Il faut entrer les conditions limites, c'est-à-dire les déplacements et les contraintes
imposées aux limites extérieures de la géométrie. Si aucune condition limite n'est fixée
sur un tronçon, par défaut le logiciel considère que l'élément n'est soumis à aucune force
extérieure et est libre de se déplacer dans toutes les directions. Les conditions limites
pouvant être imposées sont celles qui fixent un déplacement dans une direction donnée ou
celles qui imposent une force dans une direction donnée. Plusieurs outils permettent de
créer une large gamme de conditions limites (force repartie, force ponctuelle,
encastrement, glissement, etc.). De même en ce qui concerne les conditions aux limites en
déplacements : si elles sont complexes, l’ingénieur géotechnicien devra en spécifier les
subtilités d’une manière précise, face de bloc par face de bloc. Par contre, si elles ont un
caractère standard (vecteur déplacement nul à la base du domaine étudié et vecteur
déplacement horizontal nul sur ses faces latérales), l’application peut les réaliser
automatiquement (par défaut) à partir du menu avec contrôle immédiat du résultat à
l’écran. L’application des contraintes initiales dues au poids des terres peut être
matérialisée de manière exacte par activation du multiplicateur de chargement relatif au
poids propre. Par contre, si comme bien souvent en géotechnique on connaît ou on sait
estimer un état K0 (coefficient des terres au repos) donné, celui-ci peut être spécifié
directement. Dans ce cas, le massif est souvent en léger déséquilibre (incompatibilité
entre K0 et les autres caractéristiques mécaniques). Le menu permet alors, par un
changement fictif nul, de rééquilibrer le massif, puis de réinitialiser à zéro le champ de
déplacement de manière à prendre comme nouvelle origine l’état du matériau après
application de la gravité. L’option K0 est particulièrement intéressante-et réaliste- dans
le cas d’un modèle hétérogène de surface libre presque horizontale (paroi moulée dans un
sol mou par exemple).
IV.2.4.3 Les pressions interstitielles

Les pressions interstitielles ont été l’objet d’un soin particulier dans Plaxis. Pour qui
souhaite faire un calcul précis du champ de pressions interstitielles en régimes permanent
ou transitoire, c’est possible grâce au module d’écoulements en milieu poreux. Mais bien
sûr, cette opération demande du temps (d’opérateur et de machine). Si la nappe
phréatique n’est pas trop éloignée de l’horizontale, dans ses états initial et final, on sait

que la pression diffère peu de la pression hydrostatique ; si l’on adopte ce champ de

pression approchée, les calculs deviennent très simples puisqu’il s’agit seulement de manier
les variations de la poussée d’Archimède ; Plaxis offre cette possibilité qui est souvent
très appréciable.
IV.2.4.4 La conduite des calculs non linéaires

La conduite des calculs non linéaires constitue un autre exemple de la souplesse
d’utilisation que procure ce logiciel : l’opérateur peut évidemment faire lui-même ses choix
de la taille d’étape de chargement, du nombre d’étapes, de rigidité d’interface, du nombre
de nœud, de la méthode de résolution, … ; s’il ne désire pas assumer ces choix, le logiciel
peut les décider à sa place, compte tenu de l’expérience des numériciens en la matière.
Pour les calculs de consolidation, réalisés en différence finies explicites sur le temps, le
choix du pas de temps peut également être décidé par l’utilisateur, ou bien calculé dans
l’option par défaut, selon les critères numériques connus.
IV.2.4.5 Le coefficient de sécurité

Le coefficient de sécurité résume en une seule information une quantité considérable de
données expérimentales ou non. L’approche classique évalue généralement ce nombre selon
la théorie de l’équilibre limite, supposant une réduction proportionnelle généralisée de la
résistance mécanique des matériaux impliqués, ce qui ne constitue pas un scénario réel de
rupture. C’est la même approche, adaptée aux éléments finis élasto-plastiques, qui conduit
à l’évaluation du coefficient de sécurité dans Plaxis. Le critère de « rupture » est ici
qualitatif, et laissé à l’appréciation de l’observateur ; en tout état de cause, il est fondé
sur le niveau de déplacement d’un point de contrôle lié à l’ouvrage étudié. Le champ de
déplacement obtenu est évidemment tout à fait fictif.
IV.2.4.6 Les modèles de comportement utilisés dans Plaxis

Les modèles de comportement de sols sont très nombreux : depuis le modèle élastique
plastique de Mohr-Coulomb jusqu’aux lois de comportement les plus sophistiquées
permettant de décrire presque tous les aspects du comportement élasto-visco-plastique
des sols, aussi bien sous sollicitation monotone que cyclique. Ces modèles ont été
développés dans le but d’être intégrés dans les calculs par éléments finis. La modélisation
par éléments finis permet de résoudre le problème aux limites en tenant compte, par une
loi de comportement réaliste, du comportement réel du sol.
Du fait que l’un des objectifs de Plaxis est de fournir à l’utilisateur un code d’éléments
finis qui soit à la fois robuste et convivial, permettant de traiter des problèmes
géotechniques réels, dans un délai raisonnable en utilisant des modèles de comportement
de sols dont les paramètres puissent être déterminés à partir d’une étude géotechnique
normale. La démarche suivie dans le développement du code Plaxis est différente ; en ce
sens que Plaxis peut apparaître comme une « règle à calcul » de l’ingénieur géotechnicien,

ou le micro-ordinateur a remplacé la règle. C’est pourquoi les différents modèles de

comportement utilisés dans Plaxis sont des modèles qui peuvent apparaître simples, voire
simplistes, mais qui sont efficients quand ils sont utilisés dans des cas adaptés.
Pour traiter un problème de soutènement (paroi moulée, palplanche, etc.…), il est tout à
fait adapté de considérer le sol comme élastoplastique et le modèle de Mohr-Coulomb
sera bien adapté dans ce cas : on rejoint ici le calcul des soutènements par les méthodes
élasto-plastiques de coefficient de raideur. Mais pour traiter une construction de remblai
sur sols mous, avec chargement par étapes et consolidation, il faut tenir compte de
l’écrouissage. Le matériau se consolide et il est plus adapté d’utiliser le Soft-Soil-Model
qui prend en compte cette évolution du matériau. Pour un calcul d’écoulement, il suffit de
prendre un matériau élastique, mais on peut avoir à coupler écoulement et déformation :
dans ce cas un modèle élasto-plastique peut-être justifier.
IV.2.4.7 Comportement élasto-plastique

Un matériau élasto-plastique est caractérisé par un seuil d’écoulement au-dessous duquel
le comportement est purement élastique, alors qu’une fois celui-ci atteint la déformation
n’est plus réversible. En toute rigueur, cette dernière est indépendante du temps : elle
n’évolue pas si la charge est maintenue constante, elle ne disparaît pas progressivement
après décharge. De même, le seuil n’est pas fonction de la vitesse de déformation. Le
comportement élasto-plastique peut être représenté par un modèle monodimensionnel, en
série un ressort de raideur K, pour symboliser l’élasticité du matériau, à un patin de seuil
S0 (Figure IV-1).
Figure IV-1 : Modèle monodimensionnel du comportement élastoplastique.
La courbe effort-déplacement ou contrainte-déformation que l’on trouve est présentée

sur la Figure IV-2

Figure IV-2: Représentation du comportement élastique-parfaitement plastique.
Lors d’une décharge, le comportement est élastique et réversible. La longueur de la

déformation plastique est a priori indéterminée. Le type de comportement représenté par
les Figure IV-1 et Figure IV-2 est un comportement élastique- plastique sans
écrouissage.
IV.2.4.8 Modèle élastique linéaire

Le modèle élastique utilisé est classique. Les tableaux de rentrée des données dans Plaxis
(Figure 5-8) demandent le module de Young E, le coefficient de poisson v, le module de
cisaillement G et le module œdométriques. L’avantage de G est d’être indépendant des
conditions de drainage du matériau (Gu=G’), ce qui n’est pas le cas des modules d’Young :
le module de Young non drainé est supérieur au module de Young drainé. Il aurait pu
sembler logique, si G est utilisé comme paramètre élastique, d’utiliser K comme second
paramètre. D’une part Ku est infini (correspondant à νur =0.5) et il est moins courant
d’emploi. Le module de cisaillement G est en fait le module mesuré dans les essais
pressiométriques (Brinkgereve R.B.J. & Vermeer P.A). D’autre part il existe également
des relations qui lient les modules : on passe de G à E par la relation : E = 2 G (1+ν) et du
module œdométrique au module d’élasticité drainé E’ par la relation :
(1 + 𝑣 ′ )(1 − 2𝑣 ′ )
𝐸 = 𝐸𝑜𝑒𝑑 ( IV-1)
(1 − 𝑣 ′ )
Le modèle élastique de Plaxis peut être employé surtout pour modéliser les éléments de
structures béton ou métal en interaction avec le sol. Il peut aussi être intéressant pour
certains problèmes de mécanique des roches.

Figure IV-3: Interface de saisie des paramètres du sol dans le modèle

élastique linéaire
IV.2.4.9 Modèle de Mohr-Coulomb

Le comportement de Mohr-Coulomb présente un comportement élastique parfaitement
plastique sans écrouissage. Il a une grande utilisation dans la géotechnique vu les résultats
obtenus dans les calculs. Dans le plan de Mohr, la droite intrinsèque est représentée par:
𝜏 = 𝜎𝑛 tan(φ) + 𝐶 ( IV-2)
Où 𝜎𝑛 et τ sont respectivement les contraintes normales et de cisaillement, et C et φ

respectivement la cohésion et l’angle de frottement du matériau.

Figure IV-4:Courbe intrinsèque du modèle de Mohr-Coulomb
Le modèle de Mohr-Coulomb demande la détermination de cinq paramètres (Figure IV-5).

Les deux premiers sont E et ν (paramètres d’élasticité). Les trois autres sont c, φ et ψ
respectivement, la cohésion, l’angle de frottement et l’angle de dilatance. Ce sont des
paramètres classiques de la géotechnique, certes sont souvent fournis par des essais de
laboratoires, mais nécessaires aux calculs de déformation ou de stabilité.

Mohr-Coulomb

IV.2.4.10 Modèle pour sols mous « Soft Soil Model (S.S.M) »

Le modèle pour sol mou (en abrégé SSM) est un modèle dérivé du Cam-Clay.
Historiquement le modèle Cam-Clay a été développé à Cambridge dans les années 60 par
Roscoe, Schoffield et al. L’idée de base de ce modèle est de prendre en compte l’effet
d’écrouissage que provoque sur les argiles la pression moyenne. Sous l’effet d’une pression
moyenne, la teneur en eau diminue et l’argile devient plus résistante. Il s’agit d’un modèle
élasto-plastique avec une surface de charge. Sous la surface de charge, le matériau reste
élastique, tandis que si le point représentatif de l’état de contrainte effectif atteint non
réversible. Une surface de plasticité, associée, limite l’espace entre les états admissibles
et non admissibles.
NB : Le modèle Soft Soil Model (SSM) n’est pas valable pour des interfaces : il est
nécessaire d’utiliser un modèle de Mohr-Coulomb. De plus ce modèle ne permet pas des
analyses de stabilité par réduction de c et φ.
Figure IV-6:Interface de saisie des paramètres du sol dans le modèle

Soft Soil Model

IV.2.4.11 Modèle de sol avec écrouissage « Hardening Soil Model (H.S.M) »

Le modèle Hardening Soil (HSM) a pour objet d’améliorer le modèle de Mohr-Coulomb sur
différents points, il s’agit essentiellement :
 De prendre en compte l’évolution du module de déformation lorsque la contrainte

augmente: les courbes œdométriques en contraintes-déformation ne sont pas des
droites.
 De prendre en compte l’évolution non linéaire du module lorsque le cisaillement
augmente : le module E50 n’est pas réaliste : il y a une courbure des courbes effort-
déformation avant d’atteindre la plasticité.
 De tenir compte de la dilatance qui n’est pas indéfinie.
On pourrait dire que ce modèle est un dérivé du modèle hyperbolique de Duncan-Chang

car il en reprend et améliore les formulations hyperboliques des courbes effort-
déformation.

Hardening Soil

IV.2.4.12Modèle défini par l’utilisateur

User-Defined Model permet de définir et d'utiliser des lois de comportement autres que
les modèles standard de PLAXIS. Pour une description plus détaillée de cette option,
l’utilisateur de Plaxis peut consulter le manuel des modèles de matériaux (Material Models
Manual).
IV.2.5 Les paramètres fondamentaux de Plaxis

IV.2.5.1 Le module de Young
Le module d'Young est utilisé dans Plaxis comme module de déformation de référence
dans le modèle élastique et le modèle de Mohr-Coulomb, mais d’autres modules de
déformation sont également considérés. Le choix d’un module de déformation est un des
problèmes les plus difficiles en géotechnique. Le module de déformation varie en fonction
de la déformation et en fonction de la contrainte moyenne. Dans le modèle de Mohr-
Coulomb, le module est constant. Il parait peu réaliste de considérer un module tangent à
l’origine (ce qui correspondait au Gmax mesuré dans des essais dynamiques ou en très
faibles déformations). Il est conseillé de prendre un module « moyen », par exemple celui
correspondant à un niveau de 50% du déviateur de rupture car en mécanique des sols, la
pente initiale est appelée Eo et le module sécant à 50% de la résistance en compression
est noté Eo (Figure IV-8). Pour des argiles très consolidées et quelques roches avec un
large domaine élastique, il est réaliste d’utiliser Eo alors que pour les sables et les argiles
normalement consolidées, il est préférable de prendre E50. L’utilisateur doit rester
conscient de l’importance du choix du module qu’il prendra en compte.
Figure IV-8: Définition du module initiale Eo et du module à 50% de la

rupture E50.
Le module initial et le module sécant tendent à augmenter avec la pression de confinement
pour les sols. Par conséquent, les couches de sol en profondeur ont souvent une raideur
plus élevée que les couches en surface. De plus, la raideur apparente dépend du chemin
de contrainte suivi. La raideur durant les cycles charge-décharge est supérieure à celle
durant un chargement primaire. En outre, la raideur observée pour un sol en termes de
module d’Young est généralement plus petite pour des compressions drainées que pour du
cisaillement. Donc, en choisissant un module de déformation constant pour représenter le
comportement du sol, l’utilisateur devra retenir une valeur qui prenne à la fois en compte

le niveau de contrainte et le chemin de contrainte. Remarquons que des comportements

de sols dépendant des contraintes sont pris en compte dans les modèles avancés de Plaxis.
Dans les boites de dialogues, on peut aussi rentrer un gradient donnant la variation du
module avec la profondeur.
IV.2.5.2 Le coefficient de poisson

Dans des cas, comme des problèmes de déchargement particuliers, il peut être réaliste
d’utiliser une valeur initiale aussi faible du coefficient de poisson (ν), mais en général,
pour le modèle de Mohr-Coulomb, l’utilisation d’une valeur plus élevée est recommandée.
Le choix d’une valeur pour le coefficient de Poisson est particulièrement simple dans le
cas du modèle élastique ou du modèle de Mohr-Coulomb avec un chargement gravitaire.
Dans la plupart des cas, les valeurs de ν seront comprises entre 0,3 et 0,4. En général, de
telles valeurs peuvent être aussi utilisées pour des conditions de chargement autres que
la compression unidimensionnelle. On conseille une valeur de 0,2 à 0,4 pour le coefficient
de poisson. Celle-ci est réaliste pour l’application du poids propre (procédure K0 ou
chargement gravitaires). Pour certains problèmes, notamment en décharge, on peut
utiliser des valeurs plus faibles. Pour des sols incompressibles, le coefficient de poisson
s’approche de 0.5 sans cette valeur soit utilisable. Dans le cas d’un comportement non
drainé, il est conseillé d’entrer une valeur effective pour le coefficient de Poisson et de
sélectionner Undrained comme type de comportement. De cette façon, PLAXIS ajoutera
automatiquement un module de compressibilité pour le fluide interstitiel, basé sur un
coefficient de Poisson non drainé, implicite de 0,495. Dans ce cas, le coefficient de
Poisson effectif devra être inférieur à 0,35. L’utilisation des valeurs plus élevées du
coefficient de Poisson impliquerait que l’eau ne soit pas suffisamment raide par rapport
au squelette du sol pour simuler un comportement non drainé.
IV.2.5.3 L’angle de frottement

Plaxis ne prend pas en compte une variation d’angle de frottement avec la contrainte
moyenne. L’angle de frottement à introduire est soit l’angle de frottement « de pic » soit
l’angle de frottement de palier. On attire l’attention sur le fait que des angles de
frottement supérieurs à 45° peuvent considérablement allonger les temps de calcul. Il
peut être avisé de commencer des calculs avec des valeurs raisonnables d’angle de
frottement, quitte à les augmenter dans la suite. Cette valeur de 45° est compatible avec
les angles de frottement φcv (à volume constant, au palier). On peut déterminer l’angle
de frottement à partir de la courbe intrinsèque du modèle de Mohr-Coulomb
IV.2.5.4 La cohésion
Dans Plaxis, il est utile d’attribuer, même à des matériaux purement flottants, une très
faible cohésion (0.2 à 1 KPa ou prendre c>0,02 kPa) pour des questions numériques afin

d’éviter des complications aux utilisateurs peu expérimentés car PLAXIS peut manipuler
des sables sans cohésion (c=0), mais certaines options ne fonctionneront pas bien.
De plus pour les analyses en comportement non drainé avec υu=0, Plaxis offre l’option de
faire varier la cohésion non drainée avec la profondeur : ceci correspond à la croissance
linéaire de la cohésion en fonction de la profondeur observée dans des profils au
scissomètre ou en résistance de pointe de pénétromètre. Cette option est réalisée avec
le paramètre c-depth. Une valeur nulle donne une cohésion constante. Les unités doivent
être homogènes avec ce qui a été choisi dans le problème (typiquement en KPa/m).
IV.2.5.5 L’angle de dilatance (Ψ)

L’angle de « dilatance » noté Ψ est donné en degré. Les sols argileux ne présentent aucune
dilatance (ψ=0) sauf pour des couches très surconsolidées. La dilatance d’un sable dépend
de sa densité et de son angle de frottement. Pour des sables siliceux, un ordre de
grandeur est ψ = φ-30 ; c’est le paramètre le moins courant. Il peut cependant être
facilement évalué par la règle (grossière) suivante :
ψ = φ − 30 pour φ > 30°

{ Équation IV-3
ψ = 0 pour φ < 30°
Les cas où Ψ correspond à des sables très lâches (état souvent dit métastable, ou
liquéfaction statique). La valeur Ψ=0 correspond à un matériau élastique parfaitement
plastique, ou il n’y a donc pas de dilatance lorsque le matériau atteint la plasticité. C’est
souvent le cas pour les argiles ou pour les sables de densité faible ou moyenne sous
contraintes assez fortes.
IV.2.5.6 Les paramètres de compressibilité

L’indice de compressibilité Cc et l’indice de gonflement Cs décrivent le comportement
œdométrique ou isotrope observé dans des essais de laboratoires : ce seront les deux
paramètres de base réglant la position des lignes de consolidation vierge ou des lignes de
gonflement. L’axe des contraintes est tracé en logarithme, ce qui conduit à modifier la
définition de Cc et Cs en λ et K dans ce cas, l’axe des ordonnées est l’indice des vides. Il
peut être judicieux de remplacer l’indice des vides par la déformation volumique (identique
à la déformation axiale dans l’essai œdométrique). On utilise alors λ * et k *. Dans ce cas,
l’indice des vides, variable, apparaît dans les expressions du Tableau IV-1Il peut dans la
majorité des cas être pris constant, et égal à la valeur initiale. Les relations avec les
paramètres du modèle de Cam-Clay sont indiquées dans le Tableau IV-1

Tableau IV-1: Valeurs des paramètres de compressibilité et de gonflement λ et k
IV.2.5.7
IV.2.5.8 La Prise en compte de la surconsolidation

Plaxis calcule la pression de préconsolidation ơ’p à partir des données. Si un matériau est
surconsolidé, il est possible de prendre en compte le degré de surconsolidation (OCR) qui
est le rapport de la pression de préconsolidation et la contrainte effective de surcharge
ou le poids des terres de préconsolidation (POP) (Pre Overburden Pressure).
IV.2.5.9 Les paramètres avancés

Pour tenir compte des angles beta, alpha, coefficient de perméabilité et des variations de
l’indice des vides, on utilise les paramètres avancés
Figure IV-9:Interface de saisie des paramètres avancé du sol

IV.3 La modélisation numérique de la chaussée

IV.3.1 Géométrie du modèle
Le modèle géométrique représente la structure de la demie chaussée. Les dimensions des
couches de structure sont les suivantes :
Epaisseurs Matériaux
0,07 m Béton bitumineux (BB)
0,15 m Gravier (GNT)
0,20 m Grave ciment (GC3)
0,20 m Grave latéritique (GL)
La structure de la demie chaussée a été modélisée par un modèle géométrique plan (2D).
Pour mieux observer la distribution des déformations tout au long de l’emprise de la
chaussée et en profondeur, nous avons choisir 10 m de largeur sur 10 m de profondeur
pour notre modèle géométrique. Ce modèle est présenté sur la Figure IV-10.
Figure IV-10:Modèle géométrique de la structure
IV.3.2 Caractéristique du sol

Le chapitre 3 nous révèle que le phénomène de retrait et gonflement influence les
paramètres du sol expansif. Le but de cette partie de notre étude est d’observer le

comportement de la chaussée sous la variation des paramètres du sol. Notre étude

consiste donc à faire varier les paramètres du sol en maintenant fixe les caractéristiques
de la chaussée. Pour cela, nous allons utiliser les résultats trouvés des essais pour chaque
cycle comme valeurs d’entrées dans Plaxis et les caractéristiques de la nouvelle structure
de la section 2 décrit dans l’historique du tronçon Pobè-Kétou. Une des hypothèses de ce
travail est de considérer que le sol support est régi par une loi élastoplastique
représentée par le modèle Mohr Coulomb drainé. Les paramètres de ces modèles utilisés
dans la simulation sont rassemblés dans le Tableau IV-2:

Tableau IV-2:Paramètres caractérisant l’argile gonflant utilisée pour le modèle Mohr-coulomb non drainé
Cycle
Eprouvette
0- 0.5- 1- 1.5- 0- 0.5- 1- 1.5- 0- 0.5- 1- 0- 0.5- 1- 1.5-
1.5-2
0.5m 1m 1.5m 2m 0.5m 1m 1.5m 2m 0.5m 1m 1.5m 0.5m 1m 1.5m 2m
Paramètre
Poids volumique
saturé (KN /m3 17.08 16.70 16.94 17.25 14.13 14.30 14.90 15.05 13.4 13.2 12.5 14 14.95 14.95 14.95 14.95
)
Poids volumique
non saturé (KN 11.27 11.17 11.60 12.11 12.27 12.57 13.01 13.1 11.11 11.2 10 12.3 11.90 11.90 11.90 11.90
/m3 )
Module de
Young E (KN 5110 4360 4300 6230 6230 6920 4240 6400 4112 3800 3863 4299 4050 4610 4299 4112
2
/m )
Coefficient de
0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0.35
Poisson ν
CU (kPa) 16.4 21.9 22 26 24.5 15.7 9 20.3 26.6 27.7 18.9 25.4 17.5 9.5 16 7.8
𝝋U (dégrés) 1.6 0.7 3.1 5.5 21.5 26.1 24.9 24.4 10.1 14.5 12.2 15.3 5.34 9.7 12.17 10.75
Angle de
0
dilatation (°)
Facteur de
rigidité de 0.6
l’interface Rinte

Coefficient de
perméabilité Kx 0.012
(m/jour)
Coefficient de
perméabilité Ky 0.012
(m/jour)

IV.3.3 Caractéristiques des matériaux utilisés en couche de chaussée

Dans cette étude les couches de chaussée sont supposées suivre un comportement
linéaire. Le modèle utilisé pour le calcul de déformation pour ces couches est le modèle
Linear Elastic de Plaxis.
o Caractéristiques des graves traitées au ciment utilisées en couche de fondation
Les propriétés de grave ciment utilisées dans les calculs sont rassemblées dans le tableau
ci-dessous :
Tableau IV-3: Caractéristiques de graves latéritiques traitées au ciment pour Plaxis
Paramètre Valeur
Type de comportement Elastique linéaire
2
Module de Young E (kN/m ) 13800000
Coefficient de Poisson ν 0.25
o Caractéristiques des graves non traitées utilisées en couche de base
Tableau IV-4:Caractéristiques de graves non traitées pour Plaxis
Paramètre Paramètre
2
o Caractéristiques du béton bitumineux utilisé en couche de base
Tableau IV-5: Caractéristiques du béton bitumineux utilisé en couche de base
Paramètre Valeur
2
IV.3.3.1 Génération du maillage

Après la définition de la géométrie de la structure avec toutes les propriétés des
matériaux constituants les différentes couches, la structure est subdivisée en des
éléments fins afin de procéder au calcul par éléments finis. L’ensemble des couches de la
structure a été discrétisé dans son intégralité par des éléments finis de type triangulaire
à 15 nœuds.

Figure IV-11:Maillage du modèle
IV.3.3.2 Définition des conditions initiales

Les conditions initiales sont constituées de trois modes différents, l'un pour générer les
pressions interstitielles initiales (mode des conditions hydrauliques) et l'autre pour
spécifier la configuration géométrique initiale et générer le champ des contraintes
effectives initiales (mode de configuration géométrique).
IV.3.3.2.1 Conditions hydrauliques

La nappe phréatique est supposée située à 2 m de la surface de la plate-forme (Figure
IV-13).
IV.3.3.2.2 Pressions interstitielles initiales

La distribution de pressions interstitielles est très spécifique, elle est générée
automatiquement avec l’option (Phreatic level) (Figure IV-14)

Figure IV-13 : Génération des pressions

Figure IV-14: Conditions hydraulique initiale
interstitielles initiales
IV.3.3.2.3 Contrainte initiale

Pour le calcul des contraintes initiales, on génère les contraintes initiales tout en nous
assurant que les valeurs de Ko générées par le logiciel sont supérieures à 0.5 vu que notre
sol est de l’argile surconsolidée. Le Tableau IV-6 présente quelques valeurs de
coefficient Ko en fonction du type de sol. On garde le poids du sol à 1, ce qui correspond
à une application totale de la gravité.
Tableau IV-6: Valeurs de Ko pour quelques sols[27]
Type de sol Valeur de Ko

Sable lâche 0.45 à 0.50
Sable compact 0.40 à 0.45
Argile normalement consolidée 0.50
Argile surconsolidée  0.50

Figure IV-15 Contrainte effective

IV.3.4 Procédure de calcul

Pour chaque cycle, le calcul du modèle se fait en 3 phases dans l’ordre, passant :
o Phase 0 : Initiation des contraintes (procédure K0) ; on détermine les contraintes

effectives initiales.
o Phase 1 : Désactivation de la charge et des couches de chaussée (Figure IV-16)
Figure IV-16:Désactivation de la charge
o Phase 2 : Activation des couches de chaussée (Figure IV-17)
Dans la phase 2 de calcul on choisit d’appliquer la charge de référence de la méthode

française de dimensionnement.
o Phase 3 : Activation de la charge (Figure IV-18)
En effet, elle est décrite à l’aide de deux disques de 0,125 m de rayon, d’entre-axe 0,375m
et exerçant en surface de chaussée une pression verticale et statique uniformément
répartie de 0,662 MPa.

Figure IV-17: Activation des couches de chaussée
Figure IV-18: Activation de la charge

IV.4 Les résultats de calcul

Comme notre travail consiste à observer le comportement de la chaussée sous la variation
des paramètres du sol au fil des cycles, nous allons procéder cycle par cycle.
 Cycle 0
Dans ce modèle nous étudions la déformation dans la structure de chaussée sous un

chargement. Le déplacement maximum est de 57,74.10-3m. Les résultats trouvés pour
la structure sont présentés par les figures suivantes :

Figure IV-20: Déformation de maillage Figure IV-19: Déplacement total
Figure IV-22:Déplacement horizontal total
Figure IV-21: Contrainte totale

 Cycle 1
Nous utilisons ici les caractéristiques du sol au cycle 1. Le déplacement maximum est de
55,51.10-3m. Les résultats trouvés pour la structure sont présentés par les figures
suivantes :

Figure IV-24 : Déplacement total

Figure IV-23: Déformation de maillage
Figure IV-26:Contrainte totale Figure IV-25: Déplacement horizontal total

 Cycle 2
Nous utilisons ici les caractéristiques du sol au cycle 2. Les résultats trouvés pour la
structure sont présentés par. Le déplacement maximum est 70,86.10-3m. Les résultats
trouvés pour la structure sont présentés par les figures suivantes :

Figure IV-28 : Déformation de maillage Figure IV-27 : Déplacement total
Figure IV-30 : Déplacement horizontal total

Figure IV-29: Contrainte totale

 Cycle 3
Nous utilisons ici les caractéristiques du sol au cycle 3. Les résultats trouvés pour la
structure sont présentés par. Le déplacement maximum est 88,37.10-3m. Les résultats
trouvés pour la structure sont présentés par les figures suivantes :

Figure IV-31 : Déplacement total
Figure IV-32: Déformation du maillage
Figure IV-33 : Déplacement horizontal total

Figure IV-34 : Contrainte totale

Le Tableau IV-7 présente le récapitulatif des résultats trouvés :
Tableau IV-7: Récapitulatif des sollicitations de la simulation de la structure de

chaussée
Cycle 0 1 2 3
Déplacement
57.74 55.51 70.86 88.37
maximal (10−3𝑚)
Déplacement
vertical 57.74 55.51 70.86 88.37
−3
(10 𝑚)
Déplacement
6.67 7.05 10.20 13.96
horizontal (10−3𝑚)
Contrainte totale
8.73 8.43 9.94 11
(103KN/m2)
Du Tableau IV-7 nous constatons que la valeur de la déformation a diminué entre le cycle
0 et 1 et a augmentée entre du cycle 1 à 3. La diminution constatée entre le cycle 0 et 1
serait peut-être dû à l’état initial du matériau avant les essais. Par contre l’augmentation
de la déformation remarquée du cycle 1 à 3 nous montre que la déformation augmente
lorsque le nombre de cycle augmente. Cette augmentation s’explique par la variation de
certains paramètres comme l’angle de frottement et le module œdométrique constatée
au niveau du cycle 2 et 3. L’augmentation des déformations des sols fins est la cause
majeure des détériorations des chaussées. Les dégradations constatées au niveau de la
chaussée Pobè-Kétou du PK 2+480 au PK 18+860, une année après sa construction seraient
dues à l’évolution de la déformation au fil des cycles.

CONCLUSION ET PERSPECTIVES

CONCLUSION ET PERSPECTIVES
La route Pobè-Kétou contribue au désenclavement des localités et à l’accroissement des
échanges entre le Bénin et le Nigeria. Elle améliore l’accessibilité de la zone du projet en
favorisant la libre circulation des personnes et des biens dans des conditions de sécurité
acceptables. Elle accroit également la production agricole, améliore la compétitivité du
secteur routier et stimule l’utilisation des ressources locales, humaines et matérielles. La
route construite repose sur un sol argileux, très fin, très plastique et dont le pouvoir de
gonflement est élevé. Ces sols en saison pluvieuse au contact de l’eau se gonflent mais en
saison sèche se rétractent. Ceci cause ainsi des dégradations précoces sur les
infrastructures routières construites, entravant le confort de la circulation et des pertes
dues aux accidents.
Afin d’étudier l’influence du cycle de retrait-gonflement du sol sur le comportement des

chaussée, une étude a été effectuée sur le tronçon de route Pobè – Kétou précisément du
PK 2+480 au PK 18+860 soit 16,380 km. Ce mémoire de fin de formation s’est organisé en
quatre chapitres.
 Dans un premier temps, la problématique a été posée et le point a été fait sur la
bibliographie. Dans cette partie également le comportement des chaussées sur sol
expansif et les pathologies observées sur les constructions ont été abordés
 Le deuxième chapitre fait la présentation du milieu d’étude, de la méthodologie, le
but des essais géotechniques réalisés sur ces échantillons dans le laboratoire
béninois CNERTP et la procédure suivie pour le cycle d’humidification et séchage
effectué sur les matériaux.
 Le troisième chapitre présente les résultats essais (essais d’identification et
essais mécaniques). Ces résultats illustrent qu’il s’agit de l’argile très plastique, non
organique et très compressible avec un gonflement élevé. De plus les résultats des
essais de cisaillement et œdométriques effectués sur les carottes ayant subi les
deux cycles d’humidification et séchage révèlent que :
 La cohésion du sol augmente lorsque le nombre de cycle passe de 1 à 2 et
diminue du cycle 2 à 3 tandis que L’angle de frottement diminue
 On note une croissance de l’indice de compressibilité pour le sol lors du
passage du cycle 1 au cycle 2 et une décroissance du cycle 2 au cycle 3. En
ce qui concerne l’indice de gonflement, on note une augmentation de
l'indice de gonflement du cycle 1 au cycle 2 et une diminution de cet indice
du cycle 2 au cycle 3.
 Le module œdométrique décroit du cycle 1 au cycle 3 à l’exception du sol à
la profondeur 0.5 à 1.50m qui décroit du cycle 1 à 2 et croit légèrement du
cycle 2 à 3.

 Enfin dans le quatrième chapitre, le logiciel géotechnique Plaxis a permis d’étudier

le comportement de la chaussée en question sous la variation des paramètres du
sol au fil des cycles. D’après les résultats obtenus suite au calcul de la déformation
pour chaque cycle, il révèle que la déformation de la chaussée augmente lorsque le
nombre de cycle augmente.
Aux termes de ces travaux, des perspectives s’ouvrent à nous. Afin de compléter cette
étude nous proposons :
 D’augmenter le nombre de cycle jusqu’à 20 afin d’étudier la variation des

paramètres mécaniques du sol expansif après avoir subi plusieurs cycles de retrait-
gonflement dus à l'alternance des saisons humides et sèches.
 De faire l’effet réciproque c’est à dire le séchage avant l’humidification et
comparer les résultats avec les résultats obtenus avec l’humidification avant
séchage.
 De déterminer la contrainte maximale de ses sols qui aurait neutralisé le
gonflement sans créer un tassement excessif dans le sol lors du retrait.
 D’étudier la microstructure des échantillons après les cycles.
Nous suggérons aussi pour des constructions futures, que les paramètres mécaniques
variant sous les phénomènes cycliques d'assèchement et d'humidification dus à
l'alternance des saisons humides et sèches, soient pris en compte dans les études
géotechniques pour mieux garantir la stabilité des ouvrages de génie civil dont les
paramètres visiblement décroît sous le joug des cycles hydriques.

RÉFÉRENCES BIBLIOGRAPHIQUES

RÉFÉRENCES BIBLIOGRAPHIQUES
[1] Crespin and P. YABI, “Construction de route sur sols mouvants ; quelle structure pour quel
trafic ? Application au tronçon Séhouè-Massi de la RNIE N°2.,” 2011.
[2] M. A. Djellali, “Modélisation des chaussées sur sols supports expansibles et estimation de
la pression de gonflement : Cas de la région de Tébessa ( N . E Algérie ),” pp. 2017–2018,
2018.
[3] Y. Sedjro, T. Kiki, Y. “Caractérisation minéralogique, thermique et microscopique des sols

fins en technique routière Yvette,” p. 211, 2017.
[4] T. Kiki and S. Yvette, “INFLUENCE OF SWELLING ON THE SETTLEMENT OF CLAY

SOILS OF,” 2009.
[5] R. D. E. Stage, “RAPPORT DE STAGE THEME : ETUDE DU PHENOMENE DE

GONFLEMENT LIE A LA PERMEABILITE DES SOLS ARGILEUX ”.
[6] M. S. M. Shah, “Rédacteur en chef,” Bull. Cancer, vol. 103, no. 6, pp. I–II, 2016, doi:
10.1016/s0007-4551(16)30133-3.
[7] Y. T. Kiki, N. Saiyouri, V. Gbaguidi, Y. Anguy, C. Gaborieau, and R. Fabre, “Microstructural,

Chemical and Mineralogical Analyses for Understanding the Geotechnical Properties of
Clayey Soils,” J. Miner. Mater. Charact. Eng., vol. 04, no. 06, pp. 305–319, 2016, doi:
10.4236/jmmce.2016.46027.
[8] Z. Virtus “ INFLUENCE DU GONFLEMENT SUR LE PEUPLEMENT DES SOLS ARGILEUX

DE LA COMMUNE HOUEYOGBE DANS LA DÉPRESSION DE LAMA“ 2009.
[9] H. Donald, “ETUDE EXPLORATOIRE DU COMPORTEMENT DES FONDATIONS

SUPERFICIELLES SUR SOLS GONFLANTS : CAS D’ISSABA Présenté,” 2020.
[10] Y. T. K. V. T. V. S. G. Bertrand GBAFFONOU, “INFLUENCE OF DRYING-WETTING

CYCLES ON THE COMPRESSIBILITY OF CLAY SOILS IN THE COMMUNE OF
HOUEYOGBE,” Academia.Edu, vol. 6, no. 5, pp. 125–131, 2017, [Online]. Available:
https://www.academia.edu/download/33659771/148.pdf
[11] M. E. HAZOUME, “Etude des performances des techniques de construction de route sur
sol gonflant au Bénin Présenté,” pp. 2020–2021.
[12] S. Géochimie and G. Millot, “Géologie des argiles. Alterations. Sédimentologie. Géochimie
par G. Millot, 1964,” pp. 165–166, 2018.
[13] B. Mohamed, “Étude expérimentale sur la stabilisation des sols gonflants de la région de
Mostaganem et de Relizane,” 2011.
[14] N. S. M. Vincent, J. Bouchut, J.-M. Fleureau, F. Masrouri, E. Oppenheim, J.-V. Heck, N.

Ruaux, S. Le Roy, I. Dubus, “Étude des mécanismes de déclenchement du phénomène de
retrait-gonflement des sols argileux et de ses interactions avec le bâti Rapport final
Étude des mécanismes de déclenchement du phénomène de retrait-gonflement des sols
argileux et de ses interactions,” 2006.
[15] M. CHRETIEN, “Compréhension des mécanismes de retrait-gonflement des sols argileux :

approche sur site expérimental et analyse de sinistres sur constructions individuelles,”
2010.

[16] P. Magnan, “Panorama des sols gonflants en géotechnique,” pp. 85–103, 2013.
[17] B. M. NOUHOLIDJI, “Etude du phénomène de gonflement des sols argileux dans le bassin
sédimentaire côtier au Bénin: Cas de la commune de Dogbo,” 2016,
[18] Bala Wenceslas SANOU, “Monographie de la commune de Bilanga,” 2006.
[19] AFNOR, “NF P 94-050 Détermination de la teneur en eau pondérale des matériaux.”
[20] AFNOR, “NF P 94-051 - Détermination des limites d’Atterberg.” p. 16, 1993.
[21] A. Djedid, A. Bekkouche, and S. M. Aissa Mamoune, “Identification et prévision du

gonflement de quelques sols de la région de Tlemcen (Algérie),” Bull. des Lab. des Ponts
Chaussees, no. 233, pp. 67–75, 2001.
[22] AFNOR, “NF P 94-054 Détermination de la masse volumique des grains solides des sols.”
[23] AFNOR, “NF P 94-055 Détermination de la teneur pondérale en matière organique d’un
sol.” 1993.
[24] AFNOR, “N F P 94-056 S o l s : r e c o n n a i s s a n c e et e s s a i s Analyse

granulométrique Méthode par tamisage à sec après lavage,” 1996.
[25] AFNOR, “NF P94-071 - Sols : reconnaissance et essais - Essai de cisaillement rectiligne à
la boîte - Partie 1 : cisaillement direct.” p. 16, 1994.
[26] “S o l s : r e c o n n a i s s a n c e et e s s a i s,” 1997.
[27] M. B. M. ZMALI., “ETATS D’EQUILIBRE LIMITE,” pp. 16–21.
[28] H. Ejjaaouani, “Interactions of foundations and expansive soils : pathology , calculations

and experimental studies Houssine Ejjaaouani To cite this version : HAL Id : pastel-
00005300,” 2009.
[29] N. O. T. HOUANKPO, “Etude d ’ une route revêtue sur sol compressible : Cas du tronçon
de route Sèhouè-Massi de la route,” 2010.

Annexes

Annexe A : Tableau récapitulatif des résultats des essais de laboratoire sur

Puits Manuels : site d’Issaba
Analyse Granulométrique par tamisage (%Passant au tamis mm) Analyse Granulométrique par sédimentation (%
N° Type de Profondeur D
général sondage (m) MAX
10 8 6.3 5 4 3.15 2 1 0.5 0.4 0.315 0.25 0.08 0.063 0.055 0.038 0.025 0.017 0.012
1 0,00/0,50 5 - - 100 99 99 99 98 97 97 96 96 92 92 92 88 85 82 82
-
2 0,50/1,00 5 - - 100 99 99 99 98 98 98 97 97 93 93 93 89 89 86 83
-
3 1,00/1,50 4 - - 100 99 99 99 99 99 99 98 95 94 94 91 91 88 85
-
PM
4 1,50/2,00 0.4 - - - - - - - - 100 99 99 96 96 96 92 92 89 89

-

Annexe B : Courbe contrainte-déformation de l’argile d’Issaba cycle 0
40 Profondeur 0-0.5
Contrainte de cisaillement
30
essai1
20
(kPa)
essai2
10 essai3
essai4
0
0 1 Déplacement
2 3 (mm)
4 5 6
Profondeur 0.5-1
30
25
20
essai1
(kPa)
15
essai2
10
essai3
5
essai4
0
0 1 Déplacement
2 3 (mm)
4 5 6
Profondeur 1.5-2
80
60
essai1
(kPa)
40
essai2
20 essai3
essai4
0
0 1 Déplacement
2 3 (mm)
4 5 6
Profondeur
50 1-1.5
40
30 essai1
(kPa)
20 essai2
essai3
10
essai4
0
0 1 Déplacement
2 3 (mm)
4 5 6

Annexe C : Courbe contrainte-déformation de l’argile d’Issaba cycle 1

Profondeur 0-0.5
200
150
essai1
100
(kPa)
essai2
50 essai3
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6
Profondeur 0.5-1
250
200
150 essai1
(kPa)
100 essai2
essai3
50
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6
250 Profondeur 1-1.5

200
150 essai1
(kPa)
100 essai2
essai3
50
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6
Profondeur 1.5-2
250
cisaillement (kPa)
200
Contrainte de
150 essai1
100 essai2
essai3
50
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6

Annexe D : Courbe contrainte-déformation de l’argile d’Issaba cycle 2

Profondeur 0-0.5
150
cisaillement (kPa)
Contrainte de
100 essai1
essai2
50
essai3
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6
Profondeur 0.5-1
150
100
essai1
(kPa)
essai2
50
essai3
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6
Profondeur 1-1.5
120
100
80
essai1
60
(kPa)
essai2
40
essai3
20
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6
Profondeur 1.5-2
150
100
essai1
(kPa)
essai2
50 essai3
essai4
0
0 1 2 3 (mm)
Déplacement 4 5 6

Annexe D : Courbe contrainte-indice de vide de l’argile d’Issaba cycle 1
10.0 100.0 1000.0 10000.0

Contraintes normales 1.0 10.0 100.0 1000.0 10000.0
0.44
0.80 Contraintes normales
0.42
0.70
Indices des vides
0.40
Indices des vides

0.60
0.38 0.50
0.36 0.40
0.34 0.30
0.32 0.20
Courbe oedométrique 0-0.5 Courbe oedométrique 0.5-1
Contraintes normales
1.0 10.0 100.0 1000.0 1.0 10.0 100.0 1000.0 10000.0
1.00 0.70 Contraintes normales
0.95
0.60
0.90
Indices des vides
Indices des vides
0.85 0.50
0.80
0.75 0.40
0.70
0.30
0.65
0.60 0.20
0.55
0.50 0.10
Courbe oedométrique 1.5-1 Courbe oedométrique 1.5-2

Annexe E : Courbe contrainte-indice de vide de l’argile d’Issaba cycle 2
1.0 10.0 100.0 1000.0 1.0 10.0 100.0 1000.0

1.20 Contraintes normales Contraintes normales
1.20
1.10
1.10
1.00
Indices des vides
Indices des vides

1.00
0.90
0.80 0.90
0.70 0.80
0.60 0.70
0.50
0.60
0.40
0.30 0.50
Courbe oedométrique 0.5-1
Courbe oedométrique 0-0.5
1.0 10.0 100.0 1000.0 1.0 10.0 100.0 1000.0

1.20 Contraintes normales 1.20 Contraintes normales
1.10 1.15
1.10
Indices des vides
Indices des vides
1.00
1.05
0.90 1.00
0.80 0.95
0.90
0.70
0.85
0.60 0.80
0.50 0.75
Courbe oedométrique 1-1.5 Courbe oedométrique 1.5-2

Annexe F : Photos
Puit pour prélèvement Prélèvement des échantillons intact
Prélèvement remanié Echantillon intact

Bac pour humidification Humidification
Carotte
Séchage

Appareil de cisaillement
Le bâti œdométrique

Table des Matières

Sommaire ............................................................................................................................................................... i
Dédicaces .............................................................................................................................................................. ii
Remerciements ................................................................................................................................................... iii
Liste des sigles et abréviations ....................................................................................................................... v
Liste des figures ............................................................................................................................................... vii
Liste des tableaux ........................................................................................................................................... viii
Liste des photos ................................................................................................................................................. ix
Résumé .................................................................................................................................................................. x
Abstract .............................................................................................................................................................. xi
Introduction ..........................................................................................................................................................1
Chapitre I : Généralités sur les argiles et le phénomène de retrait-gonflement. ............................... 7
I.1 Les argiles .................................................................................................................................................. 8
Introduction ................................................................................................................................................ 8
I.1.1 Définition ............................................................................................................................................ 8
I.1.2 Aspect microstructure des argiles ............................................................................................... 9
I.1.2.1 Minéraux argileux ...................................................................................................................... 9
I.1.2.2 Les caractéristiques minéralogiques des argiles ............................................................... 9
I.1.3 Les différentes familles de minéraux argileux ........................................................................ 11
I.1.3.1 Les Kaolinites ............................................................................................................................ 12
I.1.3.2 Les Montmorillonites .............................................................................................................. 13
I.1.3.3 L’illite ......................................................................................................................................... 14
I.1.4 Les phénomènes de gonflement et de retrait des argiles ..................................................... 15
I.1.4.1 Gonflement au sens physico-chimique ................................................................................. 15
I.1.4.2 Gonflement au sens mécanique ............................................................................................ 17
I.1.5 Le retrait de l’argile ....................................................................................................................... 18
I.1.6 Localisation des sols gonflants au Bénin .................................................................................... 18
I.2 Généralités sur les chaussées revêtues ...........................................................................................20
I.3 Comportement et pathologie des structures de chaussées sur sols gonflants ....................... 21
I.3.1 Une route sur sol gonflant ............................................................................................................ 21
I.3.2 Pathologie des structures de chaussées sur sols gonflants .................................................23
Chapitre II : Matériels, matériaux et méthodologie ................................................................................24
II.1 Milieu d’étude ........................................................................................................................................25

II.1.1 Situation géographique ................................................................................................................25
II.1.2 Population et structure administrative ...................................................................................25
II.1.3 Climat ..............................................................................................................................................25
II.1.4 Relief [9] ........................................................................................................................................26
II.1.5 Sols [9]............................................................................................................................................26
II.1.6 Présentation des sites de prélèvement ...................................................................................27
II.2 Méthodologie de recherche ..............................................................................................................29
II.2.1 Recherche et localisation de site de prélèvement ................................................................29
II.2.2 Prélèvements .................................................................................................................................29
II.2.3 Description des essais ................................................................................................................29
II.2.3.1 Essais d’identification .........................................................................................................29
II.2.3.1.1 Essai de détermination de la teneur en eau pondérale (NFP 94 -050) [19] .....29
II.2.3.1.2 Limites d’Atterberg (NF P 94-051)[20] ..................................................................30
II.2.3.1.3 Détermination du Poids spécifique des grains solides (NF P 94-054)[22] ......32
II.2.3.1.4 Essai de détermination de la teneur en matière organique e (XP P 94-055)[23]

.............................................................................................................................................................32
II.2.3.1.5 Analyse granulométrique par tamisage (NF P 94-056) [24] ...............................34
II.2.3.1.6 Analyse granulométrique par sédimentométrie (NF P 94-057)[24] ..................34
II.2.3.2 Classifications utilisées pour la caractérisation des matériaux ................................34
II.2.3.3 Essais mécaniques ................................................................................................................36
II.2.3.3.1 Essai de cisaillement direct (NF P 94-071-1) [25] ................................................36
II.2.3.3.2 Essai œdométrique (XP P 94-090-1)[26] ................................................................36
II.3 Description du mode opératoire de l’humidification et de séchage .........................................38
II.3.1 Mode opératoire de l’humidification ........................................................................................38
II.3.1.1 Cinétique de l’humidification ...............................................................................................40
II.3.2 Mode opératoire du séchage ..................................................................................................... 41
II.3.2.1 Cinétique du séchage ............................................................................................................ 41
Chapitre III : Caractérisation au laboratoire des sols naturels ...........................................................42
Introduction ..................................................................................................................................................43
III.1 Essais d’identification .......................................................................................................................43
III.1.1 Résultats des paramètres d’identification ............................................................................43
III.1.2 Analyse et interprétation des résultats des essais d’identification ..............................44
III.1.2.1 Teneur en eau in situ ..........................................................................................................44
III.1.2.2 Analyse granulométrique par tamisage et par sédimentométrie .............................44

III.1.2.3 Matière organique ..............................................................................................................46
III.1.2.4 Limite d’Atterberg .............................................................................................................46
III.2 Essai mécanique .................................................................................................................................48
III.2.1 Essai de cisaillement ..................................................................................................................48
III.2.2 Essai œdométrique ....................................................................................................................50
III.3 Estimation du gonflement par la méthode indirecte [21] ......................................................... 51
III.4 Présentation des résultats des essais sur les échantillons ayant subi le cycle
d’humidification et de séchage ..................................................................................................................53
III.4.1 Déroulement de l’essai ..............................................................................................................53
III.4.2 Texture et photos des échantillons en fonction du nombre de cycles ..........................53
III.4.3 ........................................................................................................................................................54
III.4.4 Description des textures des échantillons au cours des cycles .....................................55
III.4.5 Les résultats des essais de cisaillement ..............................................................................55
III.4.6 Les résultats des essais œdométriques ...............................................................................56
Conclusion .......................................................................................................................................................59
Chapitre IV : Modélisation et prédiction du comportement des sols expansifs sous léffet des
sollicitations cycliques : humidification et séchage. .................................................................................60
Introduction .................................................................................................................................................. 61
IV.1 Historique de la structure ................................................................................................................. 61
IV.2 Présentation du logiciel Plaxis ..........................................................................................................62
IV.2.1 Généralité sur Plaxis ...................................................................................................................62
IV.2.2 Fonctionnalités principales du logiciel sont : .........................................................................63
IV.2.3 Les points forts de Plaxis 2D ...................................................................................................63
IV.2.4 Les options par défaut et les solutions approchés ..............................................................63
IV.2.4.1 Le maillage .............................................................................................................................63
IV.2.4.2 Les conditions aux limites ..................................................................................................64
IV.2.4.3 Les pressions interstitielles ..............................................................................................64
IV.2.4.4 La conduite des calculs non linéaires ...............................................................................65
IV.2.4.5 Le coefficient de sécurité .................................................................................................65
IV.2.4.6 Les modèles de comportement utilisés dans Plaxis......................................................65
IV.2.4.7 Comportement élasto-plastique ........................................................................................66
IV.2.4.8 Modèle élastique linéaire ...................................................................................................67
IV.2.4.9 Modèle de Mohr-Coulomb ..................................................................................................68
IV.2.4.10 Modèle pour sols mous « Soft Soil Model (S.S.M) »..................................................70
IV.2.4.11 Modèle de sol avec écrouissage « Hardening Soil Model (H.S.M) »........................ 71

IV.2.4.12 Modèle défini par l’utilisateur.........................................................................................72
IV.2.5 Les paramètres fondamentaux de Plaxis ...............................................................................72
IV.2.5.1 Le module de Young ..............................................................................................................72
IV.2.5.2 Le coefficient de poisson ...................................................................................................73
IV.2.5.3 L’angle de frottement .........................................................................................................73
IV.2.5.4 La cohésion ............................................................................................................................73
IV.2.5.5 L’angle de dilatance (Ψ) ......................................................................................................74
IV.2.5.6 Les paramètres de compressibilité .................................................................................74
IV.2.5.7 ..................................................................................................................................................75
IV.2.5.8 La Prise en compte de la surconsolidation......................................................................75
IV.2.5.9 Les paramètres avancés .....................................................................................................75
IV.3 La modélisation numérique de la chaussée.....................................................................................76
IV.3.1 Géométrie du modèle...................................................................................................................76
IV.3.2 Caractéristique du sol ................................................................................................................76
IV.3.3 Caractéristiques des matériaux utilisés en couche de chaussée ......................................80
IV.3.3.1 Génération du maillage ........................................................................................................80
IV.3.3.2 Définition des conditions initiales ................................................................................... 81
IV.3.3.2.1 Conditions hydrauliques ............................................................................................... 81
IV.3.3.2.2 Pressions interstitielles initiales .............................................................................. 81
IV.3.3.2.3 Contrainte initiale ........................................................................................................82
IV.3.4 Procédure de calcul .....................................................................................................................84
IV.4 Les résultats de calcul .......................................................................................................................86
CONCLUSION ET PERSPECTIVES ..............................................................................................................95
RÉFÉRENCES BIBLIOGRAPHIQUES ..........................................................................................................98
Annexes ............................................................................................................................................................. 101
Annexe A : Tableau récapitulatif des résultats des essais de laboratoire sur Puits Manuels :
site d’Issaba ................................................................................................................................................ 102
Annexe B : Courbe contrainte-déformation de l’argile d’Issaba cycle 0 ....................................... 103
Annexe C : Courbe contrainte-déformation de l’argile d’Issaba cycle 1 ........................................ 104
Annexe D : Courbe contrainte-déformation de l’argile d’Issaba cycle 2 ....................................... 105
Annexe D : Courbe contrainte-indice de vide de l’argile d’Issaba cycle 1 ..................................... 106
Annexe E : Courbe contrainte-indice de vide de l’argile d’Issaba cycle 2 ..................................... 107
Annexe F : Photos ....................................................................................................................................... 108
Table des Matières.......................................................................................................................................... 111


Memoire 10

Transféré par

Droits d'auteur :

Formats disponibles

Memoire 10

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Memoire 10

Transféré par

Droits d'auteur :

Formats disponibles

RÉPUBLIQUE DU BÉNIN

MEMOIRE DE FIN DE FORMATION

POUR L’OBTENTION DU DIPLOME D’INGENIEUR DE CONCEPTION

Option : Ingénierie Routière

Influence des cycles de retrait-gonflement des

Dr. Ing. AHLINHAN Marx Ferdinand

Maître de conférence des Universités Soutenu publiquement le 16 décembre

Maître-Assistant des Universités

HOUEMAVO Donald Membres

Maître de conférence des Universités

Année Académique : 2021-2022

Liste des sigles et abréviations .................................................................................................................................................. v

Liste des figures ............................................................................................................................................................................ vii

Liste des tableaux ........................................................................................................................................................................viii

Liste des photos ............................................................................................................................................................................. ix

Chapitre I : Généralités sur les argiles et le phénomène de retrait-gonflement. ......................................................7

I.1 Les argiles ..............................................................................................................................................................................8

I.2 Généralités sur les chaussées revêtues .................................................................................................................... 20

Chapitre II : Matériels, matériaux et méthodologie ........................................................................................................ 24

II.1 Milieu d’étude ................................................................................................................................................................... 25

II.2 Matériels d’étude ........................................................................................................................................................... 27

II.3 Méthodologie de recherche ........................................................................................................................................ 29

II.4 Description du mode opératoire de l’humidification et de séchage ............................................................... 38

Chapitre III : Caractérisation au laboratoire des sols naturels ................................................................................... 42

III.1 Essais d’identification ................................................................................................................................................. 43

III.2 Essai mécanique ............................................................................................................................................................ 48

III.3 Estimation du gonflement par la méthode indirecte [21] .................................................................................51

IV.1 Historique de la structure ............................................................................................................................................61

IV.2 Présentation du logiciel Plaxis ................................................................................................................................... 62

IV.3 La modélisation numérique de la chaussée ............................................................................................................. 76

IV.4 Les résultats de calcul ................................................................................................................................................. 86

CONCLUSION ET PERSPECTIVES ........................................................................................................................................ 95

RÉFÉRENCES BIBLIOGRAPHIQUES .................................................................................................................................... 98

Annexes .......................................................................................................................................................................................... 101

Table des Matières ..................................................................................................................................................................... 111

LUDOVIC METOGNISSE AGASSOUSSI II

La réalisation de ce mémoire a été possible grâce au concours de plusieurs personnes à

J’adresse également toute ma gratitude à l’endroit de Mr HOUEMAVO Donald, Ingénieur

Le Directeur de l’ENSTP Dr ALAMOU Adéchina Éric, le Directeur Adjoint de

Pr Edmond ADJOVI, Professeur Titulaire des Universités du CAMES ;

LUDOVIC METOGNISSE AGASSOUSSI III

Pr Aristide HOUNGAN, Professeur Titulaire des Universités du CAMES ;

LUDOVIC METOGNISSE AGASSOUSSI IV

Liste des sigles et abréviations

HRB : Highway Research Board

GTR : Guide de terrassement routier

MO : Matière Organique γs : Poids Spécifique

γsat : Poids volumique saturé

γd : Poids volumique sec CU : Cohésion non drainé

φU : Angle de flottement non drainé

eo : Indice des vides du sol en place

ei : Indice des vides initial

Cv : Coefficient de consolidation verticale

ơ'vo : Contrainte effective verticale

ơ'p : Contrainte de préconsolidation

𝜈𝑢𝑟 : Coefficient de poisson en charge-décharge (par défaut 0,15)

M : pente de la ligne de l'état critique

LUDOVIC METOGNISSE AGASSOUSSI V

LUDOVIC METOGNISSE AGASSOUSSI VI

Liste des figures