COURS ET EXERCICES

DE MECANIQUE DES SOLS

KHALED MEFTAH
Maître assistant

Première version : Septembre 2008

Cours Mecanique des sol 0 Khaled MEFTAH

 SOMMAIRE
Chapitres Intitulés Pages

Les sols : Structure – Identification

1 4
et classification
2 Hydraulique des sols 20
3 Contraintes dans les sols 28

4 Tassement et consolidation des sols 45

5 Résistance au cisaillement des sols 56

6 Etude de la portance des fondations superficielles 70

à partir d’essais de laboratoire

7 Poussée et butée des terres 81

Références bibliographiques

Cours Mecanique des sol 1 Khaled MEFTAH

 TABLE DES MATIERES

Chapitre 1 : Les Sols : Structure. Identification et classification 4

1- Eléments constitutifs d’un sol 4
2- Paramètres de définition d’un sol 4
3- Identification des sols 7
4- Classification des sols 14
Exercices d’application 17

Chapitre 2: Hydraulique des sols 20

1- Ecoulement linéaire 20
2- Ecoulements plans 20
Exercices d’application 24

Chapitre 3 : Les contraintes dans le sol 28

1- Notions de contraintes 28
2- Cercle de Mohr 29
3- Les contraintes dues au poids propre des sols 32
4- Les contraintes dues aux surcharges 34
Exercices d’application 38

Chapitre 4 : Tassement et consolidation des sols 45

1- Notions de déformations 45
2- Relation contraintes - deformations 45
3- Tassements des sols- consolidation 46
3- Evolution du tassement en cours du temps 49
Exercices d’application 53

Cours Mecanique des sol 2 Khaled MEFTAH

Chapitre 5 : Résistance au cisaillement des sols 56
1- Comportement élasto-plastique des sols 56
2- Comportement à cours et à long terme des sols 57
3- Détermination des paramètres de cisaillement des sols 58
au laboratoire
Exercices d’application 65

Chapitre 6 : Portance des fondations superficielles 70

1- Introduction 70
2- Calculs de la capacité portante 72
Exercices d’application 77

Chapitre 7 : Poussée et butée des terres 81

1- Introduction 81
2- La théorie de Rankine 81
3- Calcul des forces de poussée et de butée 86
5- Stabilité des murs de soutènement 90
6- Les rideaux de palplanches 92
Exercices d’application 95
Références bibliographiques 97

Cours Mecanique des sol 3 Khaled MEFTAH

 CHAPITRE 1
LES SOLS : STRUCTURE – IDENTIFICATION
ET CLASSIFICATION

1- Eléments constitutifs d’un sol

Un sol est un mélange :
d’éléments solides : Provienant de la désagrégation mécanique et/ou
chimique d’une roche mère.
On distingue les minéraux non argileux (Æ>2mm et ayant le même
comportement que la roche mère : Sols pulvérulents), les minéraux
argileux ( kaolinite, illite et montmorillonite) et le sols organiques (vases et
tourbes)
d’eau : Existe sous plusieurs formes (eau de constitution, interfeuillets, liée
et libre).

de gaz : Contenu dans les vides,c’est l’air pour un sol sec ou mélange d’air
et de vapeur d’eau pour un sol humide.

2- Paramètres de définition des sols

2-1 Modèle élémentaire d’un sol

Un sol étant composé de grains solides, d’eau et d’air , on peut rassembler chaque
phase en un volume partiel unique de section unit. Les notations suivantes sont
utilisées : Volumes Poids

Va Air Wa=0

Vw Eau Ww

Vs Grains solides Ws

Cours Mecanique des sol 4 Khaled MEFTAH

Va : volume de l’air. Vw : volume de l’eau.
Vs : volume des grains solides. Vv = Va + Vw : volume des vides.
V = Vv + Vs : volume total du sol
W w: poids de l’eau Ws : poids des grains solides
W = Ww + Ws : poids total

2-2 Les poids volumiques

Le poids volumique (spécifique) total ou humide :

W
g = V
Le poids volumique des grains solides :
g =
W
s
s
Vs
Le poids volumique du sol sec :
Ws
gd=
V

Le poids volumique de l’eau :

Ww 3 3
gw = = 10kN / m = 1t / m
Vw

Poids volumique du sol saturé

W
g sat = V

Poids volumique déjaugé

'
g =gsat-gw
La gravité spécifique :

Gs =gs
g w

Cours Mecanique des sol 5 Khaled MEFTAH

2.3 Les paramètres d’état

Ils indiquent dans quelles proportions existent les différentes phases d’un sol.
On définit :
La teneur en eau :
Ww
w%= x100
Ws

L’indice des vides :

Vv
e=
Vs

Le degré de saturation :
Vw
Sr % = x 100
V
v
La porosité :
Vv
h = V

Tous ces paramètres ne sont pas indépendants. Ils sont reliés par des relations que l’on
peut retrouver à l’aide du modéle élémetaire. Exemple de formules :

g d = 1+gw

g d =1g+se

gs
e =g -
1d

Sr = wGse

Cours Mecanique des sol 6 Khaled MEFTAH

3. Identification des sols

Pour caractériser un sol, il faut déterminer les paramètres de nature et les paramètres
d’état.

Les paramètres de nature indiquent les caracteristiques intrinsèques du sol. IIs ne

varient pas au cours du temps (poids volumique des grains solides, granularité,
argilosité, limites d’Atterberg, teneur en matières organiques,…).

Les paramètres d’état sont fonction de l’état du sol et caractérisent le comporetement du

sol sous l’effet d’un chargement donné (teneur en eau, indice des vides, porosité,
Equivalent de sable,...).

Nous regroupons dans ce paragraphe les essais géotechniques de laboratoire

classiques qui permettent de caractériser un sol.

3.1 La masse volumique des particules solides gs

Sa détermination se fait à l’aide d’un pycnomètre. Une masse de sol sec ms est
introduite dans un pycnomètre conteneant de l’eau distillée. Aprés avoir éliminé toutes

les bulles d’air, on mesure le volume d’eau déplacé par les grains solides vs.
3 3
N.B : Pour les sols (à part les sols organiques) : 26 kN/m £ gS £ 28 kN/m

3.2 Les essais granulométriques

Ils permettent d’obtenir la répartition en pourcentage des grains solides selon leurs
dimensions. Deux types d’essais sont envisageables selon le sol à tester :
- Par tamisage (par voie humide ou sèche) pour les élements de diamétre
Æ ³ 80mm.
- Par sédimentométrie pour les élements de diamétreÆ < 80mm.

Cours Mecanique des sol 7 Khaled MEFTAH

Les résultats sont traduits sous forme d’une courbe granulométrique, tracee dans des
axes semi-logarithmiques, à partir de laquelle on peut déterminer :

- Le coéfficient d’uniformité de Hazen :

d60
=
Cu

- Le coéfficient de courbure :
2
d
Cc= 30
d xd
10 60
N.B : di : diamètre correspondant à i% de pourcentage de tamisat cumulé.

Fig 1.1 :Exepmle de détermination des di :

- d10 =0.17
- d30 = 0.58
- d60= 1.80

Cours Mecanique des sol 8 Khaled MEFTAH

3.3 Essais sur sols pulvérulents

Le comportement de ces sols dépend des paramètres qui caractérisent le squelette

solide, à savoir les dimensions des grains et l’indice des vides. Les essais les plus
courants sont :
a) Equivalent de sable (ES%) : Permet de caractériser la propreté des
sables et le type de sol analysé.
.
Tableau 1.1 : Caractérisation des sols à partir de la valeur de E.S
ES Type de sol

0 Argile pure
20 Sol plastique
40 Sol non plastique
100 Sable pur et propre

b) Densité relative (ou indice de densité ): Permet de caractériser la

compacité d’un sol grénu et son aptitude à supporter des charges.

ID emax -e (1.1)

= emax -emin
Avec :
- e : indice des vides du sol en place.
- emax : indice des vides du sol à l’état le plus lâche.
emin : indice des vides du sol à l’état le plus dense.

Tableau 1.2 : Compacité d’un sol en fonction de l’indice de densité

I Compacité du sol
D
0 Tres lâche
<0.5 lâche
0.5 moyennement dense
>0.5 très compact
1 très bien compact

Cours Mecanique des sol 9 Khaled MEFTAH

 3.4 Essais sur les sols fins
Le comportement de ces sols dépend :
- de la composition minéralogique (types de minéreaux argileux)
Les argiles sont composées d’alumino-silicates hydratés. Les grains solides ont une
forme de plaquette. Ils sont formés par un empilement de feuillets (composés d’une
superposition de couches octaédriques et tétraédriques constituées par un maillage
d’ions Si, O, OH, Al et Mg) :

- les feuillets 1/1 sont formés d’une couche tétraédrique et d’une couche
octaédrique (kaolinite),

- les feuillets 2/1 sont formés d’une couche octaédrique entourée de deux
couches tétraédriques (illite, smectite).

Si
Al Al 0,96 nm
Si 0,72 nm
Si
Feuillet 1/1 Feuillet 2/1

Le tableau 1.3 présente les familles d’argile et leurs caractéristiques.

Tableau 1.3 Familles d’argile et leurs caractéristiques

Noms Feuillets Nombre de Dimension d’une Surface

élémentaires feuillets par particule l x e spécifique
particule (µm x µm) (m²/g)

Kaolinite 100 – 150 1 x 0,1 20-70

Illite 10 0,3 x 0,01 10-40

Montmorillonite 1 0,1 x 0,01 100

Cours Mecanique des sol 10 Khaled MEFTAH

 - de la structure
dispersée : contact face-face entre les particules
floculée : contact bord-face entre les particules

- de la teneur en eau w%.

Elle est obtenue par passage a l’étuve a 105°C d’un e quantité de sol. C’est le rapport de la
masse d’eau évaporée a la masse du sol sec (grains solides).

a) Les limites d’Atterberg

Suivant la consistance d’un sol remanié, qui est fonction de sa teneur en eau, on distingue 4
états schématisés comme suit :

Solide plastique liquide

ETAT
Sans retrait Avec retrait
w ws wl wp
Limite de
retrait plasticité liquidité

ws, wl et wp sont les limites d’Atterberg déterminées en laboratoire sur la fraction du sol
passant au tamis 0.40mm (méthode de la coupelle de Casagrande et du rouleau et appareil
de retrait).
A partir de ces limites, on peut déterminer :

L’indice de plasticité « Ip »
(1.2)
I P =wL -wP

Tableau 1.4 : Type de sol en fonction de Ip

Indice de plasticité (%) Type de sol

< 1% Pulvérulent
1%<Ip<7% Sable argileux
7%<Ip<17% Argile sableuse
Ip>17% Argile

Cours Mecanique des sol 11 Khaled MEFTAH

 L’indice de consistance: « Ic »
Ic
= wL -w
(1.3)
IP
Tableau 1.5 : Etat de consistance du sol en fonction de Ic

Indice de consistance consistance du sol

Ic £ 0 Liquide
0 < Ic < 1 Plastique
Ic = 1 Solide plastique
Ic >1 Solide ou semi solide.

L’indice de liquidité « IL »

w-wP
IL= (1.4)
IP
Tableau 1.6 : Etat de consistance du sol en fonction de « I L »
Indice de liquidité Consistance du sol

IL< 0 Très dure

0<IL£1 Dure à très plastique

IL>1 Fluide

b) Valeur du bleu de méthylène : « VBS »

Elle représente la quantité de bleu pouvant être adsorbée sur les surfaces internes et
externes des particules du sol. La valeur VBS s’exprime en masse de bleu pour 100g de sol.

Cours Mecanique des sol 12 Khaled MEFTAH

 Tableau 1.7 : Type de sol en fonction de la valeur «VBS »
VBS Type de sol
0,1 Sol insensible à l’eau
0,2 Seuil au-dessus duquel le sol est sensible
à l’eau.

1,5 Seuil distinguant les sols sablo- limoneux

des sols sablo- argileux.

2,5 Seuil distinguant les sols limoneux peu

plastiques de sols limoneux moyennement
plastiques

6 Seuil distinguant les sols limoneux et les

sols argileux

8 Seuil distinguant les sols argileux des sols

très argileux.

b) Teneur en carbonate : % de CaCo3

L’essai est réalisé au calcimètre Dietrich-Fruhlingafin de déterminer la teneur pondérale en

carbonates d’un sol qui est le rapport entre la masse de carbonate contenue dans le sol à sa
masse sèche totale. La détermination se fait par décomposition du carbonate de calcium

CaCo3 contenu dans le sol par l’acide chlorhydrique.

Tableau 1.8 : Type de sol en fonction du % en CaCo3

