A piecewise linear Drucker-Prager strength criterion and an isotropic continuum damage model with the damage scalar depending on an equivalent tensile strain are suggested to model rock mass behavior under blast loading. A rate-dependent constitutive relation is employed to model the energy dissipation caused by two sources, namely irreversible degradation of damage and permanent deformation caused by plasticity. The suggested model is incorporated with a commercially available software AUTODYN through its user's subroutine function. Coupling of Euler and Lagrange processors are used to include all the materials under consideration such as explosive, air and rock mass, in the calculation. Using AUTODYN and the suggested model, shock wave propagation in rock mass induced by an underground explosion is simulated. Numerical results obtained agree favorably well with those obtained from an independently conducted field test. It demonstrates that the suggested model can be used to predict the damage area, plastic zone and ground motions generated by underground explosions. (C) 1998 Elsevier Science Ltd. All rights reserved.