Icing accretion on wind turbine will degrade its performance, resulting in reduction of output power and even leading to accidents. For solving this problem, it is necessary to predict the icing type and shape on wind turbine blade, and evaluate the variation of aerodynamic characteristics. In this paper the icing types and shapes in presence of airfoil, selected from blade of 1.5 MW horizontal-axis wind turbine, are simulated under different ambient temperatures and icing time lengths. Based on the icing simulation results, the aerodynamic characteristics of icing airfoils are simulated, including lift and drag coefficient, lift-drag ratio, etc. The simulation results show that the glaze ice with two horns presents on airfoil under high ambient temperature such as -5°C. When ambient temperatures are low, such as -10°C and -15°C, the rime ices with streamline profiles present on the airfoil. With increase in icing time the lift forces and coefficients decrease, and the drag ones increase. According to the variations of lift-drag ratios of icing airfoil, the aerodynamic performance of airfoil deteriorates in the presence of icing. The glaze ice has great effect on aerodynamic characteristics of airfoil. The research findings lay theoretical foundation for icing wind tunnel experiment.