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.