Teneur en Carbonate en % Type de sol
0-10 Non marneux
10 - 30 Faiblement marneux
30 - 70 Marneux
70 - 90 Calco - marneux
90 - 100 Calcaireux – crayeux

Cours Mecanique des sol 13 Khaled MEFTAH

 c) La teneur en matières organiques : « MO »

C’est le quotient de la masse de matières organiques contenues dans un échantillon de sol

par la masse totale des particules solides minèrales et organiques. Sa détermination se fait
parcalcination.

Tableau 1.9 : Type de sol en fonction du % en MO

Teneur en matières organiques (MO%) Type de sol
MO< 3 Non organique
3 < MO < 10 Faiblement organique
3 < MO < 30 Moyennement organique
MO > o Très organique

4 - Classification des sols

Elle consiste à regrouper les sols qui ont une nature, un état et un comportement similaires
par rapport à une application géotechnique particulière (routes, fondations, etc..)
En première approximation, on peut adopter, lorsque les dimensions des grains sont peu
différentes, la classification suivante selon le diamètre moyen des grains

SOLS GRENUS SOLS FINS

cailloux grviers gros sable sable fin limon argile
20mm Ægrains
20mm 2mm 0.2mm 2mm

4.1 Classification des sols non organiques ( MO < 3%)

On distingue :
- Les sols grénus (plus de 50% des éléments solides ont unÆ>80mm)
- Les sols fins (plus de 50% des éléments solides ont un Æ>80mm).
a) Les sols grénus
On adopte la classification des laboratoires des ponts et chaussées (LPC).
(tableau 1.10 et figure 1.10 : diagramme de plasticite pour les éléments fins ).
b) Les sols fins :
Utiliser le diagramme de Casagrande (Figure 1.2).
Cours Mecanique des sol 14 Khaled MEFTAH
4.2 Classification des sols organiques (MO>3%)
Les caractéristiques utilisées pour la classification de ces sont :
- La teneur en matières organiques : % MO
· sol « fo » pour 3% < MO < 10%
· Sol « mo » pour 10% < MO < 30%
· Sol « to » pour MO > 30%

- Les limites d’Atterberg pour le sols « fo » (utiliser le diagramme de

Casagrande en rajoutant le terme « fo »).
- Les résultats du test d’humification Von Post pour les sols « mo » et « to ».
On obtient 10 classes de sols organiques de H 1 à H10 .

Cours Mecanique des sol 15 Khaled MEFTAH

 Tableau 1.10 Classification des sols grenus
(plus de 50% des éléments >
80µm)
Définitions Symboles Conditions Appellations
L.P.C.
CU = D 60 <4
D
Moins de 5% Gb 10 Grave propre bien
et
graduée
Plus de 50% d’éléments < 1<CC= (D30)²<3
D .D
Graves des éléments > 80µm 10 60
80µm ont un
diamètre > 2mm Gm Une des conditions de Gb Grave propre mal
non satisfaite graduée
Plus de 12% GL Limites d’Atterberg au-dessous de A Grave limoneuse
d’éléments <
80µm GA Limites d’Atterberg au-dessus de A Grave argileuse

CU = D 60 >6
D
Moins de 5% Sb 10 Sable propre bien
Et
gradué
Plus de 50% d’éléments <
Sables des éléments > 80µm 1<CC= (D30)²<3
D .D
80µm ont un 10 60

diamètre > 2mm

Sm Une des conditions de Sb Sable propre mal
non satisfaite gradué
Plus de 12% SL Limites d’Atterberg au-dessous de A Sable limoneux
d’éléments <
80µm

Si 5% d’éléments < 80µm < 12%, on utilise un double symbole

Figure 1.2. Classification des sols fins

Diagramme de plasticité

Cours Mecanique des sol 16 Khaled MEFTAH

 QUESTIONS A DEBATTRE
1) Quelle est la différence entre la masse volumique et le poids volumique d’un
sol? Citer la relation qui unit ces deux caractéristiques ?
2) Etant donné qu’il n’ y a pas d’essai qui mesure le degré de saturation d’un sol,
de quelle façon peut-on le quantifier ?
3) Sur quels types de sols les essais de limites d’Atterberg sont effectués ?
4) En plus de l’appareil de Casagrande, on parle aussi du pénétromètre à cône : à
quoi ca sert?
5) Peut-on réaliser l’essai de bleu sur un sable ?
Exercice 1
Montrer les égalités suivantes :
g
gd= g = s
1) 1+w 1+e

1 + w
2) g = 1 + egs

wGs
Sr =
3) e

Exercice 2
Des essais réalisés sur un échantillon de sol remanié ayant une teneur en eau à
l’état naturel de 21.5%, ont donné les résultats suivants :
- Analyse granulométrique( par voie humide et sédimentométrie)
Tamis(mm) 2,5 1,25 0,63 0,315 0,160 0 ,080 0,050 0,020 0,005 0,002
T(%) 100 99.90 99.80 99.30 98.90 98.60 85.30 65.30 43.50 31.00
- Limites d’Atterberg :

Limite de liquidité = 31.00 % et Limité de plasticité = 24.80 %.

1) Tracer la courbe granulométrique de ce sol en utilisant la fiche jointe en

annexe.
Calculer les coéfficients d’uniformité et de courbure. Commenter.
2) Déterminer les indices de plasticité, de liquidité et consistance. Commenter
3) Classer ce sol d’après la classification LPC.

Cours Mecanique des sol 17 Khaled MEFTAH

Exercice3
Les échantillons provenant d’un sondage carotté dans la région de la Soukra, nous
ont fournis les résultats suivants :
g = 19.1kN/m3 ; w = 33.56% ; gs=26.8KN/m3 ; wL= 42.2% wp= 18.3%
1) Déterminer le poids volumique sec, l’indice des vides et le degré de saturation
du sol.
2) Calculer les indices : de plasticité, de liquidité et de consistance. En déduire
l’état de consistance du sol.
3) Classer ce sol d’après la classification LPC.
Exercice4
On a effectué sur 4 échantillons de soms différents les essais d ’identification
dont les résultats sont :

Sol Tamisat Tamisat D10 D30(mm) D60(mm) Wl(%) WP(%)

2mm (%) 0,08mm(%) (mm)

1 93 14 0,06 0,16 0,35 45 65

2 70 50 - - - 38 25
3 56 3 0,2 0,75 2,4 - -
4 100 90 - - - 32 12

Appliquer à ces sols la classifications LPC.

Exercice5
Un échantillon de sol a un indice des vides égal à 0,6 et une teneur en eau de 15%.
Sachant que la gravité spécifique vaut 2,7 ; déterminer :
1) Le poids volumique sec
2) Le poids volumique total
3) La teneur en eau et son poids volumique à l’état saturé.
Exercice 6
Un échantillon d’argile saturée pesait 35.4g à l’état naturel et 24.2g après séchage à
3
l ‘étuve. Si le poid volumique des grains solides vaut 26,2 kN/m ; déterminer la
teneur en eau, l’indice des vides, la porosité, le poids volumique total, le poids
volumique sec et le poids volumique déjaugé.

Cours Mecanique des sol 18 Khaled MEFTAH

 ELEMENTS DE CORRECTION

Exercice3
gd = g = 19.1 =14.30KN /m3
1) 1+w 1.3356

g d =1g+se
gs
e =g
-1=0.87
d
wGs
Sr
e =0.3356x2.368
0.87
=1.
= 0

2) IP =wL -wP = 42.2-18.3=23.9

Ic =0.36 I = =0.6
= w L -w w -wP 4
L
IP IP

3) Il s’agit d’un sol argileux de consistance plastique

4) d’après la classification LPC : le sol est une argile peu plastique

Cours Mecanique des sol 19 Khaled MEFTAH
 CHAPITRE 2

HYDRAULIQUE DES SOLS

1- Ecoulement linéaire
1-1 Hypothèses
- Le sol sujet d’un écoulement est supposé saturé (Sr=1)
- Le régime d’écoulement est permanent et laminaire.
1-2 Mouvement de l’eau
Une molécule suit un trajet appelé « ligne de courant », son vecteur vitesse est tangeant
à cette ligne. Les lignes de courant s’appuyant sur le contour fermé d’une surface « S »
3
forment un tube de courant. Le débit « Q »en m /s, pour une vitesse « V » constante
est :
Q =VxS (2.1)

Pour effectuer des calculs,on est ramené à définir des lignes de courants fictives et des
vitesses apparentes « v »
1-3 Charge et pression hydraulique
Par sa position, la pression qu’elle subit et la vitesse à la quelle elle s’écoule, l’eau en un
point donné du sol porte une quantité d’énergie « h » en mètres d’eau (charge
hydraulique), donnée par l’équation de Bernoulli :
h v² + u +z
= 2g g (2.2)
w
v : Vitesse de l’eau.
g : Accélération de la pesanteur.
u : Pression de l’eau
z : Cote du point considéré par rapport à une surface de référence, peut être
négatif ou positif
Pour les sols, « v » est très faible, on aura alors :
u
h= +z
gw (2.3)

Cours Mecanique des sol 20 Khaled MEFTAH

La charge hydraulique est mesurée en un point donnépar l’altitude du niveau atteint par
l’eau dans un tube piézométrique placé au point considéré par rapport au plan de
référence.

1-4 Perte de Charge

Entre deux points A et B, Dh représente la variation de la charge hydraulique subie par

l’eau lors de son mouvement de A vers B. C’est une perte d ’énergie (perte de charge).

Dh=h A -h B
(2.4)

1-5 Gradient hydraulique

C’est la perte de charge par unité de longueur en un point donné.

i Dh (2.5)
= dl

Le gardient hydraulique critique (ic), est celui qui va provoquer un état de boulance
appelé phénomène de renard.
= G s - 1 = g' (2.6)
ic = Dhc
L 1 + g w
e

1-6 Loi de Darcy

Pour un sol donné, la vitesse « v » reste proportionnelle au gradient hydraulique « i »

selon la loi de DARCY :
v=k.i (2.7)

k : étant le coefficient de perméabilité du sol qui varie en fonction de la nature du sol et

qui peut être déterminé soit à partir des essais de laboratoire ou à partir d’essais en
place.

Cours Mecanique des sol 21 Khaled MEFTAH

2- Ecoulement plan
Pour résoudre un problème d’écoulement plan dans un sol saturé, il faut connaître en tout
point du sol la charge hydraulique. En se basant sur le principe de continuité du débit et
en supposant le sol homogène et isotrope vis-à-vis de la perméabilitéK, on obtient
l’équation de conservation du débit :

¶ ²h + ¶ ²h = 0
¶x² ¶z² (2.8)

Qui peut s’écrire sous la forme Dh=0 : Equation de Laplace

Cette équation admet une solution lorsque les conditions limites et initiales sont définies
pour l’écoulement. L’integration de cette équation nous donne deux familles de courbes
orthogonales. Par construction de ces courbes, on obtient un réseau d’écoulement
orthogonal constitué de lignes équipotentielles j (même charge hydraulique sur une
même ligne) et des lignes de courant y (tangeantes au gradient hydraulique). La
connaissance de ce réseau nous fournit en tout point la vitesse de l’eau « v », la charge
hydraulique « h », la pression interstitielle « u », et le débit « q ». La résolution de
l’équation (2.8) peut se faire soit par la méthode graphique, soit par la méthode
analytique par traitement numérique ou bien par la méthode par analogie électrique.

Résolution graphique : on se propose d’étudier l’exemple suivant

z
Coté amont Coté aval

A B C D P.Réf
O
I 3 4
2 VIII
ligne de courant
II VII
1 E ligne
III VI
IV V
équipotentielle
M
F G

Substratum imperméable

Fig 2.1 - Réseau d’écoulement horizontal

Cours Mecanique des sol 22 Khaled MEFTAH

 Conditions aux limites
- BEC: ligne d’écoulement.
- FG : ligne de courant
- AB : ligne equipotentielle hA=hB=h
- CD : ligne équipotentielle hC=hD=0
Pour tracer le réseau d’écoulement, certaines conditions doivent être satisfaites :
- lignes de courant orthogonales aux lignes équipotentielles.
- les quadrilataires curvilignes doivent être semblables.
- les conditions aux limites satisfaites.
- même dédit et même débit et même perte de charge entre deux lignes de
courant.
Calcul du débit
Le débit traversant un quadrilatère est donné par :
a
Dq=K.Dh.
b (2.9)

ligne
équipotentielle
b ligne de courant

C’est le débit traversant un canal « i »

q = Dq = K. h. a
b (2.10)
i
Dh étant la perte de charge élémentaire.
Si on appelle :
nh : nombre d’intervalles entre les lignes équipotentielles
nc : nombre de tubes d’écoulement (de canaux)

