Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 1/62

Organisme(s) : EDF-DIS/SEPTEN, EDF-R&D/MMC, ENS Cachan
















Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
Document : U2.04.06




Comment creuser un tunnel : méthodologie
d'excavation



Résumé :

Cette note propose une méthodologie pour simuler le creusement d'une galerie souterraine avec Code_Aster.
La méthode de base est une méthode couramment utilisée dans ce genre d'études : la méthode « convergence
­ confinement ».

Après un rappel sur le principe de la méthode, les principales étapes du fichier de commandes Code_Aster sont
décrites. Différents exemples numériques permettent de valider la procédure.
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 2/62


Table
des
matières

1 Comment simuler le creusement d'un tunnel avec Code_Aster ? ........................................................3
1.1 Principe de la méthode, mise en oeuvre et validation.....................................................................3
2 Introduction ............................................................................................................................................4
3 Une méthode pour simuler le creusement d'une galerie à partir d'un modèle 2D : la méthode
convergence-confinement......................................................................................................................5
3.1 Principe général...............................................................................................................................5
3.2 Application de la méthode pour un calcul numérique par éléments finis ........................................7
4 Avant de s'attaquer au ficher de commande Code_Aster... .................................................................8
4.1 Comment définir les modèles à partir d'un maillage simple ? .........................................................9
4.2 Comment initialiser les contraintes ?.............................................................................................11
4.3 Comment calculer les réactions nodales au bord de la « future » galerie ? .................................12
4.4 Comment simuler la création d'un « vide » dans le massif et la pose du béton ? ........................13
4.4.1 Méthode A ............................................................................................................................14
4.4.2 Méthode B ............................................................................................................................14

4.5 Résumé des méthodes proposées................................................................................................14
5 Exemples de fichiers de commandes ..................................................................................................15
5.1 Le problème traité..........................................................................................................................15
5.2 Cas n° 1 : excavation sans soutènement avec initialisation des contraintes par un calcul et
« ramollissement » des éléments « excavés »..............................................................................16
5.3 Cas n°2 : excavation avec soutènement avec initialisation des contraintes par appel à
CREA_CHAMP et déconfinement suivant la méthode A.................................................................16
5.4 Cas n°3 : excavation avec soutènement avec initialisation des contraintes par appel à
CREA_CHAMP et déconfinement suivant la méthode B..................................................................16
6 Validation du Code_Aster sur un exemple d'excavation en milieu élastique linéaire .........................17
7 En guise de conclusion : conseils et perspectives ..............................................................................18
8 Bibliographie ........................................................................................................................................19
Annexe 1
Formules analytiques pour appliquer la méthode convergence-confinement au cas
d'un massif rocheux et d'un soutènement élastiques et linéaires...................................................20
Annexe 2
Organigramme de synthèse sur les méthodes permettant de simuler une
excavation dans le Code_Aster ......................................................................................................22
Annexe 3
Fichier de maillage réalisé avec GIBI ......................................................................23
Annexe 4
Excavation sans soutènement, sur la base d'un seul modèle (cas n°1). Fichier de
commandes Code_Aster.................................................................................................................26
Annexe 5
Excavation avec soutènement, méthode A (cas n°2). Fichier de commandes
Code_Aster .................................................................................................................................37
Annexe 6
Excavation avec soutènement, méthode B (cas n°3). Fichier de commandes
Code_Aster .................................................................................................................................48
Annexe 7
Comparaison des contraintes obtenues par le calcul numérique et par la solution
analytique .................................................................................................................................61

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 3/62


1 Comment simuler le creusement d'un tunnel avec
Code_Aster ?

1.1
Principe de la méthode, mise en oeuvre et validation

Contexte

Les études de géomécanique sont généralement basées sur une simulation de creusement de galerie
souterraine. Des exemples d'application peuvent être cités :

· évaluer la zone endommagée par l'excavation (EDZ) autour d'une galerie de stockage ;
· étudier la resaturation d'une alvéole de stockage par les eaux du site.

Un certain nombre d'études ont déjà été menées par le département AMA sur ce sujet, avec
Code_Aster. Cependant, peu d'éléments pratiques sont disponibles dans les documentations pour
reproduire ce type de calcul. Le département MMC a entrepris une telle modélisation avec
Code_Aster, afin de s'approprier la procédure d'application de la méthode classiquement utilisée pour
ce genre de calcul : la méthode « convergence ­confinement ». Il ressort de cette expérience que
cette application n'est pas complètement triviale qu'il est nécessaire de se poser quelques questions
techniques pratiques de mise en oeuvre. Capitaliser cette expérience pour les futurs utilisateurs est
apparue comme assez important, dans l'intérêt collectif des études sur le stockage notamment.

Objectif

Cette note a pour objectif principal de fournir quelques conseils techniques préliminaires aux
utilisateurs de Code_Aster souhaitant modéliser une excavation souterraine.

Méthodologie

Cette note présente une application à un fichier de commande de Code_Aster de la méthode
convergence ­ confinement. Après un bref rappel sur le principe de la méthode, une description
pratique et opérationnelle des commandes à utiliser est donnée. La méthode est illustrée par des
calculs de validation de Code_Aster, dont les fichiers de commande sont fournis en annexe.

Résultat

Grâce à la mise en application du protocole proposé, deux calculs de validation de Code_Aster ont été
mis en oeuvre. L'écart relatif entre résultats numériques et solution analytique est inférieur à 2%.

Perspectives

La méthode peut être élargie aux calculs non linéaires (plasticité, endommagement) et couplés en
THM, notamment dans le cadre d'études destinées au stockage de déchets nucléaires.

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 4/62


2 Introduction

Depuis plusieurs années, des études sont réalisées avec Code_Aster afin de modéliser le
comportement d'ouvrages géotechniques (barrages en terre, tunnels, barrières ouvragées pour le
stockage de déchets...).

Code_Aster a déjà été utilisé notamment pour simuler le creusement de galeries ou de puits, dans le
cadre du projet Stockage géologique de déchets nucléaires HAVL (T4-01-10) ou lors d'études
antérieures sur le stockage profond. Les rapports rédigés jusqu'à présent (par exemple [bib6], [bib7] ou
[bib4]) se focalisent naturellement sur les résultats, afin de répondre à la question technique précise
qui a motivé l'étude. Or simuler une excavation à l'aide d'un code éléments finis n'est pas forcément
une chose facile, et même si les principes généraux sont rappelés dans les documents cités plus haut,
on trouve finalement peu d'éléments sur la structure des fichiers de commandes qui ont servi de
support aux calculs.

Afin d'aider les ingénieurs en charge des futures études d'excavation souterraine avec Code_Aster,
cette note indique quelques conseils pratiques pour débuter dans la réalisation de ce type de calcul.
En effet, dans le cadre du projet Stockage, MMC a décidé de s'approprier complètement la démarche
mise en oeuvre par AMA en 2000 et en 2001. Pour y parvenir, toute la démarche a été reproduite avec
la version 6 de Code_Aster, sur la base d'un nouveau maillage et en explorant quelques variantes.
MMC a également bénéficié de l'assistance des agents d'AMA. Par ailleurs, ce travail a abouti à une
validation de Code_Aster d'après des formules analytiques classiques en élasticité linéaire (formules
de Kirsch et méthode convergence-confinement, [bib5]).

Ce rapport présente donc :

· la méthode classique de simulation d'une excavation souterraine en 2D au moyen d'un code
éléments finis ;
· les différentes option disponibles pour appliquer cette méthode avec Code_Aster ;
· deux cas-tests de validation de Code_Aster pour les problèmes d'excavations souterraines.

