J.P.

GIROUD

Ingnieur E . C P .
Licenci s - S c i e n c e s
Docteur de Spcialit

H. W A T I S S E E
Ingnieur E . C P .
Docteur de Spcialit

A. RABATEL
Technicien

tassements et contraintes
dans une couche de sol lastique
supportant une charge
uniformment rpartie

Laboratoire de mcanique des sols

Universit de Grenoble

prsentation
F. S C H L O S S E R
Ingnieur des Ponts et C h a u s s e s
Chef de la section " mcanique des s o l s "
du Laboratoire Central

L'valuation
des tassements des ouvrages fonds
superficiellement
est un problme
important mais dlicat.
La difficult
rsulte,
pour une grande part, de la complexit
du comportement
des sols, comportement qu'il n'est pas possible de reprsenter
par des lois simples.
Les mthodes
d'valuation
sont donc nombreuses
mais toutes imparfaites, depuis
les
mthodes
utilisant la thorie
de l'lasticit
linaire
jusqu'aux mthodes
plus ou moins
empiriques
extrapolant au cas rel des rsultats
d'essais de chargement
simplifis
: essais
domtriques,
essais sur modles
ou essais de plaques, essais en semi-grandeur,
etc. Dans la mesure o il
n'existe pas de mthode
gnrale,
il s'agit dans chaque cas de savoir utiliser non pas celles qui
donneront coup sr les meilleurs rsultats,
mais celles qui sont les mieux adaptes
la fols au
cot de l'ouvrage et la prcision
que l'on demande dans la connaissance
des tassements. Ainsi,
il sera possible pour certains ouvrages de se contenter d'un rsultat
100 % prs, lorsque le tassement sera de l'ordre de quelques centimtres,
ce qui ne sera pas le cas pour un remblai construit sur un sol trs compressible,
o le tassement atteindra le mtre ou plus.
L'article de MM. Giroud, Watissee et Rabatel donne la rsolution
complte
d'un cas de
chargement souvent rencontr
dans la pratique : celui d'une couche de sol reposant sur un substratum indformable
et supportant une charge uniformment
rpartie,
un remblai, par exemple. Le
sol est suppos
obir une loi de comportement lastique
linaire.
Les rsultats
concernent la
fois les contraintes et les tassements et sont donns
sous forme
d'abaques.
S'il est bien certain que le comportement d'un sol n'est jamais lastique
et linaire,
article prsente
pour le mcanicien
des sols beaucoup d'intrt
et plusieurs titres :

cet

- la connaissance
de l'tat des contraintes dans le sol de fondation sous un remblai permet, en
effet, de s'orienter vers des mthodes
plus labores
que la mthode
domtrlque
et utilisant par
exemple l'essai triaxial. Ce processus
est justifi
par le fait que les contraintes sont moins sensibles que les tassements aux carts
qui peuvent exister entre le sol rel et celui de la
thorie
lastique
;

97

Bull.

Liaison Labo.

Routiers P. et C h . no 48 - N o v . 1970 - Rf. 9 2 5

b) D f i n i t i o n de la charge
Le calcul est fait pour une couche compressible supportant sa surface une charge normale uniforme de
longueur infinie et de largeur 2a ( f l g . 1). Bien qu'en ralit une longueur ne puisse tre infinie, on peut
reprsenter ainsi les charges dont le rapport longueur/largeur est grand pour calculer le tassement ou les
contraintes en un point loign de l'extrmit de la charge : il en rsulte des calculs plus simples
(dformation
plane) et des rsultats plus commodes. Pour avoir un ordre de grandeur de l'approximation
ainsi faite, distinguons deux cas :
pour les contraintes, nous avons calcul [3] [4] l'erreur commise dans le c a s d'un massif semi-infini
de sol homogne (autrement dit : couche compressible d'paisseur infinie). S i l'paisseur de la couche
est finie, l'erreur est encore plus faible car la prsence du substratum diminue l'influence des charges loignes. Il en rsulte que, pour les charges telles que les digues, remblais routiers et semelles filantes,
l'erreur faite sur les contraintes en supposant la longueur infinie est trs faible ;
pour les tassements, le problme est plus dlicat. En effet, il est bien connu que, sur un massif semiinfini, le tassement d'une charge de longueur infinie est infini. Le raisonnement fait pour les contraintes
n'est donc plus applicable. O n peut seulement affirmer que l'erreur est d'autant plus faible que le rapport
longueur/largeur est plus grand et que l'paisseur de la couche est plus faible ; et if semble logique de
prsumer que cette erreur est plus grande pour les tassements que pour les contraintes. Une tude en
cours G r e n o b l e [5] permettra plus tard de prciser ce point.

3 - N O T A T I O N S ET C O N V E N T I O N S
L e s notations sont rpertories comme suit :
- a : demi-largeur de la charge uniforme,
- E : module de Y o u n g du sol de la couche compressible,
- E' : module domtrique du sol de la couche compressible,
- H : paisseur de la couche de sol compressible,
- (I, J) : indices de ligne et de colonne dans la double grille,
- k m , k m , k H : coefficients sans dimensions pour le calcul de o- , T
2

z x

et o- ,
x

- p : pression uniforme exerce la surface du massif,

- P

p' :

P PP
h

coefficients sans dimensions pour le calcul du tassement,

'

explicite approche de p

valable si B est grand,

- : forces reprsentant la charge p, dans le schma de la double grille,

- t : variable

d'intgration,

- w, u : composantes du vecteur dplacement,

- w', u' : composantes du vecteur dplacement en axes inclins 4 5 ,
99

- x, z : coordonnes du point o l'on calcule les contraintes et le dplacement,

-/?,/?'

: paramtres sans dimensions lis l'paisseui de la couche,

- A : ct des mailles de la double grille,

-

v : coefficient de P o i s s o n du sol de la couche compressible,

- : paramtre de sur-relaxation.
L e s conventions suivantes sont appliques :
- le tassement, w, est positif v e r s le bas comme l'axe Oz,
ff

- les contraintes normales (

,
, a ) sont positives si ce sont des c o m p r e s s i o n s ,
- une contrainte tangentielle, exerce
sur une coupe de normale intrieure parallle un axe de
coordonnes et de mme sens, est positive si elle est de mme sens que l'autre axe de coordonnes (fig- 2).
x

% >o
Fig. 2 - Conventions de signes
pour les contraintes.

4 - EQUATIONS DU

PROBLEME

Les quations qui permettent

de rsoudre

un problme

- les deux quations de l'quilibre indfini indpendantes

o-

2 X

ox

dx

d'lasticit

en dformation

- les trois quations de H o o k e qui indiquent que le sol est lastique

1+v
E

1
+v r
1-f-v
(2)

Ex

1(1

v) a *

VOz

1+v
T

2 X

Tz:

les relations de dfinition des dformations

100

/dw

o u\

\x

dz

plane

des proprits du matriau :

ou
ox

sont

C e s trois composantes de la dformation, drivant d'un vecteur dplacement (w, u), ne sont pas indpendantes ; on en dduit par drivation la relation de compatibilit qui existe entre elles :
0
!

(4)

2 * f

0 E*
!

2X

x*

zbx

Cette relation traduit la continuit de la dformation du milieu

La solution d e s systmes prcdents doit satisfaire a u x conditions

l'expression est :
Tzx=

(5)

z= 0 {

=
T z

(6)

z= H :

aux limites

indiques

( 2a) et

dont

t p si | x [ < a
f 0 si | x | > a

w= u= 0

5 - PLAN DE L'ARTICLE
N o u s avons utilis dans cette tude deux mthodes de calcul tout fait diffrentes :
- la premire, faisant appel aux sries d e Fourier, n'a t utilise que pour les tassements. Elle a permis
d'explorer tout le champ d e s valeurs de H (paisseur de la couche) et de v (coefficient de P o i s s o n du sol) ;
- la seconde, entirement numrique, a fourni simultanment les tassements et les contraintes, mais les
temps de calcul tant trs longs, nous avons d limiter le champ d e s valeurs de H et de v examines. D e
plus, cette dernire mthode ne nous a pas permis d'tudier le c a s du s o l volumtriquement incompressible ( v = 0,5).

N o u s avons donc partag cet article en deux :

- une premire partie intitule Calcul des tassements o sera brivement expose la mthode d e s sries
de Fourier ;
- une deuxime partie consacre au Calcul des contraintes o sera mentionne la mthode numrique.

C A L C U L DES TASSEMENTS

I - DETERMINATION THEORIQUE DES TASSEMENTS

a) R s o l u t i o n du p r o b l m e par les s r i e s de Fourier
II est classique de montrer, en combinant les quations (1), (2) et (4), que l'on peut ramener un problme
de dformation plane lastique une seule quation biharmonique :
(7)

V'> = 0

avec <P (z, x) : fonction de contrainte d'Airy telle que :

s

o <>
(8)

=
x

d <>

t =

<>

a* =
z o i

dz

Il suffit donc de trouver une fonction <t> solution de (7) et satisfaisant aux conditions aux limites (5) et (6).
loi

La charge tant symtrique et finie, on peut choisir $> de la forme :

* - [-5-

o h

(-fr

') + ~jr " (~jr ' )

h

< * ( ' ) + |(-fr')]<== (TT ' )

Les constantes A , B, C et D sont dtermines d'aprs les conditions aux limites (5) et (6) qui s'expriment
respectivement :
!

>
(10)

z =

z d x

<>
x

\ \
(11)

z = H

p si | x | < a

si |

0 v)

I >a

Jo L

) Jo

dz

ftx'

Le calcul de A , B, C et D a t fait par Holl [6] en reprsentant la charge par une intgrale de Fourier :
2 i*
fx
P (x) = \
cos
^ Jo

(*

dt \

00

(12)

p (XH) c o s

Jo

(xt) dx

La fonction q> tant dtermine, l e s contraintes s'en dduisent par (8), les dformations par (2) et les dplacements par intgration de (3). En particulier, pour le tassement de la surface du s o l (z = 0), Holl a obtenu :
2

(13)

w =

4 (1 v )
f"

^ p H \

(3 4 v) cht sht t

Jo 1(3 4 v) ch t
2

itE

(1 2 v ) + t ] t
2

/ tx \
/ ta \
cos (
) sin ( - = - ) dt
V

Nous nous limiterons ici au calcul du tassement en surface l'aide de la formule (13), mais on voit que
cette mthode permet d'obtenir contraintes, dformations et dplacements en tout point de la couche de
sol c o m p r e s s i b l e .
b) Calcul des valeurs n u m r i q u e s du tassement
L'intgrale (13) fournit le tassement de tous les points de la surface du massif. Cependant, il suffit de faire
le calcul pour le point situ au bord de la charge, le principe de superposition permettant d'en dduire le
tassement des autres points. Le tassement au bord de la charge (z = 0, x = a) s e calcule l'aide de
l'intgrale suivante dduite de (13) :
(14)

