CHAPITRE 4: Calcul des ouvrages souterrains

1.Introduction :
La vérification de la stabilité d’un ouvrage souterrain et la détermination des caractéristiques de son
soutènement et de son revêtement est un problème particulièrement complexe .Le comportement des
matériaux naturels sont extrêmement varié et difficile à prévenir.
L’existence de solutions explicites au problème de l’équilibre d’un massif excavé et de son soutènement
est limitée et est basée sur des hypothèses simples telle que :
- La géométrie : le tunnel est supposé à section simple .
- La stratigraphie : on ne prend en compte qu’une seule couche de terrain supposée homogène.
- La loi de comportement du sol : supposée élastique linéaire ou élasto-plastique.
- L’état initial de contraintes est supposé isotrope et homogène.
- Les formules sont exprimées dans le plan et dans le cadre des petites déformations.

2. Méthode des réactions hyperstatiques

Cette méthode , économique , relativement simple à comprendre et à utiliser, reste une
approximation dans le cas des ouvrages construit en terrain meuble ou en rocher fracturer sous faible
couverture et avec des soutènement traditionnel lourd.

2.1. Principe de la méthode

Le principe est d’étudier le comportement du soutènement (ou du revêtement) sous l’action de
charges extérieures (Fig. 1 ). On réalise donc un calcul de structure classique .La géométrie du soutènement
est rentrée précisément pour un mètre linéaire de galerie sous forme de poutres 2D, puis on vient lui
appliquer un chargement.

Fig .1 – Schéma classique d’un modèle aux réactions hyperstatiques,

avec les forces actives et passives(ressorts).
Donc ,On distingue deux types de charges :
• Charges actives : se sont des charges indépendantes de l’état de déformation :
- pression appliquée par le poids du terrain (verticale et horizontale) ;
- pression hydrostatique si le tunnel traverse une nappe ;
- gonflement éventuel ;
- détachement d’un bloc ;
- Poids propre d’un revêtement ;
- Circulation routière à faible profondeur

1
 • Charges passives : se sont les réactions hyperstatiques issues de la déformation du soutènement.
Ces charges sont considérées comme linéairement liées aux déplacements, ce qui permet de les
modéliser par une série de ressorts, dont la rigidité K est issue des propriétés mécaniques de la roche
ou du sol environnant sont considérées comme linéairement liées aux déplacements, ce qui permet
de les modéliser par une série de ressorts,

L’équilibre de la structure établi, il est alors possible d’accéder aux efforts dans le soutènement (M, N
et T) ainsi qu’aux convergences maximales

2.2. Détermination des charges actives

Il existe de nombreuses théories permettant d’évaluer les pressions appliquées sur une structure de
soutènement. Dans un milieu continu, La méthode des réactions hyperstatiques étant plutôt adaptée aux
tunnels peu profonds dans des sols, on utilise principalement les formules de Terzaghi ou de Caquot qui
sont basées sur une rupture en coin du terrain en piédroits et l’effet du poids d’une voûte de décharge. La
formule de Terzaghi s’exprime ainsi

Section rectangulaire ou en fer à cheval section circulaire

2𝐶
𝑏(𝛾− ) −2𝐻𝑡𝑎𝑛𝑔𝜑
𝑏
𝑃𝑣 = (1−𝑒 𝑏 )
2𝑡𝑎𝑛𝑔𝜑

𝑃ℎ = 𝐾0 𝑃𝑣 ou 𝑃ℎ = 𝐾𝑎 𝑃𝑣
Pv : la pression verticale
Ph : la pression horizontale (pression sur les parois latérales)
H et b : profondeur du tunnel et la largeur du cône d’effondrement estimé en clef de voûte
C, 𝜑 et 𝛾 : les paramètres de Coulomb et le poids volumique du terrain.
Ka : coefficient de poussée
K0 : coefficient des terres au repos

N.b : Dans certains cas, il peut être intéressant de modéliser la chute d’un bloc rocheux
sur un soutènement/revêtement. Il faut alors vérifier la stabilité de la structure sous deux types de charge
active
– Chute d’un bloc en voûte. On ne considère que le poids propre du bloc, situé en clef de voûte ;

2
 – Chute d’un bloc en piédroit. Ce cas de charge dissymétrique est particulièrement défavorable. Il
faut considérer le glissement du bloc sur une ou deux discontinuités

2.3 Détermination des charges passives

