Coupled FDM–DEM Method for Analyzing EPBS Machine Tunneling Performance in Boulders

This work focuses on the development of a three-dimensional shield tunneling model coupling the finite-difference method (FDM) with the discrete-element method (DEM) and evaluates its tunneling performance in boulders. The tunneling process of cutterhead rotation, cutting the formation, soil entering the chamber, and discharging of the screw conveyor is simulated to explain the characteristics of particle movement in front of the tunnel face. The range of the soil failure zone under different shield discharge ratios is determined by the correlation between the discharge ratio and the tunnel surface stability. The evolution process of the soil arching effect in the boulders along with shield tunneling is investigated to understand the mechanisms of the delay settlement in the boulders. Subsequently, the shield discharge ratio is optimized by investigating its influence on the surface deformation. Finally, the average stress distribution on the shield cutterhead in the tunneling process is analyzed to provide the scheme of the cutterhead structure and arrangement. This work can be referred for evaluating the shield tunneling performance and optimizing the excavation design scheme of shield machines during tunneling construction under similar geological conditions.