Researchers often get the characteristics of EDZ with the experiment and the engineering detection. But the method of full size experiment is very expensive and time-consuming. The engineering detection is after the event in the engineering, which is not good for the control of EDZ. Numerical method, derived from sound mechanical principles and validated against experimental data, indicates a promising approach to reveal the EDZ formation process. For this case, various damage models have been suggested to study the damage of rock mass under blasting. For instance, Grady & Kipp (1980), Taylor (1986) and his coworkers, Yang et al. (1996) and Liu & Katsabanis (1997) developed several widely used continuum damage models to describe the dynamic damage process of rock mass in tension. The damage variable is introduced to describe the changing state of a material. However, the damage variable is established based on the tensile strain and compressive damage is neglected in these damage models. Thus far, application of an appropriate damage model to simulate the different damage pattern is rarely found.