Plus la difficulté de savoir jusqu’où placer les ressorts, il est souvent difficile, sans essais sur le
matériau, de déterminer la valeur du module K.
L’expression analytique du déplacement 𝑢 ⃑⃑⃑ en paroi d’un tunnel circulaire pour un massif de roche
élastique, isotrope et linéaire permet d’accéder à une expression approximative de ce module :

𝐸
𝑞 = 𝐾𝑢
⃑ 𝑘 = (1+𝜈) 𝑅
⃑q : pression appliquée en paroi.
K :module de réaction est issue des propriétés mécaniques de la roche ou de sol environnant .

3.Méthodes analytiques
3.1. Méthode analytique élastique
Cette méthode traite l’entourage de l’excavation comme un milieu élastique. Généralement, Les
équations d’équilibre de la mécanique sont résolues en déformation plane (𝜀𝑧 = 0).
En coordonnées cylindrique r, θ ces équations se mettent en forme :

𝜕𝜎𝑟 1 𝜕𝜏𝑟𝜃 𝜎𝑟 −𝜎𝜃

+𝑟 + =0
𝜕𝑟 𝜕𝜃 𝑟

𝜕𝜎𝑟𝜃 1 𝜕𝜎𝜃 2 𝜏𝑟𝜃

+𝑟 +𝑟 =0
𝜕𝑟 𝜕𝜃 𝜃

Les conditions aux limites sont:

Pour r =R : 𝜎𝑟 = 0 et 𝜏𝑟𝜃 = 0
Pour 𝑟 → ∞ :(𝜎𝑟 , 𝜏𝑟𝜃 , 𝜎𝜃 ) = (𝜎𝑟0 , τ0𝑟𝜃 , 𝜎𝜃0 ) état de contrainte initiale

Dans le cas d’un tunnel non revêtu à section circulaire (rayon R) creusé dans un massif homogène isotrope
élastique on a :
𝑑𝜎𝑟 𝜎𝑟 −𝜎𝜃
L’équation d’équilibre
𝑑𝑟
+ 𝑟
=0

