poinçonnement d'un sol élastique anisotrope

par

M. Dahan
Ingénieur, Attaché de Recherches C.N.R.S.
Laboratoire de Mécanique des solides*

1 Introduction 2 Chargements axisymétriques.

Pour la prévision des tassements, il est courant en Méca- Nous nous proposons de généraliser au cas d'un charge-
nique des Sols de simplifier le comportement des terres et ment quelconque la méthode que nous avons utilisée
de les considérer comme des corps élastiques pour obtenir (Dahan, 1975; Dahan et zarka, 1977]t pour la résolution
des résultats approchés. du contact élastique d'une sphère avec un massif semi-in-
Jusqu'à présent, la plupart des calculs ont été faits en utili- fini.
sant comme modèle mathématique le massif semi-infini, Considérons un massif serni-infini transversalement iso-
homogène, élastique et isotrope. Des solutions à ce pro- trope dont la surface est délimitée par le plan horizontal (r,
blème ont été données par Boussinesq dès 1885. En fait, 0) et dont l'axe de symétrie élastique (zl est vertical. Ce
le sol n'est pas parfaitement isotrope, et cela peut expli- matériau est défini par cinq coefficients élastiques aij .

quer certa ins désaccords entre la théorie et l'expérience.

Les relations contraintes-déformations associées aux pro-
Le simple dépôt, la sédimentation, la consolidation du sol blèmes à symétrie axiale sont :

sous I'effet du poids propre, peuvent entraîner l'apparition

d'une anisotropie. Les propriétés du sol dans les directions
trr - ârrorr + arzogg + ârt6ZZ,
horizontales et verticales ont alors des différences signifi-
catives. En supposant la symétrie des propriétés élastiques
teg
-- ârz6rr + arrogg + an677,
(1) Ezz: â,rorr + a*ogg + â.,,6.,..,
du sol autour de la verticale, on peut le considérer comme
un massif semi-infini transversalement isotrope. Nous Trz
avons vérifié cette propriété sur divers échantillons de sols
par !a mesure des coefficients élastiques. où les déformations sont définies à partir des compo-
santes (ur, u0 o, w) du déplacement par les relations
Le but de cet étude est de donner la distribution des - :

contraintes et des déplacements à I'intérieur d'un sol à iso- âr'.

