Code_Aster TUNNEL Manuel utilisateur

Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 1/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 2/64
Clé : U2.04.06
Révision : 1283
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 »........................................................................... 17
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........................................................... 17
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........................................................... 17
6 Validation du Code_Aster sur un exemple d’excavation en milieu élastique linéaire........................ 18
7 En guise de conclusion : conseils et perspectives............................................................................ 19
8 Bibliographie..................................................................................................................................... 20
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........................................................ 21
Annexe 2 Organigramme de synthèse sur les méthodes permettant de simuler une excavation dans le
Code_Aster.................................................................................................................................. 23
Annexe 3 Fichier de maillage réalisé avec GIBI.................................................................................. 24
Annexe 4 Excavation sans soutènement, sur la base d’un seul modèle (cas n°1). Fichier de
commandes Code_Aster............................................................................................................. 27
Annexe 5 Excavation avec soutènement, méthode A (cas n°2). Fichier de commandes Code_Aster 38
Annexe 6 Excavation avec soutènement, méthode B (cas n°3). Fichier de commandes Code_Aster 49
Annexe 7 Comparaison des contraintes obtenues par le calcul numérique et par la solution analytique
62
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Date : 07/05/2009 Page : 3/64
Clé : U2.04.06
Révision : 1283
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 4/64
Clé : U2.04.06
Révision : 1283
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)
2)
é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 œ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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
2
Date : 07/05/2009 Page : 5/64
Clé : U2.04.06
Révision : 1283
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 œ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, [bib5]).
Ce rapport présente donc :
1)
2)
3)
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 6/64
Clé : U2.04.06
Révision : 1283
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 [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 :
1)
2)
3)
4)
5)
6)
7)
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
une fraction de la contrainte initiale  0 :
 d ans le terrain, qui est
= 1−  . 0 avec 0≤≤1
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
La [Figure 3.1-a] ci-dessous illustre l’évolution de
non soutenu.
λ=1
σR=0
Date : 07/05/2009 Page : 7/64
Clé : U2.04.06
Révision : 1283
 et de la contrainte radiale  R pour un tunnel
0>λ>1
σR=(1-λ).σ0
λ=0
σR=σ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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 8/64
Clé : U2.04.06
Révision : 1283
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)
2)
3)
4)
é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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 9/64
Clé : U2.04.06
Révision : 1283
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
Excavation de la
galerie
Pose du
revêtement béton
Calcul des réactions
nodales
1
2
3
4
Figure 3.2-a.2-a : Exemple d’enchaînement typique d’un calcul d’excavation avec un code de calcul
La première configuration sert à :
1)
2)
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 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
4.1
Date : 07/05/2009 Page : 10/64
Clé : U2.04.06
Révision : 1283
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 4.1-a.1-a : Maillage utilisé et matériaux
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 11/64
Clé : U2.04.06
Révision : 1283
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)
2)
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) :
1)
2)
3)
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
Tunnel
Partie maillée
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 12/64
Clé : U2.04.06
Révision : 1283
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 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 le Code_Aster pour initialiser les contraintes dans le terrain
encaissant :
1)
2)
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 0 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,
K 0  1 (cas de l’élasticité linéaire). Il existe de nombreux cas où K 0 
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 œ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é (le
paragraphe [§3.5.3.1] de la documentation [U4.72.04] indice B1 donne un exemple adaptable à notre
problème).
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
4.3
Date : 07/05/2009 Page : 13/64
Clé : U2.04.06
Révision : 1283
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 (Doc. Aster [U4.51.03] indice F4,
paragraphe [§3.2.2]). 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 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
4.4
Date : 07/05/2009 Page : 14/64
Clé : U2.04.06
Révision : 1283
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)
2)
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 nœuds doubles aux frontières entre matériaux qui permettent de lier ou
pas les deux structures).
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
4.4.1
Date : 07/05/2009 Page : 15/64
Clé : U2.04.06
Révision : 1283
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)
2)
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 :
1)
2)
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 :
1)
2)
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 :
1)
2)
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
5
Date : 07/05/2009 Page : 16/64
Clé : U2.04.06
Révision : 1283
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 :
1)
2)
3)
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 œuvre 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 nœuds 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
Roche
K0
 v = h
Béton
Valeur
1
5 MPa
E
4 GPa

