Under the influence of engineering disturbances, the loading rate of surrounding rock is in a state of continuous adjustment. This study conducts experimental investigations on the mechanical response characteristics under different strain rates (10⁻⁵ s⁻¹, 10⁻⁴ s⁻¹, and 10⁻³ s⁻¹). During the uniaxial loading process of coal–rock composite specimens, multi-parameter monitoring was implemented, and a systematic study was carried out on the ring-down count induced by microcracks, the energy values of acoustic emission (AE) events, the stage-dependent strain characteristics on the specimen surface, and the surface temperature variation characteristics. Additionally, the stress–strain curve characteristics under different strain rates were comparatively analyzed in stages. The loading process of the coal–rock composite specimens was reproduced using the Particle Flow Code (PFC3D 6.0) simulation software. The simulation results indicate that the stress–strain results obtained from the simulation are in good agreement with the laboratory test results; based on these simulation results, the energy accumulation and dissipation characteristics of the coal–rock composite specimens under the influence of strain rate were revealed. Furthermore, a microscopic damage model considering strain rate was constructed based on the Weibull probability statistics theory. The results show that strain rate has a significant impact on the strength, elastic modulus, and failure mode of the coal–rock composite specimens. At low strain rates, the specimens exhibit obvious progressive failure characteristics and strain localization phenomena, while at higher strain rates, they show brittle sudden failure characteristics. Meanwhile, the thermal imaging results reveal that at high strain rates, the overall temperature rise in the composite specimens is rapid, whereas at low strain rates, the overall temperature rise is slow—but the temperature rise in the coal portion is faster than that in the rock portion. The peak temperature at high strain rates is approximately 2 °C higher than that at low strain rates. The PFC simulation results demonstrate that the larger the strain rate, the faster the growth rate of plastic energy in the post-peak stage and the faster the release rate of elastic energy.

Crack initiation, propagation, and slippage serve as the key mesoscopic mechanisms contributing to the deterioration of deep tunnel surrounding rocks. In this study, a secondary anisotropy of deep tunnels surrounding rocks was proposed: The axial-displacement constraint of deep tunnels forces cracks in the surrounding rock to initiate, propagate, and slip in planes parallel to the tunnel axial direction. These cracks have no significant effect on the axial strength of the surrounding rock but significantly reduce the tangential strength, resulting in the secondary anisotropy. First, the secondary anisotropy was verified by a hybrid stress–strain controlled true triaxial test of sandstone specimens, a CT 3D (computed tomography three-dimensional) reconstruction of a fractured sandstone specimen, a numerical simulation of heterogeneous rock specimens, and field borehole TV (television) images. Subsequently, a novel SSA (strain-softening and secondary anisotropy) constitutive model was developed to characterise the secondary anisotropy of the surrounding rock and developed using C++ into a numerical form that can be called by FLAC^3D (Fast Lagrangian Analysis of Continua in 3 Dimensions). Finally, effects of secondary anisotropy on a deep tunnel surrounding rock were analysed by comparing the results calculated by the SSA model and a uniform strain-softening model. The results show that considering the secondary anisotropy, the extent of strain-softening of the surrounding rock was mitigated, particularly the axial strain-softening. Moreover, it reduced the surface displacement, plastic zone, and dissipated plastic strain energy of the surrounding rock. The proposed SSA model can precisely characterise the objectively existent secondary anisotropy, enhancing the accuracy of numerical simulations for tunnels, particularly for deep tunnels.