Axisymmetric structures have applications in various fields such as engineering and architecture. This type of structure exhibits complex stress distribution and failure modes when subjected to ultimate loads, and the importance of the analysis method for its upper limit in the engineering field is self-evident. Its upper limit analysis is mainly used to evaluate the stability and load-bearing capacity of axisymmetric structures under load. The intention of the fully edge smooth finite interpolation method is to improve the accuracy and efficiency of analysis. It improves the interpolation function to maintain smoothness at the boundary, and uses adaptive mesh partitioning and local refinement to make the analysis more accurate and efficient. Therefore, the purpose of this article's in-depth study on the perfect edge smooth finite interpolation method is to introduce smooth functions and finite interpolation techniques, accurately simulate structures, avoid risks, and reduce losses. This article mainly applies numerical simulation and experimental comparison to compare the axisymmetric structures of thick walled cylinders, frustums, and spheres. The ultimate load of each structure is obtained by changing the distortion coefficient and radius ratio. The experimental results show that in cylindrical testing, the performance difference between the two methods is greatest when the radius ratio is 3.5. The error between the perfectly smooth finite interpolation method and the analytical value is smaller, while the direct iteration method has a larger error.