Investigation of the excavation damaged zone around deep TBM tunnel using a Voronoi-element based explicit numerical manifold method

Abstract

A Voronoi-element based explicit numerical manifold method is adopted to investigate the evolution process and failure characteristics of the excavation damaged zone around a deeply buried TBM tunnel. An explicit time integration scheme of numerical manifold method is adopted as the numerical platform, based on which a Voronoi polygon assemblage is generated to approximate the blocky structure of the surrounding rock mass resulting from discrete fracture networks. A new mathematical mesh generation method is developed to better adapt to the Voronoi polygons and an excavation algorithm including both in-situ stress equilibrium and TBM excavation modelling techniques is proposed to capture the mechanical response of the TBM tunnelling process. With the newly developed VE-ENMM approach, the progressive failure process of the surrounding rock mass in a headrace tunnel of Yindajihuang Project is simulated and the formation mechanism of the excavation damaged zone is explored. Simulation results show that a large extent of excavation damaged zone is easily formed in the highly fractured rock mass due to the high in-situ stress in deep grounds. Besides, relatively large convergence displacements are predicted around the tunnel surface, leading to a high risk of TBM shield jamming as in the practical situation. However, the development of damaged zone is effectively restrained and the convergence displacements are well controlled by an early installation of the tunnel lining segments, thus preventing the occurance of TBM shield jamming. Additional simulations are conducted to study the effects of the lateral in-situ stress ratio on the excavation damaged zone. It is inferred that the increase of lateral stress ratio facilitates the fracturing of the surrounding rock mass, which results in a larger extent of the damaged zone and more severe squeezing deformation compared with those under lower lateral stresses.