𝑟∞ 𝜎𝜃 = 𝜎 0
Les conditions aux limites {
𝑟=𝑅 𝜎𝑟 = 0

𝜎𝑟 = 𝜎 0 (1 − 𝛼 2 ) 𝑅
Les contraintes { /𝛼 =
𝜎𝜃 = 𝜎 0 (1 + 𝛼 2 ) 𝑟

1+𝜈 𝑅2
Déplacement radiale ur = − 𝐸
𝜎0 𝑟

3
3.2. Méthode convergence-confinement
Les idées et théories sont liées à ces deux termes : convergence et confinement.
-La convergence est liée à un déplacement.
-Le confinement est la pression radiale qui s’applique sur le pourtour de l’excavation, en présence d’un
soutènement. Il constitue en quelque sorte le chargement du soutènement. On parle aussi de déconfinement,
mais pour le terrain. Il s’agit de la décompression causée par la présence du tunnel. Ce déconfinement
s’amorce bien en avant du front de taille
3.2. 1.Hypothèses
L’hypothèse forte est la considération unidimensionnelle du problème :
– La cavité étudiée a une forme cylindrique.
– Hypothèse d’isotropie des contraintes initiales (K0 = 1) et d’isotropie du massif ;
– Hypothèse des déformations planes ;
L’état initial est défini par l’état de contrainte isotrope. H est la hauteur de couverture et le poids
volumique des terrains sus-jacents. La contrainte initiale dans le massif est donc :
𝜎0 = 𝛾𝐻
3.2. 2.Courbe de convergence
On introduit une pression fictive 𝑃𝑖 en paroi. Cette pression, uniformément répartie sur le pourtour
de l’excavation, a une valeur qui décroît avec l’éloignement au front. L’évolution de Pi est donc gouvernée
par la distance x (Fig.2), qui permet de se situer par rapport au front de taille (où x = 0). On écrit :
Pi = (1 − λ(x))σ0 0 ≤ λ(x) ≤ 1 , 0 ≤ Pi ≪ σ0
λ(x) est appelé taux de déconfinement car il caractérise l’état du massif à l’endroit x considéré.
𝜆(𝑥) = 0, Pi =𝜎0 est l’état de contrainte initial, en avant du front de taille.
λ(x) = 1 , Pi =0 est l’état complètement déconfiné, loin en arrière du front.

Fig.2 : Notion de pression fictive et de déconfinement autour du front de taille

Formules la plus simple de λ(x) , lorsque le terrain reste en élasticité est :

4
 𝐦𝟎 𝐑 𝟐
𝛌(𝐱) = 𝛂 + (𝟏 − 𝛂)(𝟏 − ( ) )
𝐦𝟎 𝐑 + 𝐱
α et m0 sont deux constantes (on prend respectivement 0.25 et 0.75) et R le rayon d’excavation.

La courbe de convergence (Fig.3 ) est la courbe donnant la valeur du déplacement en paroi u en

fonction de la pression fictive Pi , et en l’absence de soutènement. Il s’agit d’une courbe paramétrique (de
paramètre x ) qui se représente sur un graphe (Pi, u).

Fig. 3– Courbes de convergence, de confinement et assemblage.

La théorie de l’élasto-plasticité permet d’obtenir l’équation de cette courbe pour des critères simples
(Mohr-Coulomb par exemple). Le détail des calculs sort du cadre de ce cours, nous nous contenterons
simplement de distinguer deux phases et de donner les principaux résultats :
• Phase élastique ( de u = 0 à u = uic) : Le terrain se déconfine progressivement de manière linéaire.
Par prolongement de la droite sur l’axe des abscisses, on obtient le pseudo-déplacement élastique ue.
L’expression simple de ue donne une première approximation des déplacements en tunnels :
(𝟏 + 𝝂)
𝒖𝒆 = 𝑹𝝈𝟎
𝑬

• Phase plastique(de u = uic à u = uinf ) : Le terrain en paroi passe dans un état de déformations
irréversible. Il y a rupture par excès de compression, par écrasement. Parfois la courbe ne recoupe
pas l’axe des abscisses, et la paroi se referme sur elle-même (très grandes déformations). L’état
plastique est à éviter, c’est un des rôles du soutènement. Pour le critère de Mohr-Coulomb.
L’équation de la courbe plastique est :
(𝟏 + 𝝂) 𝑹 𝑹𝒑
𝒖 =𝑹 [𝑪𝟏 + 𝑪𝟐 ( )𝑲𝒑 −𝟏 + 𝑪𝟑 ( )𝜷+𝟏 ]
𝑬 𝑹𝒑 𝑹

(1−𝜈)(1+𝛽𝐾𝑝) 2(𝜎0 +𝐻)

𝐶1 = −(1 − 2𝜈)(𝜎0 + 𝐻) / 𝐶2 = ( − 𝜈)
𝐾𝑝 +𝛽 𝐾𝑝+1
1
(𝐾𝑝 −1)(𝜎0 +𝐻) 2(𝜎0 +𝐻)𝑅𝐾𝑝 −1 𝐾𝑝−1
𝐶3 = 2(1 − 𝜈) 𝐾𝑝+𝛽
/ 𝑅𝑃 = [ (𝐾 ]
𝑝 +1)(𝑃𝑖 +𝐻)

Les constantes du modèle de Mohr-Coulomb non-associé (angle de dilatance 𝜓 ) sont exprimées

ainsi :

5
 𝜋 𝜑 𝐶 1+𝑠𝑖𝑛𝜓
𝐾𝑝 = 𝑡𝑎𝑛 2 (4 + 2 ) , 𝐻 = tan 𝜑 , 𝛽 = 1−𝑠𝑖𝑛𝜓
Pour tracer la courbe il faut calculer la pression Pic d’apparition de la plasticité :

𝟐𝝈𝟎 − 𝑯(𝑲𝑷 − 𝟏)
𝑷𝒊𝒄 =
𝑲𝑷 + 𝟏
Lorsque le terrain est de bonne qualité, dans les roches dures par exemple, il se peut que la paroi reste en
élasticité durant tout le déconfinement. Un critère très utilisé en travaux souterrains pour déterminer si le
massif encaissant risque de rentrer en plasticité est le facteur de stabilité :
2𝜎0 𝑐𝑜𝑠𝜑
𝐹= ; 𝑅𝑐 = 2𝑐 1−𝑠𝑖𝑛𝜑
𝑅𝑐
Rc :résistance à la compression .
F > 1 : risque d’instabilité.
3.2. 2.Courbe de confinement
Une deuxième courbe est requise pour la méthode. Il s’agit de la courbe de confinement (Fig.2),qui
va permettre de caractériser le comportement du soutènement sous son chargement .Le chargement
considéré est purement radial, il s’agit d’une pression appliquée sur tout le pourtour extérieur de la structure.
Le calcul du déplacement radial us en fonction de la pression appliquée Ps permet de tracer la courbe de
confinement sur un graphe identique à celui de la courbe de convergence.
On distingue également deux phases dans les modèles classiques :
• Une phase élastique, (de us = 0 à us = umax) Le soutènement se comporte linéairement.
• Une phase plastique, après umax. Cette zone correspond à la rupture du soutènement, elle est
donc interdite.

Fig. 2 – Courbes de confinement

2.2. 3. Optimisation du soutènement

Après l’obtention les deux courbes (l’une caractérisant le comportement du terrain et l’autre celui
du soutènement), on va pouvoir les coupler pour étudier l’interaction roche-structure. La pression fictive de
l’un correspond parfaitement au chargement du second. Mais le couplage va nécessiter l’introduction d’un
nouveau paramètre : le déplacement à la pose du soutènement.

6
 Fig. 2 – assemblage les de convergence et de confinement

En effet, le soutènement n’est pas posé immédiatement au front de taille. Il est posé à quelques
décimètres en arrière, alors que le terrain s’est déjà partiellement déconfiné. On ajoute ainsi un paramètre ud,
qui est bien entendu étroitement lié à λd, taux de déconfinement à la pose.
Le point d’intersection des deux courbes correspond ainsi au point d’équilibre entre terrain et
soutènement.C’est ce point (Péq, uéq) qui donne l’état mécanique de la structure "à l’infini", loin du front de
taille. Toute la puissance de la méthode convergence-confinement réside donc dans cette simplicité de
représentation. En jouant sur chacun des paramètres du problème, on optimise le soutènement : pas ou peu
de plasticité pour le terrain, et chargement à 70 ou 80 % de la rupture pour le soutènement
A titre d’exemple, en jouant sur le paramètre ud : un soutènement placé trop près du front de taille sera
chargé prématurément et arrivera donc plus rapidement à la rupture. A l’opposé, un soutènement placé
trop loin du front n’aura aucun effet, car le terrain se sera déjà presque entièrement déconfiné, voire
effondré,et le chargement sera pratiquement nul. Rappelons que le soutènement est aussi là pour limiter la
convergence.

4. Comportement mécanique des soutènements

Le soutènement s’oppose à la convergence des parois du tunnel en exerçant une pression qu’on appelle
couramment pression de soutènement (Ps). Cette dernière croit avec la rigidité du soutènement et est limitée
par sa résistance. Elle a l’expression :
𝑘𝑠
𝑃𝑠 = 𝑢
𝑅
𝑘𝑠 raideur de soutènement u : déplacement radial
P :pression exercée sur le soutènement R :rayon de l’excavation
3.1. Raideur et pression maximale d’un soutènement ou d’un revêtement
Pour les différents type de soutènement le raideur( module de rigidité ) et la pression maximale sont
déterminés par les formules suivantes :
a) Soutènement (ou revêtement) mince (e << R) en béton projeté
𝐸 .𝑒
- Raideur 𝑘𝑏 = (1−𝜈𝑏 2 )𝑅
𝑏
- Pression maximale du soutènement
Pression maximale qui peut développer l’anneau en béton
𝜎𝑏 𝑚𝑎𝑥
𝑃𝑏 𝑚𝑎𝑥 = .𝑒
𝑅
R : rayon de l’excavation e : l’épaisseur du soutènement (ou revêtement)
𝐸𝑏 , 𝜈𝑏 : le module et le coefficient de poisson du béton

7
𝜎𝑏 𝑚𝑎𝑥 : la contrainte limite admissible du béton
Il habituel de dimensionner le soutènement ou le revêtement de telle sorte que pour un béton à 25MPa la
sollicitation soit en fait limitée aux environs de 6 MPa

b) Soutènement par cintres calés à la paroi