Le pré-requis à une lecture profitable de cette note est la formation de base à l'utilisation de
Code_Aster ainsi qu'un minimum de familiarisation au progiciel. Le détail des différentes commandes
utilisées est donné par la documentation Utilisateur de Code_Aster (http://www.code-aster.org).


Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 5/62


3
Une méthode pour simuler le creusement d'une galerie à
partir d'un modèle 2D

: la méthode convergence-
confinement

3.1 Principe
général

Cette partie s'inspire largement de [bib5]. Signalons que le CIH et TEGG ont également mené un
certain nombre d'études avec cette méthode (par exemple, [bib2]). Il est conseillé au lecteur de se
reporter à ces documents pour plus d'informations sur le principe de la méthode. Les paragraphes qui
suivent ne résument que l'essentiel de la démarche.

La méthode convergence-confinement est couramment utilisée dans l'ingénierie des ouvrages
souterrains. Son objectif est d'obtenir un ordre de grandeur des déplacements des parois du tunnel
ainsi que les efforts repris par la roche et le soutènement. Cette méthode permet de simplifier le calcul
d'un ouvrage tridimensionnel par un calcul bidimensionnel, par l'introduction d'un paramètre
adimensionnel appelé `taux de déconfinement'. Elle repose sur les hypothèses suivantes :

· déformations planes avec hypothèse de petites perturbations;
· le tunnel est supposé de section circulaire et d'axe horizontal ;
· terrain homogène d'extension infinie ;
· massif suivant un comportement élastique linéaire ou élasto-plastique ;
· état initial des contraintes supposé isotrope et homogène ;
· tunnel profond : pas de variation de contraintes significative sur la hauteur de la galerie. En
pratique, si H est la profondeur moyenne de l'ouvrage et R son rayon, cette hypothèse est
supposée satisfaite si H/R>10 ;
· équilibre quasi-statique (pas de termes d'accélération).

On s'intéresse à une section située dans un plan perpendiculaire à l'axe du tunnel et on souhaite
mener un calcul bidimensionnel. Le paramètre est censé prendre en compte l'influence mécanique
de la proximité du front de taille à cette section, c'est-à-dire d'un phénomène dont l'origine se situe
hors du plan considéré par le calcul. dépend de plusieurs paramètres (roche, soutènement, longueur
de tunnel non soutenu derrière le front de taille...) et sa détermination n'est pas forcément immédiate
(nombreuses publications sur le sujet, par exemple [bib1]). Ce problème de détermination analytique
du taux de déconfinement sort du cadre de ce document.

En fait, on introduit pour considérer un tenseur des contraintes fictif dans le terrain, qui est une
fraction de la contrainte initiale 0
:

= (1- ) 0
. avec 0 1
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 6/62


La [Figure 3.1-a] ci-dessous illustre l'évolution de et de la contrainte radiale R pour un tunnel non
soutenu.


=0
=1
0>>1


R=0
R=0
R=(1-).0

Figure 3.1-a : Evolution du taux de déconfinement et de la contrainte radiale R

dans le cas d'un tunnel non soutenu

Remarquons que = 1 correspond au déconfinement total de la roche : l'influence du front de taille
sur le comportement de la tranche de tunnel a disparu et le tunnel est assimilable à un tube très épais.

Puisqu'une partie, voire la totalité des contraintes initialement présentes au sein du massif
disparaissent (c'est précisément le phénomène de déconfinement), les parois de l'excavation vont
avoir tendance à se rapprocher pour atteindre un nouvel équilibre mécanique. C'est le phénomène de
« convergence ». Ce phénomène peut aboutir à la ruine de l'ouvrage si la structure ne parvient pas à
retrouver un état d'équilibre stable suite à l'excavation.

Si, pour des raisons de sécurité ou de stabilité, on décide de poser un soutènement ou un revêtement
à la paroi du tunnel, ceux-ci vont, de par leur raideur mécanique, s'opposer au phénomène naturel de
convergence. Dans ce cas, l'équilibre final dépend donc de l'interaction mécanique entre la roche et le
revêtement. D'une façon générale, cet équilibre ne permet pas aux contraintes dans le massif rocheux
de s'annuler comme dans le cas du tunnel non soutenu. On dit alors que le terrain est confiné, d'où le
nom de la méthode « convergence-confinement ».
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 7/62


Graphiquement, l'application de cette méthode revient à rechercher le point d'intersection de la courbe
de convergence, déduite du comportement du massif, et de la courbe de confinement, déduite du
comportement du soutènement [Figure 3.1-b].


Figure 3.1-b : Exemple de courbes de convergence et de confinement

Les équations de la méthode « convergence-confinement » dans le cas d'un massif élastique linéaire
sont fournies en [§Annexe 1].

Que ce soit pour des calculs analytiques ou numériques, cette méthode permet, à l'aide d'un simple
modèle 2D, de traiter le problème 3D que constitue la simulation d'une excavation.


3.2 Application de la méthode pour un calcul numérique par éléments
finis

Une particularité des calculs d'excavation par éléments finis est la nécessité de mettre en oeuvre
plusieurs modèles (au sens large).

En effet, un déroulement classique de la modélisation peut se résumer par les étapes suivantes :

· étape 1: initialisation des contraintes in situ ;
· étape 2 : calcul des réactions nodales au niveau des parois de l'excavation ;
· étape 3 : déconfinement du massif pour simuler l'excavation progressive et l'éloignement du
front de taille ;
· étape 4 : éventuelle pose d'un soutènement / revêtement et fin du déconfinement.

Si l'étude le demande, l'enchaînement des étapes 2, 3 et 4 peut être répété (cas d'une excavation en
sections divisées, par exemple).
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 8/62


Dans la plupart des cas, l'enchaînement des calculs se fait donc sur la base de quatre configurations
[Figure 3.2-a].


Massif de sol
Massif de sol
Excavation de la
Pose du
Initialisation des
galerie
revêtement béton
contraintes
Calcul des réactions
nodales
1
2
3
4

Figure 1.2-a : Exemple d'enchaînement typique d'un calcul d'excavation avec un code de calcul

La première configuration sert à :

· initialiser les contraintes d'origine géostatique ;
· initialiser la pression hydrostatique due à la présence éventuelle d'eau et la température (la
présente note ne traite pas ce point précis en détail) ;

La deuxième configuration permet de calculer les réactions aux noeuds représentant le bord de
l'excavation.

A ces stades de la modélisation, tous les éléments du maillage correspondent donc à un matériau de
type sol ou roche. On obtient donc un massif de sol dans lequel règne un état de contraintes
correspond à l'état de contraintes in situ dans le plan perpendiculaire à l'axe de la galerie. On connaît
également les réactions nodales au bord de l'excavation, ce qui va permettre un déconfinement partiel
ou total du massif dans les étapes qui suivent.

La troisième configuration est dédiée au déconfinement : on diminue les réactions nodales au bord de
l'excavation afin de simuler le creusement du tunnel. Lors de la réalisation de cette étape, les éléments
finis dans la région correspondant à l'intérieur de la galerie ne doivent plus participer à la rigidité du
modèle. Comme on le verra plus loin, ceci peut être pris en compte de plusieurs façons en pratique.

On passe éventuellement à une quatrième étape si l'on veut simuler la pose d'un soutènement béton
en cours de déconfinement par exemple. Dans ce cas, on rajoute des éléments avec des
caractéristiques de béton et on poursuit la diminution des réactions nodales calculées dans l'étape n°1
pour achever le calcul.

On remarque donc que certaines parties du modèle initial vont se voir affecter successivement des
propriétés de sol, de « vide » puis de béton. Dans cet enchaînement se situe la source de quelques
difficultés intrinsèques à ce genre de calcul.

L'application de cette démarche au moyen de Code_Aster fait l'objet des chapitres suivants. Elle est
basée sur un cas simple.


4
Avant de s'attaquer au ficher de commande Code_Aster...

Ce chapitre concerne quelques points de modélisation particuliers qu'il semble important de
commenter avant de s'intéresser aux fichiers de commandes proprement dits. Elle est composée d'une
suite de paragraphes traitant chacun une question qu'un ingénieur peut se poser lorsqu'il mène un
calcul classique d'excavation à l'aide d'un code type éléments finis comme Code_Aster.
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 9/62


4.1
Comment définir les modèles à partir d'un maillage simple ?

Le maillage choisi dans cette étude représente un quart de modèle représentant une galerie
cylindrique en milieu infini. Le rayon de la galerie est de 1,50 mètre, l'épaisseur de béton est 0,30
mètre et le maillage est un carré de 20 mètres de côté. D'après les règles usuelles de modélisation, le
rapport entre le rayon excavé R et la dimension caractéristique du maillage L est suffisant pour
considérer que les conditions aux limites ne perturbent pas le comportement de l'excavation
(L 10 x R).
Sol

Sol, vide
ou béton
Sol ou vide

Figure 2.1-a : Maillage utilisé et matériaux
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 10/62


Du point de vue des modèles (au sens Code_Aster), il faut distinguer quelques zones particulières du
maillage (en plus des autres zones plus classiques, comme les bords du maillage) et créer les objets
suivants (les noms font référence aux fichiers de commande présentés en Annexe) :

· le bord excavé où seront appliquées les réactions nodales pour simuler le déconfinement
(appelé BORD) ;
· les deux points qui sont situés aux extrémités de cette courbe, qui sont concernés à la fois par
le chargement de déconfinement et par les conditions aux limites au bord du massif.


Figure 4.1-b : Points et ensemble de points particuliers à identifier dans les modèles Code_Aster

On peut donc définir (par exemple, car plusieurs configurations sont possibles) :

· un modèle SOL, dans lequel tout le maillage est affecté d'éléments finis ;
· un modèle SOL_REST qui ne comprend pas les mailles qui correspondent à la partie excavée
(elles ne sont pas affectées d'éléments finis) ;
· un modèle SOL_REST0 qui comprend SOL_REST et les mailles correspondant au revêtement
en béton affectées d'éléments finis.

Remarque :

L'utilisation d'une telle géométrie pour faire un calcul réel d'excavation est partiellement
critiquable, car la symétrie proposée risque d'engendrer un chargement non physique. En cas
d'application du poids propre par exemple, celui-ci serait dirigé vers le haut dans la partie
inférieure du tunnel !



Pesanteur
Partie maillée
Tunnel
Pesanteur induite par les
conditions aux limites
Partie non maillée mais
simulée par symétrie

Figure 4.1-c : Exemple d'aberration que peut engendrer l'utilisation d'un quart de modèle
dans la simulation d'un tunnel
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 11/62


Pour une étude réaliste où l'on souhaiterait initialiser les contraintes à l'aide d'un chargement en poids
propre, il serait donc nécessaire de « mailler le terrain » jusqu'au substratum rigide (roche considérée
indéformable), ou tout du moins jusqu'à une profondeur suffisante pour s'affranchir du problème
évoqué ci-dessus. On maille donc un demi-modèle dans ce cas là.

Cependant, dans le cadre de notre étude, cette contrainte de maillage ne nous a pas gênés, puisque
nous n'avons pris en compte ni le poids propre, ni les couplages THM. Le chargement simulé est tout à
fait compatible avec les solutions analytiques testées.

4.2
Comment initialiser les contraintes ?

Les contraintes in situ sont généralement représentées par un tenseur d'ordre 2 dont les directions
principales correspondent à la verticale et à l'horizontale. La contrainte verticale est généralement
égale au poids des différentes formations situées au-dessus du point considéré et la contrainte
horizontale est proportionnelle à la contrainte verticale :

v = .z


h = K 0
. v

avec le poids volumique du terrain sus-jacent (en kN/m3 par exemple) et K0 un coefficient sans
dimension. K0 peut être déterminé par des mesures in situ ou estimé par des relations plus ou moins
empiriques. Dans le cas d'un massif semi-infini soumis à une contrainte externe sur son bord supérieur
ou à son poids propre, la théorie de l'élasticité linéaire fournit une valeur de K0 en fonction du
coefficient de Poisson :

K0 =

1-

Deux méthodes ont été testées avec le Code_Aster pour initialiser les contraintes dans le terrain
encaissant :

· réalisation d'un calcul (commande STAT_NON_LINE) avec un matériau fictif doté d'un
coefficient de Poisson permettant d'obtenir le rapport K0 souhaité. Ce calcul est réalisé sur le
modèle qui reprend tout le maillage de l'étude (par exemple, le modèle dénommé SOL dans le
chapitre précédent). Dans ce cas, K0 1 (cas de l'élasticité linéaire). Il existe de nombreux
cas où K0 1 (si le sol est soumis à des contraintes tectoniques, par exemple). Dans ce cas,
la méthode suivante devient obligatoire ;
· affecter directement les contraintes à tous les éléments du maillage par la commande
CREA_CHAMP (option : OPERATION = `AFFE') ;

La première solution nécessite de définir un matériau fictif et de mettre en oeuvre un calcul de plus.
Cependant, si le chargement est le poids propre (ce qui n'est pas le cas du cas-test que nous
proposons), cette méthode nous a paru à la fois intuitive et simple. Dans le cas d'un champ de
contraintes uniforme, l'utilisation de CREA_CHAMP est indéniablement la méthode la plus intéressante :
elle économise du temps de calcul et son appel est encore plus simple. Pour des distributions de
contraintes plus complexes, CREA_CHAMP fonctionne également mais nous ne l'avons pas utilisé (le
paragraphe [§3.5.3.1] de la documentation [U4.72.04] indice B1 donne un exemple adaptable à notre
problème).
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 12/62


4.3 Comment calculer les réactions nodales au bord de la « future »
galerie ?

Pour calculer les réactions nodales au bord de la galerie, il est nécessaire de bloquer cette partie du
maillage. Cette opération ne doit pas engendrer de contraintes non physiques par incompatibilité avec
le chargement appliqué lors de la phase d'initialisation. des contraintes. Une possibilité offerte consiste
à imposer le même chargement que lors de l'initialisation des contraintes en bloquant les noeuds du
bord de la galerie seulement pendant cette étape de calcul. Cette opération est sans effet sur le
résultat global, qui reste identique à celui de l'étape précédente, mais le blocage « temporaire » des
noeuds du bord de la galerie permet d'y évaluer les réactions nodales.


Même chargement
qu'à l'étape
précédente
Noeuds bloqués
u = 0

Figure 4.3-a : Blocage des noeuds du bord de la galerie pour y calculer
les réactions nodales

Concrètement, ce blocage relatif du bord de la galerie est possible grâce à l'option DIDI (pour DIrichlet
DIfférentiel) du mot-clé EXCIT de l'opérateur STAT_NON_LINE (Doc. Aster [U4.51.03] indice F4,
paragraphe [§3.2.2]). Le blocage de ces noeuds ne s'applique que sur l'incrément de déplacement
considéré et non sur le déplacement total (on impose u = 0 et non u = 0).

L'état initial de ce calcul (mot-clé ETAT_INIT de l'opérateur STAT_NON_LINE) est défini par le champ
de contraintes obtenu à l'issue de l'étape précédente.

