Preexisting fracture network (PFN) is ubiquitous in natural rock mass and is one of the key factors affecting the mechanical properties of rock mass. This research employed a combined approach of theoretical analysis and finite element numerical simulation to investigate the mechanical behavior of PFN-containing rock mass. Based on the sliding crack model and fracture mechanics theory, the initiation mechanism of closed and open cracks driven by compressive and shear stress was theoretically analyzed. Rock mass models with different inclination angles (β) of PFN were built in Fast Lagrangian Analysis of Continua in 3 Dimensions. The complete stress–strain curves of rock mass and the secondary crack propagation were obtained. Results indicated that the increase of β reduced the strength and elastic modulus of rock mass. When β is less than 45°, this influence is more prominent. Based on the load-bearing capacity of PFN-containing rock mass, the stress–strain curves can be classified into two types: type I curve has multipeak and is strain-softening in the postpeak stage when β is 15°, 30°, 45°, 60°, and 90°; type II curve is unimodal and with brittle failure when β is 60° and 75°. Secondary crack propagation can be classified into four main types: type I is the coplanar crack, type II is the wing crack around PFN, type III is the near-field wing crack, and type IV is the macroscopic shear crack.