The dynamic response is significant to the elastic wave problem in soil caused by the external load. This paper proposes a solution to calculate dynamic stress responses of an arbitrary point in a transverse isotropic multilayered soil subjected to a time-harmonic load. The generalized plane-strain equation is transformed from frequency-spatial domain into frequency-wave number domain by Fourier transformation in this algorithm. Combined with the introduction of the dual vector

the state equation is solved by the precise integration method. Based on the displacement response of the soil in the frequency-wave number domain

the dynamic stress response of any point is obtained by the inverse Fourier transformation. The time-harmonic load can be applied at the surface of the soil or under ground. The accuracy of the algorithm in this paper is verified by a comparison with an existing solution. An extensive parametric analysis on the influence of anisotropy

excited frequency and damping ratio on the dynamic stress response provides reliable numerical basis for engineering practice.