A multi-field coupled model contained volumetric strain for unsaturated frozen soil and thermal-hydro-mechanical evolution characteristics of permafrost tunnel
The entrance section of permafrost tunnels in cold regions is particularly vulnerable to frost damage caused by complex thermal-hydro-mechanical (THM) interactions in unsaturated frozen soils. The effects of temperaturedependent volumetric strain variations across different stratum materials on heat and moisture transport are often neglected in existing THM coupling models. In this study, a novel THM coupled model for unsaturated frozen soil integrating volumetric strain correction is proposed, which addresses bidirectional interactions between thermal-hydraulic processes and mechanical responses. The model was validated through laboratory experiments and subsequently applied to the analysis of the Yuximolegai Tunnel. The results indicate that distinct layered ice-water distribution patterns are formed in shallow permafrost under freeze-thaw cycles, driven by bidirectional freezing and water migration. Critical mechanical responses were observed, including a shift in maximum principal stress from the invert (1.40 MPa, frozen state) to the crown (5.76 MPa, thawed state), and periodic lining displacements (crown > invert > sidewalls). Frost damage risks are further quantified by the spatial-temporal zoning of ice-water content-sensitive regions. These findings advance unsaturated frozen soil modeling and provide theoretical guidance for frost-resistant tunnel design in cold regions.