A fully coupled thermo-hydro-mechanical (THM) finite element (FE) formulation is presented that considers freezing and thawing in water-saturated soils. The formulation considers each thermal, hydraulic and mechanical process, and their various interactions, through fundamental physical laws and models. By employing a combination of ice pressure, liquid pressure and total stress as state variables, a new mechanical model has been developed that encompasses frozen and unfrozen behaviour within a unified effective-stress-based framework. Important frozen soil features such as temperature and porosity dependence of shear strength are captured inherently by the model. Potential applications to geotechnics include analysis of frost heave, foundation stability or mass movements in cold regions. The model's performance is demonstrated with reference to the in situ pipeline frost heave tests conducted by Slusarchuk et al. Detailed consideration is given to FE mesh design, the influence of hydraulic parameters, and the treatment of air/ground interface boundary conditions. The THM simulation is shown to reproduce, with fair accuracy, the observed pipeline heave and the porosity growth driven by water migration.KEYWORDS: ground freezing; numerical modelling; temperature effects; water flow Dans la présente communication, on présente une formulation aux éléments finis thermo-hydro-mécanique (THM) à accouplement intégral, pour l'examen du gel et du dégel dans des sols saturés d'eau. La formulation examine chaque processus thermique, hydraulique et mécani-que, et leurs différentes interactions, par le biais de lois et modèles physiques fondamentaux. En employant une combinaison de pression de la glace, de pression liquide, et de contrainte totale comme variables d'état, on dével-oppe un nouveau modèle mécanique incorporant un comportement congelé et non congelé, dans un cadre effectif unifié à base de contrainte. Le modèle permet de capturer d'importantes propriétés des terrains congelés, comme la mesure dans laquelle la résistance au cisaillement est tributaire de la température et de la porosité. Parmi les applications potentielles dans la géotechnique, indiquons l'analyse du foisonnement par le gel, la stabilité des fondations, ou les mouvements de masse dans des régions froides. Les performances du modèle sont démon-trées relativement aux tests de foisonnement par le gel in situ sur le pipeline, effectués par Slusarchuk et al. On examine en détail l'étude maillée aux EF, l'influence de paramètres hydrauliques, ainsi que le traitement des conditions limites de l'interface air/sol. On démontre, avec un bon degré de précision, que la simulation THM reproduit le foisonnement relevé sur le pipeline ainsi que l'augmentation de la porosité déterminée par la migration de l'eau.