AbstractThe conventional transfer matrix method, adopted formerly in a semianalytical solution for layered elastic or viscoelastic asphalt pavement structures, has inherent deficiencies, such as ill-condition matrix, numerical overflow, and error accumulation, due to exponential items. This phenomenon is more evident in multilayered dynamic analysis with imperfect interfaces. Moreover, several studies revealed that pavement materials exhibit transverse isotropy in service. Consequently, a novel semianalytical solution methodology, wave propagation approach, was proposed herein to calculate the dynamic responses of asphalt pavement under an impact load considering the transversely isotropic material model and the imperfect interfaces. First, the transfer matrix was established based on matrix theory and wave propagation approach, while the relation between the state vector and wave vector in the transformed domain was constructed simultaneously. Then, combined with the boundary conditions and interface contact conditions, the solution of the wave vector in the transformed domain was derived. Finally, based on Laplace–Hankel inverse transform, the state vector in the time domain was obtained, followed by numerical computation with programming. The accuracy and efficiency of the proposed semianalytical solution, together with the influence regularities of several variables, were discussed. Results showed that due to the absence of positive exponential functions and a large-dimensional matrix, accuracy and efficiency requirements were satisfied during calculation. Moreover, the variation induced by the transversely isotropic properties and interface conditions, presented in the dynamic responses, reiterated that these factors should be considered during the design and analysis of asphalt pavement structures.