𝐸𝑎 .𝑆
- Raideur 𝑘𝑐 = 𝑒.𝑅
- Pression maximale du soutènement
𝑆. 𝜎𝑎
𝑃𝑐 𝑚𝑎𝑥 =
𝑅. 𝑒
𝜎𝑎 :résistance de l’acier 𝐸𝑎 : le module du matériau
e : espacement longitudinal des cintre S : section du cintre rayon de l’excavation

c)-Soutènement par cintres boulon à ancrage ponctuel

1 𝑒𝑐 .𝑒𝑙 4𝑙
- Raideur = [𝜋.𝑑2 .𝐸 + 𝑄]
𝑘𝑠 𝑅 𝑎

- Pression maximale du soutènement

𝑇𝑏𝑟
𝑃𝑆 𝑚𝑎𝑥 =
𝑒𝑐 . 𝑒𝑙
𝑒𝑐 : espacement circonférentiel des boulons l :longueur des boulons.
𝑒𝑙 : espacement longitudinal des boulons d : diamètre des boulons
𝐸𝑎 : module d’Young du matériau constituant la tige du boulon..
Q : paramètre qui permet de prendre en compte les déformations qui se produisent au niveau de l’ancrage
et de tête des boulons. Il peut être déterminé à partir de la courbe charge –déformation résultant d’un
essai d’arrachage
𝑇𝑏𝑟 : charge résistante ultime du boulon
Exercice 02
On se place dans le cas d’une portion de tunnel routier à projeter de longueur 300m pour laquelle la
hauteur de couverture h varie selon le schéma indiqué Fig.1.

