Aircraft icing can significantly degrade aerodynamic performance and compromise flight safety, highlighting the need for efficient prediction of icing characteristics and critical influencing factors. In this study, a structured set of icing simulation parameters was defined by identifying key variables and constraining their ranges in accordance with airworthiness requirements of intermittent maximum icing conditions. According to the simulation results, it concludes that heavier ice accretion leads to increased surface roughness and consequently more severe aerodynamic degradation under the intermittent maximum icing conditions. Within this admissible domain, representative samples were generated and applied to FENSAP-ICE to simulate ice accretion on the reference aircraft airfoil, thereby constructing a comprehensive icing dataset with the iced mass as the output response. Based on this dataset, a surrogate model employing the Random Forest algorithm was developed, enabling rapid and reliable prediction of icing outcomes. Furthermore, sensitivity analysis was conducted to identify the dominant parameters governing the icing process, offering valuable insights into the core factors that critically influence aircraft icing behavior.