Cours Mecanique des sol 23 Khaled MEFTAH

 on aura le débit total :
q n q n K. h. b
= c i= c D

si la perte de charges totale entre la 1ére et la dernière ligne équipotentielle est :

DH= nh . Dh
le débit total de fuites du coté amont vers le coté aval est donné par la relation :
n
Q K. . c . a
b
n
= h (2.11)

Dans le cas d’un réseau à mailles carrées (a/b =1)

Dans le cas de l’exemple de la figure 2.1, on a :

nh = 8 ; nc= 4 et DH= hA – hD= h ;

Q 1
Kh.
= 2
le débit total de fuite est :
Calcul des charges hydrauliques et des pressions :

Pour le point « M » représenté sur l’exemple de la fig 2.1

h 2
hM =hA-6Dh
8 = 8h
; hM =hA-6
u =g .h sachant que h : hauteur piézométrique
et M w wM wM

uM =gw. hM -z ( )
M zM : mezurée à partir du plan de référence (zM <0)
2
uM =g w . h -zM

8
Calcul de la force d’écoulement :

La force de l’écoulement est égale à γw i. Elle est tangente a la ligne de courant.

Cours Mecanique des sol 24 Khaled MEFTAH

 QUESTIONS A DEBATTRE

1- Pourquoi néglige-t-on la charge de vitesse en géotechnique ?

2- A quoi sert le gradient hydraulique critique ? Quels sont les facteurs qui
l’influencent ?
3- Dans un écoulement plan, où la vitesse de l’eau est-elle la plus rapide ?
Où le phénomène peut-t-il se produire ?
4- Coment peut-on diminier le risque dapparition du phénomène de renard ?

EXERCICE1
Calculer le gradient hydraulique critique ‘un sable dont la porosité est de 40%
et dont la gravité spécifique est de 2.12.

EXERCICE2
Ondoit effectuer une excavation dans un dépôt d’argile imperméable ayant
une épaisseur de 10m et reposant sur une couche de sable compact.

Sachant que la nappe phréatique se situe à 3m sous la surface du sol,

déterminer la profondeur maximale de l’excavation juste avant l’apparition du
renard dans la couche d’argile.
La masse volumique de l’argile saturée est de 18.20 kN/m3, tandis que celle
du sable est de 21.5 kN/m3.

Cours Mecanique des sol 25 Khaled MEFTAH

 EXERCICE3
Un barrage doit être fondé sur une couche d’alluvions peméables limité à 20
m de profondeur par un substratum horizontal imperméable. La largeur de ce
barrage est de 25m.
La différence du niveau d’eau entre l’amont et l’aval est de 7.50m.
Le réseau d’écoulement à mailles carrées est tracé sur la figure suivante :

Coté aval
Coté amont
A P.R

V E
IV C 15
1 III D
II
14
I
2 13
3
4 56 7 8 9 10 11 12

Substratum imperméable

1) Calculer la pression interstitielle au point C du contact barrage alluvions

situé à mi-distance du parement amont et du pied aval du barrage.

2) Evaluer le gradient hydraulique de sortie au contact du pied aval du barrage

entre les points D et E (DE=2m).
En déduire le coéfficient de sécurité vis-à-vis du phénomène de renard
i
= cr
Fs ; ( icr : gradient hydraulique critique)
i

3) Calculer le débit traversant le sol

Cours Mecanique des sol 26 Khaled MEFTAH

 ELEMENTS DE CORRECTION

Exercice1
Gs - 1
= =
ic 1 . 07

1+e

Exercice2
Profondeur maximale d’éxcavation = 6.15m

Exercice4

1)Pression interstitielle au point C :

On sait que, pour le point C , la charge hydraulique est :
u c
hc =g w +zc

zc = -2m ;
d’après le réseau d’écoulement

h c = h A - 5 . 75 D h

h =7.5m Dh= DH n 15
Avec A et On ; D h = 7 . 5 m et h=
nh
h = 4 . 62 m
aura : c , et
u c = ( 4 . 62 - ( - 2 )). 10 = 66 . 2 kN / m 2
2)Gradient hydraulique de sortie
i Dh = DH =0.25
=
DE
nhDE
i = g' Fs = g' =4.4
cr D’ou.
g 0 . 25 g w
w
La sécurité au phénomène de renard est assurée.
3) Débit traversant le sol
Q=k.Dh.nC et -5
Avec, n h = 15n c = 5 k = 4 . 10 m/s
nh

-4
Q = 10 m3/s

Cours Mecanique des sol 27 Khaled MEFTAH

 CHAPITRE 3

CONTRAINTES DANS LES SOLS

1- Notions de contraintes
Soit un solide à la surface duquel s’exercent des forces.

n d
F
II ds
M S
P
I I

s normale

dF
tangeante

0 t
fig 3.1
En coupant ce solide par un plan fictif (P), l’élément de surface « ds », autour
du point « M » sur la surface « S », est soumis à une force R
d F(fig 3.1).
R R
dF
f =
la contrainte au point « M »est le vecteur ds
Cette contrainte se décompose en une contrainte normale s et une contrainte

tangentielle t
En mécanique des sols, pour déterminer l’état de contraintes autour d’un point
« M » dans le sol, il suffit de connaître les composantes des forces s’éxerçant
sur les faces d’un paraléllipipède centré autour du point « M » et dont les arêtes
sont parallèles aux axes Ox, Oy, Oz.
L’état de contraintes au point M est défini par une matrice symétrique appelée
tenseur de contraintes :

Cours Mecanique des sol 28 Khaled MEFTAH

 s xt yxt zx
s: t xys yt zy
t xzt yz s z

Parmi les facettes autour du point M, il existe 3 plans priviligés pour les quels
la contrainte tangentielle est nulle (t = 0). Ces 3 plans sont appelés plans
principaux,
Leurs directions normales, directions principales et les contraintes
correspondantes, contraintes principales, notées
s1 : Contrainte principale majeure.
s2 : Contrainte principale intermédiaire.
s3 : Contrainte principale mineure.
Avec : s1 ³s2 ³ s3

2- Cercle de Mohr
Pour étudier l’état de contraintes autour d’un point, on utilise une
représentation appelée diagramme de Mohr qui consiste à représenter le
R
f
vecteur contrainte dans un système d’axes (s,t).
Dans le cas bidimentionnel, cas très fréquent en géotechnique, le cercle de
Mohr est le lieu des extrémités des vecteurs contraintes et les contraintes
principales se réduisent à deux.

2-1 Méthode analytique

Dans le système de repère (Ox, Oy) le tenseur de contraintes s’écrit :

s t
s x xz
t
s
xz z
La condition de nulleté du moment résultant : tij = tji càd txz=tzx

Cours Mecanique des sol 29 Khaled MEFTAH

Connaissant les contraintes sur les facettes de normales ox et oz, on peut
déterminer les contraintes sur n’importe qu’elle autre facette inclinée d’un
angle « q »

q
t
sq
s

t
q

s
t

fig 3.2

Si l’on écrit la première condition d’équlibre ( somme des forces est nulle), on
aura l’état de contrainte sur le plan incliné de « q »
sq s x +s z + s z -s x cos2q -t xz sin2q (3.1)
=
2 2
tq s z -s x sin 2q +t xz cos2q (3.2)
=
2

Le lieu de contraintes dans le plan (s,t) est défini par la relation :

s x +s z s z -s x (3.3)
s
q - ²+t²
q
= ²+t²xz

2 2

C’est l’équation d’un cercle (cercle de Mohr):

- de centre de coordonnées ((sx+sz)/2, 0)
s z -s x 2 2
- de rayon R= ( ) +t xz
2

L’orientation des plans principaux est obtenue pour tq=0, soit :

q =- 1arctg 2t xz et q 2 =q +p /2 (3.4)
1 2 s z -s x 1

Il existe donc deux plans principaux dont l’orientation est donnée par q1 et q2
Les contraintes principales majeure et mineure sont déterminées à partir de
l’équation du cercle
Cours Mecanique des sol 30 Khaled MEFTAH
 ss x + z s z -s x 2 2 (3.5)
= s + ) +t
(
1 2 2 xz

s x +s z (3.6)
s3 = (s z -s x )2 +t xz 2
-
2 2
t
t

tq
q q
0 s s sq s
-t

fig3.3 : Cercle de Mohr

A noter que, si les directions x et z sont principales (sx = s3 ; sz = s1 et

txz=0) on trouve :
sq s1+s3 + s1-s3 cos2q (3.7)
=
2 2
t =s1-s3 sin2q (3.8)
q 2
2- 2 Méthode graphique
Il s’agit de déterminer l’état de contraintes sur le plan incliné d’un angleq et dont
les valeurs des contraintes principales s1 et s3 sont connues (fig 3.4)

sq
tq
s

s
fig3.4

Cours Mecanique des sol 31 Khaled MEFTAH

 La démarche utilisée pour résoudre ce problème est la suivante :
- De s1, on trace une parallèle au plan de s1
- De s3, on trace une parallèle au plan de s3
- L’intersection des deux plans donne le pôle « P »
- Du pôle « P », on trace la parallèle à la facette sur laquelle on veut
trouver l’état de contraintes (sq et tq)
- L’intersection de cette droite avec le cercle donne sq et tq

t
q

e
n d
a
l

tq
Plan sde3

q s
0 P s3 sq s1
Plan de s1
(s s

fig 3.5
On doit souligner enfin que, en mécanique des sols, on adopte la convention
de signes suivante :
- s >0 en compression
- s <0 en traction
3- Contraintes dues au poids propre du sol
Le poids du sol augmente avec la profondeur ; réparti sur une unit é de surface
horizontale à une profondeur donnée, il correspond à la pression ou contrainte
due au poids propre.
Pour un sol de poids volumique g (en kN/m3), et à une profondeur z (en m), la
contrainte verticale est :

sV = g .z (3.9)

3-1 Cas d’un sol sec

Cours Mecanique des sol 32 Khaled MEFTAH

 Le poids volumique intervenant dans le calcul de la contrainte est gd.
Dans le cas d’un sol stratifié en plusieurs couches de différents poids
volumiques et différentes hauteurs :
sV =∑gd i.hi (3.10)
3-2 Cas d’un sol saturé
Noyé dans l’eau interstitielle, l’élément de sol est allégé d’un poids équivalent au
poids volumique de l’eau multiplié par son propre volume. Ainsi la contrainte

totale sera (sv ) supportée par le squelette solide et l’eau.

D’après la loi de Terzaghi :

sv = s’ v + u et t’ = t (3.11)

u: pression interstitielle, u = gw.h w

s’ v : contrainte effective transmise au squelette solide.
NB : Le poids volumique intervenant dans le calcul de la contrainte totale est
.
gsat
- s’v = s v – u = g’.z (g’: poids volumique déjaugé)
Exemple : Traçons les diagrammes de variation de s v , s’v et u en
fonction de la profondeur

s s

g
s g

s g

s g

fig 3.6 Diagrammes de Variation des contraintes totales, effectives et

interstitielles en fonction de la profondeur.

Cours Mecanique des sol 33 Khaled MEFTAH

 4- Contraintes dues aux surcharges

Les dépôts de sol sont normalement stables, à moins que des circonstances
naturelles ou un chargement artificiel ne contribuent à y accroître les
contraintes effectives et qu’un tassement s’en résulte. On sait qu’un
abaissement de la nappe augmente la contrainte effective, mais divers types
de surcharges induisent également des contraintes ( Ds) dans le sol. Il s’agit :
- Des charges ponctuelles.
- Des charges uniformément réparties sur les surfaces rectangulaires et
circulaires.
- Des charges en forme de remblai de longueur supposée infinie

4-1 Cas d’une surcharge uniformement répartie sur toute la surface q

Dans ce cas et quelle que soit la profondeur z, on a :
Ds =q (3.12)

4-2 Cas d’une surcharge ponctuelle Q

En considérant le sol comme milieu semi-infini élastique non pesant, la
contrainte verticale due à la force ponctuelle Q est calculée d ’après la formule de
Boussinesq :

Ds 3 . 3 (3.13)
z
= Q
2p r ²+z² 5/2

Ds

fig 3.7

Cours Mecanique des sol 34 Khaled MEFTAH

 Cette èquation peut s’écrire sous une forme plus simple :
2
Ds = Iz (Q/z ).

(3.14)
où Iz est un facteur d’influence de contrainte verticale déterminé en fonction
de r/z par des abaques.

4-3 Cas d’une surcharge circulaire uniforme q

Ds

fig 3.8
3
z )
Ds = q(1 -
(
+ ) (3.15)
r ² z²
Ou bien : Ds = Iz q (Iz : facteur d’influence fonction r/R et z/R ).

4-4 Cas d’une surcharge rectangulaire uniforme q

Sous l’effet d’une charge rectangulaire de largeur « b » et de longueur « l »,

la contrainte induite Ds sous l’un des coins de cette charge, est donnée par :
Ds = Iz q (3.16)
Iz : facteur d’influence fonction de b/z et l/z. donn par le tableau 3.1.

Cours Mecanique des sol 35 Khaled MEFTAH

 Tableau 3.1 valeurs de Iz pour une semelle rectangulaire b.l
b/z 0.1 0.3 0.5 1 2 >3

l/z

0.1 0.01 0.01 0.02 0.03 0.03 0.03

0.3 0.01 0.04 0.06 0.08 0.09 0.09
0.5 0.02 0.06 0.08 0.12 0.13 0.14
0.7 0.02 0.07 0.10 0.14 0.17 0.17
1 0.03 0.08 0.12 0.18 0.20 0.20
1.2 0.03 0.09 0.13 0.18 0.21 0.22
1.4 0.03 0.09 0.13 0.19 0.22 0.22
1.6 0.03 0.09 0.13 0.20 0.22 0.23
2 0.03 0.09 0.13 0.20 0.24 0.25

La contrainte à la verticale d’un point quelconque s’obtient en construisant à

partir du rectangle et du point, des rectangles ayant chacun un sommet au
point considéré. La contrainte cherchée est la somme algébrique des
contraintes produites par les rectangles.

exemple :
Pour calculer Ds à la verticale du point A sous l’effet de la pression « q » de la semelle
EHDC, on utilise la méthode de découpage

Ds =[ Iz1 – (Iz 2 +Iz3)+ Iz4]q

D C B

Semelle
(EHDC)
I1(DFAB) I3(ABCG)

E H I

I2(AIEF) I4(AGHI)
F G A

fig 3.9

Cours Mecanique des sol 36 Khaled MEFTAH

 4-5 Cas d’un remblai semi-infini
Pour un remblai de hauteur Hr et de poids volumique gr, la contrainte verticale
est : Ds = Iz q (3.17)
Avec : q = gr . Hr
Iz : facteur d’influence donné par le tableau 3.2 en fonction de a/z et b/z.

a b

Hr

Ds

fig 3.10
Tableau 3.2 valeurs de Iz pour un remblai semi-infini

a/z 0.01 0.05 0.1 0.3 0.5 1

b/z

0.0 0.00 0.01 0.03 0.10 0.15 0.26

0.2 0.13 0.14 0.16 0.22 0.25 0.33
0.4 0.23 0.24 0.25 0.30 0.33 0.38
0.6 0.32 0.32 0.33 0.36 0.38 0.41
0.8 0.37 0.37 0.38 0.40 0.41 0.45
1.0 0.41 0.41 0.42 0.43 0.44 0.45
1.2 0.44 0.44 0.44 0.45 0.46 0.47
1.4 0.45 0.45 0.46 0.46 0.47 0.48
1.6 0.47 0.47 0.47 0.47 0.48 0.49
2.0 0.48 0.48 0.48 0.48 0.49 0.49
3.0 0.49 0.49 0.49 0.49 0.49 0.50

Cours Mecanique des sol 37 Khaled MEFTAH

 QUESTIONS A DEBATTRE
1. Pourquoi les contraintes verticales et horizontales sont-elles la plupart du
temps associées aux contraintes principales ?
2. Que représente le cercle de Mohr ?
3. Pourqoui un sol tasse suite à un rabattement de la nappe ?
4. Une nouvelle fondation, peut-elle transmettre des charges à une fondation
proche existante ?
5. Si une surcharge est appliquée à la surface d’un sol, comment peut on
évaluer son influence à une profondeur donnée ?
Exercice1
Détrminer analytiquement et graphiquement les contraintes qui se
développent sur le plan incliné de 60°:

? ?

Exercice2
On donne l’état de contraintes illustré sur la figure ci-dessous.
On demande de déterminer :
1) les contraintes principales et leurs directions
2) L’état de contraintes sur les plans horizantal et vertical
3) La contrainte de cisaillement maximale positive et le plan sur lequel elle
est appliquée. En déduire l’angle entre ce plan et le plan de s3 ( ous1 )

