Réponse Dynamique Des Pieux en Zone Sismique

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 146

Table des Matires

Table des matires

AVANT PROPOS......1
TABLES DES MATIERES....... ..2
RESUME ..........6
ABSTRACT..........7
..8
LISTE DES FIGURES ...... ..............9
LISTE DES TABLEAUX . ..... 12
PRINCIPALES NOTATIONS......15

2
Table des matires

CHAPITRE 01 : INTRODUCTION
1.1Gnralits .......18
.2Problmatique.......18
1.3plan, du travail .........19

CHAPITRE 02 : GENERALITE
2.1introduction
2.2-Dfinition Et erminologie........20
2-2-1 dfinitions de la fondation profonde........20
2-3 historique.....22
2.4 Principe De fonctionnement ........25
2.5 Classification Des Pieux .....26
2.5.1Les pieux mis en place par refoulement du sol ....26
2.5.2- Les pieux mis en place sans refoulement du sol ...26
2.5.3- Les fondations injectes .......27
2.5.4 Diffrentes catgories de pieux ...28
2-6 Dtermination De La Charge Limite Dun Pieu Isole ....................30
2-7-Dimensionnement Des Pieux ...31
2-7-1- Dimensionnement par la mthode pressiomtrique.31
2-8- Ralisation Des Pieux ..33
2-9 comportements mcaniques des fondations profondes isoles sous chargement
monotone :
2.9.1 Paramtres influenant le comportement d'une fondation profonde isole ..34
2 .9.2 Indice de densit et contrainte de consolidation ...38
Table des matires

CHAPITRE 03 : METHODES DANALYSE DES PIEUX

3. pieu sous charge latrale statique


3-1Introduction 39
3-2 Mthode Au Module De Raction 39
3-1-2-1Principe gnral .......39
3-2-2 Rsolution dans le cas dun sol homogne ...41
3.2.3 Rsolution dans le cas dun sol non-homogne .....42
3.3-Expression Du Modle De Raction .42
3.4 Les Courbes de raction P-Y ....44
3.5 Paramtres ncessaires la construction des courbes rglementaires...57
3.6 Mthodes du continuum lastique..57
3-7 conclusion .58

CHAPITRE 04: METHODE SIMPLIFIEE DE CALCUL LA REPONSE


DYNAMIQUE DES PIEUX
4-1-Gnralit sur les sismes ..59
4-1-1 Naissance Du Sisme .59
4-1-2 Caractrisation des sismes ......................................................................................59
4-1-3 Effets particuliers des sismes ............................ ..60
4-1-4- Propagation Des Ondes ..61
4-2-Mthode Simplifie De Calcul La Rponse Dynamique Des Pieux ..... 63
4-2-1-Principe de la mthode ...63
4.3 Calcul sous dformation sismique impose . 65
4.3.1 Dtermination des efforts dans un pieu colonne ... 68
4.3.2 Dtermination des efforts dans un pieu flottant ..... 69
4.3.3-Application numrique et commentaires .................................................................70
4.4 Profil De Sol Constitu De Plusieurs Couches Homognes ......72
Table des matires

4.4.1 Calcul de la priode fondamentale quivalente........72


4.4.2 Dtermination de la dforme dun profil multicouche .......74
4.4.3 Exemples de priode pour la rgion de SAN FRANCISCO .......75
4.5 Dtermination des efforts dans un pieu en tenu compte la charge axiale .... 76

CHAPITRE 05 : SIMULATION NUMERIQUE


5.1 Simulation umrique....90
5.1 Gnralits :
5.1.1 ntroduction ....81
5.1.2 Choix dlment fini ..81
5.1.3 Dfinition de la gomtrie et les proprits des matriaux .. . 81
5.1.4 Caractristiques des matriaux .... 82
5.1.5 Maillage ....... 83
5.1.6 Conditions initiales ......... 83
5.2 Mthode simplifie .......... 84
5.2.1 Calcul les efforts internes par la mthode propose . .84
5.3 Calcul par la mthode numrique .........84
5.3.1 tape de calcul ...84
5.3.2 Modle du calcul ...85
5.4 Analyse numrique par logiciel ....91
5.4Etude aramtrique ....93
5.4.1 Module de dformation E=30.000kn/m2 ....93
5.4.2 Module de dformation E=50.000kn/m2 ...95
5.4.3 Module de dformation E=70.000kn/m2 97
5.4.4 Module de dformation E=90.000kn/m2 99
5.4.5 Variation du la vitesse V1=25.81m/s et E=3000 Kn/m2108
5.5.2 Variation du la vitesse V2=50m/s et E=11250 Kn/m2 .110
5.4.6 Variation du la vitesse V3=100m/s et E=45.000 Kn/m2112
5.4.7 Variation du la vitesse V4=150m/s et E=101250 Kn/m2 .114
Table des matires

5.4.8 Variation du la vitesse V5=200m/s et E=180000 Kn/m2 .116


5.5 Formulation propose .125

Conclusion Gnrale..128
Rfrences bibliographiques .....130

Annexe..132
Rsum

Rsum :
Les fondations profondes (pieux) sont trs utilises pour les ouvrages est devenu
une opration primordiale en Algrie vu les problmes causs par le phnomne de
sisme.
A ce jour, plusieurs tudes sont en cours de ralisation par les universits
algriennes pour trouver une solution ce phnomne, parmi ces tudes, la rponse
dynamique des pieux en zone sismique.
La prsente tude pour objectif de mieux apprhender le comportement des pieux
tout en prsentent dabord les divers mthodes disponibles dans la littrature pour
ltude les pieux bases sur les mthodes de dimensionnement et la mthode
simplifie.
Pour la partie simulation numrique on prsente un bref aperu sur le code des
lments finis PLAXIS, suivi dune tude de la rponse des pieux en zone sismique
Une tude paramtrique sur un modle choisi comme rfrence est a la fin ralis,
on fait varier chaque fois le paramtre comparer. La comparaison est prsente
par des courbes pour chaque paramtre

Mots cls : pieux mthode simplifie propos gotechnique plaxis rponse


dynamique
Rsum

Abstract:

The deep foundations (piles) are used for books has become a primary operation in
Algeria given the problems caused by the earthquake phenomenon.
To date, several studies are being carried out by Algerian universities to find a
solution to this phenomenon, among these studies, the dynamic response of piles in
seismic zone.

This study aims to better understand the behavior of piles while first present the
various methods available in the literature to study the piles based on the design
methods and simplified method.
For the numerical simulation part presents a brief overview of the finite element
code PLAXIS, followed by a study of the response of piles in earthquake zone
A parametric study on a model is chosen as reference was made late, we vary each
parameter to compare. The comparison is presented by curves for each parameter

Key words: Piles - Simplified Method Proposed - Geotechnical - PLAXIS -


Dynamic Response
Rsum

! "#$%& ' ( )*+ , -./012 345 6/ 7 8 ! 9


-6 : ; 5. < 0 => ? @ A. BC < ;? BCD E F 5G 5 @ !; $H ; @ AF . C I 2 J 5) ? =K ,0E
8+ ;
+1$;) L M2 N & ,O F 0 -6 & A. PF E0 N Q 0/; R )& < - 6 & S6 1 T & LD , O A. BC " .
P6
6; # -.8 #; U $1 <PLAXIS 8+ ; -6 & < 0 AF D 2 ? F).#; M ;A
!; : @ # A 8; ; D 5) A8 !; Q 0/2 # ? 5 V BCD ( BA 0W L2 JC XY6 , 5 0; A < AF U O L05

