Les documents d’information

LA STABILITÉ DU FRONT
DE TAILLE DES TUNNELS
Mécanismes de rupture, justification
et principes constructifs

MINISTÈRE
CHARGÉ
DES TRANSPORTS
décembre 2022 Centre d’Études des Tunnels
www.cetu.developpement-durable.gouv.fr
AVERTISSEMENT
Les documents d'information ont pour but de fournir une information sur une technique ou un problème nouveau ou insuffisamment traité
par ailleurs. Le lecteur pourra y trouver des repères susceptibles de l'aider dans son activité. Le contenu et les éventuelles conclusions
présentés ne doivent pas être considérés comme des recommandations du CETU. Bien que le maximum soit fait pour s'assurer de la
fiabilité des sources utilisées, la responsabilité du CETU ou des auteurs du document ne saurait être engagée.

2
 Les documents d’information

LA STABILITÉ DU FRONT
DE TAILLE DES TUNNELS
Mécanismes de rupture, justification
et principes constructifs

décembre 2022

Centre d’Études des Tunnels

25, avenue François Mitterrand
69500 BRON – France
Tél. 33 (0)4 72 14 34 00
Fax. 33 (0)4 72 14 34 30
[email protected]
www.cetu.developpement-durable.gouv.fr

3
 TABLE DES MATIÈRES
1 CONTEXTE ET OBJET DU DOCUMENT 5

2 PRINCIPES GÉNÉRAUX POUR L’ÉTUDE DE LA STABILITÉ 7

2.1 Position du problème 7
2.2 Hypothèses et notations 8
2.3 Configurations de calcul 9
2.4 Approche de la sécurité 10
2.5 Approches de modélisation 11

3 ÉVALUATION DES CONDITIONS DE STABILITÉ 13

3.1 Milieu discontinu 13
3.2 Milieu continu non drainé 14
3.3 Milieu continu drainé sans écoulement 17
3.4 Milieu continu drainé avec écoulement 22

4 OBSERVATIONS ET MESURES CONSTRUCTIVES SUR CHANTIER 26

4.1 Procédés constructifs en méthode conventionnelle 26
4.2 Méthodes de suivi en méthode conventionnelle 28
4.3 Procédés constructifs au tunnelier à front pressurisé 28
4.4 Méthodes de suivi au tunnelier à front pressurisé 29

5 SYNTHÈSE 30

6 GLOSSAIRE 32

7 RÉFÉRENCES 33

8 ANNEXE : RETOURS D’EXPÉRIENCES D’INSTABILITÉS DE FRONT 36

4
1 CONTEXTE ET OBJET
DU DOCUMENT
Le front de taille exerce son influence immédiate sur la zone En conséquence, une bonne compréhension des mécanismes
« active » du creusement d’un tunnel en avant et en arrière de déformations du terrain est nécessaire afin d’éviter l’initiation
du front de taille (AFTES GT16, 2018). Repoussée à chaque de mécanismes de rupture ou d’en prévoir les conséquences
pas d’avancement, cette zone est le siège de la majorité des quand ceux-ci n’ont pu être évités, en réponse à quatre enjeux :
réorganisations de contraintes dans le terrain.
• assurer la sécurité du chantier (compagnons en premier lieu)
Pour les tunnels profonds, la réorganisation des contraintes et de ses abords (constructions avoisinantes en surface
permet, sans ou avec interventions au front en fonction de la notamment) ;
qualité du terrain, la mobilisation d’effets de voûte qui favorisent • permettre la mise en place dans de bonnes conditions
le développement d’une nouvelle situation d’équilibre. Pour les des structures de pré-soutènement à l’avant et de
tunnels à faible profondeur, ces effets de voûte sont plus difficiles soutènement à l’arrière du front de taille ;
à mobiliser et ainsi les déformations, si elles sont mal maîtrisées, • participer à la maîtrise du coût et du planning du projet,
peuvent donner lieu à des tassements en surface sources en limitant les hors-profils et en évitant les situations de
d’endommagement des avoisinants. L’amplitude des déformations déformations ni anticipées ni maîtrisées ;
dans des terrains médiocres à grande profondeur peut s’accom- • contribuer à la pérennité de l’ouvrage réalisé en limitant
pagner d’une zone décomprimée qui se développe au pourtour la dégradation du terrain encaissant en périphérie de
de l’excavation et au sein de laquelle peuvent se manifester des l’excavation et donc les efforts qui s’exercent sur les
phénomènes gravitaires tout comme à proximité de la surface. structures de revêtement définitif.

De nombreuses instabilités ont été observées dans des contextes Ce document d’information a pour objectif de proposer une
géotechniques variés. CEDD (2015) dresse par exemple une synthèse générale des mécanismes de rupture, des approches
liste (non exhaustive) de plus de soixante instabilités survenues de modélisation disponibles et de leurs domaines d’application,
entre 1977 et 2014, au front ou en section courante. Les fiches ainsi que des outils de calcul les plus adaptés à chaque projet.
jointes en annexe complètent cette liste à travers quelques En effet, ce sujet n’a pour l’instant pas fait l’objet de norme
exemples assez récents. Les constats suivants en découlent : (norme d’application de l’Eurocode 7 par exemple) ni de recom-
mandation française (AFTES) ou internationale (ITA-AITES)
• les instabilités de front sont possibles dans des contextes pour les ingénieurs. Le document qui s’en rapproche le plus
géotechniques variés (fronts homogènes ou mixtes, terrains actuellement est la recommandation de la DAUB (homologue
rocheux fracturés, terrains meubles de type sols, etc.), allemand de l’AFTES) publiée en 2016, mais qui reste spécifique
et pour toutes les méthodes de creusement (méthode au cas des tunneliers, et qui introduit des spécificités nationales.
conventionnelle, tunnelier ouvert, tunnelier à pression de
terre ou à pression de boue, etc.) ; Il est principalement destiné aux ingénieurs de bureaux d’études,
• en souterrain, les volumes mobilisés peuvent aller de en maîtrise d’oeuvre et en entreprise. En proposant une large
quelques décimètres cubes pour des instabilités locales à synthèse bibliographique, les différents mécanismes d’instabilité
plusieurs centaines de mètres cubes pour des instabilités sont présentés, et structurés en fonction de la nature des terrains
globales intéressant le front dans son ensemble voire (terrains rocheux plus ou moins fracturés, terrains meubles de
les terrains de couverture ; type sols grenus ou sols fins, etc.) et de la méthode de creuse-
• l’allure des mécanismes dépend de la nature des terrains : ment (méthode conventionnelle avec ou sans renforcement du
blocs délimités par le réseau préexistant en massif rocheux, front ou présoutènement, tunnelier pressurisé ou non). Il peut
mécanismes localisés à la proximité du front en terrain ainsi être utile également aux ingénieurs chargés du suivi de
pulvérulent (avec évolution régressive vers la surface) et travaux (chargé de soutènement côté entreprise, supervision
mécanismes plus volumineux en terrain argileux cohérents ; géotechnique côté maîtrise d’œuvre) pour comprendre le com-
• les instabilités induites par une maîtrise partielle ou non portement du terrain à l’excavation afin de choisir et adapter les
appropriée de la stabilité du front de taille peuvent se meilleures dispositions constructives en réponse aux mécan-
manifester à la surface avec un effet retardé dans le ismes rencontrés.
temps et dans l’espace ;
• les conséquences des instabilités de front sont très Pour chaque configuration, les différentes méthodes de calcul
variables, qualifiées de « quasiment négligeables » à pour évaluer la stabilité du front sont décrites, et une formulation
« très importantes », avec des retards (jusqu’à plusieurs empirique explicite à vocation « opérationnelle » est proposée.
mois) ou des surcoûts importants (jusqu’à plusieurs Ce document a ainsi pour vocation d’être utilisable facilement par
millions d’euros) et des décès (les ouvriers en souterrain les ingénieurs en donnant une vision synthétique des différentes
étant les premiers exposés). méthodes existantes et de leur domaine d’application.

5
1
Ce document n’aborde pas la question essentielle que constitue Si les études conduisent à retenir des conditions d’instabilité
le choix des paramètres de calcul, en particulier les paramètres du front de taille ou si bien sûr une instabilité est rencontrée
de résistance au cisaillement des terrains. Conformément à pendant les travaux, différents procédés constructifs sont
l’Eurocode 7-1 (AFNOR, 2005), ce choix doit s’appuyer sur envisageables pour améliorer la stabilité du front ou la
les valeurs mesurées, les valeurs dérivées des essais en sécurité des travaux. Les grands principes des procédés
place et en laboratoire, complétées par les enseignements de utilisables en méthode conventionnelle et au tunnelier sont
l’expérience. Pour les problématiques de stabilité, une estimation ainsi présentés, en parallèle du suivi nécessaire pendant
raisonnablement prudente doit être effectuée, notamment les travaux afin de conforter les hypothèses géotechniques
concernant la cohésion effective. En cas d’incertitude importante, effectuées lors des études de stabilité et d’évaluer les conditions
des calculs « en fourchette » doivent être réalisés. réelles de stabilité.

6
2 PRINCIPES GÉNÉRAUX
POUR L’ÉTUDE DE LA STABILITÉ

2.1 POSITION DU PROBLÈME

On considère le cas idéal d’un tunnel parfaitement soutenu sur L’Illustration 1 en donne une vision schématique et exprime
sa section courante, avec un décousu nul tel que seul le front le lien avec les mouvements en surface dans le cas d’un tunnel
de taille, sur lequel s’exerce une certaine pression, est libre à faible profondeur hors nappe :
de se déplacer. • si la pression exercée sur le front est inférieure à la
contrainte horizontale géostatique (σT < σh0), un mouvement
d’extrusion du front apparaît, et peut conduire à
Dans le cas d’un tunnelier à front pressurisé, cette pression
l’apparition de tassements en surface. Le front peut
« réelle » est exercée par de l’air comprimé (pression d’air), éventuellement entrer en rupture par effondrement,
le matériau excavé (pression de terre) ou de la bentonite avec apparition d’un fontis dans un cas extrême ;
(pression de boue, avec boue conventionnelle ou « lourde »). • a contrario, si la pression exercée sur le front est supérieure
En méthode conventionnelle, une pression de confinement à la contrainte horizontale géostatique (σT > σh0), un
« équivalente » peut être considérée pour prendre en compte mouvement de refoulement apparaît par mise en butée
l’apport d’un boulonnage du front. des terrains à l’avant du front, ce qui peut conduire à un
soulèvement en surface (ce cas n’apparaît toutefois possible
qu’au tunnelier à front pressurisé sous très faible couverture).
En fonction de l’amplitude de cette pression, des mouvements
d’extrusion (vers l’arrière du front de taille) ou de refoulement Dans une approche globale d’évaluation du domaine de pression
(vers l’avant du front quand on le regarde depuis l’intérieur acceptable sur un projet, on identifie un « domaine de stabilité »
(en rouge sur l’Illustration 1) correspondant à la non-atteinte
du tunnel dans le sens de l’avancement) peuvent apparaître.
de mécanismes de rupture par effondrement ou refoulement
Leur amplitude dépend notamment de la géométrie de du front de taille. À cela s’ajoute l’évaluation d’un « domaine
l’ouvrage, de l’état initial des contraintes et des caractéristiques de déplacements admissibles » (en vert sur l’Illustration 1)
géotechniques des terrains. conditionné par la sensibilité des avoisinants.

Illustration 1 : Liens entre pression frontale, mécanismes d’instabilité au front de taille et déplacements générés au droit des avoisinants
en cas de tunnel à faible profondeur.

7
2
Si par ailleurs la pression d’air, de terre ou de boue exercée sur le numériques. Ce point n’est pas traité ici dans la mesure où il
front ne permet pas de contrebalancer les pressions interstitielles s’agit d’une problématique à part entière. Les recommandations
dans le terrain, deux effets liés peuvent apparaître : (i) des forces du groupe de travail n°16 de l’AFTES regroupent des éléments
d’écoulement augmentant l’instabilité du front de taille, et (ii) des sur ce sujet : évaluation de la sensibilité intrinsèque des
volumes d’eau potentiellement importants pouvant conduire à avoisinants, évaluation des déplacements induits par les travaux
l’ennoiement du chantier. souterrains et évaluation du niveau de dommage associé.

Ce document d’information est principalement consacré Le cas de pressions frontales supérieures aux contraintes
à la limite inférieure du domaine de stabilité, c’est-à-dire à horizontales en place s’avère plus simple à traiter d’un point de
la pression minimale à exercer pour empêcher la rupture par vue opérationnel. En effet, les approches théoriques (équilibre
effondrement du front de taille. limite ou calcul à la rupture) consacrées à la rupture par
refoulement du front de taille font l’hypothèse d’une mobilisation
La limite inférieure du domaine « acceptable » vis-à-vis de totale des contraintes de cisaillement le long des surfaces de
la sensibilité des avoisinants dans le cas d’un tunnel à faible rupture, ce qui conduit à des pressions frontales maximales
profondeur repose sur des approches semi-empiriques largement supérieures à celles conduisant à des mouvements
(méthode dite des « pertes de volume ») ou des calculs de soulèvement significatifs en surface (Berthoz et al., 2012).

2.2 HYPOTHÈSES ET NOTATIONS

L’étude de la stabilité du front de taille repose généralement • Méthode de creusement et pression de confinement :
sur les hypothèses suivantes avec les notations synthétisées – le creusement peut être réalisé au tunnelier ou en méthode
en Illustration 2 : conventionnelle. Dans les deux cas, on note vavct la vitesse
moyenne d’avancement du front de taille (intégrant les
• Géometrie : le front de taille est assimilé à un disque de
phases successives de creusement et soutènement) ;
diamètre D. En méthode conventionnelle, l’excavation réelle
– une pression σT constante est considérée afin de
n’est pas circulaire et un disque (voire un rectangle) équiva-
représenter les conditions aux limites sur la surface du front
lent à la surface du front est souvent retenu. Si l’excavation
et les parois latérales le long du décousu. En pratique,
est très élancée (horizontalement ou verticalement), le choix,
on retient la valeur s’exerçant dans l’axe du tunnel. En
plus difficile, s’effectue au regard de la géométrie du
cas de tunnels de très grands diamètres (supérieurs à
mécanisme en jeu. On note C la couverture, H la profondeur
une dizaine de mètres), il convient de tenir compte du
à l’axe du tunnel, et d la longueur du décousu (zone non
gradient de pression frontale sur la hauteur du front, ou
soutenue à l’arrière du front).
a minima, de s’assurer que la pression retenue dans l’axe
• Hydrogéologie : on considère le cas d’une nappe libre
est représentative des conditions réelles de creusement.
dont le toit est situé à une hauteur Hw au-dessus de l’axe du
tunnel. En cas de nappe captive, une hauteur Hw équivalente On recherche par le calcul la valeur minimale de cette pression σT
déduite des pressions hydrostatiques peut être estimée afin telle que la stabilité du front est assurée. Si cette valeur est
de se ramener au même cas d’étude. négative (avec les conventions de signe de la mécanique des sols et
• Géologie et géotechnique : le front, homogène ou composé des roches, c’est-à-dire des contraintes positives en compression),
de différents faciès, peut comporter : le front de taille est considéré théoriquement auto-stable. Si
– des terrains rocheux « peu fracturés », présentant peu de elle est positive, le front n’est pas auto-stable et des procédés de
discontinuités à l’échelle du front et se comportant comme construction adaptés doivent être mis en œuvre pour en assurer
des milieux discontinus à l’échelle de l’ouvrage. Dans ce la stabilité. Ces procédés sont dimensionnés sur la base de cette
cas, la stabilité est conditionnée par le poids volumique γ valeur de pression moyennant une marge de sécurité à définir.
des dièdres instables, l’orientation des discontinuités θ et
la résistance au cisaillement le long des discontinuités,
notamment l’angle de frottement φ’disc ;
– des terrains « meubles » (sols, roches très fracturées),
pouvant être considérés comme continus à l’échelle de
l’ouvrage, ou encore de la roche intacte, continue par
définition. Dans ce cas, on assimile généralement chaque
faciès à un milieu homogène de poids volumique γ, de C
perméabilité k, de comportement rigide parfaitement Terrains
Soutènement H
plastique caractérisé par un critère de Mohr-Coulomb )
(c’, φ’) en conditions drainées, Tresca (Cu) en conditions H OU
non drainées, ou éventuellement Hoek & Brown pour )
tunnelier OU
les terrains rocheux. OU T
• Charge en surface : on note σs la pression constante méthode OU
exercée en surface, qui prend en compte le poids des conventionnelle
bâtiments en surface (en considérant 12 kPa par étage
par exemple), des surcharges de chantier (généralement
de l’ordre de 20 kPa) ou des surcharges d’exploitation
(environ 10 kPa pour une chaussée). Illustration 2 : Position du problème et notations.

