COMPRESSIBILITE DES SOLS

1. - CALCUL DES CONTRAINTES AU SEIN D'UN MASSIF

1.1. - Cas des charges ponctuelles

Pour déterminer les contraintes au sein du massif, on

utilise les résultats de Boussinesq sur la déformation
élastique d'un milieu non pesant, isotrope, homogène, limité
à sa partie supérieure par un plan horizontal illimité et
soumis à l'action d'une force verticale isolée P Boussinesq
a démontré que, dans ce cas, la contrainte qui s'exerce sur
une facette horizontale centrée en M passe par le point
d'application 0 de la force P. La composante normale de

cette contrainte a pour valeur :

En posant  = OM et alors

On remarque donc que les courbes d'égale pression verticale ont pour équation . Ces

courbes ont l'allure de bulles. D'autre part, ces expressions sont indépendantes des constantes
élastiques du milieu : module d'Young et coefficient de Poisson. A une certaine profondeur, la
distribution de la composante normale  sur un plan horizontal se fait suivant une courbe de Gauss,
alors qu'auparavant on estimant cette répartition comme étant uniforme sur un cône d'angle au
sommet 90°.

1.2. - Cas des charges uniformément réparties

Dans ce cas, il faut utiliser les procédés d'intégration.

Chaque élément d'aire dA porte une charge équivalente à
une force concentrée q.dA (q : densité de chargement) et
engendre en M une pression verticale d.

alors

Nous verrons dans quelques cas particuliers, la valeur de cette contrainte .

2.- ETUDE DU TASSEMENT

2.1.- Généralités et définitions
Lorsqu'un massif de sol est chargé, les contraintes qui se développent en son sein provoquent des
déformations. La somme de ces déformations permet d'avoir le déplacement des limites du massif. La
plupart du temps, les charges appliquées sont verticales et la surface du sol proche de l'horizontale.
Les déplacements verticaux à la surface du sol sont prépondérants et ont des répercussions possibles
sur les structures reposant sur le sol. Ces déplacements verticaux sont appelés tassements. Si on
connaît la loi de comportement du sol, on peut connaître les contraintes et les déformations en tout
point. On peut en déduire les tassements. Mais cette loi est inconnue. On peut sans trop d'erreur,
utiliser la théorie de l'élasticité, comme nous venons de le voir pour le calcul des contraintes
verticales, puisqu'elles ne dépendent pas des caractéristiques élastiques ; module d'Young et
coefficient de poisson. Mais cette théorie n'est pas applicable pour le calcul des tassements. On peut
alors utiliser plusieurs méthodes :

a) Détermination en laboratoire du tassement relatif d'un échantillon "intact" du sol : (h : hauteur,

∆h : tassement). Le tassement relatif du sol étant supposé identique à celui de l'échantillon.

b) Détermination d'un module de déformation par un essai de chargement en place : essai
pressiométrique.
c) Applications de formule empiriques, à partir de certains essais en place : essai de chargement à la
plaque, pénétromètre.

2.2.- Généralités sur le processus du tassement dans le cas des sols grenus et des sols fins :
a) Sol grenu : Prenons comme exemple un sable propre perméable. S'il est sec ou partiellement
saturé, les contraintes appliquées provoqueront un réarrangement des grains tendant à diminuer
l'indice des vides. Il s'ensuivra un tassement qui se produira immédiatement après l'application de la
charge. S'il est saturé, la pression de l'eau interstitielle créée par les charges appliquées se dissipe
immédiatement et le même processus s'établit. Notons que le phénomène ne fait pas intervenir le
temps.
b) Sol fin : Choisissons pour exemple, une argile imperméable saturé. D'après la théorie de la
consolidation (voir les propriétés hydrauliques), les contraintes appliquées sont entièrement encaissées
par l'eau interstitielle. Ce n'est que lentement que ces contraintes se transmettent au squelette
solide, produisant petit à petit un tassement. Dans ce cas seulement, le temps intervient.
On voit donc que tout système de contraintes appliqué en un point d'un massif de sol se décompose :
 en un tenseur sphérique : les tassements correspondant seront dus au phénomène de