w = -A-!?-

(15)

p =

avec
2

v )

p\ * (v,

P\

(3
(16)

q? (v. S, t) =
v

102

S, t) dt

* (v, (3.

t)

4 v) cht sht t

(3 4 v) ch t
2

+ (1 2 v ) + t ] t

t
sin
B

en posant

L'intgrale (15)
ne peut pas tre exprime explicitement : il faut donc la calculer numriquement pour
toutes les valeurs intressantes de v et de /3. L'aire de la fonction <J> (t) (fig. 3) est une srie alterne
dcroissante ; la valeur de l'intgrale est donc finie :
(18)

A, A

0 (v,/3,

A* A

.. ..(

1)n+1 An

....

t)

Calcul

Calcul
explicite
approch

numrique

Fig 3 - Allure de la courbe <t> ( v , /3, i)

pour v et P donns.

Le nombre de boucles de la courbe tant trop grand, un planimtrage manuel tait pratiquement impossible.
Nous avons donc dtermin ces aires la calculatrice lectronique avec un pas dt = 0,01 en nous arrtant
lorsque A
est infrieur 0,001 sans toutefois dpasser t = 10. En effet, pour t > 10, on peut crire,
avec une erreur infrieure 10~ :
n

(19)

2 (1

v )

p =

l
J0

(3
| (3

4v)

4 v) cht sht t

Chr't +

(1

2 v)

t] t
2

t
sin
B

dt

f= 1
t
T
sin dt
+\

.lin t

Notons que ce dernier terme ne dpend pas de v.

On peut alors crire

2 (1 (20)

Pu

(3 4 v ) cht sht t

v )
2

B
H

[(3

4 v) ch t

(1 2 v )

t]

sin
10

Nous avons donc la somme de trois termes : le premier a t calcul numriquement comme il est dit
plus haut, le second est explicite et la valeur du troisime ( cosinus intgral ) est donne dans toutes les
tables de fonctions classiques.
N o u s avons fait ce calcul pour un grand nombre de valeurs de v et de B et les rsultats sont donns dans
le tableau 1 et la figure k.
103

Coefficient

de

POISSON

Coefficient

0,2

0,3

0,4

de POISSON

f>

0,1

0,2

0,3

0,4

0,5

1.7

0,777

0,753

0,706

0,626

0,501

0,317

1,8

0,809

0,785

0,738

0,655

0,527

0,339

1.9

0,840

0,817

0 , 7 68

0,682

0,551

0,360

/$

0
0,05

0
0,025

0,024

0,023

0,019

0,012

0,1

0,05

0,049

0,045

0,037

0,023

o 000
o, 000

0,15

0,075

0,073

0,068

0,056

0,035

000

2,0

0,870

0,848

0,796

0,709

0,575

0,381

0,2

0,100

0,098

0,090

0,074

0,047

001

2,5

1,002

0,978

0,925

0,826

0,682

0,468

0,125

0,125

0,122

0,113

0,093

0,059

004

3,0

1,112

1,085

1,029

0,926

0,773

0,548

0,3

0,150

0,147

0,135

0,112

0,072

007

3,5

1,208

1,180

1,115

1,012

0,852

0,619

0,35

0,175

0,171

0,158

0,131

0,085

011

4,0

1,289

1,263

1,197

1,086

0,920

0,683

0,4

0,201

0,196

0,181

0,150

0,099

017

5,0

1,429

1,400

1,330

1,212

1,036

0,786

0,45

0,227

0,

0,204

0,169

0,113

023

6,0

1,545

1,518

1,445

1,318

1,133

0,870

0,5

0,253

0,245

0,227

0,189

0,127

0 031

7,0

1,641

1,615

1,540

1,405

1,213

0,943

0,55

0,279

0,269

0,250

0,209

0,143

039

8,0

1,722

1,699

1,621

1,479

1,281

1,007

0,6

0,305

0,294

0,272

0,229

0,158

049

9,0

1,800

1,770

1,689

1,549

1,347

1,062

0,65

0,331

0,318

0,295

0,249

0,174

059

10,0

1,866

1,834

1,750

1,609

1,402

1,112

0,7

0,357

0,342

0,318

0,269

0,191

070

11,0

1,926

1,890

1,805

1,663

1,452

1,157

0,75

0,382

0,367

0,341

0,290

0,207

1,978

1,942

1,857

1,712

1,500

1,198

0,407

0,391

0,364

0,310

0,224

081
093

12,0

0,8

13,0

2,026

1,991

1,906

1,759

1,544

1,236

0,85

0,432

0,415

0,386

0,330

0,241

106

14,0

2,072

2,038