Une fois ce calcul intermédiaire effectué, le calcul des réactions nodales est effectué simplement par
l'appel à une commande CALC_NO munie de l'option OPTION = `REAC_NODA'. Il convient alors de
fournir à la commande CALC_NO tous les chargements ayant produit le résultat à partir duquel on
calcule les réactions nodales, sans omettre les chargement volumiques s'ils existent (non pris en
compte dans les exemples traités ici).

On construit alors un vecteur de chargement par la récupération des réactions nodales (CREA_CHAMP
avec les mots-clés TYPE_CHAM = 'NOEU_DEPL_R' , NOM_CHAM = 'REAC_NODA' et OPERATION =
'EXTR'). Il convient de noter que d'après le paragraphe [§3.1.1] de la documentation Utilisateur de
Code_Aster [U4.72.04] indice B1, l'option TYPE_CHAM = `NOEU_DEPL_R' de la commande
CREA_CHAMP est en fait sans effet ici (mais néanmoins obligatoire du point de vue syntaxique),
puisqu'on réalise une extraction. Ce vecteur est alors défini par la commande AFFE_CHAR_MECA avec
le mot-clé VECT_ASSE en tant que chargement pour l'appel suivant à la commande STAT_NON_LINE
(correspondant à l'excavation progressive de la galerie). Ce chargement est associé à une fonction
(opérateur DEFI_FONCTION) décrivant l'évolution du taux de déconfinement au fur et à mesure de la
progression du creusement.

Remarquons également que toutes les réactions nodales sont extraites : celles qui agissent sur le bord
de la galerie comme celles qui agissent sur les autres bords du modèle. Etant donné que ces dernières
agissent sur des points bloqués à toutes les étapes du calcul d'excavation, leur injection en tant que
chargement dans le STAT_NON_LINE suivant est sans effet sur les contraintes et les déformations au
sein de la structure.
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 13/62


4.4
Comment simuler la création d'un « vide » dans le massif et la pose
du béton ?


Une fois effectué le calcul des réactions nodales se pose la question de l' « élimination » de la partie
excavée du modèle numérique afin que sa rigidité n'entrave pas la convergence du tunnel. Pour y
parvenir, nous avons retenu deux méthodes [Figure 4.4-a] :

· méthode A : quasi-annulation des propriétés mécaniques des éléments situés dans la zone
excavée (exemple : E = 0,0001 Pa), puis introduction de propriétés plus réalistes lors de la
pose du soutènement ou du revêtement. Cette méthode permet de simplifier le fichier de
commande Code_Aster et donne des résultats corrects pour le cas simple que nous avons
étudié (petite galerie circulaire, excavée en une seule section dans un massif élastique). Pour
mener des études plus élaborées où le traitement numérique pourrait être affecté par la
présence d'élément à très faible rigidité, il nous semble néanmoins préférable de s'appuyer sur
la méthode suivante ;
· méthode B : initialisation des contraintes directement par création de champs aux points de
Gauss issus d'un calcul concernant une étape précédente.



Méthode A
Méthode B
Terrain
« Vide »
Béton
1
2
3

Figure 4.4-a : Différents principes de modélisation pour simuler le déconfinement d'un massif

D'autres méthodes que nous n'avons pas testées peuvent sans doute être appliquées au problème
étudié (comme la création de noeuds doubles aux frontières entre matériaux qui permettent de lier ou
pas les deux structures).
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 14/62


4.4.1 Méthode
A

La méthode A n'appelle pas de commentaire particulier : il suffit d'affecter des caractéristiques
mécaniques très faibles aux mailles devenant « vides » lors du calcul du déconfinement.

On procède en deux temps :
· une première commande STAT_NON_LINE, qui permet de réinjecter le chargement composé
du vecteur des réactions nodales et des conditions aux limites. Les mailles « vides »
correspondent donc à un matériau très mou ;
· un deuxième appel à STAT_NON_LINE qui introduit le soutènement ou le revêtement en béton
en affectant aux mailles correspondantes des caractéristiques réalistes pour un tel matériau.

A chaque appel, l'initialisation du calcul reprend l'intégralité des champs issus des calculs précédents
(opérande EVOL_NOLI pour le mot-clé ETAT_INIT).

4.4.2 Méthode
B

Cette procédure est basée sur l'enchaînement de plusieurs modèles (au sens Code_Aster). Le calcul
s'effectue en copiant certains champs d'un modèle à l'autre.

Les champs à affecter au modèle correspondant à l'étape B-3 de la [Figure 4.4-a] sont formellement la
combinaison linéaire de deux champs :

· les champs issus de l'étape de calcul précédente (B-2) et qui ne concerne que le modèle
correspondant au massif de sol privé de la zone excavée ;
· les champs affectés aux éléments du groupe de maille qui représentent les voussoirs en
béton, dans le modèle qui comprend le massif et le revêtement de la galerie. Dans notre cas,
ces champs doivent être initialisés à 0 dans B-3. Pour cela, on peut par exemple affecter un
poids nul à leur contribution dans la combinaison linéaire. Ainsi ces champs peuvent en fait
être obtenu par un calcul intermédiaire sans réelle signification physique, par exemple
l'application simple des conditions aux limites.

On utilise la commande CREA_CHAMP avec l'option ASSE pour affecter aux points de Gauss du
troisième modèle la combinaison linéaire de champs issus de calculs précédents.

4.5
Résumé des méthodes proposées

Pour initialiser les contraintes, on peut faire appel à deux méthodes :

· Méthode I : faire un calcul (appel à STAT_NON_LINE) sur matériau fictif ;
· Méthode II : créer le champ de contraintes souhaité par CREA_CHAMP.

Pour simuler le creusement et la pose des voussoirs, on a le choix entre :

· Méthode A : qui consiste à affecter des caractéristiques mécaniquement très « souples » dans
la zone excavée ;
· Méthode B : qui recourt à l'utilisation de plusieurs modèles qui s'enchaînent et qui est plus
proche de la réalité physique de la structure modélisée, les matériaux apparaissant et
disparaissant par activation d'un modèle à l'autre.

Un organigramme synthétique est proposé en [§Annexe2].
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 15/62


5
Exemples de fichiers de commandes

Cette partie présente des exemples de structures de fichier de commandes Code_Aster concernant
une excavation circulaire en milieu infini et élastique linéaire, dans le cadre d'une étude purement
mécanique (pas de couplage THM).

Trois cas de calcul sont présentés dans cette partie :

· une excavation sans soutènement avec initialisation par un calcul portant sur un matériau fictif
pour obtenir le champ de contraintes souhaité (méthode I) ;
· une excavation avec soutènement, initialisation des contraintes par un appel à CREA_CHAMP et
suivi de la méthode A pour le déconfinement et la pose des voussoirs (méthodes II + A) ;
· une excavation avec soutènement, initialisation des contraintes par un appel à CREA_CHAMP et
suivi de la méthode B pour le déconfinement et la pose des voussoirs (méthodes II + B).

Pour les cas 2 et 3, le scénario de creusement est le suivant : excavation, déconfinement à 50%
( = 0,5), pose des voussoirs de 30 cm d'épaisseur et fin du déconfinement. Ces deux cas font l'objet
de cas-test de validation de Code_Aster (mise en oeuvre prévue pour début 2003).

5.1
Le problème traité

La géométrie du maillage est rappelé au paragraphe [§4.1]. Il contient 8477 noeuds et 3304 éléments.
Le rayon de la galerie est de 1,50 mètre, l'épaisseur de béton est 0,30 mètre et le maillage est un carré
de 20 mètres de coté. Les autres données sont résumées dans le tableau suivant.

Matériau Paramètre
Valeur

K0
1
Roche
v = h
5 MPa
E
4
GPa


0,3
Béton
E 20
GPa


0,2
Tableau 5.1-1 :Données des cas tests proposés

Les conditions aux limites et le chargement sont illustrés par la figure suivante :


Pression des terrains sus jacents (5 MPa)
Ux = 0
Ux = 0
Déconfinement
Uy = 0

Figure 5.1-a : Conditions aux limites et chargement imposés
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 16/62


.
2 G
A la fin du processus de déconfinement, =
.U R = 69
,
0
.
.
0 R

Un exemple de fichier de maillage en langage gibiane (mail.dgib) est présenté en [§Annexe 3].

5.2 Cas n° 1 : excavation sans soutènement avec initialisation des
contraintes par un calcul et « ramollissement » des éléments
« excavés »


Cet exemple est relativement simple : il s'agit de simuler une excavation sans pose du soutènement,
avec déconfinement total au bord de la galerie. On n'utilise donc qu'un seul modèle pour tout le calcul.

L'état initial est engendré par un calcul (STAT_NON_LINE) qui porte sur l'ensemble du maillage. Les
propriétés des éléments sont affectées en fonction de l'état de contraintes que l'on veut atteindre (ici
K0 = 1 donc = 0,4999, la valeur de 0,5 signifiant l'incompressibilité de la roche ne pouvant être
utilisée).

Le calcul suivant concerne les réactions nodales au bord de la future galerie. Il est initialisé à partir des
contraintes issues du premier appel à STAT_NON_LINE.

Le dernier appel à STAT_NON_LINE sert à réinjecter les réactions nodales dans un modèle où les
propriétés mécaniques des éléments excavés ont été très fortement affaiblies (E tend vers 0.). On
déconfine alors complètement le terrain en faisant tendre ces réactions vers 0.

Le fichier de commande correspondant est présenté en [§Annexe 4].

5.3 Cas n°2 : excavation avec soutènement avec initialisation des
contraintes par appel à CREA_CHAMP et déconfinement suivant la
méthode A


On suit le scénario d'excavation décrit plus haut. On utilise qu'un seul modèle pour tout le calcul. Une
commande STAT_NON_LINE supplémentaire permet d'introduire les voussoirs avec une rigidité
réaliste après un déconfinement de 50%.

Le fichier de commande correspondant est présenté en [§Annexe 5].

5.4 Cas n°3 : excavation avec soutènement avec initialisation des
contraintes par appel à CREA_CHAMP et déconfinement suivant la
méthode B


On suit toujours le scénario d'excavation décrit plus haut. Cette fois, trois modèles sont utilisés et un
calcul intermédiaire (sans réalité physique, appelé « bidon ») est nécessaire pour transférer les
champs de variables d'un modèle à l'autre au moment de la mise en place des voussoirs après
déconfinement de 50%.

Le fichier de commande correspondant est présenté en [§Annexe 6].
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 17/62


6 Validation
du
Code_Aster sur un exemple d'excavation en
milieu élastique linéaire

La validation de Code_Aster repose sur la comparaison des résultats numériques issus des cas n°1, 2
et 3 listés ci-dessus à la solution analytique de [§Annexe 1]. Pour chaque calcul, on présente les
résultats obtenus au niveau de la clé de voûte et du piédroit de la galerie, en terme de contrainte
radiale R , contrainte orthoradiale et déplacement radial U R ([Tableau 6-1], [Tableau 6-2] et
Tableau 6-3]). L' [§Annexe 7] présente deux graphiques décrivant l'évolution spatiale des contraintes le
long de l'axe vertical du modèle, au droit de l'excavation. Le bon accord entre solution analytique et
résultats numériques fait que l'écart entre ces courbes est à peine visible.



R(y)

A

R
U R
B

Figure 6-a : Grandeurs comparées pour la validation du Code_Aster



Point A
Point B
Variable Analytique
Aster Ecart
relatif Analytique
Aster
Ecart
relatif
(Pa)
0.
-8.411 E3
On vérifie
0. -1.625
E4
On vérifie que
R
que
| | << | |
R

| | << | |
R

(Pa)
- 1. E7
-9.883 E6 1,2 %
- 1. E7
-1.011 E7
1,1 %
Ur (m)
- 0.0024375 - 0.0024772
1,7 %
- 0.0024375 - 0.0023982
1,6 %
Tableau 6-1 : Cas n°1, comparaison solution analytique / résultats Code_Aster pour
les contraintes radiale et orthoradiale et pour le déplacement radial en A et en B



Point A
Point B
Variable Analytique
Aster Ecart
relatif Analytique
Aster
Ecart
relatif
(Pa)
- 1.52821
-1.52974 E6
0,1%
- 1.52821
-1.52652 E6
0,1 %
R
E6
E6

-8.40987 E6
0,7 %
- 8.47179
-8.52586 E6 0,6 %
(Pa)
- 8.47179
E6
E6
Ur (m)
- 0.0016925 - 0.0017218
1,7 %
- 0.0016925 - 0.0016664
1,5 %
Tableau 6-2 : Cas n°2, comparaison solution analytique / résultats Code_Aster pour
les contraintes radiale et orthoradiale et pour le déplacement radial en A et en B
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 18/62



Point A
Point B
Variable Analytique
Aster Ecart
relatif Analytique
Aster
Ecart
relatif
(Pa)
- 1.52821
- 1.52943
0,1 %
- 1.52821
- 1.53171
0,2 %
R
E6
E6
E6
E6

-8.40822 E6
0,8 %
- 8.47179
- 8.52418
0,6 %
(Pa)
- 8.47179
E6
E6
E6
Ur (m)
- 0.0016925 - 0.0017211
1,7 %
- 0.0016925 - 0.0016658
1,6 %
Tableau 6-3 : Cas n°3, comparaison solution analytique / résultats Code_Aster pour
les contraintes radiale et orthoradiale et pour le déplacement radial en A et en B


L'écart maximal entre résultats analytiques et numériques est inférieur à 2%, à part pour la contrainte
radiale au bord de la galerie excavée dans le cas n°1, où la valeur théorique est 0. La validité du calcul
est vérifiée en considérant que la contrainte radiale est bien négligeable devant la contrainte
orthoradiale.

Bien entendu, tous ces écarts peuvent être réduits si on raffine encore le maillage.



7
En guise de conclusion : conseils et perspectives

Cette note propose une méthodologie qui permet de mener des calculs d'excavation à l'aide de
Code_Aster. Plusieurs scénarii d'excavation sont passés en revue et plusieurs méthodes sont
proposées.

La méthode et le progiciel sont validés dans le cas d'une galerie circulaire, creusée dans un massif
infini constitué par un matériau élastique linéaire. Code_Aster reproduit de façon tout à fait
satisfaisante le comportement d'une telle structure souterraine, avec ou sans prise en compte du
soutènement et/ou du revêtement.

Du point de vue de l'utilisateur, il semble plus pratique et plus rapide d'initialiser les contraintes par un
appel à la commande CREA_CHAMP plutôt que par un calcul sur matériau fictif.

Si on cherche à modéliser un comportement purement mécanique et si le phasage de l'excavation est
relativement simple, travailler avec un seul modèle paraît être la méthode la plus aisée. Il suffit
d'affecter des propriétés de matériaux très faibles aux mailles devenant « vides ». Dans les cas plus
compliqués, la mise en oeuvre de plusieurs modèles utilisés successivement peut s'avérer plus fiable
du point de vue de la mise en oeuvre pratique (erreur de modélisation) et du point de vue numérique
(erreur de calcul), malgré les procédures intermédiaires de transfert des champs (contraintes,
déplacements, pressions, températures, variables internes...) d'un modèle à l'autre.

