The aeroelasticity of the wind turbine blade has been emphasized by the related fields as the size of blade increased dramatically. The eigenvalue approach and the time domain method are applied to analyze the aeroelastic responses of wind turbine blade to determine the flutter region respectively. In order to clarify the difference of the flutter analysis for different blade, two different airfoils are used. The flutter region will be obtained directly by judging the sign of the real part of the eigenvalue of the blade system using the eigenvalue approach. Then the time domain analysis of flutter of wind turbine blade will be carried out through the use of the four-order Runge-Kutta numerical method, so the flutter region will be acquired in another way. The time domain analysis can give the changing tread of the aeroelastic responses in great detail than that of the eigenvalue method. For the two different airfoils, the flutter region given by the eigenvalue approach coincides with that of the time domain analysis method accurately. There are two critical tip speed ratios for the two airfoils, the lower tip speed ratio and the higher tip speed ratio. The flap displacement of these two different airfoils will change from convergence to divergence, and change from divergence to convergence. But the extent of flutter differs with the different blade airfoil. The flutter of airfoil NACA63-418 diverges much more dramatically than that of the airfoil FX77-W-153. So the latter is better for the wind turbine blade. The eigenvalue approach combined with the time domain method can be applied to choose the blade airfoil and to determine the flutter region in order to avoid the flutter of wind turbine blade