1,952

1,802

1,585

1,271

0,9

0,456

0,439

0,408

0,350

0,258

118

15,0

2,115

2,081

1,995

1,843

1,623

1,304

0,95

0,480

0,464

0,430

0,370

0,275

131

20,0

2,307

2,269

2,173

2,010

1,772

1,442

1,0

0,503

0,452

0,390

0,291

o
0

144

25,0

2,449

2,410

2,309

2,139

1,891

1,548

1.1

0,548

0,530

0,496

0,428

0,324

170

30,0

2,565

2,524

2,420

2,244

1,988

1,635

1.2

0,591

0,570

0,535

0,464

0,356

196

40,0

2,748

2,705

2,596

2,411

2,142

1,772

1,3

0,632

0,611

0,574

0,499

0,387

222

50,0

2,890

2,847

2,732

2,540

2,261

1,879

0,1
0

220

0,488

0.5
0

1,4

0,671

0,648

0,610

0,533

0,417

246

100,0

3,331

3,284

3,156

2,942

2,632

2,210

1.5

0,708

0,685

0,643

0;566

0,446

270

1000,0

4,797

4,735

4,563

4,276

3,863

3,309

1,6

0,743

0,720

0,676

0,597

0,474

0 294

00

CO

00

OS

T A B L E A U 1 - V a l e u r s du coefficient p

00

00

CO

pour le calcul du t a s s e m e n t d'aprs les formules du tableau 2. C e s mmes valeurs

sont reprsentes graphiquement sur la figure 4.
H

c) Tassement en un point quelconque de la surface du sol

Nous venons d'indiquer comment a t obtenu le tassement du bord
peut en dduire celui d'un point quelconque.

de la charge. V o y o n s comment on

Pour cela, considrons le point M de la figure 5. O n voit que s o n tassement s'obtient partir du tassement
du bord par superposition : ceci nous conduit la formule (24) du tableau 2. S i le point M tait l'extrieur de la charge, l'application du principe de superposition s e traduirait par une soustraction au lieu d'une
addition mais conduirait la mme formule (24).

2 - RESULTATS NUMERIQUES
D'aprs l'tude ( 1), le tassement s'exprime par les formules du tableau 2. Notons que, dans ce tableau,
les expressions (22) et (23) ne sont que d e s c a s particuliers de l'expression gnrale (24).
A partir d e s formules du tableau 2, on en dduit celles du tableau
(21)

p'

= B' p

avec : B =

= _

2a

3 en posant :

a+x

a-x

F i g. 5 - Tassement d'un p o i n t quelconque de l a surface d u s o l dduit

d u tassement au b o r d par application d u p r i n c i p e de superposition.

TABLEAU 2

Avec P =

Le tassement est exprim par :

bord

centre

2 an
(22) w =

(23) w =

- y

H
p

2 a

2 an

Calcul

du tassement par les formules

du premier type.

O n peut

aussi

utiliser celles du s e c o n d type

point
(24) w =
d'abscisse x

1 [(a + x) p

H |

+ ( a x) p

|a
H ]

x\

H -

|a -

3). L e s valeurs numriques

cient

x|

bien

(tableau

du coeffi-

sans dimensions, p , sont lues,

en fonction

de fi, dans

le tableau 1

ou la figure 4.

105

Les formules du tableau 3 sont plus c o m m o d e s que c e l l e s du tableau 2 dans le c a s d'une c o u c h e de faible
paisseur (H petit, 3' grand).

TABLEAU 3
Le tassement est exprim par :
bord

(25) w = ^

centre

(26) w -

p'

P
E

(27) w =

p'

" H
a

" |
si |

Calcul

p ,
+

/?'

2 a
H

point
d'abscisse x
Signe

Avec

p'

ja + x|
H j

1 <
x 1>

H
|a x|

du tassement par les formules

du deuxime type. O n peut aussi bien

utiliser c e l l e s du premier type
pour p '

H l

pour p '

H l

2). L e s valeurs numriques

(tableau

du coeffi-

cient sans dimensions, p ' , sont lues,

H

en fonction

de / ? ' , dans le tableau 4

ou la figure 6.

?>

F i g . 6 Valeurs numriques d u coefficient P ' H p o u r le c a l c u l d u tassement

d'aprs les formules d u tableau 3. Ces
valeurs numriques figurent galement
dans le tableau 4.

106

ife

Coefficient
0

EXEMPLE

Considrons
L

0,01
0,02

0,2

0,1

0,3

0,4

0,5

0,033

0
0,033

0,032

0,057

0,055

0,029
0,051

0,026
0,045

0,022

0,058

0,101

0,089
0,140

0,072

0,178

0,139

0,038

0,05

0,115

0,10

0,187

0,114
0,183

0,109
0,175

0,15

0,241

0,237

0,225

0,161
0,210

0,20

0,286

0,280

0,266

0,242

0,207

0,157

0,25

0,322

0,316

0,299

0,271

0,230

0,171

0,30

0,353

0,344

0,326

0,296

0,248

0,179

0,35

0,378

0,369

0,350

0,315

0,262

0,186

0,40

0,401

0,391

0,370

0,330

0,273

0,188

0,45

0,420

0,410

0,386

0,343

0,280

0,190

0,50

0,435

0,424

0,398

0,355

0,288

0,191

0,60

0,461

0,446

0,419

0,369

0,294

0,187

0,70

0,476

0,461

0,43 4

0,378

0,298

0,179

0,111

0,80

0,488

0,472

0,444

0,384

0,90

0,498

0,481

0,449

0,388

0,298
0,292

0,167
0,153

1,00

0,503

0,488

0,452

0,390

0,291

0,144

1,5
2,0
2,5

0,508

0,493

0,455

0,387

0,270

0,095

0,506
0,503

0,492

0,452

0,254

0,062

0,489

0,450

0,378
0,375

0,247

0,047

0,5

0,489

0,450

0,5

0,489

0,450

0,374
0,372

0,242

0,236

0,029
0,014

0,5

0,489

0,450

0,371

0,233

0,006

10

0,5

0,489

0,450

0,371

0,233

100

0,5

0,489

0,450

0,371

0,233

00

0,5

0,489

0,450

0,371

0,233

T A B L E A U 4 - V a l e u r s du coefficient p '
pour le c a l c u l du t a s s e m e n t d'aprs les
f o r m u l e s du tableau 3. C e s mmes v a leurs sont reprsentes graphiquement
s u r la figure 6.

un

remblai

10 m (33 ft)

parfaitement
Quel

de POISSON,

et

lastique,

tassement
p

final
=

de

grande

de pentes
de

module

peut-on

2 070

kg/m

longueur,

latrales
E

estimer
3

129

10
8

4"

de

hauteur

2/1, reposant
40

bars

sachant
Ib/cu.

(58 psi)

que

(13

ft),

de

largeur

sur une c o u c h e de 45 m d ' p a i s s e u r

et

de

coefficient

la m a s s e v o l u m i q u e

de

Poisson,

du s o l en remblai

de

plate-forme

d'un sol
=

0,3

suppos
(fia.

7).

est :

ft

8 m

= 45m

- 40

V --

bars

0,3
F i g . 7 - R e m b l a i de l'exemple 1 dont o n
c a l c u l e le tassement au centre 0 et aux
points B et B ' situs sous l a mi-pente.

107

O n peut estimer avec une bonne approximation que ce remblai exerce

uniforme s u r une largeur B B ' = 2a = 18 m (59 ft) de valeur :
p = p gh = 2 070 x 9,81 x 4 =

80 000 P a =

0,8 bar =

C a l c u l o n s le tassement sous le centre et le bord de cette

rel.

s u r le sol une contrainte

129 x 13 =

charge

167,5 Ib/sq.ft. =

(c'est--dire

sous

normale

11,6 psi

la mi-pente

du talus

A titre d'exercice, employons les deux types de formules :

- tassement au centre : formule (23) du tableau 2 ou (26) du tableau 3 :
H

j8 =

- - = 5
a

ou

fi

0,2

1 et //

0,3, on lit dans les tableaux

pour v =

1,212

0,242

.
0,8 x 18
,
2 x 0,8 x 45
. . . .
. . .
d ou : w =
x 1,212 =
x 0,242 = 0,44 m
40
40
11,6 x 59
, ,
2 x 11,6 x 147
. ...
.
. ,
x 1,212 =
^x 0,242 = 1 ft 5 mches
58
58
0

- tassement du bord : formule (22) du tableau 2 ou (25) du tableau 3 :

3 =

= 2,5

fS- = ^

ou

= 0,4

pour v = 0,3, on lit dans les tableaux 1 et 4 :

d ou : w =

x 0,826 =
.U

= 0,826

p'

0,330

x 0,330 = 0,30 m

40
=

40
x 0,826 =

58

VAJLA !
58

0,330 =

1 ft

3 - REMARQUES SUR LES RESULTATS PRECEDENTS

Les commentaires suivants nous semblent utiles la fois pour le chercheur et pour

a) Tangente l'origine des courbes de p

l'ingnieur.

Les petites valeurs de 3 reprsentent une couche d'paisseur faible ou une charge de grande largeur. Le
champ des contraintes et dformations qui rgne dans la couche compressible est alors analogue celui
d'un chantillon plac dans un domtre et le tassement du centre de la bande s'crit :
E (1

DH

(28)

w = t

v)

avec E' =

E'

: module
(1

v)

(1 2

domtrique

v)

En comparant les formules (23) et (28), on voit qu'il en rsulte :

(1

+v)

(1 - 2

v)

* - nr=- "
C e c i n'est valable que pour les petites valeurs de 3 et fournit la tangente l'origine d e s courbes de p .
N o u s avons contrl, en agrandissant la figure 4 au voisinage de l'origine, que les courbes calcules
numriquement admettent parfaitement cette tangente.
b) Branche infinie des courbes de p
H

Les grandes valeurs de 3 reprsentent d e s c o u c h e s de grande paisseur et, lorsque [3 tend vers l'infini, p
doit tendre vers l'infini puisque, le tassement d'une charge de longueur infinie s u r un sol d'paisseur
infinie est infini, comme nous l'avons rappel dans l'introduction.

108

Pour prciser la
qui tend vers la
intgr, sur une
sous le bord de

branche infinie d e s courbes de p , nous avons calcul une valeur approche du tassement
valeur exacte lorsque l'paisseur de la couche tend vers l'infini. P o u r cela, nous avons
paisseur H, la valeur de la dformation e qui existe, dans une couche d'paisseur infinie,
la charge :
H

(1 +

(30)

2a

v) p
2 v) arctg

w =

2az

h
z

%E

dz

z* +

4a J

d'o, aprs calculs et, avec B =

2a
2ap 1 + v

(31)

(1

v) Ig (1 + P )
2

soit une valeur approche de p

1 + v

(32)

'Happ

(1 -

+ (1 2 v) B arctg

P.

pour les grandes valeurs de 3 :

v) Ig (1 +

B*) +

(1

2 v) B arctg

Le calcul numrique montre que, lorsque B augmente, la diffrence p

p pp
tend vers une valeur
peu prs constante qui ne dpend que de v. C e c i est illustr par la figure
8. L'cart relatif entre p et
PHapp
donc en diminuant, comme on le prvoyait.
H a

L'intrt
triaxiale

pratique

de

de
calcul

cette

remarque

est

que

l'intgrale

des

tassements

au

cas

d'une

(31)

couche

n'est

autre

compressible

que

l'application

de

la

mthode

homogne.

F i g . 8 - C o m p a r a i s o n entre les valeurs de p

(
) et de PHapp (--).
Pour v
0,5, les deux courbes
sont pratiquement confondues.
H

109

c) Allure des courbes de p '

