Abstract：Rock engineering hazards are closely related to unstable failure of rocks. Numerical analysis of rock failure process of large-scale rock engineering needs effective，accurate and powerful computation，while traditional serial computation becomes incapable to solve these large-scale rock failure problems；and it is necessary to employ large-scale parallel computation technology. Rock failure process analysis(RFPA) code is one of the important numerical tools that can be used to investigate rock failure process. Based on the serial code of RFPA3D，a parallel computation model of rock failure process analysis is proplsed. We complete parallel stress analysis module of RFPA3D using finite element method on a cluster and integrate it with pre-processing and post-processing of RFPA3D installed on windows pc. The parallel program is performed using a distributed memory sparse linear iterative solver with preconditioning based on MPI(message passing interface). The linear solver is fully parallel，working with sub-domains determined by domain decomposition parallel algorithms. Numerical tests demonstrate that this parallel program is robust and has very high performance as expected. The parallel RFPA3D system can be used to quickly fulfill rock failure process analysis of model with more than ten millions of degrees of freedom on a 32-node cluster with 64-CPU. With the application of the RFPA3D-Parallel system to simulate the fracture spacing phenomenon in the medium of earth crust，failure process and fracture evolution in heterogeneous medium like rock are revealed；and it indicates the broad potential application of the system.
