Abstract: In the paper, the singular solutions of axisymmetric elastodynamic problem for the multilayered halfspace due to unit dynamic load acting within layered media are derived in time domain with the aid of Laplace-Hankel mixed transforms and transfer matrix method. A numerical computation of these solurions is also proposed using a fast Hankel transform algorithm. Two examples show that the solutions can be readily evaluated and the numerical results have considerably high precision. The presented solutions can...