Le coefficient p '

tant dduit de p

l'asymptote de p '

par la formule d'inversion (21), il en rsulte que, pour v donn :

a pour valeur la pente l'origine de p ,

la tangente l'origine de p '

e s t verticale.

d) D f o r m e de la surface du massif
A l'aide de la formule (27), nous avons trac la dforme de la surface du massif pour d i v e r s e s valeurs de
v et de (fig- 9). O n notera, en particulier, le soulvement du s o l , surtout pour
2a
respond, en pratique, au chargement rapide d'une argile sature.

S o u l v e m e n t

du

v = 0,5 : c e c a s c o r -

sol

w E
F i g . 9 - Dforme de l a surface d u massif p o u r diverses valeurs de v et de H / 2 a .

9 P

I 10

Par ailleurs, la figure 10 montre comment on peut diminuer le tassement en augmentant la charge dans le cas
d'un sol incompressible (v = 0,5).

F i g . 10 - Tassement d'un s o l incompressible en fonction de l a dimension de l a charge.

e) D p l a c e m e n t s des points s i t u s l ' i n t r i e u r de la couche l a s t i q u e

Dans le c a s o le coefficient de P o i s s o n du sol compressible vaut v = 0,3, nous avons calcul le dplacement de tous les points du sol par la mthode numrique de la double grille (Calcul d e s contraintes 1).
Les rsultats ainsi obtenus, pour le tassement en surface, concordent parfaitement avec ceux obtenus par les
sries de Fourier comme on le voit en comparant les figures 4 et 11.
111

N o u s n'estimons p a s utile de reproduire ici les valeurs du dplacement horizontal que nous avons galement
obtenues, avec v = 0,3, pour tous les points du massif [2], En effet, nous avons montr [5] que les charges
tangentielles ont une grande influence s u r les dplacements horizontaux alors qu'elles n'en ont pratiquement
pas s u r les tassements : il faudrait donc, pour tudier valablement l e s dplacements horizontaux, faire intervenir l'inclinaison ds contraintes de contact entre la charge et le sol, c e qui n'a p a s t fait ici. S i g n a l o n s ,
cependant, q u e les valeurs obtenues tendent bien, pour H = , vers la limite que nous avions indique par
ailleurs [7].
f) Comparaison de nos r s u l t a t s avec ceux d'autres t u d e s analogues
C e s dernires annes, plusieurs auteurs ont abord ce problme, mais de faon moins dtaille qu'ici.
D e s rsultats trs diffrents d e s ntres ont t obtenus par B e n d e l [8], sauf pour les grandes valeurs de P.
Par ailleurs, Harr [9], tendant d e s travaux de Egorov [10], a cru pouvoir crire, semble-t-il, que le coefficient
de P o i s s o n , v, intervenait seulement sous forme d'un facteur (1 v ) dans la formule de tassement. C e c i est
vrai pour un milieu semi-infini, mais ne l'est p a s pour une couche d'paisseur finie comme on le voit dans les
formules (15) et (19). L e s valeurs numriques que nous donnons ( f l g . 6) confirment que (1 v ) est loin d'tre
un facteur.
2

A u contraire, Ueshita et Meyerhof [11] ont publi rcemment d e s rsultats qui concordent avec les ntres et,
enfin, dans une remarquable tude gnrale s u r les problmes de c o u c h e s c o m p r e s s i b l e s lastiques, Poulos
[12 - fig. 42] donne un coefficient de tassement qui permet de retrouver approximativement les valeurs
donnes i c i .
Malgr les diffrences que nous venons d'indiquer avec certains d e s rsultats publis, nous sommes en droit
d'estimer que n o s rsultats sont exacts. En effet :
dans le c a s v = 0,3, nous avons retrouv l e s mmes rsultats par la mthode
grille,
112