tropie transversale lorsqu'il est soumis à des chargernents c
-
axisymétriques sur sa surface libre, puis de déterminer la "rr A,
répartition des pressions sous un poincon rigide, ainsi
qu'une relation entre la force de poinconnement et la pro- ur
fondeur de pénétration du poincon dans le sol. t00 _
Les résultats obtenus montrent que la distribution des (21
pressions sous le poincon est la même que si le massif
était isotrope; pour le tassement à la surface du massif, il ozz
câw
suffit simplement de remplacer le module 2(1 vr) I e Au a=
cas isotrope par un certain coefficient q.Nous - générali-
sons ainsi tous les résultats connus pour le massif isotrope, â,,. âw
à partir de la formule de Boussinesq, en ce qui concerne la' Trz: a=
* a,
distribution des contraintes sous le poincon et le tasse-
ment en surface, au milieu transversalement isotrope. Les équations d'équilibre sont:
Ces résultats sont illustrés par un certain nombre de
courbes montrant l'influence de I'anisotropie sur les con- àoru orr - oOe
traintes et les déplacements. 1âorr
* *T-:0
àr a=
* Le Laboratoire de Mécanique des Solides est une équipe de Recherche (3)
associée au c.N.R.S. (Groupe Commun Ecole Polytechnique, Ecole àor= àor, 6rz

ar *E*
Nationale Supérieure des Mines de Paris, Ecole Nationale des Ponts et
Chaussées).
,. -o
REVUE FRANçA|SE DE GEOTECF|NIOUE NUMERO 7
et les équations de compatibilité transformées par les rela-
tions (1 ) s'écrivent :
Fi-1-asî'
- 1'2
(a,, arrl (orr o00) t
+ (arzorr + arroge + arrorrl 9i - s',t iarro -2arral uoo,i
or
-o, (1 1)
(41 f -(O ac) /\Æ,
a2 l" - (t b) I t ,

_ (ârrorr + arzogg + arrot

oz'"
l
q - À(an - arrl (s, + sr)
. Az(ânorr + airooo + arro=rl-âqq A2o -- O. et où pH(m) est la transformée de Hankel d'ordre zéro du
# #- chargement p(r) :

Les conditions aux limites sont telles que: (121 pH(m)

- ./gft, J.(mr) p(r) dr
a) pour r ou z infini :
Jo désigne la fonction de Bessel de première espèce et
6ff O00
- 6ZZ: 6fZ: 0, d'ordre zéro.
(5) -
ur-u0-w-0, En introduisant I'expression (10) de la fonction intermé-
diaire dans les relations (8), nous obtenons pour les
b) à la surface du massif (z- O) :
contraintes et les déplacements les expressions suivantes
valables en tout point du massif semi-infini :

(6) orr(r,O) 0
-
(71 arr(r,Ol - p(r)
orr - []"(s,PzDr - srp,Dr) - (srA, sr\rlllÆl / (tt sr)
où p(r) désigne la distribution des pressions appliquée nor-
malement à la surface du massif, fonction nulle à l'exté- o00
- t- À(srpzDr - srp,D2)
rieu r du contact.
Nous avons ainsi posé de (3) à (7) les équations mathéma- Ls,p, (b-asl)n, -szp1 (b-asllnrlltl/(s,-sr)
tiques définissant le problème de mécanique considéré. : (srA,
6zz s,Ar) I (r, sr)
Introduisons une fonction intermédiaire A telle que: (1 3)
6rz- (C, Cr) I (s, sr) /d
orr: i,#.+#.a *t t (s' s')

-*(o#.+*. a*t oee: i: *ï,;':i,;,;';?ïi;'"'

où nous posons :

(8) ,6zz: (.+.f o

* #. #, Ai :Ja-s -ffisiz Jo(*r) pH(m) m dm ,

6rz-
- +(9+*13'
ôr àrz ,, a, * a uAl
à2a

1r-
Bi -T l"-tsiz Jo(mr) pH(m) dm ,
où les coefficients a, b, c, d sont définis comme suit : ,o

a- âr: (a,, arrl l(ârrâ:: ârr), (141

b - [a,r(a,3 + aool â,rar:] l(âr,â:r ârr), a- 1,2
(9) @
f
c- [a,r(a,, arrl + arraoo] l(ârrâ:r ârr), C;-J- | e-ms;z J,(mr) pH(m) m dm
o
,

d- (a,, alrl l(â,,ârr â,r),

Nous avons montré (Dahan, 1 975 ; Dahan et Zarka, 19771
que le problème est résolu si cette fonction s'écrit sous la D; - !J
r
.-ms;z J,(mr) pH(m) dm
forme :

(1 0) q(r,zl - Nous avons ainsi déterminé pour tout chargement p(r) le

comportement du sol correspondant.
1 f' ,. dm 3 Poinçons rigides axisymétriques
tpre-strnz - pre-srmzlpH(mUo (mr),,.'r*
t1s,_srlJ Dans le cas du contact du sol avec un poinçon rigide, la
o
fonction p introduite en (71 est à déterminer; nous la
-l savons simplement nulle à I'extérieur du cercle de contact
avec les notations suivantes: (r ) tJ.Les seules données du problème sont la forme du
poinçon f(r) et la profondeur de pénétration wo ou la force
de poinçonnement P. A I'intérieur du contact, au lieu de
I'hypothèse (7l., nous avons une condition aux limites en
déplacement :

(1 5) w(r,O) : wo -f(r) r(ro.

REVUE FRANÇAISE DE GEOTECHNIOUE NUMERO 7

A la su rface d u rnassif , les conditions a ux lim ites s'é c ri- ,. lo

vent, d'après les relations (13) avec z

su iva nte :
- 0, de la facon (211wo : ,"
Jof'(r)
(r"t
- r2l -t/2 6r.

w(r,O) r(ro Nous avons aussi la répartition des pressions sous le poin-
Jo çon et sa transformée en fonction de X:
(1 6)
roX'{t)
t(2 -ttz 61 r(
or=(r,Ol P(r) - - rlr| - r2l r o
Jo

(221
En supposant pH de la forme :

(171 pH(m) X(t) cos(mt) dt

PH(m) -J"I'"' sin(rnt) #
-lt"
Jo

Pour la force de pornconnement, nous avons la relation:

la condition (1 6) pour r identiquement vérifiée.
Pour r cette même condition ,lo
l'équation intégrale d'Abel (23l, P 2n I
xftl dt
:
- Jo
(18) q [ 'xtrt kz _t2) t/'dt: wo - f(r)
Jo
ll est aussi facile d'exprimer le tassement à I'extérieur du
dont la solution s'écrit : contact:

(1 9) r(t) - &two-t Jjt'(r) (t2

- r2l -t/z6vl. (241w(r,o)
- q Jf xttt k' - ('t -tt2 61 r)ro

Pour des profils de poinçons réguliers, sans discontinuité En conclusion, les formules (1 9) et (22]l permettent de
de plans tangents, p(r) reste fini sur le bord du cercle de déterminer la répartition des pressions p(r) sous le poin-
contact. Nous montrons (Dahan, 1 979) que cette condi- çon. L'expression (211 détermine le rayon du cercle de
tion est équivalente à :
contact en fonction du déplacement imposé wo et du profil
du poincon. Nous avons ainsi ramené le problème du poin-
(2Ol O. çon rigide à celui d'un chargement axisymétrique par une
X(ro) - distribution de pression p(r) sur un cercle de rayon ro:c'est
De l'inégalité (1 9), nous déduisons pour t fe problème qui a été traité au paragraphe précédent.
- ro, l'impor-
tante relation entre la profondeur wo de pénétration du Pour des poincons particuliers, le tableau suivant donne
poinçon et le rayon ro du cercle de contact: les principaux résultats (Dahan, 1979lr.

Poinçon

Plat Sphérique Conique

r lrl o r2t2R r cotg &

Prof 1I du trrclnçon

wo Wg ro2l R ro cotg &

*
Déplacement lmposé

P 4wo ro / e 8wo@F/3q 4tge *], t on

Force de 5rclnçon-
nernent

p(rl ?*o( ,"t ,rr-t'' 4t ,o2 - ,' lt'' -"o!gg arg ch la-
Presslon Qn QnR q r
aous le trrcinçon

p H(-l gg sin (m rol 4 ro -o- Gosrrre

lransformée de p
Qnm qrlRnt't^fgin
m ro #(1-cosmrsl
(r2)

REVUE FRANçAISE DE GEOTECHNIOUE NUMERO 7

4 Etude du coefficient q. 1rr et Éz sont les modules de cisaillement pour les.plans
d'isotropie et pour les plans perpendiculaires dans les
D'après les formules (22l. et (241, on remarque qu'à la sur- directions radiales et on sait eue Fr n'est pas une valeur
face du massif, la pression de poinçonnement p(r) et les indépendante.
tassements w peuvent être déduits de la solution isotrope
(Sneddon, 1965) ; il suffit simplement de remplacer le Nous rappelons que les coefficients élastiques a" sont défi-
coefficient q par le module 2(1 v2l/t du cas isotrope. nis, pour un matériau transversalement isotrop-e, à partir
de ces ( constantes techniques > par les relations :
Par conséquent, tous les résultats obtenus pour le massif
isotrope à partir de la formule de Boussinesq en ce qui
concerne la distribution des contraintes sous le poinçon et arr : -l- ,
le tassement à I'intérieur comme à I'extérieur du contact
El
sur la surface du massif, s'étendent au milieu à isotropie ârz : -'Vr
tra nsversa le. El ,

Ces résultats apportent aussi une justification aux calculs ar"

pratiques de tassement de sols faits en utilisant un module "Ê,- -1
(26)
expérimental obtenu à partir de lienfoncement d'une
plaque rigide sous une force donnée, comme le module de ^1
q11 :
--
Westergaard (module de réaction). Ce module joue le rôle E2
du coefficient q et tient compte de l'anisotropie. I
a44: :
Bien entendu, quand on s'écarte de la surface z - O, les Vz,
résultats ne sont plus les mêmes que pour un massif iso-
trope. Pour montrer I'influence de l'anisotropie, nous don- Les résultats (25) montrent l'existence d'une certaine ani-
nons ci-dessous pour le poincon plat et à titre de curiosité E,
pour le poinçon conique, les isovaleurs de la contrainte sotropie, le rapport des modules d'YoungÇest de I'ordre
normale orret celles du déplacement vertical w pour deux de 0,66, c'est-à-dire que le limon paraît plus raide dans le
sols anisotropes et deux sols isotropes tels que, dans cha- sens vertical oue dans le sens horizontal.
cun des deux cas, le coefficient q soit égal au rnodule 2
(t vzlle. Ces courbes ont été obtenues par un calcul La marne de Marseille a été prélevée à une cinquantaine
numérique simple à partir des formules explicites (13) de mètres de profondeur, et est en grande partie com-
posée de montmorillonite et d'illite. Les essais en com-
pour le cas anisotrope et des résultats de Sneddon (1946
pression simple ont donné les résultats suivants :
et 1948). Les valeurs utilisées pour les coefficients a;; des
sols anisotropes ont été mesurées au Laboratoire de Er: 1 '2 1Os Pa'
Mécanique des Solides sur des éprouvettes de lirnon pré-
levé sur le plateau de Palaiseau et sur des éprouvetteb Ez : 3.9 1O7 Pa,
d'une marne prélevée à Marseille. Ces deux sols présen- vt:Vz:0'
tent un sens d'anisotropie différent. (271
Des résultats présentés dans le tableau I et à partir des Er
relations entre la force de poinçonnement et la profondeu r tr': Z(t + uJ : 6 10? Pa,
de pénétration du poincon, nous pouvons déduire un
moyen expérimental très commode pour la mesure du Itz: 2,2 10? Pa,
coefficient q. et ainsi : q : 3,8 1O-8 Pa-r.
Nous avons cette fois une anisotropie significativement
5 Appfication pratique. différente, le rapport des modules EtlEr est pratiquement
Six éprouvettes de limon, de dimensions standard (hau- égal à 3, c'est-à-dire que la marne paraît plus raide dans le
teur :72 mm, diamètre :36 mm) ont été prélevées à Palai- sens horizontal que dans le sens vertical.
seau à une profondeur d'environ 1,2O m. Trois de ces Sur les courbes des figures 1 et 2, nous remarquons que
éprouvettes ont été taillées dans le sens vertical, deux les isovaleurs de la contrainte normale relative au limon de
dans le sens horizontal et une à 45 degrés. Palaiseau sont au-dessus de celles obtenues pour le cas
Sur ces éprouvettes, nous avons réalisé des essais de isotrope (trait discontinu), c'est-à-dire qu'en tout point du
compression simple et avons obtenu les résultats massif la contrainte o' anisotrope est supérieure à la con-
su iva nts : trainte o.. isotrope, sauf sur la surface où elles sont égales.
Nous avons le résultat contraire pour la marne de Mar-
seille. Cette différence s'explique par le sens de l'aniso-
Er _ 3'9 106 Pa, tropie des deux sols puisque le rapport E,/E, est inférieur à
1 pour le limon, alors qu'il est supérieur à 1 pour la marne.
E2- 5'9 106 Pa,
Pour les tassements, on voit sur les figures 3 et 4 que les
vl - O'10 différences varient en fonction du rapport d'anisotropie.
(25) vz - O,1 3 Nous avons ainsi de plus giands écarts avec le cas iso-
trope pour la marne que pour le limon.
Er
Itr - 2(1 +v )
_ 1,8 106 Pa, 6 Gonclusion
I
Nous avons donné dans cette étude un calcul analytique
Vz: 2'4 106 Pa' permettant d'obtenir la distribution des contraintes et des
déplacements en tout point d'un sol semi-infini transversa-
d'où la valeur de q - 3,6 1O-7 Pa-r. lement isotrope soumis à des chargements axisymétriques
sur sa surface libre. Cette méthode conduit pour la plupart
E, et E, sont les modules d'Young pour la direction paral-
des chargements classiques à des résultats explicites. dont
lèle aux plans d'isotropie du matériau et perpendiculaire à l'utilisation est aisée.
ce ux-ci.
Nous avons montré aussi que la répartition de la pression
v1 est le coefficient de Poisson qui caractérise la dilatation sous la fondation et le tassement en surface sont iden-
transversale dans le plan d'isotropie pour une compression
parallèle à ce plan.
tiques à un coefficient multiplicatif près, aux résultats
obtenus pour les sols isotropes à partir de la formule de
v2 est le coefficient de Poisson qui caractérise la dilatation Boussinesq. C:est une propriété très importante, qui per-
transversale dans le plan d'isotropie pour une compression met de généraliser aux milieux transversalement isotropes
de direction perpendiculaire à celui-ci. tous les résultats connus pour les sols isotropes.

REVUE FRANçAISE DE GEOTECHNIOUE NUMERO 7 I

 r. ou3runN 3norNHc3ro39 3cl 3srvÔNVUJ 3nA3H
6
(\ilr/d od) aleurou atulenuoe el ap srnaleAosl anfuuoc uoôuro4
- e '41
* otlt *otlt
\
\\\
\
t..
\\
\'.
N
\'---. N
o o
t I
eur e!\l uoul ! l
of,1zzs: eleurJou alu!etluoC odlzzD : eletu.rou elu!elluoC
(?tuZ/d = od) aleutou atu!enuoc el ap srnaleAosl ']eld uoôuro4 I '6lJ
not lt *otlt
so:! -
\ -cfl- \
'l
-
\
N
o N
I o
t
"urnW ftr uoul ! 1
of,1zzo : aleurrou aluleJluoC
 \-\f.
045 I
I
t Déplacement vertical : w/ro
0,030

\.tmon
\\
ô,6.0î____*_.. \
\
\
\
\
,\
\\\\
-\---1-
_ _____\__ \\
,___

-5__.\\\r,_
----\\\
0,090
ôigî--='\.\
_- - rr\\\\

rlfg+
rlro< 2

Fig. 3 Poinçon plat. lsovaleurs du déplacement vertical

Déplacement vertical : w/ro Déplacement vertical : w/re

Limon Marne

I I
o L
o
N N

\i-
\'\i
\l
rl(O+ r lro,+

Fig. 4 Poinçon conique. lsovaleurs de w

Références bibliog raphiques

BOUSSINESQ J. (1885) Applications des potenfiels à

l'étude de l'équilibre et du mouvement des so/ides élasti-
gues. Gauthiers-Villars, Paris.
DAHAN M. (1975)- Poinçonnement élastiquê par une DAHAN M., LUONG M.P. et MANDEL J. (1978) - Sur l'ani-
sphère d'Ltn massif semi-inf ini transyersa lement isotrope. sotropie des so/s. C.R. Acad. Sc. Paris, t.287, B, 179-181 .

Thèse 3" Cycle, Paris Vt. SNEDDON l.N. (1946) - Boussinesq's problems for a f lat-
DAHAN M. (1979) Poinçonnements axisymétriques ended cylinder. Proc. Camb. Phil. Soc., 42,29-39.
rigides sur un massif semFinf ini transyersa/ement isotrope. SNEDDON l.N. (1948) - Eoussinesq's problem for a rigid
A paraître, Journal de Mécanique Appliquée. cone. Proc. Gamb. Phil. Soc., 44,492-507.
DAHAN M. et ZARKA J. (1977) - E/astic contact between a SNEDDON l.N. (1965) The relation between load and
sphere and a semi-infinite transversely isotropic body. Int. penetration in the axisymrnetric Eoussinesg's problem for a
J. Sof ids Structures, 13, 229-238. punch of arbitrary profile. Int. J. Eng. Sci., 3, 47-57.

REVUE FRANçAISE DE GEOTECHNIOUE NUMERO 7