In this paper, a new direct calculation method of frost ice shape on the blade surface of horizontal axis wind turbine is proposed. Using linear interpolation algorithm, the airfoil ice shape obtained by LEWICE 2D software or ice wind tunnel experiment was fitted with equidistant step length in the first and fourth quadrants and equidistant step length in the second and third quadrants. The key point coordinates of ice shapes on cross-sections along the span-wise were mapped into lagging and flapping surfaces through the mathematical dimension reduction, respectively. The polynomial fitting was used to deal with ice projection points of multiple sections in lagging and flapping surfaces, and then the blade’s frost ice shape was obtained. By calculating the sum of squared residuals of the polar diameter at the same polar angle, the errors between experimental and airfoil frost shape fitting methods, experimental and FENSAP, and blade frost shape formation methods and FENSAP were analyzed. The results show that the new method is in good agreement with the ice shape of FENSAP simulation results and experimental results. The residual sum of squares is small. This method makes the analysis of frost ice morphology of wind turbine blades do not need to consider interdisciplinary. The calculation process is simple and reliable.
Read full abstract