Cours Mecanique des sol 38 Khaled MEFTAH

 Exercice3
Tracer les diagrammes de variation des contraintes totales et effectives et des
pressions interstitielles, en fonction de la profondeur pour les deux cas
suivants :
a) Le niveau de la nappe est à 1m au-dessous du terrain natuel.
b) On rabat la nappe jusqu’à une profondeur de 5m à partir du terrain
naturel. Entre 1 et 5m, le sol à un degré de saturation de 50% que peut
on conclure ?.

T.N
g

Cours Mecanique des sol 39 Khaled MEFTAH

 Exercice4
On donne le profil d’un sol de fondation de 0 à 22m.

T.N

1) 1-Tracer les diagrammes de variation des contraintes totales,

effectives et les pressions interstitielles de 0 à 22m.

2) 2-Calculer le supplément de contraintes dues au rabattement de la

nappe à 6m du T.N.
3) On construit à la surface du sol un bâtiment ayant la forme ci dessous.
Déterminer les valeurs des contraintes dues à ce batiment, aux
profondeurs 6m et 18m au dessous de la base de la fondation pour
les verticales passants par A et B.

Cours Mecanique des sol 40 Khaled MEFTAH

 Exercice 5
Une semelle carrée 4mx4m porte une charge totale de1520KN. Déterminer
l’accroissement de la contrainte à 4m et à 8m de profondeur :
a) Sous un coin de la semelle.
b) Sous le centre de la semelle.
c) Sous les points milieux des deux côtés de la semelle.

Exercice 6
Soit la coupe géotechnique du sol représenté par la figure ci-après.

T.N
g

1) Tracer les diagrammes des contraintes totales, effectives et

interstitielles. On considère que l’argile sableuse est sèche entre 0 et 2m, le
sable est partiellement saturé (Sr=0.6) entre 2 et 4m.
2) On construit à la surface du sol un remblai représenté par la figure ci-
dessous. Calculer alors les excès de contraintes dûs à ce remblai aux
profondeurs z=4m et z=7m sous les verticales passant par le point A , le point
B (crête du talus) , le point C (pieds du talus) et le point D.

Cours Mecanique des sol 41 Khaled MEFTAH

 g

ELEMENTS DE CORRECTION
EXERCICE 1
Analytiquement :
s3 = 30 Kpa et s1 = 100kPa
Pour déterminer sq et tq on utilise les relations (3.7) et (3.8)
sq s1+s3 + s1-s3 cos2q et tq sin2
s 1- s 3
= = q
2 2 2
sq = 47.4kPa et tq = 30.2kPa

Graphiquement :

t
t

0 s s

Cours Mecanique des sol 42 Khaled MEFTAH

 Exercice 4
1- Entre 0 et 4m le sable est sec :
3
gd =(1-h)Gsg w =13.5kN/m
Entre 4 et 8m le sable est saturé
3
g = 1-h .Gs+h g
w =18.5kN/m
Profondeur sv (kPa) u (kPa) sv’ (kPa)
4 54 0 54
8 128 40 88
14 242 100 142
22 410 180 230

s s

Diagrammes de variation des contraintes totale, effective et de

pression interstitielle
2) L’orsque la nappe s’abaisse de 2m, on a :
- sv reste constante (puisque les poids volumiques restent inchangés)
- u diminue de gw .2= 20kPa.
- s’v augmente de gw .2= 20kPa.
4) Le supplément de contraintes pour q=200kPa :

Cours Mecanique des sol 43 Khaled MEFTAH

 Au point A : Ds = q ( Iz1 +Iz2 +Iz3)
2m 4m
6m 1 2

2m 3

Au point B: Ds = q ( Iz4 -Iz5)

4 6m

5 2m

2m 4m

Sachant que : Iz1 = f(2/z,6/z)

: Iz2 = f(4/z,6/z)
: Iz3 = f(2/z,2/z)
: Iz4 = f(6/z,8/z)
: Iz5 = f(4/z,2/z)

Cours Mecanique des sol 44 Khaled MEFTAH

 CHAPITRE4

TASSEMENT ET CONSOLIDATION DES SOLS

1- Notions de déformation
Sous l’application de charges, le sol comme tout solide se déforme. Pour
déterminer les déformations qui ont lieu dans toutes les directions autour du
point M du sol, il suffit de connaître le svaleurs des déformations dans les
directions Ox,Oy et Oz autour de ce point. On définit ainsi le tenseur de
déformations :

e x e yxe zx
e: e xy e y e zy
e xze yz e z

Les déformations sont reliées aux déplacements u,v,w par les relations :
- Elongations : ex = ¶u/¶x ; ey = ¶u/¶y et ez = ¶u/¶z
- Distorsions : gxy = 2 exy = ¶u/¶y + ¶v/¶x

gxz = 2 exz = ¶u/¶z + ¶w/¶x

gyz = 2 eyz = ¶u/¶z + ¶w/¶y

ei j sont les deformations de cisaillement
La variation de volume du petit élément autour du point M est :
DV/V = ex + ey + ez (4.1)

Il existe aussi trois directions principales pour les quelles les déformations
angulaires son nulles (gi j =0). Ces directions sont appelées directions
principales de déformation, et les déformations principales sont notées :e1 , e2
et e3.
2- Relations contraintes - déformations
Ces relations sont appelées « lois de comportement » puisqu’elles permettent
de caractériser la réponse d’un matériau sous l’effet d’un chargement.
Dans le domaine de déformations élastiques dans un solide isotrope, les
relations entre les contraintes et les déformations (loi Hooke) sont :

Cours Mecanique des sol 45 Khaled MEFTAH

 ex 1 [s x -n(s y +s z )] gxz= 1 .t xz
=
E' G
ey 1 [sy-n(s x +s z )] gyz .t yz
1 (4.2)
= = G
E'
ez 1 [s z -n(s x +s y )] gxy= 1 .t xy
=
E' G

E : module d’élasticité longitudinal.

n : coefficient de Poisson.
G : module de cisaillement transversal.
G = E/[2(1+n)]
E et n peuvent être déterminés à partir des résultats d’essais en laboratoire ou
in-situ.

3- Tassements des sols-Consolidation

Sous l’effet d’un chargement donné (fondation, remblai, etc..), le sol se
déforme. On sait que dans la plus part des cas, la surface du sol est
horizantale et les charges sont verticales; les déformations et par conséquent
les déplacements, seront dans la même direction. Ils sont appelés tassements.

Pour un sol, les tassements résultent essentiellement de sa compressibilité

(diminution de volume) qui est dû :
- à la compression du squelette solide,
- à l’évacuation de l’eau contnu dans les vides,
- et à la compression de l’eau et de l’air contenus dans les vides.
A noter que pour les contraintes courantes l’eau et le squelette solide peuvent
être considérés incompressibles.

Dans le cas des sols grénus (sable et gravier ayant un coefficient de

perméabilité élevé), saturés ou non, le tassement est immédiatDhi.

Pour les sols fins saturés (faible coefficient de perméabilité), sous l’action d’une
charge, l’eau libre ne peut s’évacuer immédiatement et supporte toutes
Cours Mecanique des sol 46 Khaled MEFTAH
 les contraintes appliquées ( suppressions interstitielles Du=Ds) pendant la
phase de construction de l’ouvrage ; on aura le tassement immédiat Dhi.
La transmission des contraintes au squelette solide se fait progressivement au
cours du drainage de l’eau et les surpressions interstitielles diminuent. Cet
écoulement s’arrête lorsque Du s’annule; on obtient donc le tassement à long

terme ou le tassement final de consolidation primaire Dhc

s s', u
s

D
D
D
s=s

fig 4.1
à tinfini : Dh = Dhi + Dhc

à un instant (t) : Dh(t)= Dhi + Dhc(t) (4.3)

3.1 Relations entre le tassement, l’indice des vides et la contrainte effective

Pour une couche de sol de hauteur « h » et d’indice des vides initial « e 0 »,

après un chargement donné et à un instant « t », on a ;
D = eo - e = De (4.4)
h
h1+ eo 1+ eo
Dh et e sont le tassement et l’indice des vides à l’instant « t »

Cours Mecanique des sol 47 Khaled MEFTAH

 s : due à
Pour déterminer la relation entre l’indice des vides et la contrainte «
l’action des charges », on doit réaliser un essai oedomètrique qui permet
d’étuier la consolidation des sols ( amplitudes et durée des tassements pour
une charge donnée)

L’ essai consiste à placer un échantillon de sol saturé dans un moule

cylindrique indéformable (module oedométrique) et de le charger verticalement
(charge constante) jusqu’à dissipation des surpressions interstitielles, tout en
mesurant les tassements.

