Abstract: According to the analytical element of a single layer and the continuity conditions between each layer, the vertical displacements of transversely isotropic multi-layered soils under arbitrary vertical uniform loads are obtained. Based on the displacement coordination and smooth contact, the analytical solutions to the interaction problem between the flexible rectangular footing and transversely isotropic multi-layered soils are derived. A Fortran module is developed to analyze the effects of the rigidity of the rectangular footing, transversely isotropic and stratified characters. The results show that the rigidity ratio of plate to soil γ, transversely isotropic parameter m and stratified parameter α exert a great influence on the vertical displacements of soils, while the influence of parameters n and β is limited.