Blasting has been widely used in coal mining as a method of reducing the integrity of rock mass. Its low controllability often leads to the unsatisfactory fragmentation effect of rock mass. The empty borehole effect has great significance in avoiding the excessive breakage of rock mass and guiding the directional propagation of blasting-induced fractures. Nevertheless, the rock mass damage type evolution induced by the empty borehole has been rarely studied. A tension-compression constitutive model of rock mass damage is established in this paper. The model is incorporated into the numerical modeling code LS-DYNA as a user-defined material model. Then, LS-DYNA is used to investigate the evolution mechanism of rock mass damage under the empty borehole effect. The damage types of rock mass caused by the empty borehole effect are also studied. The Fortran language is utilized to monitor the number variation of the tensile damaged elements and the compressive damaged elements in the rock mass around the empty borehole. The results indicate that existence of the empty borehole significantly enhances the tensile stress and the stress concentration factor in the rock mass nearby the empty borehole. Meanwhile, the rock mass nearby the empty borehole mainly damages in tension. Both the number of the tensile damaged elements and the tensile stresses in the elements increases as the empty borehole diameter increases. The number of the compressive damaged elements decreases with increasing empty borehole diameter.Energies 2020, 13, 756 2 of 21 mass. Yue et al. [4] found that the empty borehole is able to change the stress distribution in the rock mass and induce the initiation and directional propagation of the blasting-induced fractures.The existence of the empty borehole in the rock can change the stress distribution and the damage zone in the rock mass [5]. The empty borehole concept has been successfully applied to rock fracturing, tunnel excavation, and coal seam gas extraction [6,7]. The dynamic evolution of the blasting-induced stress and the directional propagation of the blasting-induced fractures have been widely studied by many researchers based on theoretical analysis and laboratory experiments. Cho et al. [8] investigated the influence of the empty borehole on fracture propagation, using polymethyl methacrylate (PMMA) and revealed the law of fracture propagation from the aspect of stress concentration. Based on the evolution characteristics of the fractures between blasting boreholes, He and Yang [9] studied the law of fracture propagation between adjacent boreholes, using a high-speed camera and digital image processing. Ma and An [10] simulated the stress evolution and fracture propagation around an empty borehole in the blasting process by numerical modeling and provided a theoretical basis for the formation of blasting-induced fractures. Chen et al. [11] used LS-DYNA (a numerical modeling code) to simulate the effective stress distribution in the rock under the superimposition effect of stress waves and refl...