consolidation ;
 en un déviateur : les tassements correspondant seront dus aux contraintes de cisaillement.
Dans le premier cas, on utilisera pour l'étude des tassements la méthodes oedométrique (ex. remblai
large sur une couche d'argile peu épaisse).

Dans le 2nd Cas, on utilisera la méthode pressiométrique (ex : semelle isolée carrée sur limon sableux,
avec ou sans nappe).

3. - ESSAI DE COMPRESSIBILLITE A L'OEDOMETRE

3.1.- Généralités
Pour étudier le tassement d'une couche d'argile théoriquement infinie en (donc pratiquement de
dimensions suffisantes) et chargée quasi-uniformément par une contrainte verticale, on commence par
extraire de la couche, par des moyens de sondages appropriés, un échantillon aussi "intact" que
possible (c'est-à-dire ayant conservé ses propriétés mécaniques). On utilise alors un appareil
comprenant un moule cylindrique (oedomètre de Terzaghi) de 7 cm de Ø intérieur et généralement
l'échantillon a une hauteur de 12 ou 24 mm. Deux pierres poreuses permettent de drainer et évacuer
l'eau qui est expulsée au cours de l'essai, car pour les argiles on opère le plus souvent sur des
échantillons saturés afin d'éliminer l'influence des forces capillaires. On applique sur le piston une
charge afin d'éliminer l'influence des forces capillaires. On applique sur le piston une charge
croissante et on observe les déformations.
3.2.- Courbes tassement - temps pour une charge donnée
L'échantillon étant mis en place, on applique sur le piston, une première charge jusqu'à consolidation
complète (généralement 24h), puis ensuite  =  + ∆ . En général les laboratoires ont l'habitude
2 1 1
de doubler la charge toutes les 2'h même si le tassement n'est pas stabilisé. Pour chaque charge, on
trace la courbe tassement-temps, le tassement étant mesuré par des comparateurs au 1/100 mm.
L'exploitation des résultats conduit à déterminer le coefficient de consolidation (Cv) par l'une des
deux méthodes suivants.

3.2.1.- Méthode de Casagrande

Si on porte les tassements

(indications du comparateur) en
fonction du logarithme du
temps, la courbe présente 2
parties linéaires dont
l'intersection marque
pratiquement la fin de la
consolidation (t100 : ∆u = 0). Pour
t > t100 le sol continue à tasser
légèrement bien que ∆u = 0 :
c'est la compression secondaire.

Il faut déterminer le début de la consolidation d 0 qui n'apparait pas directement sur la courbe puisque
la mise en place de l'appareil entraîne un tassement immédiat d 1. Or on admet que pour un tassement
inférieur à 50%, il varie linéairement avec . On peut utiliser les propriétés de la parabole. Si on
prend un temps t assez petit, le tassement entre 0 et t est le même qu'entre t et 4t. Sur la figure, BA
représente le tassement entre t et 4t. Si on veut déterminer les variations du coefficient de
consolidation Cv avec la charge, il faut tenir compte des observations suivantes :
 Pendant le premier palier de chargement  , la pression moyenne effective pendant la
1

consolidation est et le coefficient de consolidation est Cv1.

 Pendant le deuxième palier, nous aurons  =  + ∆ et la pression moyenne devient ' =  +

1 1 1 2 1

et le coefficient de consolidation est Cv2.

Si on trace la courbe de
variation de Cv en fonction de la
pression moyenne effective,
nous obtenons un arc d'aspect
parabolique : (ceci afin de
pouvoir choisir par la suite la
valeur de Cv qui correspond au
système de charges appliqué à la
couche de sol compressible dont
on veut étudier le tassement).

3.3.- Courbe tassement final/charge

 Sur un diagramme semi-
logarithmique, on porte l'indice
des vides de l'échantillon après
le tassement final obtenu pour
chaque charge (∆hf = df - d0)
en fonction de la contrainte
appliquée qui est effective en
fin de consolidation (u = 0).
Notons que ce tassement final
comporte le tassement
immédiat, le tassement de
consolidation et une partie du
tassement secondaire.

La courbe obtenue comporte une première partie droite AB légèrement au-dessous de e 0, l'indice des
vides initial, et une deuxième partie droite BC. La pression de préconsolidation ' correspond
sensiblement à l'intersection dess 2 droites. Elle représente la plus grande contrainte effective qu'a
pu subir le sol au cours de son histoire. Ce qui signifie que puisque le sol tasse pour ' < 'c, celui-ci
garde "en mémoire" une partie de son histoire. La pente BC qu'on appelle indice de compression est

donnée par :

Dans cette relation e et ' caractérisent l'état initial. Si on applique une surcharge ∆ ', l'échantillon
est alors soumis à une pression '1 = ' + ∆' et l'indice des vides e1 peut s'écrire :

e1 = e – Cc log

Cette relation permet de calculer le module oedométrique E' donné par :

On porte alors C tel que BC = BA puis de C, on mène une parallèle à l'axe des temps. On obtient d 0

correspondant à t = 0. Le tassement de consolidation est donc (d 100 - d0). . On en

déduit que qui sera le temps nécessaire pour avoir 50% de consolidation.

Dans cette formule, H représente la demi-hauteur de l'échantillon et Tv est un facteur temps fonction
du degré de consolidation. Tv (50) = 0,197.

3.2.2.- Méthode de Taylor

Si on porte les tassements en fonction de la
racine carrée du temps, la première partie
de la courbe, compte tenu de la loi ∆h = A.
, devrait être une droite. Or on obtient
cette droite qu'au bout du temps t 1. Celle-
ci coupe l'axe en dc. La différence dc - d0
entre les lectures faites entre le temps t =
0 et le début de la consolidation, constitue
le tassement immédiat, résultat de la
dissolution éventuelle des gaz occlus dans
l'eau mise en pression. Le zéro du
tassement de consolidation correspond donc
en fait à dc. Cherchons le point d90
correspondant à 90% de consolidation ;
c'est le point C de la courbe.

Soit EB la droite d'équation ∆h = A. . Terzaghi a montré que si U représente le degré de

consolidation, on peut déterminer le facteur Tv :

Si U ≤ 50% Tv = Si U > 50% Tv = 1,781 - 0,933 log (100 - U)

Sur la figure la distance AB est égale à . Or la valeur de A est . Si

on joint EC, on aura AC = k.AB. On peut remplacer AC et AB par leur valeur connaissant la valeur de
Tv(90) = 0,848. Ceci permet de calculer le coefficient k ; on en déduit que k = 1,15. Il suffit donc à
partir de E de mener une droite dont la pente est 1,15 fois celle de EB. Elle coupe la courbe
expérimentale en C qui correspond à t90.

CV =

* Pour les limons Cv est de l'ordre de 10- 2 cm2/s

* Pour les argiles limoneuses 10- 3 cm2/s
* Pour les tourbes, vases, argiles 10- 4 à 10- 5 cm2/s.
On en supposant que les grains sont incompressibles et que le tassement résulte uniquement d'un
réarrangement des grains, on écrira que le volume des grains de l'échantillon est resté constant ; ce

volume par unité de surface a pour expression : . Par suite : soit, après

simplification

et ainsi

Si ∆' est petit devant ', cette expression se simplifie : ce qui correspond à un

module instantané.
Détermination de la pression de pré-consolidation 'c par la méthode de Casagrande
Ainsi que nous l'avons vu précédemment, 'c se trouve sensiblement à l'intersection des 2 parties
droites. Casagrande propose une autre méthode : Par le point de courbure maximum B, on mène la
tangente et une horizontale, puis la bissectrice de l'angle formé par ces 2 droites. L'intersection de
cette 2ème partie droite de la courbe oedométrique avec cette bissectrice donne 'c. On voit qu'il y a
une incertitude sur 'c ; c'est la raison pour laquelle l'ordre de grandeur de 'c est donné par
l'intersection des 2 parties droites.

4. - COMPORTEMENT DU SOL EN PLACE

4.1.- Généralités : couches normalement consolidée
Considérons une épaisse couche d'argile homogène, c'est-à-dire dont la nature est constante avec la
profondeur. Les variables seront le poids spécifique s, la teneur en eau w, et l'indice des vides e qui
sont fonctions de la profondeur. Comme on va le voir, les couches de surface n'ont pas le même
comportement qu'en profondeur. Nous ferons cette étude à partir d'une profondeur z0 (de l'ordre de
3 à 4 m). Supposons que la nappe soit en surface et que la couche n'ait jamais supporté de charges
autres que celles des couches sus-jacentes. On supposera aussi que la pression interstitielle u est celle
due à la nappe. A la profondeur z la contrainte verticale effective est 'oz = ( - w).z. L'indice des
vides est e. La courbe 'oz = f(ez) est la courbe de compressibilité du sol en place. On l'appelle
courbe vierge.

Si on extrait un échantillon au niveau z, et qu'on le

soumet à un essai oedométrique, il gonflera légèrement
lors de l'extraction qui correspond à un déchargement
(courbe ADB) et prendra un indice des vices e o ; puis il
tassera légèrement sous l'effet des charges jusqu'à
'oz (courbe BCA). Au delà de 'oz la loi sera la même
que pour les couches situées à des profondeurs
supérieures à z.

On retrouve donc l'allure du diagramme oedométrique e/log '. L'indice de compression Cc correspond
donc à la pente de la courbe vierge. On dit qu'une telle couche d'argile est normalement consolidée.

4.2.- Calcul du tassement total

Si on charge en surface la couche précédente par une pression verticale q, donnant au niveau z une
pression supplémentaire ∆, la couche d'épaisseur dz, va tasser. Si on considère la courbe
oedométrique d'un échantillon au niveau z, puisque la couche est normalement consolidée : 'c = 'o =

'.z. L'augmentation de contrainte ∆ diminue l'indice des vides de . Or pour

un échantillon d'épaisseur h le tassement ∆h est : ou encore : .

 Pour la couche de sol d'épaisseur dz le
tassement élémentaire dh sera donné par :

Donc pour la couche le sol situé entre zo et zf

le tassement sera :

Dans la pratique, on distingue deux cas :

a) couche homogène : Cc est constant. On découpe alors la couche en sous couches de 3 à 4 m pour
lesquelles on suppose eo, 'o, ∆' constants et ayant pour valeurs celles déterminées au centre de la
sous-couche. On calcule alors les ∆h correspondants et on fait la somme.
b) couche hétérogène : On découpe en couches de nature homogène et on détermine Cc sur chacune
d'elles. On pratique ensuite comme pour des couches homogènes.

