Investigating the microstructural damage and mechanical properties of coal under deep mine hygrothermal conditions is essential for ensuring the safe and efficient extraction of coal resources. In this study, X-ray diffraction (XRD), scanning electron microscopy (SEM), energy-dispersive spectroscopy (EDS), and nanoindentation techniques were employed to examine the surface morphology and microscale mechanical properties of coal samples exposed to four environmental conditions, initial, humidified, heated, and coupled hygrothermal, under a peak indentation load of 70 mN. The results indicate that humidification led to the formation of dissolution pores and localized surface softening, resulting in a 15.9% increase in the peak indentation depth and reductions in the hardness and elastic modulus by 29.53% and 17.14%, respectively. Heating caused localized disintegration and the collapse of the coal surface, accompanied by surface hardening, with a slight 0.4% decrease in the peak indentation depth and increases in hardness and the elastic modulus by 1.32% and 1.56%, respectively. Under the coupled hygrothermal condition, numerous fine dissolution pores and microcracks developed on the coal surface, and the mechanical properties exhibited intermediate values between those observed in the humidified and heated states. Notably, the elevated temperature suppressed the moisture penetration into the coal matrix to some extent in the hygrothermal environment. A positive correlation was found between the hardness and elastic modulus, independent of the coal sample condition. The mineralogical composition significantly influenced the microscale mechanical behavior, with hard quartz minerals corresponding to lower peak indentation depths and a higher hardness, whereas soft kaolinite showed the opposite trend.

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.