numrique de la double

 nous avons vrifi que les courbes p et p ' admettent bien la tangente l'origine et la courbe asymptotique calcules indpendamment et explicitement.
H

Reconnaissons enfin que la difficult de c e problme nous avait galement fait commettre quelques erreurs
dans les premiers calculs faits G r e n o b l e [13]. Nanmoins, c'est en continuant dans la voie ouverte par Rouget
que le calcul d e s contraintes indiqu dans la deuxime partie de cet article a pu tre men bien.

C A L C U L DES CONTRAINTES

1 - DETERMINATION THEORIQUE
DOUBLE GRILLE

DES CONTRAINTES

ET DEPLACEMENTS PAR

LA

METHODE

DE

LA

a) Principe de la m t h o d e
N o u s avons vu (Calcul d e s tassements 1 a), qu'une mthode pour rsoudre le problme consistait ramener
le systme (1), (2) et (4) une seule quation biharmonique. Une autre mthode consiste liminer les (s , 7)
et les (o-, T ) dans les quations (1), (2) et (3). O n obtient ainsi les quations de Navier-Lam :
2 (1 v)
1 - 2

(33)

w
2

2 (1 v) u
1 2 v

1 - 2

z 5 x

1
1 - 2 v

b z x

S i l'on parvient rsoudre c e systme de deux quations aux drives partielles deux inconnues, u et w,
on en dduit les dformations ( s , 7), par les formules (3) et les contraintes (o-, t ) , par les formules (2). Dans
le c a s qui nous intresse, les conditions aux limites sont telles que le systme (33) n'admet p a s de solution
explicite. O n utilise alors une mthode numrique qui consiste reprsenter le milieu lastique par un double
rseau rgulier de points masse et de points contrainte (fig. 12). Le principe de cette mthode de la double
grille est donn en dtail par ailleurs [2], [13], [14].
L'criture de l'quilibre du point masse P (I, J) sous l'action d e s quatre points contrainte environnant et sa
traduction en termes de dplacements par l'utilisation d e s relations (2) et (3) crites en diffrences finies c o n duit :
2 (1 - v)

w' (I + 1, J -

1 2v

1, J + 1) 2 w' (I,

J)j

w ' (I + 1, J + 1) + w' (I - 1, J - 1) 2 w ' (I, J)

[| w

(I + 1, J) u' (I, J + 1) + u' (l 1, J) -

+
(34)

1) + w' (I -

u' (I, J - 1)

1 2 v L

(
2_(1 - v)
1

V )

|fu' (I + 1, J + 1) + u' (I -

| u - (| + 1, j _ 1) + ' (I U

1, J -

1, J + 1) -

1) -

2 u' (I,

[ w ' (I + 1, J) + w ' (I, J + 1) + w' (I -

2 u' (I,

J)J

J)j
1, J) w' (I, J -

1)

1=

avec : w ' , u' composantes

1 - 2 v du dplacement s u r les axes inclins 45. O n en dduit immdiatement w et u.
113

Vecteur
dplacement

r
0

I =
-1

1 = 1

* \

\t

J-

'

J-

%*

r
J-

1
'

1w

r
J e

-3

-i

1
' t

Substratum

rigide

F i g . 1 2 - Schma de p r i n c i p e de l a double grille

C e s dernires relations ne sont autres que l'criture en diffrences finies des quations (33) de Navier-Lam.
La double grille permet donc de reprsenter discrtement un milieu continu lastique. D e la mme faon, on
aurait pu utiliser, au lieu de la loi d'lasticit linaire de Hooke (2), une autre loi, plus ou moins complique
(lasticit non linaire, lasto-plasticit). D'une manire gnrale, on peut donc dire que la mthode de la
double grille, plutt qu'une reprsentation physique du milieu, fournit une schmatisation pratique qui permet
de retrouver les quations d'quilibre et de comportement exprimes en diffrences finies et d'crire facilement
les conditions aux limites. L'approximation est d'autant meilleure que la maille est plus serre.

114

b) Application au cas d'une charge uniforme reposant sur une couche l a s t i q u e

V o y o n s maintenant comment la double grille a t utilise pour l'tude de notre problme.
Conditions aux limites :

A la surface, la charge est reprsente par un ensemble de forces concentres Q :

(35)

Q= X p

pour 1 = 0, J < J < Je

Q = A p/2

pour I = 0, | J | = J

Q= 0

pour I = 0, | J | > J

avec : ~k ct de la maille.

Sur le fond, la rigidit et la rugosit du substratum font que le dplacement

j

w (I = N, J) = 0

u (I = N,

de tous

les points est nul :

(36)

J)

= 0

L a t r a l e m e n t , la grille devrait tre infinie mais pratiquement, il faut que le nombre de mailles soit infrieur
1 000 pour que la relaxation ne soit pas trop longue. O n doit donc imposer cette limite artificielle une condition aussi voisine que possible de celle qui existerait si la grille tait effectivement infinie. Pour cela, nous avons
crit que, s u r la colonne J = M , le dplacement horizontal est nul et le dplacement vertical gal celui de la
colonne J = M 1 -. nous avons constat en faisant varier M entre 7 et 10a/), que ceci n'influait pratiquement pas s u r les rsultats. Enfin, la symtrie du champ de contraintes par rapport l'axe 0z permet de limiter
la grille la partie J > 0.

Rsolution :
L'ensemble d e s relations (34) forme un systme linaire de 2 N M quations 2 N M inconnues de la forme
[A] [X] = [ B ] . La matrice [A], de rang 2 N M , est tridiagonale par blocs et s a rsolution s e fait par la mthode
itrative de G a u s s - S e i d e l avec utilisation du paramtre de sur-relaxation, M . L a rapidit de la convergence est
fonction de la valeur du coefficient de P o i s s o n , v, de la dimension de la grille et de la valeur de w. N o u s
nous sommes limits au cas v = 0,3 et, pour chaque grille, nous avons dtermin la valeur optimale de M
comprise entre 1,5 et 1,8 en utilisant, notamment, les tudes de V a r g a [15]. L a prcision des rsultats est fonction de la finesse du maillage (X/a) et de la grandeur du test de convergence. Nous avons utilis d e s grilles
de N M = 200 1 000 points et obtenu une prcision s u r les rsidus d e s quations (34) de l'ordre de 10" .
3

Exploitation d e s rsultats :
Une fois le calcul men bien, nous d i s p o s o n s , pour chaque valeur de H / a , d e s valeurs numriques d e s
contraintes, d e s dplacements et des dformations en tout point de la couche compressible. L e s dplacements
ont dj t donns dans la premire partie de cette tude et les dformations ne sont gnralement pas d'un
grand intrt pratique. Nous ne donnerons donc ici que les contraintes et, grce au principe de superposition,
nous pourrons condenser les rsultats numriques en trois tableaux seulement. En effet, la contrainte en un
point quelconque de la couche compressible s e dduit de la contrainte la verticale du coin par une addition
analogue celle indique pour les tassements sur la figure 5. (La marche suivre pour faire c e s calculs est
prcise dans le C a l c u l d e s contraintes 2.)

2 - RESULTATS NUMERIQUES
D'aprs l'tude que nous venons d'exposer (Chapitre C a l c u l d e s contraintes 1), les contraintes s'expriment par les formules du tableau 5. Notons que, dans ce tableau, les expressions (37) (42) ne sont que
des c a s particuliers d e s expressions gnrales (43) (45) (Voir plus loin, chapitre C a l c u l d e s contraintes ,
3 b, pour l'influence du coefficient de Poisson).
115

.10

B
( 4 9 f t ) 15m

E X E M P L E 2 : Reprenons le remblai de l'exemple 1, mais

en considrant que la couche de 45 m est constitue de
trois s o u s - c o u c h e s supposes parfaitement
lastiques,
mais d e modules diffrents (fig. 13).

E'-" 50 bars

On demande le tassement. Faisons le calcul en passant

par l'intermdiaire de la contrainte a ( 1 de l'introduction). Il faut alors tracer la courbe o- la verticale d e s
points O et B en utilisant les formules (37) et (40) du
tableau 5 : on obtient les deux courbes reprsentes sur
la figure
14. O n en tire la valeur moyenne de o- dans
chaque couche, respectivement s o u s le centre et s o u s le
bord, d'o le calcul du tassement :

