“…Based on the discrete element method (DEM) proposed by Cundall [20], the universal distinct element code (UDEC) has been widely used to calculate the propagation problems of stress waves in a jointed rock mass [21][22][23]. Furthermore, other numerical methods and software have been adopted to solve the problems involving the stress wave propagation in a rock mass, e.g., the particle manifold method (PMM) [24,25], the numerical manifold method (NMM) [26,27], the particle flow code (PFC) [28,29], and the three-dimensional element code (3DEC) [30]. However, the above theoretical, experimental, and numerical methods have mainly focused on the effect of the parameters of the joints, e.g., joint stiffness, joint spacing, joint number, and the parameters of the stress wave, e.g., waveform, amplitude, frequency, and the incident angle of stress wave on the stress wave propagation pattern, and have proposed that the attenuation of the stress wave only occurs at the joints, while it has been assumed that the intact rock is elastic.…”