Abstract: This Program adopts nonlinear variable-modular, uses uniaxial strain load-unload curve of rocks and dynamic tensile strength. In calculation, adopted load-unload criterion of stress wave propagation was to decide load-unload state of element and to adjusts and transfers internal stress for tensile failure element as well. The program uses incremental method to obtain displacements, velocity acceleration, stresses and strains for any instant. The results show that stresses and strains of calculation are well coincident with those from measurements.