E
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 17/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
A la fin du processus de déconfinement,
=
Date : 07/05/2009 Page : 18/64
Clé : U2.04.06
Révision : 1283
2. G
. U =0, 69 .
0 . R 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
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 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
6
Date : 07/05/2009 Page : 19/64
Clé : U2.04.06
Révision : 1283
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
A
UR
B
Figure 6-a : Grandeurs comparées pour la validation du Code_Aster
Variable
σ R (Pa)
Analytique
0.
Point A
Aster
-8.411 E3
Ecart relatif
On vérifie
que
Analytique
0.
Point B
Aster
-1.625 E4
∣ R∣≪∣ ∣
σ θ (Pa)
- 1. E7
-9.883 E6
U r (m)
- 0.0024375
- 0.0024772
1,2 %
1,7 %
Ecart relatif
On
vérifie
que
∣ R∣≪∣ ∣
- 1. E7
-1.011 E7
- 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
Point A
Point B
Analytique
Aster
Ecart relatif Analytique
Aster
Ecart relatif
- 1.52821
-1.52974 E6
0,1%
- 1.52821 -1.52652 E6
0,1 %
E6
E6
- 8.47179
-8.40987 E6
0,7 %
- 8.47179 -8.52586 E6 0,6 %
  (Pa)
E6
E6
U r (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
Variable
σ R (Pa)
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 20/64
Clé : U2.04.06
Révision : 1283
Point A
Point B
Analytique
Aster
Ecart relatif Analytique
Aster
Ecart relatif
- 1.52821
- 1.52943
0,1 %
- 1.52821
- 1.53171
0,2 %
E6
E6
E6
E6
- 8.47179
-8.40822 E6
0,8 %
- 8.47179
- 8.52418 0,6 %
  (Pa)
E6
E6
E6
U r (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
Variable
 R  Pa)
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 œ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 [bib3]).
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
8
Date : 07/05/2009 Page : 21/64
Clé : U2.04.06
Révision : 1283
Bibliographie
•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
Code_Aster, rapport Ajilon Technologies Cénergys 01-A – 2002
avec le
•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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 22/64
Clé : U2.04.06
Révision : 1283
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 λ
{


 R= 1−


 . R2
. 0
2
r
. R2
. 0
2
r
R2  0
U R =⋅ ⋅
r 2G
 = 1
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 ( v b est le coefficient de Poisson du béton) :
K s=
Soit
{
E b⋅e
 1−2b ⋅R
E b⋅ R2e −R2i 
 1b ⋅[  1−2⋅b ⋅Re2R2i ]
si R10⋅e
si R≤10⋅e
P s la pression de confinement définie sur la figure suivante
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
On a donc :
Si
k s=
Ks
2⋅G
Date : 07/05/2009 Page : 23/64
Clé : U2.04.06
Révision : 1283
P s⋅R= b⋅e
représente la rigidité relative et d
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 :
{
ks
⋅ 1−d ⋅ 0
1k s
1d⋅k s  0
U R=
⋅
⋅R
1k s
2⋅G
P s=
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 24/64
Clé : U2.04.06
Révision : 1283
Annexe 2 Organigramme de synthèse sur
permettant de simuler une excavation
les
méthodes
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
SNL1
l’objet BORD en DIDI
SNL 2 avec CL sur l’objet BORD en DIDI
sur modèle SOL
Récupération des réactions
Récupération
des
réactions
SNL 1 avec CL sur l’objet
BORD en DIDI sur
modèle SOL_REST
Récupération
des
réactions
SOL_REST
BORD
Etape 3 : Déconfinement
SNL 3 avec le
chargement
du
vecteur
des
réactions nodales et
un matériau «mou»
à
la
place
du
« vide »
SOL_REST
BORD
SNL 2 (modèle SOL)
avec le chargement
du
vecteur
des
réactions nodales et
un matériau «mou »
à la place du « vide »
SNL 2 (un seul matériau
et modèle SOL_REST )
avec le chargement du
vecteur des réactions
nodales
SOL_REST
BORD
Eléments mous
Etape 4 : Pose du soutènement
SNL 4 avec 3
matériaux :
roche,
béton
et
vide
(méthode A) pour
achever
le
déconfinement
SOL_REST
BORD
SNL 3 avec 3
matériaux :
roche,
béton
et
vide
(méthode A) pour
achever
le
déconfinement
Eléments mous
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
CC pour extraire les
résultats de SNL 2
SNL 3 sur modèle
SOL_REST + BETON )
pour calcul intermédiare
Combianson des champs
CC
SNL 4 pour achever le
déconfinement
SOL_REST
BORD
BETON
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 25/64
Clé : U2.04.06
Révision : 1283
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
E2
E3
E4
=
=
=
=
0. 0. ;
0.9 0. ;
0.7 0.7 ;
0. 0.9 ;
B1 = 1.2 0. ;
B2 = 0. 1.2 ;
S1
S2
S3
S4
S5
S6
=
=
=
=
=
=
1.5 0. ;
20. 0. ;
20. 20. ;
0. 20. ;
0. 1.5 ;
(1.5 * (SIN 45)) (1.5 * (COS 45)) ;
********************************************************
*
DROITES
*
********************************************************
E1E2
E2E3
E3E4
E4E1
=
=
=
=
E1
E2
E3
E4
DROI
DROI
DROI
DROI
16
16
16
16
B1S1
S5B2
E2B1
B2E4
=
=
=
=
B1
S5
E2
B2
DROI
DROI
DROI
DROI
4
4
4
4
S1S2 =
S2S3 =
S3S4 =
S2S3S4
S4S5 =
S3S6 =
E2
E3
E4
E1
S1
B2
B1
E4
;
;
;
;
;
;
;
;
S1 DROI -70 S2 DINI 0.01 DFIN 0.50 ;
S2 DROI 16 S3 ;
S3 DROI 16 S4 ;
= S2S3 ET S3S4 ;
S4 DROI -70 S5 DINI 0.50 DFIN 0.01 ;
S3 DROI -70 S6 DINI 0.70 DFIN 0.001 ;
********************************************************
*
ARCS
*
********************************************************
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
S1S5
S5S6
S6S1
B2B1
=
=
=
=
32
16
16
32
CERC
CERC
CERC
CERC
S1
S5
S6
B2
E1
E1
E1
E1
S5
S6
S1
B1
Date : 07/05/2009 Page : 26/64
Clé : U2.04.06
Révision : 1283
;
;
;
;
********************************************************
*
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Version
default
Date : 07/05/2009 Page : 27/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 28/64
Clé : U2.04.06
Révision : 1283
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'),
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 29/64
Clé : U2.04.06
Révision : 1283
_F(GROUP_MA='NO_LEFT2'),
_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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Version
default
Date : 07/05/2009 Page : 30/64
Clé : U2.04.06
Révision : 1283
#########################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 31/64
Clé : U2.04.06
Révision : 1283
#############################################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 32/64
Clé : U2.04.06
Révision : 1283
##############################################################
# 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,
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 33/64
Clé : U2.04.06
Révision : 1283
ITER_INTE_MAXI=50,
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 34/64
Clé : U2.04.06
Révision : 1283
###########################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 35/64
Clé : U2.04.06
Révision : 1283
#############################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 36/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 37/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 38/64
Clé : U2.04.06
Révision : 1283
#-------------------------------------------# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 39/64
Clé : U2.04.06
Révision : 1283
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'),
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 40/64
Clé : U2.04.06
Révision : 1283
_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 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,
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 41/64
Clé : U2.04.06
Révision : 1283
ALPHA=0.0,),);
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Version
default
Date : 07/05/2009 Page : 42/64
Clé : U2.04.06
Révision : 1283
#########################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 43/64
Clé : U2.04.06
Révision : 1283
######################################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 44/64
Clé : U2.04.06
Révision : 1283
##########################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 45/64
Clé : U2.04.06
Révision : 1283
########################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 46/64
Clé : U2.04.06
Révision : 1283
############################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 47/64
Clé : U2.04.06
Révision : 1283
#############################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 48/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 49/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 50/64
Clé : U2.04.06
Révision : 1283
#-------------------------------------------# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 51/64
Clé : U2.04.06
Révision : 1283
Annexe 6Excavation 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'),
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 52/64
Clé : U2.04.06
Révision : 1283
_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',),),);
####################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 53/64
Clé : U2.04.06
Révision : 1283
#########################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Version
default
Date : 07/05/2009 Page : 54/64
Clé : U2.04.06
Révision : 1283
########################################################
#
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 55/64
Clé : U2.04.06
Révision : 1283
##############################################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 56/64
Clé : U2.04.06
Révision : 1283
################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 57/64
Clé : U2.04.06
Révision : 1283
########################################################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 58/64
Clé : U2.04.06
Révision : 1283
###########################
# 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',
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Version
default
Date : 07/05/2009 Page : 59/64
Clé : U2.04.06
Révision : 1283
COEF_R=0.),
_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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 60/64
Clé : U2.04.06
Révision : 1283
##############################################
# 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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 61/64
Clé : U2.04.06
Révision : 1283
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 62/64
Clé : U2.04.06
Révision : 1283
#-------------------------------------------# 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),
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 63/64
Clé : U2.04.06
Révision : 1283
TOUT_CMP='OUI',
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
Version
default
Code_Aster
Titre : Comment creuser un tunnel : méthodologie d’excavat[...]
Responsable : Sylvie GRANET
Date : 07/05/2009 Page : 64/64
Clé : U2.04.06
Révision : 1283
Annexe 7Comparaison 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
0
2
4
r/R
6
8
10
12
14
Contrainte (MPa)
0
-2
-4
-6
Solution analytique Contrainte radiale
-8
-10
-12
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
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.04 : Mécanique non linéaire
">