Abstract: The analytical layer element solution for axisymetrically transversely isotropic multilayered soils is adopted to solve the interaction between the rigid circular plate and the foundation. Then based on the displacement continuity conditions between the soil surface and the rigid circular disc, the integral equation of an axisymmetrically loaded rigid circular disc on transversely isotropic multilayered soils is developed. The integral equation is solved numerically, and the contact stresses between the rigid circular disc and the subsoil are obtained when the rigid circular disc is given a value of displacement. This solution is implemented into a computer code for numerical analysis, and the influences of parameters of transversely isotropic foundation, thickness and stratification of the soils on the contact stresses are analyzed.