This computational study focuses on the thermo-hydro-mechanical simulations of the behaviors of freezing soils used for artificial ground freezing (AGF) in a metro project. Leveraging the experimental and field data available in the literature, we simulate the sequential freezing and excavation of a twin tunneling that occurred in months during the actual construction of the tunnel. A thermo-hydro-mechanical model is developed to capture the multi-physical rate-dependent behaviors triggered by phase transitions, as well as the creeping and secondary consolidation of the soil skeleton and the ice crystals. We then calibrate the material models and establish the THM finite element model coupled with the rate-dependent multi-physical models, which may accurately predict the surface heave induced by ground freezing throughout the project. To showcase the potential of using simulations to guide the AGF, we simulate the scenario where a simultaneous freezing scheme is employed as an alternative to the actual sequential scheme design. We then compared the simulated performance with the recorded results obtained from the sequential scheme. Finally, parametric studies on the effect of ground temperature, the porosity of the frozen soil, and the intrinsic elastic modulus of the solid skeleton are conducted. The maximum surface heave is inferred from finite element simulations to quantify the sensitivity and the impact on the safety of AGF operations.

This study presents a hierarchical multiscale approach that combines the finite-element method (FEM) and the discrete-element method (DEM) to investigate tunneling-induced ground responses in coarse-grained soils. The approach considers both particle-scale physical characteristics and engineering-scale boundary value problems (BVPs) simultaneously, accurately reproducing typical tunneling-induced mechanical responses in coarsegrained soils, including soil arching and ground movement characteristics observed in laboratory tests and engineering practice. The study also unveils particle-scale mechanisms responsible for the evolution of soil arching through the underlying DEM-based RVEs. The results show that the rearrangement of microstructures and the deflection of strong contact force chains drive the rotation of macroscopic principal stress and the formation of soil arch. The microscopic fabric anisotropy direction can serve as a quantitative indicator for characterizing soil arching zones. Moreover, the effects of particle size distributions (PSD) and soil densities on ground deformation patterns are interpreted based on the stress-strain responses and contact network characteristics of DEM RVEs. These multiscale insights enrich the knowledge of tunneling-induced ground responses and the same approach can be applied to other geotechnical engineering analyses in coarse-grained soils.

In urban subway construction, shield tunneling near pile groups is common, where additional loads may threaten existing structures. This study establishes multiple 3D nonlinear FDM models with fluid-solid coupling to investigate how tunnel-pile clearances (Hc) affect the mechanical response of low-cap pile groups (2 x2) during side-by-side twin tunneling in composite strata. The advanced CYSoil model, incorporating nonlinearity, strain path dependency, and small strain behavior, is employed to simulate soil response. Results show that tunneling induces up to a similar to 66.7 % reduction in pore water pressure, forming a funnel-shaped seepage pattern. As Hc increases from 0.8D to 2.6D, the low-pressure zone shifts from sidewalls to vault and invert, while maximum displacements reduce by up to 14.04 mm (lateral), 5.28 mm (transverse), and 19.68 mm (vertical). Axial force evolution in piles follows a three-stage decline, i.e., rapid, slow, and moderate, with peak shaft resistance concentrated near the tunnel axis. These findings aid in optimizing tunnel-pile configurations and mitigating geotechnical risks.

Tunneling-induced horizontal strains for buildings with discontinuous foundations are notable and may pose significant risks to the integrity of nearby structures. This paper presents results from a series of numerical models investigating the response of framed buildings on separated footings to tunnel construction in sand. The study examines a two-story, elastic frame with varying building transverse width, eccentricity, and first story height, subjected to tunneling-induced displacements; footing embedment depth and tunnel cover depth are also varied. Results show that tunneling-induced horizontal displacements for separated footings are significant, with greater footing horizontal displacements occurring at deeper footing embedment depths. Building width and eccentricity also influence soil-footing interaction, particularly in determining the values of footing displacements and the distribution of horizontal strains. An increase in footing embedment depth slightly increases shear distortion but significantly increases horizontal strains. The presented modification factors for angular distortion and horizontal strains align well with empirical envelopes, with the horizontal strain modification factor being sensitive to the relative soil-footing stiffness. This research highlights the importance of considering horizontal strains and realistic foundation embedment depth in the damage assessment for buildings with discontinuous foundations due to tunnel construction.