4.3. - Cas des couches surconsolidées : ' c > ' o

Une couche est sur-consolidée, si la pression de pré-consolidation trouvée à l'oedomètre ( 'c) est
supérieure à 'o, pression effective actuelle des terres au niveau où est prélevé l'échantillon. Cela
peut être dû à une érosion des couches supérieures ou à une modification des liaisons inter-granulaires
(dessèchement périodique des couches de surface que nous avons laissées de côté précédemment).
Pour une telle couche d'épaisseur h, on voit que la pression passe de 'o à 'o + ∆'. On suppose que le
tassement est négligeable au-dessous de ' c .

4.4. - Cas des couches sous-consolidées : ' c < ' o

Certains dépôts récents peuvent ne pas être totalement consolidés sous leur propre poids. A un certain
niveau 'c < 'o. Dans ce cas le tassement sera :

Donc on voit que quelque soit

le cas le tassement est
toujours de la forme :

 'c = 'o : sol normalement consolidé

  'c < 'o sol sur-consolidé
 'c > 'o sol sous-consolidé

5. - CALCUL DE L'ACCROISSEMENT DE CONTRAINTE VERTICALE RESULTANT DE

L'APPLICATION D'UNE CHARGE EN SURFACE
Reprenons les résultats du paragraphe 1 qui nous permettent de déterminer dans quelques cas simples,
l'accroissement de contrainte ∆. Pour cela on utilisera :
 la géométrie de la surcharge, uniforme ou non,
 la géométrie de la surface chargée.
