L'utilisation d'argile gonflante pour les barrières ouvragées des ouvrages de stockage des déchets nucléaires, nécessite la maîtrise de son comportement thermohydromécanique. Le travail présenté ici concerne la caractérisation et la modélisation du comportement de lune de ces argiles : largile FoCa 7.<br />Les caractéristiques du matériau étudiées sont : la courbe de sorption (/ désorption) et la perméabilité à leau. Pour chacune delles, de nouveaux protocoles expérimentaux ont permis lacquisition de données dans des domaines encore peu explorés : en température (entre 20oC et 80oC) pour la courbe de sorption et dans le domaine non saturé pour la perméabilité à leau.<br />Lanalyse de ces résultats et de données bibliographiques a permis dune part dévaluer limportance du phénomène dhystérésis et de linfluence de la température sur la courbe de sorption, et dautre part, détablir la nécessité dintroduire au niveau de la modélisation de la courbe de sorption, un paramètre plastique lié aux déformations irréversibles se produisant lors du compactage.<br />De plus, les essais réalisés pour lacquisition des données ont également pu être utilisés afin de fournir des éléments de validation aux lois de comportement non linéaires proposées par O. Coussy et P. Dangla pour les milieux poreux non saturés. La particularité de ces lois est de supposer lexistence dune contrainte effective dans le domaine non saturé, doù limportance des éléments de validation présentés ici.

L'étude porte sur la recherche systématique de solutions analytiques et numériques de l'équation de l'infiltration monodimensionnelle de l'eau dans le sol. Dans une premièere partie les différents approches quasianalytiques sont analysées. La solution en série PHILIP, retenue dans cette étude est améliorée et complétée. Elles est ensuite utilisée comme base de référence pour tester les solutions numériques. Dans une seconde partie, un grand nombre de modèles numériques sont comparés afin d'essayer de dégager un ou plusieurs schémas présentant un optimum vis à vis des critères de convergence , précision et coût économique.

Dans les concepts actuels de stockage des déchets radioactifs en formations géologiques, le conteneur métallique emballant les colis se retrouvera en contact d'un milieu argileux totalement ou partiellement saturé. Après consommation totale de l'oxygène piégée, la corrosion du conteneur entame une phase anoxique en produisant de l'hydrogène. Dans la logique d'une évaluation des performances du stockage, deux problématiques essentielles sont abordées par la présente étude: la prédiction de la cinétique de corrosion anoxique à long terme en fonction de l'état hydrique de l'argile de voisinage et l'évaluation du risque d'endommagement hydro-mécanique de l'argile lié à la production d'hydrogène. Les mécanismes qui contrôlent la cinétique de la corrosion se limitent à la croissance de la couche d'oxyde et au transport à travers la porosité de ce film. Une description macroscopique de ces mécanismes couplée avec une loi d'interface a permis de proposer un modèle mécanistique de la corrosion anoxique dans un milieu poreux non saturé. Cette démarche est validée à partir des essais réalisés dans des argiles saturées sur une durée décennale. La validation est ensuite étendue vers les argiles partiellement saturées à partir des essais réalisés en laboratoire. Ces derniers ont confirmé la croissance de la vitesse de corrosion anoxique avec l'humidité relative comme dans le cas de la corrosion aérée ou atmosphérique. Concernant la problématique gaz, des simulations de l'essai oedométrique et triaxial MEGAS sont menées moyennant des choix ciblés sur les paramètres de transport dans l'argile. Une évaluation de l'impact de la corrosion est ensuite appliquée à un concept de stockage géologique. La maîtrise du terme source gaz a permis de minimiser le risque de surpressions gazeuses, de désaturation locale et des déformations mécaniques à long terme.

Afin de mieux estimer les tassements et la stabilité des ouvrages géotechniques au cours du temps, une modélisation du comportement couplé (mécanique, écoulement de l'eau et écoulement de l'air) des milieux poreux non saturés est développée. A partir des équations d'équilibre, de continuité et de comportement des trois phases du sol, nous avons d'abord établi un système non linéaire tant du point de vue mécanique que du point de vue des écoulements. Ce modèle est discrétisé dans l'espace, par éléments finis, et dans le temps, par un schéma implicite. Le traitement numérique de ce problème conduit, d'une part, à la création d'une nouvelle famille d'éléments de consolidation bidimensionnelle des sols non saturés dans le code de calcul par éléments finis CESAR-LCPC et, d'autre part, à la mise au point d'un algorithme itératif combinant la méthode des rigidités variables pour le traitement des non linéarités liées à l'écoulement des fluides et la méthode des contraintes initiales pour le traitement des non linéarités mécaniques. La programmation de cet algorithme a été effectuée dans le module CSNS (Consolidation des Sols Non Saturés) de CESAR-LCPC. La modélisation effectuée et sa programmation, qui tient compte, entre autres, du phasage de la construction du remblai et des conditions aux limites de diverses natures, ont été testées sur un ouvrage réel (remblai autoroutier). La confrontation des résultats numériques et expérimentaux sur cet ouvrage a montré que cette approche de couplage des écoulements de l'eau et de l'air et des déformations du sol produit des résultats encourageants tant du point de vue mécanique que du point de vue des écoulements.