Une étape ultérieure de validation du Code_Aster pourra se faire sur les problèmes couplés linéaires
(THM en milieu saturé et élastique) ou couplés et/ou non linéaires (modèle CJS 1 s'apparentant au
modèle de Mohr Coulomb, excavation à court terme en non drainé à comparer avec [bib3]).
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 19/62


8 Bibliographie

[1]
D. BERNAUD et G. ROUSSET : La « nouvelle méthode implicite » pour l'étude du
dimensionnement des tunnels, Revue Française de Géotechnique n°60, pp 5-26, - 1992
[2]
P. CATEL : Aval du Cycle ­ site de Bure ­ Fiche 13 ­ Méthode convergence-confinement,
note EDF TEGG EFT GG/00 168 A ­ 2000
[3]
A. GIRAUD
: Couplages Thermo-Hydro-Mécaniques dans les milieux poreux peu
perméables : application aux argiles profondes, thèse de l'ENPC ­ 1993
[4]
D. LE BOULCH : Comparaison des modélisations THM 3D et 2D d'un ouvrage de stockage
avec le Code_Aster, rapport Ajilon Technologies Cénergys 01-A ­ 2002
[5]
M. PANET : Le calcul des tunnels par la méthode convergence-confinement, Presses de
l'ENPC ­ 1995
[6]
N. SELLALI, C. CHAVANT et G. DEBRUYNE : Modélisation hydroplastique de l'excavation
d'une galerie souterraine avec le Code_Aster, note EDF MMN HI-74/00/009/A ­ 2000
[7]
N. SELLALI, C. CHAVANT et G. DEBRUYNE : Modélisation THM d'un ouvrage souterrain de
stockage avec le Code_Aster, note EDF MMN HI-74/01/014/A ­ 2001

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 20/62


Annexe
1 Formules analytiques pour appliquer la méthode
convergence-confinement au cas d'un massif rocheux
et d'un soutènement élastiques et linéaires


Le milieu est supposé élastique linéaire isotrope et soumis à un champ de contraintes initial également isotrope
( K0 = 1).

Contrainte radiale, contrainte orthoradiale et déplacement radial à la paroi du tunnel en milieu élastique
soumis à un taux de déconfinement




R
. 2

0
R = 1-
.


2

r




R
. 2
= 1+
0


2 .


r



R2
0

U R =


r
G
2

G est donné par la relation suivante :
=
E
G

2 1
( + )

Comportement du soutènement :

Soit K s la raideur du soutènement, elle est donnée par la relation suivante si on considère que le soutènement
est assimilable à un tube épais ou mince ( vb est le coefficient de Poisson du béton) :

Eb e

si R > 10 e
(1- 2
b ) R
Ks =
2
2


Eb (Re - Ri )
si R 10 e

2
2
(1+
1 2
b ) (
[ - b)Re +Ri ]

Soit s
P la pression de confinement définie sur la figure suivante


Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 21/62


On a donc :
P R = e
s
b

K
Si k
s
s =
représente la rigidité relative et le taux de déconfinement à la mise en place du
2 G
d
soutènement, alors la pression de soutènement et le déplacement radial en paroi sont donnés par :


ks
0
Ps =
(1- d )

1+ ks



1+
0
d ks
U R =

R

1+ ks
2 G
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 22/62


Annexe
2
Organigramme de synthèse sur les méthodes
permettant de simuler une excavation dans le
Code_Aster


Notations

Les noms des objets sont ceux des fichiers de commande présentées dans les annexes suivantes.

SNL signifie STAT_NON_LINE ; CC signifie CREA_CHAMP ; CL signifie conditions aux limites

Etape 1 : Initialisation des contraintes
I : SNL1 avec le chargement de poids propre II : affectation par commande CC du champ

ou de pression souhaité et un matériau doté souhaité
d'un coefficient de Poisson éventuellement
fictif
SOL


Etape 2 : Récupération des réactions nodales au bord de la future galerie
CC pour extraire les contraintes issus de SNL 1 avec CL sur SNL 1 avec CL sur l'objet

SNL1
l'objet BORD en DIDI BORD en DIDI sur
SNL 2 avec CL sur l'objet BORD en DIDI
sur modèle SOL
modèle SOL_REST
SOL_REST
Récupération des réactions
Récupération des
Récupération des
réactions
réactions
BORD


Etape 3 : Déconfinement
SNL 3 avec le
SNL 2 (modèle SOL) SNL 2 (un seul matériau

chargement du
avec le chargement et modèle SOL_REST)
vecteur des SOL_REST
du vecteur des
avec le chargement du
SOL_REST
réactions nodales et
réactions nodales et vecteur des réactions
un matériau «mou»
un matériau «mou » nodales
à la place du
à la place du « vide »
« vide »
BORD
BORD

Eléments mous


Etape 4 : Pose du soutènement
SNL 4 avec 3
SNL 3 avec 3 CC pour extraire les
matériaux :
roche,
matériaux :
roche, résultats de SNL 2
béton et vide
SOL_REST
béton et vide
SNL 3 sur modèle
SOL_REST
(méthode A) pour
(méthode A) pour SOL_REST + BETON)
achever le
achever le
pour calcul intermédiare
déconfinement
déconfinement
Combianson des champs
BORD
CC
BORD
SNL 4 pour achever le
déconfinement
Eléments mous

BETON


Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 23/62


Annexe 3 Fichier de maillage réalisé avec GIBI

********************************************************
* SAUVEGARDE DANS UN FICHIER .mgib APRES COMPILATION *
********************************************************

OPTION SAUV FORM 'mail.mgib' ;

********************************************************
* OPTION DE LA MODELISATION *
********************************************************

OPTI DIME 2 ELEM QUA8;

********************************************************
* POINTS *
********************************************************

E1 = 0. 0. ;
E2 = 0.9 0. ;
E3 = 0.7 0.7 ;
E4 = 0. 0.9 ;

B1 = 1.2 0. ;
B2 = 0. 1.2 ;

S1 = 1.5 0. ;
S2 = 20. 0. ;
S3 = 20. 20. ;
S4 = 0. 20. ;
S5 = 0. 1.5 ;
S6 = (1.5 * (SIN 45)) (1.5 * (COS 45)) ;

********************************************************
* DROITES *
********************************************************

E1E2 = E1 DROI 16 E2 ;
E2E3 = E2 DROI 16 E3 ;
E3E4 = E3 DROI 16 E4 ;
E4E1 = E4 DROI 16 E1 ;

B1S1 = B1 DROI 4 S1 ;
S5B2 = S5 DROI 4 B2 ;
E2B1 = E2 DROI 4 B1 ;
B2E4 = B2 DROI 4 E4 ;

S1S2 = S1 DROI -70 S2 DINI 0.01 DFIN 0.50 ;
S2S3 = S2 DROI 16 S3 ;
S3S4 = S3 DROI 16 S4 ;
S2S3S4 = S2S3 ET S3S4 ;
S4S5 = S4 DROI -70 S5 DINI 0.50 DFIN 0.01 ;
S3S6 = S3 DROI -70 S6 DINI 0.70 DFIN 0.001 ;

********************************************************
* ARCS *
********************************************************
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 24/62


S1S5 = 32 CERC S1 E1 S5 ;
S5S6 = 16 CERC S5 E1 S6 ;
S6S1 = 16 CERC S6 E1 S1 ;
B2B1 = 32 CERC B2 E1 B1 ;

********************************************************
* DEFINITION DES GROUPES DE MAILLE *
********************************************************

BORD = S1S5 ;
MA_HAUT = S3S4 ;
BAS_BETO = B1S1 ;
LEFT_BET = S5B2 ;
NO_DROIT = S2S3 ;
NO_LEFT2 = S4S5 ;
NO_BAS2 = S1S2 ;
NO_LEFT3 = NO_LEFT2 ET LEFT_BET ;
NO_BAS3 = BAS_BETO ET NO_BAS2 ;
NO_LEFT1 = NO_LEFT3 ET B2E4 ET E4E1 ;
NO_BAS1 = E1E2 ET E2B1 ET NO_BAS3 ;


********************************************************
* SURFACES *
********************************************************

*----------------------*
* PARTIE EXCAVATION *
*----------------------*

EXCAV1 = DALL E1E2 E2E3 E3E4 E4E1 ;
TRAC EXCAV1 ;

E4E3E2 = (INVE E3E4) ET (INVE E2E3) ;
EXCAV2 = DALL E2B1 (INVE B2B1) B2E4 E4E3E2 ;
TRAC EXCAV2;

EXCAV = EXCAV1 ET EXCAV2 ;
ELIM .005 EXCAV ;
TRAC EXCAV ;

*----------------------*
* PARTIE BETON *
*----------------------*

BETON = DALL BAS_BETO BORD LEFT_BET B2B1 'PLAN';
ELIM .005 BETON ;
TRAC BETON ;

*----------------------*
* PARTIE SOL_REST *
*----------------------*

SOL1 = DALL NO_BAS2 NO_DROIT S3S6 S6S1 'PLAN';
TRAC SOL1 ;

SOL2 = DALL MA_HAUT NO_LEFT2 S5S6 (INVE S3S6)'PLAN';
TRAC SOL2 ;
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 25/62


SOL_REST = SOL1 ET SOL2 ;
ELIM .025 SOL_REST ;
TRACE SOL_REST ;

*----------------------*
* PARTIE SOL_RES0 *
*----------------------*

SOL_RES0 = BETON ET SOL_REST ;
ELIM .005 SOL_RES0 ;
TRAC SOL_RES0 ;

*----------------------*
* TOTALITE = SOL *
*----------------------*

SOL = SOL_REST ET BETON ET EXCAV ;
ELIM 0.015 SOL ;
TRACE SOL ;

********************************************************
* SAUVEGARDE DU FORMAT *
********************************************************

SAUV FORMAT SOL ;

FIN;

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 26/62


Annexe 4 Excavation sans soutènement, sur la base d'un seul
modèle (cas n°1). Fichier de commandes Code_Aster


DEBUT( );


##########################
# LECTURE MAILLAGE GIBI
##########################

PRE_GIBI();

MAIL=LIRE_MAILLAGE();


#########################################################################
# MODELISATION D'UNE EXCAVATION SANS SOUTENEMENT D'UN TUNNEL EN D.P
#########################################################################
# DEFINITION DES GROUPES DE NOEUDS POUR LESQUELS IL Y AURA
# DES DEPLACEMENTS IMPOSES
#
# NO_BAS1 : GROUPE DE NOEUDS DU BORD INFERIEUR DE TOUT LE MASSIF
# AVANT EXCAVATION.
# NO_BAS2 : GROUPE DE NOEUDS DU BORD INFERIEUR APRES EXCAVATION,
# MAIS AVANT POSE DES VOUSSOIRS.
# NO_BAS3 : GROUPE DE NOEUDS DU BORD INFERIEUR APRES EXCAVATION,
# ET POSE DES VOUSSOIRS.
#
# NO_DROIT : GROUPE DE NOEUDS DU BORD DROIT.
#
# NO_HAUT : GROUPE DE NOEUDS DU BORD SUPERIEUR.
#

# NO_LEFT1 : GROUPE DE NOEUDS DU BORD GAUCHE DE TOUT LE MASSIF
# AVANT EXCAVATION.
# NO_LEFT2 : GROUPE DE NOEUDS DU BORD GAUCHE DE TOUT LE MASSIF
# APRES EXCAVATION, MAIS AVANT POSE DES VOUSSOIRS.
# NO_LEFT2 : GROUPE DE NOEUDS DU BORD GAUCHE APRES EXCAVATION
# ET POSE DES VOUSSOIRS.
#########################################################################


#####################################
# DEFINITION DES GROUPES DE NEOUDS
#####################################
# L'OPTION "DIFFE" PERMET D'ISOLER
# DU BORD_SOL LES NEOUDS N1 ET N8359
#####################################

