Simulation of the fracture of heterogeneous rock masses based on the enriched numerical manifold method

The destruction and fracture of rock masses are crucial components in engineering and there is an increasing demand for the study of the influence of rock mass heterogeneity on the safety of engineering projects. The numerical manifold method (NMM) has a unified solution format for continuous and discontinuous problems. In most NMM studies, material homogeneity has been assumed and despite this simplification, fracture mechanics remain complex and simulations are inefficient because of the complicated topology updating operations that are needed after crack propagation. These operations become computationally expensive especially in the cases of heterogeneous materials. In this study, a heterogeneous model algorithm based on stochastic theory was developed and introduced into the NMM. A new fracture algorithm was developed to simulate the rupture zone. The algorithm was validated for the examples of the four-point shear beam and semi-circular bend. Results show that the algorithm can efficiently simulate the rupture zone of heterogeneous rock masses. Heterogeneity has a powerful effect on the macroscopic failure characteristics and uniaxial compressive strength of rock masses. The peak strength of homogeneous material (with heterogeneity or standard deviation of 0) is 2.4 times that of heterogeneous material (with heterogeneity of 11.0). Moreover, the local distribution of parameter values can affect the configuration of rupture zones in rock masses. The local distribution also influences the peak value on the stress-strain curve and the residual strength. The post-peak stress-strain curve envelope from 60 random calculations can be used as an estimate of the strength of engineering rock masses.