: ; 5. < 0 ^ 826 [[ \ 08 +1$; 85 ] [[ ;6 Z


Liste des figures

Liste des figures

Figure Titre
2.1 : Le domaine de fondation profonde

2.2 : Dfinition de la hauteur dencastrement D et De

2.3 : Schma typique de lapplication des pieux racines (Lizzi, 1982b)

2.4 : Frottement ngatif et le frottement positif

2.5 : Dfinition de A et P pour les pieux tubulaires ouvertes les palplanches et profils H

2.6 : Liaison Poteau-Pieu

2.7 : Courbe de chargement axial d'un pieu

2.8 : Dplacements (a) et dformations (b) du sol autour d'un pieu fonc dans le cas d'un
massifsableux (Robinski and Morrison (1964) et Vesic (1965))

2.9 : Zones de dformations du sol lors du fonage de modles de pieux dans du sable ; a)
Dplacements verticaux observs (b) zones de sol compact et de sol refoul autour
des pieux et (c) zones de dplacements horizontaux du sol (Shakirev et al, 1996)

2.10 (a) : Evolution du coefficient de pression latrale avec l'arrachement statique : sable
lche et sable dense (Puech et al. 1979)
3.1 : Reprsentation de modle de Winkler.
3.2 Courbe P-Y dans le cas de sollicitations de courte dure en tte dominantes
3.3 Courbe P-Y dans le cas de sollicitations accidentelles trs brves en tte dominantes
3.4 Loi de mobilisation de la raction tangentielle.
3.5Coefficients C1, C2 et C3 en fonction de langle de frottement interne daprs A.P.L.
3. 6 : Valeurs de K en fonction de la densite relative Dr et de langle de frottement interne
daprs A.P.L 1993
3.7 : Forme caractristique de la courbe P-Y daprs Norske Veritas (1977).
Liste des figures

4.1 : dformation dun pieu sous charges latrales


4.2: diffrents formes de spectres de rponse (Seed et al 1974)
4.3: la dformation dun pieu
4.4: diagramme du moment et effort tranchant
4.5: pieu supporte une masse en tte.
5.1 : la gomtrie du modle
5-2 : Maillage du massif
5-3 : Conditions initiales
5.4 : Variation des dplacements en fonction de l'acclration et de la vitesse
5.5 : Variation des dplacements en fonction de la zone
5.6 : dformation du pieu
5.7 : Dplacements totaux
5.8 : moment flchissent
5.9 : Variation des dplacements en fonction du la zone V=81.65m/s E =30.000 Kn/m2
5.10 : V=105.41m/s E =50.000 Kn/m2
5.11 V=124.72m/sE =70.000 Kn/m2
5.12 : Variation des dplacements en fonction du la zone V=141.42m/S E =90.000 Kn/m2
5.13 : Variation des dplacements en fonction du la vitesse a=0.15g
5.14 : Variation des dplacements en fonction du la vitesse a=0.20g
5.15 : Variation des dplacements en fonction du la vitesse a=0.25g
5.16 : Variation des dplacements en fonction du la vitesse a=0.30g
5.17 : Variation des moments en fonction du E a=0.15g
5.18 : Variation des moments en fonction du E a=0.20g
5.19 : Variation des moments en fonction du E a=0.25g
5.20 : Variation des moments en fonction du E a=0.30g
5.21: Variation des moments en fonction du E Et toutes les zones
5.22 : Variation des moments en fonction du la zone E=30.000 Kn/m2
5.23 : E=50.000 Kn/m2
5.24 : E=70.000 Kn/m2
5.25 : E=90.000 Kn/m2
Liste des figures

5.26: V=25.81m/S E =3.000 Kn/m2


5.27 : Variation des dplacements en fonction du la zone V=100m/S E =45.000 Kn/m2
5.28 : Variation des dplacements en fonction du la zone V=150m/S E =101250Kn/m2
5.29 : Variation des dplacements en fonction du la zone V=200m/S E =180000Kn/m2
5.30 : Variation des dplacements en fonction du la vitesse a=0.15g
5.31 : Variation des dplacements en fonction du la vitesse a=0.20g
5.32 : Variation des dplacements en fonction du la vitesse a=0.25g
5.33 : Variation des dplacements en fonction du la vitesse a=0.30g
5.34 : Variation moments en fonction du la zone E=3.000 Kn/m2
5.35 : Variation des moments en fonction du la zone V=50m/S E=11250 Kn/m2
5.36 : Variation des moments en fonction du la zone V=100m/S E=45.000
5.37 : Variation des moments en fonction du la zone V=150m/S E=101250 Kn/m2
5.38 : Variation des moments en fonction du la zone V=200m/S E=180.000 Kn/m2
5.39 : Variation des moments en fonction du la vitesse a=0.15g
5.40 : Variation des moments en fonction du la vitesse a=0.20g
5.41: Variation des moments en fonction du la vitesse a=0.25g
5.42 : Variation des moments en fonction du la vitesse a=0.30g
5.43 : Variation des moments en fonction du la vitesse a=0.15g (par la formule
proposes)
5.44 : Variation des moments en fonction du la vitesse a=0.20g (par la formule proposes
5.45 : Variation des moments en fonction du la vitesse a=0.25g (par la formule
proposes)
5.46: Variation des moments en fonction du la vitesse a=0.30g (par la formule proposes)
Liste des tableaux

Liste des tableaux


Tableau Titre
1 : valeur du coefficient A en fonction de la densit du sable
2 : valeur moyenne de ES pour des sols sans cohsion.
3 : facteur rhologique pour divers types de sols daprs Baguelin &al (1978
4 : facteur rhologique pour divers types de sols daprs le Fascicule 62
5 : Valeurs recommandes pour nh et pour des chargements statiques et cycliques dans
le sable
6: Valeurs recommandes pour, , et dans largile
7 : Paramtres ncessaires la construction des courbes P-y
5-1 proprits des couches de sol
5-2 proprits des pieux
5-3: les valeurs de dplacement pour V1=111.8 m/s et a =0.2g
5-4 : les valeurs de dplacement pour V1=111.8 m/s et a =0.3g

5-5 : les valeurs de dplacement pour V1=111.8 m/s et a =0.15g


5-6 : les valeurs de dplacement pour V1=111.8 m/s et a =0.25g
5-7 : les valeurs de dplacement pour V2=108.5 m/s et a =0.2g
5-8 : les valeurs de dplacement pour V2=108.5 m/s et a =0.3g
5-9 : les valeurs de dplacement pour V2=108.5 m/s et a =0.15g
5-10 : les valeurs de dplacement pour V2=108.5 m/s et a =0.25g
5-11 : les valeurs de dplacement pour V3=105.4 m/s et a =0.2g
5-12 : les valeurs de dplacement pour V3=105.4 m/s et a =0.3g
5-13 : les valeurs de dplacement pour V3=105.4 m/s et a =0.15g
5-14 : les valeurs de dplacement pour V3=105.4 m/s et a =0.25g
15- Rsultats Rcapitulatif
5.16 : variation du moment flchissant MS ; Mp
5.17 : les valeurs de dplacement pour V1=81.65 m/s et a =0.15g
Liste des tableaux

5.18 : les valeurs de dplacement pour V1=81.65 m/s et a =0.20g


5.19 : les valeurs de dplacement pour V1=81.65 m/s et a =0.25g
5.20 : les valeurs de dplacement pour V1=81.65 m/s et a =0.30g
5.21 : les valeurs de dplacement pour V1=105.41 m/s et a =0.15g
5.22 : les valeurs de dplacement pour V1=105.41 m/s et a =0.20g
5.23 : les valeurs de dplacement pour V1=105.41 m/s et a =0.25g
5.24 : les valeurs de dplacement pour V1=105.41 m/s et a =0.30g
5.25 : les valeurs de dplacement pour V3=124.72 m/s et a =0.15g
5.26 : les valeurs de dplacement pour V3=124.72 m/s et a =0.20g
5.27 : les valeurs de dplacement pour V3=124.72 m/s et a =0.25g
5.28 : les valeurs de dplacement pour V3=124.72 m/s et a =0.30g
5.29 : les valeurs de dplacement pour V4=141.42 m/s et a =0.15g
5.30 : les valeurs de dplacement pour V4=141.42 m/s et a =0.20g
5.31 : les valeurs de dplacement pour V4=141.42 m/s et a =0.25g
5.32 : les valeurs de dplacement pour V4=141.42 m/s et a =0.30g
5.33 : variation du moment flchissant MS ; Mp
5.34 : les valeurs de dplacement pour V1=25.81 m/s et a =0.15g
5.35 : les valeurs de dplacement pour V1=25.81 m/s et a =0.20g
5.36 : les valeurs de dplacement pour V1=25.81 m/s et a =0.25g
5.37 : les valeurs de dplacement pour V1=25.81 m/s et a =0.30g
5.38 : les valeurs de dplacement pour V2=50 m/s et a =0.15g
5.39 : les valeurs de dplacement pour V2=50 m/s et a =0.20g
5.40 : les valeurs de dplacement pour V2=50 m/s et a =0.25g
5.41 : les valeurs de dplacement pour V2=50 m/s et a =0.30g
5.42 : les valeurs de dplacement pour V3=100 m/s et a =0.15g
5.43 : les valeurs de dplacement pour V3=100 m/s et a =0.20g
5.44 : les valeurs de dplacement pour V3=100 m/s et a =0.25g
5.45 : les valeurs de dplacement pour V3=100 m/s et a =0.30g
5.46 : les valeurs de dplacement pour V4=150 m/s et a =0.15g
Liste des tableaux

5.47: les valeurs de dplacement pour V4=150 m/s et a =0.20g


5.48 : les valeurs de dplacement pour V4=150 m/s et a =0.25g
5.49 : les valeurs de dplacement pour V4=150 m/s et a =0.30g
5.50 : les valeurs de dplacement pour V5=200 m/s et a =0.15g
5.51 : les valeurs de dplacement pour V5=200 m/s et a =0.20g
5.52 : les valeurs de dplacement pour V5=200 m/s et a =0.25g
5.53 : les valeurs de dplacement pour V5=200 m/s et a =0.30g
5.54 : variation du moment flchissant MS ; Mp
5.55 : les valeurs de moments flchissant corriges par les formules proposes
Liste des tableaux
Principales notations

Principales notations
Lettres latines :

pp : coefficient rducteur de section de leffort de pointe


ps coefficient rducteur de section de leffort de frottement latral
A : aire de la section droite
P : primtre de la section du pieu
qpu : rsistance limite de pointe
qSU : frottement latral unitaire limite dans couche i
ei : paisseur de la couche i
h : hauteur dancrage

P1e*: pression limite nette quivalente

P1*(Z) : pression limite nette la profondeur z

d : pris en gnral gal 0

D : hauteur dencastrement relle du pieu


QS lim : le frottement latral limite (avec fs ( qs) le frottement latral unitaire limite)
QP lim : la rsistance en pointe limite (avec qp la rsistance en pointe unitaire limite).
P : pression (N/m2).
Kh (z) : le coefficient de la raction une profondeur z (N/m3).
Y (z) : le dplacement du pieu pour une profondeur z (m).
P : la raction d sol (N/m)
ES : le module de raction du sol pour une profondeur z (N/m2).
B : le diamtre de pieu (m).
EP : module dlasticit du pieu.
IP : moment dinertie de la section transversale du pieu
Z : profondeur dans le sol.
A : un coefficient adimensionnel fonction de la densit du massif de sable (tableau 1.2)
Z : la profondeur (m).
Principales notations

Y : le dplacement latral (m).


PUS : raction ultime du sol en surface ;
PUS : raction ultime du sol pour les couches profondes ;
B : le diamtre du pieu ;
CU : cohsion non draine du sol (Kpa)
J : constante empirique dont les valeurs comprises entre 0.25 et 0.5 ont t dtermines par des
essais sur site
ZC : profondeur critique
KS : module de raction latrale dans un sable
KC : module de raction latrale dans une argile
Pd : rsistance latrale de dimensionnement
Plim : valeur asymptotique de lhyperbole pour y >
k1 : pente initiale de la courbe
P0 : pression effective (aprs excavation)
CU : contrainte de cisaillement caractristique non draine reprsentative des conditions de
chargements considres.
Pd : PdS pour un chargement statique
Pd : Pdc pour un chargement cyclique
triaxial statique non drain sur un chantillon de sol non remani.
K : la raideur initiale de la courbe P-y ;
Pu : la rsistance ultime du sol ;
Y : le dplacement latral
K0 : coefficient des terres au repos ;
Ka : coefficient des terres actives
kh : module initial de raction la profondeur z,
Pmax: pression limite la profondeur z,
Y : poids volumique du sol,
D : diamtre du pieu,
n : acclration centrifuge.

-
Principales notations

Lettres grecques

: Le poids volumique du sable.


: Un coefficient rhologique caractrisant le sol (tableau 4)
C : Valeur moiti de la dformation correspondant au maximum de la contrainte dviatorique
dans un essai triaxial non drain.
mc : Coefficient de scurit caractristique du matriau

: Plin/Pd = rapport de rupture, toujours suprieur 1


: coefficient dpendant du sol et des conditions de chargement (tableau 6)
d : Angle de frottement pour le projet
: Angle de frottement rel
: Coefficient empirique (tableau 5)
C : Dformation verticale la moiti de la variation des contraintes principales dun essai
: Langle de frottement (degr) ;
Principales notations
CHAPITRE 01 :

Introduction
Chapitre 01 : Introduction

1-Introduction :
1.1 Gnralit :
On fait appel des fondations sur pieux lorsque le sol de fondation n'a pas les caractristiques
suffisantes pour supporter le poids de la structure et les charges transmises par la structure au sol
de fondation.
Les pieux et micro pieux sont utiliss depuis longtemps dans le domaine du gnie civil, leurs
mises en place est trs varie et va de la fabrication en place a l'inclusion ou l'insertion dans le
sol par battage, verinage ou vibrofonage d'lments prfabriqus (tubes, pieux en bois, bton
..ect)
En Algrie, l'ingnieur gotechnicien conoit ce type de fondation en se basent sur les rgles du
fascicule 62(titre .V. cahier des clauses technique gnrales)
Le calcul et la vrification sous chargement sismique sont peu utiliss comme tout autre type de
fondation, les fondations profondes peuvent tre soumises a des chargements trs varies et
complexe : monotones et cycliques, verticaux et horizontaux.
Alors se pose le problme de leur dimensionnement afin de rsister et de bien se comporter par
rapport ces sollicitations diverses.
1-2 Problmatique:
Les fondations sur pieux sont trs utilises pour les ouvrages tels que les btiments, tours les
ponts.
En Algrie cependant leur utilisation est restreinte pour les ponts, peu les btiments sont fondes
sur pieux
Pourtant dans les sites moyenne et forte sismicit, les pieux peuvent apporter la solution
certains problmes tels que la liqufaction des sols et les sols ayant un comportement dynamique
faible. Dautre part il y a en l'Algrie:
Absence dune rglementation pour le dimensionnement et la vrification de ce type de
fondation
L'Algrie se trouvant dans une zone de sismicit importante, il est impratif de savoir
calculer et dimensionner les pieux sous chargement sismique.
Chapitre 01 : Introduction

De ce fait le choix du thme concernant la rponse dynamique des pieux en zone sismique vise
proposer un outil de calcul et de vrification de ce type de fondation pour les ingnieurs, ainsi
une contribution dans ce domaine en vue dune large vulgarisation en attendant de voir une
rglementation apparaitre. Il faut se noter que ce mmoire est le premier a luniversit de Batna et
parmi les rares mmoires proposes a travers les universits algriennes.
1-3 Plan de travail :
Pour atteindre l'objectif vis, le travail est divis en quatre chapitres en plus dune introduction et
des conclusions.
Le chapitre 1 concerne une introduction gnrale qui situe le contexte du sujet, son importance
et le plan de travail mis en uvre
Le chapitre 2 concerne ltude bibliographique qui est compos de deux parties
La 1ere partie concerne des gnralits sur les pieux et leurs utilisation ; (dfinitions
classifications dimensionnements ..)
La 2eme partie concerne tude du comportement des pieux sous charges verticales et
latrales de type monotones.

Le chapitre 3 concerne les mthodes danalyses des pieux


Quant au chapitre 4 : il prsente une mthode simplifie de calcul de la rponse dynamique des
pieux en se basant sur la dforme produite par un mouvement de sol
Cette mthode simplifie consiste dfinir l'quation de la dforme du profil de sol homogne
ou multicouche correspondant au mode fondamental d'oscillations libres de la colonne de sol et
dterminer efforts M,T sous une telle dformation sismique unipose, en faisant l'hypothse que
le pieu est articul en pied et encastr en tte.
Annexe est relatif a la modlisation dun pieu isol soumis a un mouvement de sol dform.
Le Chapitre 05: simulation numrique : Ce chapitre a pour objectif de faire une comparaison
entre les rsultats de la mthode simplifie avec ceux donns par la simulation numrique. Pour
pouvoir tirer des conclusions adquats on fait varier plusieurs paramtres : sismicit de la zone,
caractristiques du sol, (E, densit) vitesse de cisaillement.
Chapitre 02 : Gnralit

CHAPITRE 02 :

Gnralit
Chapitre 02 : Gnralit

2-Gnralits :

2.1 Introduction :

Les fondations profondes permettent de reporter les charges dun ouvrage au niveau des couches
situes entre profondeur. Elles sont en gnral utilises quand la rsistance des couches des
terrains superficielles nest pas suffisante pour supporter les charges transmises par une fondation
superficielle ou que les tassements induits par ce type sont trop importants
Les pieux et les micropieux (ou fondations profondes) permettent de profiter des couches
rsistantes, sur des sites prsentant des caractristiques mcaniques insuffisantes, pour la
reprise des efforts transmis par les ouvrages. Ces efforts sont repris par linclusion sous forme
de frottement latral, mobilis le long du ft du pieu, et de la rsistance en pointe.

Les pieux et les micropieux se diffrencient essentiellement par leur dimension et leur utilisation.

2.2 Dfinition Et Terminologie :


2.2.1 Dfinitions de la fondation profonde :
Une fondation profonde est caractrise par la manire dont le sol est sollicit pour rsister aux
charges appliques.

- rsistance en pointe

- par frottement latral

- rsistance de pointe et frottement latral (cas courant)

Ses dimensions sont dfinies par:

- D : Longueur de fondation enterre dans le sol

- B : largeur de la fondation ou diamtre

Au-del de D/B> 6, et D> 3, nous sommes dans le domaine des fondations profondes fig (2.1).

Figure2.1 : Le domaine de fondation profonde


Chapitre 02 : Gnralit

Dune manire gnrale, les fondations profondes sont souvent dsignes par le terme de pieu
Un pieu est une fondation lance qui reporte les charges de la structure sur des couches de
terrain de caractristiques mcaniques suffisantes pour viter la rupture du sol et limiter les
dplacements des valeurs trs faibles. Le mot pieu dsigne aussi bien les pieux, les puits et
les barrettes.

On dsigne par pieu, une fondation profonde ralise mcaniquement et par puits une fondation
profonde creuse la main sous la protection dun blindage.
Une barrette est un pieu for de section allonge ou composite (en T ou en croix par exemple)
Les trois parties principales dun pieu sont

La tte,

La pointe,

Le ft compris entre la tte et la pointe.


La longueur dancrage h est la longueur de pntration du pieu dans les couches de terrain
rsistantes.
Dun point de vue mcanique on distingue la longueur D du pieu de la hauteur dencastrement
mcanique De.

Cette valeur de De tient compte du fait que les caractristiques mcaniques de la couche
dancrage sont nettement suprieures celles des sols de couverture traverss par le pieu fig(2.2).

Figure2.2 : Dfinition de la hauteur dencastrement D et De

On considre (Fascicule 62 titre V) quun lment de fondation est de type profond lorsque sa
hauteur dencastrement relatif De/B est suprieure 5.
Chapitre 02 : Gnralit

2.3 Historiques:

La conduite de porte des pieux de structure dappui est lun des premiers exemples de lart et la
science de lingnieur. Dans grande-bretabne il y a de nombreux exemples de bois dans
lempilage de travaux et de pont au bord de la rivire tablissements construits par les Romains.

Dans lpoque mdivale, des pieux de chne et laulne one t utiliss dans les fondations des
grands monastres dans le fenlands. DEast anglia. En chine, lempilage de bois a t utilis par
les btisseurs de ponts de la dynastie des han (200avant JC lan 200)

La capacit de charge des pieux de bois est limite par la circonfrence du bois et de la capacit
du matriau rsister conduite par le marteau, sans avoir subir des dommages en raison de
scission ou de fragmentation .ainsi primitive rgles doivent avoir t tablies ds les premiers
jours de lempilage par laquelle la charge admissible sur un tas a t dtermine partir de sa
rsistance la conduite par un marteau de connatre le poids et connu avec une hauteur de chute.

De connaissances a t accumule en ce qui concerne galement la durabilit des pieux de


diffrentes espces de bois, et les mesures prises pour prvenir la dgradation par la carbonisation
du bois ou en maonnerie de construction des radeaux sur pilotis ttes coupes au dessous du
niveau de leau

Bois, en raison de sa force combine la lgret et la facilit de coupe et de manutention, reste


le seul matriau utilis pour lempilage jusqu relativement rcemment. Il a t remplac par le
bton et lacier seulement parce que ces nouveaux matriaux pourraient tre fabriqus dans les
units qui ont t capable de soutenir la compression, flexion et traction.

Bton, en particulier, est adaptable aux in-situ des modes de construction qui facilit
linstallation de bases entasss dans des trous percs dans les situations ou le bruit, les vibrations
et le haut-sol a du tre vit.

Bton arm, qui a t labor comme un moyen structurel la fin du XIXe et du dbut du XXe
sicle, largement remplac le bois de grande capacit pour les travaux de battage sur la terre. il
pourrait tre prfabriqu dans diverses formes structurelle de faon convenir aux imposes
chargement et de conditions de terrain, et sa durabilit a t satisfaisante pour la plupart des sols
et des conditions dimmersion. Le partielle le remplacement des conduits en bton par des pieux
de nombreuses formes de coul in situ des pieux a t du davantage llaboration de hautement
efficace des machines pour le forage des trous de forage.
Chapitre 02 : Gnralit

Lacier : a t utilis dans une mesure croissante, pour lempilage en raison de sa facilit de
fabrication et de manutention et de sa capacit rsister difficile de conduire. Problmes de
corrosion des structures maritimes ont t surmonts par lintroduction de revtement et durable
protection cathodique.
La technique des micropieux a t dveloppe en 1952 par l'entreprise Fondedile sous la direction
technique de F. Lizzi (1982a). Ces premiers pieux de petits diamtres (B = 100 mm) taient fors
et scells au terrain par un mortier (fig2.3). Ils ont t initialement utiliss en Italie pour les
travaux de reprise en sous-uvre de btiments anciens et taient connus sous lappellation de
'pieux racines' (ou "pali radice").
En l'espace d'un demi-sicle, la technique des micropieux est pratique travers le monde et
connue sous multiples termes locaux : pali radice ou micropali (en italien), pieux racines, pieux
aiguilles, minipieux ou micropieux (en franais), minipiles, micropiles, pin piles, root piles ou
needle piles (en anglais), mikropfahl, verpresspfahle ou wurzelpfahle (en allemand) ou estaca raiz
(en portugais).
L'essor de cette technique dans les annes 50 en Italie rsulte du besoin de la restauration des
btiments historiques endommags par la deuxime guerre mondiale. Les micropieux prsentent
ainsi une solution technique simple et avantageuse pour les travaux de reprise en sous-uvre :
La lgret et le faible encombrement des matriels utiliss s'adaptent bien aux conditions
d'intervention. En outre, leur mise en place n'entrane que de faibles et minimales perturbations
l'quilibre de la structure, du sol et de l'environnement.
Le dveloppement de ce systme de fondations profondes dcoule de la diversit de ses
applications et de ses multiples avantages. En effet, la fin des annes 50, la pratique des
micropieux s'tend la stabilisation des pentes, aux renforcements des sols et la protection des
ouvrages souterrains. Les micropieux y sont utiliss en rseaux enchevtrs ("rseaux de pali
radice"). Au niveau technique, on observe l'apparition des micropieux injects sous forte pression
qui ont permis d'obtenir des portances plus leves (Lizzi, 1982b).

Figure2.3 : Schma typique de lapplication des pieux racines (Lizzi, 1982b)


Chapitre 02 : Gnralit

L'entreprise Fondedile introduit la technique du pieu racine en Grande Bretagne en 1962 pour la
restauration des btiments anciens. En 1965, des micropieux similaires aux 'pieux racines' sont
utiliss en Allemagne pour la ralisation des ouvrages souterrains en milieu urbain. Les pieux
racines sont galement utiliss lors de l'excution des travaux du mtro de Milan en Italie.
Les entrepreneurs europens commencent dvelopper leurs propres techniques au dbut des
annes 60. Leurs filiales mondiales participent, par la suite, l'exportation de la technique vers
les autres continents.
En France, la pratique des systmes de fondations par micropieux concerne plusieurs domaines
d'application. Si l'usage de cette technique est moins rpandu en France qu'en Italie, les
recherches menes en France, quant elles, ont donn un lan la matrise de cette technique sur
une chelle internationale. Ds le dbut des annes soixante, le dveloppement de la mise en
place des micropieux sinspire de la mthode dinjection slective du coulis des tirants
dancrages. Cette nouvelle mthode dinjection de coulis acquiert aux micropieux une portance
plus importante et une meilleure adaptation aux contraintes du sol.
Ce n'est qu'en 1973 que l'entreprise Fondedile introduit, en Amrique du Nord, le pieu racine
pour l'excution de travaux de reprise en sous-oeuvre dans les villes de Boston et New York. Le
premier rseau de 'pieux racines' est ralis aux Etats-Unis en 1975 pour la stabilisation de la
fondation d'une pile d'un pont sur l'autoroute 55 au Missouri. L'expansion de l'application de cette
technique aux Etats-Unis dbute vers le milieu des annes 80.
2.4 Principe De Fonctionnement :

Un pieu transmet au sol les charges quil supporte:

- par lappui de sa base sur le sol rsistant (effort de pointe not QP)

- par le frottement latral entre le sol et le pieu (effort de frottement latral not QS)

Leffort de pointe est proportionnel :

- section de la base du pieu

- la rsistance du substratum

Leffort de frottement latral est proportionnel :

- la surface de contact entre le pieu et le sol

- au coefficient de frottement pieu-sol (rugosit du pieu, pression latrale, coefficient de


frottement interne du sol)
Chapitre 02 : Gnralit

Le frottement latral du pieu nest mobilisable que sil y a dplacement relatif entre le pieu et le
sol. Si le pieu a tendance senfoncer dans un sol stable, le frottement sol-pieu gnre un effort
vertical ascendant (frottement positif)

Si au contraire, le pieu tant immobile, le sol tendance tasser, le frottement sol-pieu est
ngatif. Cela pour consquence de surcharger le pieu. Pour remdier ce problme (couches
compressibles, remblais rcents non stabiliss), on chemisera le pieu par un tubage afin de
diminuer leffet du frottement ngatif fig(2.4).

Figure2.4 : Frottement ngatif et le frottement positif

2.5 Classification Des Pieux :

Le pieu est fonction de son mode dexcution, donc on distingue trois grandes classes de
fondations profondes:

2.5.1 Les pieux mis en place par refoulement du sol :

Ils concernent les pieux battus (catgories 1 et 2) et les pieux foncs (catgorie 5) :

Ce sont des pieux soit faonns lavance soit tube battu excuts en place. Pour les premiers il sagit
essentiellement de pieux en mtal et de pieux prfabriqus en bton arm, pour les seconds de pieux battus
mouls.

Les pieux mtalliques sont gnralement sous forme de tube ou en forme de H. Les tubes peuvent
tre ouverts ou ferms leur base. Les pieux mtalliques sont mis en oeuvre par battage ou par
vibration.
Les pieux en bton arm sont fabriqus sur des aires proches du chantier. Ils sont mis en uvre
par battage ou par vibration.
Chapitre 02 : Gnralit

Lexcution des pieux battus mouls consiste dabord battre un tube muni sa base dune plaque
mtallique dans le sol, mettre, si ncessaire, en place la cage darmatures, puis remplir le tube de bton
pendant son extraction.

Leur mise en place se fait par vrinage battage ou vibro-fonage. Ils repoussent le sol et le
compriment ce qui gnre un bon frottement latral.

2.5.2- Les pieux mis en place sans refoulement du sol :

Ils concernent les pieux fors et les puits (catgories 3 et 4). Leur mise en place se fait par
substitution. Ce qui pour effet de remanier le sol et de le dcomprimer. Le frottement latral est
donc diminu, sauf pour certains types de mise en uvre (pieux excuts la tarire creuse, ou
visss mouls)

Leur excution ncessite un forage pralable excut dans le sol avec les outils appropris avec
ou sans protection dun tubage ou de boue permettant dassurer la stabilit des parois du forage.
Aprs mise en place, si ncessaire, de la cage darmatures, le pieu est btonn en utilisant une
colonne de btonnage, selon la technique du tube plongeur qui descend jusqu la base du pieu.

On fera une place part aux pieux la tarire creuse qui ont connu un trs fort dveloppement
ces dix dernires annes et qui sont trs utilises dans les fondations de btiment. Le principe
consiste visser dans le sol une tarire axe creux sur une longueur totale au moins gale la
longueur du pieu raliser, puis lextraire du sol sans dvisser pendant que, simultanment, le
bton est inject par laxe creux de la tarire. On distingue trois types de tarire creuse : les
tarires creuses sans enregistrement des paramtres de forage et de btonnage (type 2), les
tarires creuses (type 2) avec enregistrement des paramtres de forage et de btonnage
(profondeur, pression de bton, quantit de bton), et ces mmes tarires quipes en plus dun
tube de btonnage tlescopable rtract pendant la perforation et plongeant dans le bton pendant
le btonnage (type 3) Ces pieux ne peuvent, videmment, tre arms quaprs lopration de
btonnage. La mise en place des armatures sous leur propre poids ou par vibration devient
dlicate pour des pieux dune longueur suprieure 22 -25m.
Au-del, des techniques de bton arm de fibres en acier ont t dveloppes et permettent
dinjecter directement dans le sol les fibres mlanges au bton. Grce lamlioration
spectaculaire du matriel de forage on peut excuter des pieux de 600 2200mm de diamtre
Chapitre 02 : Gnralit

jusqu des profondeurs allant jusqu 30-35m. Actuellement les pieux la tarire creuse sont
parmi les moins chers sur le march franais.

2.5.3- Les fondations injectes :

Elles concernent les parois moules, les barrettes.

Dautres types de fondations existent, mais ils sont plutt classs comme renforcement de sol:

- Les colonnes de Col-mix (mlange de sol en place et de liant laide de deux tarires
creuses)

- Les colonnes dejet-grouting (forage puis injection de liant haute pression en remontant)

- Les colonnes ballastes (fonage dun vibreur puis injection lair comprim de ballast
(40/80 mm) en remontant)

- Les picots de sable (fonage dun tube de petit diamtre, puis introduction de sable vibro-
compact

- Les plots pilonns (excavation la pelle mcanique, puis introduction de matriaux


granulaires sains subissant un pilonnage intensif)

2.5.4 Diffrentes catgories de pieux :

Les documents rglementaires classent les pieux selon les catgories ci-dessous

1 - Pieux faonns lavance

Battu prfabriqu

Mtal battu

Tubulaire prcontraint

Battu enrob

Battu ou vibrofonc, inject haute pression

2 - Pieux tube battu excut en place

Battu pilonn

Battu moul

3 - Pieux fors

For simple
Chapitre 02 : Gnralit

For tub

For boue

Tarire creuse (type 1 type3 selon la technologie utilise)

Viss moul

Inject haute pression

4 Puits :

Les puits sont des fondations creuses la main. De section circulaire, ils doivent avoir un
diamtre suprieur 1,20m ; de section quelconque ils doivent avoir une largeur minimale de
0,80m et une section minimale de 1,1m2.
Les parois du puits sont soutenues par un blindage qui peut-tre rcupr ou abandonn. Aprs
curage du fond du puits, le forage est btonn sec.

5 Pieux foncs

Bton fonc

Mtal fonc

6 - Micro-pieux de diamtre infrieur 250 mm

Un micropieu est un pieu for de diamtre infrieur 250 mm, gnralement entre 76 et 200 mm,

qui comporte des armatures centrales scelles dans un coulis de ciment. Ils travaillent seulement

en compression ou en traction. A l'origine ils taient destins reprendre des charges de 150

250 kN, actuellement ils peuvent atteindre 1000 kN. Lemploi des micropieux a t introduit en

France par lentreprise FONDEDILE (entreprise internationale italienne) dans les annes 1960,

puis repris et dvelopp avec les techniques propres aux entreprises franaises et europennes.

On classe en France, (DTU 13.12 et fasc. 62-titre V), suivant le matriel de forage et les

techniques dinjection 4 types de micropieux.


Chapitre 02 : Gnralit

Type I: pieu for tub, rempli de mortier (micropieux de type FONDEDILE), nest plus
actuellement utilis en France.

Type II: pieu for, quip d'armatures, scell au coulis de ciment ou au mortier par
graviter au moyen dun tube plongeur.

Type III : pieu for, quip d'armatures et d'un systme d'injection qui est un tube
manchettes mis en place dans un coulis de gaine. L'injection du coulis de ciment est faite
en tte une pression gale ou suprieure 1MPa, elle est globale et unitaire (IGU).

Type IV : pieu for, quip d'armatures et d'un systme d'injection qui est un tube
manchettes mis en place dans un coulis de gaine. L'injection du coulis de ciment est faite,
chaque niveau de manchettes, avec un obturateur simple ou double une pression gale
ou suprieure 1 MPa. L'injection est rptitive et slective (IRS)

A ct de ces micropieux, dont lexcution est "normalise" on utilise dautres types de micropieux, en

particulier les micropieux autoforeurs (de type Ischebeck par exemple)

2.6 Dtermination De La Charge Limite Dun Pieu Isole :


Considrons un pieu isol soumis une charge verticale. Le pieu traverse diffrentes couches de
sol de qualit plus ou moins bonnes pour sancrer dans une couche de sol aux caractristiques
mcaniques favorables.
Cette couche sappelle couche dancrage ou substratum rsistant.
La charge limite du pieu QU, est obtenue en additionnant la charge limite de pointe QPU, qui
correspond au poinonnement du sol sous la base du pieu et la charge limite QSU mobilisable par
le frottement latral entre le sol et le pieu :

QU = Q PU + Q SU (2.1)

La charge limite de pointe est donne par :

QPU = p P . A.q PU (2.2)

La charge limite de frottement est donne par :

QSU = p S . A.q SU .ei (2.3)

Avec:
- pp : coefficient rducteur de section de leffort de pointe
Chapitre 02 : Gnralit

- ps coefficient rducteur de section de leffort de frottement latral


- A : aire de la section droite
- P : primtre de la section du pieu
- qpu : rsistance limite de pointe
- qSU : frottement latral unitaire limite dans couche i
- ei : paisseur de la couche i
- h : hauteur dancrage
Nota : La dtermination de A et de P ne pose pas de problme particulier pour les pieux section
pleine ou pour les pieux tubulaires fermes. Pour les autres sections, on se rfrera au tableau ci-
dessous fig (2.5) :

Figure2.5 : Dfinition de A et P pour les pieux tubulaires ouvertes les palplanches et profils H

2.7-Dimensionnement Des Pieux :

Comme pour les fondations superficielles, le dimensionnement des pieux se ralise partir des
essais de laboratoires, de lessai pntromtrique ou de lessai pressiomtrique.

Actuellement la mthode pressiomtrique donne de bons rsultats quel que soit le type de sol.
Elle est prsente ci-aprs.

2.7.1- Dimensionnement par la mthode pressiomtrique

Dtermination de la contrainte limite de pointe qpu :

Le fascicule 62 titre V dfinit une hauteur dencastrement quivalente De, selon la formule
suivante:

1 D *
De = P1 (Z ).dZ (2.4)
P *1e d
Chapitre 02 : Gnralit

Avec :

- P1e*: pression limite nette quivalente

- P1*(Z) : pression limite nette la profondeur z

- d : pris en gnral gal 0

- D : hauteur dencastrement relle du pieu

La valeur de la pression limite nette quivalente est donne par la formule :

* 1 D +3a *
P1e = P1 (Z ).dZ (2.5)
b + 3a d b

Avec:

B
a: max ( , 0.50m ).
2

b: min(a, h)

h: ancrage dans la couche o se situe la pointe du pieu

La contrainte limite de pointe qpu est donne par la formule:

qpu =kp.P*1e (2.6)

KP : est appel coefficient de portance. Il traduit la proportionnalit entre la contrainte limite de


pointe et la pression limite nette quivalente mesure au pressiomtre.
Chapitre 02 : Gnralit

2.8 Ralisation Des Pieux :

Principes de conception :

La jonction entre un pieu isole et un poteau est toujours assure par une semelle de liaison en
bton arm.

Figure2.6 : Liaison Poteau-Pieu


Chapitre 02 : Gnralit

2.9 COMPORTEMENT MECANIQUE DES FONDATIONS PROFONDES ISOLEES


SOUS CHARGEMENT MONOTONE :
Lorsqu'un pieu est charg axialement en tte, la rupture, la charge limite (QT lim) est quilibre
par les ractions limites du sol (fig 2.8) soit :
Ql lim ( QT lim) = Qs lim + QP lim (2.7)
Avec :
QS lim : le frottement latral limite (avec fs ( qs) le frottement latral unitaire limite)
QP lim : la rsistance en pointe limite (avec qp la rsistance en pointe unitaire limite).
Lexpression (2.7) requiert une bonne connaissance du comportement sol-structure pour la
dtermination du frottement latral limite QS lim et de la rsistance en pointe limite QP lim et donc la
prise en compte de trs nombreux paramtres influenant le comportement de ces fondations,
notamment le comportement du sol au voisinage du pieu au cours de l'installation et du
chargement. La fig(2.8) montre galement une courbe typique de chargement d'un pieu isol. La
charge limite, sur la courbe effort-dplacement, est gnralement dfinie comme celle
correspondant un dplacement en tte gale d/10 (d tant le diamtre du pieu) (AFNOR,
1994).

Figure 2.7 : Courbe de chargement axial d'un pieu

Paramtres influenant le comportement d'une fondation profonde isole :

2 .9.1 Mode d'installation des inclusions :


Il existe diffrentes mthodes d'installation des pieux. En France, on peut distinguer cinq
catgories de fondations profondes (Fascicule n62 Titre V, 1993) : pieux faonns l'avance (le
plus souvent battus), pieux tube battu excuts en place, pieux fors, puits et micropieux.
On prsente tout dabord linfluence des diffrentes mthodes d'installation sur le sol au
voisinage du pieu (dformations, dplacements, densification) ainsi que sur la capacit portante
du pieu.
Chapitre 02 : Gnralit

Robinski et Morrison (1964) ont effectu une srie d'essais sur des modles de pieux foncs dans
un sable lche et ont visualis les dplacements du sol autour du modle par radiographie. Ils ont
constat des dplacements importants et une densification du sol sous la pointe qui sont suivis par
un dplacement du sol vers le bas (direction de fonage) proximit du ft (fig 2.9(a)). Le
mouvement des grains est quantifiable jusqu' une distance de 3 ou 4 diamtres de pieu dans la
direction latrale et de 2,5 ou 3,5 diamtres de pieu sous de la pointe.
Vesic (1965) a transform ces dplacements en dformations et on peut ainsi observer des
dformations de compression pour le sol au dessous de la pointe et des dformations de traction
pour le sol situ au dessus (fig 2.9 (b)).
Dans le but de mieux comprendre la cinmatique des mouvements du sol, Shakhirev et al. (1996)
ont prsent une tude sur le comportement d'un massif sableux lors du fonage d'un pieu, qui a
permis darriver aux conclusions suivantes :
1. proximit du pieu de section constante, les dplacements verticaux du sol sont toujours dirigs
vers le bas et, une certaine distance du ft, le mouvement est plutt ascendant (fig 2.10 (a)) ;
2. le phnomne (1) se traduit par la formation de deux zones (fig 2.10 (b)). La zone (1) est une zone
compacte proximit immdiate du pieu qui augmente avec la profondeur. Elle est entoure par
la zone (2), qui est la zone d'inversion des dplacements verticaux et de refoulement du sol. De
plus, les dplacements horizontaux du sol (fig 2.10(c)) entranent galement la formation d'une
zone comprime qui est semblable la zone comprime dans la direction verticale.
3. l'tat de contrainte est tudi par l'intermdiaire de l'examen des zones de dformation, qui permet
d'observer que lors du fonage, des contraintes verticales et horizontales de compression
apparaissent tant au niveau du ft qu'au dessous de la pointe.
De plus, dans le cas d'un massif sableux, Mokrani (1991) a mesur lors du passage de la pointe
pendant la phase de battage de pieux l'volution des contraintes verticales ( v), radiales ( r) et
orthoradiales ( ) 100 cm de profondeur 3,5 diamtres du pieu (figure 2.6). Les contraintes
( v
et ) augmentent l'approche de la pointe en passant par un maximum environ 15 20 cm
r

avant que la pointe n'atteigne le niveau des capteurs. Quant la contrainte orthoradiale ( ), on ne
constate pas dvolution sensible.
Quant l'interface sol-pieu, son comportement est li l'volution de la contrainte normale suite
l'installation et au chargement et plus particulirement l'volution du coefficient de pression
des terres K (Boulon et al. (1986), Lehanne et al. (1993)). Kraft (1991) a dcrit de manire
schmatique l'volution de ce coefficient dans le cas de sables lches et denses (figure 2.12).
L'tat de contrainte initial proche de K0 peut voluer jusqu' la condition de pousse latrale
Chapitre 02 : Gnralit

passive Kp (trajet 0-1). Le passage de la pointe entrane une diminution de K en raison des
phnomnes de broyage des grains et du comportement contractant du sol sous la pointe lors de
son passage (en raison de la forte compression ainsi provoque) (trajet 1-2). A la fin du processus
d'installation, on constate une augmentation Ki et, par consquent, une augmentation de la
contrainte normale agissant sur le ft du pieu. A partir du point 3, la phase de chargement peut
entraner une augmentation de la valeur de K dans les sables denses (comportement de l'interface
sol-pieu dilatant) ou une diminution dans le cas des sables lches (comportement de l'interface
sol-pieu contractant).
Les tudes prcdentes ont mis en vidence l'importance de la prise en compte du mode
d'installation (notamment refoulant ou non refoulant) sur la dtermination de la capacit portante
des pieux, et plus particulirement la prise en compte des contraintes rsiduelles. De nombreuses
tudes ont ainsi t entreprises (Briaud and Tucker (1984), Poulos (1987), Foray et al. (1989),
Lehanne (1992), Chow (1995)).

Figure 2.8 : Dplacements (a) et dformations (b) du sol autour d'un pieu fonc dans le cas d'un massifsableux
(Robinski and Morrison (1964) et Vesic (1965))

Figure 2.9 : Zones de dformations du sol lors du fonage de modles de pieux dans du sable ; (a)
Dplacements verticaux observs (b) zones de sol compact et de sol refoul autour des
pieux et (c) zones de dplacements horizontaux du sol (Shakirev et al, 1996)
Chapitre 02 : Gnralit

Lehanne (1992) a prsent les mesures ralises sur un pieu instrument lImperial College
(diamtre 102 mm et longueur 6 m), fonc dans un massif de sable moyennement dense. Ces
mesures ont montr que, dans ce cas, le fait de ne pas prendre en compte les contraintes
rsiduelles entranait une sous estimation de 25% de la charge reprise par le ft.
Quantitativement, les mesures faites aprs installation, ont permis de mettre en vidence un effort
de compression de 14 kN au niveau de la pointe et une mise en traction quivalente 7,5 kPa sur
le ft (frottement latral ngatif).

En ce qui concerne le forage, peu de donnes exprimentales existent. Le comportement des


pieux fors gravitaires est souvent considr comme tant quivalent celui des pieux mouls en
laboratoire (mode d'installation non refoulant), c'est dire mis en place avant la mise en place du
sol. Nanmoins, de nombreuses techniques de forage (perforateur R-SOL, forage avec pression
d'injection, etc) produisent des phnomnes de refoulement analogues ceux produits lors du
fonage.
La mthode de fabrication d'un pieu ou d'un micropieu peut donc avoir une influence
significative sur sa rponse mcanique. En effet, dans le cas de modes d'installation entranant un
refoulement du sol (battage, fonage), les grandes dformations du sol occasionnes peuvent
entraner des modifications de l'tat de contraintes du sol au voisinage de la pointe et du ft.
Foray et al. (1989) ont effectu des essais en laboratoire afin de mettre en vidence l'influence du
mode d'installation sur la capacit portante des pieux dans les sables. L'tude a t faite dans une
chambre d'talonnage sur un massif de sable d'Hostun dense (ID = 0,80) et sous conditions K0.
Les inclusions, instrumentes, ont t installes par moulage, fonage et battage, puis charges.
Les rsultats en termes de frottement latral moyen (fs) et de rsistance en pointe (qp) sont
prsents sur la figure 2.13 (a) et (b) . L'influence de la procdure d'installation sur les rsultats
est dterminante. En effet, on constate que le frottement latral et la rsistance en pointe les plus
faibles sont obtenus pour l'inclusion moule. Ces figures mettent aussi en vidence l'influence de
la contrainte verticale. On note galement les diffrences significatives observes par rapport aux
recommandations API.

2 .9.2Indice de densit et contrainte de consolidation :

Puech et al. (1979) ont ralis des essais darrachement monotone sur un pieu modle moul dans
un sable et ceci pour deux indices de densit (0,20 et 0,70). La fig 2.11 (a) prsente lvolution du
coefficient K au cours de larrachement. On observe une augmentation de ce coefficient K (et
Chapitre 02 : Gnralit

donc de la contrainte normale), par rapport au coefficient K0 au cours de larrachement, de lordre


de 300 % et 800 % pour les densits faible et forte respectivement. Cette augmentation de la
contrainte horizontale en cours darrachement est explique par le fait que le sol a tendance
augmenter de volume, tendance qui va se trouver empche par le massif environnant (notion de
dilatance empche, Schlosser et Guilloux (1981)). Le phnomne de dilatance du sable est
dautant plus important que la densit est leve et que la contrainte est faible. L'volution du
frottement latral avec la contrainte de consolidation a dj t prsente sur la fig 2.10. Francis
(1997) a galement tudi, grce une tude ralise en chambre d'talonnage, l'influence de la
contrainte de consolidation ainsi que l'influence de l'indice de densit sur le frottement latral, et
a montr des rsultats similaires ceux prsents ci-dessus en ce qui concerne la rsistance en
pointe, de nombreux essais ont t effectus en chambre d'talonnage (Schmertmann (1978),
Baldi et al. (1981), Houlsby and Hitchman (1988), Canou (1989), etc).

Figure 2.10 (a) : Evolution du coefficient de pression latrale avec l'arrachement statique : sable lche et
sable dense (Puech et al. 1979)
Chapitre03 : Mthodes danalyse des pieux

CHAPITRE 03 :

Mthode danalyse Des Pieux


Chapitre03 : Mthodes danalyse des pieux

Mthodes danalyse des pieux

3- Pieu sous charge latrale statique

3-1Introduction :
La mthode base sur les dformations est celle qui est le plus couramment utilise pour le
dimensionnement dun pieu sous charge latrale statique. Dans le cadre de cette mthode, deux
hypothses distinctes peuvent tre dgages : la mthode au module de raction et la mthode du
continuum lastique.
3.2Mthode Au Module De Raction :
Cette mthode est la plus couramment utilise .elle est base sur le model de Winkler (1867).
Elle consiste modliser linteraction entre le sol et le pieu par une srie de ressorts indpendants
entre eux et de raideur variable. La raideur permet de relier directement la raction latrale do sol
(P) et de dplacement du pieu (Y). Cette mthode est la base des courbes P-Y, les ressorts
ont alors un comportement non linaire.
3.2.1Principe gnral :
Le model de Winkler nonc en 1867, dfinit le sol comme tant un empilement de tranches
indpendantes. Chaque tranche de sol est modlise par un ressort horizontal fig3-1 sur lequel
sappuie le pieu.

Figure 3.1 : Reprsentation de modle de Winkler.


Chapitre03 : Mthodes danalyse des pieux

La pression sur une tranche de sol ne dpend que du dplacement horizontal de cette dernire
et dun coefficient de raction du sol (lq 3-1)
P = k h ( z)Y ( z) (3.1)
Avec :
P : pression (N/m2).
Kh (z) : le coefficient de la raction une profondeur z (N/m3).
Y (z) : le dplacement du pieu pour une profondeur z (m).
O,
P=ES Y (z ) (3.2)

En posant ES= k h ( z ) B
Avec :
P : la raction d sol (N/m).
ES : le module de raction du sol pour une profondeur z (N/m2).
B : le diamtre de pieu (m).
Le pieu est idalis comme une poutre lastique charge latralement. Le sol modlis par des
ressorts horizontaux indpendants les uns des autres, et de raideur ES.
Le comportement de la poutre, en flexion dans le plan (Y-Z) rsume :
2Y
M = EP I p (3.3)
Z 2
Leffort tranchant est gal :
M T
T= et = P (3.4)
Z Z
A partir des quations (3.2), (3.3) et (3.4) on peut crire lquation dquilibre statique sur un
tronon de pieu

4Y 4Y
EP I P = P EP I P + ES Y (Z ) = 0 (3.5)
Z 4 Z 4
Avec :
EP : module dlasticit du pieu.
IP : moment dinertie de la section transversale du pieu
Z : profondeur dans le sol.
Le premier terme de cette quation est un terme damplification et le second dattnuation.
Chapitre03 : Mthodes danalyse des pieux

Les solutions de cette quation peuvent tre obtenues soit par une mthode analytique soit par
une mthode numrique.
Le principal avantage de cette mthode est quen tout point le long de pieu, linteraction sol pieu
peut tre dfinie. Mais cette dfinition est restreinte par lhypothse que la pression en un point
est fonction linaire du dplacement en ce point et par sa dpendance au choix du profil des
valeurs de kh caractrisent le sol.
3.2.2 Rsolution dans le cas dun sol homogne :
Dans le cas dun sol sec et homogne ou le module de raction du sol est constant Kh(z) = Kh
quelle que soit la profondeur considre, il est possible de rsoudre lquation (3.5)
La solution gnrale de cette quation diffrentielle est (avoir par exemple Frank 1999) :
z z
z z z z
Y ( z) = e l0
cos + sin + e l0 cos + sin (3.6)
l0 l0 l0 l0

Avec , , , les constantes dintgration dtermines partir des conditions aux limites (en
tte et en pied) ;
I0 : la longueur de transfert du pieu ou longueur lastique, elle peut tre dfinie en premire
approche comme la longueur de pieu minimale pour laquelle le chargement latral en tte existe.
Lexpression de la longueur de transfert I0 est la suivante :

4E P I P
I0 = 4 (3.7)
Es

Cette solution permet bobtenir a tout niveau de sol, les expressions du moment flchissant M ;
de leffort tranchant T et la pression du sol P
M = E P I P Y ''
T = E P I P Y ( 3)
E P I PY ( 4) P
P= =

Remarque :
le pieu peut tre considre ( Frank ,1999 ) comme :
D
Souple (ou long) si 3.
I0
D
Rigide (ou court) si 1.
I0
Chapitre03 : Mthodes danalyse des pieux

Avec D : longueur de fiche

3.2.3 Rsolution dans le cas dun sol non-homogne :


Dans le cas dun sol homogne, pour lequel la raction du sol nest plus constante et est une
fonction complexe sur toute la profondeur, lq. (3.6) ne peut plus tre rsolue de manire
analytique.
Le recours au calcul numrique est ncessaire. Plusieurs solutions sont explicites dans la
littrature. En effet, le problme principal est de dterminer le module de raction du sol ES
Cest la mthode des courbes P-Y ou courbes de raction.
3-3-Expression Du Modle De Raction :
La dfinition du profil de module de raction est la difficult principale de la mthode au module
de raction il dpend en effet de nombreux paramtres tels la rigidit du pieu, le niveau du
chargement, et la nature du sol etc.
Les rsultats dessais pressiomtriques sont dusage courant pour le dimensionnement des
fondations.
Hadjadji (1993) recense la plupart des formulations publies dans la littrature. il conclut que le
module de raction ES du sol peut tre dtermin si on a obtenu le module dYoung E par des
essais en laboratoire ou le module pressiomtriques EM par des essais en place .
On retiendra :
a) Terzaghi :
Terzaghi (1955) dtermine le module de raction du sol ES partir du module dYoung du
matriau E constituant le sol
Terzaghi propose comme relation pour du sable :
ES 1
= = 0.74 (3.8)
E 1.35
Avec :
E = A Z
A : un coefficient adimensionnel fonction de la densit du massif de sable (tableau 1)
: Le poids volumique du sable.
Z : la profondeur.
Tableau 1 : valeur du coefficient A en fonction de la densit du sable.
Chapitre03 : Mthodes danalyse des pieux

Densit Du Lche Moyen Dense


Sable
A 100-300 300-1000 1000-2000
b) Poulos :
Poulos (1971) propose pour des sables (sols sans cohsion) une valeur moyenne du module de
raction en fonction du type de sol (tableau 2). Ce module est dtermin partir dessais sur des
pieux rels raliss par Broms (1964) pour des sols non cohrents
Tableau 2 : valeur moyenne de ES pour des sols sans cohsion.
Densit Du Intervalle de la valeur ES
Sable (kg/m2)
Lche 91400-210920
Moyen 210920-421840
Dense 421840-984300
Poulos a tabli que le rapport entre le module de raction du sol ES et le module dYoung du
matriau E constituant le sol est gal :
ES
= 0.82 (3.9)
E
c) Mnard, Bourdon et Ganbin :
Ces auteurs (Mnard & al ,1969) proposent de calculer le rapport entre le module de la raction
du sol ES et le module pressiomtrique EM en fonction dun coefficient rhologique , du
diamtre du pieu B et dun diamtre de rfrence B0 qui gale 0.6 m (q.3.10)
Le coefficient rhologique nous est donn en fonction de la nature du sol (tableau 3).

18
pourB B0
4(2.65 ) + 3

ES
= 3 (3.10)
EM
pourB B0
2 B0 B
2,65 +
3 B B0 2
Chapitre03 : Mthodes danalyse des pieux

Tableau 3 : facteur rhologique pour divers types de sols daprs Baguelin &al (1978)
Type
Tourbe 1
Argile 2/3
Limon
Sable 1/3
Sable et gravier 1/3

d) Matlock et Reese :
La mthode de Matlock et Reese (1960) permet de dterminer le module de raction du sol ES
partir dune analyse non linaire des courbes P-Y exprimentales qui fait intervenir la notion de
module scant en chaque point de la courbe.
Mc Clelland et fochat (1958) sont lorigine de cette loi. Ils essayrent initialement de corrler
des courbes P-Y avec des essais triaxiaux. Lexpression du module de raction est gnralement
donne en fonction de la profondeur par une loi puissante.
Remarque :
Toutes ces mthodes sont cependant limites. En effet, linteraction sol- pieu est rduite
lhypothse que la pression ou la raction du sol pour un tronon est une fonction linaire du
dplacement. La reprsentation dun sol par une loi de comportement lastique semble illusoire.
Il est ncessaire de modifier les hypothses de la mthode du module de raction et de supposer
que la raction du sol en tout point du pieu est une fonction non linaire du dplacement.
Lintroduction des courbes P-Y est alors indispensable pour la bonne modlisation du systme
sol / pieu.
3-4 Les Courbes de raction P-Y :
Ces mthodes ont t labores et dveloppes par des nombreux chercheurs par diverses
approches : essais in situ, essais de laboratoires, modlisation physique ou modlisation
numrique.
La justification des pieux soumis des sollicitations transversales se fait le plus souvent partir
des mthodes qui ncessitent de modliser linteraction sol - pieu c'est--dire la loi de raction du
sol en fonction du dplacement horizontale du pieu communment appele courbe P-Y .
Chapitre03 : Mthodes danalyse des pieux

On retrouve ainsi ce type dapproche dans diffrents rglements (M.E. L.T, Fascicule 62, 1993)
et les codes (A.P.I, 1993 ; P.H.R.I, 1980 ; D.N.V., 1992) lesquels proposent chacun leur
procdure pour dterminer les courbes P-Y partir dessais in situ ou en laboratoire.
a)-Fascicule 62 :
les articles du Fascicule 62 concernant les courbes P-Y ont t rdigs partir de linterprtation
faites des essais in situ au pressiomtre ( C.T.R.E.N4 Pressiomtres ,1991) . Mnard (1969)
fait lanalogie entre le tassement dune fondation superficielle uniformment charge et un pieu
charg latralement.
On dfinit un module KF gale deux fois le module de raction ES.
Il est calcul par la formule suivante :
12 E M