This paper presents an experimental investigation into the interaction mechanism between aqueous foam and unsaturated granite residual soil during conditioning. Contact filter paper tests and undrained shear tests were used to analyze foam's effects on soil water retention and shear behavior, while surface tension tests, capillary rise tests, and microscopic observations examined the role of soil particles in foam stability. The findings demonstrate that foam-conditioned granite residual soils exhibit three distinct saturation- dependent phases (soil-only, transition, and soil-foam mixture) governed by foam's gas-liquid biphasic nature, with foam injection effectively reducing matric suction in unsaturated conditions. Increasing foam injection ratio reduces shear stress while enhancing pore water pressure, with vertical displacement transitioning from contractive to expansive behavior at low shearing rate. Effective cohesion stress varies with gravimetric water content via a rational function, while other effective cohesion stress and friction angles with respect to foam injection ratio, shearing rate, and gravimetric water content obey exponential relationships. The probability distribution function, cumulative distribution function, and decay pattern of bubbles in foam-only systems and soil-foam mixtures all exhibit exponential relationships with elapsed time. Furthermore, a new water-meniscus interaction model was established to characterize rupture and stabilization mechanisms of foam in unsaturated granite residual soils, with particular emphasis on capillary-dominated behavior. Saturation-dependent particle contact modes were identified for foam-conditioned unsaturated granite residual soils, offering valuable guidance for enhancing soil conditioning protocols in earth pressure balance shield tunneling operations.

Shield tunneling inevitably disturbs the surrounding soil, primarily resulting in changes in stress state, stress path, and strain. Modifications to certain parameters, such as shield thrust, shield friction, and soil loss, are made based on the elastic mechanics Mindlin solution and the mirror method, and a calculation expression for additional soil stresses induced by tunneling was derived. Additional soil stresses are calculated using the parameters of the Hangzhou Metro Kanji section. 3D principal stress paths and deviations of the principal stress axes near the tunnel crown, waist, and invert during shield tunneling were obtained by applying a transition matrix orthogonal transformation. These results are compared with experimental data to validate the theoretical solution's accuracy. The stress distribution along the tunneling direction and the 3D principal stress paths and deviations of the principal stress axes in the surrounding soil are determined. The results are as follows: The additional soil stresses along the tunneling direction follow a normal distribution and an S-shape. Under the combined influence of three construction mechanics factors, the shear stress component is approximately 1/3 to 1/2 of the normal stress and should not be neglected. During shield tunneling, the deviation angle of the principal stress axis at the tunnel crown changes from 90 degrees to 180 degrees, with little change in the magnitude of the principal stress. At the invert, the magnitude of the principal stress rapidly increases from 0.25 kPa to 8 kPa, with minimal deviation in the principal stress axis. At the shoulder, the principal stress variation and axis deviation are small. At the foot of the arch, the deviation angle of the major and minor principal stress axes is larger, while the magnitude of the principal stress slightly changes. At the waist, the deviation angle of the major principal stress is larger, and the magnitude of the minor principal stress significantly changes. A strategy for addressing changes in soil stress paths during shield tunnel construction is also proposed.

Shield tunneling in soft soil is prone to losing control of posture and seriously deviating from the designed tunnel axis, which can cause significant ground settlement and tunnel misalignment. This paper develops a large-scale indoor shield tunneling model test platform and conducts shield posture adjustment model tests to examine the interaction between the shield and soft soil. Based on these tests, the distribution and evolution of earth pressure and pore water pressure around the shield, as well as the surface displacement under in-situ posture adjustment and tunneling posture adjustment modes are investigated. The results indicate that the earth pressure around the shield and ground displacement are significantly related to the shield tunneling posture mode. When the shield is in the constant posture angle tunneling mode, the sudden changes in earth pressure and ground displacement are mainly due to the in-situ posture adjustment at a large angle before excavation. However, under the subsequent dynamic tunneling effects of the shield shell approaching and leaving the soil in the excavation, the mutation amount gradually dissipates. When the small-range dynamic posture adjustment mode is adopted, the earth pressure around the shield and surface displacement will change synchronously with the shield posture mode and reach their peak value when the posture angle of the shield is maximum. Moreover, there is a plastic limit in the earth pressure around the shield during the posture adjustment, and under the plastic deformation and flow of the soil, the formation will cause soil loss and ground settlement. These findings offer important information on how the shield interacts with the surrounding soil as well as recommendations for constructing shield tunnels in soft soil.

