In view of the strong constraint zones of the concrete structure on the pile foundation, there are some differences between the calculation results of the isotropic equivalent pile foundation by the volume replacement ratio method and the actual engineering. In this paper, referring to the relevant algorithm of rock mass with anchor, the anchor and rock mass are, respectively, compared to pile and surrounding soil foundation. Eshelby equivalent inclusion theory is introduced into the equivalent mechanical model of soil foundation with pile, and a new equivalent pile foundation algorithm considering anisotropic elastic constant is compiled by Fortran. Three kinds of calculation methods are used to calculate the stress field of the concrete structure of the large pump station on the pile foundation during the construction period, and the stress in the strong constraint zones of the concrete structure are mainly analyzed. It is found that the calculation accuracy of Algorithm 3 is the highest, and the calculation results of Algorithm 2 can be modified by the coefficients to achieve the calculation accuracy of Algorithm 3 and the calculation efficiency is actually improved. Finally, the accuracy of the proposed method is verified by the engineering measured data.