RakotondrabeTojonirinaJeanC ESPA MAST 19

UNIVERSITE D’ANTANANARIVO
-------------------
ECOLE SUPERIEURE POLYTECHNIQUE

D’ANTANANARIVO
-------------------
DOMAINE : SCIENCES DE L’INGENIEUR
-------------------
MENTION : INGENIERIE MINIERE
Mémoire de fin d’étude pour l’obtention du diplôme d’ingénieur Polytechnicien
Grade Master
Parcours : Sciences et Techniques Minières
Intitulé :
MODELISATION NUMERIQUE DU
CREUSEMENT D’UN TUNNEL ROUTIER
RELIANT AMPANTSAKANA ET ANALAKELY
DANS LA COMMUNE URBAINE
D’ANTANANARIVO
Soutenu le 27 Avril 2019
par
RAKOTONDRABE TOJONIRINA Jean Claudio
Encadreur : Monsieur Rémi Roger FABIEN
Enseignant Chercheur à l’Ecole Supérieure Polytechnique d’Antananarivo
ANNEE UNIVERSITAIRE : 2017-2018

UNIVERSITE D’ANTANANARIVO
-------------------
ECOLE SUPERIEURE POLYTECHNIQUE

D’ANTANANARIVO
-------------------
DOMAINE : SCIENCES DE L’INGENIEUR
-------------------
MENTION : INGENIERIE MINIERE
Mémoire de fin d’étude pour l’obtention du diplôme d’ingénieur Polytechnicien
Grade Master
Parcours : Sciences et Techniques Minières
Intitulé :
MODELISATION NUMERIQUE DU
CREUSEMENT D’UN TUNNEL ROUTIER
RELIANT AMPANTSAKANA ET ANALAKELY
DANS LA COMMUNE URBAINE
D’ANTANANARIVO
Soutenu le 27 Avril 2019
par
RAKOTONDRABE TOJONIRINA Jean Claudio
Président du jury : Monsieur Léon Félix RANAIVOSON

Rapporteur : Monsieur Rémi FABIEN
Examinateurs : Monsieur Joachin Clotaire DAMY
Monsieur Rasamoelina Naina ANDRIATSITOMANARIVOMANJAKA
ANNEE UNIVERSITAIRE : 2017-2018

FISAORANA
Voaloahany indrindra dia isaorako i Jesosy Kristy Tompo noho ny nanafiany hery, sy ny
namonjeny ary ny nanolorany saina ahy tamin’ny fanatanterahana ity asa ity.
Manaraka izany dia tsy hay ny tsy hisaotra an’Andriamatoa isany izay voalaza anarana eto
ambany :
 Andriamatoa Rijalalaina RAKOTOSAONA, talen’ny Sekoly Ambony

Polyteknikan’Antananarivo ;
 Andriamatoa Léon Felix RANAIVOSON, mpampianatra mpikaroka sy lehiben’ny “Mention
Ingénierie Minière”, eto anivon’ny Sekoly Ambony Politeknikan’Antananarivo;
 Andriamatoa Rémi Roger FABIEN, mpampianatra mpikaroka eto anivon’ny Sekoly Ambony
Politeknikan’Antananarivo, tamin’ny torohevitra rehetra nomeny ary indrindra tamin’ny
fotoana natokany ho ahy na dia be aza ny andraikitra nanoloana azy;
 Andriamatoa Joachin Clotaire DAMY, mpampianatra mpikaroka eto anivon’ny Sekoly
Ambony Politeknikan’Antananarivo;
 Andriamatoa Rasamoelina Naina ANDRIATSITOMANARIVOMANJAKA,
mpampianatra mpikaroka eto anivon’ny Sekoly Ambony Politeknikan’Antananarivo;
 Andriamatoa Ando Malala RAZEFASON, Chef de Service Ouvrage d’Art ao amin’ny
Ministeran’ny Fanajariana ny Tany sy ny Asa Vaventy;
 ireo mpampianatra rehetra ao amin’ny « Mention Ingénierie Minière » eto amn’ny Sekoly
Ambony Politeknikan’Antananarivo;
Hotolorana ny fankasitrahana manokana koa ireto olona manaraka ireto:
 ireo namana aman-tsakaiza sy mpiara-mianatra izay nanampy sy nifanakalo hevitra tamin’ny

fanatanterahana ity asa ity ;
 Faniry RAZAFIMAHATRATRA ;
 Dada sy Neny ary Gabriel izay tsy mitsahatra ny manohana ahy.
˜ Mankasitraka anareo rehetra tompoko˜

