The pavements are subjected to the moving vehicle, pore pressure and temperature loadings, the analytical solutions of thermo-hydro-mechanical response for the pavements are paid more attentions by the road engineers while the pavements are maintained. The analytical solutions of the multi-layered pavement model are derived by the novel DVP (duel variable and position) method in this paper. Starting with the basic equations of the pavement system, obeying Biot's dynamic poroelastic theory and generalized thermoelastic theory, the original problem is divided into two sub-systems and solved twice to avoid tangling two different loading frequencies. With the aid of moving coordinate and vector functions, the decomposition expressions of stress, displacement, temperature and pore pressure variations are deduced and propagated by using DVP method, and the transformed solutions are numerically inverted. In the numerical examples, the obtained solutions are firstly verified. Then the difference among the results of several pavement models are compared and analyzed. Finally, the effects of the surface temperature, moving load speed and interface conditions on the displacement, stress, and strain distributions are discussed by the solutions of thermo-hydro-mechanical pavement models in the numerical results.