8
2.3 CONFIGURATIONS DE CALCUL
Différentes configurations de calcul découlent des hypothèses CAS C
précédentes. Ces 4 configurations sont examinées en détail Milieu continu drainé sans écoulement
dans le chapitre 4.
Dans le cas de terrains plus perméables ou d’un creusement plus
lent (avec les valeurs d’exemple retenues ci-dessus, k > 10–6 m/s
CAS A et vavct < 0,1 m/h), le terrain est sollicité dans un état « drainé »
Milieu discontinu (c’est-à-dire sans surpressions interstitielles). La stabilité du
front de taille est alors conditionnée par les contraintes effectives
Dans les massifs rocheux peu fracturés, le réseau des discontinuités
exercées dans le terrain.
découpe des dièdres potentiellement instables. Des chutes de
blocs ou des glissements peuvent avoir lieu au front de taille Il n’y a pas d’écoulement de l’eau interstitielle dans le terrain
et/ou dans la zone du décousu. La géométrie des dièdres risquant (les pressions d’eau restant égales aux pressions hydrostatiques
de se décrocher dépend de la géométrie du front de taille et de initiales) dans les cas suivants :
l’orientation des familles principales de discontinuités. La stabilité de
• creusement hors nappe, naturellement ou induit par les travaux
ces dièdres en l’absence de renforcement dépend de leur géométrie,
(rabattement préalable de la nappe, drainage à l’avancement
du poids volumique de la matrice rocheuse et de la résistance au
conduisant à l’absence de pressions hydrostatiques jusqu’à
cisaillement mobilisable sur les surfaces de discontinuités.
au moins un diamètre à l’avant du front de taille) ;
• creusement au tunnelier pressurisé sous nappe dans des
Dans les massifs rocheux globalement peu fracturés, des zones
sols relativement grossiers, en dehors de conditions non
plus fracturées ou des zones de failles présentant des matériaux
maîtrisées correspondant au cas « drainé avec écoulement ».
dégradés, altérés peuvent exister, dans lesquels il peut être
pertinent de retenir des mécanismes de type B/C/D ci-dessous. La pression frontale à exercer est alors la somme de la pression
exercée par le terrain (calculée en contraintes effectives, avec
l’hypothèse d’une isotropie de résistance au cisaillement) et
CAS B de la pression (hydrostatique) de l’eau, soit :
Milieu continu non drainé
Dans le cas d’un creusement au tunnelier ou en méthode conven-
tionnelle à une vitesse d’avancement importante dans des terrains
de faible perméabilité, il peut être raisonnable de considérer que [2]
le creusement se fait en conditions « non drainées ».
Hors nappe, les contraintes effectives sont égales aux contraintes
Dans ce contexte, le calcul est effectué en contraintes totales totales et la solution se réduit sous la forme :
avec le poids volumique total γ des terrains et la résistance au
cisaillement non drainée des terrains (paramétrée par une
« cohésion non drainée » Cu et un angle de frottement φu nul).
La solution du problème (pression frontale minimale empêchant la [3]
rupture par effondrement σT-eff) s’écrit sous la forme suivante :
CAS D
Milieu continu drainé avec écoulement
[1] En méthode conventionnelle (ou éventuellement au tunnelier
fonctionnant en mode ouvert) sous nappe, le déséquilibre
Le choix entre conditions « non drainées » et « drainées » doit hydraulique entre l’intérieur du tunnel (pression nulle) et
résulter d’une réflexion spécifique à chaque contexte géotechnique le terrain (pression interstitielle initiale égale à la pression
et à chaque phasage réel de construction. Sur la base de calculs hydrostatique) conduit à un écoulement vers le front de taille.
hydro-mécaniques, Anagnostou & Kovari (1996) donnent
par exemple la limite suivante aux conditions non drainées : Au tunnelier à front pressurisé sous nappe, l’écoulement
perméabilité k inférieure à 10–7 à 10–6 m/s et vitesse vers le front est empêché : (a) grâce à la faible perméabilité
d’avancement du front vavct supérieure à 0,1 à 1 m/h. Attention du cake, le contrôle du débit de marinage et le contrôle de la
par exemple aux périodes d’arrêts de l’avancement sur chantier pression de bulle d’air en calotte de chambre au tunnelier à
(opérations hyperbares, week-ends, incidents de chantier) qui pression de boue, (b) grâce au bon dimensionnement de la vis
pénalisent l’avancement et peuvent conduire à rendre l’hypothèse d’extraction (longueur et diamètre) couplé à l’ajout d’éventuel
« non drainée » irréaliste. En cas de doute quant à la réalité adjuvants réduisant la perméabilité dans la chambre d’abattage
de conditions non drainées, les calculs drainés et non drainés au tunnelier à pression de terre. Des surpressions interstitielles
peuvent être effectués en parallèle afin de juger de la sensibilité peuvent toutefois être générées au front si ces conditions sont
des conditions de stabilité du front à l’égard de cette hypothèse. insuffisamment maîtrisées.

9
2
Dans l’ensemble de ces conditions, la distribution des pressions à partir de laquelle le front devient instable peut être évaluée.
interstitielles (donc également des contraintes effectives) dépend La forme générale de la pression minimale empêchant la rupture
des conditions aux limites hydrauliques du modèle (fonction par effondrement du front est donc la suivante :
de Hw) et de la perméabilité des terrains (k). La pression à exercer
au front de taille pour assurer sa stabilité dépend alors des
forces d’écoulement. En régime permanent, celles-ci dépendent
seulement du gradient hydraulique (fonction de Hw et k). En
régime transitoire, le temps (t) intervient également, et la durée [4]

2.4 APPROCHE DE LA SÉCURITÉ

Trois grandes approches sont possibles dans l’absolu : En terrains meubles, la recommandation allemande (DAUB, 2016)
qui traite plus particulièrement des tunneliers à front pressurisé,
• avec des coefficients de sécurité partiels sur les effets préconise de pondérer les effets des actions, avec des
des actions : dans ce cas, les calculs sont réalisés sans coefficients partiels de sécurité de 1,5 sur la pression σT-eff
coefficients de sécurité, puis les efforts obtenus liés à la
et de 1,05 sur la pression d’eau. Cette démarche peut être
stabilité des terres (ou la contrainte σT-eff si le front est
pertinente lorsque σT-eff >> 0 mais reste discutable lorsque le
homogène et totalement sous nappe) et à la pression de
front est en limite de stabilité, la condition σT-eff = 0 restant vraie
l’eau sont multipliés par des coefficients de sécurité ;
quel que soit le coefficient de sécurité considéré.
• avec des coefficients de sécurité partiels appliqués avant
calcul, directement sur les paramètres de résistance au
Par analogie avec les problématiques de stabilité de pente,
cisaillement du terrain ;
un facteur global de sécurité peut être estimé en rapportant
• en recherchant un coefficient de sécurité global, grâce au
les efforts moteurs aux efforts résistants. Une approche de
calcul du rapport entre les efforts moteurs et les efforts
calcul avec coefficients partiels sur la résistance du terrain peut
résistants sans pondération des paramètres de terrains
également être considérée en terrains meubles (Paternesi
ni des effets des actions.
et al., 2017) puisque la dégradation des paramètres du terrain
En tant « qu’ouvrages géotechniques de catégorie 3 » (§2.1 (en particulier l’angle de frottement) majore la géométrie du
de l’Eurocode 7-1), les tunnels sortent actuellement du champ mécanisme mais ne change pas le comportement d’ensemble
d’application de l’Eurocode 7, toutefois les approches utilisées qui reste de type milieu continu. C’est la démarche retenue dans
en pratique sont compatibles avec ses principes. Pour le dimen- ce document pour les terrains meubles avec les coefficients
sionnement des structures de soutènement et de revêtement des de sécurité partiels suivants :
tunnels par des approches « en déplacements » (type éléments
• 1,25 sur c’ et tan(φ’) en conditions drainées ;
finis), l’AFTES a pris position en recommandant la pondération
• 1,4 sur Cu en conditions non drainées.
des effets des actions : efforts normaux et moments de flexion
dans les structures (AFTES GT29, 2007). Ce choix est motivé
Dans le cas du creusement au tunnelier pressurisé, les pressions
par la volonté de ne pas induire de comportement plastique par
minimales σT calculées peuvent être majorées de 10 kPa dans
la prise en compte de paramètres dégradés pour le terrain qui
serait susceptible de perturber artificiellement la distribution des le cas d’un tunnelier à pression de boue (SS) et 30 kPa dans le
efforts entre le terrain et le revêtement. cas d’un tunnelier à pression de terre (EPB) afin de prendre en
compte les difficultés de maintien et de contrôle de la pression
Pour l’étude de la stabilité du front de taille en massif rocheux peu frontale appliquée (DAUB, 2016). Il est de même raisonnable
fracturé (chutes et glissements de dièdres) par des méthodes de majorer de 5 % les pressions d’eau dans le cas d’un milieu
« à la rupture » (type équilibre limite, calcul à la rupture), une continu drainé sans écoulement (cas C, §.2.3).
approche de type « facteur global de sécurité » est recommandée.
En effet, une dégradation des paramètres du massif rocheux Attention toutefois, l’approche préconisée consistant à appliquer
(niveau de fracturation, état d’altération, résistance au cisaillement une pondération partielle sur les paramètres de résistance au
sur les discontinuités, etc.) pourrait conduire à une modification cisaillement en terrains meubles est clairement proscrite pour le
significative de la géométrie des blocs, voire à une modification dimensionnement des structures de soutènement et revêtement,
radicale des mécanismes en jeu, le terrain pouvant passer d’un en raison de la modification engendrée sur les phénomènes
comportement de type milieu discontinu à celui d’un milieu continu. d’interaction sol / structure.
Par analogie avec le boulonnage radial en section courante
(AFTES GT30, 2021), un facteur global de sécurité de l’ordre
de 1,6 à 1,8 peut par exemple être recherché au front de taille.

10
2.5 APPROCHES DE MODÉLISATION
L’ensemble des modèles existants pour l’étude de la stabilité du ainsi que les procédures de mise en place des matériaux et leur
front de taille peut être regroupé en quatre grandes catégories caractérisation. La validation des modèles de calcul, idéalement
dont les quelques lignes ci-dessous donnent les grands principes, par le retour d’expérience d’instabilités sur des ouvrages réels en
avec quelques exemples en Illustration 3. vraie grandeur, souffrent d’observations bien souvent lacunaires
en particulier sur les mécanismes de déformation et de rupture
ainsi que sur les propriétés des matériaux. Quelques informations
Les modèles physiques et expérimentations (Exp)
issues d’instabilités de front observées in situ lors de travaux de
Plusieurs campagnes d’essais ont été réalisées sur des modèles creusement en méthode conventionnelle et au tunnelier sont
réduits centrifugés ou non (en Angleterre puis en France, plus reportées en annexe.
récemment en Asie), dans des sols purement cohérents, purement
frottants, puis cohérents frottants et stratifiés. La rupture du front Les modèles développés dans le cadre
de taille est généralement provoquée par le déplacement d’une
paroi rigide ou le dégonflement d’une membrane souple au
d’une approche de type équilibre limite (EL)
front de taille (Illustration 3a). Cette approche est la plus couramment utilisée dans la pratique
courante : elle consiste en un bilan direct des efforts en jeu
Ces modèles sont précieux puisqu’ils permettent, d’une part, (Illustration 3b). Sa mise en œuvre dans les massifs de sols,
d’identifier la cinématique des mécanismes de rupture à considérer ou de matériaux que l’on peut considérer comme continus à
dans les modèles semi-empiriques et analytiques (Equilibre Limite l’échelle de l’ouvrage, nécessite des hypothèses fortes quant à :
(EL) et Calcul à la Rupture (CR) décrits ci-après), et d’autre part,
• la géométrie du mécanisme : le front étant supposé
d’évaluer les pressions frontales à la rupture dans des conditions
rectangulaire, le mécanisme le plus classique (Horn, 1961)
contrôlées sans toutefois couvrir l’ensemble de la gamme de
est constitué d’un coin, délimité par un plan incliné à l’avant
variation des différents paramètres.
du front et deux faces triangulaires latérales en retour, ainsi
Le respect des conditions de similitude et la minimisation des que d’un prisme droit vertical à base rectangulaire remontant
effets d’échelle rendent complexes les dispositifs expérimentaux en surface (extension 3D du mécanisme de Terzaghi) ;

(a) Exemple de modèle réduit (b) Bilan des efforts sur le coin, mobilisé à l’avant du front
sur matériau pulvérulent en centrifugeuse dans les approches de type équilibre limite
(Chambon & Corté, 1994) (d’après Anagnostou & Kovari, 1996)

(c) Exemple de mécanisme 3D considéré dans l’approche (d) Exemple de modèle numérique 3D en éléments finis
cinématique du calcul à la rupture (Wong & Subrin, 2006) (Vermeer et al., 2002)

Illustration 3 : Quelques représentations de modèles consacrés à l’étude de la stabilité du front de taille.

