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.

This paper presents a novel analytical solution for the consolidation behavior of viscoelastic saturated soft soil subjected to large-scale ground loading. The rheological properties of clay are described using the general Voigt model. Based on the Terzaghi-Rendulic theory, the governing equations for the dissipation of excess pore water pressure in the surrounding soil mass of a tunnel are established under the first and second boundary conditions. The governing equations are solved using the complex variable method. The obtained solutions are verified by reducing them to the forms of three traditional rheological models, demonstrating the reliability of the proposed approach. Finally, based on the established solutions, the dissipation characteristics of excess pore water pressure around the tunnel are analyzed. The case results indicate that the soil permeability coefficient (k), the independent Newtonian viscosity coefficient (K-0) and Hooke's spring modulus (E-0) in the general Voigt model have significant influences on the dissipation of excess pore water pressure and the degree of consolidation. A larger k, K-0, and E-0 lead to faster dissipation and consolidation development, while a greater tunnel buried depth (b) results in slower consolidation process. The influence of k on excess pore water pressure dissipation is more significant than that of K-0 and E-0. For the first boundary type, consolidation is instantaneous when k > 0.1 m/d. For the second type, when k < 0.0001 m/d, excess pore pressure remains unchanged within 100 d. The permeability condition of the tunnel has a considerable impact on the distribution of excess pore water pressure in the soil layer directly above the crown. When the tunnel is fully permeable, the effects of k, K-0, and E-0 on the dissipation of excess pore water pressure are more pronounced in the early stage, with almost complete dissipation of excess pore pressure above the tunnel before 100 d. However, when the tunnel is completely impermeable, their effects are more prominent in the later stage, and by 100 d, the maximum excess pore pressure within the depth range is 25 % of the initial.