i
Sommaire
FISAORANA……………………………………………………………………………….i
SOMMAIRE……………………………………………………………………………….ii
LISTE DES TABLEAUX………………………………………………………………...iii
LISTE DES FIGURES…………………………………………………………………..iv
LISTE DES ABREVIATIONS……………………………………………………….viii
LISTE DES UNITES…………………………………………………………………..x
LISTE DES SYMBOLES……………………………………………………………….xi
INTRODUCTION………………………………………………………………………..01
CHAPITRE I GENERALITES SUR LES TUNNELS…………………………………03
CHAPITRE II METHODE DE CREUSEMENT DES TUNNELS……………………19
CHAPITRE III LA METHODE DES ELEMENTS FINIS……………………………38
CHAPITRE IV LOCALISATION DE LA ZONE D’ETUDE…………………………53
CHAPITRE V UTILISATATION DE PLAXIS 2D 8.6……………………………...69
CONCLUSION…………………………………...………………………………………85
BIBLIOGRAPHIE……………………………………………………………………….87
WEBOGRAPHIE………………………………………………………………………...90
TABLE DES MATIERES…………………………………………………………….…91
ANNEXE…………………………………………………………………………………...I
ii
Liste des tableaux
Tableau 1 Quelques exemples de classification des roches et la hauteur de terrain décomprimé au-
dessus d’une cavité [12] ............................................................................................................... 5
Tableau 2 Classification des roches suivant le R.Q.D [12] .................................................................. 5
Tableau 3 Classification des roches en fonction de la résistance à la compression[3] ........................ 7
Tableau 4 Choix du type de soutènement 1974) [3] ............................................................................ 7
Tableau 5 Les modules de rigidité de différents types de soutènement [12] ...................................... 13
Tableau 6 Essais en laboratoire appliqué aux sols et rochers [12] ..................................................... 16
Tableau 7 Les caractéristiques des éléments massifs bidimensionnels ............................................ 44
Tableau 8 Exemples illustrant les significations de chaque matrice [38] ........................................... 51
Tableau 9 Les dimensions du tunnel ................................................................................................ 68
Tableau 10 Propriétés géotechniques des matériaux [17],[31],[21],[27], [11] [39] .......................................... 76
Tableau 11 Caractéristiques des infrastructures sus-jacentes [21],[22]................................................. 77
Tableau 12 Résumé de quelques valeurs des paramètres importants ............................................... 80
Tableau 13 Résumé des valeurs de quelques paramètres ................................................................. 81
Tableau 14 Résumé des valeurs de quelques paramètres ................................................................. 82
Tableau 15 Valeur du tassement maximal la méthode de Sagaseta (1980) ..................................... 83
Tableau 16 Valeur du tassement maximal la méthode de Poupelloz (1984) ................................... 83
iii
Liste des figures
Figure 1 Zone décomprimée au-dessus d’une cavité [12] .................................................................... 4
Figure 2 Configuration géologique illustrant les limites du RQD [10] ................................................ 6
Figure 3 Type de soutènement en fonction du RQD d’après LAUFFER [03] ..................................... 6
Figure 4 Les contraintes et les déplacements [12] ................................................................................ 8
Figure 5 Courbes des contraintes et du déplacement radial en fonction de la distance r [17] .............. 9
Figure 6 Portée non soutenue d'un tunnel [11] ................................................................................... 10
Figure 7 Variation du taux de déconfinement, de la pression et du déplacement radial en fonction
de la distance du front de taille [11]........................................................................................... 11
Figure 8 Front de taille stable (extrusion négligeable) [11]................................................................ 12
Figure 9 Front de taille stable, convergence et extrusion forte [11] ................................................... 12
Figure 10 Figure 10 Front de taille instable [11] ................................................................................ 12
Figure 11 Ligne caractéristique de la cavité [11] ............................................................................... 14
Figure 12 Profil gaussien de la cuvette de surface ........................................................................... 15
Figure 13 Une galerie de reconnaissance [12] .................................................................................... 18
Figure 14 Tunnel creusé en pleine section [21] .................................................................................. 20
Figure 15 Méthode de creusement en section divisée- Cloison double (à gauche), Cloison simple (à
droite) [14] ................................................................................................................................. 20
Figure 16 Creusement en demi-section [14] ....................................................................................... 21
Figure 17 Cycle d’opération de l’abattage à l’explosif [28] ............................................................... 22
Figure 18 Le bouchon pyramidal ..................................................................................................... 23
Figure 19 Exemple d'un plan de tir séquentiel - Pleine section (Tunnel de la Chamoise) [17] .......... 23
Figure 20 Schéma d'un trou de mine ................................................................................................ 25
Figure 21 Détonateurs électriques à microretard de 25 ms .............................................................. 25
Figure 22 Coupe d'un détonateur non électrique [06] ........................................................................ 26
Figure 23 Détonateur électronique [06] .............................................................................................. 26
iv
Figure 24 Bras d'une machine à attaque radiale [20] .......................................................................... 27
Figure 25 Bras d'une machine à attaque axiale [20] ........................................................................... 27
Figure 26 Machines à attaque global [35 ........................................................................................... 28
Figure 27 Graphique d'organisation ................................................................................................. 28
Figure 28 Tête en étoile [20] .............................................................................................................. 30
Figure 29 Tête à ouverture en périphérique [20] ................................................................................ 31
Figure 30 Tête à ouverture radiale [20] .............................................................................................. 31
Figure 31 Coupe d'un tunnelier à confinement par pression de boue .............................................. 32
Figure 32 Vis extraction (en rouge) [35] ............................................................................................ 33
Figure 33 Erecteur à voussoir [20] ..................................................................................................... 33
Figure 34 Exemple de boulons d’ancrage [17] ................................................................................... 34
Figure 35 Tubes perforés [17] ............................................................................................................ 35
Figure 36 Treillis soudés et de béton projeté [17] .............................................................................. 35
Figure 37 Cintres lourds [12].............................................................................................................. 36
Figure 38 Cintres coulissants [12] ...................................................................................................... 36
Figure 39 Figure 36 Discrétisation du domaine D en éléments triangulaires [01] ............................. 42
Figure 40 Discrétisation du domaine D............................................................................................ 44
Figure 41 Eléments de massif bidimensionnel [37]............................................................................ 47
Figure 42 Les éléments de massif tridimensionnels les plus utilisés [18] .......................................... 48
Figure 43 Exemple de maillage en 2D et en 3D [12] ......................................................................... 48
Figure 44 Les trois étapes du calcul : prétraitement, résolution et post-traitement [13]..................... 51
Figure 46 Tunnel Hubert Grabit en 1930 ......................................................................................... 53
Figure 46 Tunnel d’Ambanidia ........................................................................................................ 53
Figure 47 Série stratoïde dans la partie Nord d’Antananarivo (d’après Nédélec et Al.,1994) [36].... 54
Figure 48 Schéma géologique d’Antananarivo[02] (numérisé et modifié par TOJONIRINA Jean
Claudio) ................................................................................................................................... 55
Figure 49 Carte géologique améliorée ............................................................................................. 56
v
Figure 50 Coupe géologique suivant B-B’ [02] ............................................................................ 57
Figure 51 Carte de structure tectonique d’Antananarivo ................................................................. 57
Figure 52 Profil géoélectrique du profil II ....................................................................................... 59
Figure 53 Profil géoélectrique du profil I......................................................................................... 60
Figure 54 Modèle lithologique suivant le tracé du tunnel (Source : RAHAMBININTSOA Sedra,
2018) ........................................................................................................................................ 62
Figure 55 Coupe géologique longitudinale prévisionnelle du tracé ................................................. 63
Figure 56 Les deux tracés du tunnel (Fond : Google Earth) ............................................................ 66
Figure 57 Tracé 1 (Axe Ampantsakàna-Antsahavola) ..................................................................... 66
Figure 58 Tracé 2 (Axe Ampantsakàna - Totoha-baton'Antaninarenina) ........................................ 67
Figure 59 Section visible en forme de fer à cheval du tunnel routier [05] ................................... 67
Figure 60 Les contenus de l’onglet « project » [26] ........................................................................... 71
Figure 61 Fenêtre de la création géométrique [26] ............................................................................. 71
Figure 62 Relation entre la pression sur le contour du tunnel et le tassement de surface [12] ........... 75
Figure 63 Géométrie du model de la section 1................................................................................. 78
Figure 66 Déplacement horizontal (Valeur maximale : 1,26 . 10-3 m) ............................................. 79
Figure 67 Déplacement vertical (déplacement maximal : -8,09. 10-3 m) ......................................... 80
Figure 68 Maille déformée (déplacement maximal : -8,09. 10-3 m)................................................. 80
Figure 69 Maille déformée (déformation maximale : 305,6.10-6 m) ................................................ 81
Figure 70 Déformation des mailles (Déplacement maximal : 140,88 .10-3 m.................................. 81
vi
Liste des abréviations
2D : dimension 2
3D : dimension 3
AFTES : Association Française des Tunnels et de l’Espace Souterrain
B : larguer de la galerie
d : diamètre d’un trou de mines
D : Diamètre du tunnel
𝑫𝒊 : sous domaine pour la création d’un nombre fini d’éléments
E : Module d’Young
E : Est
EA : Rigidité normale
EI : Rigidité de flexion
E.D.P : Equation aux Dérivée Partielles
EPBS : Earth Pressure Balance Shield
Hp : hauteur de la surcharge du terrain uniformément répartie suivant l’horizontal
Ht : la hauteur de la galerie
H.S.M: Hardening Soil Model
K : coefficient de non uniformité d’émanation de gaz
Ks : Coefficient de rigidité du soutènement
KNO3 : nitrate de potassium ou salpêtre
L : Longueur non soutenue du tunnel
LCPC : Laboratoire Central des Ponts et Chaussées
M.E.F : Méthode des Eléments Finis
M.E.D : méthode des éléments distincts
M.M.C : mécanique des milieux continus
Ms: Moment agissant sur les soutènements
MC : Modèle de Mohr-Coulomb
N : effort axial agissant sur le soutènement ;
pi : pression fictive
𝑷𝒋 : production moyenne journalière
ℚ : débit d’air nécessaire au niveau du front d’abattage
𝒒𝒓 : débit relatif de CH4 ou de CO2
vii
R.Q.D : Rock Quality Designation
R : rayon du tunnel
Rc : résistance à la compression de la roche
Rt : résistance à la traction de la roche
S : Sud
Sfront: Tassement au front de taille
Sfinal : Tassement final
𝑺 : surface de la section visible d’une galerie
S.S.M : Modèle pour sols mous ou Soft Soil Model
S.S.C.M : Modèle pour sols mous avec effet du temps ou Soft Soil Creep Model
τ : résistance au cisaillement
Ts: force de tangentielle agissant sur les soutènements
u :déplacement radial
Ui : déplacement radial initial
𝑼(𝒓) ∶ déplacement radial
v : déplacement tangentiel
𝒗𝒂 ∶ la vitesse admissible de l’air à travers l’ouvrage
Vi : déplacement tangentiel initial
𝑽𝒕 : déplacement tangentiel
viii
Liste des unités
cm : centimètre
Ga : milliard d’années
km2 : kilomètre carré
kN/m : kilonewton par mètre
kN/m2 : kilonewton par mètre carré
kN/m3 : kilonewton par mètre cube
kNm : kilonewton mètre
kPa :kilopascal
m2 : mètre carré
m/jour : mètre par jour
m3/s : mètre cube par seconde
m/s : mètre par seconde
mm : millimètre
ms : milliseconde
Ma : million d’années
MPa : mégapascal
[°] : degré
°C : degré Celsius
ix
Liste des symboles
% : pourcentage
c : cohésion
𝛔𝟎 : contrainte initiale dans le massif
𝝈𝒄 : la résistance à la compression simple
γ : poids volumique (non saturé)
 𝒔𝒂𝒕 : poids volumique (saturé)
kx :perméabilité horizontale
ky: perméabilité horizontale
𝝂 : coefficient de Poisson
π : pi
𝝁 : rendement de l’explosif
𝛌 ∶ le taux de déconfinement
 : angle de frottement interne
ρ : masse volumique
ψ : l’angle de dilatance
𝝈𝒌𝒋 : composante de tenseur de contrainte
τ :la résistance au cisaillement
[𝒌𝒆 ] : matrice de rigidité élémentaire
KO : coefficient des terres au repos
[K] : matrice de rigidité globale
{F} : vecteur contenant les vecteurs représentant les actions extérieures sur l’élément fini
{𝐟𝐞 } ∶ vecteur élémentaire des forces appliquées aux nœuds
{𝐮𝐞 }: vecteur élémentaire des déplacements nodaux
U∞ : déplacement radial final du tunnel
𝝏Ω𝟎 : morceau de frontière du système à étudier
x
INTRODUCTION
Antananarivo, la capitale de Madagascar, est la partie que l’on appelle aussi les
hautes terres par rapport à l’ensemble de l’île. Elle ne cesse de s’améliorer surtout en termes
de construction des infrastructures. Cet avancement a remarquablement commencé depuis
la période de la colonisation pendant laquelle des voies ferrées et des routes ont été
construites et des tunnels, comme le tunnel d’Ambanidia et d’Ambohidahy, ont été creusés.
La capitale de la Grande Ile connait actuellement de nombreux problèmes sur la
circulation à cause des embouteillages. Les routes ont besoin d’entretiens surtout dans le
centre-ville et ses périphéries. Pour pallier à cela, nombreux sont les solutions déjà proposées
par les techniciens, à savoir l’élargissement des routes, la construction des nouvelles routes
liant Analakely et autres endroits de la ville.
La solution que nous allons préconiser dans cette étude est le creusement d’un
tunnel routier liant Analakely et Ampantsakàna où l’intensité de la circulation s’avère
croissante ces dernières années. Vue l’étroitesse de la Ville à cause des infrastructures déjà
existantes, la réalisation de cet ouvrage souterrain serait une solution efficace pour remédier
à l’encombrement qui est un fléau dans cette zone d’influence. La réalisation d’un tel projet
engendre, en revanche, des perturbations plus ou moins importantes dans l’équilibre interne
du massif environnant et surtout dans les ouvrages sus-jacents. D’ailleurs, pendant la
réalisation des travaux, les ingénieurs doivent assurer la stabilité de l’ouvrage et contrôler
les déformations du massif encaissant afin d’éviter tout interaction avec les constructions
avoisinantes.
D’où le problématique : « Quels seraient les effets du creusement de ce tunnel sur
les infrastructures existantes ? » Ce travail s’intitule « Modélisation numérique du
creusement d’un tunnel routier reliant Ampantsakàna et Analakely dans la commune urbaine
d’Antananarivo ». L’hypothèse de notre recherche consiste à évaluer, en faisant des calculs
se basant sur les méthodes empiriques et une modélisation numérique au moyen du logiciel
Plaxis 2D version 8.6 (version d’essai), le tassement dû au creusement de ce tunnel. En outre,
cette recherche a pour but de faire des études bibliographiques afin d’avoir plus
d’information sur les théories et les méthodes permettant d’aboutir aux résultats attendus,
particulièrement la géologie de la zone d’étude.
Cette recherche comporte cinq (05) grands chapitres. Dans un premier temps, le
chapitre I sera consacré sur les généralités sur les tunnels, puis nous allons nous focaliser sur
1
les rappels théoriques des méthodes de dimensionnement des tunnels. Ensuite, nous allons
élucider dans le troisième chapitre en ce qui concerne la Méthode des Eléments Finis suivi
de la présentation et la localisation de la zone d’étude pour le creusement du tunnel dans le
chapitre se focalisant sur la méthodologie. Enfin, nous parlerons dans les deux derniers
chapitres la présentation du logiciel Plaxis 8.5 2D et son utilisation pour la modélisation de
l’effet du creusement d’un tunnel routier, entre Analakely et Ampantsakàna, sur les
infrastructures existantes.
2
CHAPITRE I GENERALITES SUR LES TUNNELS
CHAPITRE I
GENERALITES SUR LES TUNNELS
Chapitre I. GENERALITES SUR LES TUNNELS
I.1 Introduction
Un tunnel désigne toute voie de communication ou chaussée situé sous un ouvrage
de couverture qui, quel que soit son mode de construction, crée un espace confiné. Dans le
domaine minier, les tunnels permettent le passage d’un endroit en un autre, alors que dans le
domaine géotechnique les tunnels sont construits afin d’assurer la circulation des véhicules et
des piétons dans de bonnes conditions de sécurité et de confort. Le choix de la forme d’un
tunnel dépend de plusieurs facteurs, à savoir la dureté de la roche encaissante, la durée
d’utilisation de l’ouvrage, etc.
I.2 Classification des massifs rocheux et de dimensionnement des

tunnels
On peut classifier les méthodes de calcul de dimensionnement des tunnels par les
trois (03) méthodes suivantes :
 les méthodes empiriques ;
 les méthodes analytiques ;
 les méthodes numériques.
I.2.1 Méthodes empiriques
Pour les études préliminaires, le dimensionnement de tunnels se base sur la

classification des massifs rocheux. Puisque pendant des années, l’art des tunnels reste une
science essentiellement empirique et c’est pour cela que les experts tirent des conclusions
basées sur leurs expériences passées en faisant des modélisations mathématiques. Nombreuses
sont les méthodes empiriques utilisées pour le dimensionnement des tunnels mais on ne cite
que :
 la méthode de Terzaghi ;
 le R.Q.D ;
 la classification de l’AFTES;
 la méthode de Bieniawski .
3
Ce qu’on va élucider dans les sous titres suivants sont la méthode de Terzaghi, le
R.Q.D et la classification de l’AFTES.
I.2.1.1 La méthode de Terzaghi
En vue d’estimer les charges agissant sur les cintres en acier dans les tunnels, une
classification simple des massifs rocheux a été proposée par Terzaghi en 1946. Cette
estimation de charge, donnée par la relation ci-dessous, suppose qu’une certaine quantité de
roche ou de sol se décomprime et pèsent sur soutènement :
𝑯𝒑 = 𝑲 (𝑩 + 𝑯𝒕 ) (3)
Avec
𝑯𝒑 est la hauteur de la surcharge du terrain uniformément répartie suivant l’horizontal ;
B désigne larguer de la galerie ;
Ht est la hauteur de la galerie ;
K est un coefficient multiplicateur dont la valeur dépend de la nature de terrain rencontré.
Figure 1 Zone décomprimée au-dessus d’une cavité [12]
Il est à noter que cette méthode est applicable aux tunnels à moyenne profondeur
pour lesquels la hauteur de couverture H est supérieur à et 1,5 B + Ht. Le tableau suivant nous
montre quelques exemples de classification des roches faite par Terzaghi, qui n’est valable
que pour les tunnels de section rectangulaire [12].
4
Tableau 1 Quelques exemples de classification des roches et la hauteur de terrain

décomprimé au-dessus d’une cavité [12]
Nature de la roche Charge Hp Remarque

Dure et intacte 0 à 0,25 B Quelques ancrages s’il y a
chute de pierres
Dure et stratifiée 0 à 0,5 B Soutènement léger
Massive avec quelques joints 0 à 0,25B La charge peut changer
brusquement d’un point à
l’autre
Très ébouleuse 0,25à0,35(H+Ht) Peu ou plus de pression
latérale
I.2.1.2 Rock Quality Designation ou R.Q.D
Le R.Q. D est une méthode proposée par D. Deer en 1964, basant sur l’observation
des échantillons prélevés dans des sondages carottés et qui est calculé par la longueur de passe
de sondage.
𝟏𝟎𝟎 × 𝒔𝒐𝒎𝒎𝒆 𝒅𝒆𝒔 𝒍𝒐𝒏𝒈𝒆𝒖𝒆𝒓𝒔 𝒅𝒆 𝒄𝒂𝒓𝒐𝒕𝒕𝒆𝒔 𝒔𝒖𝒑é𝒓𝒊𝒆𝒖𝒓𝒔 à 𝟏𝟎𝒄𝒎
𝑹𝑸𝑫 = (𝟒)
𝑳𝒐𝒏𝒈𝒖𝒆𝒖𝒓 𝒅𝒆 𝒍𝒂 𝒑𝒂𝒔𝒔𝒆 𝒅𝒆 𝒄𝒂𝒓𝒐𝒕𝒕𝒂𝒈𝒆
Pour que les morceaux de roches saines ne soient pas cassés lors du prélèvement il
est essentiel de choisir un carottier de diamètre supérieur à 50mm. La classification des roches
suivant sa qualité par le R.Q.D est donnée par le tableau suivant :
Tableau 2 Classification des roches suivant le R.Q.D [12]
R.Q. D Qualité de roche

R.Q. D<25% Très mauvaise
25%< R.Q.D <50% Mauvaise
50%< R.Q.D <75% Moyenne
75%< R.Q.D <90% Bonne
90%< R.Q.D <100% Très bonne
5
Il faut faire attention sur le fait que cet indicateur puisse être trompeur : pour un
même massif stratifié, on pourra conclure à un RQD de 0% ou de 100% au même endroit
selon l’inclinaison du forage. Pour pallier à ce problème, il est important de repérer la direction
et l’inclinaison des forages, puis de corréler les résultats entre eux.
Figure 2 Configuration géologique illustrant les limites du RQD [10]
Le type de soutènement en fonction des valeurs du RQD est donné par la figure
suivante :
Figure 3 Type de soutènement en fonction du RQD d’après LAUFFER [03]
I.2.1.3 Classification de l’AFTES
A.F.T.E.S ou Association Française des Tunnels et de l’Espace Souterrain a été créée

en 1972 pour répondre aux recommandations de la Conférence International sur les travaux
6
souterrains tenue à Washington en 1970 [12]

. La principale activité de l’A.F.T.E.S consiste à
diffuser des recommandations techniques sur les travaux souterrains. Ces recommandations
se reposent sur le choix du soutènement et du revêtement en fonction de la résistance à la
compression 𝜎𝑐 en Méga Pascal [MPa].
Tableau 3 Classification des roches en fonction de la résistance à la
compression[3]
Le tableau suivant présente les choix de soutènement en fonction de la résistance à la

compression de la roche, selon l’AFTES.
Tableau 4 Choix du type de soutènement 1974) [3]
7
I.2.2 Méthodes analytiques
Les méthode analytiques se basent sur des hypothèses visant à simplifier la

modélisation du problème rencontré.Ces hypothèses de simplifacation concernent :
 la géométrie : la section du tunnel est supposée circulaire et d’axe horizontal ;
 la stratigraphie : on ne prend en considération qu’une seule couche de terrain supposé
homogène ;
 la loi de comportement du sol est supposée élastique ou élasto-plastique ;
 l’état initial des contraintes est supposé isotrope et homogène ;
 les formules sont éxprimées dans le plan et dans le cadre des petites déformations [04].
I.2.2.1 Les méthodes analytiques élastiques
Pour cette méthode, l’entourage de l’excavation est traité comme un milieu élastique.
Les étapes susmentionnées sont utilisées pour la plupart des méthodes proposées par des
experts comme Lamé (1852), Kirsch (1898) et Goudman-shi (1982). Le calcul se fait en
suivant les étapes suivantes :
 calcul des contraintes (σr, σθ, τrθ) et les déformations dans le sol dus aux champs
des contraintes pour pouvoir passer au calcul des déplacements initiaux Ui et Vi ;
Figure 4 Les contraintes et les déplacements [12]
 calcul de la contrainte totale et le champ de déplacement supplémentaire dans le sol et

des contraintes de contact à l’interface sol-soutènement ;
 calcul des forces (Ms,Ts) internes dans le soutènement.
A travers ces étapes, on a abouti à l’obtention des équations finales donnant le

déplacement radial u, le déplacement tangentiel v, effort axial N et moment M dans le
8
soutènement. En effet, un tunnel de section circulaire non soutenu dans le cas linéaire, le
déplacement radial en fonction de la distance r est donné par la relation suivante :
𝟏+𝝂 𝑹𝟐
𝑼(𝒓) = − 𝒑 (5)
𝑬 𝒓
et
𝟏+𝝂
𝑽𝒕 = 𝟐 𝒑𝝅𝑹𝟐 (6)
𝑬
Avec 𝑝 = σv = σh
Figure 5 Courbes des contraintes et du déplacement radial en fonction de la

distance r [17]
I.2.2.2 Les méthodes analytiques élasto-plastiques
Ces sont des méthodes qui prennent en compte du développement d’une zone
plastique autour de la cavité souterraine due à l’excavation du tunnel. Parmi ces méthodes, on
peut citer :
a. La Théorème limite
C’est une méthode proposée par Caquot, Atkinson et Mühlhaus à partir des années
1950. L’objectif est d'étudier les comportements des tunnels à basse profondeur. Par exemple,
Mühlhaus (1985) a proposé une solution statiquement admissible par la longueur non soutenue
L et de diamètre D du front de taille, définissant une sphère de surface S1 comprise dans une
autre sphère de surface S2.
9
Figure 6 Portée non soutenue d'un tunnel [11]
En considérant les contraintes à l'intérieur de ce volume comme à l'état plastique, les

conditions limitées à la surface (σs= σr), Mühlhaus a pu obtenir la relation suivante :
𝑯
(𝟏+𝟐 )𝟐
𝑫
𝑳 = 𝑫√ 𝟏 (7)
𝛔
[𝟏+(𝛌𝒑 −𝟏) 𝒔 ]𝛌𝒆−𝟏
𝛔𝒖
Avec
λ𝑒 = 1 + 𝑠𝑖𝑛𝜑 ⁄1 − 𝑠𝑖𝑛𝜑
𝛔𝒔 est la surcharge au niveau de la surface
𝛔𝒖 est la résistance à la compression non confinée du sol
Le principal défaut de ces premières approches est qu’elles ne tiennent pas en compte
des déformations du terrain qui interviennent avant la mise en place du soutènement. Alors,
pour pallier à cela, on recommande l’utilisation de la méthode de convergence-confinement.
b. La méthode de convergence-confinement
C’est une méthode conçue par Panet, pour simplifier le problème 3D du creusement
du tunnel en un problème 2D de déformation plane de l’interaction terrain-soutènement. Pour
l’estimation des déformations du massif (convergence, tassement) et la validation du
soutènement, cette méthode est couramment introduite dans les calculs numériques 2D. La
convergence est liée à un déplacement.
Le confinement est la pression radiale qui s’applique sur le pourtour de l’excavation,
en présence d’un soutènement. En effet, elle consiste à simuler l’avancement du front par la
diminution d’une pression fictive (pi), appelée pression de confinement, appliquée aux parois
de l’excavation [11]. Cette pression est donnée par la relation suivante :
𝒑𝒊 = (𝟏 − 𝛌)𝛔𝟎 (8)
10
Avec 𝛔𝟎 la contrainte initiale dans le massif

𝛌 le taux de déconfinement variant de 0 (pour une section loin à l’avant du front du taille)
à 1 (pour une section loin à l’arrière). 𝛌=0 correspond à l’état initial en avant du front à une
distance suffisante pour négliger son influence.
Figure 7 Variation du taux de déconfinement, de la pression et du déplacement

radial en fonction de la distance du front de taille [11]
Il est à noter que la convergence des tunnels est la variation de la distance entre deux
points situés sur la paroi des tunnels. Elle est en fonction de la distance x (figure ci-dessus),
entre la section de mesure et le front de taille, et du temps t qui s’est écoulé depuis le passage
du front de taille au niveau de la section de mesure, de la distance non soutenue d derrière le
front de taille et de la rigidité du soutènement Ks [12].
11
Pour la mesure de l’extrusion du terrain en avant du front de taille, une méthode appelée
méthode de ADECO-RS a été proposée par Lunardi (1993). Trois situations (présentées par
les figures ci-dessous) peuvent exister :
 le front de taille est stable et l’extrusion au front de taille est faible ;
 le front de taille est stable mais une extrusion notable est liée aux déformations
en avant du front de taille ;
 le front de taille est instable et s’éboule.
Le domaine d’application de la méthode de convergence-confinement est constitué

par les deux premières situations.
Figure 8 Front de taille stable (extrusion négligeable) [11]
Figure 9 Front de taille stable, convergence et extrusion forte [11]
Figure 10 Figure 10 Front de taille instable [11]
12
Dans le cas classique, le soutènement est mis en place à une certaine distance d du
front de taille, appelée distance non soutenue, pour laquelle un déplacement du tunnel
ud(d) s’est déjà produit. La valeur λ du taux de déconfinement correspondant au déplacement
ud est λd. Ce paramètre (λd) dépend essentiellement de la distance non soutenue à l’arrière du
front de taille et aussi de la loi de comportement du terrain et la rigidité du soutènement. Pour
un massif dont le comportement est élastique, il peut être calculé par la relation suivante :
𝟎,𝟕𝟓𝑹 𝟐
𝛌𝒅 = 𝟏 − 𝟎, 𝟕𝟓 [𝟎,𝟕𝟓+𝒅] (9)
La pression de soutènement, notée Ps, croit avec la rigidité du soutènement mais elle
est limitée par la résistance de ce dernier. Ainsi, dans une galerie circulaire de rayon R, on
définit le module de rigidité normale au soutènement KSN par l’expression :
𝒖𝒓
𝑷𝒔 = 𝑲𝑺𝑵 (10)
𝑹
Avec ur est le déplacement radial.

Puisque dans le cas axisymétrique, ce module définit la rigidité de soutènement mais
lorsque ses conditions ne sont pas satisfaites, on passe à la détermination du module en flexion
KSF. Les modules de rigidité pour chaque type de soutènement sont donnés par le tableau
suivant :
Tableau 5 Les modules de rigidité de différents types de soutènement [12]
Type de soutènement Module de flexion

Anneau circulaire d’épaisseur constante e 𝑬 𝒆
𝑲𝑺𝑵 = 𝟏−𝛎𝑺 𝟐 × 𝑹 et 𝑲𝑺𝑭 = 𝟏−𝛎𝑺 𝟐 × 𝑹
𝑬 𝟏
𝑺 𝑺
(e<< R)
Cintres métalliques circulaire espacés 𝑲𝑺𝑵 = 𝑬𝒂 𝑨 et 𝑲𝑺𝑭 = 𝒔𝑹𝒂𝟑 𝑰
𝑬
𝒔𝑹
appliqués contre le terrain
Soutènement par boulons 𝟏 𝑺𝒕 𝑺𝒊 𝟒𝑳
= [𝛑𝚽𝟐 𝑬 + 𝑸] et 𝑲𝑺𝑭 = 𝟎
𝑲𝑺𝑵 𝑹 𝒃
Avec
L : est la longueur des boulons.

A est l’aire de la section du cintre
Φ : est le diamètre des boulons.
𝑬𝒂 le module d’Young de l’acier 𝑬𝒃 : est le module d’Young du matériau constituant la tige du boulon.
I : est le moment d’inertie du profilé. Q est un paramètre permettant de prendre en compte des déformations
qui se produisent au niveau de l’ancrage et de tête de boulon.
13
La méthode convergence-confinement peut être aussi représentée par la ligne

caractéristique de la cavité, décrivant la relation existante entre la pression fictive pi et le
déplacement radial ur.
Figure 11 Ligne caractéristique de la cavité [11]
La courbe suit le segment [AC] si le massif reste dans un domaine élastique. Par
contre, si le taux de déconfinement est augmenté, la frontière du domaine élastique est atteinte
à partir du point B et d’un λ > 𝛌é𝒍𝒂𝒔𝒕𝒊𝒒𝒖𝒆 .Ainsi, la courbe n’est plus linéaire et le déplacement
final sera élevé (cas du tunnel non soutenu). On peut ajouter dans ce même graphe la ligne de
confinement du soutènement et l’intersection des deux courbes représente l’équilibre du
cavité-soutènement [11]. En correspondance de ce point, on trouve le déplacement radial final
du tunnel (U∞) et la charge finale reprise par le soutènement (péquilibre).
I.2.3 Méthodes numériques
Ces méthodes se sont apparues dans les années 70. Elles font appel aux modèles de
comportement du sol pour la simulation des chemins de sollicitation que l’on rencontre dans
le cas des problèmes réels. Son objectif est de modéliser le déconfinement provoqué par le
creusement des tunnels. Elles peuvent être utilisées en 2D ou 3D. Ces méthodes se sont basées
sur la méthode des éléments finis qu’on va élucider dans le chapitre 03. Grâce à l’existence
de l’ordinateur performant qui facilite les calculs dans un temps très court, cette méthode est
très utilisée actuellement.
I.2.4 Méthode de prédiction du tassement
Il existe trois méthodes permettant la prédiction du tassement engendré par le

creusement du tunnel :
14
 les méthode empirique, qui se basent sur de nombreuses mesures détaillées

(modèle réelle, chatiers,etc)
Ces méthodes permettent, à partir d’un faible nombre de paramètres, d’obtenir les tassements
en tenant compte de la dimension géométrique de l’ouvrage, du type de sol et de la localisation
du niveau des aquifères.
Prenons comme exemple la méthode de Marcos (1958) a appliqué une fonction gaussienne
pour décrire le profil de la cuvette de tassement [40]. Des mesures de tassement ont été faites
par Peck (1968) dans 14 tunnels creusés dans des sols plastiques saturés et ont démontré
l’allure de la courbe de Gauss.
Figure 12 Profil gaussien de la cuvette de surface
La formulation analytique se présente comme suit :

−𝒙𝟐
𝑺(𝒙) = 𝑺𝒎𝒂𝒙 𝒆𝒙𝒑[ ] (11)
𝟐𝒊
Avec
x est la distance au centre de la cuvette
S(x) représente le tassement vertical à l’abscisse x
Smax est le tassement maximum en surface
i désigne l’abscisse du point à partir de l’axe de symétrie vertical du tunnel.
L’éstimation de i et de Smax est résumé dans un tableau à l’annexe 4-f.
 les méthodes analytiques qui utilisent les formulations analytiques basant

sur des équations ;
15
 les méthodes numériques.
I.3 Les étapes d’investigations avant le creusement d’un tunnel

Il est indispensable d’effectuer une reconnaissance géologique et géotechnique du
site hôte de l'ouvrage pour avoir une bonne connaissance des formations existantes. Cela a
pour but d’éviter l’effondrement du tunnel pendant ou après son creusement. Ainsi, pour que
l’ouvrage creusé réponde aux besoins des usagers en toute sécurité, on doit faire :
 des études géologiques et géotechniques qui permettent de connaitre les roches,
les évènements géologiques qui se sont apparues dans la zone d’étude. Cela
permet aussi de bien connaitre l’hydrogéologie du site.
 des essais in situ, à travers lesquels on peut avoir une idée générale des
paramètres tels que la puissance des couches, l’état des roches en profondeur,
etc.
On peut regrouper les essais in situ comme suit :

 les méthodes géophysiques, comme le sondage électrique, la méthode
sismique ;
 les sondages carottés, permettant d’avoir des échantillons des roches qu’on
peut analyser en laboratoire, et nous donnent les puissances des formations
sous-jacentes ;
 essais d’eau (déterminant les perméabilités, le niveau statique le de l’eau,…) ,
tels que l’essai Dupuit, Lefranc, Lugeon ;
Le tableau suivant nous résume quelques essais en laboratoires avec leurs buts
respectifs :
Tableau 6 Essais en laboratoire appliqué aux sols et rochers [12]
Nom de l’essai A analyser Objectif de l’essai

Analyse minéralogique Sols et roches Déterminer la composition
minéralogique de la roche
Essai au bleu de méthylène Sol Déterminer la nature minéralogique des

particules
argileuses et leur influence sur le
16
comportement
global du sol.
Essai brésilien Roches Déterminer la résistance à la traction 𝜎𝑡

Essai de cisaillement à la Sols, roches Déterminer la résistance au cisaillement
τ
boite de Casagrande
Essai de compression Sols, roches Déterminer la résistance à la

compression simple 𝜎𝑐
simple
Essai œdométrique Sols Déterminer le coefficient de pression des
terres au
repos K0
Essai triaxial Roches Déterminer la résistance au cisaillement

τ
Limite d’Atterberg Sols Déterminer l'aptitude de l'argile à être
plus ou moins
plastique
Mesure de densité Sols, Roches Déterminer le poids volumique γ
A part les essais susmentionnés, on utilise aussi les galeries de reconnaissances

(galerie de dimension réduite), qui permettent d’évaluer aux mieux les difficultés sur la
réalisation de l’ouvrage projeté. Ce dernier est utilisé surtout dans le cas où le creusement se
fait par la méthode conventionnelle ou par des machines à attaque ponctuelle.
Une galerie de reconnaissance (figure 11) a pour but de :
 contribuer directement au creusement du tunnel un traitement de terrain depuis
la galerie pilote, drainage des eaux d’infiltration, alésage successif pour des
tunnels de grande section ;
17
Figure 13 Une galerie de reconnaissance [12]
 faciliter l’exécution des travaux, comme la réduction des brusques venues
d’eau ;
 avoir des informations sur les terrains traversés : reconnaissance géologique et
géotechnique permettant d’optimiser les méthodes de creusement .
CONCLUSION
Nous avons pu voir dans ce chapitre qu’il existe trois méthodes de dimensionnement
des soutènements des tunnels. La première se base sur des formules empiriques proposées par
des experts. L’utilisation de ces dernières permet de classer les roches dans lesquelles le tunnel
est creusé. Il est vital de mentionner qu’il faut les utiliser avec précaution quand on est dans
l’étude de tassement causé par le creusement. Ensuite, la deuxième méthode qui est la méthode
analytique permet de simplifier la modélisation du problème rencontré, concernant la
géométrie, la stratigraphie (une seule couche homogène), la loi de comportement ainsi que
l’état initial des contraintes. Enfin, la méthode numérique, qui se base sur la méthode des
éléments finis, permet de faciliter le calcul parce que son utilisation fait appel à des ordinateurs
performants.
En guise de conclusion, les méthodes de calcul utilisées pour le dimensionnement
des ouvrages souterrains reposent sur des considérations théoriques expérimentales ou
empiriques. Dans le cas théorique, la théorie du calcul à la rupture est l’outil de base pour
l’étude de la stabilité du front de taille. Par contre, les méthodes en déplacement sont les plus
utilisées pour le dimensionnement du soutènement. Dans le chapitre suivant, nous allons
parler des différents types de méthodes de creusement des tunnels.
18
CHAPITRE II METHODES DE CREUSEMENT DES TUNNELS
CHAPITRE II
METHODES DE CREUSEMENT
DES TUNNELS
Chapitre II. METHODES DE CREUSEMENT DES TUNNELS
Ce chapitre se focalise sur les différents types de méthodes de creusement des

tunnels. Nous allons les expliquer une à une avec les paramètres à prendre en compte ainsi
que les outils respectifs à utiliser, permettant d’effectuer le creusement.
Le creusement d’un tunnel nécessite la prise en compte des paramètres suivants, qui
interviennent sur le choix des procédés des constructions et la définition du profil en travers :
 le volume utile du tunnel (partie fonctionnel) , qui conditionne l’utilisation du tunnel
à savoir la circulation des moyens de transports ( wagon, train, voitures,…) , les
espaces pour les conduites transitant l’ouvrage[32] ;
 l’environnement, qui concerne l’existence des infrastructures, des sites naturelles,…
qui sont sensibles aux vibrations lors du creusement et sans oublier la présence ou non
d’un aquifère sur la zone d’exécution du creusement ;
 le terrain encaissant l’ouvrage : cela englobe toutes les propriétés physiques et
chimiques du massif, qui sont connues à travers des essais en laboratoire après des
prélèvement des échantillons des roches sur terrain. Ces paramètres conditionnent les
méthodes de creusement, la section du tunnel ainsi que le type de soutènement
approprié.
II.1 Les procédées de creusement

Voici les différents types de procédés utilisés pour le creusement des tunnels selon
des paramètres à prendre en considération.
II.1.1 Creusement à pleine section
C’est une méthode consistant à excaver la totalité de la section du tunnel en une seule
fois [14]. Elle nécessite l’emploi d’un matériel ayant une grande hauteur de l’excavation. Elle
est couramment utilisée pour le creusement des tunnels dans des roches dures (en utilisant des
tunneliers ou Tunnel Boring Machine) ou dans des sols meubles (en utilisant des tunneliers à
boucliers). L’avantage de cette méthode est qu’elle est performante et elle permet des
avancements très importants, à peu près 12 m par jour (qui varie en fonction du modèle du
tunnelier et de la nature des roches).
19
Figure 14 Tunnel creusé en pleine section [21]
II.1.2 Creusement à section divisée
Cette méthode consiste en l’ouverture préalable d’une ou plusieurs galeries de petite

section dont une partie de soutènement participe au soutènement de
l’excavation totale à réaliser [14]
. Avec cette technique chaque phase des travaux comprend
l’excavation des terrains en sections réduites.
Figure 15 Méthode de creusement en section divisée- Cloison double (à gauche),

Cloison simple (à droite) [14]
II.1.3 Creusement par demi-section
Avec cette méthode, la première phase consiste à commencer le creusement par la

partie supérieure (ou stross) de la section avec une longueur variant de 5 à 6 m selon le cas,
mais elle peut être utilisée très largement en faisant varier la longueur de l'avancement
unitaire. Ensuite, on procède à l’excavation de la demi-section inférieure. Après chaque phase
d’avancement, la mise en place de soutènement est effectuée. Le creusement par demi-section
20
est utilisé dans le cas où le terrain est hétérogène et nécessite l’emploi d’un soutènement
important [28]
. L’inconvénient de cette méthode est qu’elle est lente dans le cas où le terrain
est médiocre car la pose des soutènements nécessite une durée importante.
Figure 16 Creusement en demi-section [14]
II.2 Abattage à l’explosif

C’est une méthode utilisée dans les roches dures et très dures et s’effectue pour
chaque volée d’avancement de manière cyclique selon les opérations suivantes [17] :
 perforation des trous de mines par le biais des jumbo perforatrice;
 bourrage ;
 le tir ;
 l’aérage et le purgeage;
 le marinage ;
 le soutènement.
Le traçage du plan de tir galerie est exécuté sur la base d'un plan de tir théorique
adapté aux conditions particulières du front de taille (fracturation, hétérogénéité, etc.), et à la
qualité des roches rencontrées. Pour l’amélioration de la qualité du traçage et la réduction des
délais nécessaires à cette opération, le plan de tir sera élaboré dans un logiciel de perforation
[14]
automatique au moyen d'un ordinateur embarqué sur le jumbo de foration . La figure
suivante nous montre le cycle d’opération d’abattage à l’explosif.
21
Figure 17 Cycle d’opération de l’abattage à l’explosif [28]
Pour le creusement d’un tunnel à l’aide de l’abattage à l’explosif, il est primordial de

choisir un bon plan de tir pour bien contrôler la fissuration du massif et donc de découper le
rocher conformément à la section prédéfinie du tunnel. La section d’un tunnel à abattre
comporte en général trois zones :
 la zone de bouchon, qui est localisée au centre du front de taille et permettant
d’obtenir une fragmentation des roches. Le type de bouchon recommandé dans
cette zone est le bouchon pyramidal (figure ci-desous). A part ce type de bouchon,
il existe également les bouchons de cisaillement qui sont utilisés pour l’abattage
des roches friables (qui se désagrègent facilement).
 la zone d’abattage (tout autour de la zone de bouchon), qui permet de dégager
la majorité de la section à abattre. Les trous de mine dans cette
zone sont soit horizontaux ou légèrement inclinés.
 la zone de contour, avec des trous inclinés dont la disposition permet de bien
cadrer le profil du tunnel.
22
Figure 18 Le bouchon pyramidal
Figure 19 Exemple d'un plan de tir séquentiel - Pleine section (Tunnel de la

Chamoise) [17]
Il est à noter qu’il y a une liaison entre l’efficacité de l’explosifs et sa capacité de

transférer de l’énergie en une courte durée. Toutefois, cette capacité de transfert d’énergie
vers la roche dépend des paramètres suivants [07] :
 les paramètres de la roche ;
 les paramètres de la charge ;
 les paramètres de l’explosif.
23
II.2.1 Les paramètres de la roche
Ce sont ses propriétés physiques (densité), mécaniques et élastiques (vitesse de

propagation des ondes, la résistance à la traction et à la compression). Plus la roche est dense,
plus la consommation en explosif augmente. D’ailleurs, de leur possibilité d’abattage à
l’explosif les roches sont classées parfois par sa résistance à la compression et à la traction.
En effet, la susceptibilité S à l’abattage à l’explosif de la roche est donnée par la relation
suivante :
𝐑𝐜
𝑺= (12)
𝐑𝐭
Où 𝐑 𝐜 désigne la résistance à la compression de la roche ;

𝐑 𝐭 désigne la résistance à la traction de la roche ;
La valeur de ce rapport est comprise entre 10 et 100.
II.2.2 Les paramètres de la charge
Ce sont le type de détonateur à utiliser, le diamètre et la longueur du trou, la distance

entre l’explosif et la paroi du trou. Le diamètre d du trou est calculé en fonction du nombre
total de trou et du volume de la charge explosive. Il est donné par la formule empirique
suivante :
𝑽
d = 𝟏, 𝟏𝟐𝟖√𝑵 (13)
𝒕 𝒍𝒕𝒆
Avec
V désigne le volume total de la charge explosif et donné par la relation :
V = 𝑵𝒕 . 𝒍𝒕𝒆 . 𝑺𝒕 (14)
𝒍𝒕𝒆 désigne la longueur effective du trou, elle est donnée par la formule suivante :
𝒍𝒕𝒆 = 𝒍𝒕 − 𝒍𝒃 (15)
𝒍𝒕 désigne la longueur totale du trou. Elle peut être calculée par la relation :
𝑽𝒊𝒎𝒑
𝒍𝒕 = 𝒋 (16)
𝝁. 𝒏𝒑 . 𝒏𝒎
𝒋
𝒍𝒗
Avec 𝝁 est le rendement de l’explosif. 𝝁 = 𝒍𝒕
𝒍𝒃 est la longueur de bourrage

𝑺𝒕 désigne la section du trou, avec
𝑺𝒕 = 𝟎, 𝟕𝟖𝟓 × 𝒅𝟐 (17)
24
Figure 20 Schéma d'un trou de mine
En ce qui concerne le détonateur, les deux catégories qui existent sont :

 les types normaux : l’explosion suivante ne se fait qu’après la précédente ;
 les types à microretard : l’explosion suivante a lieu après 20 à 30 millisecondes de
l’explosion de la précédente, par exemple le détonateur électrique.
A part cela, voici les différents types de détonateurs les plus utilisés :
 le détonateur électrique : C’est un tube en aluminium avec une tête d’allumage
électrique.
Figure 21 Détonateurs électriques à microretard de 25 ms
 le détonateur non électrique : contrairement aux détonateurs électriques, le détonateur

non-électrique n’a pas besoin d’un exploseur pour produire du courant pour faire
exploser le détonateur.
25
Figure 22 Coupe d'un détonateur non électrique [06]
 le détonateur électronique : c’est un détonateur constituant des composantes

électroniques et utilisé pour pallier aux problèmes des détonateurs électriques
(explosion en cas de courant vagabond ou surtout l’orage). Ils disposent d’une mise à
feu programmable.
Figure 23 Détonateur électronique [06]
II.2.3 Les paramètres de l’explosif
Les paramètres à prendre en compte ici sont la vitesse de détonation, la densité de

l’explosif, la pression de l’explosion, l’énergie produit.
II.3 Méthodes de creusement mécanisées
II.3.1 Les machines à attaque ponctuelle
Ce sont des machines de creusement des tunnels utilisées dans le cas où le terrain
encaissant est composé de roches tendres (marnes, craie, schiste…), dont la résistance à la
compression simple prend une valeur de quelques dizaines de MPa. Son principe de
fonctionnement est basé sur le fait qu’un bras articulé et télescopique muni d’une tête (fraise)
gratte et abatte les roches du front d’attaque. Ces machines se divisent en deux variantes :
26
 les machines à attaque radiale : la tête fraiseuse tourne autour d’un axe en
[16]
prolongement du bras (figure 23) . Pour chaque volée, le creusement
commence par un forage du front de taille pour que la fraise puisse aisément
pénétrer dans les roches. Puis on passe au fraisage progressif de la paroi du
front ;
 les machines à attaque axiale : la tête appelée tambours tourne autour d’un axe
perpendiculaire au bras (figure 24). Le creusement se fait de façon tangentielle
à la surface du front d’attaque.
En ce qui concerne le marinage, le déblai est ramassé par une chaine métallique à
raclettes qui les évacue vers les wagonnettes à travers un tapis à chaine.
Figure 24 Bras d'une machine à Figure 25 Bras d'une

attaque radiale [20] machine à attaque axiale [20]
II.3.2 Les machines à attaque globale
Ce sont des machines qui creusent la totalité de la section de l’ouvrage en une seule
opération à travers une tête circulaire rotative (figure 24). L’effort de poussée sur cette
machine est mobilisé par des « grippeurs » ou patins de vérins latéraux qui prennent appui
radicalement sur le parement naturel de l’excavation [9] ou sur les voussoirs déjà mis en place.
La tête de forage sera immobile quand le grippeur se replie. Cette technique est rapide et sûre,
parce que ces tunneliers avancent en moyenne de 12 à 15 m par jour. De plus, ils permettent
d’effectuer le creusement sans nuisance en surface et ils sont munis d’un système d’orientation
Laser permettant d’enregistrer le mouvement de la machine et de transmettre ces données à
l’ordinateur à travers lesquelles le conducteur sait où aller (vers le haut ou vers le bas) en
suivant les formations (roches) où le creusement doit être fait.
27
Il est à noter que toutes les taches à exécuter en utilisant ces machines (le creusement,
le marinage ainsi que le pose des voussoirs, etc.) se font de manière continue, pas comme dans
le cas du creusement des tunnels par la méthode classique. On peut alors présenter cette
continuité à travers le graphique d’organisation suivant :
Figure 26 Machines à attaque

global [35 Figure 27 Graphique d'organisation
II.3.2.1 Les types de boucliers mécanisés

]
Les boucliers sont des structures cylindriques rigide permettant d’assurer la stabilité
des parois près du front [20]. Il existe deux grandes classes de tunnelier se distinguant par leur
bouclier : les boucliers à front ouvert et les boucliers à front pressurisé.
a. Les boucliers mécanisés à front ouvert
Ces machines ne possèdent pas de système de stabilisation du front de taille.

L’excavation en tête du bouclier peut être partiellement ou totalement mécanisé. Donc, elles
sont utilisées pour le creusement en terrains en bonne tenue, telles que des zones rocheuses
où la venue brusque d’eau n’est pas à craindre. L’abattage s’effectue à la pression
atmosphérique et l’évacuation du déblai se fait par convoyeur à bande.
b. Les boucliers mécanisés à front pressurisé [20]
L’utilisation des boucliers mécanisés précédents sont limités seulement dans des
roches tendres et hors nappes, alors que ce n’est pas toujours le cas. D’où la conception des
boucliers à front confiné pour remédier à ce problème. Ainsi, ces nouvelles variantes de
bouclier rendent possibles au creusement d’un tunnel dans des terrains meubles et dans des
aquifères où le front de taille n’est pas stable. En effet, leur particularité est d’exercer une
pression stabilisatrice sur le front d’attaque. Cette stabilisation peut être assurée par une
28
pression mécanique ou une pression de fluide (air comprimé ou boue bentonitique) ou une
pression de terre du terrain lui-même.
 Confinement à air comprimé
C’est une mode de confinement utilisée dans la première moitié du vingtième siècle
pour le havage puis étendu dans les tunneliers. L’ensemble du tunnel est mis en pression alors
les ouvriers travaillent en atmosphère sous pressions. Ce type de bouclier a été écarté dans le
cas où la perméabilité du terrain traversé est élevée (de l’ordre de 10 -3), dans les alluvions
sablo-graveleuse [30].
 Confinement par pression de terre ou EPBS « Earth Pressure Balance Shield »
Ce type de tunnelier est utilisé dans les terrains meubles et perméables. Les roches
abattues tombent dans la chambre de travaille où elles sont maintenues sous pression. C’est
donc le terrain excavé lui-même qui contrebalance la pression du terrain au front. C’est à
travers la vis sans fin d’Archimède que les déblais se sont évacués vers l’extérieur après avoir
passé par un convoyeur à bande.
 Confinement par pression de boue ou « Slurry Shield »
Le confinement est réalisé en mettant à pression une boue de forage à base de
bentonite qui est amenée par des conduites depuis une station de fabrication et de traitement
à l’extérieur du tunnel. Cette boue, mise en pression de 20 à 30 kPa par rapport à la pression
hydrostatique, assure la stabilisation du front de taille même en terrain meuble et aquifère. La
boue bentonitique est avantageuse parce qu’elle possède la propriété ce qu’on appelle « cake »
(une peau d’argile peu perméable) qui s’essore au contact d’une paroi poreuse et diminue
localement la perméabilité du terrain. De plus, cette dernière facilite le transport et
l’évacuation du déblai jusqu’au centrale de traitement et permet de diminuer les frottements
entre la jupe du tunnelier et le terrain [30].
La roue excavatrice tourne alors dans la chambre remplie de boue et se mélange avec
les déblais. La boue chargée de déblai est extraite au moyen d’une pompe afin de la remplacer
par de boue neuve. Cette machine est munie d’un concasseur afin qu’on ait des dimensions
des matériaux convenables au transport hydraulique.
II.3.2.2 Les principaux organes des tunneliers
En général, les tunneliers sont composés des composantes suivantes :

 les organes de creusement ;
 les organes d’évacuation des déblais ;
29
 les autres équipements.
a. Les organes de creusement
 La tête de coupe
Ses rôles sont de découper les roches au front d’attaque et de ramasser les déblais.
Elle possède trois configurations d’ouverture [20] :
o ouverture périphérique : utilisée dans des roches homogènes. Elle est équipée des
godets qui ramassent les matériaux et qui se vident en partie supérieure à chaque
tour de roue. Ainsi, les matériaux abattus tombent dans une trémie qui alimente en
partie centrale un tapis distributeur.
o ouverture radiale : utilisée dans les terrains meubles qui nécessitent certain
confinement. Les matériaux abattus tombent vers l’arrière dans la chambre et sont
repris soit par une vis soit par les pompes d’extraction.
o tête en étoile : constituée de 6 à 8 bras de faible largeurs.
Figure 28 Tête en étoile [20]
30
Figure 29 Tête à ouverture Figure 30 Tête à ouverture

en périphérique [20] radiale [20]
 Les outils de coupe
Ce sont des outils qui peuvent être installés sur les têtes d’un tunnelier. Dans le cas
des machines à pleine section, les pics (ou dents) ou les disques (ou mollettes) sont souvent
utilisés. En effet, les mollettes sont une solution très efficace dans le cas où le creusement se
fait dans des roches extrêmement dures (Rc dépasse 250 MPa). Par contre, les pics conviennent
pour les terrains meubles ou cohérents. Ces outils de coupe doivent être changés depuis
l’arrière de la tête.
 Appui arrière
Pour vaincre les forces longitudinales (forces exercées sur le front, forces de
frottement du terrain le long du tunnelier à jupe, forces de confinement) qui s’opposent à la
progression du tunnelier, celui-ci doit disposer vers l’arrière d’une réaction d’appuis
suffisante. Les grippeurs sont utilisés dans le cas où l’on voudrait transmettre ces forces au
terrain encaissant, c’est-à-dire que pendant le creusement, la partie arrière du tunnelier est
bloquée au terrain alors que la partie avant se déplace longitudinalement. En revanche, il y a
une autre méthode de transmission des forces susmentionnées : au moyen d’une série de vérin
longitudinale, le tunnelier prend appui à la dernière tranche de revêtement réalisée. Cette
méthode est utilisée pour le cas des roches meuble en présence d’un aquifère.
b. Les organes d’évacuation des déblais
L’utilisation de ces organes dépendent de l’état physique des déblais qui est lié
directement au mode de creusement ainsi que le type de confinement.
31
 Tapis transporteuse
Puisque les machines à front ouvert sont utilisées dans le creusement du tunnel dans
les roches cohérent et non aquifère, donc les déblais ramassés par les godets et déversés dans
la trémie peuvent être amenés vers l’arrière par l’intermédiaire d’une série de tapis
transporteurs. Ensuite, ces tapis transporteurs distribuent les déblais dans les berlines.
 Marinage hydraulique
Le marinage hydraulique est pratiqué dans le cas des tunneliers à confinement de

boue qui travaillent dans des terrains meubles et aquifères où le transport par tapis
transporteuse est impossible. Le transport est assuré par pompage des déblais dilués de boues
dans la chambre afin de les évacuer vers l’extérieur là où elles subissent des traitements.
Figure 31 Coupe d'un tunnelier à confinement par pression de boue
 Marinage par vis d’extraction
Ce type de marinage est utilisé dans les tunneliers à pression de terre ou pression
d’air. Pour réguler le débit en fonction de l’avancement du tunnelier, la solution trouvée par
les constructeurs est l’installation au débouché de la vis à l’air libre d’un appareil régulateur
du volume de matériaux rejetés.
32
Figure 32 Vis extraction (en rouge) [35]
c. Les autres équipements
 Erecteur à voussoirs
Les voussoirs sont des éléments préfabriqués généralement en béton armé, à haute
densité de ferraillage, assurant le rôle de revêtement du tunnel après le creusement. Après le
transport des voussoirs par un convoyeur, l’érecteur (figure 33) assure alors la mise en place
de ces derniers sur la paroi du tunnel.
Figure 33 Erecteur à voussoir [20]
II.4 Les différents types de soutènement

Les objectifs essentiels du soutènement sont de contrôler et d’arrêter les déformations
engendrées par le creusement de l’ouvrage souterrain. Le procédé de soutènement peut se
subdiviser en deux étapes à savoir :
33
 le pré-soutènement ou soutènement provisoire, qui permet soit d’améliorer le terrain

avant l’excavation, donc il joue un rôle de confinement de façon provisoire [17] ;
 le soutènement définitif.
Les équipements les plus souvent rencontrés dans le soutènement ou pré-

soutènement sont :
 les boulons d’ancrage, qui fonctionne soit en mode actif (pour arrêter le
développement des instabilités), soit en mode passif (qui préserve l’ouvrage contre
tout déplacement ou perturbation au sein de l’ouvrage et le massif environnant). On
distingue les boulons à ancrage ponctuel ou à ancrage reparti.
Figure 34 Exemple de boulons d’ancrage [17]
Selon les cas qui peuvent se présenter dans les ouvrages souterrains, les
ancrages ont comme :
o rôle porteur, dans le cas où les massifs présentent des fractures et que des
blocs rochers sont susceptibles de se détacher, on utilise les boulons de
suspension pour empêcher les chutes ;
o rôle confinant, c’est-à-dire limiter l’extension de la zone décomprimée ainsi
que l’amplitude des déformations dans le cas des roches homogènes en
utilisant l’orientation radiale ;
o rôle restructurant, assuré par les boulons à ancrage repartie afin de limiter
tout mouvements susceptibles de désolidariser les blocs rochers dans le cas
où le massif est fracturé.
 les tubes perforés , utilisés dans le pré-soutènement (voûte parapluie) lors des
excavations dans des terrains en mauvaise tenue.
34
Figure 35 Tubes perforés [17]
 les treillis soudés, qui sont utilisés pour renforcer le béton en augmentant sa résistance
au cisaillement.
Figure 36 Treillis soudés et de béton projeté [17]
 les cintres métalliques, qui agissent par supportage. Ils sont en acier et peut être lourds
(possédant une forte inertie) ou légers. Ils jouent à la fois le rôle de protection, de
soutènement et renforcement. En effet, les cintres assurent la protection contre la chute
des blocs isolés, le ralentissement de la convergence des parois et de consolidation de
l’ouvrage [12].
35
Figure 37 Cintres lourds [12]
Les cintres métalliques coulissants sont présentés par la figure 35 ci-dessous.

 le béton projeté, qui contribue dans l’action de confinement.
Figure 38 Cintres coulissants [12]
CONCLUSION
Le creusement d’un tunnel est une opération difficile à exécuter parce qu’il dépend
non seulement de l’utilisation de celui-ci mais également de nombreux paramètres surtout le
type de roches dans lesquelles le creusement sera fait. L’idéal est que les formations à traverser
soient homogènes, assez compactes et imperméables où la venue d’eau n’est pas à craindre.
En revanche, dans la réalité des telles conditions sont rarement rencontrées, voire impossibles
à trouver, d’où l’existence des différents types de méthode de creusement. La méthode
d’abattage par explosif, appelée aussi méthode classique, est utilisée la plus souvent dans les
roches dures. Pourtant, cette dernière n’est pas recommandée si le creusement doit être fait
dans des endroits où les nuisances en surface ne sont pas souhaitables. Pour pallier à la
36
limitation de cette méthode, les ingénieures ont recours, actuellement, à des méthodes
d’abattage mécanique permettant d’avancer plus de dizaine de mètres par jour et pouvant être
utilisées dans des nombreuses conditions, à savoir la traversée dans les roches dures ou dans
un aquifère. A part cela, après avoir fait des études géotechniques et géologiques de l’endroit
où l’on va faire le creusement, il est essentiel d’effectuer une modélisation numérique afin
d’avoir des idées sur les phénomènes réels sur terrain. Cette modélisation est basée sur la
Méthode des Eléments Finis qu’on va entrer à expliquer dans le chapitre suivant.
37
CHAPITRE III LA METHODE DES ELEMENTS FINIS
CHAPITRE III
LA METHODE DES ELEMENTS FINIS
Chapitre III. LA METHODE DES ELEMENTS FINIS
Dans ce chapitre, on va élucider un bref aperçu de la Méthodes des Eléments Finis,

en parlant de ses principes généraux, des différentes phases de calcul ainsi que l’application
de cette dernière au calcul des ouvrages souterrains.
III.1 Historique
L’histoire de la méthode des éléments finis (MEF) a débuté en 1819, quand NAVIER
définit une méthode d’étude des systèmes hyperstatiques basés sur l’application des conditions
d’équilibre et de compatibilité, puis MAXWELL en 1864 et CASTIGLIANO en 1878 établissent
d’une façon complète les théorèmes de l’énergie. Au début du XXième siècle, des résultats
fondamentaux dans le domaine des méthodes d’approximation ont été établis en effectuant une
discrétisation spatiale du domaine tout en utilisant les méthodes d’approximation variationnelles.
En 1955, ARGYRIS présente une approche unifiée des méthodes des déplacements et des forces,
l’année suivante TURNER et CLOUGH publient une présentation systématique de la méthode
des déplacements. Ces deux publications sont particulièrement importantes et représentent
[26]
véritablement le début de la MEF comme technique de calcul des structures complexes .
D’ailleurs, il est vital de savoir que la résolution d’un problème par la MEF conduit à résoudre
un système d’équation aux dérivées partielles (E.D.P). Heureusement, l’automatisation des
calculs est concrétisée par l’avènement du premier ordinateur vers les années 50 et cela a facilité
les taches qui sont vraiment laborieuses si on les exécute manuellement.
III.2 Quelques rappels théoriques
III.2.1 Principe de la modélisation
La modélisation est une méthode qui consiste à créer un modèle physique qui reflète
le fait réel afin d’obtenir un modèle informatique exploitable.
Les étapes à suivre pour la modélisation d’un problème physique sont les suivantes :
 la modélisation physique ;
 la modélisation mathématique ;
 la modélisation numérique ;
38
 la modélisation informatique.
III.2.2 La modélisation physique
L’objectif de cette étape consiste à la reconstruction simplifiée de l’état initial du

système à étudier, en tenant compte de ses propriétés (physique et géométrique) ainsi que les
contraintes qui lui sont appliquées. Les critères suivants sont dument satisfaites pour que le
modèle physique soit précis :
 la géométrie du modèle doit vraiment imiter la forme du système à étudier, sans oublier
ses dimensions ( en tenant compte des longueurs, largeurs, la courbure,..) ;
 les propriétés physiques initiales du système à étudier doit être en fonction des lois de
comportement des matériaux (module d’élasticité E d’Young, coefficient ν de Poisson,
masse volumique ρ ) ;
 les conditions aux limites (contrainte initiale, déplacement initial,...) doivent être
connues ;
 la nature et l’emplacement des sollicitations extérieures doivent être précises.
Il existe deux méthodes de modélisation physiques possibles pour l’étude et la
prévision du comportement d’un ouvrage souterrain : les essais expérimentaux et les modèles
réduits. La deuxième approche permet de multiplier les essais qui sont plus rapides,
reproductibles, moins chères et qui peuvent aboutir à la rupture.
En effet, une campagne expérimentale a été réalisée par Mair (1979) sur un modèle réduit en
centrifugeuse. L’objet de cette expérience est d’évaluer la stabilité du front de taille et les
déplacements induits par une excavation. D’ailleurs, Chambon (1990) et Skiker (1995) ont
mené des essais sur des modèles réduits dans le Laboratoire Central des Ponts et Chaussées
(LCPC). Chambon a conclu que dans un sol pulvérulent et sec, le diamètre de la galerie est un
paramètre qui détermine la stabilité de l’ouvrage et que plus le matériau excavé est dense, plus
le front de taille est stable [11].
III.2.3 La modélisation mathématique
C’est une étape qui consiste à mettre en équation le problème à étudier après avoir
créé son modèle physique. En effet, elle assure la liaison entre les inconnues et les variables
connues à travers les hypothèses. Dans la plupart des cas, ces équations sont écrites sous-
forme des équations aux dérivées partielles (E.D.P). Cette étape est basée sur des théorèmes
39
et des hypothèses qui sont utilisées surtout dans le domaine de mécaniques des milieux
continus (M.M.C).
Ces dernières peuvent être des équations d’équilibre du système ou des équations qui
introduisent les lois de comportement des matériaux ou des équations qui lient les
déformations et les conditions aux limites.
III.2.4 La modélisation numérique
La résolution de ces équations aux dérivées partielles est un travail laborieux et

demande beaucoup de temps. Alors, pour pallier à cette difficulté, il faut transformer ces E.D.P
par des équations matricielles (transformer le problème mathématique continue en un
problème mathématique discret) dont sa dimension est finie. C’est alors une étape de
discrétisation du problème.
III.2.5 La modélisation informatique
C’est une étape qui consiste à écrire des algorithmes afin de créer un programme
permettant de résoudre rapidement et avec précision le problème discrétisé. Les méthodes
numériques sont efficaces pour la résolution des problèmes de calcul qui demandent beaucoup
de temps ainsi que d’effort. D’ailleurs, elles peuvent être utilisées dans le cas où le milieu
serait continu et isotrope et dans le milieu discret.
III.2.6 Milieu continu et isotrope
Un milieu est dit isotrope si ses propriétés physiques restent invariantes dans toutes
les directions. Quant aux milieux continus, ils peuvent être utilisés afin de simplifier la
modélisation numérique de milieux complexes comme les massifs rocheux. Il est à noter que ce
modèle exige que le matériau reflète à la fois les propriétés de la roche intacte, du litage et des
discontinuités [23]
. L’utilisation de ce modèle est généralement recommandée aux massifs
homogènes ou comportant peu de fracture. Il existe deux méthodes de calcul qui se sont basées
sur la discrétisation spatiale du massif rocheux [15] :
 la méthode des différences finies : qui est une méthode simple en termes d’écriture
et présente un faible coût de calcul mais elle se limite sur la prise en compte des formes
géométriques simples. Elle consiste à remplacer les dérivées partielles par des
40
différences divisées (combinaison de valeurs ponctuelles de la fonction) en un nombre

fini de point discret.
 la méthode des éléments finis : qui consiste à approcher dans une sous-espace
quelconque de dimension finie un problème variationnel posé dans un espace de
dimension infinie[19]. Elle permet de traiter des géométries très complexes.
III.2.7 Milieu discret
Les milieux continus et isotropes ne sont pas vraiment une généralité en travaux
souterrain et mécanique des roches. En effet, les rochers sont constitués des discontinuités
(fracture, fente, …) et surtout composés des divers éléments comme les minerais, des liquides,
des gaz, etc. Ils sont discontinus et hétérogènes donc anisotropes. C’est la raison pour laquelle
les mécaniciens des roches ont élaboré la méthode de calcul adaptée au milieu discret appelée
méthode des éléments distincts (M.E.D).
Cette dernière est utilisée en fonction du degré de fracturation des roches. A travers
la M.E.D, les milieux discrets, allant de milieux faiblement fracturés (en utilisant la méthode
des dièdres) à des milieux fortement fracturés, peuvent être représentés. D’ailleurs, elle permet
de prendre en compte les grands déplacements et les rotations.
III.3 La discrétisation
La discrétisation consiste à diviser le domaine d’étude en sous domaine de petits
éléments. Chaque élément est défini géométriquement par un nombre de nœuds bien
déterminés qui constituent en général ses sommets. Cette opération est connue également sous
le nom du maillage. Les nœuds sont les points particuliers qui caractérisent la forme
géométrique de l’élément.
La figure ci-dessous montre une discrétisation d’un domaine en éléments
triangulaires.
41
Figure 39 Figure 36 Discrétisation du domaine D en éléments triangulaires [01]
Soit D le domaine à étudier, la création d’un nombre fini d’éléments notés 𝒏𝒊

représenté par le sous domaine 𝑫𝒊 explique le principe général de la discrétisation.
𝒏
𝒊
𝐷 = ∑𝑖=1 𝑫𝒊 (18)
Les règles à respecter de la discrétisation d’un domaine sont les suivantes :
Table 1Les règles à respecter de la discrétisation
Un nœud d'un élément ne doit pas

être intérieur à un côté d'un autre du même
type.
Aucun élément bidimensionnel ne

doit être plat, éviter les angles proches de
180°
ou de 0°.
Deux éléments distincts ne peuvent
avoir en commun que des points situés dans
leurs frontières communes ; le recouvrement
est exclu.
42
L'ensemble de tout éléments doit

constituer un domaine aussi proche que
possible du domaine donné ; les trous entre
éléments sont exclus.
Parmi les méthodes de discrétisation susmentionnées, c’est la M.E.F qui présente un

grand avantage à cause de sa souplesse de prise en compte des conditions aux limites et son
adaptation aux diverses formes géométriques.
III.4 Les différentes phases de calcul

Les phases de calcul pour la modélisation par la MEF sont les suivantes :
 la phase de prétraitement ;
 la phase de calcul ;
 la phase de post-traitement.
III.4.1 La phase de prétraitement
Cette phase correspond à la création du modèle physique et du modèle mathématique

correspondant au problème. D’ailleurs, pendant cette phase, on commence par la discrétisation
du géométrique de l’élément. Le principe est effectivement basé sur le choix du type d’élément
à générer et la création du maillage.
III.4.1.1 Le maillage général
Par définition, le maillage est un assemblage de ponts isolés appelés nœuds situés dans
[23]
un domaine de définition des fonctions assujetties aux équations aux dérivées partielles
(E.D.P). Sa création relève des bases de calcul matriciel assez important et est caractérisée par
le dressage des deux matrices suivantes :
 la matrice de coordonnées, contenant les coordonnées de chaque nœud dans
le repère utilisé et dont la taille est égale au produit du nombre de nœud et la
dimension de l’espace d’étude. Celle-ci caractérise la forme de l’objet à mailler
parce qu’elle consiste à créer un grand nombre de point donnant cette forme.
43
 la matrice de connectivité, qui caractérise la connexion des nœuds à chaque

élément. En effet, il s’agit de créer une liaison entre les points qui donnent la
forme de l’élément.
Il est vital de préciser que la création d’une forme volumique (élément 3D) est
générée par les trois variables de x, y et z. Voici un exemple montrant la détermination de
la matrice de connectivité et la matrice de rigidité.
Le domaine D qu’on va discrétiser est divisé en n éléments (segment) de taille 1/n.
Chaque élément contient deux nœuds sur lesquels la fonction u est interpolée.
Figure 40 Discrétisation du domaine D
Le maillage est la division du domaine D en plusieurs éléments. Sa matrice de connectivité

𝐶𝑐𝑛 et matrice de coordonnées 𝐶𝑐𝑑 respectives sont les suivantes :
0,0
1 2
2 3 0,5
𝐶𝑐𝑛 =( ) et 𝐶𝑐𝑑 = 1
3 4
1,5
4 5
(2)
III.4.1.2 Le maillage d’élément fini pour les ouvrages souterrains
Les éléments finis de massif bidimensionnels (Tableau 7) sont constitués par des
éléments de forme triangulaire et quadrangulaire [18] dont les caractéristiques sont données par
le tableau ci-dessous :
Tableau 7 Les caractéristiques des éléments massifs bidimensionnels
Type d’élément Fonction Avantage Limite

d’interpolation
Triangle à Linéaire  En déformation  Résultats non
trois nœuds plane : satisfaisants pour
déformation et les ouvrages
contrainte sont complexes (avec
prédominance de
44
constantes sur flexion, lois de

l’élément ; comportement
 Les matrices de non linéaire)
rigidités les
forces
élémentaires
peuvent être
exprimés
analytiquement
Quadrilatère Forme bilinéaire  En déformation plane :  Un calcul qui
à quatre nœuds de degré 1 les déformations et les n’est pas
contraintes ont suffisamment
une variation linéaire rigoureux de la
sur l'élément ; flexion ;
 Performant par rapport  Résultats non
au triangle à trois satisfaisants pour
nœuds ; les problèmes
 Le maillage peut être réels (faisant
raffiné sans accroitre intervenir à des
la largeur de la bande géométries très
dans le cas des complexes)
problèmes en
comportement linéaire
et homogène.
Triangle à Forme  En déformation  Résultats pas
six nœuds quadratique plane, les contraintes vraiment
complète et les déformations satisfaisants
varient linéairement
sur l'élément ;
 Meilleure description
des variations des
contraintes ;
45
 En autorisant la
déformation des
cotés initialement
droits ou courbes, il y
a possibilité de la
description de
contours
géométriques courbes
(tunnel) grâce aux
nœuds milieux de
chaque élément ;
Quadrilatère à huit Forme  Très performant pour
nœuds quadratique le calcul des structures
incomplète des massifs
bidimensionnels ;
 Description
satisfaisante des
champs de
déformations et de
contraintes
La figure suivante nous montre les éléments bidimensionnels que constituent les
éléments finis du massif.
46
Figure 41 Eléments de massif bidimensionnel [37]
A part les éléments bidimensionnels susmentionnés, il existe aussi des éléments de

massifs tridimensionnels dont les performances relatives sont assez proches de leur
homologue bidimensionnel. L’élément hexaédrique à vingt nœuds est l’élément le plus utilisé
actuellement dans le maillage tridimensionnel. Les règles générales pour construire un
maillage d’éléments finis sont les suivantes :
 un maillage d'éléments finis de massif doit tenir compte de l'essentiel de la géométrie
de la structure (courbure et cassure) ;
 le maillage doit être suffisamment fin dans les zones où les plus fortes variations des
déplacements et des contraintes sont attendues ;
 le rapport de la plus grande dimension de l'élément sur sa plus petite dimension
(rapport déforme ou aspect ratio, en anglais) soit proche de l'unité.
Pour plus de précision, on peut créer un raffinage, qui nécessite la manipulation des
deux matrices caractéristiques, dans les zones présentant un changement brusque de géométrie
en prenant en compte :
 d’un maillage de type structuré ( les éléments sont identiques et de même taille) ;
 d’un maillage de type non structuré, qui consiste à augmenter le nombre de nœud
dans la zone à raffiner et déplacer les nœuds déjà en place dans la zone à raffiner.
47
Les éléments de massifs tridimensionnels les plus utilisés sont montrés par la
figure 39 suivante :
Figure 42 Les éléments de massif tridimensionnels les plus utilisés [18]
La figure suivante montre un exemple de maillage 2D et 3D :
Figure 43 Exemple de maillage en 2D et en 3D [12]
III.4.2 La phase de calcul
C’est la phase pendant laquelle on effectue les opérations suivantes :

 l’élaboration de la matrice de rigidité élémentaire ;
 l’assemblage des matrices de connectivité et la matrice de rigidité élémentaire afin
d’obtenir la matrice de rigidité globale [K] ;
 la prise en compte des conditions aux limites ;
 la résolution du système d’équation linéaire.
48
III.4.2.1 Création de la matrice de rigidité élémentaire
C’est une matrice qui relie le déplacement des nœuds et les sollicitations nodales. On
l’obtient à partir de la discrétisation du travail intérieur due à une déformation virtuelle. Cette
dernière est de la forme :
[𝒌𝒆 ] {𝒖𝒆 } = {𝒇𝒆 } (19)

Où
{fe } désigne le vecteur élémentaire des forces appliquées aux noeuds;
{ue } est le vecteur élémentaire des déplacements nodaux;
[𝒌𝒆 ] est la matrice de rigidité élémentaire, qui dépend de la fonction d’interpolation choisie et
des propriétés élastiques du matériau. Sa dimension est égale au produit de l’espace d’étude
et le nombre de nœuds de l’élément.
III.4.2.2 Assemblage
C’est une étape consistant à associer les matrices [𝒌𝒆 ] et les vecteurs élémentaires
{fe } et {ue } pour former un système global {F} = [K] {U} de manière à satisfaire les
conditions d’équilibres aux nœuds.
Les règles d’assemblages sont définies par les relations suivantes :
[K] = ∑𝒆[𝒌𝒆 ] et {F} = ∑𝒆[𝒇𝒆 ] (20)
e désigne le nombre du degré de liberté ;

[K] est matrice de rigidité, qui représente la loi de comportement de l’élément fini ;
{F} est le vecteur contenant les vecteurs représentant les actions extérieures sur
l’élément fini.
III.4.2.3 Les conditions aux limites
Les conditions aux limites ont pour rôle de modifier le système global. En effet, sans
la prise en compte de ces dernières, la matrice de rigidité [K] est encore singulière, ce qui
signifie que l’équilibre n’est pas encore atteint. Il est vital de noter que les conditions aux
limites se présentent sous différentes formes :
49
 une forme régulière : il existe au moins des contacts parfaits sur un morceau de la
frontière avec un autre solide indéformable :
o premier cas : si on considère 𝜕Ω0 un morceau de frontière du système. Alors, à
chaque fois qu’une composante 𝐹𝑘 des efforts extérieurs surfaciques est connue
sur 𝜕Ω0 , on peut écrire les conditions suivantes :
𝝈𝒌𝒋 𝒏𝒋 = 𝑭𝒌 𝒔𝒖𝒓 𝝏Ω𝟎 (21)

avec
𝜎𝑘𝑗 désigne une composante de tenseur de contrainte ;

𝑛𝑗 désigne le normal extérieur de la surface 𝜕Ω0 .
C’est le type de condition aux limites connus sous le nom de condition aux limites
en contrainte.
o deuxième cas : le contact rigide annonce la connaissance des déplacements sur le
morceau de 𝝏Ω𝟎 . Par conséquent, l’équilibre impose des valeurs particulières
aux translations ou rotation sur 𝝏Ω𝟎 .C’est ce qu’on appelle condition aux limites
en déplacement ou condition aux limites essentielles.
 une forme imprécise ou incomplète (les efforts sont appliqués d’une manière
imprécise) ;
 une forme unilatérale, qui conduit à des conditions de contact exprimées par des
inégalités conditionnelles. (Comme le cas des frottements indiquées par la loi de
Coulomb).
III.4.2.4 Résolution du système linéaire
La mise en équation matricielle conduit à la résolution d’un système d’équation

linéaire de la forme suivante.
𝑭 = 𝑲𝑼 (22)
.
Si l’élément fini a e degré de liberté, on peut écrire :
𝑭𝟏 𝑲𝟏𝟏 ⋯ 𝑲𝟏𝒆 𝑼𝟏
( ⋮ )= ( ⋮ ⋱ ⋮ ).( ⋮ )
𝑭𝒆 𝑲𝒆 … 𝑲𝒆𝒆 𝑼𝒆
50
Pour des problèmes en 1D, le tableau ci-dessous montre quelques exemples de

signification de ces matrices :
Tableau 8 Exemples illustrant les significations de chaque matrice [38]
Problème F K U
Ressort Force Raideur k déplacement
Barre Force Rigidité EA/L déplacement
Thermique Flux de chaleur Conductivité αA/L température
Electricité Intensité du courant Résistivité 1/R La tension du
courant
Ecoulement Débit Perméabilité Charge hydraulique
III.4.3 La phase de post-traitement
Cette phase correspond à la présentation des résultats sous forme numérique et sous
forme graphique. D’ailleurs, on effectue les fonctions complémentaires comme l’interpolation
et l’interprétation des résultats obtenus par le calcul.
Le principe de calcul par la méthode des éléments finis se résume par la figure ci-
dessus.
Figure 44 Les trois étapes du calcul : prétraitement, résolution et post-

traitement [13]
51
CONCLUSION
La Méthode des Eléments Finis est une branche des mathématiques appliqués
permettant de résoudre des problèmes en ingénierie à travers des équations aux dérivées
partielles (EDP). Des logiciels qui se basent sur le principe de cette méthode ont été
développés pour faciliter les calculs afin de résoudre des problèmes 2D ou 3D. La MEF est
une méthode de discrétisation du problème à résoudre. Pour aboutir aux résultats attendus,
elle possède trois (03) grandes phases de calcul à suivre, qui demandent beaucoup de temps
si on les fait manuellement. Heureusement, des informaticiens ont déjà développé des logiciels
qui se basent sur les concepts de la méthode des éléments finis pour faciliter le calcul. C’est
l’un de ces logiciels que nous allons utiliser pour la modélisation numérique du creusement
du tunnel dans le chapitre 5.
52
CHAPITRE IV METHODOLOGIE
Chapitre IV. METHODOLOGIE
IV.1 Introduction
Un des problèmes majeurs que la population de la capitale de Madagascar doit faire
face est l’embouteillage. Ce fléau est dû non seulement à l’étroitesse des routes mais aussi à
ses états actuels qu’elles ont vraiment besoin d’entretiens. La solution que nous proposons dans
cette étude est le creusement d’un tunnel. La construction d’un tel ouvrage nécessite des études
préliminaires se basant sur une méthodologie adéquate à savoir :
 la localisation de la zone d’influence ;
 la collection des données à travers des études bibliographique ;
 le tracé du tunnel ;
 la modélisation numérique qui consiste à :
o créer un modèle conceptuel à travers les données obtenues ;
o la justifier les paramètres d’entrée dans le logiciel Plaxis 2D ;
o afficher les résultats obtenus ;
o comparer les résultats obtenus par calcul direct (méthode analytique)
IV.2 Les tunnels routiers à Antananarivo

Avant les années 1910, pour passer très vite d’Analankely à Mahamasina ou à
Ambanidia, les populations d’Antananarivo devaient escalader les collines.
Figure 46 Tunnel Hubert Grabit en Figure 46 Tunnel d’Ambanidia

1930
Connaissant les difficultés causées par cela, l’administration coloniale décidait de
(Source : Z.R 2013)
creuser le premier tunnel routier Hubert Garbit (actuellement tunnel Ralaimongo), reliant
(
53
Analakely et Anosy, en 1914. Ce tunnel creusé en forme de fer à cheval à trois (03) voies a
116m de long et 11m de diamètre. Ensuite, un nouveau tunnel connu sous le nom Tunnel Cayla
(actuellement Tunnel d’Ambanidia) fut creusé sous le col d’Ambohijatovo. Ce dernier a les
mêmes caractéristiques que le précédent à part sa longueur qui mesure 130m.
IV.3 Localisation administrative de la zone d’étude

Le nouveau tunnel routier que nous allons modéliser dans cet ouvrage se trouve dans
le premier arrondissement de la Commune Urbaine Antananarivo, District d’Antananarivo
Renivohitra, Région Analamanga. L’entrée de ce dernier se trouve à Ampantsakàna, dans le
Fokontany Isoraka et la sortie se situe à Analakely, Fokontany Antanimalalaka.
IV.4 Etudes bibliographiques
IV.4.1 La géologie de la zone d’étude
IV.4.1.1 Géologie régionale
Le travail de synthèse cartographique de Madagascar fait par Besairie et Hottin,

respectivement en 1968 et en 1976, a montré que le socle malagasy contient des granites
stratoïdes d’épaisseur hectométrique. Ces derniers caractérisent tous les hauts plateaux
malagasy. D’ailleurs, ce sont des granites de type A, et de couleur grise ou rose. En effet,
l’étude structurale du secteur situé à Ankazobe région Analamanga a montré que les granites
[36]
stratoïdes et leur encaissant (migmatite et gneiss) sont coostructurés . Un faible pendage
vers l’ouest a été observé. (Voir figure ci-contre)
Figure 47 Série stratoïde dans la partie Nord d’Antananarivo (d’après Nédélec

et Al.,1994) [36]
A part cela, d’après des recherches faites par ARISOA Rivah Kathy dans son ouvrage
qui s’intitule « Définition des paramètres géologiques en vue de la modélisation géophysique
54
du socle cristallin à la latitude d’Antananarivo », le domaine d’Antananarivo est le domaine le

plus complexe de Madagascar.
Figure 48 Schéma géologique d’Antananarivo[02] (numérisé et modifié par

TOJONIRINA Jean Claudio)
Elle a affirmé dans son ouvrage qu’en 1996 Hottin a proposé qu’en ce qui concerne
la géochronologie du socle cristallin, Antananarivo se situe dans les terrains archéens (2,8 à
3,0 milliards d’années). Ensuite il est, d’après les nomenclatures de Bésairie (1951, 1965,
1973), constitué par les formations du groupe d’Ambatolampy du système du graphite [2]
. Le
socle cristallin d’Antananarivo, représenté par la carte établie par ARISOA Rivah Kathy ci-
contre, est constitué par des paragneiss associés à l’orogenèse Shamwaïenne (2,6 Ga).
D’ailleurs, la granitisation récente a individualisé dans la région d’Antananarivo des unités de
granitoïde (granite migmatitique, migmatite granitoïde).
En revanche, après avoir superposé la carte originaire avec l’image satellitaire
capturée par Google Earth, nous avons constaté que celle-ci n’est pas précise en termes
d’échelle. C’est pour cette raison que nous avons décidé de créer une nouvelle carte
55
géologique, en partant de celle qui a été numérisée, ayant une échelle bien précise et localisant
les endroits exacts où l’on peut trouver ces différentes formations géologiques.
Figure 49 Carte géologique améliorée
(Numérisée par TOJONIRINA Jean Claudio - Fond : Image Google Earth)

Ensuite, madame ARISOA Rivah Kathy a pu tracer des coupes géologiques montrant
les formations géologiques dominantes, qui sont les alluvions, gneiss à graphite, gneiss à
pyroxène, gneiss et migmatite ainsi que des gneiss migmatitiques.
56
Figure 50 Coupe géologique suivant B-B’ [02]

En outre cela, nous avons utilisé une carte de structure tectonique d’Antananarivo déjà tracée
par l’auteur susmentionné afin d’avoir une autre (figure ci-contre).
Figure 51 Carte de structure tectonique d’Antananarivo
(Numérisée par TOJONIRINA Jean Claudio - Fond : Image Google Earth)
IV.4.1.2 Géologie locale
D’après ARISOA Rivah Kathy, les formations qui sont en majeur partie dans notre
zone d’étude (Ampantsakàna et Analakely) sont les gneiss et les migmatites. En effet, elles
correspondent, aux anciennes formations sédimentaires métamorphisées lors de l’orogenèse
majeure Shamwaïenne de 2600 M.a et sont migmatitisées [02].
En prenant les études faites par Besairie et Hottin, notre zone d’étude est formée par
de gneiss et migmatite dans lesquelles s’enfouissent le granite de type A. C’est la même chose
pour les études faites par ARISOA Rivah Kathy mais ce qui différencie son rapport d’études
avec celui de Besairie est le type de granite parce qu’ici, elle parle du granite migmatitique.
Alors, en analysant les deux cartes géologiques de la région Analamanga élaborées par ces
auteurs, on peut constater l’alternance du granite avec du gneiss et migmatique d’Est en Ouest.
57
Il est à noter que le creusement d’un tunnel nécessite la connaissance de nature de

chaque formation géologique de la zone d’étude et surtout des paramètres géotechniques
obtenus à travers des prélèvements d’échantillons sur terrain et des analyses de celles-ci en
laboratoire. Donc, nous allons prendre en compte le fait que les granites migmatitiques
proposés par Rivah Kathy soient remplacés par des granites stratoïdes encaissés par du gneiss
et de migmatitique (proposé par Besairie et Hottin). En effet, Paquette et Nedelec ont affirmé
en 1998 que ces granites stratoïdes affleurent dans le domaine d’Antananarivo sous forme de
corps planaires de 1m à 1km d’épaisseur et pouvant étendre jusqu’à plusieurs kilomètres.
Ensuite, ils ont affirmé également que la morphologie de ces derniers est en sills concordants
dans les gneiss et les schistes encaissants. D’ailleurs, ils peuvent se présenter des veines de
charnockites [34]. De ce fait, nous allons déduire que dans notre zone d’étude pourraient se
trouver les formations géologiques suivantes :
 du sédiment détritique comme le silt ou le pélite;
 de l’argile péliteuse ou silteuse ;
 du granite stratoïde de type A (de couleur grise) : selon Besairie et Hottin ;
 du migmatite et de gneiss ;
 du charnockite, se présentant en veine ;
 des roches à structures gneissiques.
IV.4.2 Résultats des études faites par le LNTPB [39]
Les résultats que nous allons présenter dans cette partie est le fruit des études
bibliographiques.
IV.4.2.1 Résultats géophysiques
Après avoir effectué une prospection électrique, le LNTPB a pu élaborer deux profils
à travers lesquels il peut conclure que les formations suivantes sont trouvées sur le site :
 couche conductrice de nature gneissique, ayant comme puissance égale à 8 m et
une résistivité inférieure à 150Ωm;
 couche intermédiaire composée des roches dures de nature gneisso-
orthogneissique dont l’épaisseur est égale à 12m et de résistivité comprise entre
150Ωm et 800Ωm;
58
 couche résistante composée des roches dures de nature orthogneissique (granites

migmatitiques, migmatite granitique) localisé entre 12m et 20m avec une
résistivité atteignant jusqu’à 1500Ωm.
Les figures suivantes nous présentent les deux profils géoélectriques sus-mentionnés.
Figure 52 Profil géoélectrique du profil II
(Source : RAHAMBININTSOA Sedra Avojanahary, 2018)
59
Figure 53 Profil géoélectrique du profil I
IV.4.2.2 Résultats de forages
La localisation de 16 points de forages faits par le LNTPB en 2017 et l’Entreprise SIF Bachy daté de 1962, 1963 et 1969 est présentée par
la figure ci-contre et ces résultats sont présentés dans le tableau ci-contre..
60
Table 2 Les points de forages
Points X(m) Y(m) Altitude (m) Profondeur (m) Mois/année

S1 514 088,382 798 548,459 1262 26.42 Aug-17
S2 514 140,505 798 563,276 1274 40.7 Aug-17
S3 514 311,162 798 599,656 1286 40.3 Aug-17
S4 514 391,86 798 616,193 1288 39.55 Aug-17
S5 514 508,277 798 640,534 1292 40.03 Aug-17
S6 514 569,793 798 711,31 1289 39.7 Aug-17
S7 514 620,063 798 805,238 1280 29.1 Aug-17
S8 514 057,823 798 604,286 1266 20 Sep-62
S9 514 080,974 798 617,515 1270 30 Sep-62
S10 514 137,198 798 657,203 1278 15 Jan-69
S11 514 075,021 798 634,713 1276 60 Feb-69
S12 514 028,057 798 591,057 1261 60 Feb-69
S13 514 549,287 798 567,112 1296 14.3 Jan-63
S14 514 575,746 798 622,013 1294 15.75 Jan-63
S15 514 523,490 798 584,31 1295 17.4 Jan-63
S16 514 549,949 798 631,274 1293 15.15 Jan-63
En faisant l’interpolation des points de forages nous avons, d’après

RAHAMBININTSOA Sedra (2018), le modèle lihologique suivant le tracé du tunnel est
présentée sur la figure suivante.
61
Figure 54 Modèle lithologique suivant le tracé du tunnel (Source : RAHAMBININTSOA Sedra, 2018)
62
En faisant une coupe longitudinale du modèle ci-dessus, l’auteur susmentionné a

affirmé que le tracé du tunnel traverse trois (03) types de sections différentes telles qu’une
section en sol, une section variable et une section entièrement en roche. En effet, les formations
pouvant être rencontrées lors du creusement du tunnel sont le limonx argileux micacée, le
limon sableux micacée, la migmatite, le charnockite et le granite migmatitique (Source :
RAHAMBININTSOA Sedra,2018).
Figure 55 Coupe géologique longitudinale prévisionnelle du tracé
(Source : RAHAMBININTSOA Sedra, 2018)
IV.4.2.3 Résultats Hydrogéologiques
L’étude hydrogéologique est à ne pas négliger dans l’étude de creusement d’un tunnel
parce qu’elle nous permet de connaitre l’existence de l’eau ainsi que sa circulation dans les
massifs, ce qui conditionne le choix du mode de creusement du tunnel. En effet, l’eau peut
avoir une grande influence sur la résistance de la roche
D’après RAKOTOVAHOAKA E. (2008) dans son ouvrage qui s’intitule « Cartographie de la
vulnérabilité des nappes d’eaux souterraines de la plaine d’Antananarivo et de ses environs »,
trois types de nappe existent sous la formation cristalline d’Antananarivo, à savoir :
 la nappe de socle fissuré ;
 la nappe de la complexe arène micacé ;
 la nappe d’altérité.
En outre, l’hydrogéologie d’Antananarivo est caractérisée par :

 une perméabilité K ayant des valeurs comprises entre 10-6 et 10-5 m/s ;
 un coefficient d’emmagasinement variant de 0,1 à 0,001% pour les nappes captives et
1 - 20% pour les nappes libres [38] .
63
La nappe rencontrée sur le site est surtout celle de socle fissuré. En effet, l’eau souterraine
circule dans un réseau de fracture pouvant être interconnecté avec un débit spécifique de l’ordre
de 0,5l/s/m.
A part cela, des essais piézométriques (présenté par le tableau suivant) ont été faits par
l’Entreprise SIF Bachy (1962, 1963, 1969) et LNTPB (2017) dans différentes saisons [39] :
 la saison d’hiver (frais et sec) : Mai à Octobre ;
 la saison d’été (doux et pluvieux) : novembre à avril.
A travers ces essais, les formations hydrogéologiques trouvées sur le site sont les dépôts
meubles (composées de limon argileux et sable limoneux) et des roches comme des granites.
Il est alors vital de mentionner qu’à travers ces essaies, le tunnel à creuser est situé sous le
niveau de la nappe.
IV.4.2.4 Résultats géotechniques
D’après les essais en laboratoires des échantillons prélevés au niveau des forages
susmentionnés, nous avons les résultats suivants [39] :
IV.4.2.4.1 Propriétés des sols
Les résultats sur les propriétés du sol sont résumés dans le tableau suivant.
Table 3 Les propriétés du sol
Classification Teneur en eau Indice de Matière

moyenne W [%] plasticité Ip [%] Organique [%]
LPC GTR
Argile peu plastique A2 20,65 13,82 1,4
Argile très plastique A4 27,35 18,23 0,2
Limon peu plastique A2 27,20 12,77 0,5
Limon très plastique A3 58,30 26,70 0,16
Sable argileux A2 15,15 9,33 0,4
LPC : Laboratoire des Ponts et Chaussées (voir annexe 4a)

GTR : Guide de Terrassement Routier (voir annexe 4b)
64
Selon la LNTPB (2017), la teneur en eau du sol varie de 15,15 % à 58,30 % et il est
peu plastique vu que son indice de plasticité est entre 9,3% et 26,7%. D’ailleurs, l’indice de
plasticité du sol affecté par le tracé 2 (voir figure 56) est 4.6% à 16.5% (peu plastique à
plastique) [39]. Quant aux matières organiques, la LNTPB (2017) affirme que leur teneur est
comprise entre 0,1% à 1,4%, qui n’ont aucune influence sur le tunnel à creuser [39].
IV.4.2.4.2 Propriétés des roches
Les résultats concernant les propriétés des roches sont résumés dans le tableau
suivant :
Table 4 Résultats sur les propriétés des roches sur le site
Paramètre Valeur Interprétation

Porosité 2,3% à 3,53% Porosité faible
Résistance à la traction 4 MPa à 6,9 MPa Les roches sont résistantes par
indirecte (Rt) rapport à la traction (voir annexe 3a)
Résistance à la 41 MPa à 68 MPa Roche résistante (voir annexe 1)
compression (Rc)
Classification RMR 62 à 77 Massif rocheux en bonne qualité
(voir annexe 3b)
Espacement des 0, 09 m à 1,50 m Espacements modérés en majeure
discontinuités partie (voir annexe 4c)
IV.4.3 Sismicité
Antananarivo et ses environs connaissent une activité séismique peu importante et

dont la magnitude dépasse très rarement 5 dans l’échelle de Richter. Cela peut être expliqué
par le fait que la lithosphère de la partie centrale de Madagascar est disloquée en des blocs en
quête de leur position d’équilibre [02].
IV.5 Tracé du tunnel

Les figures ci-après nous montrent les tracés du tunnel à creuser. Le choix de
l’orientation de ces tracés est expliqué par le fait de faciliter le trajet pour aller d’Ampantsakàna
à Analakely et surtout les structures géologiques. D’ailleurs, on sait très bien que la réalisation
65
d’un tel projet nécessite des expropriations. Ainsi, d’après des analyses, on a abouti au choix
de ces deux axes suivants, pour éviter la démolition d’un grand nombre d’infrastructures
existantes :
 tracé 1 : c’est un axe liant Ampantsakana, sur le Lalana Rakotomahefa (Latitude :
18°54'39.67"S et Longitude : 47°31'11.71"E), et Antsahavola, au Lalana Rabearivelo
(Latitude : 18°54'33.34"S et Longitude : 47°31'20.48"E) pour satisfaire aux critères
susmentionnées. La longueur de ce tracé est à peu près 300 m.
 tracé 2 : c’est un axe liant Ampantsakana et Analakely (Latitude :18°54'30.24"S et
Longitude :47°31'37.01"S). Sa longueur est approximativement égale à 700 m.
Figure 56 Les deux tracés du tunnel (Fond : Google Earth)
Les figures suivantes nous présentent les coupes longitudinales de ces deux tracés.
Figure 57 Tracé 1 (Axe Ampantsakàna-Antsahavola)
66
Figure 58 Tracé 2 (Axe Ampantsakàna - Totoha-baton'Antaninarenina)
IV.6 Description du tunnel

En tenant compte de la géologie, de l’existence des infrastructures (bâtiments) ainsi
que son utilisation, on préconise les caractéristiques (dimensions, méthode de creusement, etc.)
suivantes du tunnel à creuser.
IV.6.1 Les dimensions du tunnel
Le tunnel à creuser est un tunnel routier bidirectionnel et nous allons choisir la forme
de section en fer à cheval (tunnel vouté).
Figure 59 Section visible en forme de fer à cheval du tunnel routier [05]

Le tableau suivant nous résume les dimensions du tunnel routier bidirectionnel à creuser :
67
Tableau 9 Les dimensions du tunnel
Section en forme de fer à cheval

Diamètre de l’intrados 11,2
Hauteur [m] 8
Aire de surface de la section visible du 75,5
tunnel [m2]
Profondeur de l’axe du tunnel [m] Atteint jusqu’à 30 m
Pente moyenne [°] 3°
IV.6.2 Méthodes de creusement
Compte tenu de l’existence des maisons et d’autres infrastructures sur et aux alentours
de l’endroit où l’on va creuser le tunnel, et compte tenu de la qualité des roches rencontrées,
l’utilisation de l’explosif comme moyen d’abattage n’est pas recommandé parce qu’ils
entrainent des vibrations qui peuvent endommager ces infrastructures sus-jacentes. D’ailleurs,
puisque le massif est hétérogène et que la forme du tunnel est en fer à cheval alors, on
préconise, l’utilisation d’une machine à attaque ponctuelle combiné avec la technique de
creusement en demi-section.
CONCLUSION
On a parlé, dans ce chapitre, des points importants concernant la zone d’étude, qui est
située dans la région Analamanga. Plus précisément, la zone d’influence où l’on ferait le
creusement est entre Analakely et Ampantsakàna. D’après des études bibliographiques nous
avons pu connaitre les formations géologiques pouvant être rencontrées lors du creusement du
tunnel à travers des forages et des études géotechniques. D’ailleurs, d’après des essais
piézométrique faits par le LNTPB et l’Entreprise SIF Bachy le tunnel à creuser est situé au-
dessous de la nappe phréatique. Ces résultats vont être utilisés pour que la modélisation
numérique traitée dans la partie suivante soit faite.
68
CHAPITRE V MODELISATION NUMERIQUE
CHAPITRE V
MODELISATION NUMERIQUE
Chapitre V. MODELISATION NUMERIQUE
Dans ce chapitre, on va se focaliser sur la présentation du logiciel Plaxis 8.5 version

d’éssai, dans un premier temps puis son utilisation pour la modélisation du creusement du
tunnel routier liant Analakely et Ampantsakàna.
V.1 Historique du Plaxis

En 1987, les plaines de la Hollande ont subi des fréquents remontées d’eau. Pour faire
face à cela, les ingénieurs géotechniciens Hollandais ont eu l’idée de développer un programme
donnant la possibilité de réaliser des analyses concrètes sur le comportement du sol, surtout au
niveau des digues, talus et des berges.
Afin de trouver des solutions pratiques et efficaces, l’Université DELFT de
technologie du Pays Bas développa un logiciel qui se base sur la méthode des éléments finis,
nommé PLAXIS. La première version en 2D de ce dernier exploitable sous Windows a apparu
en 1998. Depuis, le logiciel est utilisé par nombreux ingénieurs géotechniciens dans le monde
entier et a connu ses succès. Ce succès encourage les développeurs d’améliorer le logiciel
jusqu’à ce qu’ils obtiennent la version 3D.
Plaxis 8.6 est un logiciel permettant d’analyser des problèmes élastiques, élastoplastique,
élasto-viscoplastiques en 2D. En effet, ces problèmes peuvent conduire à des analyses
d’ouvrages ainsi que des essais géotechniques voire des modélisations afin de suivre ses
évolutions. Le choix de Plaxis dans cette étude est expliqué par le fait qu’il peut résoudre ces
problèmes en utilisant la Méthode des Eléments Finis, et qui est à la fois robuste et convivial,
permettant de traiter des problèmes géotechniques réelles, dans un délai raisonnable en utilisant
des modèles de comportement de sol.
V.2 Présentation du logiciel [26]

L’interface d’utilisation de Plaxis se divise en quatre (04) sous-programme :
 Plaxis Input ;
 Plaxis Calculation ;
 Plaxis Output ;
 Plaxis Curves.
69
V.2.1 Plaxis Input
C’est un programme d’entrée des données qui permet à l’utilisateur de créer un modèle
numérique et de spécifier les propriétés des matériaux concernés et surtout de définir les conditions
initiales et les conditions aux limites. Pour faire cela, les étapes à suivre sont les suivantes :
 dessin des contours géométriques ;
 ajout des couches de sol et les éléments de structure ;
 définition des conditions aux limites ;
 application des charges.
V.2.1.1 Réglages généraux (General settings)
Pendant la création d’un nouveau projet, la fenêtre des réglages généraux apparait.
Elle contient deux onglets project et dimensions. L’onglet project contient le nom et la
description du projet ainsi que le type du modèle à utiliser. En revanche, l’onglet dimensions
permet à l’utilisateur de choisir les unités de longueur, de force et du temps et de dimensionner
la planche à dessin [04].
Sur la figure ci-dessous montre les contenus de l’onglet projet dans lequel
(1) désigne le nom du projet ;
(2) : le nom du modèle dont
 plane strain : est à choisir lorsque la section est plus ou moins uniforme avec
chargement uniforme ;
 axisymetric : est à choisir lorsque les structures sont circulaires (comme les
tunnels) et les charges sont reparties autour de l’axe central.
(3): éléments permet de choisir le nombre des nœuds :
 15-node, permet l’utilisation du triangle à 15 nœuds qui produit un résultat en
contrainte de haute qualité sur différents problèmes ;
 6-node, utilise un triangle à 06 nœuds qui est un élément relativement précis
donnant un résultat satisfaisant pour les analyses standard en déformation.
70
Figure 60 Les contenus de l’onglet « project » [26]
(4) acceleration, dont :

 x-acceleration et y-acceleration sont des accélérations indépendantes afin de modéliser
les efforts dynamiques ;
 earth gravity : c’est la gravité terrestre.
V.2.1.2 La création géométrique
Après que les valeurs nécessaires sont entrées dans la fenêtre general setting,
l’utilisateur peut procéder à la création géométrique.
Figure 61 Fenêtre de la création géométrique [26]
(8) origine
71
(9) –éléments géométriques, composés de :

ligne géométrique ;
plaque, servant à modéliser les structures élancées placées dans le sol,
l’influence des murs, coques ou soutènement s’étendant vers l’axe z ;
géogrille, pour modéliser des renforcements de sol ;
interface, pour la modélisation de l’interaction ( exemple : palplanche) ;
élément nœud à nœud , servant à modéliser des liaisons entre deux points ;
ancrage à tête fixe, pour modéliser une liaison à partir d’un seul point ;
tunnel, pour la création de la section d’un tunnel à section circulaire ou non.
(10) blocage, composé des
blocages standards, permettant d’imposer automatiquement un réglage des
conditions aux limites générales du modèle géométrique ;
blocage en rotation, pour la fixation de degré de liberté en rotation d’un
plaque autour de l’axe z ;
(11) chargement
déplacements imposés, placés sur certaines lignes géométriques pour
contrôler leur déplacement ;
charges reparties ;
charges ponctuelles ;
(12) drains et puits
drains ;
puits ;
(13) propriétés des matériaux
bases de données des matériaux ;
(14) génération du maillage
maille, pour la division des éléments afin de permettre le calcul par éléments
finis.
V.2.2 Plaxis Calculation
Une fois que le modèle soit défini, l’utilisateur peut commencer le calcul proprement dit, tout
en faisant en sorte que le type de calcul à réaliser dans Plaxis soit précis.
En effet, le processus de calcul dans Plaxis est divisé en étape de calcul, à savoir :
72
 l’activation de charge ;
 simulation d’étape de construction ;
 introduction d’une période de consolidation ;
 calcul d’un coefficient de sécurité.
Après avoir cliqué sur l’icône calculate depuis la fenêtre des conditions initiales du
programme Input, la fenêtre principale du programme de calcul apparait.
V.2.3 Plaxis Output
Une vaste gamme d’outils permettant d’afficher les résultats d’une analyse aux
éléments finis est offerte aux utilisateurs par le logiciel Plaxis. Les résultats d’un calcul
d’élément finis sont les déplacements aux nœuds, les efforts dans les éléments de structure et
les déformations. Les résultats de base offerts par Plaxis sont la déformation (deformation) et
les contraintes (stress). A part cela, on peut visualiser toutes les données numériques
correspondant à chaque type de résultats graphiques sous forme de tableau.
V.2.4 Plaxis Curves
Un des avantages que Plaxis offre à ses utilisateurs est qu’il possède un programme
permettant de dessiner des courbes qui relient les charges et les déplacements, contraintes et
déformations pour les points choisis dans la géométrie.
V.3 Les modèles de comportement intégrés dans Plaxis 8.6

Les modèles de comportement de sols intégré dans Plaxis :
 le modèle élastique linéaire ;
 le modèle de Mohr-Coulomb ;
 le modèle spécifique aux sols mous.
V.3.1 Le modèle élastique linéaire
Ce modèle est utilisé pour des structures rigides massives placées dans le sol. D’ailleurs, il
est très limité pour la simulation des comportements d’un sol et il comporte deux paramètres
de rigidité élastique : le module d’Young E et le coefficient de Poisson ν.
On sait que le module d’Young est lié avec le module de cisaillement et le coefficient
de Poisson ν par la relation :
73
𝑬 = 𝟐𝑮(𝟏 + 𝝂) (23)
V.3.2 Le modèle de Mohr-Coulomb
En ce qui concerne ce modèle, il demande la détermination des paramètres, fournies

généralement par des essais en laboratoire, suivante :
 le module d’Young E, qui mesure la rigidité des matériaux et est constante dans le
modèle de Mohr-Coulomb. Il varie en fonction de la déformation et de la contrainte
moyenne, et est obtenu par la courbe des contraintes en fonction des déformations après
la réalisation d’un essai en traction.
le coefficient de Poisson ν, dont les valeurs prises sont généralement comprises entre
0,2 et 0,4.Dans Plaxis, la valeur par défaut de ν est 1,15.
 la cohésion c [kN/m2], dont les valeurs à prendre même s’il s’agit des matériaux
purement frottant sont comprises entre 0,2 et 1 kPa[28] ;
 l’angle de frottement  [°], qui est déterminé par des essais de compression et varie de
15° à 45° suivant la nature des matériaux. Celui-ci est constant et indépendant de la
contrainte moyenne.
 l’angle de dilatance ψ en [°] , qui règle le comportement non associé du sol. Il est vital
de mentionner que lorsque le matériau atteint sa plasticité, la dilatance n’existe plus
comme le cas des argiles ou des sables de densité faible sous assez forte contrainte.
V.3.3 Les modèles spécifiques aux sols mous
V.3.3.1 Modèle pour sols mous ou Soft Soil Model (S.S.M)
Il s’agit d’un modèle élastoplastique avec une surface de charge sous laquelle le
matériau reste élastique. Ce modèle se base sur la prise en compte de l’effet d’une pression
moyenne et la diminution de la teneur en eau pour que l’argile devienne plus résistante.
V.3.3.2 Modèle pour sols mous avec effet du temps ou Soft Soil Creep Model (S.S.C.M)
Ce modèle simule le comportement des sols mous comme les argiles, tourbes normalement
consolidées en fonction du temps.
74
V.3.4 Le modèle de sol avec écrouissage ou Hardening Soil Model (H.S.M)
Le modèle H.S.M prend en compte l’écrouissage en compression pour simuler le

compactage irréversible d’un sol sous son premier chargement [25]. Il permet de simuler le
comportement des sables et des sols mous comme les argiles et les limons. D’ailleurs, c’est
un modèle qui a pour but d’améliorer celui de Mohr-Coulomb en prenant en compte
 de l’évolution du module de déformation lorsque la contrainte augmente ;
 de l’évolution non linéaire du module lorsque le cisaillement augmente .
V.4 UTILISATION DU PLAXIS 2D 8.6
V.4.1 Les étapes principales de la modélisation
Puisque la modélisation en 3D nécessite de temps de calcul assez important, on réalise

les calculs en 2D même si l’effet du creusement d’un tunnel sur les contraintes et les
déformations du massif est tridimensionnel. La modélisation en 2D se base sur le principe de
la méthode de convergence-confinement
Figure 62 Relation entre la pression sur le contour du tunnel et le tassement de

surface [12]
. La pression initiale qui est égale aux contraintes initiales du terrain diminue
progressivement jusqu’au passage du front à l’aplomb de la section. A ce moment, un tassement
Sfront apparait en surface et la charge se repartie, ensuite, sur le soutènement mis en place
derrière le front de taille. Ainsi, on atteint une valeur de tassement égale à Sfinal et un état
d’équilibre est atteint.
75
Par le biais d’une construction par étape du tunnel (qui s’appuie sur le concept du taux de
déconfinement λ), le code Plaxis 2D permet d’intégrer la méthodologie susmentionnée. Les
étapes principales de cette approche sont les suivantes :
 initialisation des contraintes : on calcule les contraintes initiales du massif à partir des
caractéristiques des terrains (γ, K0) et des contraintes hydrauliques ;
 excavation du tunnel ;
 mise en place des soutènements.
V.4.2 Les propriétés géotechniques des matériaux
Le tableau suivant va nous présenter les propriétés géotechniques des sols, que nous
supposons, rencontrés sur le site. Les valeurs prises dans ce tableau sont des valeurs moyennes
trouvées dans des études et recherches déjà faites.
Tableau 10 Propriétés géotechniques des matériaux [17],[31],[21],[27], [11] [39]
Paramètre Nom Granite Granite Charnockite Limon Limon Limon Unité

altéré sain sableux argileux Argileux
micacé
Modèle du MC MC MC MC MC MC MC -
matériau
Poids  13 14 12.95 14 16 16 kN/m3
volumique
(non saturé)
Poids  𝑠𝑎𝑡 24.5 25 24.35 18 19 19 kN/m3
volumique
(saturé)
Perméabilité kx 10-5 10-9 10-8 10-4 1.7*10-7 10-7 m/s
–horizontale
Perméabilité ky 10-5 10-9 10-8 10-4 10-7 10-7 m/s
verticale
Module Eref 8,250.106 12.106 12.106 104 104 104 kN/m2
d’Young
76
Coefficient Ν 0,25 0,3 0.21 0.33 0.35 0.3 -

de Poisson
Cohésion cref 420 441 439 4 4 4 kN/m2
Angle de  41 41 40 25 20 22 °
frottement
interne
Angle de ᴪ 10 12 10,5 2 0 0 °
dilatance
V.4.3 Les caractéristiques mécaniques des soutènements, et des

infrastructures sus-jacentes
Le tableau ci-après nous présente les caractéristiques (valeurs moyennes) des

infrastructures au-dessus du tunnel à creuser. Dans cette étude, on prend ces dernières comme
une charge linéairement repartie, ayant comme propriétés données par le tableau suivant.
Tableau 11 Caractéristiques des infrastructures sus-jacentes [21],[22]
Paramètre Nom Infrastructures Soutènement Unité
Modèle Modèle Elastique Elastique -
Rigidité EA 5.108 1,5.107 kN/m

normale
Rigidité de EI 9.108 1,5.105 kNm2/m
flexion
Coefficient de Ν 0 0,15 -
Poisson
Charge Q 160 8,5 kN/m2
V.4.4 Modélisation numérique
Dans cette modélisation, nous allons montrer les trois sections différentes traversées
par le tunnel, que nous allons appeler :
 section 1 : section sol composé du limon argileux et des limons argileux micacés ;
 section 2 : composée de limon argileux, argile sableuse, du granite altéré et de
granite saint ;
77
 section 3 : composé de limon argileux, de limon sableux et de charnockite.
V.4.4.1 Hypothèses de modélisation
La structure sur la section 1 est représentée par un portique de 70 m de longueur et 10 m

de hauteur. Cette dernière est soumise à une charge linéairement repartie égale à 90kN/m2.
Quant à la section 2 et la section 3, la structure est un portique de longueur respective égale à
100m et 50m et de même hauteur qui vaut 15 m et sont soumises à une charge linéairement
repartie qui vaut 130 kN/m2.
Figure 63 Géométrie du model de la section 1
78
V.4.4.2 Phases de calcul
Etant donné que le tunnel sera creusé sous des infrastructures déjà en place, alors les
phases de creusement du tunnel sont les suivantes :
 phase 0 : détermination des contraintes effectives initiales en calculant le
coefficient du sol au repos K0= 1-sinφ (formule de Jaky) ;
 phase 1 : mise en place des infrastructures sus-jacentes ;
 phase 2 : creusement du tunnel ;
 phase 3 : mise en place des soutènements.
V.4.4.3 Les résultats obtenus
V.4.4.3.1 Section 1
Figure 66 Déplacement horizontal (Valeur maximale : 1,26 . 10-3 m)
79
Figure 67 Déplacement vertical (déplacement maximal : -8,09. 10-3 m)
Figure 68 Maille déformée (déplacement maximal : -8,09. 10-3 m)
Le tableau suivant nous résume les valeurs de quelques paramètres importants.

Tableau 12 Résumé de quelques valeurs des paramètres importants
Contrainte Deplacement Deplacement Force Force de Moment

effective vertical moyen Horizontal axiale cisaillement maximal
2
[kN/m ] [m] moyen [m] [kN/m] [kN/m] [kNm/m]
-214,83 -8,09 . 10-3 1,26 . 10-3 -438,08 -1,49 -2,31
80
V.4.4.3.2 Section 2
Figure 69 Maille déformée (déformation maximale : 305,6.10-6 m)
Le tableau suivant nous résume les valeurs de quelques paramètres importants.
Tableau 13 Résumé des valeurs de quelques paramètres
Contrainte effective Deplacement vertical Deplacement Horizontal

[kN/m2] moyen moyen [m]
[m]
-383,71 305,6.10-6 -60,62.10-6
V.4.4.3.3 Section 3
Figure 70 Déformation des mailles (Déplacement maximal : 140,88 .10-3 m
81
Tableau 14 Résumé des valeurs de quelques paramètres
Contrainte effective Deplacement vertical Deplacement Horizontal

2
[kN/m ] moyen moyen [m]
[m]
-532,33 -140,88. 10-3 5,42.10-3
V.5 Observation et interprétation des résultats

Les résultats précédents nous ont donné des déplacements horizontaux, verticaux et
des contraintes effectives.
 Contraintes effectives
Les contraintes effectives moyennes dans la section 1, section 2 et section 3 ont des valeurs
respectives égales à -214,83 kN/m2, -383,71 kN/m2, -532,33 kN/m2.
 Déplacements horizontaux
Les déplacements horizontaux se produisent dans le massif pendant le creusement du tunnel.

Pour chaque section, les valeurs respectives de ces derniers sont égales à 1,26. 10-3m,
-60,62.10-6m et 5,42.10-3m. On constate que ces déplacements horizontaux sont de types
compression.
 Tassement
Le tassement en surface, qui est symétrique par rapport à l’axe du tunnel, se présente
sous forme de cuvette tridimensionnelle dont les dimensions dépendent du type du sol, de la
proximité de la surface, du diamètre du tunnel [40]. Les déplacements verticaux moyens trouvés
dans chaque section ont des valeurs respectives égales à -8.09 mm, 0.30mm, -1.40 mm. Ces
déplacements verticaux sont de type décompression. En ce qui concerne la section 2 (section
ayant de la roche au plancher et des sols en voûte), on observe un déplacement qui est quasi-
zéro.
Il est vital de mentionner que les tassements admissibles en surface sont similaires à
ceux pour les fondations : de 25 à 40mm. De plus, d’après l’AFTES en 1995, la valeur absolue
des tassements est aujourd’hui de l’ordre de 10 à 20mm. Donc, ces valeurs de tassement sont
acceptables.
82
V.6 Comparaison de ces résultats avec ceux de la méthode

empirique
Nous avons pris la méthode de Sagaseta [1980], de pouelloz [1984]. Les détails de ces
méthodes sont résumés dans l’annexe 4-f.
Les résultats obtenus par ces méthodes sont résumés dans le tableau suivant.
Tableau 15 Valeur du tassement maximal la méthode de Sagaseta (1980)
Tassement
Poids volumique Diamètre du. Mod Coefficient de maximal
saturé du son Tunel [m] Young Poisson [m]
19 5.25 10000 0.3 0.028
16 5.25 10000 0.3 0.024
Tableau 16 Valeur du tassement maximal la méthode de Poupelloz (1984)
Hauteur du centre Diamètre Coefficient Tassement

du tunnel [m] du Tunnel k Coefficient n maximal [m]
11.5 5.25 2.32 0.0574 0.037
Les détails de calcul concernant la méthode de Marcos (courbe de Gauss) est dans l’
annexe 5-a.
Les valeurs de tassement maximal obtenus à travers chaque méthode empirique sont comprises
entre 20mm et 30mm. Par contre, celles données par la méthode numérique nous donne des
valeurs comprises entre 2mm et 8mm. Ces différences de valeurs sont dues au choix des
différents paramètres de calcul pour les méthodes empiriques.
Malgré tous cela, ces valeurs de tassement sont toujours acceptables si on se réfère aux valeurs
admissibles qui sont comprises entre 20mm et 40mm. En revanche, si on se réfère à ce que
l’AFTES a annoncé en 1995, seule la méthode de Marco qui est fiable.
83
V.7 Recommandations
Les études faites dans cet ouvrage nous a permis de conclure que les résultats obtenus
par la modélisation numérique sont fiables si on se réfère à l’intervalle de valeurs de tassements
admissibles imposées par l’AFTES. Pourtant, pour qu’on obtienne des résultats beaucoup plus
proches de la réalité, on préconise les solutions suivantes :
 mettre en œuvre des modèles en trois dimensions en utilisant Plaxis 3D ou avec
FLAC3D pour une meilleure simulation numérique ;
 préciser les valeurs de charges dues aux ouvrages sus-jacents pour chaque section ;
 préciser les profondeurs des fondations des infrastructures sus-jacentes ;
 préciser les valeurs des paramètres géotechniques utilisés dans la modélisation
numérique.
CONCLUSION PARTIELLE
Les mouvements engendrés par le creusement d’un tunnel sont liés à un phénomène
complexe qui dépend de la nature du sol, des méthodes de creusement et de la géométrie de
l’ouvrage à réaliser. On peut dire que le logiciel Plaxis 2D est un outil puissant pour la
modélisation numérique. Il nous a permis d’obtenir des valeurs de tassement beaucoup plus
faible par rapport à celles données calculées à travers les méthodes empiriques. La modélisation
qu’on a faite dans ce chapitre n’est qu’une simplification du phénomène réel. Les résultats
obtenus demandent donc à être confrontés aux observations enregistrées in situ.
84
CONCLUSION
En guise de conclusion, on a pu découvrir à travers des études bibliographiques que

le dimensionnement des tunnels se basent sur des formules empiriques. Ces dernières
peuvent classifier les roches selon des critères précises pour pouvoir dimensionner l’ouvrage
souterrain ainsi que son soutènement. A part cela, pour qu’il y ait plus de précision et qu’on
puisse prendre la décision à propos des méthodes de creusement du tunnel, de son
soutènement, et surtout son tracé, des études géologiques et géotechniques doivent être
faites. D’ailleurs, en connaissant les formations existantes avec leurs caractéristiques
respectives obtenues par les études susmentionnées de l’endroit favorable au creusement du
tunnel, on peut modéliser le creusement en utilisant la méthode des éléments finis. Cette
méthode est maintenant utilisée dans des programmes informatiques afin de faciliter le
calcul. Dans notre cas, nous avons utilisé Plaxis 2D version 8.6 version d’essai, qui nous a
permis d’envisager l’effet du creusement sur les infrastructures existantes.
Le creusement d’un tunnel entraine une modification de l’équilibre préexistant dans

le terrain. Cette modification provoque des mouvements qui sont répercutés en surface et
peuvent être la source d’importants dégâts, surtout lorsqu’il s’agit d’un ouvrage réalisé en
site urbain à faible profondeur comme notre cas. En partant sur des de calculs, nous avons
pu conclure un tassement pas assez important provoqué par le creusement du tunnel. De
plus, l’utilisation d’une machine à attaque ponctuelle est préconisée pour le creusement de
ce tunnel routier.
La solution préconisée dans cette étude est le creusement d’un tunnel dans la zone
où la circulation est intense. Nous avons fait la modélisation de ce tunnel et cela nous a
permis de déduire que le tracé choisi entre Analakely et Ampatsakàna, passant
essentiellement dans des roches granitiques, convient. Des données géologiques,
géophysiques, hydrogéologiques et géotechniques trouvées à travers des études
bibliographiques nous permis d’avoir les résultats trouvés dans cet ouvrage.
Nous pouvons dire que le creusement de ce tunnel reliant Ampantsakàna et

Analakely est une solution efficace pour pallier au problème des embouteillages dans le
centre-ville. Le fait est que le problème de la circulation n’affecte non seulement le centre-
ville mais aussi la périphérie de celui-ci. Alors, au lieu d’exproprier les habitants de leurs
biens pour l’élargissement des routes dans des endroits comme Andravoahangy,
85
CONCLUSION
Ambohibao, n’est-il pas une solution efficace de creuser des tunnels pour remédier au
problème des embouteillages dans ces zones ?
86
BIBLIOGRAPHIE
[01] Abdelghani SEGHIR, Docteur en Sciences de l’université à Mira, Béjaia Algérie.

Cours de Méthode des Eléments Finis.2005-2014,106p.
[02] ARISOA Rivah Kathy. Définition des paramètres géologiques en vue de la
modélisation géophysique du socle cristallin à la latitude d’Antananarivo. Mémoire de fin
d’étude en vue de l’obtention du diplôme d’Ingénieur de l’ESPA. 14 Avril 2000,59 p.
[03] BERABEZ Amirouche et BOURBET Hassan. Dimensionnement du soutènement
d’un tunnel routier Cas du tunnel sis à SIDI AICH. Mémoire de fin de cycle en vue
d’Obtention du diplôme de Master En Génie Civil. Promotion 2015, 105 p.
[04] Catherine Lavergne. Modélisation numérique de tunnel de métro dans le massif
rocheux sédimentaire de la région de Montréal. Mémoire présenté en vue de l’obtention du
diplôme de maîtrise en sciences appliquées. Département des génies civil, géologique et des
mines école polytechnique de Montréal, Avril 2016,139 p.
[05] Centre d’Etude des tunnels. Prix des tunnels, Mars 2016, 178p.
[06] COSTECALDE Vincent. Recyclage cpt, 54 p.
[07] DAMY Clotaire. Cours de fonçage de puits et creusement des galeries, 2016
[08] Dictionnaire de Géologie. 335p.
[09] Dossier pilote des tunnels génie civil. Section 3 : conception et dimensionnement,
ministère de l’équipement, des transports et du logement direction des routes. Juillet 1998,
45 p.
[10] François MARTIN. Mécanique des Roches et Travaux Souterrains Cours et exercices
corrigés. Janvier 2012, 88p.
[11] Jean-Pierre Janin. Tunnels en milieu urbain : Prévisions des tassements avec prise en
compte des effets des pré-soutènements (renforcement du front de taille et voûte-parapluie).
Thèse pour obtenir le grade Docteur.2012, 268 p.
[12] Kitchan Fethi. Etude numérique de la stabilité d’une section du tunnel T4 de
l’autoroute Est-Ouest. Mémoire en vue d’obtention du diplôme de Magistère en Génie Civil.
28 Juin 2012, 165p.
[13] Lionel GENDRE. Introduction à la simulation par éléments finis. Edité le 10/01/13,
07 p.
87
BIBLIOGRAPHIE
[14] Ministère de l’équipement, des transports et du logement direction des routes. Dossier
pilote des tunnels, g é n i e c i v i l - Section 4 : Procédés de creusement et de soutènement.
Juillet 1996, 80p.
[15] Mme Z.Djezzar. Dimensionnement des ouvrages souterrains : Méthodes numériques,
la méthode de calcul aux éléments finis, 22 p.
[16] Module ouvrage souterrains « Les tunneliers ». Ecole nationale des ponts et chaussées
centre des hautes études de la construction, 43 p.
[17] Mr. KEITA Alpha Mamoudou Talibé. Caractérisation et modélisation d’une phase
de creusement d’un tunnel (cas du Tunnel d’Ait Yahia Moussa - wilaya de Tizi-Ouzou).
Mémoire de fin d’étude en vue d’obtention de Master en Mines.187 p.
[18] Philippe MESTAT.Maillages d'éléments finis pour les ouvrages de géotechnique
Conseils et recommandations. Directeur de recherche et chef de la section rhéologie, 26 p
[20] Pierre GESTA.Tunnelier. Ingénieur de l’Ecole Centrale de Paris. 20 p.
[21] Plaxis Version 8 Tutorial Manual, 110p.
[22] Plaxis Version 8, Dynamics Manual, 70p.
[23] RANDRIATSILANY Herilala Antonio Andrien. Modélisation des structures
élastiques par éléments finis volumiques application sous Matlab. Mémoire de Fin d’Etude
en vue d’obtention du Diplôme d’Ingénieur en Bâtiment et Travaux Publics. Promotion
2009, 243 p.
[24] Remi Roger FABIEN, cours des ouvrages souterrains, 2017
[25] RAZAFIADJAVOLA Maeva Lalantenaintsoa. Analyse et modélisation des risques
de glissements de terrain dans les Zones Hautes d’Antananarivo-Ville. Mémoire de fin
d’études en vue de l’obtention du diplôme de Master. Septembre 2015,140p.
[26] RAZAFINTSALAMA ZAFINIMANGA Raivohavana.Application du logiciel
Plaxis sur la stabilité des talus. Mémoire de fin d’étude en vue d’obtention du diplôme
d’ingénieur des Mines. Promotion 2010, 143p.
[27] TCHOUANI NANA et M. CALLAUD. Cours de mécanique des sols Tome I.
Décembre 2004, 137p.
[28] TOUAMRIA SABRI. Modélisation numérique tridimensionnelle de l’effet de
creusement d’un tunnel sur le comportement des superstructures sus-jacentes. Mémoire pour
l’obtention du diplôme de Magistère en Génie Civil Option Géotechnique,98 p.
[29] Some Useful Numbers on the Engineering Properties of Materials (Geologic and
Otherwise) GEOL 615, 6p.
88
BIBLIOGRAPHIE
[30] Sylviane BERNAT. Modélisation du creusement d’un tunnel en terrain meuble,

Qualification sur chantier expérimental. Thèse pour obtenir le titre de Docteur en Génie
Civil. 1996, 320p.
[31] Sylvie Magdeleine. La Terre dans l'Univers, la vie, l'évolution du vivant.Chap.1 :
Caractéristiques du domaine continental, 73p.
[34] USGS-BRGM. Cartes géologique et Métallogéniques de la République de
Madagascar à 1/1000 000.Projet de gouvernance des ressources minérales de Madagascar.
Notice Explicative 2012. 236p
[36] Vincent Grégoire. Caractérisation pétrostructurale du socle précambrien de
Madagascar (pays betsileo, province de Fianarantsoa) : méthodologie et implications
géodynamiques. Thèse présentée Devant l’Université Paul-Sabatier de Toulouse en vue de
l’obtention du doctorat de l’Université Paul-Sabatier. 26 Novembre 1999. 210p.
[38] Zhi-Qiang FENG et Al. MEF Méthode des Eléments Finis, 19p.
[37] Z.Djezzar, Ingénieure principale des T.P Dimensionnement des ouvrages
souterrains: Méthodes numériques la méthode de calcul aux éléments finis. 22p.
[38] RAONIMAHARIVO R ; « Apport des prospections géophysique à la caractérisation

des domaines hydrogéologiques du haut bassin versant Ikopa, zone du socle cristallin
(Hautes terres de Madagascar) », Mémoire de maitrise, Université d’Antananarivo, (2010)
82p
[39] RAHAMBININTSOA Sedra Avojanahary, Etude géologique et géotechnique de
faisabilité du projet de tunnel, reliant Ampefiloha- Analakely dans la Commune Urbaine
d’Antananarivo, Mémoire de fin d’études en vue de l’obtention du diplôme de Master, 2020,
126p.
[40] Aspect général sur la construction souterraine, Partie I, 30p
89
WEBOGRAPHIE
[19] Méthode des éléments finis. In Wikipédia, disponible sur

https://fr.wikipedia.org/wiki/Méthode des éléments finis/Présentation générale. Consulté le
27 Août 2018.
[32] Tunnel. In Larousse, disponible sur www.Larousse.com.Consulté le 13 Janvier 2017.
[33] Tunnelier. In Wikipédia, disponible sur :
https://fr.wikipedia.org/wiki/Tunnelier.Consulté le 13 Janvier 2017
[35] Tunnel. Meet Crossrail's giant tunnelling machines. In You tube. Consulté le 31 Janvier
2019
90
TABLE DES MATIERES
TABLE DES MATIERES
FISAORANA ............................................................................................................. i
Sommaire................................................................................................................. ii
Liste des tableaux ................................................................................................... iii
Liste des figures...................................................................................................... iv
Liste des abréviations ............................................................................................ vii
Liste des unités ....................................................................................................... ix
Liste des symboles................................................................................................... x
INTRODUCTION ................................................................................................... 1
Chapitre I. GENERALITES SUR LES TUNNELS ............................................. 3
I.1 INTRODUCTION ................................................................................................................ 3
I.2 CLASSIFICATION DES MASSIFS ROCHEUX ET DE DIMENSIONNEMENT DES TUNNELS..... 3
I.3 LES ETAPES D’INVESTIGATIONS AVANT LE CREUSEMENT D’UN TUNNEL .................... 16
Chapitre II. METHODES DE CREUSEMENT DES TUNNELS ..................... 19
II.1 LES PROCEDEES DE CREUSEMENT ................................................................................ 19
II.2 ABATTAGE A L’EXPLOSIF .............................................................................................. 21
II.3 METHODES DE CREUSEMENT MECANISEES .................................................................. 26
II.4 LES DIFFERENTS TYPES DE SOUTENEMENT .................................................................. 33
Chapitre III. LA METHODE DES ELEMENTS FINIS .................................... 38
III.1 HISTORIQUE.............................................................................................................. 38
III.2 QUELQUES RAPPELS THEORIQUES ........................................................................... 38
III.3 LA DISCRETISATION ................................................................................................. 41
III.4 LES DIFFERENTES PHASES DE CALCUL..................................................................... 43
Chapitre IV. METHODOLOGIE ....................................................................... 53
IV.1 INTRODUCTION ......................................................................................................... 53
IV.2 LES TUNNELS ROUTIERS A ANTANANARIVO ............................................................ 53
IV.3 LOCALISATION ADMINISTRATIVE DE LA ZONE D’ETUDE......................................... 54
IV.4 ETUDES BIBLIOGRAPHIQUES .................................................................................... 54
IV.5 TRACE DU TUNNEL.................................................................................................... 65
91
TABLE DES MATIERES
IV.6 DESCRIPTION DU TUNNEL......................................................................................... 67
Chapitre V. MODELISATION NUMERIQUE ................................................. 69
V.1 HISTORIQUE DU PLAXIS ................................................................................................ 69
V.2 PRESENTATION DU LOGICIEL [26]................................................................................... 69
V.3 LES MODELES DE COMPORTEMENT INTEGRES DANS PLAXIS 8.6 ................................. 73
V.4 UTILISATION DU PLAXIS 2D 8.6 ............................................................................ 75
V.5 OBSERVATION ET INTERPRETATION DES RESULTATS .................................................. 82
V.6 COMPARAISON DE CES RESULTATS AVEC CEUX DE LA METHODE EMPIRIQUE ............ 83
V.7 RECOMMANDATIONS..................................................................................................... 84
CONCLUSION ..................................................................................................... 85
BIBLIOGRAPHIE ................................................................................................ 87
WEBOGRAPHIE .................................................................................................. 90
TABLE DES MATIERES .................................................................................... 91
Annexes .................................................................................................................. I
92
Annexe
Annexes
I
Annexe
Annexe 1 Classification des roches selon la résistance à la compression [13]
II
Annexe
Annexe 2 Classification des massifs rocheux par la méthode de N.
Barton [12]
C’est une méthode développée Barton et Al du Norwegian Geotechnical Institute pour la

classification empirique des massifs rocheux. En faisant des analyses de plus de 200 ouvrages
souterrains, ils ont pu proposer un indice Q qui varie de 0.001 à 1000 (Q-system) permettant de
déterminer la qualité des massifs rocheux, de définir le mode de soutènement à mettre en place et
déterminer les pressions qui s’exercent sur le soutènement. L’indice de qualité Q est déterminé par la
relation suivante :
𝑹𝑸𝑫 𝑱𝒓 𝑱𝒘
𝑸=
𝑱𝒏 𝑱𝒂 𝑺𝑹𝑭
Où
RQD est le Rock Quality Designation de Deere ;
Jn désigne l'expression du nombre de familles principales de discontinuités
Jr caractérise la rugosité des faces des joints.
Ja définit le degré d'altération des joints (épaisseur du joint et nature du matériau de
remplissage).
Jw spécifie les conditions hydrogéologiques : importance des venues d'eau et
pression.
SRF ou (Stress Reduction Factor) précise l'état des contraintes dans le massif.
Classification des massifs rocheux par N. Barton

Valeur Qualité du massif
400-1000 Exceptionnellement bonne
100-400 Extrêmement bonne
40-100 Très bonne
10-40 Bonne
4-10 Moyenne
1-4 Mauvaise
0,1-1 Très mauvaise
0,01-0 ,1 Extrêmement mauvaise
0,001-0,01 Exceptionnellement
mauvaise
III
Annexe
Annexe 3 : Classification des roches en fonction de la résistance à la

traction
a. Classification des roches selon sa résistance à la traction
b. Classification des roches selon la RMR
IV
Annexe
Annexe 4 Classification LPC et GTR du sol – Espacement des

discontinuités-Calcul du tassement
a. Classification LPC du sol

 Sols fins dont 50% d’éléments en poids sont inférieurs à 80 ;
b. Classification GTR du sol
V
Annexe
c- Espacement des discontinuités
VI
Annexe
d-Indice de plasticité
e- Matière organique (MO)
VII
Annexe
f- Estimation du tassement maximal en fonction des formules empiriques
VIII
Annexe
Annexe 5 Résultats des essaies piézométriques faits par le
LNTPB et Entreprise SIF
Niveau d'eau Profondeur de la

Point de forage Date (m) piézomètre (m)
S4 2/9/2017 13,7 35,8
S7 24/08/2017 9,46 29,5
S9 18/01/1962 18,5 30
S11 6/1/1969 15,25 60
S12 10/1/1969 13,4 60
S13 23/01/1963 10,75 14,3
S15 21/01/1963 14,5 15,15
S16 11/1/1963 11,6 17,4
(Source : RAHAMBININTSOA Sedra, 2018)
a. Résultat de calcul du tassement par la méthode de Marco (courbe de

Gauss)
i x 2i -x² Smax i² S(x)
1 0.1 2 0.01 0.02 1 0.02010025
0.2 2 0.04 0.02 0.02040403
0.3 2 0.09 0.02 0.02092056
0.4 2 0.16 0.02 0.02166574
0.5 2 0.25 0.02 0.02266297
0.6 0 0.36 0.02 0.02394435
0.7 2 0.49 0.02 0.02555243
0.8 2 0.64 0.02 0.02754256
0.9 2 0.81 0.02 0.02998605
1 2 1 0.02 0.03297443
1.414213562 2 2 0.02 0.05436564
1.732050808 2 3 0.02 0.08963378
V
Titre : Modélisation numérique du creusement d’un tunnel routier reliant
Ampantsakana et Analakely dans la Commune Urbaine d’Antananarivo
Nombre de pages : 120

Nombre de figures : 70
Nombre de tableaux : 16
Nombre d’équations : 23
Résumé
Le problème des encombrements dans les périphéries et surtout dans le centre-ville de
Madagascar ne cesse de s’aggraver. Ce fléau est dû à l’étroitesse des voies de communication liant
plusieurs endroits dans la capital de la Grande Ile. Dans cet ouvrage, notre zone d’influence est à
Analakely où l’intensité de la circulation s’avère élevé et où l’extension des routes semble
impossible. Le recours au creusement d’un tunnel routier liant Analakely et Ampantsakàna est une
solution efficiente pour pallier à ce problème. Les objectifs de cette étude sont le choix du tracée du
tunnel et la modélisation numérique du creusement de ce dernier, afin de connaitre ses effets sur les
infrastructures sus-jacentes, par la Méthode des Eléments Finis en utilisant le logiciel Plaxis 2D
version 8.6. Une étude bibliographique sur la géologie de la zone d’influence a été faite et a abouti à
l’obtention d’un modèle géologique de son sous-sol. D’ailleurs, les valeurs des paramètres
géotechniques qu’on a utilisées sont tirées dans des études et recherches déjà faites. La modélisation
numérique de ce creusement nous a donné des valeurs de tassement assez faibles, et ça nous incite à
confronter ces résultats aux observations in situ.
Mots clés : encombrement, tunnel, Méthode des Eléments Finis, Plaxis 2D, tassement
Abstract
The problems about congestions at the outskirt and especially in downtown area of
Madagascar is getting worst and worst. This scourge is caused by the small wideness of the roads
which communicate many places in the Capital city of this Big Island. In this work, our influence
zone is Analakely, where the intensity of the traffic is high and where the extension of roads seems
to be impossible. Boring a tunnel which links Analakely and Ampantsakàna is an efficient solution
to sort this problem out. The objectives of this study are choosing the layout of the tunnel and making
a numerical modelling using Finite Element Method through Plaxis 2D version 8.5, to find out the
effect of the tunnel boring with the neighboring infrastructures.A bibliographical study about the
influence zone’s geology was done and that ended at having a geological model of its underground.
Moreover, the values of all geotechnical parameters used in the modelling process were taken from
researches and studies. The numerical modelling of this tunnel-boring gave us some settlements’
values and that incites us to compare these results with in situ observation.
Key words: congestion, tunnel, Finite Element Method, Plaxis 2D, settlement
Auteur : RAKOTONDRABE TOJONIRINA Jean Claudio Encadreur : Monsieur Rémi Roger FABIEN
Téléphone : (+261) 34 64 191 78
Enseignant chercheur au sein de l’Ecole
E-mail : [email protected]
Supérieure Polytechnique d’Antananarivo

RakotondrabeTojonirinaJeanC ESPA MAST 19

Transféré par

Droits d'auteur :

Formats disponibles

RakotondrabeTojonirinaJeanC ESPA MAST 19

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

RakotondrabeTojonirinaJeanC ESPA MAST 19

Transféré par

Droits d'auteur :

Formats disponibles

UNIVERSITE D’ANTANANARIVO

ECOLE SUPERIEURE POLYTECHNIQUE

DOMAINE : SCIENCES DE L’INGENIEUR

MENTION : INGENIERIE MINIERE

Mémoire de fin d’étude pour l’obtention du diplôme d’ingénieur Polytechnicien

Parcours : Sciences et Techniques Minières

RAKOTONDRABE TOJONIRINA Jean Claudio

Encadreur : Monsieur Rémi Roger FABIEN

Enseignant Chercheur à l’Ecole Supérieure Polytechnique d’Antananarivo

ANNEE UNIVERSITAIRE : 2017-2018

ECOLE SUPERIEURE POLYTECHNIQUE

DOMAINE : SCIENCES DE L’INGENIEUR

MENTION : INGENIERIE MINIERE

Mémoire de fin d’étude pour l’obtention du diplôme d’ingénieur Polytechnicien

Parcours : Sciences et Techniques Minières

RAKOTONDRABE TOJONIRINA Jean Claudio

Président du jury : Monsieur Léon Félix RANAIVOSON

ANNEE UNIVERSITAIRE : 2017-2018

 Andriamatoa Rijalalaina RAKOTOSAONA, talen’ny Sekoly Ambony

Hotolorana ny fankasitrahana manokana koa ireto olona manaraka ireto:

 ireo namana aman-tsakaiza sy mpiara-mianatra izay nanampy sy nifanakalo hevitra tamin’ny

˜ Mankasitraka anareo rehetra tompoko˜

LISTE DES TABLEAUX………………………………………………………………...iii

LISTE DES FIGURES…………………………………………………………………..iv

LISTE DES ABREVIATIONS……………………………………………………….viii

LISTE DES UNITES…………………………………………………………………..x

LISTE DES SYMBOLES……………………………………………………………….xi

CHAPITRE I GENERALITES SUR LES TUNNELS…………………………………03

CHAPITRE II METHODE DE CREUSEMENT DES TUNNELS……………………19

CHAPITRE III LA METHODE DES ELEMENTS FINIS……………………………38

CHAPITRE IV LOCALISATION DE LA ZONE D’ETUDE…………………………53

CHAPITRE V UTILISATATION DE PLAXIS 2D 8.6……………………………...69

TABLE DES MATIERES…………………………………………………………….…91

dessus d’une cavité [12] ............................................................................................................... 5

Tableau 2 Classification des roches suivant le R.Q.D [12] .................................................................. 5

Tableau 3 Classification des roches en fonction de la résistance à la compression[3] ........................ 7

Tableau 4 Choix du type de soutènement 1974) [3] ............................................................................ 7

Tableau 5 Les modules de rigidité de différents types de soutènement [12] ...................................... 13

Tableau 6 Essais en laboratoire appliqué aux sols et rochers [12] ..................................................... 16

Tableau 7 Les caractéristiques des éléments massifs bidimensionnels ............................................ 44

Tableau 8 Exemples illustrant les significations de chaque matrice [38] ........................................... 51

Tableau 9 Les dimensions du tunnel ................................................................................................ 68

Tableau 10 Propriétés géotechniques des matériaux [17],[31],[21],[27], [11] [39] .......................................... 76

Tableau 11 Caractéristiques des infrastructures sus-jacentes [21],[22]................................................. 77

Tableau 12 Résumé de quelques valeurs des paramètres importants ............................................... 80

Tableau 13 Résumé des valeurs de quelques paramètres ................................................................. 81

Tableau 14 Résumé des valeurs de quelques paramètres ................................................................. 82

Tableau 15 Valeur du tassement maximal la méthode de Sagaseta (1980) ..................................... 83

Tableau 16 Valeur du tassement maximal la méthode de Poupelloz (1984) ................................... 83

Figure 1 Zone décomprimée au-dessus d’une cavité [12] .................................................................... 4

Figure 2 Configuration géologique illustrant les limites du RQD [10] ................................................ 6

Figure 3 Type de soutènement en fonction du RQD d’après LAUFFER [03] ..................................... 6

Figure 4 Les contraintes et les déplacements [12] ................................................................................ 8

Figure 6 Portée non soutenue d'un tunnel [11] ................................................................................... 10

Figure 7 Variation du taux de déconfinement, de la pression et du déplacement radial en fonction

de la distance du front de taille [11]........................................................................................... 11

Figure 8 Front de taille stable (extrusion négligeable) [11]................................................................ 12

Figure 9 Front de taille stable, convergence et extrusion forte [11] ................................................... 12

Figure 10 Figure 10 Front de taille instable [11] ................................................................................ 12

Figure 11 Ligne caractéristique de la cavité [11] ............................................................................... 14