11
2
• le profil des contraintes appliquées sur chaque face du de terrains hétérogènes puisque les champs de contraintes et
mécanisme (ou leur résultante globale) : la résistance au de vitesse requièrent la vérification des conditions de continuité
cisaillement est définie par le critère de résistance de entre matériaux de caractéristiques différentes.
Coulomb (éventuellement Tresca) et les contraintes
normales sur les surfaces de rupture présupposent Cette approche, quoique peu pratiquée, présente le mérite de
une valeur du rapport entre les contraintes horizontales proposer des fuseaux de valeurs possibles et donc d’encadrer
et verticales dans le massif. les conditions de rupture selon une approche rigoureuse et
donc de vérifier la pertinence des approches par équilibre
Le principe consiste à comparer les efforts moteurs (poids du limite reposant sur des hypothèses plus nombreuses. Les
coin et contrainte verticale s’exerçant sur sa face supérieure, mécanismes de rupture virtuels construits « à la main » dans
éventuelle surcharge en surface, forces d’écoulement) aux le cadre de l’approche cinématique peuvent être affinées par
efforts résistants mobilisables (résistance au cisaillement du des approches numériques.
terrain sur les différentes faces du coin et pression frontale
simulant l’effet d’un confinement).
Les modèles numériques (Num)
Il est possible de prendre en compte des terrains stratifiés hori-
Les modèles numériques (Illustration 3d) considérés consistent
zontalement (par exemple en contexte de bassin sédimentaire)
à résoudre le système aux dérivées partielles issu des équations
et de considérer l’effet d’une nappe sans changer le mécanisme
d’équilibre, de la loi de comportement du terrain, des conditions
de rupture considéré.
aux limites et des conditions intiales, par l’écriture du problème
sous forme discrète (en éléments finis ou différences finies
Les modèles développés dans le cadre explicites). La précision de la solution ainsi approchée est
de la théorie du calcul à la rupture (CR) d’autant plus grande que le maillage est fin. Compte-tenu du
caractère tridimensionnel des mécanismes en jeu au niveau du
Les conditions de stabilité du front peuvent être estimées par front de taille, les temps de calcul nécessaires pour obtenir une
la théorie du calcul à la rupture (Salençon, 1990), dont les précision acceptable peuvent être significatifs. Une attention
résultats tendent vers ceux de l’analyse limite considérant particulière doit être portée à la représentativité des hypothèses
l’hypothèse d’un comportement rigide parfaitement plastique retenues : blocage radial des nœuds en périphérie du front
associé du terrain.
de taille, gradient vertical de pression frontale (en particulier
dans le cas des tunnels de grand diamètre), caractéristiques
Basée sur la détermination de champs de contraintes statiquement
mécaniques des terrains, etc.
admissibles, l'approche statique, dite par l’intérieur, conduit à
une estimation par défaut, donc dans le sens de la sécurité, du
Ce type d’approche en déformation est utilisé pour approcher les
domaine des chargements supportables. En d’autres termes,
conditions de rupture, soit en simulant la réduction progressive
s’il est possible de trouver un champ de contraintes dans le terrain
de la pression au front, soit en utilisant une méthode consistant
vérifiant les équations d’équilibre, les conditions aux limites en
à réduire les paramètres de résistance. Par rapport à l’approche
contraintes et le critère de résistance, alors la stabilité du front de
par l’extérieur du calcul à la rupture, la recherche de la cinéma-
taille est assurée. Dans le cas d’une rupture par effondrement,
tique de déformation et de rupture se fait parmi une « infinité »
la pression frontale étant un paramètre de chargement qui
de géométries possibles, avec optimisation résultant du calcul
s’oppose à l’instabilité, cette approche conduit à un majorant
numérique sans figer a priori la géométrie de la rupture.
de la pression frontale de rupture.

Partant d’un état initial de contraintes géostatiques dans le

A contrario, l'approche cinématique, dite par l’extérieur, est
dérivée de la proposition de mécanismes de rupture virtuels terrain, la rupture par effondrement du front de taille peut être
cinématiquement admissibles (Illustration 3c), c’est-à-dire recherchée en diminuant progressivement la pression frontale
vérifiant les conditions aux limites en vitesses. Pour une famille simulant la procédure d’essai habituellement pratiquée dans
de mécanisme caractérisée par un ensemble de paramètres les approches expérimentales en laboratoire. Les isovaleurs
géométriques, la puissance des efforts extérieurs est comparée des déformations plastiques de cisaillement permettent alors
à la puissance de dissipation que le terrain peut mobiliser d’exhiber la géométrie du mécanisme de rupture et notamment
compte tenu de ses capacités de résistance. L'optimisation par de mettre en évidence les surfaces de rupture dans le cadre
rapport à la géométrie du mécanisme conduit à un minorant de mécanismes par blocs rigides. La pression frontale limite
de la pression frontale de rupture : pour la valeur calculée, la de rupture par effondrement peut également être estimée en
rupture par effondrement du front est certaine. représentant l’évolution de l’extrusion de certains points du front
de taille en fonction de la réduction de la pression frontale.
Par rapport aux approches de type équilibre limite, la géométrie
optimale du mécanisme de rupture résulte du calcul par optimi- Si le front est auto-stable (c’est-à-dire pour une pression au
sation des paramètres géométriques de la famille de mécanisme front nulle), la marge de sécurité associée peut être estimée
considérée. Le résultat ne dépend pas des conditions initiales en diminuant progressivement la résistance au cisaillement du
concernant l’état de contraintes. Les modèles développés terrain jusqu’à la rupture (méthode dite « c-phi réduction ») selon
en milieu homogène ne s’appliquent pas directement en cas une approche de type pondération des paramètres du terrain.

12
3 ÉVALUATION
DES CONDITIONS DE STABILITÉ
On présente, successivement pour les quatre configurations de sous le mot clé « Tunnel face stability » dans les revues
calcul (A, B, C, D) identifiées au §.2.3, les différents modèles scientifiques internationales « Tunnelling and Underground Space
existants suivant les quatre approches (Exp, EL, CR, Num) Technology » et « Computers and Geotechnics », dont 60 sur
explicitées au §.2.5. la période 2018-2020 (3 % des articles publiés dans ces revues).
De très nombreux articles scientifiques ont été consacrés à Les références citées ci-après ne sont pas exhaustives et ne
l’étude de la stabilité du front de taille depuis une quarantaine retracent pas totalement l’historique des travaux ayant conduit à
d’années. En guise d’illustration, citons que, à la date de la construction de l’ensemble de ces modèles : tous les modèles
rédaction de cette note, près de 100 articles sont référencés « intermédiaires » n’ont pas été rappelés.

3.1 MILIEU DISCONTINU

Dans les massifs rocheux modérément fracturés, assimilables à
des milieux discontinus à l’échelle de l’ouvrage, les instabilités au
front de taille se manifestent par des chutes ou des glissements
de dièdres isolés par le réseau de discontinuités (Illustration 4).

La géométrie des mécanismes de rupture est dictée par

l’orientation et le pendage des discontinuités par rapport à l’axe
du tunnel, leur extension, leur espacement ainsi que la forme
du front de taille. La stabilité naturelle de ces dièdres dépend du
poids volumique de la roche et de la résistance au cisaillement
mobilisable sur les discontinuités (fonction de la rugosité des
épontes, de leur ouverture, d’un éventuel remplissage, de la
présence d’eau, etc.). Le phénomène de dilatance empêchée
est un facteur essentiel dans la stabilité naturelle des excavations
souterraines en milieu rocheux.
Illustration 4 : Exemple d’instabilité en massif rocheux dans le tunnel
À ce titre, la démarche de dimensionnement peut être résumée
de Chavannes (France, 70).
sous la forme des étapes suivantes :
La stabilité des blocs de plus petits volumes (stabilité locale)
• Identification des données géométriques et géo-
pourra être assurée grâce à la résistance au cisaillement du
mécaniques des discontinuités :
béton projeté, ou en dégageant ces petits blocs instables.
Les principales familles de discontinuités, leur orientation
• Calcul du coefficient de sécurité de chaque bloc et
et leur espacement moyen sont évalués grâce à l’analyse calcul du boulonnage nécessaire pour obtenir un
structurale du massif déduite d’informations issues princi- coefficient de sécurité acceptable :
palement de levés de terrain et de sondages carottés. La Le calcul du coefficient de sécurité est obtenu, dans une
résistance au cisaillement des discontinuités peut être démarche d’équilibre limite, en faisant le rapport entre les
mesurée à partir d’essais de cisaillement sur discontinuités efforts moteurs (poids du bloc) et résistants (résistance
en laboratoire ou bien évaluée à partir des travaux de au cisaillement le long des discontinuités et apport du
Barton (1977) sur l’état de surface (JRC) et la résistance boulonnage). Les boulons ont deux effets : d’une part, ils
des épontes (JCS). En l’absence de boulons ou avec reprennent directement une partie du poids du dièdre, et
des boulons non précontraints, il est recommandé par d’autre part, ils augmentent la résistance au cisaillement
AFTES GT30 (2021) de ne pas retenir de cohésion sur mobilisable le long des faces du dièdre, en augmentant
les discontinuités (résistance résiduelle). la contrainte normale s’ils sont précontraints. Ces deux
• Identification des blocs potentiellement instables : effets sont bien sûr à intégrer dans les calculs de stabilité.
La géométrie des blocs les plus volumineux susceptibles de Un coefficient de sécurité de l’ordre de 1,6 à 1,8 peut être
glisser au front de taille est évaluée par une approche basée sur recherché, par analogie avec les valeurs préconisées par
les projections stéréographiques ou par modélisation numéri- AFTES GT30 (2021) pour un boulonnage radial en section
que 3D. Un schéma de boulonnage permettant d’assurer courante. Pour des situations où l’excavation reste ouverte
la stabilité générale peut en être déduit (cf. point suivant). longtemps, ce coefficient peut monter jusqu’à 2,0.

13
3

3.2 MILIEU CONTINU NON DRAINÉ

Lors du creusement à une vitesse importante dans des terrains modélisant le front de taille. L’ensemble est disposé sur un bâti
peu perméables homogènes isotropes, une hypothèse de compor- d’essai triaxial. Dans un premier temps, une contrainte verticale σ0
tement non drainé des terrains peut être légitime quelle que soit est imposée à l’échantillon, puis la contrainte normale sur le
la méthode de creusement (méthode conventionnelle ou tunnelier), trou p0 est réduite progressivement jusqu’à la rupture du front
même si les conditions nécessaires (cf. §2.3) sont rarement réunies de taille (Illustration 6).
en pratique. Dans ce contexte, le critère de résistance des terrains
peut être considéré de type « Tresca » (sols purement cohérents) Ces analyses conduisent les auteurs à écrire la condition de
caractérisé par la résistance au cisaillement non drainée (ou stabilité du front de taille dans les terrains purement cohérents
par abus de langage « cohésion non drainée ») Cu. L’Illustration 5 sous la forme de l’équation [5], où N couramment appelé
synthétise les principales références consacrées à ce sujet. « nombre de stabilité » exprime en fait un facteur de charge par
rapport à la résistance du terrain.
Broms & Bennermark (1967) sont les premiers auteurs connus
à s’être intéressés à cette problématique. Ils réalisent des essais
en laboratoire à l’aide d’un échantillon cylindrique d’argile contenu
dans un cylindre en acier, comportant un trou de 2 cm de diamètre [5]

Cas Type Référence Description succincte

Premiers essais en laboratoire réalisés en Suède (dispositif expérimental s’apparentant à une
Broms & Bennermark (1967) cellule triaxiale avec un trou permettant l’extrusion du terrain) ; Analyse d’instabilités de front sur
14 chantiers ; Introduction du « nombre de stabilité ».
Essais réalisés en centrifugeuse à Cambridge (Angleterre), avec un tunnel modélisé sous la forme d’un
Milieu continu homogène non drainé

Exp
Kimura & Mair (1981) tube de 60 mm pressurisé. Différentes longueurs de décousu (nul à « infini », i.e. déformations planes)
et épaisseur de couvertures étudiées.
Essais réalisés en centrifugeuse à Delft (Pays-Bas), avec un tunnel modélisé sous la forme d’un tube
Bezuijen & Van Seters (2006)
de 150 mm pressurisé.
Adaptation du modèle d’Anagnostou (2012) réalisé en terrain cohérent-frottant, au cas purement cohérent.
Perazelli & Anagnostou (2017) La nouveauté réside dans le calcul de la force verticale V s’exerçant sur le dièdre, déduite des approches
statiques du calcul à la rupture de Gunn (1980), au lieu d’une voûte de décharge de type « Terzaghi ».
EL
Extension du modèle de Perazzelli & Anagnostou (2017) en considérant une géométrie circulaire
Champagne de Labriolle (2018) du front de taille (au lieu de rectangulaire) et une distribution plus réaliste des contraintes
horizontales le long des surfaces de rupture.
Approche cinématique : plusieurs cas étudiés dont mécanisme 3D à deux blocs (cylindres elliptiques).
Davis et al. (1980) Approche statique : champs de contraintes sphériques et cylindriques + comparaison avec résultats
CR
expérimentaux dans le cas d’un mécanisme de rupture en déformation plane en coupe transversale.
Mollon et al. (2013) Approche cinématique : mécanisme 3D ayant une forme de tore.
Num Ukrichton et al. (2017) Calculs aux éléments finis 3D avec prise en compte d’une augmentation linéaire avec la profondeur de Cu.

Illustration 5 : Synthèse des approches existantes en conditions « non drainées ».

(a) (b)

Illustration 6 : Essais en laboratoire de Broms & Bennermark (1967) : (a) dispositif expérimental, (b) synthèse des résultats.

14
Broms & Bennermark (1967) ont également réalisé une synthèse Bezuijen & Van Seters (2006) ont également effectué deux essais
des conditions de stabilité observées au front lors de onze fonçages en centrifugeuse à l’université de Delft, pour des faibles épaisseurs
de conduites d’assainissement (diamètres compris entre 0,8 et de couverture (C/D égal à 0,6 et 0,8). Les auteurs ne décrivent
2,3 m environ) réalisés en Suède, ainsi que trois tunnels d’environ pas la géométrie du mécanisme de rupture observé, mais donnent
7 m de diamètre en Norvège et aux États-Unis, pour des épaisseurs les valeurs des nombres de stabilité associées à la rupture
de couverture comprises entre 1,5.D et 4.D environ. La plupart (cf. Illustration 8).
présente des nombres de stabilité inférieurs à 6, sans probléma- Différentes approches théoriques ont été développées en parallèle
tiques fortes d’instabilité du front. L’expérience du tunnel de Tyholt par différents auteurs. Dans le domaine du calcul à la rupture,
(Norvège) en 1950 est en revanche intéressante avec deux citons notamment les travaux de Davis et al. (1980) et Mollon et al.
instabilités de front observées, pour des rapports C/D (couverture / (2013), concernant l’approche cinématique (approche par
l’extérieur). Davis et al. (1980) ont considéré différents mécanismes
diamètre) égaux respectivement à 2,0 et 2,2, et des nombres de
en section courante du tunnel, mais également un mécanisme au
stabilité de l’ordre de 6,2 et 7,8. La longueur du décousu lors des
front de taille (décousu nul), constitué de deux cylindres elliptiques.
effondrements n’est pas connue. Une pression d’air comprimée, Mollon et al. (2013) l’ont amélioré en considérant un mécanisme
imposée ensuite dans le tunnel, réduit les nombres de stabilité de rupture en forme de tore. Concernant l’approche statique,
de 2 unités permettant ainsi d’assurer la stabilité. Davis et al. (2013) ont étudié différents champs de contraintes,
de géométrie cylindrique ou sphérique au front de taille.
Des essais en centrifugeuse sur un modèle réduit de tunnel ont
ensuite été réalisés à Cambridge sur des argiles (Schofield, 1980 ; Parmi les approches d’équilibre limite, Perazzelli & Anagnostou
(2017) ont étendu l’approche développée dans le cas des terrains
Kimura & Mair, 1981). Le tunnel est modélisé sous la forme d’un
cohérents-frottants (cf. §3.3) au cas de terrains purement cohérents.
demi-tube pressurisé de 60 mm de diamètre, porté à une accélé-
La principale modification apportée concerne le calcul de la
ration de l’ordre de 120.g simulant un tunnel réel de l’ordre de 7 m de contrainte verticale s’exerçant sur la face supérieure du coin
diamètre. Différentes épaisseurs de couverture (1,2 < C/D < 3,3) sont considéré au front de taille. Contrairement au cas cohérent-frottant où
étudiées, ainsi que différentes longueurs de décousu d, comprise celle-ci est déduite d’un mécanisme de Terzaghi, basé sur l’hypothèse
entre 0 (déplacements autorisés uniquement au front de taille), de surfaces de glissement verticales, celle-ci est ici déduite de
et l’infini (étude en section courante d’un tunnel non soutenu). l’approche statique du calcul à la rupture conformément aux travaux
de Gunn (1980). Trois formes de champs de contraintes dans un
Ces essais mettent bien en évidence que : terrain situé au-dessus d’une trappe rectangulaire sont proposés
• le mécanisme d’effondrement (Illustration 7a) affecte la avec : un champ de contrainte de géométrie cylindrique dont l’axe
est parallèle à la longueur de la trappe (longitudinalement au tunnel
totalité du front de taille ;
dans notre cas), un champ de contrainte de géométrie cylindrique
• la surface de glissement inférieure est inclinée à moins de 45° dont l’axe est parallèle à la largeur de la trappe (c’est-à-dire transver-
par rapport à l’horizontale, ce qui conduit à un mécanisme salement au tunnel dans notre cas), et un champ de contrainte de
s’étendant assez loin à l’avant du front de taille (quasiment géométrie sphérique englobant la trappe. La valeur minimale de
un diamètre à hauteur d’axe du tunnel) ; la contrainte verticale issue de ces trois champs est retenue.
• le mécanisme initial se propage très rapidement jusqu’à la
Champagne de Labriolle (2018) a amélioré l’approche de Perazzelli
surface en cas de couverture faible (1 à 2.D) avec un fontis & Anagnostou (2017) en considérant la géométrie cylindrique du
final s’avançant jusqu’à environ 2.D à l’avant du front, et 0,5.D tunnel. L’intégration des efforts considère en effet des tranches
à l’arrière ; de hauteurs infinitésimales et de largeur variable sur la hauteur
• il n’apparaît pas clairement de surface de cisaillement à du front, contrairement à l’approche de Perazzelli & Anagnoustou
l’intérieur de ce mécanisme. (2017) où le front est considéré rectangulaire. Enfin, Ukrichton
et al. (2017) ont effectué des calculs numériques en éléments finis
Les modélisations numériques réalisées plus tard par Ukrichton et en considérant un tunnel circulaire et un décousu nul. Perazzelli &
al. (2017) confortent également ces observations expérimentales Anagnostou (2017) ont également effectué deux calculs numériques
(Illustration 7b). en différences finies avec un tunnel de section carrée.

