Disastrous rock slope failures have been posing a hazard to people’s lives and causing enormous economic losses worldwide. Numerical simulation of rock slope failure can lead to improve the degree of understand of such phenomenon so as to predict and avoid the occurrence of these disastrous events. In order to simulate the global behaviors of rock slope failure under the high seepage pressure and the local behaviors of the occurrence of hydraulic fracture in the pre-existing rock joints effectively, a powerful finite element tools F-RFPA2D, is adopted. The simulation takes into account of the growth of existing fractures and the initiation of new fractures under various of hydraulic pressure in different heterogeneities medium. The behavior of fluid flow and damage evolution, and their coupling action are studied in small specimens that are subjected to both hydraulic and biaxial compressive loadings. The influence of the ratio (the initial horizontal stress to the initial vertical stress) and the distance between the two existing cracks on the fracture propagation behaviors are investigated. Moreover, based on the fundamental study of hydraulic fracture, the progressive failure of rock slope under the influence of the increase in hydraulic pressure was also studied in the paper.