The stability of slope under complex geological conditions is one of the essential issues in geotechnical engineering. The mechanical damage caused by external force load and the hydraulic property of rock have both significant influences on slope stability. However, the two factors and the coupling process between them are not considered in the traditional slope safety calculation method. Firstly, in this paper, an elastoplastic-damage-seepage coupling model is established based on the modified Mohr-Coulomb (M-C) criterion. This model takes into account the weakening effect of damage on the strength of rock mass. Further, according to the Biot's theory and the evolution equation of permeability coefficient, a complete hydraulic-mechanics (HM) coupling model is established. Secondly, the fully implicit backward Euler algorithm and "strict-cornered" algorithm are used to integrate the stress precisely. The complete coupling analysis of HM is implemented by the step-by-step iterative method. By combining the centrifugal loading method (CLM) with the FEM program, the factor of slope safety (FOS) under the coupling influence of damage-seepage can be solved. The calculation results showed that, compared with the traditional rounded M-C model, the results obtained by the "strict-cornered" algorithm method in this paper are closer to the theoretical solution. The solution process is robust with a global second-order convergence rate. After considering the effects of damage and seepage, the FOS is obviously reduced. Finally, the model is applied to the stability evaluation of an actual engineering. The variation rules of pore water pressure, displacement field, damage zone and FOS are calculated under different head heights. The model built in this paper considers the coupling characteristics of rock seepage-stress-damage well. The proposed numerical algorithm ensures the accuracy and rationality of the calculation results which provides a theoretical foundation for engineering application.

复杂地质条件下边坡的稳定性是岩土工程中的重要问题之一。由外力载荷引起的机械损伤和岩石的水力特性均对边坡稳定性产生重大影响。但是，传统的边坡安全计算方法没有考虑这两个因素以及它们之间的耦合过程。首先，基于改进的莫尔-库仑（MC）准则，建立了弹塑性-损伤-渗流耦合模型。该模型考虑了破坏对岩体强度的减弱作用。此外，根据毕奥特的理论和渗透系数的演化方程，建立了完整的液力耦合模型。其次，完全隐式的向后欧拉算法和“ 在不同的头部高度下计算损伤区域和FOS。本文建立的模型考虑了岩石的渗流-应力-损伤耦合特征。所提出的数值算法保证了计算结果的准确性和合理性，为工程应用提供了理论基础。