(a) (b)

Illustration 7 : Géométrie du mécanisme : (a) Essais en centrifugeuse de Kimura & Mair (1981), (b) Modélisation numérique de Ukrichton et al. (2017).

15
3
L’Illustration 8 compare les nombres de stabilité N obtenus par Notons également les limites suivantes des approches précé-
ces différentes approches dans le cas d’un décousu nul et d’une demment citées :
« cohésion non drainée » constante. Cette synthèse montre tout • les approches statiques du calcul à la rupture considèrent
d’abord que les résultats expérimentaux, notamment Kimura & des champs de contraintes simples, dont il est difficile
Mair (1981), sont encadrés par : d’imaginer qu’ils soient pleinement représentatifs de l’état
de contrainte régnant autour du front de taille ;
• les approches par l’extérieur, pour lesquelles la rupture est
• les modèles d’équilibre limite sont basés sur un coin, délimité
certaine, qu’elles soient de type « approche cinématique
au niveau du front de taille par une surface de glissement
du calcul à la rupture » (Davis et al. (1980), amélioré par
plane, contraire aux observations expérimentales qui
Mollon et al. (2013)), ou numérique (Ukrichton et al. (2017)
tendent à considérer un arc de spirale ;
et Perazzelli & Anagnostou (2017)). Noter que l’écart entre
• il est surprenant qu’un écart aussi important existe entre
les deux approches numériques précédentes peut être lié
les approches numériques et les approches expérimentales.
à la géométrie considérée (front carré ou circulaire), mais
également au maillage assez grossier considéré par Compte-tenu des développements précédents, la pression
Perazzelli & Anagnostou (2017) ; frontale limite de rupture par effondrement peut ainsi être
• les approches par l’intérieur, pour lesquelles le front est calculée à partir de l’équation [7], en intégrant la réflexion du
potentiellement stable, qu’elles soient de type « approche §2.4 sur la prise en compte de la sécurité. Le choix de la valeur
statique du calcul à la rupture » (Davis et al., 1980) ou du nombre de stabilité critique Ncrit est à effectuer par l’ingénieur
« équilibre limite du dièdre en tenant compte d’un effort sur la base de l’Illustration 8, à l’intérieur du fuseau illustré. Les u
vertical sur le dièdre déduit de l’approche statique du valeurs basses du fuseau conduisent à un dimensionnement + +
calcul à la rupture » (Perazzelli & Anagnoustou (2017), plus sécuritaire.
Champagne de Labriolle (2018)). u
+ + , ,
L’écart entre ces deux grandes catégories d’approches reste
assez important, de l’ordre de 30 %. À titre d’exemple, pour u u u
++ + + sur l’Illustration 8,
+ + choisi
un tunnel de 10 m de diamètre avec une couverture de 15 m
d’épaisseur, creusé dans un terrain présentant 30 kPa de u

« cohésion non drainée », le nombre de stabilité est compris

,
entre 6,5 et 9, ce qui conduit à une pression limite au front, sans u u u
sinon [7]
coefficient de sécurité, comprise entre 130 et 200 kPa.

INSTABILITÉ

Broms & Bennermark (1967) -

chantier du tunnel de Tyholt (Norvège, 1950)
Kimura & Mair (1981) -
Expérience sur tunnel centrifugé (cas P/D=0)
Bezuijen & Van Seters (2006) -
Expérience sur tunnel centrifugé

STABILITÉ
Davis et al. (1980) - Approche cinématique
(mécanisme à deux cylindres elliptiques)
Broms & Bennermark (1967) - Mollon et al. (2013) - Approche cinématique
chantier du tunnel de Tyholt (Norvège, 1950) (mécanisme tore)
Kimura & Mair (1981) - Ukrichton et al. (2017) -
Expérience sur tunnel centrifugé (cas P/D=0) Calcul numérique (front cylindrique)
Bezuijen & Van Seters (2006) - Perazzelli & Anagnostou (2017) -
Expérience sur tunnel centrifugé Calculs numériques (front carré)
Broms & Bennermark (1967) - Davis et al. (1980) - Approche cinématique Davis et al. (1980) - Approche statique
chantier du tunnel de Tyholt (Norvège, 1950) (mécanisme à deux cylindres elliptiques) (max. des champs cylindrique et sphérique)
Kimura & Mair (1981) - Mollon et al. (2013) - Approche cinématique Perazzelli & Anagnostou (2017) -
Expérience sur tunnel centrifugé (cas P/D=0) (mécanisme tore) Équilibre limite - Force V déduite
Bezuijen & Van Seters (2006) - Ukrichton et al. (2017) - de l’approche statique de Gunn (1980)
Expérience sur tunnel centrifugé Calcul numérique (front cylindrique) Champagne de Labriolle (2018) -
Davis et al. (1980) - Approche cinématique Perazzelli & Anagnostou (2017) - Approche similaire à P & A (2017),
(mécanisme à deux cylindres elliptiques) Calculs numériques (front carré) mais avec « front cylindrique »
Mollon et al. (2013) - Approche cinématique Davis et al. (1980) - Approche statique
Illustration 8 : Comparaison
(mécanisme tore) des nombres de stabilité déduits
(max. desdes différents
champs modèles
cylindrique existants.
et sphérique)
Ukrichton et al. (2017) - Perazzelli & Anagnostou (2017) -
Calcul numérique (front cylindrique) Équilibre limite - Force V déduite
Perazzelli & Anagnostou (2017) - de l’approche statique de Gunn (1980)

16 Calculs numériques (front carré) Champagne de Labriolle (2018) -

Davis et al. (1980) - Approche statique Approche similaire à P & A (2017),
mais avec « front cylindrique »
 L’influence de la longueur du décousu sur la stabilité du front de taille 0,2.D (cf. Figure 24 de l’auteur) dans le cas d’une « cohésion non
a également été étudiée par différents auteurs. Citons notamment drainée » supérieure ou égale à 50 kPa (soit Cu / (γ.D) > 0,3).
les essais en centrifugeuse de Kimura & Mair (1981), réalisés pour Pour des « cohésions non drainées » très faibles (Cu / (γ.D) < 0,3),
un décousu important, égal à 0,5.D, 1.D, 2.D et infini (c’est-à-dire l’influence du décousu dans l’approche de Champagne de Labriolle
équivalent à un calcul de stabilité en section courante). Les résultats (2018) est plus importante qu’indiqué au paragraphe précédent.
obtenus révèlent des nombres de stabilité significativement plus
faibles en « déformations planes » (décousu infini) que dans le cas Quelques auteurs, notamment Ukrichton et al. (2017) et
d’un décousu nul (N divisé par deux environ). Pour un décousu Champagne de Labriolle (2018) ont par ailleurs intégré une
limité à 0,2.D, on peut estimer que le nombre de stabilité est réduit évolution de la « cohésion non drainée » avec la profondeur.
d’environ 0,5 d’après les conclusions obtenues par les auteurs. Les résultats obtenus ne sont pas présentés ici dans la mesure
où il est rare, en pratique, de disposer de suffisamment de
Les calculs analytiques de Champagne de Labriolle (2018) con- résultats d’essais pour mettre en évidence cette évolution de la
duisent au même ordre de grandeur pour un décousu de longueur résistance au cisaillement avec la profondeur.

3.3 MILIEU CONTINU DRAINÉ SANS ÉCOULEMENT

3.3.1 Terrain homogène

Dans le cas de terrains perméables et d’un creusement lent
(à titre indicatif, comme précédemment, avec k > 10–6 m/s et
vavct < 0,1 m/h), le terrain est sollicité dans un état « drainé ».
La stabilité du front de taille est dans ce cas conditionnée par les
contraintes effectives exercées dans le terrain. On considère ici
qu’il n’y a pas d’écoulement de l’eau interstitielle dans le terrain
(les pressions d’eau restent égales aux pressions hydrostatiques
initiales). Cette condition est rencontrée si le creusement a
lieu hors nappe ou si un tunnelier à front pressurisé est utilisé
(hors conditions décrites ultérieurement au §3.4).

Dans cette configuration, les calculs sont effectués en contraintes

effectives avec :
• le poids volumique γ* des terrains égal à γ’ = γsat – γw (poids
volumique déjaugé) si le tunnel est sous la nappe, ou à γh
(poids volumique humide) si le tunnel est situé hors nappe ;
• la résistance au cisaillement drainée, exprimée sous la forme Illustration 9a : Géométrie du mécanisme de rupture en terrain cohérent-
d’une cohésion effective c’ et d’un angle de frottement interne frottant (coupes transversales et longitudinales) : essais de Chambon
φ’ dans le cadre d’un critère de Mohr-Coulomb. & Corte (1989).

Illustration 9b : Géométrie du mécanisme de rupture en terrain cohérent-frottant (coupes transversales et longitudinales) : essais de Berthoz et al. (2012).

17
3
Le cas des sols partiellement saturés est généralement traité • Chambon & Corté (1989) réalisés en centrifugeuse avec
en considérant une cohésion effective et un angle de frottement des sables secs, où le creusement est généré via le
« apparents » déduits d’essais en laboratoire réalisés sur le dégonflage d’une membrane souple au front de taille, avec
sol pris dans un état de saturation identique à celui observé un soutien radial constant grâce à un tube métallique ;
sur site. • Berthoz et al. (2012) réalisés sur un modèle réduit de tunnelier
à pression de terre non centrifugé, avec des sables secs ou
Différentes campagnes d’essais ont été réalisées en laboratoire légèrement humides, leur conférant une cohésion de quelques
dans cette configuration. Citons notamment les essais de : kPa correspondant à une dizaine de kPa en vraie grandeur.

Cas Type Référence Description succincte

Essais centrifugés de tubes rigides à fronts flexibles sur des terrains cohérents frottants
Chambon & Corté (1989) → Forme des mécanismes + valeur de pression frontale pour quelques valeurs de H/D.

Essais non centrifugés de tubes rigides à fronts déplaçables sur des sables secs
Exp Messerli et al. (2010) → Forme des mécanismes + valeurs de pression frontale limite pour quelques valeurs de H/D.

Essais sur modèle réduit non centrifugés en terrains purement frottant et cohérents frottant
Berthoz et al. (2012) → Forme des mécanismes, valeurs de pression frontale limite pour quelques valeurs de c’.

Modèles en éléments finis 3D avec un terrain de comportement élastique parfaitement plastique régi
par un critère de Mohr-Coulomb, et un déplacement autorisé uniquement au front (déplacements
Vermeer et al. (2002)
radiaux bloqués en paroi). Différentes valeurs d’angle de dilatance et de coefficient de pression des
Num terres au repos ont été testées, sans effet significatif sur la pression de rupture par effondrement.
Un modèle similaire à celui de Vermeer et al. (2002) a été utilisé, avec quelques études paramétriques
Alagha & Chapman (2019)
complémentaires.
Mécanisme type Horn (1961). Equilibre du coin au front avec ces hypothèses :
(1) front rectangulaire (B*H),
(2) Effort vertical sur le coin déduit d’une voûte de décharge avec σH = K1.σV0 où K1 = 0,8,
Milieu continu homogène drainé sans écoulement

Anagnostou & Kovari (1996)

(3) Effort de cisaillement le long du coin déduit en supposant que la contrainte verticale augmente
linéairement avec la profondeur sur la hauteur du coin, avec la contrainte horizontale proportionnelle
via un facteur K2 = 0,4.
Reprise du modèle d’Anagnostou & Kovari (1996) en :
EL (1) supprimant l’hypothèse (3) précédente et en la remplaçant par une analyse de type « voûte de
Anagnostou (2012)
décharge » dans le plan horizontal,
(2) choisissant σH =1,0 σV (au lieu de 0,8) pour le bloc (hypothèse (2)) et le coin.
Préconise d’utiliser la méthode d’Anagnostou & Kovari (1996) mais avec un coefficient K2 = (K0 + Ka)/2
DAUB (2016)
au lieu de 0,4.
Reprise du modèle d’Anagnostou & Kovari (1996) en : (i) exprimant K à partir de l’hypothèse que le
Champagne de Labriolle (2018) critère est vérifié sur les frontières du mécanisme (donc relation entre σV et σH , qui ne sont pas forcément
contraintes principales) ; (ii) considérant le caractère circulaire du front de taille (au lieu de rectangulaire).
Leca & Dormieux (1990) Approche cinématique : modèle initial : 1 ou deux blocs 2D.
Approche cinématique : modèle 3D à un bloc défini par des spirales logarithmiques sans surface de
Wong & Subrin (2006)
discontinuité interne.
Approche cinématique : modèle 3D avec un nombre important (mais variable) de surfaces de
Mollon et al. (2011)
discontinuités découpant des blocs définis par des arcs de spirales logarithmiques.
Approche cinématique : reprise du modèle de Mollon et al. (2011) en remplaçant la partie du
mécanisme située au-dessus de la clé du tunnel par une pression déduite d’un mécanisme de
Zou et al. (2019)
type « équilibre limite » avec une hypothèse de contrainte normale sur ce mécanisme (Ks) tiré de
diverses approches (analytique, numérique, expérimentale).
CR
Approche cinématique : reprise de la géométrie du modèle de Wong & Subrin (2006), en considérant
Li et al. (2019)
un sol partiellement saturé.
Approche cinématique : modèle 3D à deux blocs (un tétraèdre au front, et un prisme à base
Quarmout et al. (2019) triangulaire au-dessus). Résolution par la méthode KEM (Kinematic Element Method), qui n’est pas
exactement du Calcul à la Rupture, mais s’y apparente.
Approche cinématique : mécanisme 3D à un bloc défini par des spirales logarithmiques sans surface
Li & Yang (2020)
discontinuité interne, et avec troncature du critère de résistance en traction.
Approche cinématique : mécanisme 2D à 3 blocs, en considérant l’existence d’un décousu. Prise en
Senent et al. (2020)
compte possible également d’une voûte parapluie.

