Abstract: Based on the theory of porous media containing porous solids, established is a two-dimensional calculation model for layered saturated frozen soil foundation, and studied is the dynamic response of layered saturated frozen soil foundation under harmonic load. Through Helmholtz theorem and Fourier integral transform, the stiffness matrix of layered saturated frozen soil foundation is derived by using one transfer matrix method and combining boundary and interlayer continuity conditions, and the dynamic response solution of layered saturated frozen soil foundation in frequency domain is obtained. The numerical solutions of displacement and stress of each phase in layered saturated frozen soil are obtained by fast Fourier transform. After comparing with the existing literatures, the effects of surface soil shear modulus, of surface soil temperature, of surface soil porosity and of load frequency on dynamic response in two typical layered foundations of upper soft and lower hard and upper hard and lower soft are studied by numerical examples. The computational results show that: the stress amplitude shared by the vertical displacement and pore water decreases with the increase of shear modulus of surface soil and, increases with the increase of surface soil temperature, of surface soil porosity, and of loading frequency. The arrangement order of soft and hard soil layers in the layered saturated frozen soil foundation has a significant effect on the vertical displacement and the stress shared by pore water.