FIG. 1 – Profil en long du tunnel.

8
 On désire dimensionner le soutènement lors du creusement du tunnel. L’excavation est supposée
circulaire de rayon R égal à 5 mètres et creusé dans le sens d’une augmentation de la hauteur de couverture.
Le terrain est homogène dans la portion considérée et possède un comportement élastoplastique
parfait dont les paramètres sont donnés en Fig 1 On supposera que la méthode convergence-confinement
est valable et la pression initiale P est calculée de la façon suivante pour chaque section :P = γ.h
avec :
γ: poids volumique du terrain pris égal à 0,025 MN/m3 ;
h : hauteur de couverture de la section considérée ; la valeur de h sera prise égale à la distance
entre la côte du terrain naturel au droit de la section et le centre de l’excavation.

I. Dimensionnement de la section de hauteur de couverture minimale (Si)

On s’intéresse dans cette partie à la section (Si) (voir Fig. 1) de la portion dont la hauteur de couverture
est minimale. On cherche à dimensionner l’épaisseur de soutènement à mettre en œuvre.
1. Montrer que la pression initiale P du terrain avant creusement pour la section (Si) est égale à
1.25 MPa. Y a-t-il apparition de la plasticité dans le terrain pour la section (Si) ?
2. Quelle est la valeur finale du déplacement de la paroi au cas où aucun soutènement n’est mis en place ?
3. Tracer la courbe de convergence du massif pour la section (Si).
4. On décide de mettre en œuvre en (Si) un soutènement à base de béton projeté dont les caractéristiques
sont les suivantes :
– Module d’Young Eb = 10 000 MPa – Coefficient de Poisson νb = 0,2
– Contrainte maximale de compression admissible σmax = 20 MPa
L’épaisseur choisie (e) est de 20 centimètres et la distance de pose (d) du soutènement par rapport
au front de taille est de 1,5m.
(a) Calculer le taux de déconfinement _p à la pose du soutènement. Calculer la pression fictive Pid
à la pose du soutènement.
(b) Quel est le déplacement ud du terrain à la pose du soutènement ?
(c) Tracer la courbe de confinement.
5. On recherche les valeurs atteintes à l’équilibre :
(a) Quelle est la valeur du déplacement de la paroi à l’équilibre ?
(b) Quelle est la valeur du taux de déconfinement à l’équilibre ?
(c) Quel pourcentage de la contrainte maximale admissible _s a-t-on mobilisé à l’équilibre ?
(d) Conclusions sur l’épaisseur de soutènement choisie pour la section (Si).
II. Dimensionnement du soutènement lorsque la hauteur h varie
On s’intéresse dans cette partie à l’effet de la variation de hauteur h sur le dimensionnement du
soutènement.
1. A partir de quelle distance (x) du début de la portion étudiée (voir Figure 3.20) y aura-t-il apparition
de la plasticité dans le terrain ?
2. Le dimensionnement est réalisé en étudiant dans cette question la section (Sf) de hauteur h égale à
300m (voir Fig 1).
(a) Montrer que la pression initiale P du terrain est de 7,5 MPa. Calculer la pression d’apparition
de la plasticité à la paroi dans le cas où la section étudiée est (Sf).
(b) Tracer la courbe de convergence du terrain dans le cas où la section étudiée est (Sf).
(c) On suppose que la distance de pose (d) est telle que le taux de déconfinement à la pose est de
0,6. Déterminer graphiquement le déplacement de la paroi à la pose à partir de la courbe de
convergence. Est-ce que cette valeur de déplacement est admissible ?
(d) On propose de mettre en oeuvre un béton projeté don les caractéristiques sont identiques à la
question 2 de la partie précédente. On garde la même épaisseur e de soutènement égale à 20cm.
Tracer la courbe de confinement correspondante. Conclusion sur l’épaisseur de soutènement
choisie pour la section (Sf).