Illustration 10 : Synthèse des approches existantes en milieu continu homogène drainé sans écoulement.

18
Ces essais, dont les résultats sont présentés en Illustration 9, In fine, l’ensemble des approches citées ici convergent. Vermeer
révèlent que le mécanisme : et al. (2002) ont également effectué des calculs numériques en
considérant une longueur non soutenue d non nulle. Ils ont montré
• intercepte tout le front de taille ; que Nγ augmente assez significativement (c’est-à-dire que la
• est peu étendu à l’avant du front de taille (moins de 0,5.D) pression minimale nécessaire à exercer au front de taille et au droit
et remonte peu au-dessus de la clé (0,1.D environ pour du décousu augmente), lorsque la longueur non soutenue devient
un angle de frottement interne de l’ordre de 35°) ; supérieure à environ 0,3.D. Senent et al. (2020) l’ont également
• peut se propager par rupture progressive jusqu’à la surface montré via une approche analytique en calcul à la rupture (mécanisme
dans le cas d’une faible épaisseur de couverture. longitudinal 2D à 3 blocs définis par des spirales logarithmiques).

De nombreux modèles numériques et analytiques ont également Compte-tenu de la cohérence entre ces différents résultats,
été consacrés à cette configuration. L’Illustration 10 en donne l’expression [10], déduite de Vermeer et al. (2002) peut être
les principales références. retenue pour la plupart des cas rencontrés en pratique (d/D ≤ 0,5,
C/D ≥ 0,8 et φ’ ≥ 20°).
Les modélisations numériques réalisées par Vermeer et al. (2002)
ou Berthoz et al. (2012), confirment ces grands traits des
mécanismes observés en milieu cohérent-frottant, sur lesquels
se sont appuyés différents modèles analytiques récents (Wong
& Subrin (2006), Mollon et al. (2011), etc.). avec

La pression frontale minimale de rupture par effondrement associée

à ce type de mécanisme peut s’exprimer sous la forme suivante :
où : ,

avec , ,

[9] (sous nappe)

,
(hors nappe)
sinon
Différents auteurs ont exprimé la fonction Nγ de l’équation [9] via
différentes approches. L’Illustration 11 en synthétise les résultats. [10]

Illustration 11 : Comparaison entre les différents modèles existants.

19
3
Dans le cadre de l’approche cinématique du calcul à la rupture,
le mécanisme peut intercepter la surface du terrain naturel en
cas de tunnel à très faible profondeur dans des terrains peu
frottants (typiquement pour C/D ≤ 0.3 (cot φ’ – 1)) environ. La
pression minimale de stabilité du front est alors fonction de la
surcharge en surface σs et peut s’exprimer sous la forme de
l’équation [11] (Wong & Subrin 2006), analogue aux approches
classiques à la rupture pour le dimensionnement des fondations
superficielles ou profondes.

[11]

Les calculs numériques réalisés par Alagha & Chapman

(2019) montrent également que l’influence de l’épaisseur de
couverture (rapport C/D) devient significative lorsque l’angle
de frottement interne du terrain est très faible (φ’ < 20°). Cette Illustration 12: Influence de la résistance au traction sur les pressions
hypothèse de tunnel très peu profond et de terrains très peu limite de rupture par effondrement (Li & Yang, 2020). Ici, d est le diamètre
frottants est toutefois peu fréquente dans le cadre des projets. du tunnel (D) et ξ = 0 signifie que le critère n’est pas du tout tronqué.

L’ensemble des modèles numériques et analytiques présentés En terrain purement frottant (sables secs ou saturés), un
considère un critère de Mohr-Coulomb avec une résistance à mécanisme d’écoulement prenant naissance en clé et remontant
la traction non nulle des terrains. Li & Yang (2020) ont réalisé jusqu’en surface, à l’image d’un « sablier » (Illustration 13), peut être
des calculs analytiques (par l’approche cinématique du calcul à la mobilisé d’après les essais sur modèle réduit réalisés par Berthoz
rupture) et numériques pour cerner l’influence de cette hypothèse et al. (2012). En pratique, ce type de mécanisme d’écoulement
discutable dans les sols. Les calculs réalisés montrent que de grains solides ne pourrait être rencontré que dans le cas d’un
les approches classiques (sans tronquer le critère en traction) ont creusement au tunnelier à pression de terre hors nappe dans
tendance à sous-estimer les pressions de confinement à la rup- les sables. En effet, en méthode conventionnelle, des traitements
ture. Cette sous-estimation reste faible pour φ’ ≥ 20° et c’ ≤ 20 kPa de terrain seraient nécessaires et de la bentonite serait injectée au
(Illustration 12) mais devient significative en cas de terrains peu tunnelier à pression de boue. Ces deux derniers cas conduiraient
frottants et fortement cohérents (φ’ ≤ 20° et c’ ≥ 20 kPa). En à l’existence d’une cohésion sur une faible longueur de terrain
pratique cependant, les terrains ayant de telles résistances au à l’avant du front, empêchant l’apparition de ce mécanisme
cisaillement peuvent vraisemblablement être considérés comme d’écoulement. Dans l’ensemble de ces cas, utiliser l’approche
« non drainés » donc traités par l’approche décrite au §3.2. de l’équation [10] avec une cohésion nulle est sécuritaire.

Illustration 13 : Géométrie du mécanisme de rupture en terrain purement frottant (Berthoz et al., 2012).

20
3.3.2 Terrain stratifié partie supérieure du front (de hauteur β.D) à une résistance
au cisaillement inférieure à celle de la partie inférieure.
La diminution du facteur N γ associée, c’est-à-dire l’aug-
Quelques auteurs ont étendu les approches précédentes au
mentation de la stabilité liée à la taille du front réduite, est
cas de terrains stratifiés horizontalement, avec des couches
donnée en Illustration 16b.
de résistances au cisaillement différentes, mais toujours
avec l’hypothèse qu’elles sont individuellement isotropes.
L’Illustration 14 en dresse un rapide panorama. Les approches d’équilibre limite (EL) permettent assez faci-
lement de tenir compte de plusieurs couches de résistances
Retenons notamment les modélisations physiques en terrains au cisaillement différentes, sur la hauteur du front et dans
stratifiés de Berthoz et al. (2012). Ces essais ont permis d’exhiber l’épaisseur de la couverture, tels les développements
les mécanismes de rupture associés, tel que celui présenté en présentés dans Broere (2001), repris par Vu et al. (2015)
Illustration 15 dans le cas d’un terrain tricouche (partie inférieure puis Champagne de Labriolle (2018). Ces deux modèles
du front cohérente frottante autostable, partie supérieure du considèrent un mécanisme de type « coin et prisme vertical »
front purement frottante, et couverture cohérente frottante). similaire à celui de Horn (1961), repris dans Anagnostou &
Kovari (1996), avec des spécificités relatives notamment à la
Berthoz (2012) a également étendu le mécanisme de Wong prise en compte des effets de voûte, via le rapport K entre
& Subrin (2006) au cas d’un massif stratifié bicouche où la les contraintes horizontale et verticale.

Cas Type Référence Description succincte

Essais sur modèle réduit non centrifugés en terrains bicouche et tricouche
Exp Berthoz et al. (2012) → Forme des mécanismes.

Front découpé en 10 tranches dans la hauteur. Pour chaque tranche, pression au front calculée
Piakowski & Kowalewski (1965) dans chaque tranche en fonction du coefficient de poussée du terrain (théorie de Rankine) et d’un
Milieu continu stratfié drainé

coefficient de forme permettant la prise en compte des effets de voûte dans la direction horizontale.
Approche identique à celle d’Anagnostou & Kovari (1996) mais avec :
sans écoulement

EL Broere (2001) (1) un découpage par tranches du front, permettant l’application à des fronts hétérogènes ;
et Vu et al. (2015) (2) la prise en compte des effets de voûte horizontaux via le choix de K1 = K2 = K0 et le remplacement de
R par R/(1+tanθ).
Approche en sol homogène transposable en raisonnant par couches horizontales de résistance au
Champagne de Labriolle (2018)
cisaillement différentes.
Berthoz (2012) Extension du modèle de Wong & Subrin (2006) au cas d’un milieu stratifié bicouche.
Ibrahim et al. (2015) Extension du modèle de Mollon et al. (2011) au cas d’un milieu stratifié bicouche.

CR Extension du modèle de Mollon et al. (2011) au cas d’un tunnel non circulaire. In fine, le cas étudié
Pan & Dias (2017)
est similaire à celui étudié par Berthoz (2012) et Ibrahim et al. (2015).
Extension du modèle de Mollon et al. (2011) au cas d’une évolution linéaire avec la profondeur de
Zou et al. (2019)
la cohésion et de l’angle de frottement du terrain.
Illustration 14 : Synthèse des approches existantes en milieu continu stratifié drainé sans écoulement.

Phase 3 : bloc 2

Phase 2 : bloc 1

Chambre
d’abattage Phase 1 : extrusion

Axe du
tunnel

(massif tricouche)

Illustration 15 : Mécanisme de rupture en sol stratifié (couche purement frottante entre deux couches cohérentes frottantes), (Berthoz et al., 2012).

21
3
Ces travaux restent cependant assez difficiles à exploiter dans Dans le cas d’un fort contraste, avec une couche très résistante
la mesure où il n’existe pas d’expression directe du résultat : en partie haute (dalle rocheuse), le mécanisme peut simplement
une feuille de calcul est nécessaire. Lorsque les contrastes de être tronqué à la partie inférieure (meuble) du front, la contrainte
résistances entre les différentes couches ne sont pas trop impor- verticale s’exerçant sur le coin étant automatiquement réduite
tants, considérer un milieu homogène de résistance au cisaillement grâce à la cohésion considérée dans la couche supérieure. Si la
égale à celle de la couche réelle la moins résistante peut être couche très résistante est en partie basse du front de taille, il peut
considéré en toute première approche dans le cadre des projets. simplement être considéré une hauteur réduite du front de taille.

Nγ

(a) (b)

Illustration 16 : Terrain stratifié bicouche (Berthoz, 2012) : (a) Géométrie du mécanisme, (b) Abaque.

3.4 MILIEU CONTINU DRAINÉ AVEC ÉCOULEMENT

3.4.1 Écoulement vers le front de taille le régime permanent atteint. Cette hypothèse est légitime dès
que la perméabilité des terrains k est supérieure ou égale à
Lors du creusement en méthode conventionnelle ou avec 10–7 à 10–6 m/s et que la vitesse d’avancement du front vavct
un tunnelier à front ouvert dans un terrain saturé en eau, est inférieure ou égale à 0,1 à 1 m/h (à titre indicatif comme
un écoulement vers le front de taille apparaît avec des forces précédemment d’après Anagnostou & Kovari (1996)). Rappelons
d’écoulement diminuant la stabilité. Les principaux auteurs ayant que l’hypothèse d’une isotropie des terrains est effectuée, à la fois
étudié cette problématique sont synthétisés en Illustration 17. en ce qui concerne la résistance au cisaillement et la perméabilité.
Schuerch et al. (2019) ont également étudié la phase transitoire
La force d’écoulement atteint son maximum lorsque le régime de l’écoulement, en quantifiant la durée pendant laquelle le front
hydraulique permanent est atteint avec une vitesse d’écoulement de taille reste stable à court terme. Cette question est par exemple
maximale. De manière conservative, la plupart des calculs posée lors des opérations hyperbares de maintenance de la
effectués dans ce cas de calcul considèrent par conséquent roue de coupe des tunneliers pressurisés.

Cas Type Référence Description succincte

Essais sur modèle réduit non centrifugés dans des sables secs, saturés sans écoulement, et saturés avec
Exp Lü et al. (2018)
écoulement. Comparaison des mécanismes et des pressions limites de rupture par effondrement.
avec écoulement vers le front

Num Schuerc et al. (2019) Modélisation en éléments finis 3D avec couplage hydro-mécanique pour une analyse en régime transitoire.
Milieu continu drainé

Extension du modèle d’Anagnostou (2012) en prenant en compte les forces d’écoulement vers le tunnel.
EL Perazelli et al. (2014) Pour cela, la charge hydraulique h est supposée ne dépendre que de la distance au front x pour le coin,
avec une expression calée grâce une modélisation numérique aux éléments finis de l’écoulement.
Extension du modèle de Mollon et al. (2011) en ajoutant une composante Nw modélisant les forces
Pan & Dias (2016) d’écoulement (i.γw) au calcul classique en contraintes effectives (fonction de Nγ’). Nw est déduit
d’un calcul d’écoulement sur Flac.
CR Extension du modèle de Pan & Dias (2016) au cas d’un critère de plasticité de Hoek & Brown pris en
compte via une optimisation de la géométrie du mécanisme en fonction de l’angle de frottement
Pan & Dias (2018)
interne correspondant à la linéarisation de H&B et en gardant un couple (c,φ) cohérent avec le
critère H&B initial (si φ augmente, c diminue et Nc augmente).

Illustration 17 : Synthèse des approches consacrées à un milieu continu drainé avec écoulement vers le front.

22
Illustration 18 : Résultats expérimentaux de Lü et al. (2018) : géométrie des mécanismes en présence d’un écoulement (en violet) pour
différentes hauteurs de nappe (H) par rapport à la clé du tunnel. Les mécanismes observés pour des terrains secs (trait bleu pointillé) et
saturés sans écoulement (trait bleu plein) sont également présentés. Les axes des abscisses et ordonnées correspondent aux distances
par rapport au pied du front. L’échelle de couleur correspond aux isovaleurs de déformations plastiques.

(a) (b)

Illustration 19 : Calculs en éléments finis des forces d’écoulement de Perazzelli et al. (2014) : (a) maillage, (b) isovaleurs de charge hydraulique.