pourB B0
4 B0 B
2.65 +
Kf = 3 B B0 (3.11)
12 E M
pourB B0
4
(2,65) +

3
Avec :
: Un coefficient rhologique caractrisant le sol (tableau 4)

B0 = 0.6m.

Tableau 4 : facteur rhologique pour divers types de sols daprs le Fascicule 62.
Type Tourbe Argile Limon Sable Grave
EM EM EM EM
PI PI PI PI

Sur consolid ou trs - >16 1 >14 2/3 >12 1/2 >10 1/3
serr
Normalement 1 9-16 2/3 8-14 1/2 7-12 1/3 6-10 1/4
consolid ou
normalement serr
Sous consolid altr 1 7-9 1/2 5-8 1/2 5-7 1/3 - -
et remani ou lche
Chapitre03 : Mthodes danalyse des pieux

On admet que le sol exerce en chaque section de llment une raction perpendiculaire laxe
de celui-ci, fonction du dplacement transversal de la section considre.
Dans le cas de pieux de sections carrs ou circulaires, on considre que cette raction se compose
uniquement de pressions frontales.
La pression frontales est modlise par une pression uniforme sexerant suer la largeur de
llment perpendiculairement au sens du dplacement, note B.
La loi de mobilisation de la raction frontale en fonction du dplacement du pieu est dfinie par :
Un segment de droite passant par lorigine et de pente KF
Un palier PF gale B PF.
Cette loi est illustre par la (fig3.2) dans le cas de sollicitations de courte dure en tte
dominantes.

Figure3.2 courbe P-Y dans le cas de sollicitations de courte dure en tte dominantes
Pour certains calculs, par exemple vis--vis des sollicitations accidentelles trs brves (chocs), les
justifications peuvent tre menes partir dun diagramme dinteraction tel que celui reprsent
par la (fig 3.3)

Figure3.3 courbe P-Y dans le cas de sollicitations accidentelles trs brves en tte dominantes
Chapitre03 : Mthodes danalyse des pieux

Dans le cas de fondations allonges, de type barrettes, un frottement se dveloppe sur les surfaces
latrales des lments.
La courbe de raction globale P-Y est dcompose en deux courbes :
Une courbe de raction frontale.
Une courbe de raction tangentielle (fig 3.4) dfinie par
1. un segment de droite passant par lorigine et de pente KF
2. un palier PS.
Le module KS est pris gal au module de raction frontale KF (dfini prcdemment)
Le palier PS est pris gal :
PS = 2LS qS (3.12)
Ou LS est la longueur sur laquelle est calcul le frottement latral et qS le frottement latral
unitaire limite.

Figure3.4 loi de mobilisation de la raction tangentielle .


Dans les deux cas de fondations (pieux ou barrettes), le module de raction du sol et le palier doit
tre minors pour les zones proches de la surface
La profondeur ZC sur laquelle sapplique cette minoration, compte partir de la surface du sol
aprs travaux, est prise gale 2B pour les sols cohrents et 4B pour les sols frottants.
Pour Z Z C les lois effert-dplacement dfinie prcdemment est modifi par une affinit :
daxe y
de direction P

Z
de rapport : 0.5 1 +
ZC
Chapitre03 : Mthodes danalyse des pieux

b) - A.P.I. American Petroleum Institute, 1993:


Ce code amricain regroupe un ensemble de rgles techniques de conception et de calcul de
fondations. Il est reconnu pour le dimensionnement de plates-formes offshores.
Les essais en grandeur nature Mustang Island, Texas [Cox et al. 1974 ; Reese et al. 1974] sont
lorigine de la dtermination des P-y dans les sables tablis dans ce code.
Des modifications ont t apportes aprs les travaux de Murchison et al [1984].
A partir des essais de chargement latral dun pieu sur site, les relations semi- empiriques tablies
ont permis de donner lallure gnrale des P-y.
Celles-ci sont ensuite cales sur une banque de donnes dessais sur sites.
La raction latrale du sol pour le sable est non linaire. En labsence dinformation plus
restrictive, elle est approche, toute profondeur z par lexpression suivante:

kz
P = APu tanh y (3.13)
APu

Avec :
A : un facteur pour la prise en compte dun chargement continu cyclique ou statique ;
A = 0.9 pour un chargement cyclique
Z
A = 3 0.8 0.9 Pour un chargement statique.
B
PU : la raction ultime du sol la profondeur Z (KN/m).
PU la plus faible valeur de PUS = (C1 Z + C2 B)Z proche de la surface (3.14)

PUd = C3 BZ En profondeur
C1, C2 et C3 des coefficients dtermins partir dabaques qui sont fonction de langle
de frottement du sable (fig 3.5)
K : le module initiale de la raction du sol (KN/m3) comme une fonction de langle de
frottement interne et de la densit relative Dr ( fig 3.6)
Z : la profondeur (m).
Y : le dplacement latral (m).
PUS : raction ultime du sol en surface ;
PUS : raction ultime du sol pour les couches profondes ;
B : le diamtre du pieu ;
: Le poids volumique du sol ;
Chapitre03 : Mthodes danalyse des pieux

Pour les sols cohrents, linteraction sol pieu est base sur le modle tabli par Matlock (1970).
Celui-ci a ralis des essais sur un pieu de 0,32 in de diamtre fonc dans des argiles molles.
La seule modification apporte par lA.P.I. porte sur lexpression de la raction ultime du sol
PU. :
1
3
Y
P = 0.5PU (3.15)
yC

Ou PU : raction ultime du sol la profondeur Z (KN/m3)


YC = 2.5 C B (3.16)
Avec
C : Valeur moiti de la dformation correspondant au maximum de la contrainte dviatorique
dans un essai triaxial non drain.
Pu varie de 3CU 9CU quand Z varie de 0 ZC dtermin de la manire suivante :
CU Z
pU = 3CU + Z + J pourZ ZC
B (3.17)
PU = 9CU pourZ Z C

O
CU : cohsion non draine du sol (Kpa)
: Poids volumique du sol (Kn/m3)
J : constante empirique dont les valeurs comprises entre 0.25 et 0.5 ont t dtermines par des
essais sur site
Par exemple J =0.5 pour les argiles molles du Golfe du Mexique
ZC : profondeur critique
6B
ZC = (m)
B
+J
CU
On remarque que pour des argiles plus consistantes (CU> 96 kPa), lA.P.I. ne donne pas de rgle
particulire. Il recommande seulement de prendre en compte le caractre plus fragile de ces
argiles et par consquent la possible dtrioration rapide de leur capacit portante pour de grands
dplacements en chargement cyclique.
Chapitre03 : Mthodes danalyse des pieux

Figure3.5Coefficients C1, C2 et C3 en fonction de langle de frottement interne daprs A.P.L.

Figure3. 6 : Valeurs de K en fonction de la densite relative Dr et de langle de frottement interne daprs


A.P.L 1993

c)- PHRI [Port and Harbour Research Intitute , [1980 ]:


Suit des essais sur modle rduit 1g dun pieu soumis un chargement latral [Kubo,1965], le
rglement japonais propose une forme parabolique pour les courbes de raction:

P = K S ZY 0.5 Pour un sable (3.18)

P = K S Y 0.5 Pour une argile (3.19)

O
KS : module de raction latrale dans un sable
KC : module de raction latrale dans une argile
Chapitre03 : Mthodes danalyse des pieux

Cette relation a ensuite t confirme par des travaux raliss avec la centrifugeuse du P.H.R.I.
portant sur un pieu isol charg latralement dans du sable [Terashi et al., 1989]. Dans cette
tude, diverses configurations de pieux ont t testes. La rigidit EI du pieu a une influence sur
la rsistance du sol kS (si El diminue, kS augmente). Par contre linfluence est peu marque sur le
moment flchissant maximum ou sur les dplacements.
Ltude montre aussi que la rsistance du sol kS est inversement proportionnelle la racine carre
du diamtre B du pieu jusqu un diamtre proche de 80 cm. Terashi trouve ce rsultat
concordant avec une tude faite au P.H.R.I. par Takahashi et Kasugai (1987).
d) - Det Norske Veritas [1977]:
Cet organisme norvgien tablit, comme lA.P.I., des recommandations pour la construction de
structures offshore. Il se base sur linterprtation dessais in situ.
De nombreuses simplifications sont apportes tant donn la connaissance limite du
comportement latral dun pieu et lincertitude permanente concernant les donnes du sol.
La campagne de reconnaissance des paramtres gotechniques doit tre ralise dans les
directions horizontales et verticales sur au moins dix diamtres de pieu. Lanalyse de ces
caractristiques conduit la dcision dexcaver ou non les couches de sols suprieurs de faibles
qualits.
La courbe p-y reprsentant les caractristiques effort-dplacement dun lment de pieu de
diamtre B peut alors tre construite. Elle se dcompose en portions (fig 3.7).
La partie initiale, pour p pd et y .b , est une hyperbole dexpression:
Y 1 Y
= + (3.20)
P K 1 p d
O
Pd : rsistance latrale de dimensionnement
: Plin/Pd = rapport de rupture, toujours suprieur 1
Plim : valeur asymptotique de lhyperbole pour y >
k1 : pente initiale de la courbe
Chapitre03 : Mthodes danalyse des pieux

