Laminar flow phenomena may occur when pore water flows at low velocities across the interfaces between soils of different properties, thus causing flow contact resistance. To explore the impacts of interfacial flow contact resistance and rheological characteristics on the thermal consolidation process of layered viscoelastic saturated soil foundation featuring semi-permeable boundaries. This paper established a new thermal consolidation model by introducing a fractional order derivative model, Hagen-Poiseuille law and time-dependent loadings. The semi-analytical solutions for the proposed thermal consolidation model are derived through the Laplace transform and its inverse transform. The reliability and correctness of the solutions are verified with the experimental data in literatures. The influence of constitutive parameters, flow contact resistance model parameters on thermal consolidation process and the interfacial flow contact resistance on foundation settlement, is further explored. The results indicate that the impact of the constitutive parameters and permeability coefficient on the thermal consolidation of viscoelastic saturated soil is related to the flow contact resistance. The enhanced flow contact resistance effect leads to a significant increase in pore water pressure and displacement during the consolidation process.

This study investigates the rheological properties of saturated soft clay surrounding a tunnel using the generalized Voigt viscoelastic model. The model incorporates linear semi-permeability boundary conditions to describe the behavior of the clay. Furthermore, two-dimensional rheological consolidation control equations are derived based on the Terzaghi-Rendulic theory, considering the excess pore water pressure as a variable. To solve the equations, conformal transformation and separation of variables methods are employed, resulting in two independent equations representing the excess pore pressure in terms of time and space variables. The Laplace transformation and partial fractional summation method are then utilized to obtain the solution for excess pore pressure dissipation in the time domain. The reliability of the solution is verified by comparing it with the existing four-element Burgers and five-element model, both of which are derived from the generalized Voigt model. Furthermore, the influence of liner permeability, Kelvin body number, independent Newtonian dashpot viscosity coefficient, and tunnel depth on the dissipation and distribution of excess pore pressure is analyzed based on the established solutions. The findings indicate that a higher relative permeability of the liner and soil leads to an earlier onset of excess pore pressure dissipation and a faster dissipation rate. Increasing the number of Kelvin bodies results in slower dissipation rate. Moreover, larger independent viscous coefficients lead to smaller viscous deformation and faster dissipation rates. Additionally, greater tunnel depth prolongs soil percolation path, slowing down the dissipation of excess pore pressure. When the relative permeability coefficient is 0.01, the excess pore pressure gradually decreases with distance from the outer wall of the tunnel. However, when the relative permeability coefficient is 1, the excess pore pressure initially increases and then decreases with distance. As the relative permeability coefficient increases, the influence of the number of Kelvin bodies on the dissipation of super pore pressure diminishes, the variation in super pore pressure dissipation caused by different independent Newtonian dashpot viscosity coefficients gradually decreases, and the role of tunnel liners as new permeable boundaries within the soil layer is becoming increasingly prominent.