MAIL = DEFI_GROUP(reuse=MAIL,


MAILLAGE=MAIL,


CREA_GROUP_NO=(_F(GROUP_MA='SOL'),


_F(GROUP_MA='SOL_REST'),


_F(GROUP_MA='EXCAV'),

_F(NOM='NO_HAUT',




GROUP_MA='MA_HAUT'),

_F(GROUP_MA='NO_DROIT'),

_F(GROUP_MA='NO_LEFT1'),

_F(GROUP_MA='NO_LEFT2'),
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 27/62


_F(GROUP_MA='NO_LEFT3'),

_F(GROUP_MA='NO_BAS1'),

_F(GROUP_MA='NO_BAS2'),

_F(GROUP_MA='NO_BAS3'),

_F(GROUP_MA='LEFT_BET'),

_F(GROUP_MA='BAS_BETO'),

_F(GROUP_MA='BORD'),

_F(NOM='NOEUD1',




NOEUD='N1'),

_F(NOM='NOEUD8359',



NOEUD='N8359'),

_F(NOM='BORD_SOL',


DIFFE=('BORD','NOEUD1','NOEUD8359'),),),),

##################################################
# MODELE DU SOL AVANT L EXCAVATION POUR L'ETAPE
# D'INITIIALISATION DU CHAMP DE CONTRAINTES
##################################################
#'MA_HAUT' FIGURE DANS LE MODELE VU QU'ON
# APPLIQUE LA DESSUS UN CHARGEMENT
##################################################

MO=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA=('SOL','MA_HAUT'),
PHENOMENE='MECANIQUE',
MODELISATION='D_PLAN',),),);



##############################################
# SOL POUR INITIALISER LES CONTRAINTES
##############################################

SOL0=DEFI_MATERIAU(ELAS=_F(E=4.0E9,
NU=0.4999,
RHO=2000.0,
ALPHA=0.0,),);

########################################
# MATERIAU SOL REEL (DONNEE DU CALCUL)
########################################

SOL=DEFI_MATERIAU(ELAS=_F(E=4.0E9,
NU=0.3000,
RHO=2000.0,
ALPHA=0.0,),);

################################################
# PROPRIETES MECANIQUES DES ELEMENTS EXCAVES
################################################

VIDE=DEFI_MATERIAU(ELAS=_F(E=0.0001,
NU=0.2,
RHO=0.0,
ALPHA=0.0,),);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 28/62


#########################################
# MATERIAU CALCUL ELASTIQUE ===> CHMAT0
#########################################

CHMAT0=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL',
MATER=SOL0,),),);

#####################################################
# MATERIAU AVEC LES DONNEES DE L ETUDE ===> CHMAT1
#####################################################

CHMAT1=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL_REST',
MATER=SOL,),
_F(GROUP_MA='EXCAV',
MATER=VIDE,),
_F(GROUP_MA='BETON',
MATER=VIDE,),),);

########################################################
# LISTE DES INSTANTS DE CALCUL
########################################################
# DE 0 A 1 ==> POUR LA PHASE D'INITIALISATION
# DE 1 A 10 ==> POUR LE BLOCAGE DU BORD DE LA GALERIE
# 10 CORRESPOND A UN TEMPS DE DECONFINEMENT = 0
# 500 CORRESPOND A UN TEMPS DE DECONFINEMENT = 50 %
# 1000 CORRESPOND A UN TEMPS DE DECONFINEMENT = 100 %
########################################################

LI=DEFI_LIST_REEL(DEBUT=0,
INTERVALLE=(_F(JUSQU_A=1.0,
NOMBRE=1,),
_F(JUSQU_A=10.0,
NOMBRE=1,),
_F(JUSQU_A=500.0,
NOMBRE=1,),
_F(JUSQU_A=1000,
NOMBRE=1,),),);

######################################################
# FO FONCTION MULTIPLICATRICE POUR LE DECONFINEMENT
######################################################

F0=DEFI_FONCTION(NOM_PARA='INST',
VALE=(10.0,1.0,
500.0,0.5,
1000.0,0.0,),);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 29/62


#############################################################
# CONDITIONS AUX LIMITES EN DEPLACEMENTS :
# SYMETRIE SUR LES COTES LATERAUX => DX=0
# CONTINUITE SUR LA PARTIE INFERIEURE => DY=0
# POIDS DES TERRES SUR LA FACE SUPERIEURE => PRES
#############################################################

CH0=AFFE_CHAR_MECA(MODELE=MO,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT1',
DX=0.0,),
_F(GROUP_NO='NO_BAS1',
DY=0.0,),),
PRES_REP=_F(GROUP_MA='MA_HAUT',
PRES=5.0E6,),);


##########################
# PREMIERE STAT NON LINE #
########################################################
# 1ERE PHASE : INITIALISATION DU CHAMP DES CONTRAINTES
########################################################

RESU1=STAT_NON_LINE(MODELE=MO,
CHAM_MATER=CHMAT0,
EXCIT=(_F(CHARGE=CH0,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL',),),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=0.,
INST_FIN =1.,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=10,),
CONVERGENCE=_F(ITER_GLOB_MAXI=10,
ITER_INTE_MAXI=5,),
PARM_THETA=0.57,);


###############################
# EXTRACTION DES CONTRAINTES
###############################

RES1=CREA_CHAMP(TYPE_CHAM='ELGA_SIEF_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELGA',
INST=1,);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 30/62


##############################################################
# BLOCAGE DES NOEUDS AU BORD DE LA PARTIE EXCAVEE => DX+DY=0
##############################################################

CH210=AFFE_CHAR_MECA(MODELE=MO,
DDL_IMPO=(_F(GROUP_NO='BORD_SOL',
DX=0.0,
DY=0.0,),
_F(NOEUD=('N1'),
DX=0.0,),
_F(NOEUD=('N8359'),
DY=0.0,),),);


######################################################################
# CONDITIONS AUX LIMITES EN DEPLACEMENTS :
# SYMETRIE SUR LES COTES LATERAUX => DX=0
# CONTINUITE SUR LA PARTIE INFERIEURE => DY=0
# POIDS DES TERRES SUR LA FACE SUPERIEURE => PRES
######################################################################

CH220=AFFE_CHAR_MECA(MODELE=MO,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT2',
DX=0.0,),
_F(GROUP_NO='NO_BAS2',
DY=0.0,),
_F(GROUP_NO=('BAS_BETO'),
DY=0.0,),
_F(GROUP_NO=('LEFT_BET'),
DX=0.0,),),
PRES_REP=_F(GROUP_MA='MA_HAUT',
PRES=5.0E6,),),


##########################
# DEUXIEME STAT NON LINE #
###########################################################
# 2EME PHASE BLOCAGE DU BORD DE LA GALERIE EN DIDI
###########################################################
# RMQ : DIDI ===> DELTA U = 0
###########################################################

RESU1=STAT_NON_LINE(reuse =RESU1,
MODELE=MO,
CHAM_MATER=CHMAT0,
EXCIT=(_F(CHARGE=CH210,
TYPE_CHARGE='DIDI'),
_F(CHARGE=CH220,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL',),),
ETAT_INIT=_F(SIGM=RES1,),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=1,
INST_FIN=10,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=200,
ITER_INTE_MAXI=50,
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 31/62

ITER_INTE_PAS=-40,),
PARM_THETA=0.57,);

#####################################################################
# CONDITIONS AUX LIMITES EN DEPLACEMENTS :
# SYMETRIE SUR LES COTES LATERAUX => DX=0
# CONTINUITE SUR LA PARTIE INFERIEURE => DY=0
# POIDS DES TERRES SUR LA FACE SUPERIEURE => PRES
######################################################################

CH51=AFFE_CHAR_MECA(MODELE=MO,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT3',
DX=0.0,),
_F(NOEUD='N1',
DY=0.0,),
_F(GROUP_NO=('NO_BAS3'),
DY=0.0,),
_F(NOEUD='N8359',
DX=0.0,),),
PRES_REP=_F(GROUP_MA='MA_HAUT',
PRES=5.0E6,),);


########################
# CALCUL DES REACTIONS
########################

RESU1=CALC_NO(reuse =RESU1,
RESULTAT=RESU1,
INST=10.,
OPTION='REAC_NODA',
MODELE=MO,
CHAM_MATER=CHMAT0,
EXCIT=_F(CHARGE=CH220,),);


#######################################
# RECUPERATION DES REACTIONS NODALES
#######################################

REANODA=CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='REAC_NODA',
INST=10.,);


########################################################################
# CONSTITUTION D UN VECTEUR CHARGEMENT OBTENU CONSTITUE DES REACTIONS
########################################################################

CH3=AFFE_CHAR_MECA(MODELE=MO,
VECT_ASSE=REANODA,);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 32/62


###########################
# TROISIEME STAT NON LINE #
##########################################################
# 3EME PHASE : REINJECTION DE LA REACTION
##########################################################

RESU1=STAT_NON_LINE(reuse =RESU1,
MODELE=MO,
CHAM_MATER=CHMAT1,
EXCIT=(_F(CHARGE=CH3,
FONC_MULT=F0,),
_F(CHARGE=CH51,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL_REST',),
_F(RELATION='ELAS',
GROUP_MA='EXCAV',),
_F(RELATION='ELAS',
GROUP_MA=('BETON'),),),
ETAT_INIT=_F(EVOL_NOLI=RESU1,),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=10,
INST_FIN=1000,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=500,
ITER_INTE_MAXI=100,
ITER_INTE_PAS=-10,),
PARM_THETA=0.57,);

#############################
# CALCULS ET POST TRAITEMENT
#############################

RESU1=CALC_ELEM(reuse =RESU1,
MODELE=MO,
GROUP_MA='SOL_REST',
CHAM_MATER=CHMAT1,
OPTION=('SIEF_ELNO_ELGA',),
RESULTAT=RESU1,);

RESU1=CALC_NO(reuse=RESU1,
CHAM_MATER=CHMAT1,
OPTION=('SIEF_NOEU_ELGA',),
RESULTAT=RESU1)

##############################################
# IMPRESSION DES RESULTATS EN FORMAT CASTEM
# POUR VISUALISATION DES ISOVALEURS
##############################################

IMPR_RESU(MODELE=MO,
RESU=_F(FORMAT='CASTEM',
MAILLAGE=MAIL,
RESULTAT=RESU1,),);


Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 33/62


#############################
# TABLES DE POST-TRAITEMENT
#############################

#-------------------------------------------------
# DEPLACEMENTS NOEUD N1 FONCTION DU DECONFINEMENT
#-------------------------------------------------

DEP_1=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_FONC_DECONF_N1',
NOEUD='N1',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
TOUT_ORDRE='OUI',
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_1,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','DX','DY',),);

#------------------------------------------------
# CONTRAINTES NOEUD N1 FONCTION DU DECONFINEMENT
#------------------------------------------------

SIG_1=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_FONC_DECONF_N1',
NOEUD='N1',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
TOUT_ORDRE='OUI',
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_1,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','SIXX','SIYY'),);

#----------------------------------------------------
# DEPLACEMENTS NOEUD N8359 FONCTION DU DECONFINEMENT
#----------------------------------------------------

DEP_8359=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_FONC_DECONF_N8359',
NOEUD='N8359',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
TOUT_ORDRE='OUI',
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_8359,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','DX','DY',),);

#---------------------------------------------------
# CONTRAINTES NOEUD N8359 FONCTION DU DECONFINEMENT
#---------------------------------------------------
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 34/62


SIG_8359=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_FONC_DECONF_N8359',
NOEUD='N8359',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
TOUT_ORDRE='OUI',
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_8359,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','SIXX','SIYY'),);

#-------------------------------------------
# DEPLACEMENTS NO_LEFT2 ===> 50 %
#-------------------------------------------

DEP_L50=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_LEFT2_50%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(500),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_L50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

#-------------------------------------------
# DEPLACEMENTS NO_LEFT2 ===> 100 %
#-------------------------------------------

DEP_L100=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_LEFT2_100%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(1000),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_L100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);


#--------------------------------------------
# CONTRAINTES NO_LEFT2 INSTANT ===> 50 %
#--------------------------------------------

SIG_L50=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_LEFT2_50%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(500),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_L50,
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 35/62

FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#--------------------------------------------
# CONTRAINTES NO_LEFT2 INSTANT ===> 100 %
#--------------------------------------------

SIG_L100=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_LEFT2_100%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(1000),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_L100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#-------------------------------------------
# DEPLACEMENTS NO_BAS2 ===> 50 %
#-------------------------------------------

DEP_B50=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_NO_BAS2_50%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(500),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_B50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

#-------------------------------------------
# DEPLACEMENTS NO_BAS2 ===> 100 %
#-------------------------------------------

DEP_B100=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_NO_BAS2_100%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(1000),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_B100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 36/62


#--------------------------------------------
# CONTRAINTES NO_BAS2 INSTANT ===> 50 %
#--------------------------------------------

SIG_B50=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_BAS2_50%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(500),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_B50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#--------------------------------------------
# CONTRAINTES NO_BAS2 INSTANT ===> 100 %
#--------------------------------------------