La fondation est toujours supposée souple et les résultats se présentent sous forme d'abaques.

5.1. – Application pratique de la théorie de Boussinesq (Théorie de l’élasticité)

On a recourt à la théorie de l’élasticité ( voir paragraphe ! 1.1.) pour calculer les contraintes dans les
sols. Boussinesq a developpé en 1885 un certain nombre de résultats qui permettent de déterminer les
contraintes induites par une charge ponctuelle dans un espace homogène, isotrope et dans lequel la
théorie de l’élasticité linéaire reste valable. La contrainte verticale est égale à :

Q : Charge ponctuelle,

z : Profondeur du point d’application de z ;

r : Distance horizontale entre le point d’application de la contrainte et l’azimut
d’application de Q.

On peut aussi écrire z de la forme suivante : (avec NB qui est un facteur de forme qui

regroupe les termes constants de z en fonction de r/z. On retrouve aussi les contraintes radiales,
tangentielles et de cisaillement à partir de ces travaux de Boussinesq.

Les termes constants sont

exprimés à la figure ci-
contre.

L’intégration des charges ponctuelles sur un axe permet d’en déduire la contrainte due à une charge

linéique. Ainsi, z s’exprime par : où : P : Charge linéique

Newmark (1935) a trouvé les solutions de pour la contrainte verticale sous le coin d’une

charge rectangulaire chargée uniformément :

 1 2mn(m 2  n 2  1) (m 2  n 2  2) 2mn(m 2  n 2  1)1/ 2
 z  qo [ 2   arctan ]
4 m  n 2  1  m 2 n 2 (m 2  n 2  1) m2  n2  1  m2n2
avec :
qo : Contrainte superficielle,
m = B/z (B : largeur de la surface uniformément chargée) ;
n = L/z (L : longueur de la surface uniformément chargée).

M et n sont inter-changéables. L’équation peut être écrite comme suit . I est le

facteur d’influence, il est fonction de m et n.

5.1.- Surface circulaire chargée uniformément

La surface de rayon r est chargée uniformément par une contrainte verticale q. L'accroissement de
contrainte verticale, sur l'axe à la profondeur z est donnée par :

5.2.- Surface rectangulaire chargée uniformément

L'accroissement de contrainte ∆ s'exerçant sur un plan horizontal en un point M à la profondeur z,
z
situé à la verticale d'un coin d'un rectangle chargé uniformément est donné par :

k est donné dans une abaque en fonction de 2 paramètres : et . Pour avoir la contrainte à

la verticale d'un point M quelconque, on applique le principe de superposition.

2 cas :

kM = KI + kII + kIII + kIV

kM = kACMH + kHMDG - kBCME - kEMDF

5.3. - Bande de largeur infinie supportant une charge trapézoidale

Ce cas se rapporte à celui des remblais de grande longueur. L'accroissement de contrainte est fonction

d'un coefficient d'influence I lui-même fonction de et :

 z = q.I
Les valeurs de I sont données par
l’abaque suivante

Dans le calcul des contraintes sous les remblais de grande longueur, la charge appliquée sur le sol en un
point est supposée égale à  • H ( : poids spécifique du sol constituant le remblai, H : hauteur du
remblai au-dessus du point considéré) ; on utilise ensuite le principe de superposition. Deux cas peuvent
se présenter :

a) Point sous la surface chargée

z = q.(I1 + I2)

Les valeurs de I1 avec a1 et b1
Les valeurs de I2 avec a2 et b2

b) Point en dehors de la surface chargée

z = q.(I1 - I2)