( 725 psi)

(33ft)IOm

E'= 35 bars

(500psi

(65ft)20m

& ? ; E ' = 55bars

(800psi)

Tassement du centre :
_
w -

15 x 0,7

S ubstratum

F i g . 13 - C o u p e d u s o l supportant le r e m b l a i de
l'exemple 2. L e s modules indiqus sont les modules
domtriques. L e tassement est calcul en 0 et B .

50

10 x 0,44
_
35

49 x 10
725

33 X 6,3
zrr^.
510

20 x 0,31

55

1+

65 x 4,5
rr^
800

_ o,45 m
.
1 ft. 6 inches

Tassement des bords :

15 x 0,39
50

49 x 5,6
725

10 X 0,34
35

, 33 x 4,8
510

20 x 0,27
55

= 0,31 m

65 X 3,9
= 1 ft.
800

Notons que c e s valeurs sont trs voisines de celles obtenues dans l'exemple 1 pour une c o u c h e homogne de
module E = 40 bars.

TABLEAU 5
0,40bar(5,8psi)

Avec P

Les contraintes sont exprimes par :

bord

centre

point
d'abscisse x

(37)

<r

H 0

(38)

(39)

<r

k 2

(40)

(41)

(42)

<r

(43)

(44)

(45)

z x

Z X

z x

i
2 a

2 p k o
0

2 p

p [kffi

[ j

p f {1)

H
a

k 2
H

k<f]
k

|a

En tout point : <r

= v (a

+ o )
x

si | x | <
.
^
si
x
>

x|

ffl]
K

$ ]

