Prise en compte de l'anisotropie dans le comportement instantané des géomatériaux pour les calculs d'ouvrages souterrains

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.