This paper proposes a fatigue damage life prediction method for wind turbine blades that incorporates aeroelastic characteristics of pre-bend blade profiles. Based on the modified blade element momentum theory and combined with the dynamic stall model, an aerodynamic load calculation model for pre-bent blades was established. The quadratic Lagrange interpolation function was used to fit the blade centerline, and the blade mass, stiffness, and energy matrices were calculated based on the geometric exact beam theory. The Galerkin method was applied to discretize the aerodynamic loads, and the Newton-Raphson iteration method was used to solve the dynamic equations, thereby establishing an aeroelastic coupling model. Based on the aeroelastic coupling model, precise aerodynamic loads under different wind speeds were calculated and applied to the blade structural model. Wind speed time histories are generated using the Weibull distribution, and the blade stress spectrum is obtained through the Rain-flow counting method. The remaining fatigue life is then predicted based on the Miner rule, using the S-N curve to accumulate fatigue damage. Comparative results indicate that traditional methods for predicting wind turbine blade life tend to overestimate the fatigue lifespan, whereas the proposed method, by considering aeroelastic coupling effects, provides a more accurate estimation of blade life. The method proposed in this paper offers a precise and reliable tool for the assessment of fatigue damage and the prediction of fatigue life in wind turbine blades.