|a .
( +
Signe
-

15m

x\

( 2 )

a
a

(46)

C a l c u l des contraintes en tout point (x, z), de la couche lastique. Les valeurs
numriques des coefficients sans dimensions kno- k m et k m sont donnes en
fonction de P et de z / H dans les tableaux 6, 7 et 8 et les figures 15, 16 et 17.

116

0,80bor(11,6 p s i )

F i g . 14 - D i s t r i b u t i o n des contraintes sous

le centre (
) et sous le point B ()
de l'exemple 2. L a premire courbe est
obtenue avec la f o r m u l e (40) et P = 5
et l a seconde avec l a f o r m u l e (37) et
P
=
2,5.

0,5

1,25

1,5

0 500
0,500 0 500
0 502 0,5015 0 501
o 5025 0,502 0 5015

0
0 ,0 5
0 .1

0 500
o 503
o 504

0 .1 5
0 .2
0 ,2 S

0 5035 0 502
o 503
0 5015
0 5025 0 501

0 ,3
0 .3 5
0 ,4

0 502
o 502
0 502

0,5015 o 501
0,501
0 500
0,5005 0 499

0 5005 0,500
0,499
o 500
o 499
0,498

1,75

0 500
0,500
0 5005 0,499
0 500
0,498

10

20

0,500 0,500
0,4985 0,498
0,497
0,496

0 500
o 500
0 4975 0 497
0 4955 0 495

0,500 0 500
0,4965 0 496
0,492 o 488

0 500
0 495
o 477

o 500
0 494
o 472

0 500
0,494
0 462

0 500
0 482
0 402

0,500
0,407
0,275

2,25

2,5

2,75

3,5

o 499
0 498
o 496

0,497
0,495
0,489

0,495
0,490
0,481

0,492
0,485
0,473

0 490
0 480
0 466

0 486
0 472
o 453

0,477
0,460
0,437

0 472
0 450
0 417

0 447
0 407
o 369

o 420
0 380
o 332

0 397
0 337
0 295

0 325
0 267
0 225

0,194
0,152
0,122

o 497
o 495
0 491

o 492
o 487
o 480

0,481
0,469
0,457

0,470
0,457
0,442

0,460
0,442
0,425

0 449
0 429
o 410

0 430
o 407
o 385

0,410
0,382
0,358

0 378
0 347
o 322

o 333
o 303
0 278

0 295
o 265
o 240

0 260
0 235
0 216

0 195
0 170
0 154

0,103
0,089
0,080

0 .5
0 ,5 5

0 5015 0 498
0 501 0 500
0 5005 0 499

0,4965 0 487
0,495
0,482
0,493
0 475

o 473
0 466
o 458

0,445
0,433
0,423

0,428
0,414
0,400

0,409
0,394
0,380

0 392
0 374
0 357

0 367
0 350
0 333

0,338
0,318
0,303

o 301
0,285
0 269

0 256
0 240
0 227

o 220
0 205
0 191

0 199
0 185
0 172

0 143
0 132
0 124

0,072
p,066
0,061

0 ,6
0 .65
0 ,7

0 500
0 500
0 499

0 498
0 497
0 495

0,491
0,487
0,481

0 467
0 460
0 453

o 449
o 440
o 430

0,412
0,401
0,389

0,387
0,366
0,375
0,354
0,3634 0,343

o 343
0 331
0 320

o 319
o 308
0,298

0,289
0,277
0,266

0 256
0 245
0 235

o 216
o 205
0 196

o 180
o 170
o 162

o 162
0 153
0 146

o 1175
o 112
0 107

0,058
0,055
0,052

0 ,75
0 ,8
0 ,85

0 498
0 497
0 496

0 492
0 489
0,485

0,475
0,469
0,462

0 445
o 437
0 428

0 420
o 410
0 400

0,379
0,369
0,360

0,352
0,342
0,333

0,332
0,321
0,312

0 310
0 301
0 292

0 289
0 280
0 271

0,255
0,246
0,238

0 225
o 217
0 210

0 188
o 181
o 175

0 157 ? o 140
o 151 o 135
0 147
0 13(0

0 1045
o 1025
0 101

0,051
0,0505
0,050

0 ,9
0 ,9 5

0 481
0 495
0 4945 0 476
0 494
0 471

0,454
0,446
0,438

0 419
0 410
0 400

0 390
0 380
0 370

0,350
0,340
0,332

0,325
0,317
0,308

0,303
0,295
0,288

0 283
0 275
0 269

0 262
0 256
0 251

0,230
0,223
0,217

o 204
o 200
o 197

0 171
o 167
o 165

o 145
o 1425
o 141

0 100
0,099
o 099

0,050
0,0495
0,049

,4

0 1265
0,125
0 124

T A B L E A U 6 - V a l e u r s du coefficient k o pour le calcul de cz d'aprs les formules du tableau 5.

C e s mmes valeurs sont reprsentes graphiquement s u r la figure 15.
H

I 17

0,5

0,75

2,5

2,75

3,5

0
0 280
0 27 5

0
0 282
0 275

0
0 284
0 275

0
0 288
0 265

0
0 289
0 259

0
0 284
0 243

0
0 278
0 226

0
0 269
0 204

259
234
207

252

227

0,148
0 112

149
0 108
0 07 7

191
156

0 172
0 127
0,093

0 219
0 17 8
0 143

221
0,191

182
158
143

0 171
0 147
0 136

131
109
091

0 116
0 094
0 077

0,086
0 069
0 055

0 07 2
0 058
0 046

0 060
0 047
0 037

117
0 101
0 087

0 107
0 091
0 07 7

0 099
0 084
0 071

07 6
065
0 054

0,063
0 052
0 043

0 044
0 035
0 029

0 036
0 028
0 023

0 029
0 023
0 019

0,084
0,082
0,065

0 07 5
0 070
0 057

0 066
0 057
0 049

0 060
0 051
0 044

045
0,038
0 032

0 035
0 029
0 023

0 023
0 017
0 013

0 019
0 015
0 012

0 015
0 013
0 010

0,057
0,051

0 044
0 038
0 036

0 037
0 032
0 030

0 027
0 024
0 022

0 019
0 017
0 016

0 011
0 010
0 009

0 009
0 008
0 007

0 007
0 006
0 005

0 036
0 037
0 038

0 029
0 029
0 030

0 022
0 023
0 023

0 015
0 015
0 016

0 009
0 009
0 007

0 007
0 007
0 005

0 005
0 005
0 006

2,25

1,25

',5

0
0 ,05
0 l1

0
0 ,234
0,274

0
0 240
0 276

0
0 246
0 27 7

0
0,252
0 278

0
0 261
0 277

0
0,267
0 276

0
0,272
0 276

0,280
0,276

0 ,15
0 ,2
0 ,25

0,266
0 255
0 240

0 267
256
0 241

0 268
0 256
0 244

0 268
0 257
0 243

267
255
240

0 266
0 251
0 234

0 264
0 247
0 237

0,262
0,242
0,217

0 >3
0 227
0 . 35 0 2 1 2
0 , 4
0 200

0 228
0 214
0 202

0 230
0,216
0 204

0 228
0 213
0 198

224
207
190

0 216
0 198
0 180

0 206
0 185
0 164

0,193
0,169
0,149

0 190
0 ,45
0 > 5 0 181
0 , 55 0 1 7 4

0 190
0 17 9
0 171

0 189
0 178
0 167

0
0
0

184
170
158

173
157
143

162

0 146
0 129
0 114

0,129
0,111
0,096

0 146
0 131

0 76
0 ,65
0

0 163
0 157
0 153

0
0

155
146
0 139

0 145
0 134
0 125

0
0

130
119
0 109

0 118
0,105
0 095

0 102
0,093
0 080

,1

0 168
0,164
0 162

0 , 75 0 , 1 6 2
0 > 8 0,163
0 , 85 0 169
0 ,9
0 17 8
0 > 9 5 0 193
1
0 209

0
0

0
0

0
0

1,75

0
0
0
0
0

151
0 150
0 153

133
130
129

0,118
0 112
0 110

0 101
0 094
0 088

0 086
0 080
0,07 5

0 07 1
0 065
0 061

0,048

0 051
0 045
0 042

0 160
171
184

0 132
0 140
0,163

0 114
0 119
0 130

0 092
0,097
0 106

0 075
0 077
0 083

0 061
0 062
0 065

0,048
0,049
0,050

0 042
0 043
0 044

256
0 227
0 197

0 161
0,139
125

0
0
0

194

10
0
0

0
0

234
126
07 6
052
037

02 7
019
0 014
0

011
009
0 007
0

005
004
0 003
0

0
0
0
0
0

0025
002
0015
0015
0015
0015

T A B L E A U 7 - V a l e u r s du coefficient k m pour le
calcul de t z x d'aprs les formules du tableau 5.
C e s mmes valeurs sont reprsentes graphiquement sur la figure 16.

F i g . 16 - G r a p h i q u e des valeurs de k n i p o u r le
calcul de T . Signalons que ces courbes sont peuttre imprcises au voisinage d u m a x i m u m . V o i r
galement le tableau 7.
2 x

I 18

0,5

0,75

1,25

1,5

1,75

2,5

3,5

10

0 500
0 195
0 188

0 500
0 210
0 192

0 500
0 221
0 198

0 500
0 237
0 212

0 500
0 252
0 222

0 500
0 272
0 233

0 500
o 280
0 234

0 500
0 288
0 224

0 500
o 291
0 216

0 500
0 293
0 215

0 500
0 293
0 196

0 500
0 265
0 162

0 500
0 255
0 139

0,500
0,251
0,109

0,500
0 195
0 028

o 187
o 188
0 189

0 183 0 188
o 17 6 0 178
0 174 0 168

0 94
0 177
0 162

0 199 0 200
0 17 6 0 168
0 154 0 140

0 197
0 159
0 128

0 178
0 133
0 097

0,162
0 111
0 073

0 147
0 091
0 054

0 121
0 074
0 040

0 091
0 048
0 025

0 064
0 029
0 014

0,032 0,006
0,004 0,0005
0,0015 -0 ,001

0,35
0,4

0 190
0 191
0 193

0 173
0 172
0 171

0 160
0 153
0 147

0 147
0 136
0 126

0 134
0 116
0 103

0 119
0 101
0 086

0 103
0 084
0 069

0 070
0 051
0 039

0 047
0 031
0 021

0 033
0 020
0 013

0 021
012
0 006

0 Oil
0 006 0,0005 -0 0015
0 005 0 003 0,000
0 0025 0 001 -0,0005

0,45
0,5
0,55

0 194
0 195
0 196

0 170
0 169
0 169

0 140
0,136
0 134

0,118
0 113
0 109

0 093
0 084
0 078

0 074
0 066
0 061

0 059
0 053
0 049

0 033
0 030
0 030

0
0

017
016
0 017

0 Oil
0 009
0 010

0 005
0 004
0 006

0,6

0,65
0,7

0 198
0 200
0 202

0 170 0 135
o 1735 0 138
0 177 0 141

0 106
0 107
0 111

0 075
0 077
0 082

0 060
0 062
0 067

0 049
0 051
0 056

0 032
0 036
0 040

0 020
0 023
0 029

0 013
0 016
o 022

0 008
0 010
0 015

0 0075 0 007
0 010 0 009
0 014 0 013

0,004
0,006
0,010

0,75
0 ,8
0 ,85

o 203
o 204
0 206

0 180
0,183
0 189

o 146
0 154
0 162

0 119
0 128
0,138

0,0905 0,075
0 102 0 087
0 116 0 100

o 064
0 075
0 088

0 047
0 057
0 070

0 034
0 043
0 057

0 029
0,036
0 048

0 021
0 028
0 037

0 020
0 027
0 033

0 019
0 022
0 028

0,014
0,018
0,024

0,9

0 208
0 211
o 214

0 194
o 200
o 206

0 173
0,183
0 195

0 152
0 167
0 186

0 132
0 149
0 167

0 103
0 122
o 143

0 086
0 105
0 125

0 072
0 091
0 111

0,063
0 082
0 103

0 049
0 062
0 080

0 042
0 052
0 063

0 034
0 043
0 051

0,029
0,036
0,043

0
0,05

0,1
0,15
0,2
0,25
0,3

,95

0 115
0 133
0 153

TABLEAU 8
km pour le
formules du
valeurs sont
ment

0015 0 000 -0,0005

0 0025 0 001 0,000
0 005 0 004 0,002

- V a l e u r s du coefficient
calcul de a* d'aprs les
tableau 5. C e s mmes
reprsentes graphiquesur la figure 17.

F i g . 17 - G r a p h i q u e des valeurs de kn2 p o u r

le c a l c u l de o-x. V o i r galement le tableau 8.

119

120

3 - REMARQUES SUR LES RESULTATS PRECEDENTS

a) Comparaison avec le milieu semi-infini
Pour faire cette comparaison, nous avons trac les figures
geant l'ordonne.

18 et 19 dduites d e s figures

15 et 17 en chan-

On voit que les contraintes o- sont infrieures aux valeurs correspondantes du milieu semi-infini dans la partie
haute de la couche et suprieures au voisinage du substratum. C e phnomne local est d essentiellement
la rugosit de contact entre la couche compressible et le substratum.
x

A u contraire, les valeurs de o- sont toujours suprieures celles qui rgneraient la mme profondeur dans
une couche lastique d'paisseur infinie. Toutefois, ceci n'est valable que pour les points situs proximit de
l'axe 0 . Pour les grandes valeurs de l'abscisse x, c'est le contraire qui s e produit, comme le montre la figure
20
est alors plus grand dans le milieu semi-infini que dans la couche d'paisseur finie. C e c i est tout
fait logique c a r la somme d e s contraintes o- s u r le substratum est gale 2pa par tranche d'paisseur unit,
comme la profondeur correspondante dans le milieu semi-infini. Donc, un excs en un point doit tre compens par un dfaut ailleurs.
z

CTz

F i g . 20 D i s t r i b u t i o n de o- a u contact de l a couche compressible et d u substratum (

).
O n a galement trac l a courbe <r dans u n m i l i e u semi-infini l a mme profondeur (---).
x

121

122

La comparaison avec le milieu semi-infini e s t galement illustre par la figure 21.

Enfin, il nous a sembl utile d'tablir
une couche d'paisseur finie partir
pas sous les yeux le tableau 6, mais
du tiers et du quart dont le principe

une rgle simple permettant d'obtenir une valeur approche de o- dans

des valeurs du milieu semi-infini. C e c i peut intresser l'ingnieur qui n'a
qui possde l e s rsultats bien connus du milieu semi-infini. C ' e s t la rgle
est le suivant (fig- 22)
z

on trace la courbe de o- sous la charge comme si le milieu tait semi-infini,

z

3H
la valeur de o- pour z = est reporte la profondeur z = H, c e qui donne le point C ,
4
la courbe de o- dans la couche c o m p r e s s i b l e comprend l'arc A B confondu avec la courbe du milieu semiinfini entre z = 0 et z = H/3 et l'arc B C obtenu par raccordement d e s points B et C en s'loignant progressivement de la courbe du milieu semi-infini,
notons que c e procd n'est valable que pour <j . , v = 0,3 et | x | < a.
z

Partie commune aux

deux

y
dans

courbes

un

milieu s e m i -

infini

couche
d'paisseur H

'4

F i g . 22 . Rgle d u tiers et d u quart. L a courbe de o dans

u n m i l i e u semi-infini tant trace, o n en dduit l e point C
en menant l a verticale C ' C . L a courbe de <r dans l a couche
compressible suit l e trac A B C .
2

b) Influence du coefficient de Poisson sur les contraintes

