Numerical Investigation on the Creep-Induced Microdamage Evolution in Rock

Rock creep, a key factor in the long-term stability of deep geotechnical engineering, remains challenging to study due to the complexity of its microscopic damage mechanisms. Laboratory creep tests are limited by long durations and scale effects, while phenomenological models cannot fully capture the underlying processes. This study employs the parallel-bonded stress corrosion (PSC) model in PFC2D to simulate sandy mudstone’s creep behavior, systematically correlating macroscopic creep deformation with microscopic damage evolution and energy conversion. The model reproduces the four stages of the idealized creep curve and quantifies the effects of axial stress level and confining pressure on creep lifetime, rate, and failure mode. Increasing axial stress shortens creep lifetime; every 10% increase raises the creep rate by a factor of 4–14, and high stress enhances nonlinear deformation, producing stair-stepping curves due to unstable microcrack propagation. In contrast, confining pressure prolongs lifetime; at 90% uniaxial compressive strength (UCS), 15 MPa extends it from 2.78 h to ~25 years. Confinement also enhances ductility by suppressing tensile stresses and delaying damage accumulation. This study reveals the coupling mechanism of stress-corrosion-induced subcritical crack propagation and energy dissipation, clarifies the microscopic origin of stepped creep curves, and provides a micromechanical framework for long-term stability evaluation in deep geotechnical engineering.