Consideration of anisotropy in the instantaneous behaviour of geomaterials for underground structure calculations
Prise en compte de l'anisotropie dans le comportement instantané des géomatériaux pour les calculs d'ouvrages souterrains
Abstract
Within the framework of radioactive waste management, for long-lived and high-level radioactive waste, France has chosen deep geological storage. This storage will be located in a rock formation, in the east of the basin Parisian, in a tectonically stable area. This Callovo-Oxfordian claystone formation was chosen for its radionuclides retention properties since it is not very permeable. An underground research laboratory was built 500 meters deep to study the properties of the rock in-situ. The Callovo-Oxfordian claystone has been shown to exhibit anisotropy related to its geological formation. The main objective of this PhD thesis is to reproduce anisotropic mechanical behaviour using an elasto-plastic constitutive model. The thermo-elasto-viscoplastic and isotropic behaviour model named LKR is the result of all EDF’s knowledge and expertise in the design of underground structures. Therefore, we seek to apply a method of taking into account structural or inherent anisotropy characterizing several types of geomaterials including Callovo-Oxfordian claystone to this LKR constitutive model. In order to achieve this objective, two methods of taking anisotropy into account and developed in the scientific literature are applied to a Drucker-Prager model with linear softening and then are compared to each other. The first method consists of introducing a microstructure or fabric tensor to define the preferential orientations of the material. The second method is the weakness plane approach which consists of describing the anisotropic behaviour of the material through two distinct mechanisms, one describing the behaviour of the isotropic rock matrix and the other describing the weakness planes. It is through this second mechanism that anisotropy is introduced. These two applications were made using the free mechanical simulation software, Code Aster, developed by EDF. This allowed to understand the numerical difficulties of each of these methods, and to choose the most relevant approach for the extension of the LKR model. Thus, it is the method with the fabric tensor which is, in our case, the most suitable. It was therefore applied to the LKR model. This new extension to the model enables the material orientation dependency of the compression strengths to be taken into account. In this work, the mechanics of continuous media are considered. Anisotropy is introduced only into the elasto-plastic mechanism of the LKR behaviour model.
Dans le cadre de la gestion des déchets radioactifs dits de moyenne activité à vie longue et de haute activité, la France a fait le choix du stockage géologique profond dans une formation rocheuse située à l’est du bassin parisien, dans une zone tectoniquement stable. Cette formation d’argilites du Callovo-Oxfordien a été choisie pour ses propriétés de rétention des radionucléides puisqu’elle est très peu perméable. Un laboratoire de recherche souterrain a été construit à environ 500 mètres de profondeur afin d’étudier in-situ les propriétés de la roche. Il a été montré que l’argilite du Callovo-Oxfordien présente une anisotropie liée à sa formation géologique. L’objectif principal de cette thèse est de pouvoir reproduire le comportement mécanique anisotrope à l’aide d’un modèle de comportement élasto-plastique. Le modèle de comportement thermo-élasto-viscoplastique et isotrope nommé LKR est la résultante de tout le savoir-faire et l’expertise d’EDF en ce qui concerne le dimensionnement d’ouvrages souterrains. On cherche donc à appliquer une méthode de prise en compte de l’anisotropie structurelle caractérisant plusieurs types de géomatériaux dont l’argilite du Callovo-Oxfordien à ce modèle de comportement. Afin d’atteindre cet objectif, deux méthodes de prise en compte de l’anisotropie développées dans la littérature scientifique sont appliquées à un modèle de Drucker-Prager à écrouissage linéaire négatif et sont comparées. La première méthode consiste à introduire un tenseur de microstructure ou de fabrique permettant de définir les orientations préférentielles du matériau. La seconde méthode est l’approche par plan de faiblesse qui consiste à décrire le comportement anisotrope du matériau via deux mécanismes distincts, l’un décrivant le comportement de la matrice rocheuse isotrope et l’autre décrivant les plans de faiblesse. C’est par ce deuxième mécanisme que l’anisotropie est introduite. Ces deux applications ont été faites dans le logiciel libre de simulation en mécanique, Code Aster, développé par EDF et ont permis d’appréhender les difficultés numériques de chacune de ces méthodes, et de choisir l’approche la plus pertinente pour l’extension du modèle LKR. Ainsi, c’est la méthode avec le tenseur de fabrique qui est, dans notre cas, la plus adaptée. Elle a donc été appliquée au modèle LKR. Cette nouvelle extension au modèle permet de prendre en compte la dépendance à l’orientation du matériau des résistances en compression. Dans ces travaux, on se place dans le cadre de la mécanique des milieux continus. L’anisotropie est introduite seulement dans le mécanisme élasto-plastique du modèle de comportement LKR.
Origin | Files produced by the author(s) |
---|
Loading...