En appliquant un chargement discontinu par paliers et en déterminant la

contrainte effective s’ et l’indice des vides « e » pour chaque palier de

chargement, on peut tracer la courbe : e= f (logs’) appelée courbe
oedométrique. Cette courbe peut être schématisée par la figure 4.3 suivante.

A Courbe : e=f(logs’)

B
C

h
a
r
g
e
m
e
D c
é
n
t
h
ar
g
e
m

nt
e

C
0 s s

fig 4.3
On déduit de cette courbe :
- L’indice de compression Cc De (4.5)
(pente de la droite BC): =-
D logs'
- Cs De (4.6)
L’indice de gonflement ( pente de la droite DC): =-
D logs'
Cours Mecanique des sol 48 Khaled MEFTAH
 - La pression de préconsolidation :sp’ C’est la plus grande
contrainte qu’a supporté le sol durant son histoire géologique.

Le rapport de surconsolidation (Roc)nous permet de connaître l’état de consolidation

du sol. Il est défini par :

Roc s'p (4.7)

= s'o
s0’ : Contrainte effective vertivale du sol en place
· Si Roc = 1, le sol est normalement consolidé, c.à.d qu’il n’a
jamais été le ite d’une érosion ou supporté des surcharges
(Glaciers..) qui en suite ont disparue. Ces sols sont sujettis
à développer des grands tassements, puisque toute
surcharge peut augmenter la contrainte effective à un
niveau que le sol n’a jamais atteint, ce qui pourra réduire
considérablement l’indice des vides.
· Si Roc > 1, le sol est surconsolidé. Dans ce cas les
tassements seront très faibles tant que la contrainte
effective finale est inférieure à la pression de
préconsolidation. Dans le cas contraire les tassements
seront plus prononcés.

· Si Roc < 1 , le sol est sousconsolidé. C’est le cas des sols

en cours de formation (vases,etc..) .

4- Evolution du tassement au cours du temps

4-2 Tassement final de consolidation primaire

Pour une couche de sol de hauteur « h », d’indice de compression « Cc »,

d’indice de gonflement «Cs » et de pression de préconsolidation « sp’ », si la

contrainte verticale initiale due aux poids propres du sol est « s0’ », calculée
au centre de la couche et si le chargement provoque un supplément de
contrainte Ds (défini dans le chapitre3), les tassements se calculent comme
suit :

Cours Mecanique des sol 49 Khaled MEFTAH

 · Roc = 1 : le sol est normalement consolidé :
Cc so'+Ds )
Dhc 1+e (4.8)
= o hlog ( s' o

· Roc > 1 : le sol est sur consolidé :

Deux cas peuvent se présenter :

- Si s0’ +Ds > sp’ alors :
Cs s p' C c o
' +D s
Dhc = hlog + hlog
s

(4.9 a)
1+eo s'0 1+eo s'p

- Si s0’ +Ds < sp’ alors :

Dhc Cc s o'+Ds
hlog (4.9 b)
= s'p
1+eo

· Roc < 1 : le sol est sous consolidé :

Dhc =
C
c hlog (s '+Ds )
o
(4.10)
1+eo s'o

On doit noter que :

- Si le sol est composé de plusieurs couches compressibles, le tassement total
sera la somme des tassements de chaque couche.
- Si la hauteur de la couche est supérieure à 3m, on subdivisera la couche en
sous couches.

4-3 Tassement à un instant « t »

le tassement à un temps donné Dhc(t) est calculé en fonction du tassement

final Dhc, à partir de la formule suivante :
Dhc(t) = U. Dhc (4.11)

U : étant le degré de consolidation moyen défini à un instant « t »par la

relation :

Cours Mecanique des sol 50 Khaled MEFTAH

 h
∫ u ( z ,t ) dz
D u (t )
U =1- 0 =1-
h D u 0 (4.12)
∫ u ( z , 0 ) dz
0
Du(t) : surpression interstitielle au temps t.
Du(o) : surpression interstitielle à l’instant initiale (Ds).
Pour déterminer le degré de consolidation « U », il suffit de déterminer
l’équation de u(z,t).
La théorie de consolidation unidimensionnelle de Terzagi, basée sur des
hypothèses simplifiées nous donne l’équation de consolidation en fonction du

coefficient de consolidation : Cv (déterminé à partir de la courbe du tassement

en fonction du temps pour un des paliers de chargement de l’essai
oedométrique).
On peut alors déterminer le degré de consolidation U pour un temps « t », ou
plutôt pour un facteur temps « Tv » fonction du temps, Cv et des conditions de
drainage :
(4.13)
T = Cv.t
v
H' 2
H’ : représente le chemin de drainage.
Pour une couche de hauteur h :
H’ = h/2 si la couche est doublement drainée.
H’ = h si couche est simplement drainée.
Le tableau suivant donne les valeurs de U en fonction deTv

U (%) Tv U (%) Tv

12 0,01 70 0,40
20 0,03 80 0,57
30 0,07 90 0,85
40 0,13 95 1,20
50 0,20 99 2,00
60 0,29 100 2,08
Notons que pour la fin de la consolidation primaire (U=100%), on prendra
Tv=2.08 ( correspondant à un degré de consolidation de 99.5%)

Cours Mecanique des sol 51 Khaled MEFTAH

 5- Contrainte effective et pression interstitielle au cours du temps

A une profondeur « z » et a un instant « t » donnés, le degré de consolidation

Uz est calculé à partir de la relation :
u(z,t)
U z =1- (4.14)
u(z,0) ,
Uz est déterminé en fonction de Tv et Z/H’.
(Z est comptée à partir du niveau supérieur de la couche compressible
étudiée).
La pression interstitielle : u(z,t)=(1-Uz)Ds(z).
La contrainte effective : s’(z,t) = Uz. Ds(z), puisque u(z,0)= Ds(z).

QUESTIONS A DEBATTRE
1- Pourquoi les tassements surviennent-ils plus rapidement dans les sols à
gros grains que dans les sols à grains fins ?
2- Comment appelle-t-on un tassement qui évolue avec le temps ?
3- Résumez les principales étapes de l’essai de consolidation.
4- Lors de l’essai de consolidation, que représente la pression de pré
consolidation déterminée graphiquement sur la courbe.
5- Quelle est l’ampleur du tassement auquel on peut s’attendre lorsque la
pression exercée par une fondation induit un contrainte effective totale
supérieure à la pression de pré consolidation ?
6- A quoi correspond le degré de consolidation ? Pour quoi on Tv=2.08 pour
U=100% ?
7- Comment le mode de drainage d’un dépôt d’argile influe-t-il sur le temps de
consolidation ?

Cours Mecanique des sol 52 Khaled MEFTAH

Exercice1

Suite à un sondage à ciel ouvert, on donne une représentation de la coupe

d’un sol.

Les résultats d’essais oedométriques réalisés sur des eprouvettes de la couche

d’argile raide sont :
-7 2
e0= 0.8; Cc= 0.26; s’p=70kPa et Cv=2.10 m /s
1) Pour une largeur de la fondation B=3m, calculer le tassement final de
consolidation primaire dû au remblai et à la pression « q » appliquée par la
semelle.
On prendra : q= ((Q/B) + 20x2.5) en kPa, et on la suppose uniformément
répartie sur une bande infinie.
2) Combien de temps après l’application des charges obtient-t-on 30%,50% et
100% de consolidation.
3) Calculer, lorsque U=50%, les valeurs du tassement et de l’indice des vides.

Exercice2

A l’aide des données de la figure suivante, calculer le tassement de

consolidation d’un dépôt d’argile sous le centre d’un remblai en longueur de
4m de hauteur, sans tenir comte du tassement immédiat d’une couche de
gravier. L’essai de consolidation effectué sur un échantillon prélevé à mi-

Cours Mecanique des sol 53 Khaled MEFTAH

 hauteur du dépôt argileux a révélé une pression de pré consolidation de 110kPa
,un coefficients de gonflement de 0.023 et de compression de 0.38.

Exercice3
Estimer le temps nécessaire pour qu’un dépôt argileux de 6m de hauteur situé
entre deux couches de sables perméable atteigne 95% de son tassement de
consolidation. L’essai de consolidation a fourni un coefficient de consolidation
-3 2
de 0.5 10 cm /s.

Exercice4
Un dépôt d’argile de 3.5m d’épaisseur repose sur une marne imperméable
très compacte.
-3 2
Si le coefficient de consolidation est de 0.625 10 cm /s, quel sera le degré
de consolidation du dépôt sous la pression provoquant le tassemnt :
1) après six mois.
2) après un an

Cours Mecanique des sol 54 Khaled MEFTAH

 ELEMENTS DE CORRECTION
Exercice 2
1)
h = 5m
s’o = 18x1 +3.5x11=56.5kPa.
s’p = 70kPa.
Dhc Cc so'+Ds
hlog
=
1+eo s'p
Ds = qr h r + 2 q Iz sachant que q = 383.33 kPa.
Iz = f (b/z ,l/z) b = 1.50m ; l = infini (semelle filante)
Z = 2.50m
D’après l’abaque 3-3, on aura Iz = 0.157
Ds = 150.36 kPa.

D’où : Dhc = 0.322m.

2)
-7 2
Tv Cv.t avec : Cv=2.10 m /s ; H’ = h/2 (couche simplement drainée)
= H' 2

U 30% 50% 100%

Tv 0.07 0.2 2.08

Temps (ans) 0.28 0.8 8.32
3)
Lorsque U = 50%, le temps correspondant est environ de 9mois et 18jours.
Dht =0.5 x 32.2 = 16.1 cm
et D/(1+eo) = Dh/h avec De = eo – e
e = eo –(Dh/h) (1+eo) = 0.74
e = 0.74

Cours Mecanique des sol 55 Khaled MEFTAH

 CHAPITRE 5

LA RESISTANCE AU CISAILLEMENT DES SOLS

1- Comportement élastoplastique des sols

Dans les chapitres précédents, on a étudié le comportement des sols sous

faibles taux de chargement (domaine élastique). Dans ce chapitre, nous
traitons les grandes déformations, c’est à dire le comportement des sols à l’état
d’écoulement plastique ou à l’état de rupture.

La loi de Hooke n’est plus valable, on utilise alors une nouvelle loi appelée :
- Critère d’écoulement plastique qui représente la frontière du domaine
d’élasticité.
- Ou le critère de rupture représenté par la courbe intrinsèque qui est
l’enveloppe des cercles de Mohr correspondant à la rupture.

Au moment de la rupture d’un sol, il y a un glissement entre les particules

solides, d’où le terme de résistance au cisaillement.

Tout sol présente une résistance au cisaillement qui est due au frottement
inter-granulaire (contact entre les grains) et aux forces d’attraction entre les
particules dans le cas des sols fins.
du à la dureté des grains et à l’état de surface de contact.

On définit aussi la rupture dans un sol à partir des courbes contraintes-

déformations dans des essais à déformation contrôlée.

Cours Mecanique des sol 56 Khaled MEFTAH

 Rupture 3
Rupture
2
Rupture
1
e
e
%

fig 5.1
1 : Comportement élastique parfaitement plastique.
2 : Comportement élasto-plastique écrouissable.
3 : Comportement élasto-plastique écrouissable .

La connaissance de la résistance au cisaillement d’un sol ( ou de sa courbe

intrinsèque modelisée généralement par le critère de rupture de Coulomb) sert
pour la résolution d’un grand nombre de problèmes de Génie Cicil :
- Capacité portante des fondations superficielles et profondes.
- Stabilité des ouvrages de soutènement.
- Stabilité des talus, pentes et barrages, etc..
2- Comportement à court terme et à long terme des sols
(cas des sols saturés)
Un sol saturé soumis à des forces, subit une variation de volume. L’eau est
incompressible, alors cette variation a lieu s’il y a mouvement de l’eau ( apport ou
export).
Dans le cas d’un sol grenus, on sait que sa perméabilité est élevée, donc
l’écoulement de l’eau est très rapide. Le comportement de ce sol, de même que sa
résistance, ne sont régis que par celui du squelette solide.
Dans le cas d’un sol fin saturé, la perméabilité est très faible, l’écoulement de
l’eau est très lent et le changement de volume se traduit par une surpression
(export de l’eau) ou une dépression (apport d’eau) dans l’eau interstitielle.

Cours Mecanique des sol 57 Khaled MEFTAH

 On considère donc deux comportements extrêmes du sol :
- Un comportement à court terme (C.T) : L’eau n’a pas eu encore le temps
de s’évacuer et le sol se déforme à volume constant. Il s’agit alors d’un
comportement non drainé (étude en contraintes totales).
- Un comportement à long terme (L.T) : Au bout d’un temps assez long
(fin de la consolidation primaire), l’eau libre s’est évacuée. Le
comportement du sol est régi uniquement par celui du squelette solide
(étude en contraintes effectives).

