Based on the two-node Euler-Bernoulli beam, the tower system is discretized by finite element method, and the cubic Hermite polynomial is taken as the shape function of the beam element, and the structural characteristic matrix of the tower system is calculated, and the wind turbine-nacelle-tower multi-degree of freedom is established Finite element numerical model. The aerodynamic load calculation formula for any nacelle attitude angle is deduced. The influence of the vibration feedback of the flexible tower on the aerodynamic load of the wind turbine is studied. The results show that when the rigidity of the tower is large, the impact of tower vibration feedback on the aeroelastic load of the wind turbine is small. For a tower system with greater flexibility, the time-varying feedback of wind-induced vibration will cause greater aeroelastic load changes, especially the overturning moment of the tower top, which will cause a greater impact on the dynamic behavior of the tower in the downwind and crosswind directions. As the flexibility of the tower system increases, the interaction between tower vibration and aerodynamic load is gradually increasing. Taking the impact of the flexible tower on the aeroelastic load of the wind turbine into account, on the one hand, helps to predict the wind more accurately. The aerodynamic load of the wind turbine improves the efficiency of wind energy utilization. On the other hand, it can more accurately analyze the dynamic behavior of the flexible structure of the wind turbine, which is extremely beneficial to the structural optimization design of the wind turbine.