The shutdown of earth pressure balance (EPB) shield tunneling in gravel stratum can easily lead to significant unexpected ground deformation. In order to study the response of gravel strata during shield shutdown and the characteristic change of soil state in the chamber, this paper establishes a coupled Eulerian-Lagrangian finite element method (CEL-FEM) coupling analysis model that reflects the interaction between the spoiled soil and gravel strata. The plastic flow parameters of CEL spoiled soil are calibrated using the slump method, and a quantitative relationship between the slump value, plastic flow parameters, equivalent coefficient of loosening, and excavation face support pressure is established. The reliability and applicability of CEL method in the simulation of shield shutdown are verified by the field measurements. Results show that: (1) The chamber's soil equivalent loose coefficient is inversely proportional to the soil slump value which is related to soil's plastic flow parameters. (2) The shield shutdown in gravel strata has a more significant impact on the deep strata displacement than on the surface. (3) During the shield shutdown stage, the chamber pressure should be dynamically adjusted based on the soil deformation characteristics, and an increase of 16% could result in a stable rebalance.

The deformation caused by tunnel excavation is quite important for safety, especially when it is adjacent to the existing tunnel. Nevertheless, the investigation of deformation characteristics in overlapped curved shield tunneling remains inadequate. The analytical solution for calculating the deformation of the ground and existing tunnel induced by overlapped curved shield tunneling is derived by the Mirror theory, Mindlin solution and Euler-Bernoulli-Pasternak model, subsequently validated through both finite element simulation and field monitoring. It is determined that the overcutting plays a crucial role in the ground settlement resulting from curved shield tunneling compared to straight shield tunneling. The longitudinal settlement distribution can be categorized into five areas, with the area near the tunnel surface experiencing the most dramatic settlement changes. The deformation of the existing tunnel varies most significantly with turning radius compared to tunnel clearance and grouting pressure, especially when the turning radius is less than 30 times the tunnel diameter. The tunnel crown exhibits larger displacement than the tunnel bottom, resulting in a distinctive 'vertical egg' shape. Furthermore, an optimized overcutting mode is proposed, involving precise control of the extension speed and angular velocity of the overcutting cutter, which effectively mitigates ground deformation, ensuring the protection of the existing tunnel during the construction. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

To investigate the mechanism by which clay shock slurry fills excavation gaps and reduces ground layer deformation during shield tunneling, we conducted a study using the project example of Beijing Metro Line 19 from Youanmenwai Station to Niujie Station, which passes through Guang'anmennei Station to CaiShiKou Station of Beijing Metro Line 7 at a close distance. We employed physical and mechanical testing, numerical simulation calculations, and other methods to examine the deformation law and mechanism of the clay shock method in shield tunneling construction. Our results indicate that (1) as the mass concentration of clay shock slurry increases, its permeability decreases significantly; at a mass concentration of 400 kg/m3, clay shock slurry can prevent synchronous grouting slurry from flowing forward, providing optimal filling and support for excavation gaps. (2) Clay shock slurry can reduce friction between the shield shell and soil body by 50%, avoiding super-consolidation, shear damage, and volumetric expansion of the surrounding soil body. (3) Radial grouting with a two-fluid slurry of cement-water glass at a 1:1 ratio within 15 rings after shield tail removal effectively reduces settlement of the existing tunnel. (4) Numerical simulations demonstrate that using clay shock slurry to fill shield tunnel gaps not only significantly reduces construction settlement but also effectively inhibits strata displacement along the tunnel axis.