3- Détermination des paramètres de résistance des sols en laboratoire

Les principaux essais de laboratoire sont :
a) Essai de cisaillement direct.
b) Essai tri axial.
c) Essai de compression simple.

o Phase1 : Remise sous l’état de contraintes en places (saturation et

application de contraintes jusqu’à u=0).

o Phase 2 : Cisaillement proprement dit.

Selon le drainage pendant l’une ou l’autre de ces deux phases, on distingue

trois types d’essais :
· Essai non consolidé non drainé (UU).

· Essai consolidé non drainé (CU).

· Essai consolidé drainé (CD).

3-1 Essai de cisaillement direct

L’échantillon de sol, placé dans deux demi-boites qui peuvent glisser l’une par
rapport à l’autre, est soumis sur l’élément supérieur à une contrainte normale

(s =N/S, S :section horizontale de la boite de cisaillement), puis à un effort de

Cours Mecanique des sol 58 Khaled MEFTAH

 cisaillement T jusqu’à la rupture tout en mesurant le déplacement horizontal de
la demi boite.

Charge verticale:N
s=N/S

Charge
horizontale: T

Plan de t=N/S'
cisaillement

fig 5.2 Boite de cisaillement

Trois échantillons identiques doivent être testés, pour trois contraintaintes
normales différentes
Les valeurs des contraintes de cisaillement à la rupture sont représentées en
fonction des contraintes normales t = f (s).
Cette représentation n’est autre que la traduction graphique de l’équation de
Coulomb : trupt = C + s rupt tg j (5.1)

a) Essai Consolidé Drainé(C.D)

j
t
g
t

t3
s
+
'
C
=
t
t2
t1

j
s s s
s

Fig 5.3

Cours Mecanique des sol 59 Khaled MEFTAH

 c’,j’ : Les paramètresrésistance au cisaillement
c’ : Cohésion non drainée.
j’ : Angle de frottement interne effectif.
b) Essai non consolidé non drainé(U.U)
Essai réalisable pour les sols fins

t
t t2 1t1 t =
3 Cu

Cu

s s s

Fig 5.4
Cu : Résistance au cisaillement non drainé.
3-2 Essai triaxial :
Pour trois éprouvettes identiques ( Æ=36 mm, h=2Æ) d’un sol donné, on
applique pour chacune les états de contraintes représentés comme suit :

PHASE1 PHASE2

s3
Ds1 Effort du
piston
s3 s3

s3 Ds1

s3 : Pression hydrostatique = s’ v0

Fig 5.5 : Champs de contraintes à l’appareil triaxial

Cours Mecanique des sol 60 Khaled MEFTAH

 A la rupture, on note les valeurs des contraintes principales :
s1 = s3 + Ds1 et s3
Ds1 = s1 - s3 étant le déviateur de contraintes.

a) Essai non consolidé non drainé(UU).

L’essai étant rapide, la représentation ne peut se faire qu’en contraintes
totales (fig 5.6)

t
t = Cu

Cu

s
a

Fig 5.6

Essai 1 : sI3 et sI1 ; Essai 1 : sII3 et sII1 ; Essai 1 : sIII 3 et sIII 1

L’orientation du plan de rupture : a =p/4
Cu= (s1 - s3)/2 (5.2)
NB : Cet essai ne permet pas la détermination des paramètres effectifs même en
mesurant la pression interstitielle.
Le critère de rupture est :
- En terme de contraintes appliquées sur le plan de rupture :
t = Cu s =(s1 + s3)/2
- En terme de contraintes principales : Cu =(s1 - s3)/2

b) Essai Consolidé Drainé (C.D)

(Essai très lent : Contraintes effectives)

Cours Mecanique des sol 61 Khaled MEFTAH

 sj j
t
g

+
'
C
=

t
t
t
t
s s s s

FIG5.7 RSULTATS DE L’ESSAI CD A L’APPAREIL

TRIAXIAL

Orientation du plan de rupture : a =p/4 + j’/2

Le critère de rupture est :
En terme de contraintes sur le plan de rupture : t = c’ + s’tg j’
En terme de contraintes principales :
s’1 = s’3 Kp +2 c’ÖKp. (5.3)

Kp = tg² (p/4 +j’/2)

c) Essai Consolidé non Drainé (C.U)

Il a pour but :
- la détermination des paramètres de résistance effectifs (c’,j’) en mesurant
la pression interstitielle à la rupture.
- La détermination des paramètres de résistance consolidés non drainés
(Ccu, jcu)
- L’étude de la variation de la cohésion non drainée en fonction de la pression
de consolidation.

La représentation des résultats peut se faire en contraintes totales et en

contraintes effectives ( fig 5.8)

Cours Mecanique des sol 62 Khaled MEFTAH

 t
c
u
j
g
t

j s

t
g
c
u +
s j= C
+ t

jcu
'

Cu3 =
c

t
Cu2
Cu1

ccu
s
Fig 5.8
s
t
'

l c
l
=a +
cu 3
c

u
cu 2
Dcu
cu 1
Ds'c
a
sI3 sII3 sIII3 s
Fig 5.9 cu= f (s ‘c)

L’orientation du plan de rupture : a =p/4 + j’/2

Le critère de rupture est :
- En terme des contraintes effectives :
s’1 = s’3Kp + c’ÖKp. Kp = tg² (p/4 +j’/2) (5.4)
t = c’ + s’tg j’
- En terme des contraintes totales :
s1 = s3Kp +2 CcuÖKp. Kp = tg² (p/4 +jcu/2) (5.5)
t = Ccu + stg jcu
L’accroissement de la cohésion non drainée est : DCu = l D s’c (5.6)
D s’c : l’accroissement de la pression de consolidation.

Cours Mecanique des sol 63 Khaled MEFTAH

 3-3 Essai de compression simple
L’essai consiste à appliquer une charge axiale sur l’échantillon du sol et
l’augmenter progressivement jusqu’à la rupture. L’essai étant rapide, les
résultats sont représentées en contraintes totales (fig 5.10)

s1

a s1

t
t = Cu
a

s s1
s
fig 5.1 Résultats de l’essai de compression simple

Orientation du plan de rupture : a =p/4

La cohésion non drainée est : cu = s1 / 2 (5.7)

Cours Mecanique des sol 64 Khaled MEFTAH

 QUESTIONS A DEBATTRE

1- Selon le critère de Coulomb, de quoi dépend la résistance au cisaillement

d’un sol à la rupture ?
2- Quelles sont les circonstances qui provoquent généralement une rupture du
sol ?
3- Pourquoi utilise-t-on l’angle de frottement interne effectif pour les sols
pulvérulents ?
4- Dans un sol cohérent, qu’arrive-t-il à la valeur des paramètres de la
résistance au cisaillement lorsqu’il n’y a aucun drainage.
5- A quoi servent les deux phases lors des essais de cisaillement direct et
triaxial ?
6- Quels sont les paramètres de résistance au cisaillement qu’on doit utiliser
lors de la construction d’un ouvrage, et lorsqu’il est en service.
7- Parmi les essais de cisaillement, on parle du scissomètre de laboratoire et
du pénétromètre à cône. Comment sont -ils réalisés?

Exercice 1
Chercher les relations entre les contraintes principales et l’angle de frottement
interne dans les conditions où la cohésion est nulle.

Exercice 2
On a effectué des essais de cisaillement direct sur un sable moyen. A l ’aide
des résultats qui suivent, déterminer les paramètres de la résistance au cisaillement de ce sable.

Essai N° s’rupt trupt

1 67.3 57.1
2 118.6 99.2
3 170.2 150.7
4 231.8 194.7

Commenter les résultats trouvés.

Cours Mecanique des sol 65 Khaled MEFTAH

 Exercice 3
Les résultats de 2 essais triaxiaux consolidés non drainés (CU) avec mesure
de la pression interstitielle sur matériau de remblai, sont les suivants :

Essai N° s3 (kPa) s1-s3 (kPa) U(kPa)

1 105 466.5 20.70
2 307.5 1056 93.01

Déterminer les caractéristiques de résistance à long terme et les paramètres

de résistance consolidés non drainés.
Exercice 4
On prélève d’un sondage carotté à 11 m de profondeur, des carottes de sol
vaseux.
Les résultats d’essais sur les différentes éprouvettes de ce sol, sont fournis
dans les tableaux suivants :
- Essai de compression simple :
Rc = 100kPa.
- Essai triaxial UU.
s1 = 248 kPa s3 = 140 kPa
Essai de cisaillement direct drainé

Essai N° s’(kPa) t(kPa)

1 30 22
2 50 36
3 80 57
Essai tiaxiaux CU :

Essai N° s3 (kPa) s1 (kPa) U(kPa)

1 100 170 70.8

2 340 580 240

On demande :
1) L a cohésion non drainée de la vase.

Cours Mecanique des sol 66 Khaled MEFTAH

 2) Les paramètres de résistance effectifs déterminés à partir des essais de
cisaillement direct et des essais CU.
3) Les paramètres de résistance consolidés non drainés ainsi que la variation de
la cohésion non drainée en fonction de la pression de consolidation (s’c)
4) L’inclinaison des plans de rupture dans les éprouvettes :
a) de l’essai de cisaillement direct.( Essai 1)
b) de l’essai tiaxial CU (Essai 2)

Cours Mecanique des sol 67 Khaled MEFTAH

 CORRIGES

Exercice 1
La cohésion est nulle, donc la courbe intrinsèque passe par l’origine (t =tgj)

j
s3 s1
s
0 C

Le triangle OBC est rectangle

OC = (s1 + s3)/2 et BC = (s1 - s3)/2
= BC/OC d’où s3/ s1 =(1- sinj) / (1+ sinj) = tg² (p/4 -j/2)
de même s1/ s3 =(1+ sinj) / (1- sinj) = tg² (p/4 +j/2)
Exercice 2
1- Paramètre de résistance à court terme :
Essai de compression simple : Cu = Rc/2 = 50kPa.

Essai triaxial UU : Cu =(s1 - s3)/2 = (248-140)/2 = 54 kPa.

2- Paramètre de résistance effectifs :

· Cisaillement direct :
t = C’ + s’tg j’
22 = 30 tg j’ + c’ (1)
36 = 50 ’tg j’ + c’ (2)
57 = 80 ’tg j’ + c’ (3)
(1) et (2)⇒ j’ = 35° ; (2) et (3) j’ = 35°. c’ = 0

Cours Mecanique des sol 68 Khaled MEFTAH

 · Cisaillement triaxial :
s’1 = s’3Kp –2 c’ ÖKp.

Kp = tg² (p/4 +j’/2) = Nf

Essai1 : 99.2 = 29.2 Nf +2c’Ö Nf
Essai2 : 340 = 100 Nf +2c’Ö Nf
⇒ j’ = 33° et c’ @0

3- Paramètre : Ccu et jcu

Essai CU : Contraintes totales
s1 = s3Kp +2 CcuÖKp.
Kp = tg² (p/4 +jcu/2)
Essai1 : 170 = 100 Nf +2c’Ö Nf
Essai2 : 580 = 340 Nf +2c’Ö Nf
⇒ jcu = 15.1° et Ccu =0
Variation de la cohésion non drainée :
Essai1 : s3 = 100 kPa Cu = 35kPa

Essai2 : s3 = 340 kPa Cu = 120kPa

Cu = a + l s’c

CU = A + l . 100
Cu = a + l . 340
Cu=0.354 s’c

⇒ l = 0.354 et a =0
4- Inclinaison des plans de rupture :
· Eprouvette1, essai de cisaillement direct : Plan horizontal.
· Eprouvette 2, essai CU : a = p/4 +j’/2 = 61.5°

Cours Mecanique des sol 69 Khaled MEFTAH

 CHAPITRE 6

ETUDE DE LA PORTANCE DES

FONDATIONS SUPERFICIELLES A
PARTIR D’ESSAIS DE LABORATOIRE

1- Introduction

La fondation superficielle, élément inferieur d’un ouvrage, permet de

transmettre à une couche de sol peu profonde, dans des conditions favorables,
les charges issues de la superstructure.

Selon le chapitre B.1 du fascicule 62 titre V : En règle générale on peut

considérer qu’une fondation est superficielle lorsque sa hauteur
d’encastrement D est inférieure à 5 fois sa largeur.

De point de vue comportement mécanique, si l’on soumet une fondation à un

chargement croissant, le sol est en équilibre élastique au début de de
l’application de la charge; il le garde pour les valeurs faibles valeurs de celle-ci.
Lorsque la charge devient importante, le sol se met progressivement à
l’équilibre plastique et au moment de la rupture :

