Abstract: The deep in-situ occurrence environment leads to obvious viscoelastic-plastic deformation behavior of deep rocks. In view of the viscoelastic-plastic deformation behavior of deep rocks in the mine, the constitutive equation of Abel visco-pot is derived based on the consensus derivative, and a three-dimensional constitutive model of viscoelastic-plastic components is established. The residual form of the local differential algebraic equation and its implicit method are used to solve the problem, and the numerical simulation is carried out by relying on MFront, MTest and FEniCS software. The analytical solution and experimental data are used to verify the reliability of the constitutive model, and it has a good convergence effect. Simulation results show that the constitutive model can well describe the triaxial compression and expansion effect of rock under the condition of confining pressure. By controlling the relevant parameters of the constitutive model, the viscoelastic deformation behavior during creep can be effectively analyzed.