: coefficient dpendant du sol et des conditions de chargement (tableau 6)

Figure3.7 : Forme caractristique de la courbe P-Y daprs Norske Veritas (1977).


Lhyperbole est limite y = .b , et la rsistance latrale de conception Pd est dpendante de la
valeur donne au coefficient . Celui-ci a alors pour expression:
1
= (3.21)
p
1 d
K 1 b
Si : .b pd/K1 lhyperbole peut tre remplace par une droite de pente k1.
Pour des dplacements suprieurs (aprs le point o de la courbe p-y, fig3.7), des effets
spcifiques peuvent conduire une dtrioration de la rsistance latrale. Une rsistance
rsiduelle Pd est dfinie sur la base dessais de laboratoires.
Cet aspect est vident pour largile dans une gamme de grands dplacements y > .b.
Les valeurs de sont tablies pour largile dans le tableau 5.
Pour un sable, la rsistance latrale du sol nest pas dgrade. Det Norske Veritas fait donc
lhypothse que, pour des dplacements suprieurs y = .B, la rsistance latrale est constante et
gale la rsistance latrale rsiduelle : Pd = p d.
Cette rsistance latrale dpend du type de chargement. En statique, Pd = PdS est dfini par:
PdS = 4K P P'0 (3.22)

1 + sin d
O KP =
1 sin d
P0 : pression effective (aprs excavation)
d : Angle de frottement pour le projet
Chapitre03 : Mthodes danalyse des pieux

: Angle de frottement rel


Pour un chargement cyclique la rsistance latrale du projet Pdc est
Pdc = 3K P P'0 (3.23)
Pour des profondeurs suprieures deux diamtres B.
Pour des profondeurs Z infrieures deux diamtres lq (3.22) est remplace par:
Z
PdC = 3 K P P' 0 (3.24)
2B
On note que ces expressions sous-estiment sans doute la rsistance du pieu pour de grandes
profondeurs. Mais du fait quelles se placent du ct de la scurit et que la mobilisation des
contraintes de cisaillement y est faible, elles peuvent tre conserves ces profondeurs.
Linteraction sable-pieu est reprsente par la courbe p-y prcdemment dfinie (3.20) avec
Pd=Pd
Z
k1 = nh
B
O :
nh : coefficient de raction du sol (tableau 5)
Z : profondeur sous le niveau du sol (aprs excavation ventuelle).
Pour une argile, la rsistance latrale de projet dun pieu peut tre prise comme tant:
Cu
Pd = N P (3.25)
mc
O
CU : contrainte de cisaillement caractristique non draine reprsentative des conditions de
chargements considres.
mc : Coefficient de scurit caractristique du matriau
NP : constante semi-empirique augmentant linairement de la valeur 1 (en statique) ou O (en
cyclique) au niveau du sol jusqu 8 la profondeur z, gale Nr.B en dessous du niveau du sol
Nr =10 pour des argiles normalement consolides
=5 pour des argiles surconsolides
Linteraction argile-pieu est reprsente par la courbe p-y prcdemment dfinie (3.25) avec pour
pente initiale
Chapitre03 : Mthodes danalyse des pieux

Pd
K1 = 0.25
(3.26)
B( c )
O :
Pd : PdS pour un chargement statique
Pd : Pdc pour un chargement cyclique
: Coefficient empirique (tableau 5)
C :dformation verticale la moiti de la variation des contraintes principales dun essai
triaxial statique non drain sur un chantillon de sol non remani.
Tableau 5 : Valeurs recommandes pour nh et pour des chargements statiques et cycliques dans
le sable
Densit relative du sable
lche moyen Dense
nh (KPa/m) 5.0 12.0 18.0
0.04 0.04 0.04

Tableau6: Valeurs recommandes pour, , et dans largile


Paramtre Type de chargement Argile normalement consolide Argiles surconsolide
Statique 10 30
Statique 20 C 5C

Statique 80 C 8C
Cyclique 10 30
Cyclique 7.5 C 2.5 C

Cyclique 20 C 5C

e) Autres expressions des courbes P-y :


1. Li Yan & Byrne :[1992]
Utilisant la modlisation par la similitude du gradient hydraulique, Li Yan et Byrne [1992] ont
men une campagne dessais de chargement latral statique sur un pieu isol dans le sable. Le but
tait dtudier linteraction sol-pieu en termes de courbes de raction P-y. Comparant leurs
Chapitre03 : Mthodes danalyse des pieux

rsultats au rglement A.P.I., ils ont propos une nouvelle expression de lissage des courbes P-y
en deux segments pour une profondeur suprieure un diamtre de pieu.
Le premier segment est une droite lorigine de pente Emax applicable pour :
1
P y 1
= = (3.26)
E max B B
O :
Emax = 2Gmax (1 + ), le module de cisaillement maximal Gmax tant mesur dans le conteneur
exprimental, tant choisi gal 0,2
- 08
= 5(Dr) , Dr tant la densit relative en pourcentage
B : ayant une valeur proche de 0,5
Le second segment est une courbe parabolique dexpression:

P y
= (3.27)
E max B B
Le diamtre du pieu na pas dinfluence sur la pente initiale des courbes P-y. Par contre, pour de
grands dplacements, le palier plastique est dpendant du diamtre.
Linfluence de lexcentricit du chargement est tudie pour diffrentes profondeurs. Une
augmentation de lexcentricit conduit des moments flchissant plus importants et des
courbes P-y moins raides. Ceci est surtout vrifi prs de la surface. En profondeur, Li Yan et
Byrne [1992] concluent quun mme jeu de courbes de raction P-y peut tre utilis, en pratique,
pour dimensionner un pieu avec diffrentes configurations dexcentricit de chargement.
Un pieu fixe en tte est la configuration la plus courante en pratique. Compare celle dun pieu
libre en tte, ltude montre que le moment maximal dvelopp est moindre, donc que le pieu
peut fournir une rsistance latrale plus importante.
Par contre, pour une mme condition de sol et de pieu, un mme jeu de courbes de raction peut
tre utilis pour une liaison en tte fixe ou libre.
Ils confrontent leurs rsultats aux rgles donnes par lA.P.I. Les conclusions renforcent la
validit de la nouvelle expression propose pour les courbes de raction P-y.
2. Georgiadis et al. [1992]
Cette quipe a ralis deux tudes exprimentales sur des pieux chargs latralement.
La premire porte sur le comportement dans largile dun groupe de pieux. La modlisation est
faite sur modle rduit en gravit normale.
Chapitre03 : Mthodes danalyse des pieux

La seconde porte sur le comportement dun unique pieu charg transversalement dans du sable.
La modlisation est ralise en centrifugeuse 50g dans du sable compact manuellement une
densit de 16.3 kN.m-3
Dans les deux cas, les pieux sont quips de jauges sur leur fiche. Le dispositif exprimental est
suffisamment instrument pour obtenir les conditions aux limites et ainsi obtenir par double
drivation et double intgration les courbes de raction P-y.
Par ces deux approches, aux conditions exprimentales trs diffrentes, Georgiadis et al
obtiennent les mmes conclusions concernant la forme des courbes P-y:
Y
PU = (3.28)
1 Y
+
K PU
O
K : la raideur initiale de la courbe P-y ;
Pu : la rsistance ultime du sol ;
Y : le dplacement latral
Pour les sables, se basant sur la thorie de Terzaghi (1955), les auteurs tablissent que les valeurs
de K augmentent proportionnellement avec la profondeur : K = z.nh o n est un coefficient
dpendant de la densit du sable.
Pour la rsistance ultime PU, la valeur est le minimum de celles tablies par Reese et al.
K 0 Z tan sin
tan( ) cos
tan
+ (B + Z tan tan )
tan ( )
PU 1 = AZ + K 0 Z tan (tan sin tan )
(3.29)
K0 B

( ( )
pU 2 = AZ K a tan 8 1 + K 0 tan tan 4 )
Avec
: Le poids volumique (N/m3) ;
: Langle de frottement (degr) ;
K0 : coefficient des terres au repos ;
Chapitre03 : Mthodes danalyse des pieux

Ka : coefficient des terres actives ;

=2 ;

= 45 + 2

A : un facteur de profondeur qui dpend du rapport Z/D. les auteurs choisissent ce coefficient
entre 1 et 2, Reese et al (1974) recommandent entre 0.9 et 3
Le choix effectu pour mieux coller aux courbes P-Y exprimentales na pas dinfluence sur
la rponse du pieu.
3. Kouda et al. [1998]
A la suite du tremblement de terre de Kobe en 1995, des lacunes dans les mthodes de
dimensionnement ont t mises en avant. Cette quipe de chercheurs a travaill sur la dfinition
de relations p-y pour de grands dplacements. Des essais sur modles en centrifugeuse sont la
base exprimentale de ltude.
Le pieu modle est constitu de 13 cylindres indpendants et joints. Le dispositif exprimental
permet dappliquer une mme charge chaque tronon indpendamment. Les dplacements de
chaque lment du pieu sont aussi mesurs. Ainsi chaque palier de chargement, les conditions
aux limites dun tronon sont connues, permettant une analyse aise des courbes de chargement.
Le principe dinterprtation des donnes exprimentales nest malheureusement pas dcrit. La
courbe de raction obtenue a pour expression:
kh y
.
p n D
= (3.30)
nD kh
n y
1+ .
Pmax D
nD
O :
kh : module initial de raction la profondeur z,
Pmax: pression limite la profondeur z,
Y : poids volumique du sol,
D : diamtre du pieu,
n : acclration centrifuge.
Chapitre03 : Mthodes danalyse des pieux

Ces essais mis en oeuvre ne correspondent pas des essais reprsentatifs dun prototype mais ils
utilisent un procd original permettant dobtenir de grands dplacements du pieu sur toute sa
fiche.
Les courbes p-y proposes sont valides en calculant la rponse dun pieu. Ces rsultats sont
confronts ceux dun essai sur pieu modle isol centrifug. Une bonne concordance est
obtenue pour de faibles et de grands dplacements.
Suite au tremblement de terre destructeur de Kob en 1995, il a t mis en avant des problmes
lors du calcul des pieux soumis de grands dplacements.
3.5 Paramtres ncessaires la construction des courbes P-y rglementaires :
Nous avons rsum dans le tableau ci -dessous les paramtres de sol s ou de pieu ncessaires la
construction des courbes P-y pour les trois rglements
Nous nous sommes placs dans le cas dun sable sec.
Tableau 7 : Paramtres ncessaires la construction des courbes P-y
Rglement Paramtre de sol Paramtre de pieu
Fascicule 62 EM le module pressiomtrique B le diamtre du pieu
A.P.L/D.N.V langle de frottement interne B le diamtre du pieu
ID lindice de densit
P.H.R.I - B le diamtre du pieu (pour B<0.8m)
Il est important de remarquer quaucun des rglements ne prend en compte les efforts appliqus
sur la structure, ni mme la mthode de mise en place des pieux.
3.6) Mthodes du continuum lastique
Les mthodes du continuum lastique supposent que le massif de sol est un milieu continu et
lastique. La plupart du temps, le sol est galement considr comme homogne et isotrope. Ces
mthodes reposent sur la solution des quations de Mindlin (1936).
3.6.1) Mthode de Poulos :
Le comportement du pieu et le sol est caractrise par deux paramtres, le module dlasticit ES
et le coefficient de Poisson.
Poulos (1971) modlise le pieu par une plaque verticale de largeur B et de longueur L.
La rigidit la flexion est constante sur toute la longueur et gale EpI p.
Le pieu est dcompos en n lments de mme longueur, chaque lment tant soumis
une contrainte horizontale constant P.
Lexpression du dplacement et de la rotation du pieu est fonction de deux paramtres :
Chapitre03 : Mthodes danalyse des pieux

Llancement du pieu L/B et le facteur de flexibilit KR (rapport entre la rigidit la flexion et le


produit du module dlasticit et de la longueur du pieu).
Des abaques permettent de calculer le dplacement et l a rotation du pieu partir des paramtres
noncs ci-dessus (ES, EpI p, L/B, KR,) et de la charge latrale applique.
Cette mthode est toutefois limite. En effet le sol ntant pas lastique isotrope il est difficile
dobtenir les coefficients ES et parti r dessais gotechniques standards.
3.6.2) Mthode de Banerjee & Davis :
Lapproche de Banerjee et Davis (1978) utilise la solution analytique des quations de Mindlin
(1936). On suppose que le module dlasticit ES varie linairement avec la profondeur.
Deux variables adimensionnelles KR et X permettent, partir dabaques de dterminer les
dplacements du pieu ainsi que les moments.
Cette mthode, trs proche de celle de Poulos (1971) prsente les mmes inconvnients
(obtention de ES).
Toutefois, la simplicit dutilisation des abaques permet un dimensionnement ais dun pieu si on
admet un comportement lastique du sol.
3.7) Conclusion :
Il existe deux mthodes pour modliser le comportement dun pieu sous une charge latrale : la
mthode du continuum lastique et la mthode au module de raction.
1. La mthode du continuum lastique repose sur deux hypothses fondamentales : le sol est
un milieu lastique, continu et adhre au pieu au cours du chargement. Cette dernire
hypothse suppose une certaine rsistance la traction du sol qui dans la ralit est trs
faible. Cette mthode est aujourdhui trs peu utilise, mise part pour le calcul aux
lments finis.
2. La mthode au module de raction est aujourdhui la base de toutes l es rglementations
nationales et internationales. Elle dcrit linteraction sol- pieu par une loi de raction du
terrain en fonction du dplacement horizontal (courbe P-y). Toutefois, la modlisation du sol
par des courbes P-y suppose que chaque couche de sol est indpendante ce qui, dans la
ralit, est inexact. Dailleurs, pour la modlisation dun pieu sous charge dynamique, Novak
(1991) ajoute une loi de comportement entre chaque couche. Qui plus est, les courbes P-y
rglementaires ne prennent pas en compte le mode dinstallation du pieu (for, refoulant, )
ou la rugosit de linterface sol
CHAPITRE 04 :

Mthode simplifie de calcul la rponse


dynamique des pieux
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

4-Mthode Simplifie De Calcul La Rponse Dynamique Des Pieux.

4-1-Gnralit sur les sismes :


4-1-1 Naissance Du Sisme :
Le tremblement de terre correspond une vibration du sol provoque par une libration soudaine
de lnergie de dformation accumule dans la croute terrestre ou dans la couche sous-jacente
appele manteau.
Lendroit ou se produit le sisme es appel foyer ou hypocentre le point de la surface du sol
directement au-dessus est appel picentre.
4-1-2 Caractrisation des sismes :
On peut caractriser la "taille" du tremblement de terre par diffrents paramtres :
! " #$%&'(# ) * +
* % , - ) * )
- ) "-
! + . /
. /
!) " 0(#%1!!$ )
* %
"- , - !
+ ,
"
Un paramtre important pour lvaluation des effets des sismes un endroit donn est :
!) , 231 2 . 3 1
4 5,6 5,7 / - 8 ,' -
5 5,9 * - : - !) )
; -
F<ma * * F=ma ,
Cest en termes de dacclration maximale ag au niveau du bedrock quon exprime les cartes de
zonation ou zonage sismique quantifiant le niveau sismique prendre en compte pour
lapplication des rgles de constructions parasismiques.
Un autre paramtre utile pour lvaluation des effets des sismes un endroit donn est :

67
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

,-
* - -
* - , - / -
! * /
(

Cette dure est au maximum de lordre de 60 s en zone trs sismique, mais n'est que de quelques
secondes en zone peu sismique.
La caractrisation la plus explicite d'un tremblement de terre est videmment constitue
d'acclrogrammes enregistrs dans la zone gographique considre, qui contiennent la fois
les aspects acclrations, dure et contenu frquentiel.
Les acclrogrammes constituent une donne d'action directement utilisable par la dynamique
des structures.
Enfin, les spectres de rponse constituent la caractrisation des tremblements de terre la plus
couramment utilise par l'ingnieur des constructions. Ils sont des sous-produits des
acclrogrammes et permettent un calcul simple des efforts internes dans une structure soumise
sisme.
4-1-3 Effets particuliers des sismes :

Des sables secs soumis vibration peuvent subir des tassements importants, qui peuvent tre
estims par des mesures des vides du sable. Ces tassements peuvent atteindre quelques dizaines
de cm.
Diffrentiels ou non, ils peuvent tre suffisants pour rendre une construction inutilisable.

Les sables fins saturs en eau peuvent tre le sige dun phnomne encore plus spcifique au
contexte sismique : la liqufaction des sols.
Lexplication du phnomne est la suivante.
La rsistance au cisaillement d'un sol est donne par la relation :
R = ( u) tg (4.1)

68
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Avec :
-
- u)
- u
Comme la mise en mouvement du sol entraine une monte de la pression interstitielle u, il y a
rduction de la rsistance R du sol, jusqu son annulation :
R=( u ) tg =0 (4.2)
Un matriau sans rsistance au cisaillement se comporte comme un liquide, do le terme
Liqufaction du sol.
Ce phnomne peut entrainer des effets catastrophiques :
- une couche non liqufiable situe au dessus dune couche liqufie peut glisser sur celle-ci, si
des pentes existent. Des dplacements de sol de grande amplitude sont possibles: ainsi, des
glissements de terrain de plusieurs km ont t observs au Prou
4-1-4- Propagation Des Ondes :
A partir du foyer, la secousse sismique se propage sous forme dondes de divers types, qui,
compte tenu de lhtrognit du sol, engendrent en surface un mouvement complexe
difficilement prdictible en un site donn.
On peut distinguer les types dondes suivants :

Elles prennent naissance dans le foyer et se propagent lintrieur de la terre sous deux formes :
> 2
Ces ondes se propagent avec une vitesse de 7 8 Km/s et saccompagnent dun changement de
volume (compression et dilatation alternes).
> 4
Ces ondes se dplacent avec une vitesse de 4 5 Km/s saccompagnent dune distorsion dans le
plan perpendiculaire la direction de propagation, provoque un cisaillement sans changement de
volume.
Al diffrence des ondes longitudinales, les ondes transversales peuvent se propager dans les
milieux liquides ou gazeux, en raison de linaptitude de ces milieux transmettre les efforts de
cisaillement.

69
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Ondes de surface : Les ondes de volume qui arrivent la surface de la terre produisent des
ondes de surface nintressant le sol que sur une profondeur extrmement faible.
On distingue principalement :
! # # + "
Ce sont des ondes pour lesquelles les points do sol dcrivant des ellipses dans le plan vertical de
propagation.
Ce mouvement est semblable au mouvement de la houle et entraine des compressions (ou des
tractions) ainsi que des cisaillements dans le sol.
! ? !
Ce sont des ondes pour lesquelles les points du sol se dplacement dans un plan tangent la
surface, perpendiculairement a la direction de propagation ; elles nengendrent que des
contraintes de cisaillement.
Les ondes de surface ont une vitesse de propagation de 1.5 5 Km/S dans les terrains compacts
ou rocheux et de 0.5 1.5 Km/S dans les terrains meubles.
Si on assimile le sol un milieu lastique et homogne caractris par le module dlasticit
longitudinal E, le coefficient de poison et la masse volumique , la vitesse VP des ondes P est
donne par :

Et la vitesse VS des ondes S par :

Avec :

Constantes de Lam (4.5)

Le rapport des vitesses de propagation des ondes P et S peut se mettre sous la forme :

70
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

4-2-Mthode De Calcul Propose :


4.2-1 Introduction :
La tendance qui se dveloppe a lheure actuelle pour dimensionner les pieux en zone sismique
consiste a dterminer la dforme dun profil de sol sous action sismique et a calculer les pieux en
admettent que leur dforme propre est identique a celle du profil de sol considre.
La mthode simplifie propose consiste a dfinir lquation de la dforme du profil de sol tudie
(profil homogne ou profil multicouche) correspondant au mode fondamental doscillations libres
de la colonne de sol.
4-2-2-Principe de la mthode :
Profil de sol homogne sur substratum rocheux :

Les paramtres gotechniques ncessaires pour la dtermination des proprits dynamiques du
terrain sont :
!
! 3
! 4

On considre un dpt sdimentaire dlimit par une base rocheuse en partie infrieure, on
considre un milieu lastique parfait, homogne, isotrope et continu de hauteur H
Dans la situation dune oscillation horizontale, un point situ la cot x, de la surface libre tel
que :
!"
#
$ #
Lquation de mouvement est :
% !" % !"
&
%" %
On introduit un vecteur propre reprsentant la fonction harmonique stationnaire qui admet pour
lquation (4.7) ;
Une solution particulire de type Fourier-Bernoulli, on considre une variable en fonction du
' " tel que : (! " ) ' " *
En remplaant (4.8) dans (4.7) on a :

71
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

) ' ++ " ) ++ ' " ,

$ #
-.. / 0 .. 1
- / 0 1
Ou

Comme les deux variables ' " et X(t) sont indpendantes (4.10) on obtient deux quations
indpendantes variables spares :
' ++ " $ ' " #
$
) ++ 2 3 ) #



La solution des quations (4.11) et (4.12) est :
' " 4567$" 89:5$" ;567 $" <
$ $
) =567 >9:5

Les constantes dintgrations M, N, C et < sont dtermines par les conditions aux limites ainsi
que les constantes A et B
les conditions aux limites :
#
?0 1
?1
4 * <

1 <&, <5
Lq (4.13) @ A=0
$ $ C
9:5 B # @ B 6 DEF9 > G #
A A

En consquence, lexpression de W et du vecteur propre :


)
H
>
6 C
$I ,
B

6 C
HI 9:5 J K #
B
i=1, 2, 3..

72
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

La priode de vibration : LI T
M PQ
NO I RS

La priode fondamentale :
B C
L LU $
B
De la mme manire :
B B B B
L ! LV ! LP ! LW
T & ,
Ou :
L L L L
! T! &! ,!
L LV LP LW
La priode est proportionnelle la hauteur H et inversement proportionnelle la vitesse de
propagation de londe de cisaillement.
La dforme pour le mode fondamental est :
C
X 9:5 &
B
Dans le cas dun dpt constitu de plusieurs couches de proprits peu diffrentes, la vitesse de
propagation de londe S peut tre value par :
]