- Dans le cas d’une semelle, il se produit un coin rigide sous la base

de la semelle qui s’enfonce dans le sol en le refoulant de part et
d’autre( fig 6.1a)
- Dans le cas d’une fondation profonde (pieu) on obtient une zone de
rupture lovcale sous la pointe du pieu (fig 6.1b)

Cours Mecanique des sol 70 Khaled MEFTAH

 t
t

Fig 6.1-a Schéma de rupture sous une semelle chargée

fig 6.1-b Rupture sous un pieu

Un projet de fondation est donc très délicat, il doit répondre à trois sortes de
préoccupations :
- une bonne reconnaissance du sol support.
- les tassements doivent être inférieurs aux
tassements admissibles définis pour l’ouvrage
- les contraintes doivent être compatibles avec la
résistance à la rupture du sol : C’est le problème de la
capacité portante.

Cours Mecanique des sol 71 Khaled MEFTAH

 La contrainte de cisaillement le long du plan de rupture est donnée par :
- Pour un sol pulvérulent : t = s . tgj
- Pour un sol cohérent :
A court terme : t = Cu
A long terme : t = c’ + s’tg j’

La capacité portante est la pression maximale que peut supporter le sol

avant la rupture.
La contrainte admissible est la pression maximale qui puisse être appliquée
par une structure sur le sol, sans qu’il y ait des tassements excessifs et des
risques de rupture du sol

2- Calcul de capacité portante

Dans le cas d’une semelle filante de largeur B, encastrée à une profondeur D

dans le sol et soumise à une charge verticale centrée (fig 6.2), l ’équation générale
de la charge limite est :

1
q
2 2
γ BN γ cN q γD N
l= + c
+ + 1 q (6.1)

Ng, Nq, Nc : facteurs de portance donnés en fonction de l’angle de frottement

interne du sol j sous la base de la fondation

g
j g
fig 6.2

Cours Mecanique des sol 72 Khaled MEFTAH

2-1 Rupture à court terme et à long terme

A court terme, Les calculs sont menés en contraintes totales en utilisant cu et

ju = 0 (Ng =o, Nq =1 et Nc =5,14)

= + + (6.2)

A long terme, les calculs sont menés en contraintes effectives avec c’ et j’

q γ'BN c' N q γ' D

l 2 2 γ c + 1 (6.3)
+ + ( ) Nq
=

2-2 Effet de la présence d’une nappe

Il y a lieu de déjauger les poids volumiques si les sols correspondants sont

immergés :

o Pour une nappe affleurant à la surface du sol (sol saturé) :

1
qu g'BN l +c 'Nc + q+g' D Nq (6.4)
=
2 2 1
o Pour une nappe à grande profondeur (sol sec ou humide) :
1
g
qu =
2 BN l +cNc + q+g D Nq (6.5)
2 1

2-3 Influence de la forme de la semelle

L’équation générale de la capacité portante est modifiée par l’introduction des

coefficients de forme : Sg, Sc, Sq

q γ BN S cN c S q γDN
l 2 2 λ γ c S (6.6)
+ + (+ 1 q) q

=
Tableau 6.1 - Coefficients de forme
Fondations Rectangulaires Carrées Circulaires
Sg 1-0.2B/L 0.8 0.6
Sc 1+0.2B/L 1.2 1.3
Sq 1 1 1
Cours Mecanique des sol 73 Khaled MEFTAH
 2-4 Influence de l’excentrement de la charge

Si l’excentricité « e » est parallèle à « B » , on remplace dans tout ce qui

précède B par une largeur réduite B’ = B - 2.e
Si l’excentricité « e’ » est parallèle à « L » , on remplace L par une lorgueur
réduite L’ = L - 2.e’.

Figure 6.3

La charge de rupture QL est :

o QL =ql .B’. L’ (semelle rectangulaire ou carrée)
o QL = ql.p.B’B /4(semelle circulaire) (6.7)

2-5 Influence de l’inclinaison de la charge

L’équation générale de la capacité portante est modifiée par l’introduction des

coefficients minorateurs : ig, ic et iq
q γ BN i cN i q γDN i
u 2 2 λ + 1 q (6.8)
γ+
c
c+ ( ) q
=
Pour une inclinaison d de la charge par rapport à la verticale, le DTU 13.12
propose les relations suivantes pour les coefficients d’inclinaison :

o ig = (1-d/j’)² (6.9)
o ic= iq = (1-2d/p)²

Cours Mecanique des sol 74 Khaled MEFTAH

 2-6 Cas d’un sol bicouche

On utilise dans la pratique la méthode de la semelle fictive qui consiste à diffuser

la pression sous un angle de répartition de l’ordre de 27° (fig 6.4 )

q1 D
1 sol compact B
q2 h g1
2 sol mou
B1=B+h g2
c,j

Figure 6.4

On calcule la portance du sol (2) sous l’effet de la pression q 1 , avec :

q2
q γ
1 1 h1
= + + (6.10)

3- Calcul de la contrainte admissible

On prend pour la pression admissible :

qadm = charge des terres enlevées + capacité portante nette
FS
q Dγ. l
adm = + - (6.11)

Fs étant le coéfficient de sécurité pris généralement égal à 3.

On vérifie ensuite que la pression moyenne appliquée sous la fondation est
inférieure à qadm.

Cours Mecanique des sol 75 Khaled MEFTAH

 N.B : Dans le cas d’une semelle posée à une profondeur d en fond de fouille
après un premier creusement de D ( cave, sous-sol, etc …) : fig6.5

fig 6.5
q Dγ. l
adm

= + -

Avec γ' cNc

ql BN+ dγ.
+ ( )
= N
q

Cours Mecanique des sol 76 Khaled MEFTAH

 QUESTIONS A DEBATTRE
1- Disntinguer les divers types de fondations suprficielles?
2- Pour un ouvrage reposant sur des fondations de différents types, quel sorte
de tassement peut se produire ?
3- Si on demande votre avis (sommaire) sur le choix du type de fondation,
quels documents vous devez consulter avant de répondre.
4- Définisser les expressions « capacité portante admissible» et « capacité
portante à la rupture »
5- Quels sont les trois parties, ou termes, de l’équation générale de la capacité
portante ?
6- Enumérer les facteurs et les paramètres d’état du sol dont on tient compte
lors du calcul de capacité portante.
Exercice1
Calculer la capacié portante sous la semelle filante dans les deux cas :
- La nappe est au niveau (a).
- La nappe est au niveau (b).
L’angle de frottement interne du sol est 30°
D=2 m

4.00

g =18.9 kN/m3
502.
b B=1.50

g =22.1 kN/m3

Semelle filante reposant sur le sable

Cours Mecanique des sol 77 Khaled MEFTAH

 Exercice2
Une semelle filante de 1 m de largeur repose sur une argile saturée. Calculer
la capacité portante admissible de l’argile à court et à long terme en utilisant
un coefficient de sécurité égal à 3.

D=2.00
g =17.3 kN/m3
1.00

Cu=65 kPa

C'=16 kPa ; j'=25°

g =23.0 kN/m3
1.00

Semelle filante sur une argile saturée.

Exercice3
Au niveau d’un mur mitoyen, on était contraint de réaliser une semelle
rectangulaire ayant comme dimensions B = 1.20 m, et L= 1.00 m. La charge
transmise par le poteau est alors excentrée de 45 cm dans la direction
parallèle à B. Le sol support est pulvérulent de poids volumique 20 kN/m3 et
d’angle de fottement interne 35°.
Déterminer la capacité portante admissible de ce sol en utilisant un coefficient
de sécurité F= 3

0.30
D=1.60

g =20kN/m3
1.20

Semelle rectangulaire excentrée.

Cours Mecanique des sol 78 Khaled MEFTAH

 Exercice 4
Sur un sol très compact, on construit à 1.60 m de profondeur une semelle
filante de largeur B = 2 m produisant une charge Q = 250 kN/ml , inclinée de
30° et excentreé de 30 cm.
Calculer les coefficients de sécurité à court et à long terme.

12°
Q

30 g =18
601. kN/m3
Cu=165kPa ; C'=10kPa ;j=30°

B=2.00

Semelle filante avec charge excentrée et inclinée

Cours Mecanique des sol 79 Khaled MEFTAH

 ELEMENTS DE CORRECTION
Exercice1
Nappe en (a) :
La nappe phréatique se trouve sous la profondeur d’influence de la semelle
Hw> (D+B) puisque 4 > (2+1.5)

Pour j = 30°⇒ abaques 6-1, on a : Ng = 18.10, Nq = 18.4, Nc = 30.00

Les facteurs de formes pour une semelle filante sont égaux à 1. Pas de
facteurs d’inclinaison, ni d’excentricité. D’où :
qadm = (qnet/ Fs ) + g . D telleque qnet = ql - g .D et Fs =3
1
( )
ql = 2gBNl + g D Nq
ql = (1/2.x18.9 x 2 x 18.1) + (18.9x2x18.4) = 1037.61 kPa.
qnet = 1037.61 – (18.9x2) = 999.81 kPa
qadm = 999.81/3+(18.9x2) = 371.07 kPa
Nappe en (a) :
La nappe phréatique se trouve à une profondeur intermédiaire entre la
profondeur de la semelle et la profondeur d’influence de la semelle. Il faut alors
déjauger le sol dans le terme de surface.
1
( )
ql = 2g'BNl + g D Nq = 856.60 kPa

qnet = 856.6 – 8.9x2= 838.8 kPa

qadm = 297.4 kPa
Exercice2
A court terme:
ju = 0⇒ abaques 6-1 on a : Ng = 0, Nq = 1, Nc = 5.14
qadm = 120.7 kPa
A Long Terme:
j’ = 25°⇒ abaques 6-1 on a : Ng = 8.10, Nq = 10.70, Nc = 20.70
qadm = 227.6 kPa
Exercice3
A Court Terme : qadm = Q/1.4 Fs = 4.20
A long terme : Fs = 4.40

Cours Mecanique des sol 80 Khaled MEFTAH

 CHAPITRE 7

POUSSEE ET BUTEE DES TERRES

1- Introduction

Les sols exercent des pressions verticales sur les couches sous-jacentes, ce
qui engendre des pressions horizontales dites « pousseé des terres ». Un
ouvrage de soutènement doit être alors conçu de manière à résister à cette
poussée.
Pour un massif de sol, on distingue trois états d’équilibre :
- état au repos : le rapport entre la contrainte horizontale effective et la
contrainte verticale effective s’exprime au moyen du coefficient de
pression des terres au repos K0.

K0=s’h /s’v
(7.1)
Alors s’ =s’ . K0
h v
- état actif ou de poussée.

- état passif ou de butéé.

2- La théorie de Rankine

L’équilibre des contraintes dans un sol est exprimé par l’enveloppe de rupture

(courbe intrinsèque) issue de l’équation de Coulomb : t = c + s’ tg j.

L’état de contrainte agissant sur un élément de sol s(’h,s’v ) est représenté

par le cercle de Mohr. Rappelons que pour garder l’équilibre de l’élément de
sol, le cercle doit être toujours à l’interieur de l’enveloppe de rupture (fig 7.1)

Cours Mecanique des sol 81 Khaled MEFTAH

 t
etat
3 supérieur(butée)
j
s
etat inferieur t
(poussée)
2

1
s

0s s s
s
s s

Fig 7.1 Les Trois états d’équilibre d’un sol

Les cercles représentent trois états différents des contraintes d’un élément de
sol.
Cercle1 : exprime l’état d’équilibre au repos.La contrainte de cisaillement
agissant sur l’élement est inférieure à la résistance au cisaillement du sol : pas
de risque de rupture.

Cercle2 et cercle 3 : Si la contrainte horizontale diminues3i ou la contrainte

verticale s1i augmente, le cercle1 ira en augmentant jusqu’à ce qu’il touche

l’enveloppe de rupture. A ce moment les contraintes seront plus fortes que la

résistance au cisaillement du sol, ce qui entraînera une rupture le long d’un plan
de cisaillemnt.

Pour qu’un relachement de la contrainte horizontale contre un ouvrage

de soutènement entraîne la rupture du sol, le mur doit nécessairement
être flexible : c’est ce sur quoi s’appuie la théorie de Rankine.

Cours Mecanique des sol 82 Khaled MEFTAH

 Expansion latérale
t
a s

plan de rupture

a p/4+/2
j/2
0s s s 01 s s

fig 7.2 Etat de poussée

s’A =KA.s’v
KA : Coéfficient de poussée

Compression latérale
t
0s s b b s

01 02 s s

Remarquons que KA.KP = 1

Fig 7.3 Etat de butée

s’P =KP.s’v
KP : Coéfficient de butée

