Abstract: An elasto-plastic damage seepage coupling model for rock based on the Mohr-Coulomb criterion is established to analyze slope stability. The stress return method in the principal stress space is given for the Mohr-Coulomb criterion under damage effect to solve the singular point problem during stress updating. The elastoplastic-damage-seepage finite element program for rock is compiled based on the stepwise iterative method. The centrifugal loading method is combined with the program which is used to calculate the factor of safety of slopes under coupling effect, and the damage evolution rules during failure process are represented. Finally, the numerical simulation is done for an actual slope, and the results may provide reference for the evaluation of engineering safety.