Le stockage en formation géologique profonde est l'une des solutions possibles pour la gestion à très long terme des déchets radioactifs de Haute Activité et à Vie Longue (HAVL). Son étude nécessite d'acquérir une bonne connaissance du comportement et des propriétés des formations géologiques potentielles ainsi que de leur évolution au cours du temps sous l'effet des sollicitations induites par une éventuelle installation de stockage.<br />La formation géologique hôte sera soumise à des chargements mécanique et thermique dus respectivement au creusement des ouvrages de stockage et au dégagement de chaleur des colis de déchets stockés. Ces chargements thermomécaniques vont engendrer une modification de la distribution des contraintes dans la couche hôte et les déformations des ouvrages ainsi que l'extension des zones endommagées (EDZ) pourraient provoquer des instabilités locales et globales. <br />Les objectifs de ce travail consistent en l'élaboration de techniques de calculs pour optimiser la modélisation numérique du comportement thermoélastique des ouvrages à grande échelle et en l'évaluation de la perturbation thermomécanique induite par le stockage sur la formation géologique hôte. Dans cette optique, après une présentation de l'état des connaissances sur les aspects thermomécaniques des roches liés au stockage souterrain, des modélisations numériques 2D et 3D du comportement thermoélastique d'ouvrage individuel et d'un réseau de galeries ont été réalisées par une approche discrète. Toutefois, cette approche classique est pénalisante pour l'étude du comportement global d'installations de stockage. Pour pallier à cela, une approche de modélisation numérique, basée sur la technique d'homogénéisation des milieux périodiques, a été proposée. Des formulations ainsi que des procédures numériques ont été élaborées pour déterminer le comportement thermoélastique effectif d'une structure hétérogène équivalente. Le modèle, obtenu par cette méthode, a été validé avec des méthodes d'homogénéisation existantes telles que le modèle auto-cohérent, ainsi que les bornes de Hashin-Shtrikman. La comparaison du comportement thermoélastique effectif avec le comportement thermoélastique réel de référence a montré une bonne cohérence des résultats. Pour une application au stockage profond, les coefficients thermoélastiques effectifs d'un réseau de galeries circulaires ont pu être déterminés en 2D pour différentes dimensions de l'entraxe des galeries.

Dans les massifs rocheux fissurés, il est très difficile de réaliser des ouvrages sans avoir fait, au préalable, une étude de sa fracturation et de son comportement mécanique et hydraulique. L'étude d'un massif rocheux fissuré doit comporter les étapes suivantes : - Relevé systématique des discontinuités du massif, - Étude statistique de ses discontinuités pour identifier les différentes familles qui découpent le massif, - Étude du comportement mécanique des discontinuités, - Étude générale de stabilité, - Étude hydraulique. Notre recherche est relative à l'étude générale de stabilité. Les méthodes classiques de calcul de stabilité, basées sur l'analyse des forces en situation d'équilibre limite et ne prenant pas en compte la déformabilité du massif et la répartition des contraintes existantes dans celui-ci, ne donnent pas entière satisfaction. La méthode que nous avons développée est une méthode dite par "éléments distincts". Elle permet de modéliser les massifs rocheux fissurés par un assemblage de blocs définis par le système de discontinuités et de calculer leurs interactions réciproques. L'interaction entre les blocs se transmet à travers leurs contacts (ou joints) suivant les lois de comportement mécanique de ces derniers. On peut exprimer les contraintes au niveau des joints en fonction de leurs rigidités et des déplacements des blocs. La minimisation de l'énergie potentielle totale de la structure étudiée permet de produire un système d'équations où sont exprimés les forces exercées sur le système et les déplacements des blocs par l'intermédiaire d'une matrice de rigidité. Il est souvent nécessaire de considérer des lois de comportement mécanique des joints non-linéaires, ce qui conduit à un système d'équations non-linéaires dont la résolution se fait d'une manière incrémentale itérative. Par cette méthode, on peut alors étudier les déformations causées par la réalisation d'une excavation dans le massif, rendre compte des ruptures locales qui peuvent évoluer suivant des mécanismes de ruptures progressives. Dans cette recherche, nous avons également modélisé les systèmes d'ancrages par boulons ou tirants actifs et par armatures passives scellées dans le massif. Le comportement mécanique des ancrages est considéré comme élasto-plastique avec écrouissage cinématique.