Numerical Simulation Study of the Permeability Evolution of Sandstone during the Loading and Failure Process Based on FDEM

Abstract Rock failure under hydromechanical coupling is a complex process that has attracted considerable attention in deep underground engineering. In this study, a hydraulic coupling analysis method for rock loading and failure is developed based on the finite–discrete-element method (FDEM) combined with the grain-based model, which integrates the dual-medium seepage–stress coupling theory to simultaneously capture pore and fracture seepage. Numerical simulations of sandstone under varying confining pressures were performed and validated against laboratory experiments, and the effects of confining pressure and weak joint content on the permeability evolution were systematically investigated. The results demonstrate that the permeability evolution during loading is governed by the competitive interplay between matrix seepage and fracture seepage, exhibiting a nonlinear pattern. Higher confining pressures prolong the stage of slow permeability growth, whereas an increase in the weak joint content reduces strength and stiffness, leading to an earlier occurrence of the permeability inflection point and an enhanced dominance of fracture flow in controlling macroscopic permeability. These findings provide new insights into the mechanisms of hydromechanical coupling in fractured porous rocks and offer theoretical support for predicting and controlling hydraulic hazards in deep rock engineering.