[ !\ B\ *
YZ?
B
^_

Avec :
!\ ` vitesse de propagation couche k ;
B\ : Hauteur de la couche k
4.3-Calcul sous dformation sismique impose :
Il y a quelques dcades, le principe de la protection parasismique des fondations profondes
consistait sur dimensionner ces ouvrages (pieux larges et raides) de manire ce que les
ruptures par cisaillement qui pourraient se produire au cours dun sisme par suite des
dformations du sol ne compromettent pas la transmission des charges verticales dans le fut des
pieux

73
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Figure4.1 : dformation dun pieu sous charges latrales

Depuis quelques annes, cette approche a t abandonne, et la tendance qui se dveloppe


lheure actuelle, comme dans les recommandations de lApplied Technology Council [2] de
1978, consiste dimensionner les pieux de fondation de manire ce quils soient suffisamment
flexibles pour sadapter aux dformations du sol induites par le sisme (pieux longs et flexibles)
Une mthode simplifie consiste dterminer la dformation dun profil de sol sous action
sismique, et calculer les pieux en admettant que leur dforme propre est identique celle du
profil de sol considr ; dans une telle approche , on nglige la raction du sol , ce qui va dans le
sens de la scurit ; on nglige aussi les forces dinertie apportes par la structure au cours du
sisme , ce qui va cette fois en sens inverse du conservatisme ; toutefois le calcul aux lments
finis de fondations profondes faisant intervenir le module de raction du sol et linteraction sol-
structure a montr que le rsultat global diffrait peu de celui donn par la mthode simplifie ,
sauf peut-tre dans le cas de sols anormalement compressibles , type MEXICO[11]
Ces dernires annes , le nombre grandissant denregistrements de sismes dstructures a permis
de procder des tudes statiques sur les caractristique des mouvements du sol observs en y
incluant les effets de la gologie locale : ainsi partir de linvestigation de 187 enregistrements
californiens , TRIFUNAC et BRADY[10] ont obtenu des moyennes de pics dacclration plus
basses dans les sites alluvionnaires que dans les sites rocheux ( rapport de lordre de 0.8) ;
paralllement , les dpouillements de SEED[8] portant sur 104 acclrogrammes enregistrs en
Californie , ou dans dautres pays , et pour les quels les conditions de sol taient connues avec
une certaine prcision , ont montr que la moyenne des pics de vitesse et de dplacement taient
nettement plus grande dans les sites alluvionnaires (fig4-2), la prsence du sol provoquant un
dplacement des amplifications dans le domaine des grandes priodes ; en normalement 1g les
spectres de rponse correspondant ces acclrogrammes , SEED a montr que les immeubles de
grande hauteur fonds sur des alluvions profondes ou lches taient soumis des forces
sismiques bien plus grandes que ceux fonds sur sols rocheux lacclration au sol tant gale par
ailleurs .

74
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Figure4.2: diffrents formes de spectres de rponse (Seed et al 1974)


Ces rsultats ont t confirms par les dpouillements effectus par NEWMARK [6] sur 49
enregistrements californiens utiliss pour confectionner le spectre de calcul propos par
lUSNRC pour les centrales nuclaires ; en normalisant les pics dacclration dplacement
maximum du sol , correspondant au domaine des grandes priodes (T > 4 secondes ) , variait en
moyenne plus cart type de 0.30m sur les sites rocheux 0.90 sur les sites alluvionnaires .
La valeur importante des dplacements du sol dans le cas de sites ou les fondations profondes
sont exiges, qui peut atteindre exceptionnellement 1.60m en champ libre pour une acclration
de 1 g sur le site dalluvions lacustres particulirement lches de MEXICO, est un argument qui
fortifie lhypothse faite de considrer les pieux de fondations comme des lments flexibles
aptes supporter les dformations sismiques imposes .
La mthode simplifie propose consiste dfinir lquation de la dformation du profil de sol
tudi correspondant au mode fondamental doscillations libres de la colonne de sol ; au passage,
il sera utile de calculer le priode de ce mode de vibration, permettant dune part de dterminer le
dplacement maximum du sol pour une acclration en champ libre donne , dautre par de

75
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

comparer cette priode de vibration celle de la structure supporte , pour dceler sil ny a pas
mise en rsonance de lensemble sol-structure .

4.3.1 Dtermination des efforts dans un pieu colonne :


On considre un pieu colonne transvasant une couche de sol homogne et sappuyant sur un
substratum rocheux ; en reprenant lhypothse de base selon laquelle le pieu est suffisamment
flexible pour suivre la dforme du sol, on considre que la dforme du pieu est une demi-arche
de sinusode ; dans ce cas la ligne de charge et la courbe des moments sont aussi des sinusodes,
affines de la dforme a un coefficient de proportionnalit prs :

Considrons une poutre sur appuis simples de porte 2H et dont la ligne de charge est la
sinusode dquation :
a b acde
fg
h
(4.19)
Equation de leffort tranchant :
i b klc
jh fg
f h
(4.20)

an ob n
Equation du moment flchissant :
m b cde a b
o n o

anP ob nP
Equation de la dform :
pqr s cde a b
oP n oP

C
En liminant P(z) entre les quations (4.15) et (4.16), il vient :
4 t uvw x
B
Le moment maximum et la flche maximum tant obtenus la surface du sol (Z=H), et le
dplacement maximum du sol tant li lacclration en champ libre a par la relation
fondamentale :
y {
r y
z o

On en dduit la valeur maximum du moment en fonction de lacclration :


L
4Y|1 uvD T
B
76
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

On peut dailleurs crire cette relation en fonction de la vitesse des ondes de cisaillement :
uvD
4Y|1

En calculant la charge maximum P partir de la valeur de la flche maximum (Z=H), soit :


C V
} ~ uvw &
B

On en dduit la valeur maximale de leffort tranchant (Z=0) :


CL
8Y|1 uvD *
BV
Que lon peut aussi crire en fonction de V :
C
8Y|1 uvD ,
B
4.3.2 Dtermination des efforts dans un pieu flottant :
Lorsque le pieu considr ne descend pas jusquau substratum rocheux, on calculera comme
prcdemment la priode de vibration du mode fondamentale de la couche complte :
B
L #

Le dplacement relatif de la tte du pieu a pour valeur :


{ on
yJ cde K
o n

Et le moment flchissant correspondant sexprime par la relation :


uvD CB
4Y|1 J 567 K
B
Ou :
L CB
4Y|1 uvD J 567 K
B B

Dans les prescriptions du SEAOC[7], SEED et DONOVAN recommandation de considrer


comme substratum rocheux (bed-rock) les matriaux pour les vitesses des ondes de cisaillement
77
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

est au mois gale 750 m/s ; ils conseillent dautre part de ne pas prendre en compte dans les
calculs de profondeur H du substratum rocheux suprieure 150 mtres.

4.3.3-Application numrique et commentaires :


On considre une couche homogne dalluvions lches surmontant un substratum rocheux situ
36 m de profondeur ; la densit des alluvions est & et le module de cisaillement vaut
G=24000 kPa.

###
#c
&

n
2 * "

{ ckler
#
Si le champ libre est soumis une acclration de 0.3g (zone de forte sismicit), le dplacement
maximum la surface du sol a pour valeur :
y { ##
r y # * k
z o #
Dplacement relatif dun pieu flottant de 24m de profondeur :
o o
#*J cde K #* cde #*#T T k

Ainsi un pieu de 36 m reposant en pointe sur le substratum rocheux devra tre calcul pour un
dplacement en tte de 10.8cm, alors quun pieu flottant de 24 m devra tre calcul pour un
dplacement relatif de 5.4cm
4-3-4 Remarques importantes :
Pour les pieux colonnes , il conviendra de limiter lencastrement du pieu dans le substratum
rocheux de faon ce que la rotation en pied puisse seffectuer effectivement ; si le pieu est
parfaitement encastr ( rotation nulle ) , la rupture par cisaillement est quasi invitable , moins
que le sol ne soit suffisamment lche pour absorber la dforme dinflexion du pieu , on

78
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

rappellera ce sujet que pour une profondeur identique au pieu articul en pied , le pieu
parfaitement encastr sa base supporte pour un dplacement en tte donne des moments de
flexion doubles et des efforts tranchant quadruples.

Figure4.3: la dformation dun pieu

Il convient de signaler que le problme du calcul dun pieu sous sisme, ou cest le sol qui
provoque la dformation du pieu : est fort diffrent du problme du calcul dun pieu sollicite en
tte par un effort horizontal, ou cest le pieu qui dforme le sol :

Figure4.4: diagramme du moment et effort tranchant

79
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

On rappellera pour mmoire que la solution de ce dernier problme peut tre apporte par la
mthode de WINCLER ; si lon suppose que le module de raction du sol K est constant avec la
profondeur, on sera conduit utiliser les solutions de lquation diffrentielle (chapitre2) :

uv > #
?
?
(4.34)

Le module de raction du sol pouvant tre dduit de lessai pressiomtre classique.

4.4 Profil De Sol Constitu De Plusieurs Couches Homognes :


La mthode consiste dterminer la priode quivalente du premier mode de vibration du profile
multicouche pour apprcier le dplacement relatif au niveau du sol, et valuer la forme de la
dforme du profil correspondant ce mode fondamental.

4.4.1 Calcul de la priode fondamentale quivalente :


Linventaire complet des mthodes approches permettant dvaluer la priode quivalente dun
profil multicouche a t dress par DOBRY, OWEIS et URZUA dans article essentiel [3] paru en
1976.
4.4.1.1-mthode de SHIMA [9] ,1962:
Cette mthode, qui t reprise par ZEEVAERT [11] en 1972, puis par ONAMOTO en 1973
consiste faire la somme des priodes lmentaires de chacune des couches rencontres :
]
BI BI
L [ DEF9 LI *
I I

4.4.1.2mthode de MADERA [4] ,1970 :


80
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Cette mthode, qui t reprise par SCHNABLE en 1973, et recommand par le SEAOC [7] en
1974, consiste pondrer les vitesses des ondes du cisaillement dans les diffrentes couches :
BV
L ,
]I I BI
4.4.1.3mthode dAMBRASEYS [1], 1959 :
Cette mthode, qui a t reprise par IDRISS, et qui sera recommande dans les futures rgles
parasismiques PS.84, consiste a pondrer les modules de cisaillement des diffrentes couches :
B
L #
I BI
I BI

4.4.1.4 Mthode de MADERA [4], 1970 :


Cette mthode qui t reprise par CHEN en 1971, puis par URZUA en 1973, et recommande
par le SEAOC [7] en 1974, consiste appliquer successivement chaque couche la mthode
exacte applicable un systme bicouche : pour chacune des deux couches 1 et 2,
On peut trouver une quation similaire lquation 4.7 avec V constante, la solution de ce
systme de deux quations diffrentielles exige quatre conditions aux limites :
% *
@ * "
% 9AB
% 9AB
Lquation donnant la priode fondamentale de vibration de la bicouche est alors :
B L CL CL
" "
B L L L
On trouvera dans labaque suivant (figure4.10) une solution graphique de la priode
fondamentale correspondant la plus grande racine de lquation (4.28)
Pour un profil multicouche, on appliquera successivement la procdure prcdente au bicouche
constitue de la couche i et des couches 1 i-1 dont on aura pralable dtermin la priode
quivalente ; il sera toutefois admissible de prendre :
B B

B B
A condition que les densits des couches contiges prsentent une faible variation.
81
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Si les densits des couches contiges sont trs diffrentes, on pourra pondrer ces densits par la
hauteur des couches correspondantes.
Ainsi, la priode quivalente des 2 premires couches tant rigoureusement donne par la
solution de lquation (4.31), se transformant comme suit aprs pondration des densits par les
hauteurs correspondantes :
V BV L C L C LV
" " D
B B LV L L
On procdera par rcurrence jusqu la dernire couche, qui combine avec les (n-1) couches
prcdentes, donne la priode quivalente des n couches par rsolution de lquation suivante :
] B] L] C L] C L]
" "
] I BI L] L L
Dans les quations prcdentes, L] dsigne la priode quivalente des i premires couches
Pour utiliser la rsolution graphique donne par labaque de la figure 4.12, il suffira de
remplacer :
B B
D(
B B
Et :
B BI B BI
D(
I
BI I BI

4.3.1.5 Formule de RAYLIGH :


La premire itration de la procdure de RAYLEIGH, donne lexpression suivante de la pulsation
pour une forme linaire du mode fondamental :
BI

I
BV
4.4.2 Dtermination de la dforme dun profil multicouche :
On pourra utiliser la procdure de RAYLEIGH, dont lalgorithme est bas sur lquation exacte
donnant la pulsation :
w)
U ~t wx
Q


U ) wx
Q

Obtenue en galant lnergie cintique totale maximum et lnergie potentielle du systme


lorsquil vibre librement dans le mode fondamental.

82
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

La procdure considre conduit une mthode de rsolution itrative convergence rapide, qui
peut amplement se simplifier si, comme le propose DOBRY [3], on admet une densit constante
t et si lon fait lhypothse dune variation linaire de la dforme entre les points
extrmes dune mme couche :
B xI
)I )I BI
I

B xI ; dsignant le milieu de la couche considre.


Dans ces conditions, et en partant du pied du profil (X1=0), on pourra estimer X chaque
interface de couche, et partir de ces valeurs, utiliser une version modifie de lquation (4.33)
pour estimer la valeur de la pulsation :
w)
] ~ t wx
O
I
O
T
] O ) wx

O

Variables ( w)wx et X sont remplaces par des valeurs moyennes dans chaque couche :

) )I )I

w) )I )I B xI
wx BI I

Lquation (4.35) se transforme donc comme suit :


B xI
]I ~ BI
I
&
] )I )I BI
La priode de vibration se dduit de la pulsation en appliquant la relation fondamentale :
C
L

Lutilisation de la procdure de RAYLEIGH simplifie est fort intressante, car pour des profils
ou les couches sont caractrises par des vitesses de cisaillement peu diffrentes, le contour
polygonal de la dforme du premier mode pourra toujours tre assimil une sinusode voisine
de celle de la couche homogne, et partir de celle-ci il sera commode dvaluer les efforts dans
un pieu transvasant un profil multicouche.

83
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Par contre, lorsquon rencontrera de grandes variations se la vitesse de cisaillement, le calcul


mettra en vidence les points anguleux de la dforme, et cest aux interfaces correspondant ces
points anguleux quil conviendra de disposer des renforts darmatures.
4.4.3 Exemples de priode pour la rgion de SAN FRANCISCO :
Le SEAOC [7] donne les valeurs suivantes pour la priode fondamentale de sols multicouches
pour la rgion de SAN FRANCISCO
H(m) Sables compacts Remblai (5m) sur argile consolide
9 0.3 0.5 s 0.5 1.0 s
30 0.6 1.2 s 1.5 2.3 s
60 1.0 1.8 s 1.8 2.8 s
90 1.5 2.3 s 2.0 3.0 s
150 3.0 3.5 s -

4.5 Dtermination des efforts dans un pieu en tenu compte la charge axiale :
Le cas dun pieu traversant une couche de sol dpaisseur H pour se ficher dans le bed-rock la
base de la couche.

fig4.5: Dtermination des efforts dans un pieu en tenu compte la charge axiale.

En ngligeant les forces dinertie dans le pieu, lquation de sa dform U(Z) quant le sol a un
dplacement US(Z) scrit :
w Q
uv 4Q B x / x + x x+ x + wx + *
wx Q

Avec :
84
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

E : tant le module dYoung du matriau constitutif le pieu,


I : linertie de flexion,
MH et VH : le moment flchissant et leffort tranchant appliqus en tte du pieu,
Kt : la raideur transversale de linteraction sol-pieu tel que :
/ / Et

Diffrentes expressions des coefficients / et , faisant intervenir le coefficient de poison , ont


t proposes [AFP 98], conduisant des rsultats numriques assez variables (de 0.5 2 ou 3)
pour ces facteurs adimensionnels, lassociation franaise du gnie parasismique recommande,
compte tenu de cette dispersion de prendre simplement :
/ 5:6" / G
Pour les calculs de premire approximation.
En drivant deux fois par rapport Z, lquation (4.38) se met sous la forme :
wP
P P ,
w P

En posant : P #
\ Q Q
Q /

Les conditions aux limites associes (4.39) sont

#! # F" # D("69D"6:7 F7 :67"F


?
?-

w wV BV Q
! # F"
w w V uv
(Encastrement dans la masse m qui exerce leffort tranchant VH)
Pour rsoudre lquation (4.39), on suppose que US ( correspond la dforme du mode
fondamental de la couche de sol c'est--dire :
C
567

Le dplacement de la surface du sol ( tant tel que lacclration correspondante


soit gale une valeur donne avec, comme expression de la pulsation :
CE

B
VS : la vitesse des ondes de cisaillement dans la couche.
On a donc pour `
B
T
C EA
85
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Dans ces conditions, la solution U ( de lquation (4.39) se dtermine par des calculs
lmentaires ; on trouve lquation (4.46) qui suit :

P C
567
P C P

BV Q
9 5 9:55 9 5 5679
V uv 5

En ayant pos :
C=cos ; s = sin ; 9 5 &
Dou pour le dplacement en tte du pieu :
P BV Q
95 *
P C P V uv 5
Si on suppose que le btiment pos sur les pieux soit suffisamment rigide, il se comporte comme
une masse qui suit le mouvement impos par la couche de sol et la masse m de la figure (4.14)
reprsente le quotient de la masse totale du btiment par le nombre de pieux.
Lacclration Q en tte du pieu est donc :
Q (4.49)
Et on a pour leffort tranchant VH exerc sur la tte du pieu
Q Q T#

Do, compte tenu de lq (4.48), lexpression suivante de rapport Q des acclrations :


Q P C 95
T
P C P
/ BV 5

Ou lon a introduit la masse volumique de sol (=G/vs2).


Le moment flchissant MH en tte de pieu se calcul partir de la drive seconde de (voir
lquation 4.46) :

BP BEQ
4Q 95 T
/
P CP 5

86
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

Do en remplaant VH par m Q ( voir 4.50et en utilisant la formule 4.51) :

B BV 95 C 95
4Q V T
P C P /
5 / BV 5

Dans les mthodes simplifies de dimensionnement des pieux sous sollicitation sismique, la
plupart des codes parasismiques admettent que les efforts et les moments peuvent tre calculs en
supposant que les pieux suivent exactement la dformation du sol ; on trouve alors en tte le
moment MS donn par :

C uv nV /
4 B T
B P
Do en divisant lq (4.53) par lq (4.54) :
4Q P V 95 C 95
TT
4 P C P
/ BV 5 / BV 5

Les formules (3.51) et (4.55) permettent de mettre en vidence les principales caractristiques de
linteraction sol-pieu.
A prsence dun dnominateur pouvant a priori sannuler pour certaines valeurs des paramtres
fait apparatre un risque de rsonance (lamplitude de la rponse devient infinie en labsence
damortissement) ; la condition de rsonance scrit :
BV C 95
T
/ 5
Pour en apprcier la signification, il convient de prciser la gamme des valeurs numriques des
facteurs qui interviennent dans cette relation. Pour des pieux cylindriques de diamtre , en
prenant / on a :
C B P P
B
v HP P 2 3 2 3 T&
C u H Cu H
En prenant E= # U
8 (bton), *## V ; VS =100m/s ou 200m/s on trouve
pour `

87
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

B
# T :( ##5
H A

B
# :( ##5
H A

variant typiquement de 5 20, on voit que varie environ de 1 10 ; pour ces


Q
H
Le rapport

valeurs du le coefficient 95 5 qui intervient dans (3.56) :


X * B u P
2 3 T*
B CV }
On introduisant la pression statique (P=4mg/ CX ) dans les pieux on trouve avec les mme
valeurs qui prcdemment pour E, et VS, g=10 m/s2 et P=T # 8 :
X
# ##T, B :( ##5
B A

X
# ## #B :( ##5
B A

Ces conditions peuvent tre dans des cas plausibles (par exemple X #T, :(
B # . Le risque de rsonance fait donc partie des critres de vrification dun systme de
pieux en zone sismique, comme cest le cas dans les rgles PS92(AFN95).
La longueur H des pieux tant fixe par les caractristiques du site, cest en agissant sur le
diamtre X que lon peut viter le risque de rsonance.
La formule (4.55) montre que, mme si on sest plac assez loin de la rsonance, les efforts et les
moments subis par les pieux peuvent tre sensiblement majors par rapport ceux que lon
calcule en tenant compte que des dformations imposes par le sol. il ne faut donc pas oublier, en
dimensionnant les pieux, la part deffort qui rsulte des forces dinertie agissant sur la structure
supporte. Par contre, pour des ouvrages linaires entirement enterrs (canalisations, tunnels),
qui ne sont pas sollicits par des ractions inertielles provenant de parties lair libre , la pratique
courante consistant les calculer sous la seule action des dformations du sol est tout fait
justifie ; elle revient ngliger les effets dinteraction sol/structure dans des conditions ou cette
omission est ncessairement dans le sens de la scurit .

88
Chapitre04 : mthode simplifie de calcul
La rponse dynamique des pieux

89
CHAPITRE 05 :

Simulation Numrique

90
Chapitre 05 : Simulation Numrique

5. SIMULATION NUMERIQUE

5.1 Gnralits :
5.1.1 Introduction :
Le but de ce chapitre est de propose un modle de rfrence pour ltude, pour cela on va raliser
une modlisation de lexemple prsente dans le chap.4.Pour cet exemple, on se propose dtudier
la rponse dynamique dun pieu en calculant les efforts par la mthode simplifie dune part et
de voir linfluence des paramtres de calcul sur la variation des efforts avec le logiciel.
Le modle de rfrence est tabli et test avec le logiciel PLAXIS 8.2 [voir Annexe]
5.1.2 Choix de llment fini :
Pour les rglages gnraux, on choisi
1) Elments 15 nuds pour le sol.

2) Elments 2 nuds pour le pieu.


3) Problme de dformation plane.

5.1.3 Dfinition de la gomtrie et les proprits des matriaux :


La gomtrie du modle tudi est prsente sur la figure (5.1)
Hauteur du massif du sol H=36 m et la largeur L=2 H
Le massif du sol est compose de deux couches de sol
1ere couche : -5,00-0,00 m de profondeur argile rigide ;
2eme couche : 0- 36.00 m de profondeur sable, zone dencastrement.

91
Chapitre 05 : Simulation Numrique

5.1.4 Caractristiques des matriaux :


a- Proprits des couches de sol :

Les proprits des couches de sol utilises dans la simulation sont rsumes dans le tabl 5.1 :
Tableau 5-1 proprits des couches de sol
paramtre Dsignations Couche 1 Couche2
Modle et type de Mohr- coulomb drain Mohr- coulomb drain
comportement
Poids volumique apparent 16 kn/m3 16, 17,18 kn/m3
Poids volumique sature 18 kn/m3 20 kn/m3
Module dYong Eref 70.000 kn/m2 30.000 - 90.000 kn/m2
Coefficient de poison 0.35 0.25
Cohsion Cref 51 kn/m2 1 kn/m2
Angle de frottement 35 25
Angle de dilatance 0.2 0.2
b-Proprits des pieux :
Les proprits des pieux sont rsumes dans le tabl5.2 :
Tableau 5-2 proprits des pieux :
paramtre Dsignations
type de comportement Matriel type Elastoplastique
Rigidit normale EA 4.000X104 Kn/m
Rigidit de flexion EI 3.333X104 Kn/m
Epaisseur quivalente D 1.000m
Poids W 2.00kn/m
Coefficient de poison 0.3

Figure 5.1 : la gomtrie du modle

92
Chapitre 05 : Simulation Numrique

5.1.5 Maillage :
Licne maillage permet deffectuer automatiquement le maillage du problme. On rgle la
finesse du maillage dans le menu Mesh. Le maillage est prsent sur la fig 5.2

Figure 5-2 : Maillage du massif


5.1.6 Conditions initiales :
Une fois le modle gomtrique cr et le maillage d'lments finis gnr, l'tat de contraintes
initiales et la configuration initiale doivent tre spcifis. Les conditions initiales sont constitues
de deux modes diffrents, l'un pour gnrer les pressions interstitielles initiales (mode des
conditions hydrauliques) et l'autre pour spcifier la configuration gomtrique initiale et gnrer
le champ des contraintes effectives initiales (mode de configuration gomtrique).

Figure 5-3 : Conditions initiales

93
Chapitre 05 : Simulation Numrique

5.2 Mthode simplifie :


5.2.1 Calcul les efforts internes par la mthode propose :

4 Y|1
|
Le moment flchissant maximal est :

R
(5.1)
5.3 Calcul par la mthode numrique :
5.3.1 tape de calcul :
La mthode numrique se fait par les tapes suivantes :
1- On devise la profondeur du massif en dix couches c'est--dire le sol est dfinit comme
un empilement de tranches chaque tranche a une paisseur gale a H/10.

2- Ensuite on applique les incrments de dplacements pour chaque tranche .les valeurs des
dplacements sont en fonction de la vitesse, la densit et de la priode

wI D
| O
PM
(5.2)
LI
PQO
R
(5.3)



(5.4)

Le calcul numrique est fait laide de logiciel Plaxis .Le calcul du modle de rfrence se fait
selon les phases suivantes :
Phases 0 : Initiation Des Contraintes :
On value les contraintes effectives initiales
Phase 01 : Application des Dplacements,
On applique les dplacements sur le massif de sol cette phase est caractris par : Calcul
plastique Activer loption ignore undrained, Entre total multiplers weight =1
Phase 02 : Calcul Les Efforts Internes.
On dtermine les moments, et effort tranchant dans le pieu : Calcul plastique Activer loption
staged construction dans le menu loading input Activer le pieu et lancrage a partir du bouton
define qui nous permet daccder au menu de gomtrie.

94
Chapitre 05 : Simulation Numrique

5.3.2 Modle du calcul :


On considre une couche de sable avec des diffrentes densits ( dans une
zone de moyenne sismicit avec une acclration en champ libre de 0.15g , 0.2g puis dans une
zone de forte sismicit avec une acclration de 0.25g et 0.3g en champ libre, le module
dformation est pris gale a 50.000Kn/m2,
On calcul les dplacements quon applique sur le pieu, afin de dduire les valeurs de des
moments Mmax. Les rsultats obtenus sont rsums dans les tabls :
Tableau 5-3: les valeurs de dplacement pour V1=111.8 m/s et a =0.2g
Profondeur Vitesse priode Acclration Dplacement
(m)
0 111,8 0 0.2g 0
3,6 111,8 0.129 0.2g 0,00082532
7,4 111,8 0,265 0.2g 0,00348722
10,8 111,8 0,386 0.2g 0,00742786
14,4 111,8 0,515 0.2g 0,01320508
18 111,8 0,644 0.2g 0,02063294
21,6 111,8 0,773 0.2g 0,02971143
25,2 111,8 0,902 0.2g 0,04044056
28,8 111,8 1,03 0.2g 0,05282032
32,4 111,8 1,159 0.2g 0,06685072
36 111,8 1,288 0.2g 0,08253176
Tableau 5-4 : les valeurs de dplacement pour V1=111.8 m/s et a =0.3g

Dplacement
Profondeur Vitesse priode Acclration
(m)
0 111,8 0 0.3g 0
3,6 111,8 0,1288 0.3g 0,001238
7,4 111,8 0,26476 0.3g 0,0052308
10,8 111,8 0,3864 0.3g 0,0111418
14,4 111,8 0,51521 0.3g 0,0198076
18 111,8 0,64401 0.3g 0,0309494
21,6 111,8 0,77281 0.3g 0,0445671
25,2 111,8 0,90161 0.3g 0,0606608
28,8 111,8 1,03041 0.3g 0,0792305
32,4 111,8 1,15921 0.3g 0,1002761
36 111,8 1,28801 0.3g 0,1237976

95
Chapitre 05 : Simulation Numrique

Profondeur Vitesse priode Acclration Dplacement (m)


0 111,8 0 0.15g 0
3,6 111,8 0,1288 0.15g 0,00061899
7,4 111,8 0,26476 0.15g 0,00261542
10,8 111,8 0,3864 0.15g 0,00557089
14,4 111,8 0,51521 0.15g 0,00990381
18 111,8 0,64401 0.15g 0,0154747
21,6 111,8 0,77281 0.15g 0,02228357
25,2 111,8 0,90161 0.15g 0,03033042
28,8 111,8 1,03041 0.15g 0,03961524
32,4 111,8 1,15921 0.15g 0,05013804
36 111,8 1,28801 0.15g 0,06189882

Tableau 5-5 : les valeurs de dplacement pour V1=111.8 m/s et a =0.15g :

Tableau 5-6 : les valeurs de dplacement pour V1=111.8 m/s et a =0.25g :

96
Chapitre 05 : Simulation Numrique

Tableau 5-7 : les valeurs de dplacement pour V2=108.5 m/s et a =0.2g

Profondeur Vitesse priode Acclration Dplacement (m)


0 111,8 0 0.25g 0
3,6 111,8 0,1288 0.25g 0,00103165
7,4 111,8 0,26476 0.25g 0,00435903
10,8 111,8 0,3864 0.25g 0,00928482
14,4 111,8 0,51521 0.25g 0,01650635
18 111,8 0,64401 0.25g 0,02579117
21,6 111,8 0,77281 0.25g 0,03713929
25,2 111,8 0,90161 0.25g 0,0505507
28,8 111,8 1,03041 0.25g 0,06602541
32,4 111,8 1,15921 0.25g 0,0835634
36 111,8 1,28801 0.25g 0,1031647

Profondeur Vitesse priode Acclration Dplacement (m)


0 108,5 0 0.2g 0
3,6 108,5 0,133 0.2g 0,00087677
7,4 108,5 0,273 0.2g 0,00370462
10,8 108,5 0,398 0.2g 0,00789093
14,4 108,5 0,531 0.2g 0,01402831
18 108,5 0,664 0.2g 0,02191924
21,6 108,5 0,797 0.2g 0,0315637
25,2 108,5 0,929 0.2g 0,0429617
28,8 108,5 1,062 0.2g 0,05611325
32,4 108,5 1,195 0.2g 0,07101833
36 108,5 1,328 0.2g 0,08767695

Tableau 5-8 : les valeurs de dplacement pour V2=108.5 m/s et a =0.3g

97
Chapitre 05 : Simulation Numrique

Profondeur Vitesse priode Acclration Dplacement (m)


0 108,5 0 0.3g 0
3,6 108,5 0,13276 0.3g 0,0013152
7,4 108,5 0,27289 0.3g 0,0055569
10,8 108,5 0,39827 0.3g 0,0118364
14,4 108,5 0,53102 0.3g 0,0210425
18 108,5 0,66378 0.3g 0,0328789
21,6 108,5 0,79653 0.3g 0,0473456
25,2 108,5 0,92929 0.3g 0,0644426
28,8 108,5 1,06204 0.3g 0,0841699
32,4 108,5 1,1948 0.3g 0,1065275
36 108,5 1,32756 0.3g 0,1315154
Tableau 5-9 : les valeurs de dplacement pour V2=108.5 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 108,5 0 0.15g 0
3,6 108,5 0,13276 0.15g 0,00065758
7,4 108,5 0,27289 0.15g 0,00277847
10,8 108,5 0,39827 0.15g 0,00591819
14,4 108,5 0,53102 0.15g 0,01052123
18 108,5 0,66378 0.15g 0,01643943
21,6 108,5 0,79653 0.15g 0,02367278
25,2 108,5 0,92929 0.15g 0,03222128
28,8 108,5 1,06204 0.15g 0,04208493
32,4 108,5 1,1948 0.15g 0,05326375
36 108,5 1,32756 0.15g 0,06575771
Tableau 5-10 : les valeurs de dplacement pour V2=108.5 m/s et a =0.25g

Profondeur Vitesse priode Acclration Dplacement (m)


0 108,5 0 0.25g 0
3,6 108,5 0,13276 0.25g 0,00109596
7,4 108,5 0,27289 0.25g 0,00463078
10,8 108,5 0,39827 0.25g 0,00986366
14,4 108,5 0,53102 0.25g 0,01753539
18 108,5 0,66378 0.25g 0,02739905
21,6 108,5 0,79653 0.25g 0,03945463
25,2 108,5 0,92929 0.25g 0,05370213
28,8 108,5 1,06204 0.25g 0,07014156
32,4 108,5 1,1948 0.25g 0,08877291
36 108,5 1,32756 0.25g 0,10959618

Tableau 5-11 : les valeurs de dplacement pour V3=105.4 m/s et a =0.2g

Profondeur Vitesse priode Acclration Dplacement (m)


98
Chapitre 05 : Simulation Numrique

0 105,4 0 0.2g 0
3,6 105,4 0,137 0.2g 0,00092841
7,4 105,4 0,281 0.2g 0,00392283
10,8 105,4 0,41 0.2g 0,00835571
14,4 105,4 0,546 0.2g 0,0148546
18 105,4 0,683 0.2g 0,02321032
21,6 105,4 0,82 0.2g 0,03342286
25,2 105,4 0,956 0.2g 0,04549222
28,8 105,4 1,093 0.2g 0,05941841
32,4 105,4 1,229 0.2g 0,07520143
36 105,4 1,366 0.2g 0,09284127
Tableau 5-12 : les valeurs de dplacement pour V3=105.4 m/s et a =0.3g

Profondeur Vitesse priode Acclration Dplacement (m)


0 105,4 0 0.3g 0
3,6 105,4 0,13661 0.3g 0,0013926
7,4 105,4 0,28081 0.3g 0,0058842
10,8 105,4 0,40983 0.3g 0,0125336
14,4 105,4 0,54644 0.3g 0,0222819
18 105,4 0,68305 0.3g 0,0348155
21,6 105,4 0,81966 0.3g 0,0501343
25,2 105,4 0,95627 0.3g 0,0682383
28,8 105,4 1,09288 0.3g 0,0891276
32,4 105,4 1,22948 0.3g 0,1128021
36 105,4 1,36609 0.3g 0,1392619

Tableau 5-13 : les valeurs de dplacement pour V3=105.4 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 105,4 0 0.15g 0
3,6 105,4 0,137 0.15g 0,00069631
7,4 105,4 0,281 0.15g 0,00294212
10,8 105,4 0,41 0.15g 0,00626679
14,4 105,4 0,546 0.15g 0,01740774
18 105,4 0,683 0.15g 0,02506714
21,6 105,4 0,82 0.15g 0,03411917
25,2 105,4 0,956 0.15g 0,04456381
28,8 105,4 1,093 0.15g 0,05640107
32,4 105,4 1,229 0.15g 0,06963095
36 105,4 1,366 0.15g 0,01740774

Tableau 5-14: les valeurs de dplacement pour V3=105.4 m/s et a =0.25g

Profondeur Vitesse priode Acclration Dplacement (m)


0 105,4 0 0.25g 0

99
Chapitre 05 : Simulation Numrique

3,6 105,4 0,137 0.25g 0,00116052


7,4 105,4 0,281 0.25g 0,00490354
10,8 105,4 0,41 0.25g 0,01044464
14,4 105,4 0,546 0.25g 0,01856825
18 105,4 0,683 0.25g 0,0290129
21,6 105,4 0,82 0.25g 0,04177857
25,2 105,4 0,956 0.25g 0,05686528
28,8 105,4 1,093 0.25g 0,07427302
32,4 105,4 1,229 0.25g 0,09400179
36 105,4 1,366 0.25g 0,11605159

a=0.2g
a=0.3g

Figure 5.4 : Variation des dplacements en fonction de


l'acclration et de la vitesse

Remarque :

On remarque que pour la mme acclration, le dplacement augmente avec la diminution


de la vitesse.
Pour la mme vitesse le dplacement augment avec laugmentation de lacclration

Tableau 5.15:rsultats rcapitulatif :


La vitesse (m/s) Lacclration Dmax
111.8 0.15g 0,06189882
0.20g 0.08253176
0.25g 0,1031647
0.30g 0.1237976
108.5 0.15g 0,06575771
0.20g 0.08767695
0.25g 0,10959618
0.30g 0.1315154
100
Chapitre 05 : Simulation Numrique

105.4 0.15g 0,01740774


0.20g 0.09284127
0.25g 0,11605159
0.30g 0.1392619

5,96
5,9B
5,9
5,5C 5B
5,57 5D
5,56 5 9E
5,5B 5 BE
5
5 95 B5 D5 65

Figure 5.5 : Variation des dplacements en fonction de la zone

Daprs ces rsultats on rsume que le dplacement diminue avec diminution du la vitesse.
Le dplacement augment avec laugmentation de lacclration en champ libre
Le dplacement maximum tant obtenu la surface de sol

5.4 Analyse numrique par logiciel :


Pour tester la mthode simplifie on fait les calculs avec le logiciel pour la comparaison
Tableau 5.16 : variation du moment flchissant MS ; Mp
Mthode Mthode
Module de Priode Dplacement
densit Acclration EI Vitesse simplifi numrique MP/MS
dformation T D
Ms Mp
50000 1,6 1,472 3333 111,80 1,29 0,061895055 0,392361 0.3885 0,99047487
50000 1,7 1,472 3333 108,47 1,33 0,065763495 0,416883 0,4110 0,98588813
50000 1,8 1,472 3333 105,41 1,37 0,069631936 0,441406 0.4395 0,99568198
50000 1,6 1,962 3333 111,80 1,29 0,082526739 0,523148 0,6449762 1,23287521
50000 1,7 1,962 3333 108,47 1,33 0,087684661 0,555844 0,4474658 0,80502047
50000 1,8 1,962 3333 105,41 1,37 0,092842582 0,588541 0,54903 0,93286619
50000 1,6 2,453 3333 111,80 1,29 0,103158424 0,653935 0.3908 0,59761291
50000 1,7 2,453 3333 108,47 1,33 0,109605826 0,694806 0.3483 0,50129101
50000 1,8 2,453 3333 105,41 1,37 0,116053227 0,735676 0,515 0,70003643
50000 1,6 2,943 3333 111,80 1,29 0,123790109 0,784722 0,4398835 0,56055966

101
Chapitre 05 : Simulation Numrique

50000 1,7 2,943 3333 108,47 1,33 0,131526991 0,833767 0,3563576 0,4274067
50000 1,8 2,943 3333 105,41 1,37 0,139263873 0,882812 0,3795545 0,42993808

Figure
5.6 :
dfor
matio
n du
pieu

Le
tabl5.16 :
montre une
bonne
corrlation des
rsultats pour
les
acclrations
0.15g et 0.20g (moyenne),
Pour les acclrations plus grandes les rsultats divergents.

102
Chapitre 05 : Simulation Numrique

Figure 5.7 : Dplacements totaux

Figure 5.8 : moment flchissent

Daprs ces rsultats on remarque que les valeurs des moments par la mthode simplifie sont
augment avec laugmentation du la zone et la densit mais par la mthode numrique on trouve
un dcalage sur les valeurs. Donc on peut prend la densit = 1.8

5.4Etude paramtrique :

103
Chapitre 05 : Simulation Numrique

Aprs lapplication de la mthode simplifie et calculer les efforts et les dformations dun pieu,
ltude paramtrique sera effectue sur le mme modle montre leffet des diffrents paramtres
sur la rponse dun pieu avec les mmes donnes en changeant chaque fois le paramtre choisi
pour comparaison (les paramtres ont t modifi un par un).
Les paramtres tudis sont :
Module de dformation
La vitesse
La zone sismique
5.4.1 Module de dformation E=30.000kn/m2 :
Tableau 5.17 : les valeurs de dplacement pour V1=81.65 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 81,65 0 1,4715 0
3,6 81,65 0,17636252 1,4715 0,00116052
7,4 81,65 0,36252296 1,4715 0,00490357
10,8 81,65 0,52908757 1,4715 0,0104447
14,4 81,65 0,70545009 1,4715 0,01856836
18 81,65 0,88181261 1,4715 0,02901306
21,6 81,65 1,05817514 1,4715 0,04177881
25,2 81,65 1,23453766 1,4715 0,05686561
28,8 81,65 1,41090018 1,4715 0,07427344
32,4 81,65 1,58726271 1,4715 0,09400233
36 81,65 1,76362523 1,4715 0,11605226
Tableau 5.18 : les valeurs de dplacement pour V1=81.65 m/s et a =0.20g
Profondeur Vitesse priode Acclration Dplacement (m)
0 81,65 0 1,962 0
3,6 81,65 0,17636252 1,962 0,00154736
7,4 81,65 0,36252296 1,962 0,00653809
10,8 81,65 0,52908757 1,962 0,01392627
14,4 81,65 0,70545009 1,962 0,02475781
18 81,65 0,88181261 1,962 0,03868409
21,6 81,65 1,05817514 1,962 0,05570508
25,2 81,65 1,23453766 1,962 0,07582081
28,8 81,65 1,41090018 1,962 0,09903126
32,4 81,65 1,58726271 1,962 0,12533644
36 81,65 1,76362523 1,962 0,15473634

104
Chapitre 05 : Simulation Numrique

Tableau 5.19 : les valeurs de dplacement pour V1=81.65 m/s et a =0.25g


Profondeur Vitesse priode Acclration Dplacement (m)
0 81,65 0 2,4525 0
3,6 81,65 0,17636252 2,4525 0,0019342
7,4 81,65 0,36252296 2,4525 0,00817261
10,8 81,65 0,52908757 2,4525 0,01740784
14,4 81,65 0,70545009 2,4525 0,03094727
18 81,65 0,88181261 2,4525 0,04835511
21,6 81,65 1,05817514 2,4525 0,06963135
25,2 81,65 1,23453766 2,4525 0,09477601
28,8 81,65 1,41090018 2,4525 0,12378907
32,4 81,65 1,58726271 2,4525 0,15667054
36 81,65 1,76362523 2,4525 0,19342043
Tableau 5.20 : les valeurs de dplacement pour V1=81.65 m/s et a =0.30g
Profondeur Vitesse priode Acclration Dplacement (m)
0 81,65 0 2,943 0
3,6 81,65 0,17636252 2,943 0,00232105
7,4 81,65 0,36252296 2,943 0,00980713
10,8 81,65 0,52908757 2,943 0,02088941
14,4 81,65 0,70545009 2,943 0,03713672
18 81,65 0,88181261 2,943 0,05802613
21,6 81,65 1,05817514 2,943 0,08355762
25,2 81,65 1,23453766 2,943 0,11373121
28,8 81,65 1,41090018 2,943 0,14854689
32,4 81,65 1,58726271 2,943 0,18800465
36 81,65 1,76362523 2,943 0,23210451

Variation des dplacements en fonction du la zone


5,BE V1=81.65m/s E =30.000 Kn/m2

5,B

5,9E
@9A <5 9E
5,9
@B A <5 B5
5,5E G@D A <5 BE G

5 G@6 A <5 D

5 D,7 F,6 95,C96,6


95,C 9C B9,7
B9,7BE,BBC,CDB,6
DB,6 D7

Figure 5.9
9 : Variation des dplacements en fonction du la zone
V1=81.65m/s E =30.000 Kn/m2

105
Chapitre 05 : Simulation Numrique

5.4.2 Module de dformation E=50.000kn/m2 :


Tableau 5.21 : les valeurs de dplacement pour V2=105.41 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 105,41 0 1,4715 0
3,6 105,41 0,13660943 1,4715 0,00069631
7,4 105,41 0,28080827 1,4715 0,00294212
10,8 105,41 0,40982829 1,4715 0,00626679
14,4 105,41 0,54643772 1,4715 0,01114095
18 105,41 0,68304715 1,4715 0,01740774
21,6 105,41 0,81965658 1,4715 0,02506714
25,2 105,41 0,95626601 1,4715 0,03411917
28,8 105,41 1,09287544 1,4715 0,04456381
32,4 105,41 1,22948487 1,4715 0,05640107
36 105,41 1,3660943 1,4715 0,06963095
Tableau 5.22 : les valeurs de dplacement pour V2=105.41 m/s et a =0.20g
Profondeur Vitesse priode Acclration Dplacement (m)
0 105,41 0 1,962 0
3,6 105,41 0,13660943 1,962 0,00092841
7,4 105,41 0,28080827 1,962 0,00392283
10,8 105,41 0,40982829 1,962 0,00835571
14,4 105,41 0,54643772 1,962 0,0148546
18 105,41 0,68304715 1,962 0,02321032
21,6 105,41 0,81965658 1,962 0,03342286
25,2 105,41 0,95626601 1,962 0,04549222
28,8 105,41 1,09287544 1,962 0,05941841
32,4 105,41 1,22948487 1,962 0,07520143
36 105,41 1,3660943 1,962 0,09284127
Tableau 5.23 : les valeurs de dplacement pour V2=105.41 m/s et a =0.25g
Profondeur Vitesse priode Acclration Dplacement (m)
0 105,41 0 2,4525 0
3,6 105,41 0,13660943 2,4525 0,00116052
7,4 105,41 0,28080827 2,4525 0,00490354
10,8 105,41 0,40982829 2,4525 0,01044464
14,4 105,41 0,54643772 2,4525 0,01856825
18 105,41 0,68304715 2,4525 0,0290129
21,6 105,41 0,81965658 2,4525 0,04177857
25,2 105,41 0,95626601 2,4525 0,05686528
28,8 105,41 1,09287544 2,4525 0,07427302
32,4 105,41 1,22948487 2,4525 0,09400179
36 105,41 1,3660943 2,4525 0,11605159

106
Chapitre 05 : Simulation Numrique

Tableau 5.24 : les valeurs de dplacement pour V2=105.41 m/s et a =0.30g

Profondeur Vitesse priode Acclration Dplacement (m)


0 105,41 0 2,943 0
3,6 105,41 0,13660943 2,943 0,00139262
7,4 105,41 0,28080827 2,943 0,00588425
10,8 105,41 0,40982829 2,943 0,01253357
14,4 105,41 0,54643772 2,943 0,0222819
18 105,41 0,68304715 2,943 0,03481548
21,6 105,41 0,81965658 2,943 0,05013429
25,2 105,41 0,95626601 2,943 0,06823833
28,8 105,41 1,09287544 2,943 0,08912762
32,4 105,41 1,22948487 2,943 0,11280214
36 105,41 1,3660943 2,943 0,13926191

!"#$%&" ' ( !$#%### ) '


5,97

5,96

5,9B

5,9 @9A <5 9E


5,5C @B A <5 B5
5,57 G@D A <5 BE G

5,56 G@6 A <5 D

5,5B

5
5 D,7 F,6 95,C96,6 9C B9,7BE,BBC,CDB,6 D7

Figure 5.10 :
V2=105.41m/s E =50.000 Kn/m2

107
Chapitre 05 : Simulation Numrique

5.4.3 Module de dformation E=70.000kn/m2 :


Tableau 5.25 : les valeurs de dplacement pour V3=124.72 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 124,72 0 1,4715 0
3,6 124,72 0,11545863 1,4715 0,00049739
7,4 124,72 0,23733162 1,4715 0,00210161
10,8 124,72 0,34637588 1,4715 0,00447648
14,4 124,72 0,46183451 1,4715 0,00795818
18 124,72 0,57729314 1,4715 0,01243466
21,6 124,72 0,69275176 1,4715 0,0179059
25,2 124,72 0,80821039 1,4715 0,02437193
28,8 124,72 0,92366902 1,4715 0,03183272
32,4 124,72 1,03912765 1,4715 0,04028828
36 124,72 1,15458627 1,4715 0,04973862
Tableau 5.26 : les valeurs de dplacement pour V3=124.72 m/s et a =0.20g
Profondeur Vitesse priode Acclration Dplacement (m)
0 124,72 0 1,962 0
3,6 124,72 0,11545863 1,962 0,00066318
7,4 124,72 0,23733162 1,962 0,00280215
10,8 124,72 0,34637588 1,962 0,00596863
14,4 124,72 0,46183451 1,962 0,01061091
18 124,72 0,57729314 1,962 0,01657954
21,6 124,72 0,69275176 1,962 0,02387454
25,2 124,72 0,80821039 1,962 0,0324959
28,8 124,72 0,92366902 1,962 0,04244363
32,4 124,72 1,03912765 1,962 0,05371771
36 124,72 1,15458627 1,962 0,06631816
Tableau 5.27 : les valeurs de dplacement pour V3=124.72 m/s et a =0.25g
Profondeur Vitesse priode Acclration Dplacement (m)
0 124,72 0 2,4525 0
3,6 124,72 0,11545863 2,4525 0,00082898
7,4 124,72 0,23733162 2,4525 0,00350268
10,8 124,72 0,34637588 2,4525 0,00746079
14,4 124,72 0,46183451 2,4525 0,01326363
18 124,72 0,57729314 2,4525 0,02072443
21,6 124,72 0,69275176 2,4525 0,02984317
25,2 124,72 0,80821039 2,4525 0,04061988
28,8 124,72 0,92366902 2,4525 0,05305453
32,4 124,72 1,03912765 2,4525 0,06714714
36 124,72 1,15458627 2,4525 0,08289771

108
Chapitre 05 : Simulation Numrique

Tableau 5.28 : les valeurs de dplacement pour V3=124.72 m/s et a =0.30g

Profondeur Vitesse priode Acclration Dplacement (m)


0 124,72 0 2,943 0
3,6 124,72 0,11545863 2,943 0,00099477
7,4 124,72 0,23733162 2,943 0,00420322
10,8 124,72 0,34637588 2,943 0,00895295
14,4 124,72 0,46183451 2,943 0,01591636
18 124,72 0,57729314 2,943 0,02486931
21,6 124,72 0,69275176 2,943 0,03581181
25,2 124,72 0,80821039 2,943 0,04874385
28,8 124,72 0,92366902 2,943 0,06366544
32,4 124,72 1,03912765 2,943 0,08057657
36 124,72 1,15458627 2,943 0,09947725

!" &%* ' ( !*#%### ) '

5,9B

5,9

5,5C
@9A <5 9E
5,57 @B A <5 B5

5,56 G@D A <5 BE G


G@6 A <5 D
5,5B

5
5 D,7 F,6 95,C 96,6 9C B9,7 BE,B BC,C DB,6 D7

Figure 5.11
V3=124.72m/sE =70.000 Kn/m2

109
Chapitre 05 : Simulation Numrique

5.4.4 Module de dformation E=90.000kn/m2 :


Tableau 5.29 : les valeurs de dplacement pour V4=141.42 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 141,42 0 1,4715 0
3,6 141,42 0,10182435 1,4715 0,00038685
7,4 141,42 0,20930561 1,4715 0,00163457
10,8 141,42 0,30547306 1,4715 0,00348166
14,4 141,42 0,40729741 1,4715 0,00618962
18 141,42 0,50912176 1,4715 0,00967129
21,6 141,42 0,61094612 1,4715 0,01392665
25,2 141,42 0,71277047 1,4715 0,01895572
28,8 141,42 0,81459482 1,4715 0,0247585
32,4 141,42 0,91641918 1,4715 0,03133497
36 141,42 1,01824353 1,4715 0,03868515

Tableau 5.30 : les valeurs de dplacement pour V4=141.42 m/s et a =0.20g


Profondeur Vitesse priode Acclration Dplacement (m)
0 141,42 0 1,962 0
3,6 141,42 0,10182435 1,962 0,0005158
7,4 141,42 0,20930561 1,962 0,00217942
10,8 141,42 0,30547306 1,962 0,00464222
14,4 141,42 0,40729741 1,962 0,00825283
18 141,42 0,50912176 1,962 0,01289505
21,6 141,42 0,61094612 1,962 0,01856887
25,2 141,42 0,71277047 1,962 0,0252743
28,8 141,42 0,81459482 1,962 0,03301133
32,4 141,42 0,91641918 1,962 0,04177996
36 141,42 1,01824353 1,962 0,0515802
Tableau 5.31 : les valeurs de dplacement pour V4=141.42 m/s et a =0.25g
Profondeur Vitesse priode Acclration Dplacement (m)
0 141,42 0 2,4525 0
3,6 141,42 0,10182435 2,4525 0,00064475
7,4 141,42 0,20930561 2,4525 0,00272428
10,8 141,42 0,30547306 2,4525 0,00580277
14,4 141,42 0,40729741 2,4525 0,01031604
18 141,42 0,50912176 2,4525 0,01611881
21,6 141,42 0,61094612 2,4525 0,02321109
110
Chapitre 05 : Simulation Numrique

25,2 141,42 0,71277047 2,4525 0,03159287


28,8 141,42 0,81459482 2,4525 0,04126416
32,4 141,42 0,91641918 2,4525 0,05222495
36 141,42 1,01824353 2,4525 0,06447525

Tableau 5.32 : les valeurs de dplacement pour V4=141.42 m/s et a =0.30g

Profondeur Vitesse priode Acclration Dplacement m)


0 141,42 0 2,943 0
3,6 141,42 0,10182435 2,943 0,0007737
7,4 141,42 0,20930561 2,943 0,00326913
10,8 141,42 0,30547306 2,943 0,00696333
14,4 141,42 0,40729741 2,943 0,01237925
18 141,42 0,50912176 2,943 0,01934258
21,6 141,42 0,61094612 2,943 0,02785331
25,2 141,42 0,71277047 2,943 0,03791145
28,8 141,42 0,81459482 2,943 0,04951699
32,4 141,42 0,91641918 2,943 0,06266994
36 141,42 1,01824353 2,943 0,0773703

!"&"%& ' ( !+#%### ) '


5,5H
5,5C
5,5F
5,57
@9A <5 9E
5,5E
5,56 @B A <5 B5

5,5D G@D A <5 BE G

5,5B G@6 A <5 D


5,59
5
95,C
96,6

B9,7
BE,B
BC,C
DB,6
5
D,7
F,6

9C

D7

Figure 5.12 : Variation des dplacements en fonction du la zone


V4=141.42m/S E =90.000 Kn/m2

111
Chapitre 05 : Simulation Numrique

!#%"$
5,96

5,9B

5,9

5,5C 7<C9 7E
5,57 F<95E 69

5,56 C<9B6 FB
H<969 6B
5,5B

5
5 95 B5 D5 65

Figure 5.13 : Variation des dplacements en fonction du la vitesse


a=0.15g

!#% #

5,9C
5,97
5,96
5,9B
5,9 7<C9 7E
5,5C
F<95E 69
5,57
5,56 C<9B6 FB
5,5B H<969 6B
5
5 95 B5 D5 65

Figure 5.14 : Variation des dplacements en fonction du la vitesse


a=0.20g

112
Chapitre 05 : Simulation Numrique

!#% $

5,BE

5,B

5,9E
7<C9 7E
5,9 F<95E 69
5,5E C<9B6 FB
H<969 6B
5
5 95 B5 D5 65

Figure 5.15 : Variation des dplacements en fonction du la vitesse


a=0.25g

!#%,#
5,BE

5,B

5,9E
7<C9 7E

5,9 F<95E 69
C<9B6 FB
5,5E H<969 6B

5
5 95 B5 D5 65

Figure 5.16 : Variation des dplacements en fonction du la vitesse


a=0.30g

113
Chapitre 05 : Simulation Numrique

Tableau 5.33 : variation du moment flchissant MS ; Mp


Module
Densit G V A EI Mc Mp MP/MS
De dformation
30000 1,8 12000 81,65 1,4715 3333 0,73568 0,862 1,172
30000 1,8 12000 81,65 1,962 3333 0,98090 1,15 1,172
30000 1,8 12000 81,65 2,4525 3333 1,22613 0,936 0,763
30000 1,8 12000 81,65 2,943 3333 1,47135 0,907 0,616
50000 1,8 20000 105,41 1,4715 3333 0,44141 0,4395 0,996
50000 1,8 20000 105,41 1,962 3333 0,58854 0,552 0,938
50000 1,8 20000 105,41 2,4525 3333 0,73568 0,515 0,7
50000 1,8 20000 105,41 2,943 3333 0,88281 0,618 0,7
70000 1,8 28000 124,72 1,4715 3333 0,31529 0,398 1,262
70000 1,8 28000 124,72 1,962 3333 0,42039 0,575 1,368
70000 1,8 28000 124,72 2,4525 3333 0,52548 0,444 0,845
70000 1,8 28000 124,72 2,943 3333 0,63058 0,533 0,845
90000 1,8 36000 141,42 1,4715 3333 0,24523 0,263 1,072
90000 1,8 36000 141,42 1,962 3333 0,32697 0,3516 1,075
90000 1,8 36000 141,42 2,4525 3333 0,40871 0,315 0,771
90000 1,8 36000 141,42 2,943 3333 0,49045 0,378 0,771

-( .!#/"$
9,55555

5,H5555

5,C5555

5,F5555
0
5,75555

5,E5555

5,65555

5,D5555

5,B5555

5,95555

5,55555
5 B5555 65555 75555 C5555 955555

Figure 5.17 : Variation des moments en fonction du E


a=0.15g

114
Chapitre 05 : Simulation Numrique

-( .!#/ #
9,65555

9,B5555

9,55555
0
5,C5555

5,75555

5,65555

5,B5555

5,55555
5 B5555 65555 75555 C5555 955555

Figure 5.18 : Variation des moments en fonction du E


a=0.20g

-( .!#/ $
9,65555

9,B5555

9,55555
0
5,C5555

5,75555

5,65555

5,B5555

5,55555
5 B5555 65555 75555 C5555 955555

Figure 5.19 : Variation des moments en fonction du E


a=0.25g

115
Chapitre 05 : Simulation Numrique

-( .!#/,#
9,75555

9,65555

9,B5555

9,55555 0

5,C5555

5,75555

5,65555

5,B5555

5,55555
5 B5555 65555 75555 C5555 955555

Figure 5.20 : Variation des moments en fonction du E


a=0.30g

0-(
9,6

1<5,9E
9,B
1<5,B
9 1<5,BE

1<5,D
5,C

5,7

5,6

5,B

5
5 B5555 65555 75555 C5555 955555

Figure 5.21: Variation des moments en fonction du E

116
Chapitre 05 : Simulation Numrique

Et toutes les zones

- !2"%3$ '

9,75555

9,65555 0

9,B5555

9,55555

5,C5555

5,75555

5,65555

5,B5555

5,55555
5 5,E 9 9,E B B,E D D,E

Figure 5.22 : Variation des moments en fonction du la zone


E=30.000

moment-zone V=105.41m/s

9,55555
5,H5555 0
5,C5555
5,F5555
5,75555
5,E5555
5,65555
5,D5555
5,B5555
5,95555
5,55555
5 5,E 9 9,E B B,E D D,E

117
Chapitre 05 : Simulation Numrique

Figure 5.23 :
E=50.000
moment-zone V=124.72m/s

5,F5555

0
5,75555

5,E5555

5,65555

5,D5555

5,B5555

5,95555

5,55555
5 5,E 9 9,E B B,E D D,E

Figure 5.24 :
E=70.000

- !"&"%& '

5,75555
0
5,E5555

5,65555

5,D5555

5,B5555

5,95555

5,55555
5 5,E 9 9,E B B,E D D,E

118
Chapitre 05 : Simulation Numrique

Figure 5.25 :
E=90.000

Pour une zone de faible sismicit de 0.15g le rapport Mp/Ms= 0,491/0.4344=1.229


Pour une zone de moyenne sismicit de 0.25g le rapport Mp/Ms=0.657/0.491=1.135
Pour une zone de forte sismicit de 0.30g le rapport Mp/Ms= 0.5525/0.724=0.701
Le moment flchissant augmente quand lacclration augmente et diminue quand le module de
dformation de sol augment

Profondeur Vitesse priode Acclration Dplacement (m)


0 25,81 0 1,4715 0
3,6 25,81 0,55792329 1,4715 0,01161422
7,4 25,81 1,14684231 1,4715 0,04907365
10,8 25,81 1,67376986 1,4715 0,10452796
14,4 25,81 2,23169314 1,4715 0,18582748
18 25,81 2,78961643 1,4715 0,29035544
21,6 25,81 3,34753971 1,4715 0,41811183
25,2 25,81 3,905463 1,4715 0,56909666
28,8 25,81 4,46338628 1,4715 0,74330992
32,4 25,81 5,02130957 1,4715 0,94075162
36 25,81 5,57923286 1,4715 1,16142175
5.4.5 Variation du la vitesse V1=25.81m/s et E=3000 Kn/m2 :
Tableau 5.34 : les valeurs de dplacement pour V1=25.81 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 25,81 0 1,962 0
3,6 25,81 0,55792329 1,962 0,01548562
7,4 25,81 1,14684231 1,962 0,06543154
10,8 25,81 1,67376986 1,962 0,13937061
14,4 25,81 2,23169314 1,962 0,24776997
18 25,81 2,78961643 1,962 0,38714058
21,6 25,81 3,34753971 1,962 0,55748244
25,2 25,81 3,905463 1,962 0,75879554
28,8 25,81 4,46338628 1,962 0,99107989
32,4 25,81 5,02130957 1,962 1,25433549
36 25,81 5,57923286 1,962 1,54856233

Tableau 5.35 : les valeurs de dplacement pour V1=25.81 m/s et a =0.20g

119
Chapitre 05 : Simulation Numrique

Tableau 5.36 : les valeurs de dplacement pour V1=25.81 m/s et a =0.25g

Profondeur Vitesse priode Acclration Dplacement (m)


0 25,81 0 2,4525 0
3,6 25,81 0,55792329 2,4525 0,01935703
7,4 25,81 1,14684231 2,4525 0,08178942
10,8 25,81 1,67376986 2,4525 0,17421326
14,4 25,81 2,23169314 2,4525 0,30971247
18 25,81 2,78961643 2,4525 0,48392573
21,6 25,81 3,34753971 2,4525 0,69685305
25,2 25,81 3,905463 2,4525 0,94849443
28,8 25,81 4,46338628 2,4525 1,23884986
32,4 25,81 5,02130957 2,4525 1,56791936
36 25,81 5,57923286 2,4525 1,93570291

Tableau 5.37 : les valeurs de dplacement pour V1=25.81 m/s et a =0.30g

Profondeur Vitesse priode Acclration Dplacement (m)


0 25,81 0 2,943 0
3,6 25,81 0,55792329 2,943 0,02322843
7,4 25,81 1,14684231 2,943 0,09814731
10,8 25,81 1,67376986 2,943 0,20905591
14,4 25,81 2,23169314 2,943 0,37165496
18 25,81 2,78961643 2,943 0,58071087
21,6 25,81 3,34753971 2,943 0,83622366
25,2 25,81 3,905463 2,943 1,13819331
28,8 25,81 4,46338628 2,943 1,48661984
32,4 25,81 5,02130957 2,943 1,88150323
36 25,81 5,57923286 2,943 2,3228435

120
Chapitre 05 : Simulation Numrique

! $%2" '5 ( !,%### ) '


B,E

9,E
@9A <5 9E

9 @B A <5 B5
G@D A <5 BE G
5,E G@6 A <5 D

5
5 D,7 F,6 95,C96,6 9C B9,7BE,BBC,CDB,6 D7
4

Figure 5.26 :
V=25.81m/S E =3.000 Kn/m2

5.5.2 Variation du la vitesse V2=50m/s et E=11250 Kn/m2 :

Tableau 5.38 : les valeurs de dplacement pour V2=50 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 50 0 1,4715 0
3,6 50 0,288 1,4715 0,00309475
7,4 50 0,592 1,4715 0,01307629
10,8 50 0,864 1,4715 0,02785277
14,4 50 1,152 1,4715 0,04951604
18 50 1,44 1,4715 0,07736882
21,6 50 1,728 1,4715 0,1114111
25,2 50 2,016 1,4715 0,15164288
28,8 50 2,304 1,4715 0,19806417
121
Chapitre 05 : Simulation Numrique

32,4 50 2,592 1,4715 0,25067497


36 50 2,88 1,4715 0,30947527

Tableau 5.39 : les valeurs de dplacement pour V2=50 m/s et a =0.20g

Profondeur Vitesse priode Acclration Dplacement (m)


0 50 0 1,962 0
3,6 50 0,288 1,962 0,00412634
7,4 50 0,592 1,962 0,01743505
10,8 50 0,864 1,962 0,03713703
14,4 50 1,152 1,962 0,06602139
18 50 1,44 1,962 0,10315842
21,6 50 1,728 1,962 0,14854813
25,2 50 2,016 1,962 0,20219051
28,8 50 2,304 1,962 0,26408557
32,4 50 2,592 1,962 0,33423329
36 50 2,88 1,962 0,4126337

Tableau 5.40 : les valeurs de dplacement pour V2=50 m/s et a =0.25g


Profondeur Vitesse priode Acclration Dplacement (m)
0 50 0 2,4525 0
3,6 50 0,288 2,4525 0,00515792
7,4 50 0,592 2,4525 0,02179381
10,8 50 0,864 2,4525 0,04642129
14,4 50 1,152 2,4525 0,08252674
18 50 1,44 2,4525 0,12894803
21,6 50 1,728 2,4525 0,18568516
25,2 50 2,016 2,4525 0,25273814
28,8 50 2,304 2,4525 0,33010696
32,4 50 2,592 2,4525 0,41779162
36 50 2,88 2,4525 0,51579212

Tableau 5.41 : les valeurs de dplacement pour V2=50 m/s et a =0.30g

122
Chapitre 05 : Simulation Numrique

Profondeur Vitesse priode Acclration Dplacement (m)


0 50 0 2,943 0
3,6 50 0,288 2,943 0,00618951
7,4 50 0,592 2,943 0,02615257
10,8 50 0,864 2,943 0,05570555
14,4 50 1,152 2,943 0,09903209
18 50 1,44 2,943 0,15473764
21,6 50 1,728 2,943 0,2228222
25,2 50 2,016 2,943 0,30328577
28,8 50 2,304 2,943 0,39612835
32,4 50 2,592 2,943 0,50134994
36 50 2,88 2,943 0,61895055

5.4.6 Variation du la vitesse V3=100m/s et E=45.000 Kn/m2 :


Tableau 5.42 : les valeurs de dplacement pour V3=100 m/s et a =0.15g

Profondeur Vitesse priode Acclration Dplacement (m)


0 100 0 1,4715 0
3,6 100 0,144 1,4715 0,00077369
7,4 100 0,296 1,4715 0,00326907
10,8 100 0,432 1,4715 0,00696319
14,4 100 0,576 1,4715 0,01237901
18 100 0,72 1,4715 0,0193422
21,6 100 0,864 1,4715 0,02785277
25,2 100 1,008 1,4715 0,03791072
28,8 100 1,152 1,4715 0,04951604
32,4 100 1,296 1,4715 0,06266874
36 100 1,44 1,4715 0,07736882
Tableau 5.43 : les valeurs de dplacement pour V3=100 m/s et a =0.20g
Profondeur Vitesse priode Acclration Dplacement (m)
0 100 0 1,962 0
3,6 100 0,144 1,962 0,00103158
7,4 100 0,296 1,962 0,00435876
10,8 100 0,432 1,962 0,00928426
14,4 100 0,576 1,962 0,01650535
18 100 0,72 1,962 0,02578961
21,6 100 0,864 1,962 0,03713703
25,2 100 1,008 1,962 0,05054763
28,8 100 1,152 1,962 0,06602139
32,4 100 1,296 1,962 0,08355832
36 100 1,44 1,962 0,10315842
Tableau 5.44 : les valeurs de dplacement pour V3=100 m/s et a =0.25g
Profondeur Vitesse priode Acclration Dplacement (m)
123
Chapitre 05 : Simulation Numrique

0 100 0 2,4525 0
3,6 100 0,144 2,4525 0,00128948
7,4 100 0,296 2,4525 0,00544845
10,8 100 0,432 2,4525 0,01160532
14,4 100 0,576 2,4525 0,02063168
18 100 0,72 2,4525 0,03223701
21,6 100 0,864 2,4525 0,04642129
25,2 100 1,008 2,4525 0,06318453
28,8 100 1,152 2,4525 0,08252674
32,4 100 1,296 2,4525 0,1044479
36 100 1,44 2,4525 0,12894803

Tableau 5.45 : les valeurs de dplacement pour V3=100 m/s et a =0.30g


Profondeur Vitesse priode Acclration Dplacement (m)
0 100 0 2,943 0
3,6 100 0,144 2,943 0,00154738
7,4 100 0,296 2,943 0,00653814
10,8 100 0,432 2,943 0,01392639
14,4 100 0,576 2,943 0,02475802
18 100 0,72 2,943 0,03868441
21,6 100 0,864 2,943 0,05570555
25,2 100 1,008 2,943 0,07582144
28,8 100 1,152 2,943 0,09903209
32,4 100 1,296 2,943 0,12533749
36 100 1,44 2,943 0,15473764

124
Chapitre 05 : Simulation Numrique

Variation des dplacements en fonction du la zone


V=100m/S E =45.000 Kn/m2
5,9C
5,97
5,96
5,9B
5,9 @9A <5 9E
5,5C @B A <5 B5
5,57 G@D A <5 BE G
5,56 G@6 A <5 D
5,5B
5
5 D,7 F,6 95,C96,6 9C B9,7BE,BBC,CDB,6 D7
4

Figure 5.27 : Variation des dplacements en fonction du la zone


V=100m/S E =45.000 Kn/m2

5.4.7 Variation du la vitesse V4=150m/s et E=101250 Kn/m2 :

Tableau 5..46 : les valeurs de dplacement pour V4=150 m/s et a =0.15g


Profondeur Vitesse priode Acclration Dplacement (m)
0 150 0 1,4715 0
3,6 150 0,096 1,4715 0,00034386
7,4 150 0,197 1,4715 0,00145292
10,8 150 0,288 1,4715 0,00309475
14,4 150 0,384 1,4715 0,00550178
18 150 0,48 1,4715 0,00859654
21,6 150 0,576 1,4715 0,01237901
25,2 150 0,672 1,4715 0,01684921
28,8 150 0,768 1,4715 0,02200713
32,4 150 0,864 1,4715 0,02785277

125
Chapitre 05 : Simulation Numrique

36 150 0,96 1,4715 0,03438614


Tableau 5.47: les valeurs de dplacement pour V4=150 m/s et a =0.20g
Profondeur Vitesse priode Acclration Dplacement (m)
0 150 0 1,962 0
3,6 150 0,096 1,962 0,00045848
7,4 150 0,197 1,962 0,00193723
10,8 150 0,288 1,962 0,00412634
14,4 150 0,384 1,962 0,00733571
18 150 0,48 1,962 0,01146205
21,6 150 0,576 1,962 0,01650535
25,2 150 0,672 1,962 0,02246561
28,8 150 0,768 1,962 0,02934284
32,4 150 0,864 1,962 0,03713703
36 150 0,96 1,962 0,04584819
Tableau 5.48 : les valeurs de dplacement pour V4=150 m/s et a =0.25g
Profondeur Vitesse priode Acclration Dplacement (m)
0 150 0 2,4525 0
3,6 150 0,096 2,4525 0,0005731
7,4 150 0,197 2,4525 0,00242153
10,8 150 0,288 2,4525 0,00515792
14,4 150 0,384 2,4525 0,00916964
18 150 0,48 2,4525 0,01432756
21,6 150 0,576 2,4525 0,02063168
25,2 150 0,672 2,4525 0,02808202
28,8 150 0,768 2,4525 0,03667855
32,4 150 0,864 2,4525 0,04642129
36 150 0,96 2,4525 0,05731024

Tableau 5.49 : les valeurs de dplacement pour V4=150 m/s et a =0.30g


Profondeur Vitesse priode Acclration Dplacement (m)
0 150 0 2,943 0
3,6 150 0,096 2,943 0,00068772
7,4 150 0,197 2,943 0,00290584
10,8 150 0,288 2,943 0,00618951
14,4 150 0,384 2,943 0,01100357
18 150 0,48 2,943 0,01719307
21,6 150 0,576 2,943 0,02475802
25,2 150 0,672 2,943 0,03369842
28,8 150 0,768 2,943 0,04401426
32,4 150 0,864 2,943 0,05570555
36 150 0,96 2,943 0,06877228
126
Chapitre 05 : Simulation Numrique

Variation des dplacements en fonction du la zone


V=150m/S E =101250 Kn/m2
5,5C

5,5F

5,57

5,5E
@9A <5 9E
5,56
@B A <5 B5
5,5D
G@D A <5 BE G
5,5B G@6 A <5 D

5,59

5
5 D,7 F,6 95,C96,6 9C B9,7BE,BBC,CDB,6 D7
4

Figure 5.28 : Variation des dplacements en fonction du la zone


V=150m/S E =101250Kn/m2

5.4.8 Variation du la vitesse V5=200m/s et E=180000 Kn/m2 :


Tableau 5.50 : les valeurs de dplacement pour V5=200 m/s et a =0.15g
Profondeur Vitesse priode Acclration Dplacement (m)
0 200 0 1,4715 0
3,6 200 0,072 1,4715 0,00019342
7,4 200 0,148 1,4715 0,00081727
10,8 200 0,216 1,4715 0,0017408
14,4 200 0,288 1,4715 0,00309475
18 200 0,36 1,4715 0,00483555
21,6 200 0,432 1,4715 0,00696319
25,2 200 0,504 1,4715 0,00947768
127
Chapitre 05 : Simulation Numrique

28,8 200 0,576 1,4715 0,01237901


32,4 200 0,648 1,4715 0,01566719
36 200 0,72 1,4715 0,0193422
Tableau 5.51 : les valeurs de dplacement pour V5=200 m/s et a =0.20g
Profondeur Vitesse priode Acclration Dplacement (m)
0 200 0 1,962 0
3,6 200 0,072 1,962 0,0002579
7,4 200 0,148 1,962 0,00108969
10,8 200 0,216 1,962 0,00232106
14,4 200 0,288 1,962 0,00412634
18 200 0,36 1,962 0,0064474
21,6 200 0,432 1,962 0,00928426
25,2 200 0,504 1,962 0,01263691
28,8 200 0,576 1,962 0,01650535
32,4 200 0,648 1,962 0,02088958
36 200 0,72 1,962 0,02578961
Tableau 5.52 : les valeurs de dplacement pour V5=200 m/s et a =0.25g
Profondeur Vitesse priode Acclration Dplacement (m)
0 200 0 2,4525 0
3,6 200 0,072 2,4525 0,00032237
7,4 200 0,148 2,4525 0,00136211
10,8 200 0,216 2,4525 0,00290133
14,4 200 0,288 2,4525 0,00515792
18 200 0,36 2,4525 0,00805925
21,6 200 0,432 2,4525 0,01160532
25,2 200 0,504 2,4525 0,01579613
28,8 200 0,576 2,4525 0,02063168
32,4 200 0,648 2,4525 0,02611198
36 200 0,72 2,4525 0,03223701

Tableau 5.53 : les valeurs de dplacement pour V5=200 m/s et a =0.30g


Profondeur Vitesse priode Acclration Dplacement (m)
0 200 0 2,943 0
3,6 200 0,072 2,943 0,00038684
7,4 200 0,148 2,943 0,00163454
10,8 200 0,216 2,943 0,0034816
14,4 200 0,288 2,943 0,00618951
18 200 0,36 2,943 0,0096711
21,6 200 0,432 2,943 0,01392639
25,2 200 0,504 2,943 0,01895536
28,8 200 0,576 2,943 0,02475802
128
Chapitre 05 : Simulation Numrique

32,4 200 0,648 2,943 0,03133437


36 200 0,72 2,943 0,03868441

Variation des dplacements en fonction du la zone


V=200m/S E =180000 Kn/m2
5,56E
5,56
5,5DE
5,5D
5,5BE
@9A <5 9E
5,5B
@B A <5 B5
5,59E G@D A <5 BE G
5,59 G@6 A <5 D
5,55E
5
95,C
96,6

B9,7
BE,B
BC,C
DB,6
5
D,7
F,6

9C

D7

Figure 5.29 : Variation des dplacements en fonction du la zone


V=200m/S E =180000Kn/m2

129
Chapitre 05 : Simulation Numrique

Variation des dplacements en fonction du la vitesse


a=0.15g
9,6

9,B

5,C G 9<BE C9
B<E5
5,7
D<955
5,6 6<9E5
5,B E<B55

5
5 95 B5 D5 65

Figure 5.30 : Variation des dplacements en fonction du la vitesse


a=0.15g

Variation des dplacements en fonction du la vitesse


a=0.20g

9,C
9,7
9,6
9,B
G 9<BE C9
9
5,C B<E5

5,7 D<955
5,6 6<9E5
5,B E<B55
5
5 95 B5 D5 65

Figure 5.31 : Variation des dplacements en fonction du la vitesse


a=0.20g

130
Chapitre 05 : Simulation Numrique

Variation des dplacements en fonction du la vitesse


a=0.25g

B,E

9,E G 9<BE C9
B<E5
9
D<955

5,E 6<9E5
E<B55
5
5 95 B5 D5 65

Figure 5.32 : Variation des dplacements en fonction du la vitesse


a=0.25g

Variation des dplacements en fonction du la vitesse


a=0.30g

B,E

9,E G 9<BE C9
B<E5
9
D<955

5,E 6<9E5
E<B55
5
5 95 B5 D5 65

Figure 5.33 : Variation des dplacements en fonction du la vitesse


a=0.30g

131
Chapitre 05 : Simulation Numrique

Tableau 5.54 : variation du moment flchissant MS ; Mp


Vitesse densit EI G E A MS Mp Mp/Ms
25,82 1,8 3333 1200 3000 1,4715 7,35676425 8,62 1,172
50 1,8 3333 4500 11250 1,4715 1,9618038 2,29 1,167
100 1,8 3333 18000 45000 1,4715 0,49045095 0,4125 0,841
150 1,8 3333 40500 101250 1,4715 0,2179782 0,346 1,587
200 1,8 3333 72000 180000 1,4715 0,12261274 0,1625 1,325
25,82 1,8 3333 1200 3000 1,962 9,809019 11,5 1,172
50 1,8 3333 4500 11250 1,962 2,6157384 3,06 1,17
100 1,8 3333 18000 45000 1,962 0,6539346 0,55 0,841
150 1,8 3333 40500 101250 1,962 0,2906376 0,2697 0,928
200 1,8 3333 72000 180000 1,962 0,16348365 0,325 1,988
25,82 1,8 3333 1200 3000 2,4525 12,2612738 14,37 1,172
50 1,8 3333 4500 11250 2,4525 3,269673 3,82 1,168
100 1,8 3333 18000 45000 2,4525 0,81741825 0,687 0,84
150 1,8 3333 40500 101250 2,4525 0,363297 0,44 1,211
200 1,8 3333 72000 180000 2,4525 0,20435456 0,406 1,987
25,82 1,8 3333 1200 3000 2,943 14,7135285 17,25 1,172
50 1,8 3333 4500 11250 2,943 3,9236076 4,59 1,17
100 1,8 3333 18000 45000 2,943 0,9809019 0,825 0,841
150 1,8 3333 40500 101250 2,943 0,4359564 0,366 0,84
200 1,8 3333 72000 180000 2,943 0,24522548 0,488 1,99

moment-zone V=25.81m/s

B5
9C 0
97
96
9B
95
C
7
6
B
5
5 5,E 9 9,E B B,E D D,E
1

Figure5.34 : Variation moments en fonction du la zone


E=3.000 Kn/m2

132
Chapitre 05 : Simulation Numrique

moment-zone V=50m/s
E

6,E
0
6

D,E

B,E

9,E

5,E

5
5 5,E 9 9,E B B,E D D,E
1

Figure 5.35 : Variation des moments en fonction du la zone


V=50m/S E=11250 Kn/m2

moment-zone V=100m/s

9,B

0
9

5,C

5,7

5,6

5,B

5
5 5,E 9 9,E B B,E D D,E

1
Figure 5.36 : Variation des moments en fonction du la zone
V=100m/S E=45.000 Kn/m2
133
Chapitre 05 : Simulation Numrique

moment-zone V=150m/s
5,E
5,6E
0
5,6
5,DE
5,D
5,BE
5,B
5,9E
5,9
5,5E
5
5 5,E 9 9,E B B,E D D,E
1

Figure 5.37 : Variation des moments en fonction du la zone


V=150m/S E=101250 Kn/m2
moment-zone V=200m/s
5,7

0
5,E

5,6

5,D

5,B

5,9

5
5 5,E 9 9,E B B,E D D,E
1
Figure 5.38 : Variation des moments en fonction du la zone
V=200m/S E=180.000 Kn/m2
Pour une zone de faible sismicit de 0.15g le rapport Mp/Ms= 2.36/2.0299=1.166
Pour une zone de moyenne sismicit de 0.25g le rapport Mp/Ms=3.1409/2.0299=1.60
134
Chapitre 05 : Simulation Numrique

Pour une zone de forte sismicit de 0.30g le rapport Mp/Ms= 4.7038/4.704=1


Le moment flchissant diminue quand la vitesse de cisaillement augmente.

moment-vitesse A=0,15g
95
H
C
F
7
E
6
D
B
9
5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.39 : Variation des moments en fonction du la vitesse a=0.15g

moment-vitesse A=0,20g
96

9B

95

5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.40 : Variation des moments en fonction du la vitesse


a=0.20g

135
Chapitre 05 : Simulation Numrique

moment-vitesse A=0,25g
97

96

9B

95

5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.41: Variation des moments en fonction du la vitesse


a=0.25g

Moment-Vitesse A=0,3g
B5

9C

97

96

9B

95

5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.42 : Variation des moments en fonction du la vitesse


a=0.30g

136
Chapitre 05 : Simulation Numrique

5.5 Formulation propose :


A partir des rsultats de la simulation numrique ou propose les formules suivantes pour corriger
les moments par la mthode simplifie :
4 4AIY ~ ~ T#
W
R U W
(5.5)

4 4AIY # ~ ~ T# ##
WU
R U W
(5.6)

4 4AIY ## (5.7)
Avec
A : coefficient dacclration de la zone
V : vitesse (m/s)
Les rsultats obtenus sont rsumes dans le tab5.55
Tableau 5.55 : les valeurs de moments flchissant corrige par les formules proposes :

Module de
Vitesse acclration Mp Mcorr
dformation
25,82 3000 1,4715 7,35670098 8,32494266
50 11250 1,4715 1,9618038 2,4618038
100 45000 1,4715 0,49045095 0,34045095
150 101250 1,4715 0,2179782 0,2179782
200 180000 1,4715 0,12261274 0,12261274
25,82 3000 1,962 9,80893464 11,0999235
50 11250 1,962 2,6157384 3,28240507
100 45000 1,962 0,6539346 0,4539346
150 101250 1,962 0,2906376 0,2906376
200 180000 1,962 0,16348365 0,16348365
25,82 3000 2,4525 12,2611683 13,8749044
50 11250 2,4525 3,269673 4,10300633
100 45000 2,4525 0,81741825 0,56741825
150 101250 2,4525 0,363297 0,363297
200 180000 2,4525 0,20435456 0,20435456
25,82 3000 2,943 14,713402 16,6498853
50 11250 2,943 3,9236076 4,9236076
100 45000 2,943 0,9809019 0,6809019
150 101250 2,943 0,4359564 0,4359564
200 180000 2,943 0,24522548 0,24522548
Figure 5.43 : Variation des moments en fonction du la vitesse
a=0.15g (par la formule proposes)

137
Chapitre 05 : Simulation Numrique

moment-vitesse A=0,15g
95

5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.43 : Variation des moments en fonction du la vitesse


a=0.15g (par la formule proposes)

moment-vitesse A=0,20g
96

9B

95

5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.44 : Variation des moments en fonction du la vitesse


a=0.20g (par la formule proposes)

138
Chapitre 05 : Simulation Numrique

moment-vitesse A=0,25g
97

96

9B

95

5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.45 : Variation des moments en fonction du la vitesse


a=0.25g (par la formule proposes)

moment-vitesse A=0,30g
B5

9C

97

96

9B

95

5
5 BE E5 FE 955 9BE 9E5 9FE B55 BBE

Figure 5.46: Variation des moments en fonction du la vitesse


a=0.30g (par la formule proposes)
La formulation propose donne de bons rsultats.

139
Conclusion gnrale

Conclusion Gnrale

1-Introduction
Lobjet de notre travail a porte sur ltude de la rponse dynamique des pieux en zone sismique,
en utilisant le code lments finis (Plaxis) et une mthode simplifie
Cette simulation concerne ltude de linfluence de diffrents paramtres sur leffort internes :
moment flchissant
Les paramtres choisit sont :
Sismicit du la zone
Caractristiques de sol (module de dformation et la densit)
Vitesse de cisaillement

2 Conclusions :
Cette comparaison permet dmettre ce qui suit
Les rsultats obtenus avec la mthode simplifie sont assez concordance avec ceux du calcul par
lment finis.
En effet obtient :
Pour une zone de faible sismicit de 0.15g le rapport Mp/Ms= 0,96375/1,10351=0.873
Pour une zone de moyenne sismicit de 0.25g le rapport Mp/Ms=0.531/0.66121=0.802
Pour une zone de forte sismicit de 0.30g le rapport Mp/Ms= 0.327/0.36781=0.889

1- Le dplacement et tte du pieu augmente quand lacclration augmente et diminue quand


la vitesse de cisaillement augmente
2- Le moment flchissant augmente quand lacclration augmente et diminue quand le
module de dformation de sol augment

140
Conclusion gnrale

En fin a partir des rsultats obtenus par la simulation numrique (limite) et une formulation
simplifie propose en fonction de la vitesse de cisaillement et la zone sismique pour corriger les
rsultats de la mthode simplifie ; cette formulation a donne des bons rsultats

3 Recommandations :
La mthode simplifie peut tre utilise a condition de corriger le moment obtenu en fonction de
la vitesse de cisaillement et de la zone sismique pour ce la les formulations propose donne de
bons rsultats et peut tre utilis.

141
Rfrences Bibliographique

REFERENCES BIBLIOGRAPHIQUES

Polous : h(1971) Beharviour of laterally loaded piles 1-single pile journal of Geotechnical
engineering ;A.S.CE, vol101 nSM5 , PP 733-751
Tarzaghi .K,(1955) Evalution of the coefficient of sub grade reaction , Geotechnique vol 5.Pp
297.326
Mnard .L , Bourdon G et Gambin .M (1969): method gnerale de pressiometrique , sols-soils -
n22-23 vol VI PP 16-29
Matlock .H & Reese LC (1960) Generalized soulotions for lateral piles journal. SMFE ASCE
.SM5 octobre 1960 PP63-91
Frank(1999) calcul foundations superficielles et profondes , editions Technique de lingenieur et
Presse de ponts et chausses 141P
Li Yan & byrne (1992) lateral pile reponse of a monotonic head loading Canadian Geotechnical
journal n29-PP955-970
Kouda & al (1998) : direct measurement of P-Y relation Ships of piles in sand , centrifufe 98-
Vol 1. PP 551-556
Hadjadji : (1993) analyse du comportement de pieux sous chargement horizontale Thse doctorat
de lcole Nationale des ponts et chausses , 307P
Wikler (1867) , Die leher vonder elastzitat and Festigkeit , verlag 182P
Novak , (1991) Ples under dynamic load , Proc, Second Geotechnical Earthquake Engineering
and soil dynamics , March 11-15 St Louis Missouri PP 2433-2456
AFNOR (1993): Essai Statique de pieu isol sous effort transversal Norme Francais P94-151
Foray. P., Genevois, J.M., Labanieh, S. & Goulois, A. (1989). Effet de la mise en place sur la
capacit portante des pieux dans les sables. Compte-rendu du douzime Congrs International
de Mcanique des sols et des travaux de Fondations. Rio de Janeiro, Vol. 2, pp. 913-914.
Lizzi, F. (1982a). The pali radice (root pile). A state-of-the-art report. Proceedings of
symposium on Recent Development in ground Improvement Techniques, Bangkok, pp. 417-432.
Lizzi, F. (1982b). The static restoration of monuments. Sagep Publisher, Genova, 146 pages.
Poulos, H.G. (1987). Analysis of residual stress effects in piles. Journal of Geotechnical
Engineering Division, ASCE, Vol. 113, N 3, pp. 216-229.
Puech, A., Foray, P., Boulon, M. & Desrues, J. (1979). Calcul des pieux l'arrachement partir
d'un modle numrique en contraintes effectives-premiers rsultats. Proceeding of the 7th
European Conference on Soils Mechanics and Foundation Engineering, Brighton, Vol. 1, pp.
227-233.
Shakhirev, V., Magnan, J.P. & Ejjaaouani, H. (1996). Etude exprimentale du comportement du
sol lors du fonage des pieux. Bulletin de liaison des laboratoires des Ponts et Chausses, N :
206, pp. 95-116.

142
Rfrences Bibliographique

Schlosser, F., et Guilloux A. (1981). Le frottement dans le renforcement des sols. Revue
franaise de gotechnique , N16, pp 65-67.
Robinsky, E.I. & Morrison, C.F. (1964). Sand displacement and compaction around model
friction piles. Revue Canadienne de Gotechnique, Vol 1, N : 2, pp. 81-93.
Mokrani, L. (1991). Simulation physique du comportement des pieux grande profondeur en
chambre dtalonnage. Thse de Doctorat de linstitut National Polytechnique de Grenoble
Lehane, B.M. (1992). Experimental investigations of displacement pile behaviour using
instrumented field piles. PhD Thesis, Imperial College of Science, Technology & medicine
Chow, F. (1995). Field measurements of stress interactions between displacement piles in sand.
Ground Engineering, July/August 1995, pp. 36-40.
Boulon, M., Plytas, C. & Foray, P. (1986). Comportement d'interface et prvision du frottement
latral le long des pieux et tirants d'ancrage. Revue franaise de gotechnique. N 35, pp. 31-
48.
Briaud, J.L., Tucker, L.M. & E. Ng. (1989). Axially loaded 5 pile group and a single pile in sand.
Comptes-rendus du douzime Congrs International de Mcanique des Sols et des travaux de
Fondations, Rio de Janeiro, Vol. 2, pp. 1121-1124.
Canou, J. (1989). Contribution ltude et lvaluation des proprits de liqufaction des
sables. Thse de doctorat de lEcole Nationale des Ponts et Chausses
PHRI : Port and Harbour research instate (1980) Technical Facilities in Japon , Bureau of ports
and Harbous Port and harbor research institute Ministry ob transport 317.P
Fasicule N62 ttre v (1993) regle technique de conception et de calcul des foundations des
ouvrages de genie civil cahier des clauses techniques grrales applicable aux marches
publics de travuax . Ministre de lEquipement du logement et de transport
[Council [2] de 1978, MEXICO[11] TRIFUNAC et BRADY[10], NEWMARK [6], SEAOC[7],
SEED et DONOVAN, DOBRY, OWEIS et URZUA dans article essentiel [3] paru en 1976.
ZEEVAERT [11] en 1972, SCHNABLE en 1973, SEAOC [7] en 1974, CHEN en 1971, URZUA
en 1973, DOBRY [3],] SOULOUNIC , mthode simplifie des pieux en zone sismique ,
sminaire international de gnie civil, universit de Constantine (1985)

143
Annexe

Annexe

144
Annexe

Prsentation de plaxis
1-Bref historique
Lhistoire de La mthode des lments finis a commenc en 1819, quand Navier dfinit
une mthode dtude des systmes hyperstatiques bass sur lapplication des conditions
dquilibre et de compatibilit, puis Maxwell en 1864 et Castigliano en 1878 tablissent dune
faon complte les thormes de lnergie. Au dbut de 20me sicle, des rsultats fondamentaux
dans le domaine des mthodes dapproximation en ont tabli les bases en effectuant une
discrtisation spatiale du domaine tout en utilisant les mthodes dapproximation variationnelles.
En 1955, Argyris prsente une mthode unifie des mthodes de dplacements et des forces,
lanne suivante Turner et Clough publient une prsentation systmatique de la mthode des
dplacements. Ces deux publications sont particulirement importantes et prsentent
vritablement le dbut de la MEF comme technique de calcul des structures complexes.

2- Concept de base
La MEF consiste remplacer la structure physique tudier par un nombre fini
dlments ou de composants discrets qui reprsentent un maillage. Ces lments sont lis entre
eux par un nombre de points appels nuds. On considre le comportement de chaque partie
indpendante, puis on assemble ces parties de telle sorte quon assure lquilibre des forces et la
compatibilit des dplacements rels de la structure en tant quobjet continu.
La MEF est extrmement puissante puisquelle permet dtudier correctement des
structures continues ayant des proprits gomtriques et des conditions de charges compliques.
Elle ncessite un grand nombre de calculs qui, cause de leur nature rptitive, sadaptent
parfaitement la programmation numrique.

3 - Calcul par la MEF


La mthode des lments finis est thoriquement, la plus satisfaisante puisquelle
permet la modlisation des problmes gotechniques complexes.
Elle ncessite :

145
Annexe

1. La dfinition de la gomtrie du problme, pour que les frontires de calcul


ninfluencent pas sur les rsultats.
2. Le choix dune loi de comportement du sol, de type Mohr-Coulomb, Ducan etc.
3. Les caractristiques des ouvrages et des lments dinterface pour introduire
linteraction sol-structures et les conditions hydrauliques.
4. Les caractristiques des boutons et des tirants dancrage.
5. Ltat initial des contraintes et des pressions interstitielles.
Elle permet
1. deffectuer les calculs dcoulement ;
2. de simuler toutes phases de travaux ;
3. de prendre en compte les variations des caractristiques : des ouvrages ; des couches de
terrain et de lits de boutons ou de tirants ;
4. de calculer un Cfficient de scurit.

Les sorties fournissent


Pour les ouvrages :
1. Les dplacements des ouvrages,
2. les sollicitations internes des ouvrages,
3. les efforts dans les boutons ou tirants.
Et pour le sol :
4. les dplacements du sol,
5. les dformations du sol,
6. les contraintes totales et effectives dans le sol,
7. les pressions interstitielles,
8. de calculer un coefficient de scurit.
4- Prsentation de PLAXIS
4-1 Introduction
PLAXIS est un programme dlments finis en deux dimensions spcialement conu pour
raliser des analyses de dformation et de stabilit pour diffrents types dapplications
gotechniques. Les situations relles peuvent tre reprsentes par un modle plan ou
axisymtrique. Le programme utilise une interface graphique pratique permettant aux utilisateurs
de gnrer rapidement un modle gomtrique et un maillage dlments finis bass sur la coupe
146
Annexe

verticale de louvrage tudier. Les utilisateurs sont supposs tre capables de travailler dans un
environnement Windows.
Lutilisation de PLAXIS consiste en quatre sous-programmes (Input, Calculations,
Output, Curves) :
a) Le programme dentre de donnes (Input)
Le programme contient tout ce qui est ncessaire pour crer et modifier un modle
gomtrique, pour gnrer le maillage dlments finis correspondant et pour gnrer les
conditions initiales.

b) Le programme de calcul (Calculations)


Ce programme contient tous les lments pour dfinir et amorcer un calcul par la mthode
des lments finis. Au dbut du programme de calcul, lutilisateur doit choisir le projet pour
lequel les calculs vont tre dfinis.

c) Le programme de rsultats (Output)


Ce programme contient tous les lments qui permettent de voir les rsultats des donnes
gnres et des calculs dlments finis. Au dbut du programme de rsultats, lutilisateur doit
choisir le modle et la phase de calcul approprie ou le numro du pas pour lequel les rsultats
seront affichs.

d) Le programme courbes (Curves)


Ce programme contient tout le ncessaire pour gnrer des courbes charges-
dplacements, des chemins de contrainte et des courbes contraintes-dformations.

4.2 Les modles de comportement utiliss dans PLAXIS


Les modles de comportement de sols sont trs nombreux : depuis le modle lastique-
plastique de Mohr-Coulomb jusquaux lois de comportement les plus sophistiques permettant de
dcrire presque tous les aspects du comportement lasto-visco-plastique des sols, aussi bien sous
sollicitation monotone que cyclique. Ces modles ont t dvelopps dans le but dtre intgrs
dans des calculs par lments finis. Dans ce schma, la modlisation par lment finis permet de
rsoudre le problme aux limites en tenant compte, par une loi de comportement raliste, du
comportement rel du sol. Deux difficults majeures ont empch la ralisation complte de ce
147
Annexe

schma : dune part les lois de comportement qui dcrivent bien le comportement des sols sont
complexes et demandent, pour la dtermination des paramtres quelles contiennent, des tudes
spcifiques lourdes sortant du cadre des projets dingnierie. La validation des lois de
comportement a fait lobjet, dans les annes 80 de plusieurs ateliers pour comparer les rponses
des diffrents modles sur diffrents chemins de sollicitation. La seconde difficult a t
lintgration de ces lois de comportement dans des codes par lments finis, bi ou
tridimensionnels. Peu de codes sont oprationnels actuellement, avec des lois sophistiques. Le
cot de ces calculs est gnralement important [BOU1, 04].
Pour traiter un problme de soutnement (paroi moule, palplanche, etc.), il est tout
fait adapt de considrer le sol comme lastoplastique et le modle de Mohr-Coulomb sera bien
adapt dans ce cas. Mais pour traiter un remblai sur sols mous, avec chargement par tapes et
consolidation, il faut tenir compte de lcrouissage. Le matriau se consolide et il est plus adapt
dutiliser le Soft Soil Model qui prend en compte cette volution du matriau. Pour un calcul
dcoulement, il suffit de prendre un matriau lastique, mais on peut avoir coupler coulement
et dformation ; dans ce cas un modle lastoplastique peut tre justifi.
Les rgles dor dans le domaine de la simulation du comportement dun ouvrage sont :
Quel est le comportement principal modliser ?
Utiliser un modle qui dcrive ce comportement ;
Interprter les rsultats, notamment en fonction des paramtres de la modlisation.
En ce sens, la modlisation numrique ne fournit sous une autre forme que les donnes du
problme pos.

148

Vous aimerez peut-être aussi