Mesoscopic Damage Evolution of Water-Bearing Mudstone Under Low-Strain-Rate Cyclic Dynamic Loading: A Particle-Flow Simulation Study

Groundwater-related weakening and repeated low-strain-rate disturbances can jointly affect the long-term stability of soft surrounding rock. In this study, a PFC3D model was developed by combining a cumulative residual-strain damage variable with a Weibull statistical damage variable. Water effects were represented indirectly through water-content-dependent mesoscopic parameters and contact-strength degradation, thereby linking macroscopic irreversible deformation to progressive mesoscopic bond degradation. The cyclic responses of mudstone specimens with water contents of 0%, 3%, 5%, and 7.04% were simulated under a roadway-like three-directional, five-face boundary condition, consisting of axial loading, lateral pressure, one laterally constrained side, and one free face. The model reproduced the main stress–strain trends and failure characteristics. From the dry to saturated state, the experimental and simulated peak strengths decreased by 68.51% and 67.44%, respectively. Increasing water content also promoted earlier crack initiation, weakened strong force-chain continuity, and shifted failure from localized shear instability to distributed tensile–shear damage. Energy dissipation became increasingly important as water content increased. These findings clarify the mesoscopic damage mechanism of water-bearing mudstone under coupled water-induced softening and low-strain-rate cyclic loading, providing a reference for stability assessment and support design in underground soft-rock engineering.