Directional fracture patterns of excavated jointed rock mass within rough discrete fractures

Abstract

The influence of the fractures on the anisotropic mechanical behavior of excavated fractured rock mass is important for the stability of rock tunnels. Rough joints were generated using fractal principles, and the Monte Carlo method was employed to create joint network faces. Biaxial compression numerical simulations were conducted using the discrete element method to analyze the behavior of the rock mass containing these networks. Quantitative analysis was carried out to evaluate the strength variability and total number of cracks in the jointed rock mass. A linear regression analysis was conducted to establish the relationship between rock strength and excavation size. The findings demonstrate significant disparities between the linear and fractal models regarding rock damage mode, stress–strain curves, and crack evolution patterns. Tensile rupture and local shear damage were found to be the primary fracture modes for the fractured rock tunnels, forming a V-shaped damage zone within the central tunnel. The shape of this zone was influenced by fracture direction and excavation ratio. Excavation at the center of the rock mass substantially impacted the variability of rock strength and the occurrence of internal cracks. An increase in the slope ratio led to a greater variability of rock mass strength, transitioning from nearly isotropic to weakly anisotropic behavior. The total number of cracks exhibited a W-shaped trend, initially increasing and later decreasing as the slope ratio decreased. Additionally, a clear linear correlation was observed between rock strength and excavation ratio at the center of the rock mass.