Abstract: Based on the theory of unsaturated porous media, a computational model of two-dimensional layered unsaturated soil foundation is established and the dynamic response of the layered unsaturated soil under moving load is studied. The dynamic response solutions of two-dimensional layered unsaturated soil foundation in frequency domain are obtained by Fourier transform and Helmholtz vector decomposition principle using the transmission-reflection matrix (TRM) method. Combined with boundary and interlayer continuity conditions derived, the numerical solutions of displacement, stress and pore pressure in layered unsaturated soil are obtained by using Fourier inverse transform. The effects of load movement velocity, soil shear modulus, saturation and load frequency on the dynamic response of three types of layered foundations, namely, homogeneous soil, soft interlayer and hard interlayer foundation, are discussed through numerical calculations. The results of this study show that when the load movement velocity is close to the Rayleigh wave velocity in the unsaturated ground foundation, the soil resonates and the vertical displacement of the surface increases rapidly to the peak. The order of soil layers in the layered unsaturated foundation has a significant influence on surface displacement. Compared with the homogeneous foundation, the vertical displacement of the surface of the hard interlayer foundation is large, while the vertical displacement of the surface of the soft interlayer foundation is small. The saturation has a significant influence on the dynamic response of the layered unsaturated foundation. With the increase of load frequency, the displacements of three types of layered foundations demonstrate obvious fluctuation along the horizontal distance direction.