Abstract: The destruction of the rock mass is the shear damage that mainly occurs along the failure surface in nature, and it is mainly characterized by shear rheological damage. Based on the analysis of experimental data of weak intercalation under different shear stress levels, the variable D which can characterize the parameter damage in the rheological process was introduced. A viscoelastic plastic nonlinear rheological model based on D, which can reflect the accelerated rheological characteristics of weak intercalation, was proposed. And a new shear rheological damage model of weak intercalation was constructed by connecting it with the Burgers model, which can fully reflect the three rheological stages. Then the rheological damage model was redeveloped with C++ language based on FLAC^3D, and compared with the indoor rheological test. The results show that: (1) The instantaneous strain of the increment value of strain under each shear stress is quite close, and both the stable strain value and the accelerated strain value are very close after the same time of stress loading. (2) When the applied shear load does not reach the yield stress, the deceleration period of the numerical simulation test is shorter than that of the indoor test, that is, it reaches the stable rheological stage faster. When the applied shear load exceeds the yield stress threshold, the displacement change rate of the accelerated rheological stage is larger than that of the indoor test, and the final displacement is slightly larger than the indoor test displacement.