A computational method based on nonlinear wake model was established for horizontal axis wind turbines aerodynamic performance prediction. This method makes use of finite difference method to solve the integral differential equation of the model, the induced velocity of wake vortex can be calculated from equations and compared with the induced velocity of wake vortex in linear model. The comparison between the calculated results of wind turbine under axis flow condition, including tip vortex geometry and aerodynamic performance, and available experimental data shows that this method is suitable for wind turbine aerodynamic performance analysis. Finally, a series of numerical calculations were made to investigate the change of wake geometry and aerodynamic performance of the wind turbine when yawing and pitch angle increasing, which provide foundations for aerodynamic optimization design of horizontal axis wind turbines.