SIG_B100=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_BAS2_100%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(1000),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_B100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);
FIN();
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 37/62


Annexe 5 Excavation avec soutènement, méthode A (cas n°2).
Fichier de commandes Code_Aster


DEBUT( );


##########################
# LECTURE MAILLAGE GIBI
##########################

PRE_GIBI();

MAIL=LIRE_MAILLAGE();


#########################################################################
# MODELISATION D'UNE EXCAVATION AVEC SOUTENEMENT D'UN TUNNEL EN D.P
#########################################################################
# DEFINITION DES GROUPES DE NOEUDS POUR LESQUELS IL Y AURA
# DES DEPLACEMENTS IMPOSES
#
# NO_BAS1 : GROUPE DE NOEUDS DU BORD INFERIEUR DE TOUT LE MASSIF
# AVANT EXCAVATION.
# NO_BAS2 : GROUPE DE NOEUDS DU BORD INFERIEUR APRES EXCAVATION,
# MAIS AVANT POSE DES VOUSSOIRS.
# NO_BAS3 : GROUPE DE NOEUDS DU BORD INFERIEUR APRES EXCAVATION,
# ET POSE DES VOUSSOIRS.
#
# NO_DROIT : GROUPE DE NOEUDS DU BORD DROIT.
#
# NO_HAUT : GROUPE DE NOEUDS DU BORD SUPERIEUR.
#
# NO_LEFT1 : GROUPE DE NOEUDS DU BORD GAUCHE DE TOUT LE MASSIF
# AVANT EXCAVATION.
# NO_LEFT2 : GROUPE DE NOEUDS DU BORD GAUCHE DE TOUT LE MASSIF
# APRES EXCAVATION, MAIS AVANT POSE DES VOUSSOIRS.
# NO_LEFT2 : GROUPE DE NOEUDS DU BORD GAUCHE APRES EXCAVATION
# ET POSE DES VOUSSOIRS.
#########################################################################

#####################################
# DEFINITION DES GROUPES DE NEOUDS
#####################################
# L'OPTION "DIFFE" PERMET D'ISOLER
# DU BORD_SOL LES NEOUDS N1 ET N8359
#####################################

MAIL = DEFI_GROUP(reuse=MAIL,
MAILLAGE=MAIL,
CREA_GROUP_NO=(_F(GROUP_MA='SOL'),
_F(GROUP_MA='SOL_REST'),
_F(GROUP_MA='EXCAV'),
_F(NOM='NO_HAUT',
GROUP_MA='MA_HAUT'),
_F(GROUP_MA='NO_DROIT'),
_F(GROUP_MA='NO_LEFT1'),
_F(GROUP_MA='NO_LEFT2'),
_F(GROUP_MA='NO_LEFT3'),
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 38/62

_F(GROUP_MA='NO_BAS1'),
_F(GROUP_MA='NO_BAS2'),
_F(GROUP_MA='NO_BAS3'),
_F(GROUP_MA='LEFT_BET'),
_F(GROUP_MA='BAS_BETO'),
_F(GROUP_MA='BORD'),
_F(NOM='NOEUD1',
NOEUD='N1'),
_F(NOM='NOEUD8359',
NOEUD='N8359'),
_F(NOM='BORD_SOL',
DIFFE=('BORD','NOEUD1','NOEUD8359'),),),),

##################################################
# MODELE DU SOL AVANT L EXCAVATION POUR L'ETAPE
# D'INITIIALISATION DU CHAMP DE CONTRAINTES
##################################################
#'MA_HAUT' FIGURE DANS LE MODELE VU QU'ON
# APPLIQUE LA DESSUS UN CHARGEMENT
##################################################

MO=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA=('SOL','MA_HAUT'),
PHENOMENE='MECANIQUE',
MODELISATION='D_PLAN',),),);

##################################
# SOL POUR CALCUL REACTIONS NODALES
##################################

SOL=DEFI_MATERIAU(ELAS=_F(E=4.0E9,
NU=0.4999,
RHO=2000.0,
ALPHA=0.0,),);

#########################################
# MATERIAU SOL REEL
#########################################

SOL2=DEFI_MATERIAU(ELAS=_F(E=4.0E9,
NU=0.3000,
RHO=2000.0,
ALPHA=0.0,),);

###################################################
# VIDE
###################################################

VIDE=DEFI_MATERIAU(ELAS=_F(E=0.001,
NU=0.2,
RHO=0.0,
ALPHA=0.0,),);

####################################################
# PROPRIETES MECANIQUES DES VOUSSOIRS EN BETON
####################################################

BETON1=DEFI_MATERIAU(ELAS=_F(E=2.E10,
NU=0.2,
RHO=0.0,
ALPHA=0.0,),);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 39/62


#########################################
# MATERIAU CALCUL ELASTIQUE ===> CHMAT0
#########################################

CHMAT0=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL',
MATER=SOL,),),);



#####################################################
# MATERIAU AVEC LES DONNEES DE L ETUDE ===> CHMAT2
#####################################################

CHMAT2=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL_REST',
MATER=SOL2,),
_F(GROUP_MA='EXCAV',
MATER=VIDE,),
_F(GROUP_MA='BETON',
MATER=VIDE,),),);

#####################################################
# MATERIAU AVEC LES DONNEES DE L ETUDE ===> CHMAT3
#####################################################

CHMAT3=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL_REST',
MATER=SOL2,),
_F(GROUP_MA='EXCAV',
MATER=VIDE,),
_F(GROUP_MA='BETON',
MATER=BETON1,),),);

########################################################
# LISTE DES INSTANTS DE CALCUL
########################################################
# DE 0 A 1 ==> POUR LA PHASE D'INITIALISATION
# DE 1 A 10 ==> POUR LE BLOCAGE DU BORD DE LA GALERIE
# 10 CORRESPOND A UN TEMPS DE DECONFINEMENT = 0
# 500 CORRESPOND A UN TEMPS DE DECONFINEMENT = 50 %
# 1000 CORRESPOND A UN TEMPS DE DECONFINEMENT = 100 %
########################################################

LI=DEFI_LIST_REEL(DEBUT=0,
INTERVALLE=(_F(JUSQU_A=1.0,
NOMBRE=1,),
_F(JUSQU_A=10.0,
NOMBRE=1,),
_F(JUSQU_A=500,
NOMBRE=1,),
_F(JUSQU_A=1000,
NOMBRE=1,),),);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 40/62


######################################################
# FO FONCTION MULTIPLICATRICE POUR LE DECONFINEMENT
######################################################

F0=DEFI_FONCTION(NOM_PARA='INST',
VALE=(10.0,1.0,
500.0,0.5,
1000.0,0.0,),);


#################################################
# CREA CHAMP PAR l'OPERATEUR 'AFFE'
#################################################
# 1ERE PHASE : INITIALISATION DES CONTRAINTES
#################################################

RES1=CREA_CHAMP(TYPE_CHAM='ELNO_SIEF_R',
OPERATION='AFFE',
MODELE=MO,
AFFE=(_F(TOUT='OUI',
NOM_CMP=('SIXX','SIYY','SIZZ','SIXY'),
VALE=(5.0E6,5.0E6,0.,0.),),),),


##############################################################
# BLOCAGE DES NOEUDS AU BORD DE LA PARTIE EXCAVEE => DX+DY=0
##############################################################

CH210=AFFE_CHAR_MECA(MODELE=MO,
DDL_IMPO=(_F(GROUP_NO='BORD_SOL',
DX=0.0,
DY=0.0,),
_F(NOEUD=('N1'),
DX=0.0,),
_F(NOEUD=('N8359'),
DY=0.0,),),);


######################################################################
# CONDITIONS AUX LIMITES EN DEPLACEMENTS :
# SYMETRIE SUR LES COTES LATERAUX => DX=0
# CONTINUITE SUR LA PARTIE INFERIEURE => DY=0
# POIDS DES TERRES SUR LA FACE SUPERIEURE => PRES
######################################################################

CH220=AFFE_CHAR_MECA(MODELE=MO,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT2',
DX=0.0,),
_F(GROUP_NO='NO_BAS2',
DY=0.0,),
_F(GROUP_NO=('BAS_BETO'),
DY=0.0,),
_F(GROUP_NO=('LEFT_BET'),
DX=0.0,),),
PRES_REP=_F(GROUP_MA='MA_HAUT',
PRES=5.0E6,),),
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 41/62



##########################
# PREMIERE STAT NON LINE #
###########################################################
# 2EME PHASE BLOCAGE DU BORD DE LA GALERIE EN DIDI
###########################################################
# RMQ : DIDI ===> DELTA U = 0
###########################################################

RESU1=STAT_NON_LINE(MODELE=MO,
CHAM_MATER=CHMAT0,
EXCIT=(_F(CHARGE=CH210,
TYPE_CHARGE='DIDI'),
_F(CHARGE=CH220,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL',),),
ETAT_INIT=_F(SIGM=RES1,),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=1,
INST_FIN=10,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=200,
ITER_INTE_MAXI=50,
ITER_INTE_PAS=-40,),
PARM_THETA=0.57,);


################################
# EXTRACTION DES CONTRAINTES
################################

RES2=CREA_CHAMP(TYPE_CHAM='ELGA_SIEF_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELGA',
INST=10,);


#####################################################################
# CONDITIONS AUX LIMITES EN DEPLACEMENTS :
# SYMETRIE SUR LES COTES LATERAUX => DX=0
# CONTINUITE SUR LA PARTIE INFERIEURE => DY=0
# POIDS DES TERRES SUR LA FACE SUPERIEURE => PRES
######################################################################

CH51=AFFE_CHAR_MECA(MODELE=MO,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT3',
DX=0.0,),
_F(NOEUD='N1',
DY=0.0,),
_F(GROUP_NO=('NO_BAS3'),
DY=0.0,),
_F(NOEUD='N8359',
DX=0.0,),),
PRES_REP=_F(GROUP_MA='MA_HAUT',
PRES=5.0E6,),);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 42/62


########################
# CALCUL DES REACTIONS
########################

RESU1=CALC_NO(reuse =RESU1,
RESULTAT=RESU1,
INST=10.,
OPTION='REAC_NODA',
MODELE=MO,
CHAM_MATER=CHMAT0,
EXCIT=_F(CHARGE=CH220,),);

#######################################
# RECUPERATION DES REACTIONS NODALES
#######################################

REANODA=CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='REAC_NODA',
INST=10.,);

########################################################################
# CONSTITUTION D UN VECTEUR CHARGEMENT OBTENU CONSTITUE DES REACTIONS
########################################################################

CH3=AFFE_CHAR_MECA(MODELE=MO,
VECT_ASSE=REANODA,);

###########################
# DEUXIEME STAT NON LINE #
##########################################################
# 3EME PHASE : REINJECTION DE LA REACTION
##########################################################

RESU1=STAT_NON_LINE(reuse =RESU1,
MODELE=MO,
CHAM_MATER=CHMAT2,
EXCIT=(_F(CHARGE=CH3,
FONC_MULT=F0,),
_F(CHARGE=CH51,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL_REST',),
_F(RELATION='ELAS',
GROUP_MA='EXCAV',),
_F(RELATION='ELAS',
GROUP_MA=('BETON'),),),
ETAT_INIT=_F(EVOL_NOLI=RESU1,),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=10,
INST_FIN=500,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=500,
ITER_INTE_MAXI=100,
ITER_INTE_PAS=-10,),
PARM_THETA=0.57,);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 43/62


############################
# TROISIEME STAT NON LINE #
##########################################################
# 4EME PHASE : ACTIVATION DU BETON
##########################################################
RESU1=STAT_NON_LINE(reuse =RESU1,
MODELE=MO,
CHAM_MATER=CHMAT3,
EXCIT=(_F(CHARGE=CH3,
FONC_MULT=F0,),
_F(CHARGE=CH51,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL_REST',),
_F(RELATION='ELAS',
GROUP_MA='EXCAV',),
_F(RELATION='ELAS',
GROUP_MA=('BETON'),),),
ETAT_INIT=_F(EVOL_NOLI=RESU1,),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=500,
INST_FIN=1000,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=500,
ITER_INTE_MAXI=100,
ITER_INTE_PAS=-10,),
PARM_THETA=0.57,);

#############################
# CALCULS ET POST TRAITEMENT
#############################

RESU1=CALC_ELEM(reuse = RESU1,
MODELE=MO,
CHAM_MATER=CHMAT2,
GROUP_MA='SOL_REST',
OPTION=('SIEF_ELNO_ELGA',),
RESULTAT=RESU1,);

