The scaled boundary finite element method coupled with the cohesive crack model is extended to investigate the hydraulic fracture at the dam-foundation interface. The concrete and rock bulk are modeled by the scaled boundary polygons. Cohesive interface elements model the fracture process zone between the crack faces. The cohesive tractions are modeled as side-face tractions in the scaled boundary polygons. The solution of the stress field around the crack tip is expressed semi-analytically as a power series. Accurate displacement field, stress field and stress intensity factors can be obtained without asymptotic enrichment or local mesh refinement. The proposed procedure is verified by the hydraulic fracture of a rectangular embankment on rigid foundation and applied to the modeling of hydraulic fracture on the dam-foundation interface of a benchmark dam. Different distributions of water pressure inside the crack are investigated. It is found that the water pressure inside the crack decreases the peak overflow to less than 20% of the case without water in the crack. Considering the water lag or not is significant to the response, while different distribution of pressure following the water lag region in the fracture process zone has negligible influence.