Most icing studies use 2D or quasi-3D models, ignoring the effects of blade span and tip icing. In this paper, a new model for predicting 3D vertical axis wind turbine blades icing by supplementing the blade span is established, the icing characteristics of blades under multiple working conditions are simulated step by step from the time scale, and the unsteady icing laws are compared and analyzed. Ice properties and surface temperature distribution determine the corresponding deicing power. The analysis results show that the unfrozen water film flows in the spanwise direction of the blade, resulting in a smaller icing limit on the pressure surface than the two-dimensional model, which indicates that the chordwise anti-icing structure can be set in a small range. At the same time, the annular ridge ice formed at the end of the blade slows down the reduction rate of the lift-drag ratio of the blade, effectively reducing the influence of the leading edge ice shape on the aerodynamic performance of the blade. Finally, the critical ice melting power density required to reach equilibrium temperature at the leading edge of the blade is reduced by about 11% compared to the trailing edge due to the increased thermal resistance of the gas-to-blade matrix interface by the ice layer.