Numerical Simulation of Fracture Failure Propagation in Water-Saturated Sandstone with Pore Defects Under Non-Uniform Loading Effects

The instability of mine roadways is significantly influenced by the coupled effects of groundwater seepage and non-uniform loading. These interactions often induce localized plastic deformation and progressive failure, particularly in the roof and sidewall regions. Seepage elevates pore water pressure and deteriorates the mechanical properties of the rock mass, while non-uniform loading leads to stress concentration. The combined effect facilitates the propagation of microcracks and the formation of shear zones, ultimately resulting in localized instability. This initial damage disrupts the mechanical equilibrium and can evolve into severe geohazards, including roof collapse, water inrush, and rockburst. Therefore, understanding the damage and failure mechanisms of mine roadways at the mesoscale, under the combined influence of stress heterogeneity and hydraulic weakening, is of critical importance based on laboratory experiments and numerical simulations. However, the large scale of in situ roadway structures imposes significant constraints on full-scale physical modeling due to limitations in laboratory space and loading capacity. To address these challenges, a straight-wall circular arch roadway was adopted as the geometric prototype, with a total height of 4 m (2 m for the straight wall and 2 m for the arch), a base width of 4 m, and an arch radius of 2 m. Scaled physical models were fabricated based on geometric similarity principles, using defect-bearing sandstone specimens with dimensions of 100 mm × 30 mm × 100 mm (length × width × height) and pore-type defects measuring 40 mm × 20 mm × 20 mm (base × wall height × arch radius), to replicate the stress distribution and deformation behavior of the prototype. Uniaxial compression tests on water-saturated sandstone specimens were performed using a TAW-2000 electro-hydraulic servo testing system. The failure process was continuously monitored through acoustic emission (AE) techniques and static strain acquisition systems. Concurrently, FLAC3D 6.0 numerical simulations were employed to analyze the evolution of internal stress fields and the spatial distribution of plastic zones in saturated sandstone containing pore defects. Experimental results indicate that under non-uniform loading, the stress–strain curves of saturated sandstone with pore-type defects typically exhibit four distinct deformation stages. The extent of crack initiation, propagation, and coalescence is strongly correlated with the magnitude and heterogeneity of localized stress concentrations. AE parameters, including ringing counts and peak frequencies, reveal pronounced spatial partitioning. The internal stress field exhibits an overall banded pattern, with localized variations induced by stress anisotropy. Numerical simulation results further show that shear failure zones tend to cluster regionally, while tensile failure zones are more evenly distributed. Additionally, the stress field configuration at the specimen crown significantly influences the dispersion characteristics of the stress–strain response. These findings offer valuable theoretical insights and practical guidance for surrounding rock control, early warning systems, and reinforcement strategies in water-infiltrated mine roadways subjected to non-uniform loading conditions.