RESU1=CALC_NO(reuse = RESU1,
CHAM_MATER=CHMAT2,
OPTION=('SIEF_NOEU_ELGA',),
RESULTAT=RESU1)

##############################################
# IMPRESSION DES RESULTATS EN FORMAT CASTEM
# POUR VISUALISATION DES ISOVALEURS
##############################################

#IMPR_RESU(MODELE=MO2,
# RESU=_F(FORMAT='CASTEM',
# MAILLAGE=MAIL,
# RESULTAT=RESU1,),);


Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 44/62


#############################
# TABLES DE POST-TRAITEMENT
#############################

#-------------------------------------------------
# DEPLACEMENTS NOEUD N1 FONCTION DU DECONFINEMENT
#-------------------------------------------------

DEP_1=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_FONC_DECONF_N1',
NOEUD='N1',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
TOUT_ORDRE='OUI',
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_1,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','DX','DY',),);

#------------------------------------------------
# CONTRAINTES NOEUD N1 FONCTION DU DECONFINEMENT
#------------------------------------------------

SIG_1=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_FONC_DECONF_N1',
NOEUD='N1',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
TOUT_ORDRE='OUI',
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_1,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','SIXX','SIYY'),);

#----------------------------------------------------
# DEPLACEMENTS NOEUD N8359 FONCTION DU DECONFINEMENT
#----------------------------------------------------

DEP_8359=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_FONC_DECONF_N8359',
NOEUD='N8359',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
TOUT_ORDRE='OUI',
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_8359,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','DX','DY',),);

#---------------------------------------------------
# CONTRAINTES NOEUD N8359 FONCTION DU DECONFINEMENT
#---------------------------------------------------
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 45/62


SIG_8359=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_FONC_DECONF_N8359',
NOEUD='N8359',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
TOUT_ORDRE='OUI',
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_8359,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','SIXX','SIYY'),);

#-------------------------------------------
# DEPLACEMENTS NO_LEFT2 ===> 50 %
#-------------------------------------------

DEP_L50=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_LEFT2_50%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(500),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_L50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

#-------------------------------------------
# DEPLACEMENTS NO_LEFT2 ===> 100 %
#-------------------------------------------

DEP_L100=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_LEFT2_100%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(1000),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_L100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);


#--------------------------------------------
# CONTRAINTES NO_LEFT2 INSTANT ===> 50 %
#--------------------------------------------

SIG_L50=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_LEFT2_50%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(500),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 46/62


IMPR_TABLE(TABLE=SIG_L50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#--------------------------------------------
# CONTRAINTES NO_LEFT2 INSTANT ===> 100 %
#--------------------------------------------

SIG_L100=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_LEFT2_100%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(1000),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_L100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#-------------------------------------------
# DEPLACEMENTS NO_BAS2 ===> 50 %
#-------------------------------------------

DEP_B50=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_NO_BAS2_50%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(500),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_B50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

#-------------------------------------------
# DEPLACEMENTS NO_BAS2 ===> 100 %
#-------------------------------------------

DEP_B100=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_NO_BAS2_100%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(1000),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_B100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 47/62


#--------------------------------------------
# CONTRAINTES NO_BAS2 INSTANT ===> 50 %
#--------------------------------------------

SIG_B50=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_BAS2_50%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(500),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_B50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#--------------------------------------------
# CONTRAINTES NO_BAS2 INSTANT ===> 100 %
#--------------------------------------------

SIG_B100=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_BAS2_100%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(1000),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_B100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);


FIN();
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 48/62


Annexe 6 Excavation avec soutènement, méthode B (cas n°3).
Fichier de commandes Code_Aster


DEBUT( );


##########################
# LECTURE MAILLAGE GIBI
##########################

PRE_GIBI();

MAIL=LIRE_MAILLAGE();

#########################################################################
# MODELISATION D'UNE EXCAVATION AVEC SOUTENEMENT D'UN TUNNEL EN D.P
#########################################################################
# DEFINITION DES GROUPES DE NOEUDS POUR LESQUELS IL Y AURA
# DES DEPLACEMENTS IMPOSES
#
# NO_BAS1 : GROUPE DE NOEUDS DU BORD INFERIEUR DE TOUT LE MASSIF
# AVANT EXCAVATION.
# NO_BAS2 : GROUPE DE NOEUDS DU BORD INFERIEUR APRES EXCAVATION,
# MAIS AVANT POSE DES VOUSSOIRS.
# NO_BAS3 : GROUPE DE NOEUDS DU BORD INFERIEUR APRES EXCAVATION,
# ET POSE DES VOUSSOIRS.
#
# NO_DROIT : GROUPE DE NOEUDS DU BORD DROIT.
#
# NO_HAUT : GROUPE DE NOEUDS DU BORD SUPERIEUR.
#
# NO_LEFT1 : GROUPE DE NOEUDS DU BORD GAUCHE DE TOUT LE MASSIF
# AVANT EXCAVATION.
# NO_LEFT2 : GROUPE DE NOEUDS DU BORD GAUCHE DE TOUT LE MASSIF
# APRES EXCAVATION, MAIS AVANT POSE DES VOUSSOIRS.
# NO_LEFT2 : GROUPE DE NOEUDS DU BORD GAUCHE APRES EXCAVATION
# ET POSE DES VOUSSOIRS.
#########################################################################


#####################################
# DEFINITION DES GROUPES DE NEOUDS
#####################################
# L'OPTION "DIFFE" PERMET D'ISOLER
# DU BORD_SOL LES NEOUDS N1 ET N8359
#####################################

MAIL = DEFI_GROUP(reuse=MAIL,
MAILLAGE=MAIL,
CREA_GROUP_NO=(_F(GROUP_MA='SOL'),
_F(GROUP_MA='SOL_REST'),
_F(NOM='NO_HAUT',
GROUP_MA='MA_HAUT'),
_F(GROUP_MA='NO_DROIT'),
_F(GROUP_MA='NO_LEFT1'),
_F(GROUP_MA='NO_LEFT2'),
_F(GROUP_MA='NO_LEFT3'),
_F(GROUP_MA='NO_BAS1'),
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 49/62

_F(GROUP_MA='NO_BAS2'),
_F(GROUP_MA='NO_BAS3'),
_F(GROUP_MA='LEFT_BET'),
_F(GROUP_MA='BAS_BETO'),
_F(GROUP_MA='BORD'),
_F(NOM='NOEUD1',
NOEUD='N1'),
_F(NOM='NOEUD8359',
NOEUD='N8359'),
_F(NOM='BORD_SOL',
DIFFE=('BORD','NOEUD1','NOEUD8359'),),),),


##################################################
# MODELE DU SOL AVANT L EXCAVATION POUR L'ETAPE
# D'INITIIALISATION DU CHAMP DE CONTRAINTES
##################################################
#'MA_HAUT' FIGURE DANS LE MODELE VU QU'ON
# APPLIQUE LA DESSUS UN CHARGEMENT
##################################################

MO=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA=('SOL','MA_HAUT'),
PHENOMENE='MECANIQUE',
MODELISATION='D_PLAN',),),);

####################
# MODELE SOL_REST
####################

MO1=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA=('SOL_REST','MA_HAUT'),
PHENOMENE='MECANIQUE',
MODELISATION='D_PLAN',),),);

#####################################
# MODELE DU SOL APRES L EXCAVATION
#####################################

MO2=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='BETON',
PHENOMENE='MECANIQUE',
MODELISATION='D_PLAN',),
_F(GROUP_MA=('SOL_REST','MA_HAUT'),
PHENOMENE='MECANIQUE',
MODELISATION='D_PLAN',),),);


################################################
# SOL POUR CALCUL REACTIONS NODALES
################################################

SOL0=DEFI_MATERIAU(ELAS=_F(E=4.0E9,
NU=0.4999,
RHO=2000.0,
ALPHA=0.0,),);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 50/62


#########################################
# MATERIAU SOL REEL (DONNEES DU CALCUL)
#########################################

SOL=DEFI_MATERIAU(ELAS=_F(E=4.0E9,
NU=0.30,
RHO=2000.0,
ALPHA=0.0,),);



####################################################
# PROPRIETES MECANIQUES DES VOUSSOIRS EN BETON
####################################################

BETON=DEFI_MATERIAU(ELAS=_F(E=2.E10,
NU=0.2,
RHO=0.0,
ALPHA=0.0,),);


#########################################
# MATERIAU CALCUL ELASTIQUE ===> CHMAT0
#########################################

CHMAT0=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL',
MATER=SOL0,),),);


##############################################
# MATERIAU PHASE D INITIALISATION ===> CHMAT
##############################################

CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL_REST',
MATER=SOL0,),),);

############################################################
# MATERIAU PHASE DE REINJECTION DE LA REACTION ===> CHMAT1
############################################################

CHMAT1=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL_REST',
MATER=SOL,),),);

#####################################################
# MATERIAU AVEC LES DONNEES DE L ETUDE ===> CHMAT2
#####################################################

CHMAT2=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='SOL_REST',
MATER=SOL,),
_F(GROUP_MA='BETON',
MATER=BETON,),),);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 51/62


########################################################
# LISTE DES INSTANTS DE CALCUL
########################################################
# DE 0 A 1 ==> POUR LA PHASE D'INITIALISATION
# DE 1 A 10 ==> POUR LE BLOCAGE DU BORD DE LA GALERIE
# 10 CORRESPOND A UN TEMPS DE DECONFINEMENT = 0
# 500 CORRESPOND A UN TEMPS DE DECONFINEMENT = 50 %
# 1000 CORRESPOND A UN TEMPS DE DECONFINEMENT = 100 %
########################################################

LI=DEFI_LIST_REEL(DEBUT=0,
INTERVALLE=(_F(JUSQU_A=1.0,
NOMBRE=1,),
_F(JUSQU_A=10.0,
NOMBRE=1,),
_F(JUSQU_A=500.0,
NOMBRE=1,),
_F(JUSQU_A=1000,
NOMBRE=1,),),);


##################################################
# LISTE DES INSTANTS POUR LE CALCUL BIDON POUR
# INITIALISER LES CHAMPS A 0 DALS LE BETON
##################################################

LI1=DEFI_LIST_REEL(DEBUT=0,
INTERVALLE=(_F(JUSQU_A=1.E6,
NOMBRE=1),),);


######################################################
# FO FONCTION MULTIPLICATRICE POUR LE DECONFINEMENT
######################################################

F0=DEFI_FONCTION(NOM_PARA='INST',
VALE=(10.0,1.0,
500.0,0.5,
1000.0,0.0,),);


#################################################
# CREA CHAMP PAR l'OPERATEUR 'AFFE'
#################################################
# 1ERE PHASE : INITIALISATION DES CONTRAINTES
#################################################

RES1=CREA_CHAMP(TYPE_CHAM='ELNO_SIEF_R',
OPERATION='AFFE',
MODELE=MO,
AFFE=(_F(TOUT='OUI',
NOM_CMP=('SIXX','SIYY','SIZZ','SIXY'),
VALE=(5.0E6,5.0E6,0.,0.),),),),

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 52/62


##############################################################
# BLOCAGE DES NOEUDS AU BORD DE LA PARTIE EXCAVEE => DX+DY=0
# RQ = ON TRAVAILLE MAINTENANT AVEC LE MODELE MO2
##############################################################

CH210=AFFE_CHAR_MECA(MODELE=MO1,
DDL_IMPO=(_F(GROUP_NO='BORD_SOL',
DX=0.0,
DY=0.0,),
_F(NOEUD=('N1'),
DX=0.0,),
_F(NOEUD=('N8359'),
DY=0.0,),),);

######################################################################
# CONDITIONS AUX LIMITES EN DEPLACEMENTS =
# SYMETRIE SUR LES COTES LATERAUX => DX=0
# CONTINUITE SUR LA PARTIE INFERIEURE => DY=0
# POIDS DES TERRES SUR LA FACE SUPERIEURE => PRES
######################################################################

CH220=AFFE_CHAR_MECA(MODELE=MO1,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT2',
DX=0.0,),
_F(GROUP_NO='NO_BAS2',
DY=0.0,),),
PRES_REP=_F(GROUP_MA='MA_HAUT',
PRES=5.0E6,),),


##########################
# PREMIERE STAT NON LINE #
###########################################################
# 2EME PHASE BLOCAGE DU BORD DE LA GALERIE EN DIDI
###########################################################