Considérons tout d’abord le régime permanent atteint. Lü et al. l’horizontale d’après l’Illustration 18) que pour le même matériau
(2018) ont effectué des essais sur modèle réduit non centrifugé, sec (63°) ou saturé sans écoulement (57°). Le mécanisme s’étend
avec un tunnel de diamètre D = 15 cm et une épaisseur de en conséquence plus loin à l’avant du front en présence des forces
couverture comprise entre 0,5 et 2.D. Le terrain est constitué de d’écoulement, et mobilise un plus grand volume de terrain à la rupture.
sables. Compte-tenu du d10 de ce sable (85 microns environ),
Différents modèles analytiques ont été développés pour évaluer
sa perméabilité est de l’ordre de 5.10–5 à 10–4 m/s d’après la
la pression frontale σ’T nécessaire en présence d’un écoulement.
relation de Hazen (AFTES GT8, 2006). L’instabilité du front est
Dans ces modèles, les forces d’écoulement sont estimées grâce à
générée en reculant une plaque rigide. Trois conditions
des calculs numériques (aux éléments finis ou aux différences
d’essais sont étudiées : sable sec, sable saturé en eau mais
finies) comme le montre l’Illustration 19. Celles-ci sont ensuite
sans écoulement (plaque au front « étanche ») et sable saturé
injectées dans les modèles analytiques de stabilité du front de taille
en eau avec écoulement (plaque au front totalement perméable).
présentés aux paragraphes précédents. Perazzelli et al. (2014)
Ces essais montrent que les forces d’écoulement conduisent à ont à ce titre adapté le modèle d’équilibre limite d’Anagnostou
une inclinaison du pied de la surface de glissement du mécanisme & Kovari (1996), et Pan & Dias (2016, 2018) ont étendu le modèle
plus faible en présence d’écoulement (environ 40° par rapport à de calcul à la rupture de Mollon et al. (2011).

23
3
Les calculs effectués par les différents auteurs montrent que : permanent (à droite) dans un cas particulier. Ces figures sont
déduites de Schuerch et al. (2019) et Perazzelli et al. (2014),
• l’instabilité du front de taille augmente significativement
après uniformisation des notations utilisées. Les points rouges
lorsque la hauteur d’eau augmente (cf. Illustration 20),
correspondent aux hypothèses suivantes avec les notations
ce qui s’explique par l’augmentation de l’amplitude
de l’Illustration 2 : D = 10 m, Hw = H = 15 m, γ = 20 kN/m³,
des forces (motrices) d’écoulement lorsque le gradient
K0 = 1, k = 10–8 m/s, E = 20 MPa, c’ = 25 kPa et φ’ = 25°. Avec
hydraulique augmente ;
ces hypothèses, les calculs montrent que la durée ts critique
• l’instabilité du front de taille augmente lorsque le poids
d’autostabilité du front est égale à 700 s, et qu’une pression σ’T
volumique des terrains diminue (traits pointillés par
égale à 15 kPa doit être appliquée sur le front de taille en régime
rapport aux traits pleins sur l’Illustration 21b) ;
permanent pour en assurer la stabilité.
• les tendances sont similaires, mais un écart existe entre
les résultats obtenus par Perazzelli et al. (2014) et Pan
& Dias (2016, 2018). Cet écart dépend notamment de
l’amplitude de la charge hydraulique (Illustration 20) ;
• une perméabilité horizontale des terrains supérieure à la
perméabilité verticale (anisotropie) améliore la stabilité du
front par rapport au cas isotrope de perméabilité minimale
(Pan & Dias, 2016).

Les analyses menées en régime transitoire par Schuerch et al.

(2019) sur la base de calculs aux éléments finis hydro-mécaniques
couplés permettent d’évaluer la durée pendant laquelle les forces
d’écoulement restent suffisamment faibles pour que le front reste
stable. Différentes valeurs de charge hydraulique, épaisseur de
couverture et coefficients de pression des terres au repos ont
été étudiées. Des abaques en ont été déduits. À titre d’exemple, Illustration 20 : Comparaison entre les solutions existantes avec
l’Illustration 21 présente la durée maximale d’auto-stabilité les hypothèses suivantes : tunnel circulaire de diamètre D, γw / γ’ = 1,0
(à gauche) et la pression de confinement nécessaire en régime (i.e. γ = 20,0 kN/m³), c’ = 0 kPa, φ’ = 30°.

Illustration 21 : (a) Évaluation de la durée limite de stabilité, d’après Schuerch et al. (2019), et (b) évaluation de la pression de confinement
nécessaire pour assurer la stabilité en régime permanent, d’après Perazzelli et al. (2014). Sur ces deux figures, les notations ont été
modifiées pour être cohérentes avec celles de l’Illustration 2. Elles correspondent au cas suivant : Hw = H = 1,5.D, K0 = 1. Les deux points
rouges correspondent au cas décrit dans le corps de texte.

24
3.4.2 Écoulement vers le terrain Broere (2001) a mené des développements analytiques pour
exprimer l’amplitude des surpressions interstitielles Δu en
Le creusement avec un tunnelier à front pressurisé peut fonction de la granulométrie (d10) et de la perméabilité (k) du
générer des surpressions interstitielles dans le terrain à l’avant terrain, de la résistance au cisaillement de la bentonite (τF) et de la
du front de taille, diminuant de fait les contraintes effectives, pression hydrostatique initiale (u0). Le champ de contraintes
donc la résistance au cisaillement sur les surfaces de rupture effectives en tout point du mécanisme de rupture peut en
du mécanisme d’instabilité du front. être déduit et appliqué aux modèles décrits dans le §3.1.

Des mesures effectuées par Xu & Bezuijen (2018) aux Pays-Bas Les études paramétriques effectuées par Broere (2001) montrent
lors du creusement avec un tunnelier à pression de boue dans cependant que la prise en compte des surpressions interstitielles
des sables (10–5 < k < 10–3 m/s) montrent que celles-ci atteignent dans l’étude de la stabilité du front de taille conduit à une
50 kPa à proximité immédiate du front et 10 kPa à 10 m à faible augmentation de la pression frontale (inférieure à 20 kPa)
l’avant du front (Illustration 22a). Lors des phases de pose des dans la grande majorité des cas étudiés.
voussoirs (arrêt du creusement), ces surpressions se dissipent
rapidement jusqu’à s’annuler (retour à la pression hydrostatique). En conséquence, une première approche dans ce cas de
Ce phénomène a également été mesuré par les mêmes auteurs figure (creusement au tunnelier à pression de boue ou au
dans le cas d’un creusement au tunnelier à pression de terre tunnelier à pression de terre avec beaucoup d’adjuvants fins
utilisant des adjuvants (Illustration 22b), avec des pics de sur- dans des terrains sableux de perméabilité comprise entre
pressions qui sont même plus importants (100 kPa à proximité 10–5 et 10–3 m/s) peut être de calculer la pression frontale
immédiate du front, et 20 à 30 kPa à 10 m à l’avant du front) minimale de rupture par effondrement sans tenir compte
dans le cas présenté. Les conditions de creusement associées des surpressions interstitielles, puis d’augmenter celle-ci
ne sont par contre pas décrites précisément par les auteurs. forfaitairement de 30 kPa.

(a) (b)

Illustration 22 : Surpressions interstitielles mesurées à l’avant d’un tunnelier à pression de boue (a) et d’un tunnelier à pression de terre (b)
par Xu & Bezuijen (2018).

25
4 OBSERVATIONS ET
MESURES CONSTRUCTIVES SUR CHANTIER
Cette partie est consacrée aux principes constructifs en méthode conventionnelle et le cas des tunneliers à front
permettant d’assurer la stabilité du front de taille, avec pressurisé. Le cas des tunneliers ouverts est en partie
quelques indications quant au suivi de ces conditions de couvert par les informations indiquées pour le creusement
stabilité. Deux cas sont distingués : le cas du creusement en méthode conventionnelle.

4.1 PROCÉDÉS CONSTRUCTIFS EN MÉTHODE CONVENTIONNELLE

Pour contrôler la stabilité du front de taille, différents procédés sont Rabattement et drainage
possibles en méthode conventionnelle. On retiendra notamment la
Les calculs présentés dans la partie 3.4 ont démontré, si
diminution de la section, le drainage, la réalisation de boulons de
c’était nécessaire, l’impact très négatif de l’eau sur la stabilité
front, puis les traitements de terrains. Les paragraphes suivants
du front de taille. Rabattre la nappe au préalable ou bien
donnent quelques grands principes pour chacun de ces procédés.
drainer le terrain à l’avancement constitue ainsi une excellente
méthode pour réduire les pressions interstitielles à l’avant du
Diminution de la section front de taille et empêcher l’apparition de forces d’écoulement
La stabilité du front est d’abord fonction de la taille de l’ouvrage : significatives défavorables à la stabilité. Il est préférable que
celle-ci est d’autant plus grande que ses dimensions sont ce drainage soit périphérique afin d’éviter les instabilités
réduites. Un creusement en section divisée (calotte puis stross) locales du noyau d’avancement.
ou le maintien d’un merlon central permettent par conséquent
d’améliorer la stabilité, avec toutefois des conséquences non Zingg & Anagnostou (2016) ont effectué des calculs analytiques
négligeables sur la mise en œuvre des soutènements en paroi, et numériques de stabilité du front de taille dans des sols
et plus globalement sur l’avancement du chantier. (0 < c’ < 400 kPa et φ’ = 30°) en faisant varier le nombre, la
longueur et la position de ces forages. Les auteurs concluent,
sur la base d’abaques permettant de quantifier l’effet du
Couche de béton projeté drainage sur la stabilité du front, qu’une très bonne efficacité
Bien qu’elle ne soit pas justifiée par le calcul en pratique, une est obtenue avec 4 à 6 forages de 10 cm de diamètre et de
purge suivie d’une fine couche de béton projeté fibré (2 à 5 cm 1,5.D de longueur situés à proximité de la paroi du tunnel, en
d’épaisseur) doivent être effectuées sur le front de taille après partie supérieure du front.
chaque pas d’excavation afin, en premier lieu, d’assurer la
sécurité des ouvriers vis-à-vis de la chute de « petits blocs ». La mise en place d’un rabattement ou d’un drainage à
l’avancement peut toutefois conduire à des tassements de
L’efficacité du béton projeté résulte en effet de deux actions. consolidation sur une zone étendue autour de l’ouvrage.
D’une part, celui-ci permet de protéger le terrain contre l’altération
superficielle (limitation de la dessication liée à la mise à l’air du
terrain et limitation des circulations d’eau source de délavage
Boulonnage de front
et de forces d’écoulement). D’autre part, cette couche de béton En cas d’insuffisance des procédés précédents, des boulons
projeté a une « action mécanique très locale (d’ordre centimétrique), de front peuvent être mis en œuvre.
mais suffisante pour empêcher les micro-déplacements, et par
suite les micro-ruptures, qui seraient susceptibles de conduire Dans le cas de milieux discontinus, où les mécanismes sont
à une dégradation des caractéristiques du massif en profondeur, de type « chutes de dièdres », quelques boulons à ancrage
par un phénomène de « décohésion en chaîne : le grain de continu de longueur modérée (jusqu’à 6 m) orientés en
sable qui s’échappe libère le caillou qui assurait le calage du fonction des discontinuités sont généralement suffisants. Leur
bloc, et ainsi de suite » (béton projeté de « type 1 » au sens de dimensionnement repose sur des méthodes d’équilibre limite
AFTES GT20, 2001). conformément au §3.1.

26
Dans le cas de milieux continus ou continus équivalents, des L’effet de la voûte parapluie sur la stabilité du front reste
boulons horizontaux longs (généralement 1,5 à 2,0.D, jamais cependant limité. À titre d’illustration, les calculs 2D numériques
moins de 0,5.D) en fibre de verre sont nécessaires. La longueur et analytiques (calcul à la rupture) réalisés par Senent et al. (2020)
des boulons et leur fréquence de renouvellement (généralement dans des terrains cohérents frottants sans écoulement estiment
par tiers) doivent être suffisants à tout instant pour assurer la que la présence de la voûte parapluie réduit la pression
stabilité et contrôler les déformations du front. En particulier, de confinement nécessaire pour assurer la stabilité du front
le calcul devra être effectué en considérant la longueur la plus d’environ 15 %.
faible des boulons, c’est-à-dire juste avant leur renouvellement.
Il convient également de rappeler ici que la réalisation d’une
La longueur, la section et le nombre de ces boulons de front voûte parapluie est un procédé constructif lourd. Compte-tenu
peuvent être calculés par différentes approches. La plus simple du diamètre important des forages à réaliser, il est indispensable
consiste à calculer la contrainte σT nécessaire pour assurer que les tubes soient scellés les uns après les autres, après
la stabilité du front avec les approches décrites au §2.5 et la foration de chaque trou, afin d’éviter un « pré-découpage »
en déduire la force de confinement totale à reprendre par le du terrain et l’apparition de tassements très importants en surface
boulonnage. Le nombre et la section des boulons permettant liée à la phase de mise en place préalable au creusement
de reprendre cet effort tient compte de la longueur d’ancrage proprement dit.
minimale compte-tenu du frottement latéral (qs) mobilisable
autour de chaque boulon, en tenant compte du fait que l’ancrage Traitements de terrain
doit être vérifié sur la longueur « active » (en tête) comme sur
la longueur « passive » (en pied) sans atteindre la rupture en Différents procédés de traitement du terrain visant à augmenter
traction des boulons (Peila, 1994). La recommandation AFTES globalement la résistance au cisaillement du terrain et à réduire
GT30 (2021) peut fournir quelques informations sur ce point sa perméabilité peuvent être réalisés. Citons notamment les
concernant les facteurs de sécurité à utiliser et la détermination injections, la congélation ou le jet-grouting.
de la valeur du qs notamment. L’évaluation des caractéristiques mécaniques des sols traités
n’est pas un sujet immédiat. Une marge de sécurité suffisante
D’autres approches plus élaborées tenant compte explicitement
doit être considérée lors de la conception, et des essais
des boulons dans l’étude de la stabilité du front sont possibles.
de contrôle à l’issue d’un plot d’essai sont recommandés.
Une synthèse (ancienne) des différentes approches de Au stade de la conception, le lecteur intéressé pourra par
dimensionnement est disponible dans Clouterre (2002). exemple se référer aux références suivantes pour obtenir des
Citons notamment les travaux de Leca (1997) et Subrin (2002) premières informations :
s’inscrivant dans l’approche cinématique du calcul à la rupture,
ou ceux plus récents de Anagnostou & Perazzelli (2015) par • pour les injections : des cohésions comprises entre
équilibre limite. Il est également possible de nos jours de quelques dizaines et quelques centaines de kPa sont
réaliser des modélisations numériques 3D avec simulation observées (Dano, 2001 ; Chang et al., 2009 ; etc). L’efficacité
explicite des boulons (Zapata Franco (2020) par exemple), du traitement à grande échelle reste un sujet difficile.
qui permettent en outre d’évaluer l’effet du boulonnage sur L’ingénieur doit rester prudent quant à la valeur de
les déplacements induits dans le terrain. cohésion retenue dans son dimensionnement ;
• pour le jet-grouting : des ordres de grandeur des
cohésions et perméabilités atteignables sont données
Voûte parapluie dans Croce et al. (2014) et Toraldo et al. (2018). Ces
La réalisation d’une voûte parapluie a principalement pour derniers proposent également une méthodologie pour
objectif d’assurer la stabilité du décousu, et de limiter la estimer la résistance au cisaillement des terrains à
pré-convergence et la convergence d’où une limitation des grande échelle à partir d’essais de compression uniaxiale
déplacements générés en surface. Celle-ci a cependant et de tomographie sonique ;
également un effet bénéfique sur la stabilité du front de taille • pour la congélation : ISGF WG2 (1992) propose une
dans la mesure où elle réduit le volume du mécanisme de première estimation des caractéristiques mécaniques
rupture en tronquant la partie supérieure. des sols gelés, notamment l’augmentation de leur
cohésion lorsque la température du sol diminue.
Différentes approches de dimensionnement des voûtes sont
possibles, allant d’un simple calcul aux réactions hyper-
statiques d’un tube isolé soumis à une charge de terrain estimé
grâce à un mécanisme de type Terzaghi, à des modélisations
numériques 3D. Quelques références peuvent être trouvées
dans Gilleron (2016) par exemple.

27
4

4.2 MÉTHODES DE SUIVI EN MÉTHODE CONVENTIONNELLE

En méthode conventionnelle, le front de taille est visible après Des paramètres caractérisant le milieu rocheux fracturé, basés
chaque pas d’avancement. Ce moment doit par conséquent sur le Rock Mass Rating de Bienawski (1989), l’indice Q de
être mis à profit pour valider les hypothèses géotechniques Barton (NGI, 2013) ou le Geological Strength Index de Marinos
(au sens large) retenues lors des études, et en particulier valider et al. (2005) peuvent aussi être déduits des observations et
les conditions de stabilité du front de taille. Les levés de front comparés aux valeurs attendues.
de taille en sont la formalisation. Des exemples de levés du front
sont disponibles en annexes 1 et 2 de AFTES GT24 (2008). Concernant la « validation des conditions
de stabilité du front de taille »
Concernant la « validation Les signes d’instabilité du front de taille sont à relever avec
des hypothèses géotechniques » précision sur les levés du front. Ces signes consistent en
Les informations importantes concernent : des chutes de blocs, l’apparition de fissures au front ou des
mouvements d’extrusion (mesurés grâce à un extrusomètre ou
• la nature des terrains ; à des cibles topographiques en cas d’arrêt prolongé).
• un ordre de grandeur de la résistance au cisaillement
Ces indices doivent permettre au chargé de soutènement
du terrain par des observations empiriques (pelable
de conclure quant à l’auto-stabilité ou non du front de taille, à la
au couteau, difficile à morceler au marteau, etc.), des
géométrie des mécanismes en jeu en cas d’instabilité, et à l’adap-
essais au pénétromètre de poche dans les sols, des
tation des procédés de soutènement du front en cas d’instabilité.
essais Franklin pour les blocs rocheux ;
• l’orientation et l’état des discontinuités dans les massifs Il est également important de vérifier la bonne répartition des
rocheux ; boulons au front (s’assurer que leur déviation reste modérée)
• le nombre et le débit des venues d’eau. ainsi que leur bon scellement au terrain.

4.3 PROCÉDÉS CONSTRUCTIFS AU TUNNELIER À FRONT PRESSURISÉ

Les tunneliers à front pressurisé permettent d’assurer la stabilité confinement est celui exercé par les tunneliers à pression
du front de taille et limiter les déplacements du terrain autour de boue ;
de l’excavation. Différents modes de confinement existent, en • un confinement gazeux, où la chambre d’abattage est
fonction de la nature du matériau contenu dans la chambre remplie d’air comprimé. Ce mode de confinement est
d’abattage : celui exercé par les tunneliers à pression d’air, mais
aussi par les tunneliers à pression de terre lorsque la
• un confinement pâteux, où la chambre d’abattage chambre d’abattage est partiellement remplie, ou par
est principalement remplie du matériau excavé, avec les tunneliers à pression de terre et pression de boue
quelques additifs (eau / mousses / polymères / argile). lors des opérations hyperbares.
Ces additifs visent à rendre le matériau de confinement
suffisamment imperméable afin de limiter les circulations Dans tous les cas, ces dispositifs nécessitent l’évaluation
d’eau internes, mais pas trop visqueux pour que le couple d’une pression de confinement frontale « de consigne »
sur la roue de coupe reste raisonnable, ni trop liquide pour permettant d’assurer la stabilité du front de taille, de contenir les
que le marin puisse être transporté facilement par bande pressions hydrostatiques et, lorsque c’est pertinent, de limiter les
transporteuse. Ce mode de confinement est celui exercé par déplacements admissibles au droit des avoisinants.
les tunneliers à pression de terre avec chambre pleine ;
• un confinement liquide, où la chambre d’abattage est La question de la pression de confinement régnant le long du
remplie de boue. Cette boue doit être assez visqueuse bouclier est également primordiale : si celle-ci est très faible,
pour permettre la création d’une membrane étanche au la stabilité du front en sera dégradée (passage au cas d’un
front (cake) permettant d’exercer réellement une pression creusement avec décousu non nul dans les approches théoriques
sur le squelette granulaire du terrain. Elle ne doit cependant précédemment décrites). Plus largement, l’expérience prouve que
pas être trop visqueuse afin de transporter facilement la pression régnant le long du bouclier conditionne directement
les matériaux excavés jusqu’à la centrale de traitement l’amplitude des déplacements induits en surface dans le cas
des boues et d’être recyclée facilement. Ce mode de d’un creusement à faible profondeur.

28
4.4 MÉTHODES DE SUIVI AU TUNNELIER À FRONT PRESSURISÉ
Le mode « fermé » des tunneliers à front pressurisé ne permet • améliorer la cadence de creusement, grâce par exemple
pas d’obtenir des informations précises relatives à la nature à la possibilité de revoir les pressions de confinement
des terrains, ni aux conditions de stabilité réelles du front en intégrant le retour d’expérience ;
de taille. La borne inférieure du domaine de stabilité a par • mieux maîtriser les impacts sur les avoisinants, lorsque
conséquent peu de chance d’être réévaluée pendant les cela est nécessaire.
travaux. Néanmoins, les interventions hyperbares peuvent
être utilisées pour faire un levé partiel du front. Pour répondre à ces objectifs, le PAT (AFTES GT16, 2018)
doit inclure :
Par contre, le domaine des « déplacements admissibles » • une synthèse des données recueillies lors des travaux :
et la pression de consigne qui en découle sont à recaler (i) les données d’auscultation des avoisinants et du terrain,
régulièrement pendant les travaux grâce à la rétro-analyse des (ii) les données relatives au modèle géotechnique, et
déplacements mesurés dans le cadre du Plan d’Avancement (iii) les données relatives aux conditions de pilotage du
Tunnelier (PAT), fait par sections successives de creusement. tunnelier (pressions dans la chambre d’abattage, effort de
Cet outil, qui se veut dynamique et itératif, a les objectifs poussée sur la roue de coupe, effort de poussée total,
suivants pour le maître d’ouvrage, le maître d’œuvre et mesures au fontimètre, masses extraites...) ;
l’entreprise : • la comparaison de ces données « réelles » avec les
hypothèses retenues lors des études ;
• synthétiser dans un document unique et partagé, et • la réalisation de retro-analyse pour recalibrer certaines
par tronçon du projet, les conditions géologiques, hypothèses des modèles et permettre la mise à jour des pré-
hydrogéologiques et géotechniques attendues, les visions des déplacements induits sur la suite du projet (visant
valeurs de consigne des paramètres-clé de la machine, une amélioration de la fiabilité du scénario de référence).
les contrôles nécessaires, l’identification spatiale des
possibles risques résiduels, etc. ; Pour répondre aux points ci-dessus et être utile pendant le
• synthétiser le retour d’expérience au fur et à mesure chantier, la démarche mise en œuvre doit être suffisamment
de l’avancement et le réutiliser pour les sections de flexible et synthétisée dans des formats relativement simples
creusement successives au sein du même projet, et (par exemple des synoptiques sur le profil en long du tunnel)
aussi pour les projets futurs ; pour être rapide et faciliter la communication.

29
5
SYNTHÈSE

30
31
6
GLOSSAIRE
Paramètres géométriques
D (m) Diamètre du tunnel
C (m) Épaisseur de couverture
H (m) Profondeur de l’axe du tunnel
d (m) Longueur du décousu

Paramètres géologiques, hydrogéologiques et géotechniques

Hw (m) Hauteur d’eau par rapport à l’axe du tunnel
k (m/s) Coefficient de perméabilité du terrain
γ (kN/m³) Poids volumique humide des terrains
γ’ (kN/m³) Poids volumique déjaugé des terrains
Cu (kPa) Cohésion « non drainée » du terrain
c’ (kPa) Cohésion effective du terrain
φ’ (°) Angle de frottement interne du terrain

Paramètres liés au creusement et aux avoisinants

EPB « Earth pressure balanced shield » : tunnelier à pression de terre
SS « Slurry shield » : tunnelier à pression de boue
σs (kPa) Pression en surface
σT (kPa) Pression exercée sur le front de taille (supposée constante)
σT-eff (kPa) Valeur de pression frontale à exercer pour empêcher sa rupture par effondrement
vavct (m/h) Vitesse d’avancement du front de taille

32
7
RÉFÉRENCES
AFNOR (2005), Eurocode 7 : Calcul géotechnique – Partie 1 : Bezuijen, A., & van Seters, A. (2006). The stability of a tunnel face
règles générales. NF EN 1997-1, 138 p. in soft clay. In H. van Lottum, & A. Bezuijen (Eds.), Tunnelling;
a decade of progress (pp. 149-155). Taylor & Francis.
AFTES WG8 (2006), La conception et la réalisation des travaux
d’injection des sols et des roches, Recommandation R2F1, Bienawski Z.T. (1989). Engineering Rock Mass Classification,
Tunnels et Espace Souterrain, n°194-195, pp. 70-154. Wiley, New York, 251 p.

AFTES WG16 (2020), Effects caused by excavation on Broere W. (2001), Tunnel face stability & new CPT application,
neighbouring structures in the design and construction of PhD Thesis, Delft university of technology, Delft university press.
underground works, Recommandation R2A1, 68 p.
Broms B.B., Bennermark H. (1967), Stability of clay at vertical
AFTES WG20 (2001), Design of sprayed concrete for openings, Journal of Soil Mechanics and Foundations Division,
underground support, Recommandation R1A1, Tunnels et vol. 93, pp. 71-94.
Ouvrages Souterrains, n°164, pp.68-108.
CEDD (2015), Catalogue of notable tunnel failures – case histories
AFTES WG24 (2008), Les reconnaissances à l’avancement, (up to April 2015), Civil Engineering and Development Department,
Recommandation R1F1, Tunnels et Espace Souterrain, n°209, Government of Hong Kong, https://www.cedd.gov.hk
pp. 302-325.
Chambon P. & Corté J.F. (1989). Stabilité du front de taille
AFTES WG29 (2007), Compatibilité des recommandations d'un tunnel faiblement enterré: modélisation en centrifugeuse,
AFTES relatives aux revêtements des tunnels en béton avec Proceedings of the Int. Conf. Tunnelling and Microtunnelling in Soft
les Eurocodes, Recommandation R2F1, Tunnels et Ouvrages Ground: from field to theory, Paris, presses ENPC, pp. 307-315.
souterrains, n°204, 18 p.
Chambon, P., Corte, J.F. (1994). Shallow tunnels in cohe-
AFTES WG30 (2021), Radial rock bolting in tunnels – Design sionless soil: stability of tunnel face, Journal of Geotechnical
and sizing guide, Recommandation R1A1, 97 p. Engineering 120 (7), pp. 1148–1165.

Alagha S.N, Chapman D.N. (2019), Numerical modelling of Champagne de Labriolle G. (2018), Détermination de la fenêtre
tunnel face stability uin homogeneous and layered soft ground, de pilotage de la pression de confinement d’un tunnelier
Tunnelling and Underground Space Technology, vol. 94, fermé dans un sol cohérent-frottant ou purement cohérent,
pp. 1-14. Revue Française de Géotechnique, vol. 155, n°3.

Anagnostou G. (2012), The contribution of horizontal arching Chang M., Chen C.C., Huang R.C., Chang J., Yang P.J. (2009),
to tunnel face stability, Geotechnik, vol 35, n°1, pp34-44. Investigation on mechanism of grouting and engineering
characteristics of in-situ grouted soils, Proceedings of the
Anagnostou G., Kovari K. (1996), Face stability conditions with 17th International Conference on Soil Mechanics and Geotechnical
Earth Pressure Balanced Shield, Tunnelling and Underground Engineering, M. Hamza et al. (Eds.), pp. 2346-2349,
Space Technology, vol 11, n°2, pp. 165-173.
Clouterre (2002), Additif 2002 aux recommandations Clouterre
Anagnostou G., Perazzelli P. (2015), Analysis method and 1991, pour la conception, le calcul, l’exécution et le contrôle
design charts for bolt reinforcement of the tunnel face in des soutènements réalisés par clouage des sols, Presses de
cohesive-frictional soils, Tunnelling and Underground Space l’École Nationale des Ponts et Chaussées, 217 p.
Technology, vol 47, pp. 162-181.
Croce P., Flora A., Modoni G. (2014), Jet-grouting technology, design
Berthoz N. (2012), Modélisation physique et théorique du and control, London, UK : CRC Press (Taylor & Francis Group).
creusement pressurisé des tunnels en terrains meubles
homogènes et stratifiés, Thèse de doctorat de l’ENTPE, 292 p. Dano C. (2001), Comportement mécanique des sols injectés,
Thèse de doctorat de l'Université de Nantes, 216 p.
Berthoz N., Branque D., Subrin D., Wong H., Humbert E. (2012),
Face failure in homogeneous and stratified soft ground : Theoretical DAUB (2016), Recommendations for Face Support Pressure
and experimental approaches on 1g EPBS reduced scale model, Calculations for Shield Tunnelling in Soft Ground, German
Tunnelling and Underground Space Technology, vol. 30, pp. 25-37. Tunnelling Committee, 64 p.

33
7
Davis E.H., Gunn M.J., Mair R.J., Seneviratne H.N. (1980). Marinos V., Marinos P., Hoek E. (2005), The geological strength
The stability of shallow tunnels and underground openings in index : applications and limitations, Bulletin of Engineering
cohesive material, Géotechnique, n°40, pp. 397-416. and Geological Environment, n°64, pp. 55-65.

Gilleron N. (2016), Méthode de prévision des tassements Messerli J., Pimentel E., Anagnostou G. (2010), Experimental
provoqués par le creusement des tunnels urbains et influence study into tunnel face collapse in sand, Proceedings of the
des présoutènements, thèse de doctorat de l’université 7th International Conference in Physical Modelling in Geotechnics,
Paris-Est, 238 p. Zurich, Switzerland.

Grasso P., Xu S., Del Fedele M., Russo G., Chiriotti E. (2003). Mollon G., Dias D., Soubra A.H. (2011), Rotational failure
Particular failure mechanisms of weathered granite observed mechanisms for the face stability analysis of tunnels driven
during construction of metro tunnels by TBM, ITA World Congress, by a pressurized shield, International Journal for Numerical and
Amsterdam. Analytical Methods in Geomechanics, vol 35, n°12, pp. 1363-1388.

Gunn MJ. (1980), Limit analysis of undrained stability problems Mollon G., Dias D., Soubra A.H. (2013), Continuous velocity fields
using a very small computer, Proc., Symp. on Computer for collapse and blowout of a pressurized tunnel face in purely
Applications to Geotechnical Problems in Highway Engineering, cohesive soil, International Journal for Numerical and Analytical
Cambridge Univ., Cambridge, UK, pp. 5–30. Methods in Geomechanics, vol 37, n°13, pp. 2061-2083.

Horn, M. (1961), Alagutak homlokbiztositstasra hato vizszintes NF EN 1997-1 (2005), Eurocode 7 : Geotechnical Design – Part 1:
földnyomasgalat néhany eredménye, Azorszdgos mélyépitoipari General Rules.
konferencia éloadasai, Közlekédési Dokumentacios Vallalat,
NGI (2013), Using the Q-system – Rock mass classification
Budapest (in Hungarian).
and support design, Handbook of the Norvegian Geotechnical
Institute, 57 p.
Ibrahim E., Soubra A.H., Mollon G., Raphael W., Dias D.,
Reda A. (2015), Three-dimensional face stability analysis of
Pan Q., Dias D. (2016), The effect of pore water pressure on
pressurized tunnels driven in a multilayered purely frictional
tunnel face stability, International Journal for Numerical and
medium, Tunnelling and Underground Space Technology,
Analytical Methods in Geomechanics, vol. 40, n°15. pp. ..-..
vol. 49, pp. 18-34.
Pan Q., Dias D. (2017), Upper-bound analysis on the face
ISGF WG 2 (1992), Frozen ground structures – Basic principles
stability of a non-circular tunnel, Tunnelling and Underground
of design, by Andersland O., Berggren A.L and Fish A.,
Space Technology, vol. 62, pp. 96-102.
Ground Freezing 91, vol 2. pp. 503-513.
Pan Q., Dias D. (2018), Three dimensional face stability of a tunnel
Kimura, T. and Mair, R. J., (1981), Centrifugal testing of in weak rock masses subjected to seepage forces, Tunnelling
model tunnels in soft clay, Proceedings of 10th International and Underground Space Technology, vol. 71, pp. 555-566.
Conference on Soil Mechanics and Foundation Engineering,
Stockholm 1981, Vol. 1, pp 319-322. Paternesi A., Schweiger H.F., Ruggeri P., Fruzzetti V.M.E.,
Scarpelli G. (2017), Comparisons of Eurocodes design
Leca E. (1997), Développement d’outils de calcul pour le approaches for numerical analysis of shallow tunnels, Tunnelling
dimensionnement des tunnels creusés en terrains meubles, and Underground Space Technology, vol. 62, pp. 115-125.
Thèse d’habilitation, Université de Lille, 83 p.
Peila D. (1994), A theoretical study of reinforcement influence
Leca E., Dormieux L. (1990), Upper and lower bound solutions on the stability of a tunnel face, Geotechnical and Geological
for the face stability of shallow circular tunnels in frictionnal Engineering, n°12, pp. 145-168.
material, Geotechnique, vol. 40, n°4, pp. 581-606.
Perazzelli P., Anagnostou G. (2017), Analysis Method and
Li T.Z., Yang X.L. (2020), Stability of plane strain tunnel Design Charts for Bolt Reinforcement of the Tunnel Face in
headings in soils with tensile strength cut-off, Tunnelling and purely Cohesive Soils, J Geotech Geoenviron Eng 143 (9).
Underground Space Technology, vol. 95, pp. 1-11.
Perazzelli P., Leone T. Anagnostou G. (2014), Tunnel face stability
Li Z.W., Yang X.L., Li T.Z. (2019), Face stability analysis of under seepage flow conditions, Tunnelling and Underground
tunnels under steady unsaturated seepage conditions, Tunnelling Space Technology, vol. 43, pp. 459-469.
and Underground Space Technology, vol. 93, pp. 1-8.
Piakowski A., Kowaleski Z. (1965), Application of thixotropic
Lü X., Zhou Y., Huang M., Zeng S. (2018), Experimental study clay suspensions for stability of vertical sides of deep trenches
of the face stability of shield tunnel in sands under seepage without strutting, Proceedings of 6th International Conference
condition, Tunnelling and Underground Space Technology, on Soil Mechanics and Foundations Engineering, Montreal,
vol. 74, pp. 195-205. vol. 111.

34
Quarmout M., König D., Gussmann P., Thewes M., Schanz T. Viana Da Fonseca A., Topa Gomes A. (2011), A tunnel collapse
(2019), Tunnel face stability analysis using Kinematical Element on the construction in metro do Porto : solutions for optimization
Method, Tunnelling and Underground Space Technology, vol. 85, of advance control parameters of a EPB TBM, Proceeding, of
pp. 354-367. the XXV technical conference, Miedzyzdroje, 24-24 mai 2011.

Salençon J. (1990), An introduction to the yield design theory Vu M.N., Broere W., Bosch J. (2015), The impact of shallow
and its application to soil mechanics, European Journal of cover on stability when tunnelling in soft soils, Tunnelling and
Mechanics A-Solids, vol 9, n°5, pp. 477-500. Underground Space Technology, vol. 50, pp. 507-515.

Schofield AN. (1980), Cambridge geotechnical centrifuge Wong H., Subrin D. (2006), Stabilité frontale d'un tunnel : mécanisme
operations, Geotechnique, vol. 30, n°3, pp. 227-268. 3D en forme de corne de rhinocéros et influence de la profondeur,
Revue européenne de génie civil, vol. 10, n°4, pp. 429-456.
Schuerch R., Poggiati R., Anagnostou G. (2019), Design charts
for estimating face stand-up time in soft ground tunnelling, Xu T., Bezuijen A. (2018), Analytical methods in predicting
Proceedings of the World Tunnel Congress, Naples, 10 p. excess pore water pressure in front of slurry shield in saturated
sandy ground, Tunnelling and Underground Space Technology,
Senent S., Yi C., Jinenez R. (2020), An upper bound solution vol. 73, pp. 203-211.
for tunnel face stability analysis considering the free span,
Tunnelling and Underground Space Technology, vol. 103, Zapata Franco D., Janin J.P., Dano C., Le Bissonnais H.,
pp. 1-14. Falconi F., Gérardin C. (2020), Méthodes 3D simplifiées pour
évaluer l’impact du boulonnage au front : REX de la L11-GC03,
Subrin D. (2002), Études théoriques sur la stabilité et le Paris, Journées Nationales de Géotechnique et de Géologie
comportement des tunnels renforcés par boulonnage, thèse de l’Ingénieur, Lyon, 8 p.
de doctorat de l’INSA de Lyon, 210 p.
Zingg S., Anagnostou G. (2016), An investigation into efficient
Toraldo C., Modoni G., Ochmanski M., Croce P. (2018), The drainage layouts for the stabilization of tunnel faces in
characteristic strength of jet-grouted material, Géotechnique, homogeneous ground, Tunnelling and Underground Space
vol. 68, n°3, pp. 262-279. Technology, vol. 58, pp. 49-73.

Ukrichton B., Yingchaloenkitkhajorn K., Keawsawasvong S. Zou J., Chen G., Qian Z. (2019a), Tunnel face stability in cohesion-
(2017), Three-dimensional undrained tunnel face stability frictional soils considering the soil arching effect by improved
in clay with a linearly increasing shear strength with depth, failure models, Computers and Geotechnics, vol. 106, pp. 1-17.
Computers and Geotechnics, vol. 88, pp. 146-151.
Zou J.F., Qian Z.H., Xiang X.H., Chen G.H. (2019b), Face stability
Vermeer P.A., Ruse N., Marcher T. (2002), Tunnel heading stability of a tunnel excavated in saturated nonhomogeneous soils,
in drained ground, Tunnelling, Felsbau, vol. 20, n°6, pp. 8-18. Tunnelling and Underground Space Technology, vol. 83, pp. 1-17.

35
8 ANNEXE : RETOURS D’EXPÉRIENCES
D’INSTABILITÉS DE FRONT
Quelques retours d’expériences sont présentés dans les tableaux suivants. Si vous disposez de données relatives à d’autres
effondrements, veuillez les envoyer dans un format similaire à [email protected] afin qu’ils puissent
être inclus dans une future mise à jour de ce document d’information.

Nom du tunnel : Schirmeck (France, 68)

Méthode de creusement : conventionnelle

Année : 2004

Nature des terrains :

Zone de transition entre (i) un faciès rocheux fortement fracturé
modérément altéré plutôt sec, et (ii) un second faciès très
fortement fracturé et complètement altéré avec des venues d’eau
assez importantes.

Observations et adaptation des méthodes de soutènement :

Au PM 218, à la fin du marinage, le front de taille s’est effondré
sur presque toute sa largeur selon un plan de fracturation avec
un pendage de l’ordre de 45° vers le tunnel. Le volume total de
l’effondrement était de l’ordre de 100 m³ (illustration a). (a) Effondrement du front au PM 218.
Le creusement a pu être poursuivi grâce à une stabilisation du front
par 28 m³ de béton projeté et 30 boulons expansifs longitudinaux
de 4 m de longueur.
Compte-tenu des conditions précaires de stabilité observées, le
boulonnage de front par des boulons expansifs de 4 m de longueur a
été systématisé jusqu’au PM 465. Cela a permis d’empêcher l’apparition
d’un nouveau mécanisme global d’instabilité, mais des instabilités
locales ont persisté, telle celle du PM 235 visible en illustration (b).

Conséquences globales pour le projet :

La durée totale des travaux de creusement-soutènement a été
de l’ordre du double de celle attendue. Cela est lié notamment à
l’augmentation du linéaire de tunnel cintré (au détriment d’un profil
boulonné), accompagnée, de fait, d’une épaisseur de béton projeté
plus importante. La mise en œuvre de davantage de boulons de front y
contribue également, sans que cela puisse être quantifié précisément. (b) Effondrement localisé
en piédroit du front au PM 235.

36
Nom du tunnel : Bois de Peu (France, 25)

Méthode de creusement : conventionnelle

Année : 2006

Nature des terrains :

Argiles humides.
Observations :
Trois ruptures de front majeures, de volume supérieur à 50 m³ se
sont produites aux PM 518, 515 et 510 de la contre-attaque du
tube descendant sous une épaisseur de couverture de l’ordre d’une
dizaine de mètres. Les boulons en fibre de verre, mis en place au
front avec une densité de 0,25 boulons par m², n’ont pas empêché
les instabilités de front qui ne se sont toutefois pas propagées en
surface en raison de la voûte parapluie initialement mise en place.
Confortement et adaptation des méthodes :
Poursuite du creusement en passant en section divisée (calotte
et stross) pendant une quinzaine de mètres, jusqu’à ce que la
résistance au cisaillement des terrains s’améliore.
Conséquences pour le projet : Vue du front après un des effondrements. Les boulons de front
Augmentation du délai de travaux d’environ 1 mois, et surcoût d’environ ainsi que la voûte parapluie sont clairement visibles.
500 k€ HT d’après le Décompte Général Définitif du marché de travaux. © Nicolas Dupriez, DDE 25.

Nom du tunnel : Métro de Porto (Portugal)

Méthode de creusement : tunnelier à front pressurisé (pression de terre)

Année : 2001

Nature des terrains :

Front mixte et hétérogène composé de granites sains en partie basse
du front de taille, et granites très altérés en partie supérieure.
Observations :
Trois instabilités locales observées lors du creusement des premiers
600 m de la ligne C du métro, dont une conduisant à un fontis
engendrant l’effondrement d’un immeuble, et causant la mort d’une
habitante. Ce fontis est apparu environ 30 m après le passage du
front de taille. Celui-ci a été généré par la sur-excavation de la partie
supérieure du front, liée au contraste de résistance entre les granites
sains et altérés et par une difficulté initiale (d’ordre technique et
organisationnel) à garantir la chambre pleine et pressurisée.
(a) Fontis apparu en surface.
Confortement et adaptation des méthodes :
Suite à cet accident, le conditionnement a été revu, les équipes de
chantier ont été renforcées, les procédures de suivi ont été revues, et le
PAT (plan d’avancement tunnelier) pour chaque tronçon de creusement
de 500 m a été introduit pour la première fois. En outre, un dispositif
d’injection de bentonite sous pression a été ajouté sur le tunnelier. Celui-ci
se déclenchait automatiquement lorsque la pression mesurée en calotte
de la chambre d’abattage était inférieure à une valeur prédéfinie.
La fin de la ligne C et la ligne S n’ont pas vu l’apparition d’autres fontis.
(b) Coupe longitudinale illustrant le front mixte rencontré
Conséquences pour le projet : et le système complémentaire d’injection de bentonite en clé
Le décès d’une tierce personne a nécessairement eu un impact très négatif ajouté suite à l’effondrement.
sur le démarrage de ce projet. Les conséquences financières associées Référence : Grasso et al. (2003),
ne sont pas publiées. Le projet a subi un retard de 6 mois environ. Viana Da Fonseca & Topa Gomes (2011).

37
8

Nom du tunnel : Métro de Rennes (France, 35)

Méthode de creusement : tunnelier à front pressurisé (pression de terre)

Année : 2016

Nature des terrains :

Schistes plus ou moins altérés.

Observations :
Plusieurs fontis en surface ont été générés par le creusement de cette
ligne de métro. Citons notamment celui d’environ 20 m² et au moins
3 m de profondeur généré dans le magasin Noz. Ces fontis sont liés à
la structure très particulière du massif rocheux, constitué de bancs très
inclinés (70°), dont les discontinuités avaient tendance à s’ouvrir et
à être lubrifiées sous l’effet de la pression de confinement.

Confortement et adaptation des méthodes :

Ajout sur le tunnelier en cours de chantier d’un système d’injection
de bentonite (bulle intégrée), et contrôle accru de la pression de
Photographie du fontis apparu dans le magasin Noz.
confinement en veillant à ce qu’elle ne soit pas trop importante. Celui-ci a une surface d’environ 20 m² et 3 m de profondeur.
Trois personnes sont tombés dans ce fontis,
Conséquences pour le projet : mais ont pu s’en extraire sans gravité.
Arrêt du tunnelier pendant 133 jours, soit un surcoût de l’ordre de © Journal Ouest-France du 19/11/2016.
10 M€. Référence : fiche tunnelier AFTES GT4 n°126, TES n°274.

Nom du tunnel : Métro de Nice (France, 06)

Méthode de creusement : tunnelier à front pressurisé (pression de boue)

Année : 2017

Nature des terrains :

Fronts mixtes composés de terrains indurés (calcaires et poudingues)
et de terrains meubles (sables et galets).

Observations :
Plusieurs fontis ont été induits en surface, dont un de 6 m de
diamètre et 1,5 m de profondeur environ au milieu de la rue de
France. Ce dernier n’a pas engendré de dommages aux constructions
avoisinantes ni d’accident de la circulation. Ce fontis a été généré
par la rencontre par le tunnelier d’une ancienne cavité remplie de
matériaux de construction, notamment des aciers coincés dans la
roue de coupe, avec une perte de pression associée dans la chambre.

Confortement et adaptation des méthodes : Photographie du fontis généré dans la rue de France.
© Patrick Allemand.
Des injections ont été réalisées sur une quinzaine de mètres autour
Référence : Batiactu.com, 05/07/2017 et 25/08/2017.
du fontis afin d’améliorer la stabilité des terrains, puis la roue de
coupe a été dégagée, permettant la reprise du creusement.

Conséquences pour le projet :

Environ 1,5 mois d’arrêt du tunnelier.

38
Auteurs : Nicolas BERTHOZ, Didier SUBRIN.

Relecteurs internes : Michel DEFFAYET, Eric PREMAT, Cédric GAILLARD, Gilles HAMAIDE, Johan KASPERSKI.

Relecteurs externes : Denis BRANQUE (ENTPE), Martin CAHN (TERRASOL-SETEC), Guillaume CHAMPAGNE DE
LABRIOLLE (ARCADIS), Elena CHIRIOTTI (INCAS PARTNERS), Fabrice EMERIAULT (Grenoble INP), Jean-Pierre JANIN
(LOMBARDI INGENIERIE), Christophe JASSIONNESSE (SPIE BATIGNOLLES), Adrien SAITTA (EGIS TUNNELS).

39
 Crédits photos : CETU / Couverture : Nicolas Dupriez (DDE 25) - Rédaction et coordination : collaborateurs CETU - Réalisation : PAO Concept - Impression : VASSEL GRAPHIQUE

Centre d’Études des Tunnels

25 avenue François Mitterrand MINISTÈRE

69500 BRON - FRANCE CHARGÉ
Tél. +33 (0)4 72 14 34 00 DES TRANSPORTS
Fax. +33 (0)4 72 14 34 30
[email protected]

www.cetu.developpement-durable.gouv.fr