Non-parallel overlapping flaws widely exist in engineering rock mass. Understanding their crack initiation and energy evolution characteristics is of great significance to ensure the stability of rock engineering. Based on the existing experiments, the influence of flaw inclination angles (β) on the crack initiation and energy evolution characteristics of rock samples with non-parallel overlapping flaws was studied by numerical simulation. The results show that (1) the uniaxial compressive strength, elastic modulus and crack initiation stress increase with the increase of flaw angle. (2) The boundary energy, strain energy and dissipated energy under peak stress increase with the increase of flaw angle; the dissipated energy increases the most. (3) With the increase of flaw angle, the tension stress zone is transferred to the flaw tip, and the zone is reduced gradually; the maximum tensile stress and the tension stress concentration decrease. (4) In the crack initiation stage, the influence of a lower flaw inclination angle (β ≤ 60°) on the lateral displacement field of the sample is higher than that of a high flaw inclination angle (β = 75°).