Hydro-mechanical modelling of damage-plasticity couplings in unsaturated geomaterials

Modélisation hydro-mécanique du couplage endommagement-plasticité dans les géomatériaux non saturés Cette thèse a pour but le développement d'un modèle constitutif hydromécanique prenant en compte le couplage de l'endommagement et de la plasticité dans les géomatériaux non-saturés. Le couplage de ce modèle avec les transferts hydrauliques est rendu possible par son implémentation dans le code aux éléments finis Θ-Stock. Un algorithme local a été développé spécifiquement dans ce but. Des problèmes hydromécaniques complètement couplés, tels que la création de la zone endommagée par excavation autour d'un tunnel ou la création d'un endommagement de dessiccation et d'humidification ont été simulés. Une contrainte doublement effective incorporant les effets de la succion et de l'endommagement (considéré isotrope) a été définie en s'appuyant sur des bases thermodynamiques. Cette approche a l'avantage de considérer une unique variable de contrainte étant thermodynamiquement conjuguée aux déformations élastiques. Une formulation hyperélastique dépendante de la pression de confinement est utilisée pour décrire le comportement dans le domaine élastique. L'évolution des rigidités élastiques avec l'endommagement est comparée en considérant deux hypothèses : le principe des déformations équivalentes et le principe de l'énergie élastique équivalente. L'hypothèse d'équivalence des déformations permet d'introduire la contrainte doublement effective dans les équations de plasticité et ainsi de coupler plasticité et endommagement. Les équations de plasticité sont dérivées du modèle de Barcelone (Alonso extit{et al.}, 1990), lui même basé sur le modèle de Cam-Clay. Deux critères distincts d'endommagement et de plasticité sont définis qui peuvent être activés aussi bien indépendamment que simultanément. Les surfaces de charges étant exprimées en fonction de la contrainte effective et de la succion, elles évoluent automatiquement en fonction de la succion et de l'endommagement dans l'espace des contraintes totales. Cela permet de représenter la transition d'un comportement ductile vers un comportement fragile lorsque le matériau est asséché. Un algorithme local explicite a été développé pour gérer le couplage des deux phénomènes dissipatifs. L'implémentation du modèle constitutif dans le code aux éléments finis Θ-Stock permet l'étude de problèmes hydromécaniques complètement couplés, les lois de transfert dépendant elles aussi de l'état de saturation du matériau. Ce modèle a d'abord été appliqué à la simulation du développement de micro-fissures lors du séchage d'une éprouvette. L'apparition de l'endommagement est expliquée principalement par le gradient de pression très important créé à la surface de l'échantillon lorsqu'on applique une forte variation de succion. Enfin, un problème à échelle réelle est simulé. L'excavation d'un tunnel, la désaturation du sol environnant dû à la ventilation, ainsi que la création de la zone endommagée par excavation sont étudiés. L'étendue des zones l'endommagement et de déformations plastiques autour du tunnel est étudiée