Abstract
In order to predict blast vibration more precisely, we attempted to use the following equation described by circular frequency, that is used in the earthquake engineering field.
O(ω) =G(χ)·H(ω)·S(ω)
where, O(ω) is spectra of the wave motion of ground surface, G(ω) is the geometrical attenuation related energy spread during wave propagation, H(ω) is the transfer function and S(ω) is the source function.
We developed a computer program for a multi-reflection method based on Haskell (1960) to get H(ω) in above equation. For S(ω), we used the equation derived by Sharpe (1942). Finally, we had reasonable results for wave motion of a point of ground surface by calculating inverse FFT of O(ω) obtained by convolution of H(ω) and S(ω) about two layers model and four layers model.