Les calculs tant trs longs, nous nous s o m m e s limits au cas v = 0,3, valeur moyenne gnralement admise
pour les sols. Q u e l q u e s essais, faits avec d'autres valeurs de v, nous ont montr que cette influence est a s s e z
faible. C e c i est confirm par P o u l o s [12], notamment pour o- o le rle de v semble tre absolument ngligeable.
z

O n pourra donc retenir les valeurs donnes ici pour les contraintes, et notamment pour o- , comme une bonne
approximation des valeurs relles quel que soit le coefficient de P o i s s o n du s o l , v.
z

c) Comparaison de nos r s u l t a t s avec ceux d'autres t u d e s analogues

C o m m e pour les tassements, les valeurs que nous avons obtenues sont a s s e z diffrentes de celles donnes
par B e n d e l [8]. C e l u i - c i prsente s e s rsultats sous la forme du produit (o- /p) (H/2a) qui serait constant pour
une a b s c i s s e x donne. C e c i nous parat conduire une absurdit lorsque H/2a est petit.
z

A u contraire, nos rsultats sont v o i s i n s 10 % prs d e s valeurs donnes par Egorov [10] pour o- au point
(a; = 0, z = H) et ils sont en bon a c c o r d avec les valeurs que l'on peut calculer pour o- et T
en tout point
partir des coefficients de P o u l o s [12 - fig. 32 et 4 0 ] .
2

2 x

N o u s avons donc tout lieu de croire la validit de nos rsultats, d'autant plus que le mme calcul de relaxation a fourni simultanment les valeurs du tassement dont nous avons contrl l'exactitude par ailleurs ( C a l c u l
des tassements 3 f).

123

CONCLUSION

N o u s esprons que cette tude, que nous avons essay de faire aussi complte que possible, pourra fournir
des informations utiles ceux qui ont s ' o c c u p e r de tassements. Toutefois, il ne s'agit que d'une approche
thorique et la possibilit d'en appliquer les rsultats des cas concrets devrait faire l'objet d'une d i s c u s sion trs dtaille. Il nous parat prfrable, pour le moment, de nous limiter l'expos prcdent qui parat
dj suffisamment long.

REMERCIEMENTS
Ce travail a t ralis dans le cadre de la thse de Doctorat
s Sciences du premier auteur et de la thse
de Doctorat de Spcialit du second. Les auteurs tiennent remercier
M. M. Rey qui s'est charg de l'excution des dessins et M. U. Ilker qui a bien voulu relire le manuscrit.
Que soient galement
remercis
M. F. Schlosser, M. E. Leflaive et M. J.M. Runacher qui ont bien voulu nous faire de trs nombreuses
et
trs intressantes
suggestions.

BIBLIOGRAPHIE

[1]
[2]

J.-P. G I R O U D , Fondation
VI T.B.T.P., SF 70, 253,
H .

rectangulaire
linairement
(janv. 1969), 79-111.

charge

tassement

et

contraintes,

Annales de

Problmes
bidimensionnels
d'lasticit
linaire ou non linaire traits par
thse de Doctorat de Spcialit, Facult des Sciences de Grenoble, (1969).

WATISSEE,

numrique,

calcul

[3]

J.-P. G I R O U D , Calcul rapide des contraintes

provoques
Labo. Routiers P. et Ch., 35, (dc. 1968), 83-88.

[4]

J.-P. G I R O U D , Calcul
( paratre).

[5]

J.-P.

[6]

D. L. H O L L , Plane-strain
148, Ames (1941).

[7]

J.-P. G I R O U D , Dplacement
horizontal d'une droite particulire
de la surface d'un massif lastique semiinfini linairement
charg, Comptes rendus l'Acadmie des Sciences, 268, Srie A , (27 janv. 1969),
252-255.

[8]

H. B E N D E L , Die Berechnung
von Spannungen
Versuchsanstalt
fr Wasserbau und Erdbau,

[9]

M. E .

[10]
[11]

GIROUD,

HARR,

pratique

des contraintes

dans un sol supportant

une charge

Bull. Liaison des

de grande

longueur

thse de Doctorat s Sciences (soutenance en 1970).

Foundations

distribution

of stress in clastic media, Iowa Engineering Experiment Station,

of theoretical

K. E . E G O R O V , Sur la dformation
troiizdat, Moscou (1958).

und Verschiebungen
55, Zurich (1959).

soil mechanics,

des couches

in

Erddmmen,

Mitteilungen

der

Mc Graw-Hill, New York (1966).

d'paisseur

and G. G. M E Y E R H O F , Surface displacement

Highway Research Record, 228 (1968).

K . UESHITA

loads,

dans le sol par un remblai,

le

finie,

(en russe),

of an elastic

Mekanika Gruntov, 34 Goss-

layer under uniformly

[12]

H . G. P O U L O S , Stresses and displacements

technique 17 (dec. 1967), 378-410.

[13]

M. R O U G E T , Integration
numrique
des champs de contraintes
et de dformations
lastiques
par la
mthode de la double grille, thse de Doctorat de Spcialit, Facult des Sciences de Grenoble (1967).

[14]

D. N . de G.

[15]

R. S.

124

VARGA,

ALLEN,

Matrix

Relaxation
iterative

methods,
analysis,

in an elastic layer underlain

distributed

by a rough rigid base,

M c Graw Hill, New York (1954).

Prentice Hall, Engelwood - Cliffs - New Jersey (1967).

Geo-