Code_Aster TUNNEL Manuel utilisateur
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
1/24
Clé :
U2.04.06
Révision :
13998
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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
2/24
Clé :
U2.04.06
Révision :
13998
Table des matières
1 Comment simuler le creusement d’un tunnel ?
...................................................................................
4.1 Comment définir les modèles à partir d’un maillage simple ?
.......................................................
4.2 Comment initialiser les contraintes ?
...........................................................................................
4.3 Comment calculer les réactions nodales au bord de la « future » galerie ?
................................
4.4 Comment simuler la création d’un « vide » dans le massif et la pose du béton ?
5.2 Cas n° 1 : excavation sans soutènement avec initialisation des contraintes par un calcul et
5.3 Cas n°2 : excavation avec soutènement avec initialisation des contraintes par appel à
5.4 Cas n°3 : excavation avec soutènement avec initialisation des contraintes par appel à
Annexe 2 Organigramme de synthèse sur les méthodes permettant de simuler une excavation .........
Annexe 3 Comparaison des contraintes obtenues par le calcul numérique et par la solution analytique
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
3/24
Clé :
U2.04.06
Révision :
13998
1 Comment simuler le creusement d’un tunnel ?
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 :
1.évaluer la zone endommagée par l’excavation (EDZ) autour d’une galerie de stockage ;
2.é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 œuvre. 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 œuvre. 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
4/24
Clé :
U2.04.06
Révision :
13998
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 [6], [7] ou [4]) 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 œuvre 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, [5]).
Ce rapport présente donc :
1) la méthode classique de simulation d’une excavation souterraine en 2D au moyen d’un code éléments finis ;
2) les différentes option disponibles pour appliquer cette méthode avec
Code_Aster ;
3) 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
5/24
Clé :
U2.04.06
Révision :
13998
3 Une méthode pour simuler le creusement d’une galerie à partir d’un modèle 2D : la méthode convergenceconfinement
3.1
Principe général
Cette partie s’inspire largement de [5]. Signalons que le CIH et TEGG ont également mené un certain nombre d’études avec cette méthode (par exemple, [2]). 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 :
1) déformations planes avec hypothèse de petites perturbations;
2) le tunnel est supposé de section circulaire et d’axe horizontal ;
3) terrain homogène d’extension infinie ;
4) massif suivant un comportement élastique linéaire ou élasto-plastique ;
5) état initial des contraintes supposé isotrope et homogène ;
6) 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 / R10
;
7) é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 [1]). 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
6/24
Clé :
U2.04.06
Révision :
13998
La [Figure 3.1-a] ci-dessous illustre l’évolution de soutenu.
et de la contrainte radiale
R
pour un tunnel non
λ
=1
σ
R
=0
σ
R
0>λ>
1
=(1-
λ).σ
0
σ
λ=0
R
=
σ
0
Figure 3.1-a : Évolution 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
7/24
Clé :
U2.04.06
Révision :
13998
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 œuvre plusieurs modèles (au sens large).
En effet, un déroulement classique de la modélisation peut se résumer par les étapes suivantes :
1) étape 1: initialisation des contraintes in situ ;
2) étape 2 : calcul des réactions nodales au niveau des parois de l’excavation ;
3) étape 3 : déconfinement du massif pour simuler l’excavation progressive et l’éloignement du front de taille ;
4) é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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
8/24
Clé :
U2.04.06
Révision :
13998
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
Initialisation des contraintes
Calcul des réactions nodales
Excavation de la galerie
Pose du revêtement béton
1 2
3
4
Figure 3.2-a : Exemple d’enchaînement typique d’un calcul d’excavation avec un code de calcul
La première configuration sert à :
1.initialiser les contraintes d’origine géostatique ;
2.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 nœuds 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
9/24
Clé :
U2.04.06
Révision :
13998
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
, l’épaisseur de béton est le maillage est un carré de 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 (
L10×R
).
0,30 m
et
20 m
de côté. D’après les règles usuelles de modélisation, le rapport entre
Sol
Sol, vide ou béton
Sol ou vide
Figure 4.1-a : Maillage utilisé et matériaux
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
10/24
Clé :
U2.04.06
Révision :
13998
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) :
1) le bord excavé où seront appliquées les réactions nodales pour simuler le déconfinement (appelé
BORD
) ;
2) 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) :
1) un modèle
2) un modèle
SOL
, dans lequel tout le maillage est affecté d’éléments finis ;
SOL_REST
qui ne comprend pas les mailles qui correspondent à la partie excavée (elles ne sont pas affectées d’éléments finis) ;
3) 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
11/24
Clé :
U2.04.06
Révision :
13998
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
h
=
. z
=
K
0
.
v
avec le poids volumique du terrain sus-jacent (en kN /m
3
par exemple) et
K
0
un coefficient sans dimension.
K
0
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
K
0
en fonction du coefficient de Poisson :
K
0
=
1−
Deux méthodes ont été testées avec Code_Aster pour initialiser les contraintes dans le terrain encaissant :
1.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
K
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é précédent). Dans ce cas,
K
0
K
0
≥
1
(si le sol est soumis à des contraintes tectoniques, par exemple). Dans ce cas, la méthode suivante devient obligatoire ;
0
SOL
dans le chapitre
≤
1
(cas de l’élasticité linéaire). Il existe de nombreux cas où
2.affecter directement les contraintes à tous les éléments du maillage par la commande
CREA_CHAMP (OPERATION = ‘AFFE’) ;
La première solution nécessite de définir un matériau fictif et de mettre en œuvre 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é (se référer à la documentation [U4.72.04] qui donne un exemple adaptable à notre problème).
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
12/24
Clé :
U2.04.06
Révision :
13998
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 nœuds 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 nœuds du bord de la galerie permet d’y évaluer les réactions nodales.
Même chargement qu’à l’étape précédente
Nœuds bloqués
∆
u = 0
Figure 4.3-a : Blocage des nœuds 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 (Document [U4.51.03]). Le blocage de ces nœuds 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 avec la commande CALC_CHAMP (FORCE = ‘REAC_NODA’). Il convient alors de fournir à la commande
CALC_CHAMP 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], 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. Étant 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
13/24
Clé :
U2.04.06
Révision :
13998
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] :
1) 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 ;
2) 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 nœuds 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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
14/24
Clé :
U2.04.06
Révision :
13998
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 :
1.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 ;
2.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 formellement la combinaison linéaire de deux champs :
B−3
de la [Figure 4.4-a] sont
1.les champs issus de l’étape de calcul précédente (
B−2
) et qui ne concernent que le modèle correspondant au massif de sol privé de la zone excavée ;
2.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 obtenus 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 :
1.Méthode
2.Méthode
I
: faire un calcul (appel à STAT_NON_LINE) sur matériau fictif ;
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 :
1.Méthode
A
: qui consiste à affecter des caractéristiques mécaniquement très « souples » dans la zone excavée ;
2.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. Cette méthode sera privilégiée.
Un organigramme synthétique est proposé en [§Annexe2].
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
15/24
Clé :
U2.04.06
Révision :
13998
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). Les cas présentés ici font l'objet des cas test de validation
SSLP105A, SSLP105B et SSLP105C.
Trois cas de calcul sont présentés dans cette partie :
1.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
) ;
2.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
) ;
3.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.
5.1
Le problème traité
La géométrie du maillage est rappelé au paragraphe [§4.1]. Il contient 8477 nœuds et 3304 éléments.
Le rayon de la galerie est de
1,50 m
, l’épaisseur de béton est
0,30 m
et le maillage est un carré de
20 m
de coté. Les autres données sont résumées dans le tableau suivant.
Matériau
Roche
Béton
Paramètre
K
0
v
=
h
E
E
Valeur
1
5 MPa
4 GPa
0,3
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 :
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Pression des terrains sus jacents (5 MPa)
Version 12
Date :
23/10/2015
Page :
16/24
Clé :
U2.04.06
Révision :
13998
Ux = 0 Ux = 0
Déconfinement
Uy = 0
Figure 5.1-a : Conditions aux limites et chargement imposés
à la fin du processus de déconfinement,
=
2. G
0
. R
. U
R
=
0, 69
.
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
K
0
=
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 correspond au cas test SSLP105C.
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 correspond au cas test SSLP105B.
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
17/24
Clé :
U2.04.06
Révision :
13998
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 correspond au cas test SSLP105A.
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
18/24
Clé :
U2.04.06
Révision :
13998
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)
σ
θ
σ
R
U
R
A
B
Figure 6-a : Grandeurs comparées pour la validation de Code_Aster
σ
Variable
U
R
θ
(
r
(
(
Pa
)
Analytique
0.
Pa
)
−
1. E7
Point
−
A
Aster Écart relatif Analytique
−
8.411 E3
On vérifie que
0.
∣
R
∣≪∣
∣
9.883 E6
1,2 % −
1. E7
Aster Écart relatif
−
1.625 E4
On vérifie que
∣
R
∣≪∣
∣
−
Point
B
1.011 E7
m
)
−
0.0024375
−
0.0024772
1.70%
−
0.0024375
−
0.0023982
1,1 %
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
Variable
R
(
Analytique
Point
A
Aster Écart relatif Analytique
Pa
)
−
1.52821 E6
−
1.52974 E6 0,1%
−
1.52821 E6
−
Point
Aster
B
1.52652 E6
U
r
(
(
Pa
)
−
8.47179 E6
−
8.40987 E6
0,7 %
m
)
−
0.0016925
−
0.0017218
1,7 %
−
8.47179 E6
−
8.52586 E6
−
0.0016925
−
0.0016664
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
Écart relatif
0,1 %
0,6 %
1,5 %
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
19/24
Clé :
U2.04.06
Révision :
13998
Variable
R
(
Analytique Aster Écart relatif Analytique
Pa
)
−
1.52821 E6
−
1.52943 E6
0,1 %
Aster
−
1.52821 E6
−
1.53171 E6
U r
(
(
Pa
)
−
8.47179 E6
−
8.40822 E6 0,8 %
m
)
−
0.0016925
−
Point
A
0.0017211
1,7 %
−
8.47179 E6
−
8.52418 E6
−
0.0016925
−
Point
B
0.0016658
Écart relatif
0,2 %
0,6 %
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 scenarii 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 œuvre de plusieurs modèles utilisés successivement peut s’avérer plus fiable du point de vue de la mise en œuvre 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 [3]).
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
20/24
Clé :
U2.04.06
Révision :
13998
8 Bibliographie
[1]
[2]
[3]
[4]
[5]
[6]
[7]
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
P. CATEL : Aval du Cycle – site de Bure – Fiche 13 – Méthode convergence-confinement, note EDF
TEGG EFT GG/00 168 A – 2000
A. GIRAUD : Couplages Thermo-Hydro-Mécaniques dans les milieux poreux peu perméables : application aux argiles profondes, thèse de l’ENPC – 1993
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
M. PANET : Le calcul des tunnels par la méthode convergence-confinement, Presses de l’ENPC –
1995
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
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
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
21/24
Clé :
U2.04.06
Révision :
13998
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
K
0
=
1
).
Contrainte radiale, contrainte orthoradiale et déplacement radial à la paroi du tunnel en milieu élastique soumis à un taux de déconfinement
{
U
=
R
=
R
=⋅
r
R
1−
r
2
. R
1
r
2
. R
2
2
2
⋅
0
2G
.
0
.
0
G
est donné par la relation suivante :
G=
E
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 (
b
est le coefficient de Poisson du béton) :
K
{
=
s
E b
⋅
e
1−
2
b
⋅
R
E b
⋅
R
2
e
−
R i
2
1
b
⋅
[
1−2⋅
b
⋅
R e
2
R i
2
] si R10⋅e si R≤10⋅e
Soit
P s
la pression de confinement définie sur la figure suivante
On a donc :
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
22/24
Clé :
U2.04.06
Révision :
13998 s
=
P s
⋅
R=
d b
⋅
e
Si
k
K s
2⋅G
représente la rigidité relative et le taux de déconfinement à la mise en place du soutènement, alors la pression de soutènement et le déplacement radial en paroi sont donnés par :
U
{
P s
=
R
=
k
1k
s
1
d
⋅
k
1k
s s
⋅
1−
d s
⋅
0
2⋅G
⋅
⋅
R
0
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
23/24
Clé :
U2.04.06
Révision :
13998
Annexe 2 Organigramme de synthèse sur les méthodes permettant de simuler une excavation
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
I
:
SNL1
avec le chargement de poids propre ou de pression souhaité et un
Étape 1 : Initialisation des contraintes
II
: affectation par commande
CC
du champ souhaité matériau doté d’un coefficient de Poisson
éventuellement fictif SOL
Étape 2 : Récupération des réactions nodales au bord de la future galerie
CC
pour extraire les contraintes issus de
SNL1
SNL2
avec
CL
sur l’objet
BORD
en DIDI
Récupération des réactions
SNL1
avec
CL
sur l’objet
BORD
en DIDI sur modèle
SOL
Récupération des réactions
SNL1 l’objet
avec
CL
BORD
en
sur
DIDI sur modèle
SOL_REST
Récupération des réactions
SOL_REST
BORD
SNL3 chargement du vecteur des réactions nodales et un matériau «mou»
à la place du
« vide »
SNL4
avec le
avec 3 matériaux : roche, béton et vide
(méthode
A
) pour achever le déconfinement
SOL_REST
BORD
Eléments mous
Étape 3 : Déconfinement
SNL2
(modèle
SOL
) avec le chargement du vecteur des réactions nodales et un matériau «mou »
à la place du « vide »
SNL2
(un seul matériau et modèle
SOL_REST
) avec le chargement du vecteur des réactions nodales
SOL_REST
BORD
SOL_REST
BORD
Eléments mous
Étape 4 : Pose du soutènement
SNL3
avec 3 matériaux : roche, béton et vide
(méthode
A
) pour achever le déconfinement
CC
pour extraire les résultats de
SNL2
SNL3
sur modèle
SOL_REST
+ BETON ) pour calcul intermédiaire
Combinaison des champs
CC
SNL4
pour achever le déconfinement
BETON
SOL_REST
BORD
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable :
Sylvie GRANET
Version 12
Date :
23/10/2015
Page :
24/24
Clé :
U2.04.06
Révision :
13998
Annexe 3 Comparaison des contraintes obtenues par le calcul numérique et par la solution analytique
Cas du tunnel non soutenu
-10
-12
-6
-8
-2
-4
0
0
Evolution des contraintes selon l'axe vertical
2 4 6
r/R
8 10 12 14
Solution analytique Contrainte radiale
Solution analytique Contrainte orthoradiale
Calcul Code_Aster Contrainte radiale
Calcul Code_Aster Contrainte orthoradiale
Cas du tunnel soutenu (à partir de 50% de déconfinement)
Manuel d'utilisation Fascicule u2.04 : Mécanique non linéaire
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
