A better understanding of the mechanical behavior of the multi-layered structure under external loading is the most important item for the structural design and the risk assessment. The objective of this study are to propose and develop an analytical solution for the mechanical behaviors of multi-layered structure generated by axisymmetric loading, and to investigate the impact of anisotropic layers and interlayer conditions on the multi-layered structure. To reach these objectives, first, according to the governing equations, the analytical solution for a single layer was formulated by adopting the spatial Hankel transform. Then the global matrix technique is applied to achieve the analytical solution of multi-layered structure in Hankel domain. The sliding and bonded interlayer conditions were considered in this process. Finally, the numerical inversion of integral transform was used to solve the components of displacement and stress in real domain. Gauss-Legendre quadrature is a key scheme in the numerical inversion process. Moreover, following by the verification of the proposed analytical solution, one typical three-layered flexible pavement was applied as the computing carrier of numerical analysis for the multi-layered structure. The results have shown that the anisotropic layers and the interlayer conditions significantly affect the mechanical behaviors of the proposed structure.