In this paper, the interaction between level ice and wind turbine tower is simulated by the explicit nonlinear code LS-DYNA. The isotropic elasto-plastic material model is used for the level ice, in which ice crushing failure is considered. The effects of ice mesh size and ice failure strain on ice forces are investigated. The results indicate that these parameters have a significant effect on the ice crushing loads. To validate and benchmark the numerical simulations, experimental data on level ice-wind turbine tower interactions are used. First, the failure strains of the ice models with different mesh sizes are calibrated using the measured maximum ice force from one test. Next, the calibrated ice models with different mesh sizes are applied for other tests, and the simulated results are compared to corresponding model test data. The effects of the impact speed and the size of wind turbine tower on the comparison between the simulated and measured results are studied. The comparison results show that the numerical simulations can capture the trend of the ice loads with the impact speed and the size of wind turbine tower. When a mesh size of ice model is 1.5 times the ice thickness, the simulations can give more accurate estimations in terms of maximum ice loads for all tests, i.e., good agreement between the simulated and measured results is achieved.