Fracture is one of the most common failure modes of materials and components and greatly restricts engineering design. Understanding of the crack propagation and evolution of rock and other engineering materials is of great significance to engineering construction. For the current numerical methods there are more or less limitations when analyzing the evolution of cracks

such as the grid dependence of the crack path

the difficulty to deal with crack bifurcation and merging by the classic fracture criterion. In recent years

the phase field method (PFM) has been widely used in simulating crack growth. A phase field numerical manifold method (PFNMM) makes use of the advantages of the phase field method in simulating crack propagation and those of the numerical manifold method (NMM)

and is proposed for crack growth in rock. The implementation details of the proposed numerical model are presented. Several benchmark examples

including notched semi-circular bend test and Brazilian disc test

are adopted to validate the proposed numerical approach. After that

the multi-crack propagation process with different rock bridge inclination angles under uniaxial compression is simulated

which is in good agreement with the results derived from laboratory and PFC. And the results indicate that the PFNMM has broad application prospects in simulating crack growth of rock.