Abstract:Starting from the governing equations of thermoelasticity for the plane strain problem, the analytical solution of the problem is obtained by using the Fourier transform and Laplace transform. And then, the exact stiffness matrix for the multi layered plane strain problem in stable temperature field is established, which is called analytical layer element. The global stiffness matrix is assembled by considering the boundary conditions and the continuity between adjacent layers, and the solution in the transform domain are obtained by solving the equations of the global stiffness matrix. The actual solutions in the physical domain are acquired by inverting the Laplace Fourier transform. Numerical calculation and analysis are carried out by corresponding computer program whose results agree well with those by finite element analysis software, and numerical calculation shows the layered characteristics has a significant effect on the thermal stress and vertical displacement in layered pavement system.