Cette étude porte sur la modélisation numérique du creusement de tunnels à l'aide de boucliers à front pressurisé et notamment sur l'adaptation de la méthode du ramollissement, initialement mise au point pour le calcul, en coupe bidimensionnelle, de tunnels creusés à l'aide de procédés non mécanisés. La thèse est composée de 8 chapitres et comprend une étude bibliographique concernant le calcul des ouvrages souterrains, une analyse de la réponse du terrain à une pressurisation du front de taille à partir des calculs axisymétriques, la mise au point d'une méthode permettant de reproduire, dans un calcul plan, en coupe transversale, les effets mis en évidence en conditions axisymétriques et l'application de cette méthode à l'analyse d'un ouvrage instrumenté dans le département du Val-de-Marne. Cet ouvrage a été équipé d'inclinomètres et d'extensomètres installés dans le terrain, ainsi que de témoins à cordes vibrantes disposés dans un anneau de revêtement. Les mesures ont permis d'analyser la réponse du massif encaissant et de l'ouvrage de soutènement au cours des différentes étapes de creusement. Le creusement de cet ouvrage a été simulé à l'aide de la méthode de calcul développée dans le cadre de ce travail et les résultats de calculs confrontés à ceux donnés par les mesures.

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.