Cours Mecanique des sol 83 Khaled MEFTAH

 2-1 Massif à surface horizontale

b) Cas d’un sol pulvérulent

Pour déterminer le coefficient de poussée, il suffit d’établir une relation entre

s’A ,s’v et l’angle de frottement interne j , on aura :

- en poussée : s’A = KA s’v (7.2)
- en butéé : s’P = KP s’v (7.3)

1- sinj' p j' 1
K = =tg² - =Nf K

avec : A 1+ sinj' 4 2 et P = (7.4)

KA

c) Cas d’un sol cohérent et frottant:

A court terme : (cu¹0 ; ju= 0)

s A=sv -2Cu (7.5)

s P=sv +2Cu (7.6)

Les résultats sont exprimés en terme de contraintes totales.

A long terme : (c’¹0 ; j’¹ 0)

s ' =K s ' - 2C' K A (7.7)

A
A v
s ' =K s ' + 2C' K P (7.8)
P P v

Cours Mecanique des sol 84 Khaled MEFTAH

 2-2 Massif à surface inclinée

s s
t

figure 7.3

La contrainte verticale appliquée à une profondeur z sur une facette parallèle

au talus est :

s v=g.z.cosb (7.9)

Par projection :

s N=sv cosb =g .z .cos²b (7.10)

t = s v sin b = g .z. cos b. sin b (7.11)

Les contraintes de poussée et de butée se calculent de l même manière que

précédemment mais en utilisant les coefficients suivants :

cosb - cos²b - cos²j K (b ) 1 (7.12)

K A(b )
=
=
cosb + cos²b - cos²j P
K A(b )

Cours Mecanique des sol 85 Khaled MEFTAH

 3- Calcul des efforts de poussée et de butée

3-1 Massif à surface horizontale

a) cas d’un sol pulvérulent

En exprimant la contrainte verticale en fonction du poids volumique du sol, on

aura : s’A = KA.g.h qui sera représentée par le diagramme suivant (fig 7.4)

s g

s O s g

fig 7.4 Forces de poussée et de butée

Poussée :
h
P = s' =1 gh²K
A ∫ A 2 A située à h/3 par rapport à O (7.12)
0
Butée :
h
P = s' =1 gh'²K
P ∫ A 2 p située à h’/3 par rapport à O (7.13)
0

Lorsqu’une surcharge uniformément répartie est appliquée à la surface du sol,

à une profondeur z, sv = g.z + q

La contrainte de poussée est : sA = KA.(g.z + q) (7.14)
La contrainte de butée est : sP = KP.(g.z + q) (7.15)

Cours Mecanique des sol 86 Khaled MEFTAH

 Puisque q est constante, sa répartition le long du mur est un rectangle de
surface (q.h) ; la poussée résultante s’applique à h/2 du mur :
PA= KA q h (7.16)

Et par analogie, la butée résultante s’applique à h’/2 du mur :

PP= K p q h (7.17)

Si le sol est saturé , l’action de l’eau est équivalente dans toutes les directions.
La poussée de l’ouvrage de soutènement est indépendante des coefficients de
poussée et de butéé. Elle a une répartition triangulaire comme d ’habitude :
P
ω γωh² ω
= (7.18)

hw : hauteur de l’eau sur l ‘ouvrage

b) Cas d’un sol cohérent et frottant

On a démontré que :

A long terme : s ' = K s ' - 2C ' K

A A v A
j'
K = tg ²
p

Avec -
A 4 2

et s ' = K P s ' + 2C' K P

P v
p j'
K p =tg² +
Avec 4 2

les forces de poussée et de butée s’écrivent :

1
P A= 2gh ²K A- 2c' K A .h (7.19)

1 (7.20)
PP= 2gh '²K p + 2c ' KP .h'

Cours Mecanique des sol 87 Khaled MEFTAH

A court terme :
s'A=g.h-2cu et s 'p =g .h'+2cu et les forces sont
:
P A =1 gh²-2cu.h (7.21)

2
PP =1 gh '²+ 2cu .h' (7.22)

3-2 Massif à surface inclinée d’un angle (b)

fig 7.5 Massif a surface incliné

Si le sol est pulvérulent, à une profondeiur z :
En poussée :

s'A=g.zK A(b )cosb

1
P A= 2g .h²K A(b )cosb (7.23)

Cette force de poussée a deux composantes :

P AH =PAcosb P AV =PAsinb (7.24)

En butée :

s'P b cosb

=g.zK
P (7.25)
1
P P= 2g .h²KP(b )cosb
Cours Mecanique des sol 88 Khaled MEFTAH
 De la même maniére que la force de poussée, la force de butée se

P =P cosb P =P cosb
décompose en deux foces. PH P et AV P
Pour un sol cohérent et frottant, l’effet de la cohésion se retranche de la
poussée et s’ajoute à la butée.
1
P A= gh ²K b - 2c' K A (7.26)
2 A

P =1gh ²K b + 2c' K

P (7.27)
p
P 2
3-3 Massif en plusieurs couches hétérogèn es

Les méthodes de calcul citées précédemment restent valables. Mais, les

contraintes de poussée ou de butée doivent être calculées à l’interface des
couches en tenant compte des paramètres de chacune d’elles. Pour l’exemple
de la figure 7.6 suivante, les contraintes de poussée aux points A,B et C sont :
A

j
sB1 sB2

g
j
sC

figure 7.6

Cours Mecanique des sol 89 Khaled MEFTAH

 p j'
A 1 A1 avec K A =tg² 4 - 2

s =-2c K

sB1=K A g1h1-2c1 K A
1
1

sB2 =K A2g1h1-2c2 K A2
sC =s B2+K A2g2h2
4- La stabilité des murs de soutènement

Deux types de murs de soutènement sont envisageables : les murs poids en

maçonnerie ou en béton et les murs cantilevers en béton armé
Pour qu’un mur de soutènement soit stable, il ne doit ni glisser, ni se
renverser, ni poinçonner le sol de fondation.
4.1 Vérification vis-àvis du glissement

fig 7.7 Forces appliquées à un mur de soutènement

Cours Mecanique des sol 90 Khaled MEFTAH

 La poussée active est la seule force qui tend à faire glisser un mur de
soutènement sur sa base. La résistance au glissement par le frottement entre
le matériau constituant le mur et le sol support (T). On calcule le coefficient de
sécurité Fs :

F = T = c .B + ∑ j
Fv.tg
s ∑F ∑F (7.28)
H H

Pour que le mur soit stable :

- Fs >1.5 si on néglige la butée.
- Fs >2 si on tient compte de la butée

tgj: coefficient de frottement entre le sol de fondation et la base du mur

c: cohésion réduite.
On prendra en général :

j £ £j
3 j
et c : = c ( à long terme) c : =cu ( à court terme).

4-2 Vérification vis-à-vis du renversement

Pour s’assurer qu’un mur ne basculera pas autour du point aval le plus éloigné
sous la semelle(pointO), il faut connaître la valeur des moments résistants
(stabilisateurs) et des moments moteurs (de renversement) par rapport au
point(O). On calcule le coefficient de sécurité Fs :
(7.29)
F s=∑Momentsrésistants
∑momentsmoteurs

Pour que le mur soit stable :

- Fs >1.5 si on néglige la butée.
- Fs >2 si on tient compte de la butée
-

Cours Mecanique des sol 91 Khaled MEFTAH

 4-3 V érification vis-à-vis de la portance du sol de fondation

La résultante des forces doit se situer à l’interieur du noyeau central de façon

que la pression entre la semelle et le sol soit positive en amont ; en aval, cette
pression de contact ne doit pas être supérieure à la capacité portante
admissible du sol. Pour déterminer le lieu de passage de cette résultante, on
doit calculer le moment total autour du point (O) :

Mrésultant =∑M (moteurs)-∑ M(résistants)=Mrn

L’éxcentricité de la force par rapport au point O est :

e = Mrn
o
∑FV
L’éxcentricité de la force par rapport au centre de la semelle est :
B
ec = 2 - eo
On calcule alors les contraintes à l’aval et à l’amontde la semelle (smax) et

(smin)
smax ∑F (1+ 6e ) c £ qadm (7.30)

= V
B B
smin ∑F (1- V
6ec ) ³0 (7.31)

=
B B

5- Les rideaux de palplanches

5-1 Les types constructifs de palplanches

Les rideaux de palplanches métallaques, sont très utilisés dans les ouvrages,
aussi bien provisoires que définitifs. Cependant, l’étude de leur stabilité dépend
de leur mode de fonctionnement. On distingue alors :

Cours Mecanique des sol 92 Khaled MEFTAH

 - Les rideaux simplement encastrés en pied, dont la stabilité est
assurée par la réaction du sol de la partie enterrée, appelée fiche.
C’est le cas essentiellement des batardeaux.
- Les rideaux encastrés en pied et ancrésdont la stabilité, outre la
réaction sur la fiche, est due à un ou plusieurs ancrages enterrés
dans le sol. C’est le cas des murs de quai, etc..
- Les rideaux simplement butés en pieds et ancrés.

5-1 Les rideaux simplement encastrés en pied:

A l’quilibre limite, le
rideau pivote, sans se
déformer, autour du point

figure 7.8 Rideau simplement encastré

Les inconnues sont : la fiche « f » et la contre butée « C ».

Le principe fondamental de la statique, nous fourni deux équations :

{∑ F = 0}et{∑ Mo = 0}
Ainsi le système sera isostatique et résolvable :
On aura la longueur de la palplanche et les éfforts qui lui sont appliqués.
En traçant les diagrammes des moments fléchissants et des efforts tranchants,
on peut
Pour dimensionner la palplanche selon son module de résistance( I / v) en
appliquant l’équation de la résistance des matériaux :

v. Í σ' a
. (7.32)

Cours Mecanique des sol 93 Khaled MEFTAH

 5-2 Les rideaux ancrés simplementbutés en pied
La méthode simplifiée de calcul consiste à ne pas tenir de la contre-butée. Le
problème est alors isostatique et se résout en appliquant le principe
fondamental de la statique :

Tirant

fig7.9 Rideau ancré simplement butée

{∑ F = 0}et{∑ M F / A = 0}

Cours Mecanique des sol 94 Khaled MEFTAH

 QUESTIONS A DEBATTRE

1- A l’aide des cercles de Mohr, expliquer le principe de base des pousseés et

des butées selon la théorie de Rankine.
2- Quelle est l’utilité d’un mur de soutènement ?
3- La stagnation d’eaux a l’amont d’un mur de soutènement est favorable ou
défavorable a sa stabilité ?
4- Quels est le role des barbacanes dans les murs de sotènement ?
5- Expliquer le role des bêches dans les murs de soutènement.

Exercice1
On veut remblayer un mur de soutènement en béton de 5 m de hauteur avec
3
un sable de poids volumique 18kN/m et ayant un angle de frottement interne
de 30°.

0.80
m5

1.50 m

1) Calculer la poussée activesur le mur en utilisant la théorie de Rankine. Le

mur supporte uniquement le remblai.

Cours Mecanique des sol 95 Khaled MEFTAH

 2) Sachant qu’il y a une surcharge de 30kPa sur le terrain situé en amont du
mur. Calculer dans ce cas la poussée active.

3) Vérifier la stabilité du mur au glissement et au renversement, lorsqu’il

supporte le remblai et la surcharge. Le poids volumique du béton est de 25kN/
3
m
4) Vérifier uniquement le renversement si en plus la hauteur d’eau à l’amont
est de 5 m

Exercice 2
Un rideau de palplanches simplement encastré est représenté par le schéma
suivant :

004.

Eau
008.

g=16
kN/m3

j=30°
C=0

004.
008.

g=16 g=16
kN/m3 kN/m3

j=20° j=20°
C'=10 C'=10

1) Tracer les diagrammes des contraintes de poussée et de butée sur le

rideau de palplanches (méthode de Rankine)
2) Calculer les valeurs des forces ainsi que leurs positions par rapport
au point O. En déduire les valeurs des moments.
Cours Mecanique des sol 96 Khaled MEFTAH
 REFERENCES BIBLIOGRAPHIQUES

*Eléments de mécanique des sols : François schlosser

*Aide memoire de mecanique des sols : Publication de l’ENGREF

*Notes de Cours de géptechnique (IUT Paul Sabatier) : Claude Legrand

*Mecanique des sols Tome1 et Tome2 : J.Costet et G.Sanglerat

*Problèmes pratiques de mécanique de mécanique des sols : G.Sanglerat

G.Olilari et B. Cambou.

*Mécanique des sols : V.Robitaille et D.Tremblay

*Fascicule 62 titreV

*DTU 13.12

Cours Mecanique des sol 97 Khaled MEFTAH