RESU1=STAT_NON_LINE(MODELE=MO1,
CHAM_MATER=CHMAT,
EXCIT=(_F(CHARGE=CH210,
TYPE_CHARGE='DIDI'),
_F(CHARGE=CH220,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL_REST',),),
ETAT_INIT=_F(SIGM=RES1,),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=1,
INST_FIN=10,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=200,
ITER_INTE_MAXI=50,
ITER_INTE_PAS=-40,),
PARM_THETA=0.57,);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 53/62


################################
# EXTRACTION DES CONTRAINTES
################################

RES2=CREA_CHAMP(TYPE_CHAM='ELGA_SIEF_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELGA',
INST=10,);


#####################################################################
# CONDITIONS AUX LIMITES EN DEPLACEMENTS =
# SYMETRIE SUR LES COTES LATERAUX => DX=0
# CONTINUITE SUR LA PARTIE INFERIEURE => DY=0
#####################################################################

CH51=AFFE_CHAR_MECA(MODELE=MO2,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT3',
DX=0.0,),
_F(NOEUD='N1',
DY=0.0,),
_F(GROUP_NO=('NO_BAS3'),
DY=0.0,),
_F(NOEUD='N8359',
DX=0.0,),),
PRES_REP=_F(GROUP_MA='MA_HAUT',
PRES=5.0E6,),);


########################
# CALCUL DES REACTIONS
########################

RESU1=CALC_NO(reuse =RESU1,
RESULTAT=RESU1,
INST=10.,
OPTION='REAC_NODA',
MODELE=MO1,
CHAM_MATER=CHMAT,
EXCIT=_F(CHARGE=CH220,),);


#######################################
# RECUPERATION DES REACTIONS NODALES
#######################################

REANODA=CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='REAC_NODA',
INST=10.,);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 54/62


########################################################################
# CONSTITUTION D UN VECTEUR CHARGEMENT OBTENU CONSTITUE DES REACTIONS
########################################################################

CH3=AFFE_CHAR_MECA(MODELE=MO1,
VECT_ASSE=REANODA,);


###########################
# DEUXIEME STAT NON LINE #
##########################################################
# 3EME PHASE = REINJECTION DE LA REACTION
##########################################################

RESU1=STAT_NON_LINE(reuse=RESU1,
MODELE=MO1,
CHAM_MATER=CHMAT1,
EXCIT=(_F(CHARGE=CH3,
FONC_MULT=F0,),
_F(CHARGE=CH220,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL_REST',),),
ETAT_INIT=_F(EVOL_NOLI=RESU1),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=10,
INST_FIN=500,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=500,
ITER_INTE_MAXI=100,
ITER_INTE_PAS=-10,),
PARM_THETA=0.57,);


###########################################################
# EXTRACTION DES CHAMPS = DEPLACEMENTS, CONTRAINTES,
# ET VARIABLES INTERNES OBTENUS LORS DU CALCUL PRECEDENT
###########################################################

SIG1=CREA_CHAMP(TYPE_CHAM='ELGA_SIEF_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELGA',
INST=500.,);

DEP1=CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='EXTR',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=500.,);


CHBID = AFFE_CHAR_MECA(MODELE=MO1,
DDL_IMPO=(_F(GROUP_NO='NO_DROIT',
DX=0.0,),
_F(GROUP_NO='NO_LEFT2',
DX=0.0,),
_F(GROUP_NO='NO_BAS2',
DY=0.0,),),),
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 55/62


###########################
# TROISIEME STAT NON LINE #
#########################################################
# 4EME PHASE = CHARGEMENT BIDON,
# CALCUL BIDON POUR PERMETTRE :
# - UNE INITIALISATION DES CHAMPS A 0 DANS LE BETON,
# - ET ENSUITE L'ASSEMBLAGE
#########################################################

BIDON=STAT_NON_LINE(MODELE=MO2,
CHAM_MATER=CHMAT2,
EXCIT=_F(CHARGE=CHBID),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL_REST',),
_F(RELATION='ELAS',
GROUP_MA='BETON',),),
INCREMENT=_F(LIST_INST=LI1,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_MAXI=1.,
ITER_GLOB_MAXI=1,
ITER_INTE_MAXI=10,
ITER_INTE_PAS=-10,),
PARM_THETA=0.57,);


DEP2=CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='EXTR',
RESULTAT=BIDON,
NOM_CHAM='DEPL',
INST=1.E6);


SIG2=CREA_CHAMP(TYPE_CHAM='ELGA_SIEF_R',
OPERATION='EXTR',
RESULTAT=BIDON,
NOM_CHAM='SIEF_ELGA',
INST=1.E6);

###################
# ASSEMBLAGE
###################

DEPINI=CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
OPERATION='ASSE',
MAILLAGE=MAIL,
ASSE=(_F( CHAM_GD=DEP2,
GROUP_MA='BETON',
CUMUL='OUI',
COEF_R=0.),
_F(CHAM_GD=DEP1,
GROUP_MA='SOL_REST',
CUMUL='OUI',),),);

SIGINI=CREA_CHAMP(TYPE_CHAM='ELGA_SIEF_R',
OPERATION='ASSE',
MODELE=MO2,
ASSE=(_F(CHAM_GD=SIG2,
GROUP_MA='BETON',
CUMUL='OUI',
COEF_R=0.),
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 56/62

_F(CHAM_GD=SIG1,
GROUP_MA='SOL_REST',
CUMUL='OUI',),),);

# VARINI=CREA_CHAMP(TYPE_CHAM='ELGA_VARI_R',
# OPERATION='ASSE',
# MAILLAGE=MAIL,
# MODELE=MO2,
# ASSE=_F(CHAM_GD=VAR2,
# GROUP_MA=(BETON,),
# CUMUL='OUI',
# COEF_R=0.),
# ASSE=_F(CHAM_GD=VAR1,
# GROUP_MA=(SOL_REST,),
# CUMUL='OUI',),);

############################
# QUATRIEMME STAT NON LINE #
#######################################################
# 5EME PHASE = PRESENCE DES VOUSSOIRS,
# DECONFINEMENT DE 50 A 100 %
#######################################################

RESU1=STAT_NON_LINE(reuse=RESU1,
MODELE=MO2,
CHAM_MATER=CHMAT2,
EXCIT=(_F(CHARGE=CH3,
FONC_MULT=F0,),
_F(CHARGE=CH51,),),
COMP_INCR=(_F(RELATION='ELAS',
GROUP_MA='SOL_REST',),
_F(RELATION='ELAS',
GROUP_MA='BETON',),),
ETAT_INIT=_F(DEPL=DEPINI,
SIGM=SIGINI,),
INCREMENT=_F(LIST_INST=LI,
INST_INIT=500.,
INST_FIN=1000.,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=5.E-6,
ITER_GLOB_MAXI=500,
ITER_INTE_MAXI=100,
ITER_INTE_PAS=-10,),
PARM_THETA=0.57,);

#############################
# CALCULS ET POST TRAITEMENT
#############################

RESU1=CALC_ELEM(reuse=RESU1,
MODELE=MO2,
CHAM_MATER=CHMAT2,
OPTION=('SIEF_ELNO_ELGA',),
RESULTAT=RESU1,);

RESU1=CALC_NO(reuse=RESU1,
CHAM_MATER=CHMAT2,
OPTION=('SIEF_NOEU_ELGA',),
RESULTAT=RESU1);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 57/62


##############################################
# IMPRESSION DES RESULTATS EN FORMAT CASTEM
# POUR VISUALISATION DES ISOVALEURS
##############################################
# IMPR_RESU(MODELE=MO2,
# RESU=_F(FORMAT='CASTEM',
# MAILLAGE=MAIL,
# RESULTAT=RESU1,),);

#############################
# TABLES DE POST-TRAITEMENT
#############################

#-------------------------------------------------
# DEPLACEMENTS NOEUD N1 FONCTION DU DECONFINEMENT
#-------------------------------------------------

DEP_1=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_FONC_DECONF_N1',
NOEUD='N1',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
TOUT_ORDRE='OUI',
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_1,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','DX','DY',),);

#------------------------------------------------
# CONTRAINTES NOEUD N1 FONCTION DU DECONFINEMENT
#------------------------------------------------

SIG_1=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_FONC_DECONF_N1',
NOEUD='N1',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
TOUT_ORDRE='OUI',
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_1,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','SIXX','SIYY'),);

#----------------------------------------------------
# DEPLACEMENTS NOEUD N8359 FONCTION DU DECONFINEMENT
#----------------------------------------------------

DEP_8359=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_FONC_DECONF_N8359',
NOEUD='N8359',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
TOUT_ORDRE='OUI',
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 58/62


IMPR_TABLE(TABLE=DEP_8359,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','DX','DY',),);

#---------------------------------------------------
# CONTRAINTES NOEUD N8359 FONCTION DU DECONFINEMENT
#---------------------------------------------------

SIG_8359=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_FONC_DECONF_N8359',
NOEUD='N8359',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
TOUT_ORDRE='OUI',
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_8359,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('INST','SIXX','SIYY'),);

#-------------------------------------------
# DEPLACEMENTS NO_LEFT2 ===> 50 %
#-------------------------------------------

DEP_L50=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_LEFT2_50%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(500),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_L50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

#-------------------------------------------
# DEPLACEMENTS NO_LEFT2 ===> 100 %
#-------------------------------------------

DEP_L100=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_LEFT2_100%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(1000),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_L100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 59/62


#--------------------------------------------
# CONTRAINTES NO_LEFT2 INSTANT ===> 50 %
#--------------------------------------------

SIG_L50=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_LEFT2_50%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(500),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_L50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#--------------------------------------------
# CONTRAINTES NO_LEFT2 INSTANT ===> 100 %
#--------------------------------------------

SIG_L100=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_LEFT2_100%',
GROUP_NO='NO_LEFT2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(1000),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_L100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#-------------------------------------------
# DEPLACEMENTS NO_BAS2 ===> 50 %
#-------------------------------------------

DEP_B50=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_NO_BAS2_50%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(500),
TOUT_CMP='OUI',
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_B50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);

#-------------------------------------------
# DEPLACEMENTS NO_BAS2 ===> 100 %
#-------------------------------------------

DEP_B100=POST_RELEVE_T(ACTION=_F(INTITULE='DEPL_NO_BAS2_100%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='DEPL',
INST=(1000),
TOUT_CMP='OUI',
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 60/62

OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=DEP_B100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','DX','DY',),);


#--------------------------------------------
# CONTRAINTES NO_BAS2 INSTANT ===> 50 %
#--------------------------------------------

SIG_B50=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_BAS2_50%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(500),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_B50,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

#--------------------------------------------
# CONTRAINTES NO_BAS2 INSTANT ===> 100 %
#--------------------------------------------

SIG_B100=POST_RELEVE_T(ACTION=_F(INTITULE='SIEF_ELNO_ELGA_NO_BAS2_100%',
GROUP_NO='NO_BAS2',
RESULTAT=RESU1,
NOM_CHAM='SIEF_ELNO_ELGA',
INST=(1000),
NOM_CMP=('SIXX','SIYY'),
OPERATION='EXTRACTION',),);

IMPR_TABLE(TABLE=SIG_B100,
FICHIER='RESULTAT',
FORMAT='AGRAF',
NOM_PARA=('NOEUD','COOR_X','COOR_Y','SIXX','SIYY'),);

FIN();
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 61/62


Annexe 7 Comparaison des contraintes obtenues par le calcul
numérique et par la solution analytique

Cas du tunnel non soutenu




Evolution des contraintes selon l'axe vertical


r/R


0
2
4
6
8
10
12
14

0


-2



-4


-6


Solution analytique Contrainte radiale

-8
Solution analytique Contrainte orthoradiale

Contrainte (MPa)

-10
Calcul Code_Aster Contrainte radiale


-12
Calcul Code_Aster Contrainte orthoradiale




Cas du tunnel soutenu (à partir de 50% de déconfinement)

Evolution des contraintes selon l'axe vertical
0
2
4
6
r/R
8
10
12
14
0
Solution analytique Contrainte radiale
Solution analytique Contrainte orthoradiale
-2
Calcul Code_Aster Contrainte radiale
-4
Calcul Code_Aster Contrainte orthoradiale
-6
-8
Contraintes (MPa)
-10


Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Code_Aster ®
Version
6.5

Titre :

Comment creuser un tunnel : méthodologie d'excavation
Date
:

11/06/04
Auteur(s) :
A. COURTOIS, P. SEMETE, A. SAIDANI Clé
:
U2.04.06-A Page
: 62/62




























Page laissée intentionnellement blanche.
Manuel d'Utilisation
Fascicule U2.04 : Mécanique non linéaire
HT-66/04/004/A

Document Outline