Numerical Simulation of Failure Behavior of Brittle Heterogeneous Rock under Uniaxial Compression Test

Rocks have formed heterogeneous characteristics after experiencing complex natural geological processes. Studying the heterogeneity of rocks is significant for rock mechanics. In this study, a linear parallel bond model with Weibull distribution in two-dimensional particle flow code (PFC2D) is adopted to study the mechanical characteristics and brittle failure mode of granite rock specimens with different heterogeneity. Firstly, we selected several combinations of key micro-parameters of the parallel bond model. Then, we subjected them to a Weibull distribution to satisfy heterogeneity, respectively. Finally, we chose one optimal combination plan after comparing the stress–strain curves of heterogeneous rock specimens. We analyzed the simulated results of heterogeneous rock specimens. The crack distribution of rock specimens under peak stress shows different characteristics: a diagonal shape in rock specimens with low heterogeneity indexes, or a rotated “y” shape in rock specimens with high heterogeneity indexes. As for failure mode, the numerical simulation results show high consistency with the laboratory experiment results. The rock specimen breaks down almost diagonally, and the whole specimen tends to form an x-shaped conjugate shear failure or the well-known “hour-glass” failure mode. With the increase of the homogeneity index of the rock specimen, the shear rupture angle becomes larger and larger. Generally, the crack number increases with time, and when the rock specimen reaches the peak failure point, the number of cracks increases sharply. The development of cracks in numerical rock specimens under compression test is a result of the coalescence of many microscopic cracks. Furthermore, tensile cracks formed initially, followed by shear behavior along the macroscopic crack plane. We also preliminarily study the mechanical characteristics of heterogeneous rock specimens with discontinuous structural planes. The discontinuous structural planes are simulated by the smooth-joint model. We can conclude that the discontinuous structural planes and the microscopic structural planes which contribute to the heterogeneity have a mutual influence on each other.