Numerical Investigation of Rock–Backfill Composite Fracture Evolution Laws Under Deep Mining and Filling Stress Paths

Abstract

Fracture evolution of rock–backfill composites (RBCs) under complex loading–unloading and dynamic disturbances is critical for stope stability in deep backfill mining. Using PFC3D, this study constructs numerical models of RBCs to investigate this process, considering burial depths (500, 1000, 1500, 2000 m), interface angles (IA: 60°, 90°), and cement–tailings ratios (CTR—1:4, 1:8), while replicating true triaxial paths and blasting impacts. Systematic analysis of mesoscopic crack quantity, spatiotemporal distribution, and multiscale fracturing reveals that shear cracks dominate damage, with crack counts evolving in stages as strain increases. With greater depth, the number of propagation stages and growth rate inflection points shift systematically. During mining–filling disturbance, crack quantity negatively correlates with depth but turns positive during late static loading beyond 70% peak stress. Spatial crack distribution is synergistically controlled by IA, CTR, and depth. For IA 60°, shear crack angles spread broadly yet concentrate at 50–70°; for IA 90°, they are near-axial, concentrated at 80–90°. The synergistic process progresses through microscopic initiation, mesoscopic accumulation, and macroscopic instability. In terms of failure modes, IA 60° exhibits shear failure along the cemented interface plus tensile fracturing in rock, while IA 90° shows combined diagonal shear and axial tension. Higher CTR yields more extensive fracture networks in backfill, indicating superior synergistic bearing capacity.