Ice and snow environment simulation test method based on multi-field coupling

Abstract

The invention relates to the technical field of ice and snow environment simulation, and discloses an ice and snow environment simulation test method based on multi-field coupling. The method comprises the following steps: collecting a dynamic multi-field data set in a target ice and snow environment, wherein the dynamic multi-field data set covers a periodic temperature change sequence, a relative humidity fluctuation sequence and surface ice layer thickness real-time monitoring data; carrying out multi-physics field coupling numerical simulation processing to generate coupling field distribution characteristics including an axial heat conduction coefficient gradient, a radial fluid velocity cumulant and a tangential ice layer stress fluctuation coefficient; calling a pre-configured ice and snow growth simulation model to execute ice layer evolution analysis to obtain an ice layer growth rate prediction value and a high-risk icing area position identifier; performing environment dynamic compensation correction on the ice layer growth rate prediction value by using the correlation characteristic of the relative humidity fluctuation sequence and the material heat capacity attribute; based on the high-risk icing area position identification, scientific guidance is provided for related protection in the ice and snow environment.

Description

Technical Field

[0001] The present invention relates to the technical field of ice and snow environment simulation, and in particular to an ice and snow environment simulation test method based on multi-field coupling. Background Art

[0002] The impact of icy and snowy environments remains a key concern in many fields, including engineering construction, transportation, and power transmission in cold regions. With the increasing frequency of extreme weather events, phenomena such as ice and snow cover and icing pose significant threats to the safe operation of infrastructure, the proper functioning of equipment, and the safety of life and property. For example, icing on transmission lines can cause overloads and tower collapses; icy roads can cause traffic accidents; and ice on aircraft surfaces can seriously impact flight safety. Therefore, accurately simulating and testing icy and snowy environments, and subsequently developing effective anti-icing and de-icing strategies, has become a key research area in this field. Traditional ice and snow environment simulation and testing methods often have certain limitations. Some methods only focus on the influence of a single physical field, such as only considering the effect of temperature changes on the icing process, while ignoring the coupling effects of other factors such as humidity and airflow. This single-factor simulation is difficult to truly reflect the complex conditions of the actual ice and snow environment, resulting in a large deviation between the simulation results and the actual situation. In addition, some simulation methods lack an effective dynamic compensation mechanism when dealing with dynamically changing environmental parameters. When environmental parameters such as temperature and humidity fluctuate, the accuracy of key indicators such as the simulated ice layer growth rate drops significantly, and it cannot provide a reliable basis for the formulation of anti-icing treatment plans. At the same time, traditional methods often rely on empirical judgment or simple numerical calculations when identifying high-risk icing areas, making it difficult to accurately locate risk areas, making the deployment of anti-icing measures lack of pertinence and reducing anti-icing efficiency. With the development of technology, multi-physics field coupling simulation has gradually attracted attention, but existing technologies are still insufficient in extracting and analyzing the distribution characteristics of coupled fields. How to effectively integrate multi-field data such as temperature, humidity, and ice thickness, accurately generate coupling characteristic parameters that reflect the interaction between various physical fields, and how to accurately predict the evolution of the ice layer based on these characteristic parameters are still urgent problems to be solved in the current field of ice and snow environment simulation and testing. In addition, when generating anti-icing treatment solutions, existing methods are mostly limited to a single technical means, such as only using heat source heating or a single surface ice melting technology, and fail to provide a diversified and targeted solution set according to different icing conditions and environmental conditions, affecting the maximization of anti-icing effects. Summary of the Invention

[0003] The purpose of the present invention is to provide a snow and ice environment simulation test method based on multi-field coupling to solve the problems raised in the above background technology.

[0004] To achieve the above objectives, the present invention provides a method for simulating ice and snow environment testing based on multi-field coupling, the method comprising: Collect dynamic multi-field data sets in the target ice and snow environment, wherein the dynamic multi-field data sets include periodic temperature change series, relative humidity fluctuation series and real-time monitoring data of surface ice thickness; Performing multi-physics field coupling numerical simulation processing on the dynamic multi-field data set to generate coupled field distribution characteristics of the target ice and snow environment, wherein the coupled field distribution characteristics include an axial heat transfer coefficient gradient, a radial fluid velocity accumulation, and a tangential ice layer stress fluctuation coefficient; Calling a preconfigured ice and snow growth simulation model to perform ice layer evolution analysis and processing on the coupled field distribution characteristics, and generating an ice layer growth rate prediction value and a location identifier of a high-risk icing area; performing an environmental dynamic compensation correction process on the ice layer growth rate prediction value to generate a corrected ice layer growth rate prediction value, wherein the environmental dynamic compensation correction process is implemented by utilizing the correlation characteristics between the relative humidity fluctuation sequence and the heat capacity property of the material; Based on the high-risk icing area location identifier, an anti-icing treatment solution set is generated, and the anti-icing treatment solution set includes a heat source deployment strategy and a surface ice melting technology solution.

[0005] Preferably, performing multi-physics field coupling numerical simulation processing on the dynamic multi-field data set to generate coupling field distribution characteristics of the target ice and snow environment includes: Dividing the periodic temperature variation sequence into a plurality of temperature subsequences according to a preset time window, each temperature subsequence corresponding to a simulation period unit; For each temperature subsequence, perform the following operations: Constructing a three-dimensional ice layer deformation topology model based on the real-time monitoring data of the surface ice layer thickness, wherein the three-dimensional ice layer deformation topology model includes spatial coordinate information of thermal expansion field distribution, fluid velocity field distribution, and stress field distribution; Performing field coupling interaction processing on the three-dimensional ice layer deformation topology model and the temperature subsequence to generate a coupled field simulation output of the current simulation period unit; Performing cumulative integration calculation on the coupled field simulation outputs of a plurality of consecutive simulation period units to obtain the axial heat transfer coefficient gradient, radial fluid velocity accumulation, and tangential ice layer stress fluctuation coefficient; Among them, the axial heat transfer coefficient gradient is the maximum rate of change of the thermal expansion field along the axial direction, the radial fluid velocity accumulation is the total integral of the fluid velocity field in the normal direction of the ice layer contact surface, and the tangential ice layer stress fluctuation coefficient is the coefficient of variation of the stress field distribution over time.

[0006] Preferably, performing field coupling interaction processing on the three-dimensional ice layer deformation topology model and the temperature subsequence to generate a coupled field simulation output of the current simulation period unit includes: Based on the corresponding relationship between the heat conduction law and the temperature subsequence, a thermal field-temperature balance equation is established, and an initial distribution function of the heat conduction component is obtained by solving the thermal field-temperature balance equation; Based on the correlation between fluid dynamics principles and ice deformation data, a fluid-ice coupling calculation framework is constructed, which integrates dynamic adjustment parameters of fluid viscosity coefficient and ice elastic modulus; Combining the surface stress change rate with the friction effect factor, an iterative ice growth prediction process is implemented, which includes a feedback calibration mechanism for the stress increment and ice deformation. A spatial field fusion operation is performed on the output results of the initial distribution function, the fluid-ice layer coupling calculation framework, and the ice layer growth iterative prediction process to generate a three-dimensional coupled field distribution data set containing heat conduction, fluid velocity, and ice layer stress components.

[0007] Preferably, the calling of a preconfigured ice and snow growth simulation model to perform ice layer evolution analysis and processing on the coupled field distribution characteristics to generate an ice layer growth rate prediction value and a high-risk icing area location identifier includes: Inputting the axial heat transfer coefficient gradient into the initial characteristic analysis layer of the ice and snow growth simulation model, and determining the distribution coordinates of the heat concentration area and the heat amplitude change trajectory through the heat concentration factor calculation unit; Inputting the radial fluid velocity accumulation into the secondary characteristic analysis layer of the ice and snow growth simulation model, performing fluid scour damage accumulation calculation, and generating the melting probability and crack growth rate prediction value of the ice surface; Inputting the tangential ice layer stress fluctuation coefficient into the final characteristic analysis layer of the ice and snow growth simulation model, and calculating the material wear depth of the friction surface and the ice layer roughness evolution data based on the surface ice layer degradation model; The thermal amplitude variation trajectory, the melting probability, and the material wear depth are integrated to generate a comprehensive evolution index of the target ice and snow environment, and a predicted value of the ice layer growth rate is determined based on a comparison result of the comprehensive evolution index with a preset threshold library; Based on the spatial field superposition results of the distribution coordinates, the crack growth rate prediction value and the ice layer roughness evolution data, the geometric position information of the heat concentration area, the crack path and the high-risk wear area is identified.

[0008] Preferably, performing environmental dynamic compensation correction processing on the ice layer growth rate prediction value to generate a corrected ice layer growth rate prediction value includes: Extracting the extreme humidity values ​​and humidity change frequencies in the relative humidity fluctuation sequence, and calculating the dynamic correction amount of the material heat capacity property as the humidity fluctuates; performing a thermal stress compensation calculation on the axial heat transfer coefficient gradient according to the dynamic correction amount to generate a compensated heat transfer coefficient gradient; Based on the correlation between the humidity change frequency and the phase change characteristics of the material, performing phase change damage correction processing on the radial fluid velocity accumulation to generate a corrected fluid velocity accumulation; performing surface strength adaptation adjustment processing on the tangential ice layer stress fluctuation coefficient according to the material thermal conductivity variation data under extreme humidity to generate an adjusted ice layer stress fluctuation coefficient; The compensated thermal conductivity gradient, the corrected fluid velocity accumulation and the adjusted ice layer stress fluctuation coefficient are input into the ice and snow growth simulation model for recalculation to generate an ice layer growth rate prediction value including environmental factor compensation.

[0009] Preferably, performing thermal stress compensation calculation on the axial thermal conductivity gradient according to the dynamic correction amount to generate a compensated thermal conductivity gradient includes: Obtaining the initial heat capacity property of the target ice and snow environment under reference humidity conditions and the dynamic correction amount, and establishing a heat capacity-humidity correlation function; Calculating an axial thermal stress compensation value according to the heat capacity-humidity correlation function, wherein the axial thermal stress compensation value is the product of the humidity change amplitude and the heat capacity change amplitude; superimposing the axial thermal stress compensation value into the calculation process of the axial thermal conductivity gradient to generate an intermediate thermal conductivity gradient value including thermal stress compensation; A thermal relaxation effect compensation process is performed on the intermediate value of the thermal conductivity gradient, wherein the thermal relaxation effect compensation process is implemented by using a multiplication factor of a material thermal relaxation curve and a humidity retention time.

[0010] Preferably, generating an anti-icing treatment solution set based on the high-risk icing area location identifier includes: Calculating an optimal heat source deployment path based on the identification information of the heat concentration area, wherein the optimal heat source deployment path is achieved by adjusting the power distribution ratio of adjacent heat sources; Constructing a surface ice melting technical solution based on the identification information of the crack path, wherein the surface ice melting technical solution includes the selection of a microwave heating area and the optimization configuration of energy parameters; generating a lubricant application strategy based on the identification information of the high-risk wear area, wherein the lubricant application strategy dynamically adjusts the application frequency and application thickness according to the wear rate prediction value; The optimal heat source deployment path, the surface ice melting technical solution and the lubricant application strategy are prioritized to generate an anti-icing solution set including execution timing and technical parameters.

[0011] Preferably, the surface ice melting technology solution comprises: Extracting geometric characteristic data of the crack path and calculating the path curvature radius and expansion direction angle; selecting a coverage density parameter of the microwave heating spot according to the curvature radius, wherein the coverage density parameter is in a nonlinear inverse proportional relationship with the curvature radius; Adjusting the direction of application of microwave heating energy based on the expansion direction angle so that the energy application direction forms a preset angle with the crack path direction; Dynamically adjust the duration of microwave heating pulses based on the test data of the material's surface thermal conductivity to ensure that the heating energy is below the critical value of the material's melting point; Generate an ice melting parameter configuration table including coverage density parameters, application direction and pulse duration.

[0012] Preferably, the method further comprises: Collecting actual ice thickness and crack propagation length data of the target ice and snow environment within a preset verification period; performing deviation analysis on the actual ice thickness and the predicted ice thickness to generate a thickness error correction coefficient; Performing a time series comparison process on the crack propagation length and the predicted propagation speed to generate a length error correction coefficient; Adjusting the parameter weights of the ice and snow growth simulation model according to the thickness error correction coefficient and the length error correction coefficient to generate an optimized ice and snow growth simulation model; The optimized ice and snow growth simulation model is applied to ice and snow environment simulation tasks in subsequent test batches.

[0013] Preferably, after generating the anti-icing treatment solution set, the method further includes: Inputting the anti-icing treatment plan set into the environmental simulation control system to perform real-time scheduling processing; updating the acquisition parameters of the dynamic multi-field data set according to the scheduling result feedback; Restarting the multi-physics field coupling numerical simulation process using the updated acquisition parameters to generate an iterative coupling field distribution feature; The anti-icing treatment scheme set is optimized based on the iterative coupling field distribution characteristics, and a final simulation test report is output.

[0014] Compared with the prior art, the present invention has the following beneficial effects: This multi-field coupled ice and snow environment simulation test method collects dynamic, multi-field data sets from the target ice and snow environment, encompassing periodic temperature change sequences, relative humidity fluctuation sequences, and real-time monitoring data on surface ice thickness. This method comprehensively captures the multi-dimensional and dynamic information of the ice and snow environment, providing rich and realistic basic data for subsequent simulation analysis. Compared with traditional methods that only consider a single physical field, this multi-field data collection can more comprehensively reflect the complex characteristics of the actual ice and snow environment and avoid simulation bias caused by missing data. Applying multi-physics coupled numerical simulations to dynamic multi-field data sets generates coupled field distribution characteristics, including axial heat transfer coefficient gradients, radial fluid velocity accumulations, and tangential ice layer stress fluctuation coefficients. These coupled characteristic parameters can deeply reveal the interaction mechanisms between various physical fields, such as temperature, humidity, and ice layer. These coupled characteristic parameters quantify the coupling effects between multiple fields from different dimensions, providing key clues for understanding the inherent laws of ice and snow formation and development. This makes the simulation process more consistent with actual physical processes and overcomes the shortcomings of traditional methods that lack consideration of multi-field coupling effects. A preconfigured ice and snow growth simulation model analyzes the coupled field distribution characteristics to analyze ice evolution, generating predicted ice growth rates and identifying high-risk icing areas. This allows for precise deduction of ice layer dynamics and accurate location of risk areas. The ice and snow growth simulation model analyzes coupled field distribution characteristics and leverages the rich information provided by multi-field coupling, making ice growth rate predictions more scientific and identifying high-risk icing areas more realistic. This addresses the inaccurate positioning issues often associated with traditional methods that rely on empirical judgment or simple calculations. The ice growth rate predictions are corrected using dynamic environmental compensation, leveraging the correlation between relative humidity fluctuations and material heat capacity. This correction effectively improves the accuracy of the predictions in dynamic environments. This dynamic compensation mechanism allows for timely adjustments to predictions when environmental parameters fluctuate, preventing bias caused by environmental changes. This makes ice growth rate predictions more adaptable to dynamic changes in the actual environment and provides a more reliable reference for developing anti-icing strategies. Based on the location of high-risk icing areas, a collection of anti-icing solutions, including heat source deployment strategies and surface ice melting technology solutions, is generated. This provides diverse and targeted solutions for different high-risk areas and icing conditions. This collection of solutions is no longer limited to a single anti-icing method. Instead, it flexibly selects appropriate anti-icing technologies based on specific icing characteristics and environmental conditions. This improves the relevance and effectiveness of anti-icing measures, overcoming the limitations of traditional anti-icing methods, which often rely on a single method and have limited effectiveness. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1This is a working principle diagram of the ice and snow environment simulation test method based on multi-field coupling according to the present invention; Figure 2 Flowchart for multi-physics field coupling numerical simulation processing; Figure 3 Flowchart for ice layer evolution analysis and processing; Figure 4 Construct a flow chart for the surface ice melting technology solution. DETAILED DESCRIPTION

[0016] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0017] See also Figure 1 The present invention provides a method for simulating and testing an ice and snow environment based on multi-field coupling, the method comprising: This system collects dynamic multi-field data sets from the target ice and snow environment, including periodic temperature change sequences, relative humidity fluctuation sequences, and real-time monitoring data on surface ice thickness. Multi-physics field coupling numerical simulation is then performed on this dynamic multi-field data set to generate the coupled field distribution characteristics of the target ice and snow environment, including the axial heat transfer coefficient gradient, radial fluid velocity accumulation, and tangential ice stress fluctuation coefficient. A pre-configured ice and snow growth simulation model is then used to perform ice layer evolution analysis on the coupled field distribution characteristics, generating ice layer growth rate predictions and identification of high-risk icing areas. The predicted ice layer growth rate is then corrected for environmental dynamic compensation to generate a corrected ice layer growth rate prediction. Based on the identification of high-risk icing areas, a set of anti-icing solutions is generated, including heat source deployment strategies and surface ice melting technology solutions.

[0018] Example 1: See Figure 2 , a multi-physics coupled numerical simulation process for dynamic multi-field data sets. The process begins with the segmentation of the periodic temperature variation sequence, dividing the continuously collected temperature data into several discrete temperature subsequences according to preset time windows. Each temperature subsequence represents a complete temperature fluctuation cycle and serves as the basic unit for subsequent simulation calculations. The division of the time window is determined by the actual temperature variation characteristics of the target ice and snow environment, and the minimum time span that can reflect the complete freeze-thaw cycle is usually selected.

[0019] For each temperature subsequence, the system first uses real-time monitoring data of the surface ice thickness to construct a three-dimensional ice deformation topology model. This model divides the ice structure into a finite number of calculation units through a spatial discretization method. Each unit records the spatial coordinate information of the thermal expansion field distribution, fluid velocity field distribution, and stress field distribution. The thermal expansion field data is derived from the product relationship between the material thermal expansion coefficient and the local temperature change. The fluid velocity field data is obtained by interpolation of the boundary layer velocity measurement results, and the stress field data is inverted based on the strain monitoring results inside the ice layer. The grid density of the three-dimensional model is automatically adjusted according to the gradient of the ice thickness change, and a finer grid division strategy is adopted in the area of ​​thickness mutation.

[0020] During the field coupling interaction processing phase, the three-dimensional ice layer deformation topology model is coupled with the current temperature subsequence for multi-physics field coupling calculations. The thermal field-temperature balance equation is established based on Fourier's heat conduction law, taking into account the constitutive relationship between temperature gradient and heat flux density. The equation solution process adopts an implicit difference format to increase the time step while ensuring the stability of the calculation. The initial distribution function of the heat conduction component is obtained by iterative solution, and its spatial distribution characteristics reflect the initial thermal influence of the temperature field on the ice layer structure. The fluid-ice layer coupling calculation framework integrates the Navier-Stokes equations and the ice layer elastic deformation equation, and realizes fluid-solid coupling calculations through an alternating solution strategy. The fluid viscosity coefficient adopts a temperature-dependent dynamic adjustment algorithm, and the ice layer elastic modulus is updated in real time according to the measured stress-strain curve.

[0021] The iterative prediction process for ice layer growth introduces the surface stress change rate as a key control parameter and combines it with the friction effect factor to construct a closed-loop feedback system. The stress increment calculation takes into account the influence of ice crystal growth orientation and is described using an anisotropic material model. Feedback calibration of ice layer deformation is achieved by comparing the deviation between the simulated and measured values. When the deviation exceeds a threshold, a parameter adjustment mechanism is automatically triggered. The spatial field fusion operation vector-superimposes the heat conduction component, fluid velocity field, and ice layer stress field, and uses a weighted average algorithm to address the coupling effects of different physical fields. Weight coefficients are dynamically assigned based on the degree of influence of each field on the evolution of the ice layer. The heat conduction field has a higher weight during the freezing phase, while the fluid field dominates during the melting phase.

[0022] The cumulative integration calculation of continuous simulation cycle units uses a time series analysis method. The calculation of the axial heat transfer coefficient gradient selects the extreme value of the thermal expansion field along the main axis, and the maximum value of the spatial change rate is obtained through differential operation. The radial fluid velocity accumulation is obtained by surface integration of the normal velocity component, and the integration area covers the entire ice layer contact surface. The calculation of the tangential ice layer stress fluctuation coefficient is based on the statistical characteristics of the stress field time series, and the ratio of the standard deviation to the mean is used to represent the stress fluctuation intensity. The integration process of multi-cycle data introduces an attenuation factor, and the contribution of data from earlier cycles gradually decreases with simulation time.

[0023] Post-processing of coupled-field simulation outputs includes data normalization and outlier correction. Normalization converts the dimensions of different physical fields into dimensionless parameters, facilitating subsequent model invocation. An outlier correction algorithm detects and removes simulation results that significantly deviate from the normal range, and missing data are supplemented by interpolation from adjacent nodes. The resulting coupled-field distribution feature dataset contains complete spatial distribution information and temporal evolution characteristics, providing accurate input parameters for the snow and ice growth simulation model.

[0024] The update mechanism of the three-dimensional ice layer deformation topology model adopts a dynamic adaptive strategy. When the monitored ice layer thickness change exceeds the preset threshold, the system automatically triggers the model reconstruction process. The reconstruction process retains the topological structure of the original model and only updates the node coordinates and physical field parameters. The correction of the thermal expansion field data takes into account the influence of the latent heat of the material phase change, and adopts a special constitutive relationship in the ice-water phase transition temperature range. The update of the fluid velocity field introduces a turbulence model and considers the influence of eddy effects in high flow velocity areas. The stress field data is recalculated using an incremental method, superimposing the current strain increment on the stress state at the previous moment.

[0025] The adaptive time step adjustment algorithm dynamically optimizes the time step based on computational convergence. When increasing nonlinear effects are detected, the time step is automatically reduced to ensure computational accuracy; during periods of gradual field changes, the step size is appropriately increased to improve computational efficiency. The step size adjustment strategy comprehensively considers multiple factors, including temperature change rate, flow velocity gradient, and stress fluctuation amplitude, and implements intelligent control through a fuzzy logic algorithm. This dynamic time management approach significantly improves overall simulation efficiency while maintaining computational accuracy.

[0026] The core of multi-physics field coupling numerical simulation lies in the precise description of the interaction mechanisms between the various physical fields. The coupling of the thermal field and the flow field is achieved through the Businesq approximation, taking into account the buoyancy effect caused by temperature changes. The coupling of the flow field and the stress field uses the arbitrary Lagrangian-Euler method to accurately capture the mechanical effects of fluid loads on the ice structure. The coupling of the thermal field and the stress field introduces the theory of thermoelasticity to analyze the thermal stress distribution caused by temperature gradients. The solution to the fully coupled three-field problem adopts a step-by-step iterative strategy, gradually approaching the true solution through alternating calculations between the fields.

[0027] Example 2: See Figure 3, the ice and snow growth simulation model analyzes and processes the evolution of the ice layer based on the coupled field distribution characteristics. This process uses the coupled field distribution characteristics generated in Example 1 as input data, and realizes the accurate prediction of the dynamic evolution of the ice layer through a multi-level feature analysis mechanism. The axial thermal conductivity coefficient gradient is first input into the initial feature analysis layer for processing. The analysis layer contains a heat concentration factor calculation unit, and a spatial gradient analysis method is used to identify abnormal heat conduction areas. The determination of the heat concentration area is based on the local extreme value characteristics of the heat conductivity coefficient gradient. The entire calculation domain is scanned by a sliding window algorithm, and the spatial coordinates where the gradient change exceeds the threshold are marked. The tracking of the thermal amplitude change trajectory adopts a time correlation algorithm to correlate and match the coordinates of the heat concentration area in consecutive time steps to construct the development path of the heat affected area.

[0028] The secondary feature analysis layer processes radial fluid velocity accumulation data and includes a module for calculating fluid scour damage accumulation. The module's operating mechanism considers the interaction between fluid velocity and the ice surface, with the normal component of the velocity vector used to calculate the local scour intensity. The calculation of the melt probability incorporates a phase transition dynamics model to analyze the microstructural changes on the ice surface under fluid scour. The crack propagation rate prediction utilizes the principles of fracture mechanics, based on the interaction between fluid pressure fluctuations and inherent defects in the ice layer. The calculation considers the stress intensity factor at the crack tip and, combined with the material's fracture toughness parameters, derives the dynamic crack propagation law. The output of the secondary analysis layer includes a spatially distributed melt probability field and a crack propagation velocity field, reflecting the erosion effect of fluid action on the ice structure.

[0029] The final feature analysis layer processes the stress fluctuation coefficient of the tangential ice layer and integrates the surface ice degradation model. The calculation framework of the degradation model includes a wear depth prediction algorithm and a roughness evolution analysis module. The calculation of wear depth is based on the cumulative damage theory, which correlates the historical data of stress fluctuations with the fatigue properties of the material. The roughness evolution analysis uses surface morphology reconstruction technology to calculate the surface contour changes through the micro-plastic deformation caused by stress fluctuations. The analysis process considers the influence of ice crystal orientation on wear characteristics and adopts an anisotropic material response model to improve prediction accuracy. The output of the final analysis layer includes a wear depth distribution map and a roughness evolution curve, which characterizes the long-term effect of mechanical stress on the ice surface.

[0030] The generation of the comprehensive evolution index utilizes a multi-field coupling weighted algorithm, spatially superimposing the thermal amplitude variation trajectory of the initial characteristic resolution layer, the melting probability of the secondary characteristic resolution layer, and the material wear depth of the final characteristic resolution layer. The weighting mechanism considers the relative contributions of different physical fields to the evolution of the ice layer. The thermal field weight automatically increases in low-temperature environments, while the fluid field weight increases accordingly during the melting phase. A preset threshold library contains reference data for a variety of typical ice and snow environmental conditions. A pattern matching algorithm is used to determine the threshold combination that best matches the current environmental state. The comparison process utilizes fuzzy logic reasoning to address the nonlinear relationship between the index and the threshold. Probabilistic statistical methods are used to determine the final predicted value of the ice growth rate, calculating the probability distribution of different growth rates.

[0031] The identification process of high-risk areas uses spatial clustering analysis technology. The distribution coordinates of the heat concentration area are first converted into three-dimensional point cloud data, and the critical heat-affected zone is identified through a density clustering algorithm. The predicted results of the crack path are spatially correlated with the fluid velocity field to mark the dangerous sections where the flow velocity is consistent with the crack propagation direction. The determination of high-risk areas for wear is based on the gradient change of wear depth, and an edge detection algorithm is used to identify the boundary areas of sudden wear changes. The integration of geometric position information uses a spatial overlay analysis method to project areas of different risk types into a unified coordinate system. The output of the identification results includes a risk level distribution map and three-dimensional spatial position data, which supports the precise positioning of subsequent anti-icing measures.

[0032] An adaptive learning algorithm is used to optimize the parameters of the feature resolution layer. The heat concentration threshold of the initial feature resolution layer is dynamically adjusted based on the ambient temperature, with a more sensitive threshold setting used under low-temperature conditions. The fluid erosion model parameters of the secondary feature resolution layer are automatically updated as the flow rate changes, accounting for turbulence effects at high flow rates. The wear calculation coefficient of the final feature resolution layer is modified in real time based on the material state, taking into account the impact of changes in ice density on wear characteristics. The optimization process utilizes a closed-loop feedback mechanism, comparing predicted results with actual observed data and automatically adjusting model parameters to reduce prediction bias.

[0033] The time-marching mechanism for ice layer evolution analysis utilizes a variable step-size control strategy. During periods of drastic environmental parameter changes, a smaller time step is used to ensure accurate analysis of transient processes. During steady-state evolution, the time step is appropriately increased to improve computational efficiency. A step-size adjustment algorithm monitors the rate of change of each physical field and automatically triggers a step-size reduction when a sudden change in field intensity gradient is detected. A hybrid implicit-explicit time discretization scheme is employed, with implicit discretization used for rapidly changing quantities to ensure stability and explicit discretization used for slowly changing quantities to improve computational efficiency.

[0034] Adaptive meshing technology is used to manage spatial resolution. The computational grid is automatically refined in areas of heat concentration and near crack propagation paths to improve local feature resolution. A coarser mesh is used in low-gradient regions to optimize computational resource allocation. The mesh refinement criterion is based on the gradient change rate of the field variables, triggering local refinement when the spatial derivative exceeds a threshold. Dynamic adjustment of the mesh topology utilizes an octree data structure, supporting fast spatial queries and neighborhood searches.

[0035] Multi-source data fusion utilizes a Bayesian inference framework. Numerical simulation results are probabilistically fused with real-time monitoring data to address measurement errors and model uncertainties. The data assimilation process considers the reliability differences between different data sources, assigning greater weight to high-precision measurement data. The confidence level of the fusion results is expressed through a posterior probability distribution, supporting reliability assessment for risk-based decision-making.

[0036] The risk visualization system uses augmented reality technology to present analysis results. Heat concentration areas are displayed as a red heat map overlay, crack paths are marked with blue wireframes, and high-risk wear areas are rendered with yellow semi-transparent volumetric rendering. Visualization parameters can be adjusted interactively, supporting multi-level observation from overall distribution to local details. The roaming function of the 3D scene allows users to view the three-dimensional distribution of risk areas from any perspective, assisting in the planning and deployment of anti-icing measures.

[0037] Example 3: See Figure 4 This process uses dynamic environmental compensation to correct ice growth rate predictions. This process establishes a compensation mechanism for the effects of humidity fluctuations on material properties in icy environments, improving prediction accuracy through a multi-stage correction algorithm. Preprocessing of relative humidity fluctuation sequences includes extreme value detection and frequency analysis. A sliding time window statistical method is used to extract extreme humidity values ​​from the sequence, and Fourier transform decomposition is used to obtain the characteristic frequencies of humidity variations. The dynamic correction calculation for the material's heat capacity properties considers the nonlinear relationship between humidity and heat capacity, establishing the following compensation relationship:

[0038] in, represents the heat capacity correction, is the humidity change amplitude, is the basic heat capacity coefficient of the material, is the humidity sensitivity factor, These parameters are obtained through material property experimental calibration and have different values ​​for different types of ice and snow covered surfaces.

[0039] Thermal stress compensation calculations utilize an incremental superposition method to convert heat capacity corrections into compensation values ​​for the axial thermal conductivity. The calculation considers the coupled effects of temperature gradients and thermal expansion coefficients, introducing an additional source term into the heat conduction equation to account for the effects of humidity changes. The compensated thermal conductivity gradient is solved iteratively, with each iteration updating the humidity-related material parameters. A time delay factor is introduced to compensate for thermal relaxation effects, accounting for the hysteresis in the material's response to humidity changes. This delay time is calculated based on the material's thermal diffusion characteristics and multiplied with the humidity hold time to correct for thermal conductivity.

[0040] Phase change damage correction compensates for the effects of humidity on fluid velocity accumulation. A model correlating humidity change frequency with material phase change properties is constructed based on phase change kinetics theory, analyzing the degree of completeness of the ice-water phase change under different humidity fluctuation cycles. A correction algorithm detects regions of incomplete phase change and adjusts the release rate of latent heat during fluid calculations. The phase change front is tracked using the level set method, with the signed distance function describing the position of the phase change interface. The distribution of corrections to the fluid velocity field is correlated with the geometric characteristics of the phase change interface, with a greater correction weight applied to regions near the interface.

[0041] Surface strength adaptation and adjustment analyzes changes in material thermal conductivity under extreme humidity conditions. The adjustment algorithm establishes a humidity-thermal conductivity correlation database and uses interpolation to obtain equivalent thermal conductivity parameters under the current humidity conditions. The recalculation of the stress fluctuation coefficient considers the thermal stress redistribution caused by uneven thermal conductivity and employs an equivalent homogenization method to address local parameter changes. The adjusted stress field utilizes virtual node technology to achieve continuous parameter transition, avoiding numerical calculation instability.

[0042] The prediction of ice growth rate with environmental compensation utilizes a multi-field coupled iterative computational framework. The compensated thermal conductivity gradient, the corrected fluid velocity accumulation, and the adjusted ice stress fluctuation coefficient are used as input parameters, and the ice and snow growth simulation model is re-solved. The iterative process monitors the convergence of key physical quantities and terminates the calculation when the difference between the prediction results of two consecutive iterations falls below a threshold. The final output of the ice growth rate prediction includes humidity compensation effects, reflecting the actual impact of dynamic changes in ambient humidity on ice evolution.

[0043] The adaptive adjustment mechanism for the humidity sensitivity factor is dynamically optimized based on environmental monitoring data. The adjustment algorithm analyzes the relationship between historical humidity data and material response and updates the sensitivity factor parameters through regression analysis. The optimization process also factors in seasonal variations, establishing differentiated models for the impact of humidity during different time periods. The sensitivity factor is lower in dry winter conditions, while a more sensitive adjustment strategy is employed during periods of high humidity fluctuations in spring and autumn.

[0044] The spatial and temporal distribution of the latent heat correction for phase change is addressed using a nonuniform compensation approach. In the spatial dimension, different correction intensities are assigned based on the magnitude of local humidity variations, while in the temporal dimension, the duration of the correction is adjusted based on the frequency of humidity fluctuations. This four-dimensional compensation strategy accurately captures the spatial heterogeneity and temporal dynamics of the phase change process, avoiding the error accumulation associated with traditional uniform compensation methods.

[0045] The material parameter database is constructed using a hierarchical storage structure. Basic material parameters, including inherent properties such as heat capacity and thermal conductivity, are stored in the core layer. Environmental compensation parameters are stored in the middle layer, recording parameter correction rules for different humidity conditions. Temporary calculation parameters are stored in the cache layer, preserving intermediate results of the current simulation state. The database update mechanism uses version control, retaining historical parameter sets for retrospective analysis.

[0046] Stability control for coupled thermomechanical calculations utilizes a step-by-step relaxation technique. The humidity compensation is applied incrementally, with each incremental step performing a separate thermodynamic balance calculation. The relaxation factor is dynamically adjusted based on the convergence rate, automatically reducing the incremental step size when the risk of divergence is high. This progressive compensation method effectively handles strong nonlinearities and ensures computational stability under complex humidity conditions.

[0047] The Extreme Conditions Identification and Processing module monitors for abnormal humidity fluctuations. When humidity fluctuations exceeding safety thresholds are detected, the system automatically switches to enhanced compensation mode. This mode employs a finer time step and spatial grid, and utilizes higher-order numerical formats in critical areas. The duration of the abnormal event is recorded in a log for subsequent optimization of the compensation model.

[0048] The effectiveness of humidity compensation is verified using a dual validation method. Theoretical validation uses energy conservation analysis to check the energy balance in the compensation calculation. Numerical validation uses a grid convergence test to analyze the consistency of compensation results at different spatial resolutions. Deviations discovered during the validation process trigger automatic fine-tuning of compensation parameters, forming an adaptive error correction loop.

[0049] Assimilation of real-time monitoring data incorporates the latest humidity observations into compensation calculations. The assimilation algorithm uses an optimal interpolation method to update field variables while ensuring computational continuity. The data assimilation interval automatically adjusts based on the rate of humidity change, with a higher assimilation frequency during periods of rapid fluctuation. This dynamic data fusion mechanism ensures that compensation calculations consistently track actual environmental conditions.

[0050] The output of compensation calculations contains complete field distribution information and statistical characteristics. Spatially distributed data records the physical field state after compensation, while statistical characteristics extract the variation range and trend indicators of key parameters. The output format supports direct access to various post-processing tools, facilitating further analysis and visualization.

[0051] Example 4: This article focuses on the generation process of a set of anti-icing solutions, using specific examples to illustrate the complete process from identifying high-risk icing areas to finalizing the solution. Using the snow and ice disaster prevention and control of an overhead high-voltage transmission line as an application scenario, the system first receives the risk area identification data from Example 2, which includes three typical risk areas: the heat concentration area at the upper end of the insulator string hardware, the crack path between the composite insulator sheds, and the high-risk wear area at the conductor suspension point.

[0052] The optimal heat source deployment path for the heat-concentrated area at the upper end of the insulator string hardware was calculated based on the following parameters: the hardware was made of cast aluminum alloy, the surface temperature fluctuated periodically between -12°C and -5°C, and the spacing between adjacent heating elements was constrained to within 15 cm. The system-generated heat source configuration adopted a staggered layout. Table 1 shows the heat source deployment parameters for this area.

[0053] Table 1: Parameters of heat source deployment at the top of insulator string.

[0054] Heat source number Installation position (cm from the top) Rated power (W) Working cycle (min) Temperature trigger threshold (°C) HS-01 3.2 120 8 -8.5 HS-02 18.7 80 12 -9.2 HS-03 34.5 150 6 -7.8 HS-04 50.0 100 10 -8.0 Microwave ice melting technology is used to treat the crack path between the sheds of composite insulators. Systematic analysis of the crack geometry shows that the main crack path extends in a spiral shape, with a curvature radius varying in the range of 2.1-5.3 cm, and the extension direction forms an average angle of 22° with the shed axis. The coverage density of the microwave heating spots is dynamically adjusted according to the curvature radius. At the minimum curvature radius, 3 heating points are set per square centimeter, and the larger curvature area is reduced to 1.5 heating points per square centimeter. The energy application direction is set to a 15° deflection angle with the crack path. This angle design can effectively block the crack extension and avoid excessive concentration of energy in a single direction. The pulse duration is controlled between 200-350 milliseconds. The specific value is dynamically adjusted according to the real-time surface temperature of the silicone rubber material to always maintain the heating temperature below the critical softening point of the material.

[0055] A lubricant application strategy is adopted in the high-risk wear area of ​​the conductor suspension point. System analysis shows that there are two wear modes in this area: periodic vibration wear in the vertical direction and breeze vibration wear in the horizontal direction. The lubricant formula uses low-temperature silicon-based composite materials, and a correlation model is established between the application frequency and the conductor vibration amplitude. When the vibration acceleration exceeds 0.5g, the application interval is automatically shortened from the standard 4 hours to 2 hours. The application thickness is controlled within the range of 0.1-0.3mm, and the upper limit value is used in extremely cold weather conditions to enhance the protection effect. The movement trajectory of the lubricant nozzle is optimized according to the wear distribution map, and a spiral covering path is used in key areas to ensure a uniform and complete protective layer.

[0056] Prioritizing options utilizes a multi-criteria decision-making approach, taking into account key factors: the hazard level of the risk area, the immediate effectiveness of the treatment measures, implementation costs, and ease of maintenance. In the case of transmission lines, the system prioritizes heat concentration in insulator strings, as they could potentially trigger insulation failure. Crack path treatment was ranked second, and conductor suspension point lubrication, a routine protective measure, ranked third. The implementation schedule is as follows: heat source deployment is completed 24 hours before the onset of the cold wave, the microwave de-icing system is activated when the temperature drops to -5°C, and lubricant application is automatically applied according to a pre-set cycle.

[0057] Optimal configuration of technical parameters is achieved through a feedback loop. Heat source power allocation uses an adaptive algorithm, with initial settings calculated based on a standard heat conduction model. After actual operation, adjustments are made dynamically based on infrared temperature measurement results. Monitoring revealed a 5°C temperature deviation in the area covered by the HS-03 heat source. The system automatically increased its power from the initial 150W to 170W, while simultaneously reducing the HS-02 power by 10W to maintain overall balance. Adjustment of microwave heating parameters is more refined, with the energy output of each heating spot independently controllable. By monitoring changes in the crack tip strain rate in real time, the pulse parameters at each point are dynamically optimized.

[0058] Material compatibility testing was conducted throughout the design process. Contact corrosion testing was performed before heat source installation to confirm chemical compatibility between the heating element and the cast aluminum alloy. The microwave frequency was selected at 2.45 GHz to avoid overlap with the molecular resonance frequency of the composite insulator. The lubricant formulation underwent a 48-hour low-temperature aging test to confirm good adhesion at -30°C. All material selection was based on a database of over 200 commonly used materials in power equipment.

[0059] Differentiated treatment of microenvironmental factors reflects the refined nature of the solution. Different heat source configurations are used for different orientations of the same tower: the sun-facing side reduces the power setting by 15% to account for the auxiliary heating effect of solar radiation; the shady side increases the power margin by 10% to cope with persistent low temperatures. An altitude correction factor is applied to all heating parameters, with the power increasing by 1.2% for every 100-meter increase in altitude. The influence of local wind speed is incorporated into the microwave heating direction adjustment. When the wind speed exceeds 5 m / s, the energy application direction is deflected 8° downwind to compensate for the wind chill effect.

[0060] The layout of the condition monitoring network was designed in parallel with the treatment plan. Temperature sensors were placed at the boundaries of the heat source's radius, using wireless networking technology for real-time data transmission. A high-frequency acoustic emission sensor array was used for crack monitoring to capture stress wave signals generated by crack propagation. Wear monitoring employed resistive wear sensors, quantifying the degree of wear through impedance changes. Monitoring data was sampled at one-minute intervals and analyzed and processed in real time by edge computing nodes, providing data support for solution adjustments.

[0061] A cost-benefit analysis model evaluated the economic viability of different configuration options. The heat source deployment scheme compared three technical routes: resistance heating, induction heating, and photovoltaic heating, ultimately selecting the resistance heating option with the best operating and maintenance costs. A quantitative relationship was established between the energy consumption of the microwave de-icing system and its anti-icing effectiveness to determine the optimal energy consumption control curve. A lubricant usage prediction model comprehensively considers weather trends and equipment operating history to minimize material consumption. All economic analyses are based on full lifecycle cost calculations, including factors such as initial investment, operating energy consumption, and maintenance costs.

[0062] Emergency response plans for abnormal situations are embedded in the master plan. When a heat source failure is detected, compensation mode is automatically activated for adjacent heat sources, and a maintenance alarm is issued. In the event of a microwave system failure, the system switches to backup heating mode, using localized hot air compensation. In the event of a lubricant supply interruption, a mechanical ice scraper is activated as a temporary alternative. These plans, pre-designed through fault tree analysis, cover over 90% of typical failure scenarios.

[0063] Solution version management records all adjustments and decision-making processes. Each parameter modification generates a new solution version, documenting the reason for the change and the expected results. A version comparison tool visually displays technical parameter differences between solutions and enables quick rollback to historical versions. This strict version control ensures transparency and traceability throughout the solution evolution process.

[0064] The operator interface is designed to align with engineering practices. Heat source installation steps are broken down into detailed process cards, including key parameters such as torque values ​​and wiring sequences. The microwave equipment commissioning interface offers a simple wizard mode to guide operators through the calibration process. The lubricant filling system features voice prompts, providing real-time guidance on filling volume and position control. These user-friendly features reduce the need for specialized operator skills and improve the accuracy and efficiency of solution implementation.

[0065] Example 5: Closed-loop optimization and iterative update mechanism for the ice and snow environment simulation test system. This process is based on cross-validation of previous simulation test results and actual monitoring data, and achieves continuous improvement of model parameters through dynamic feedback adjustment. The system collects measured ice data in the target area during a typical cold wave weather cycle, including ice thickness distribution on the insulator surface, cross-sectional morphology of ice-covered conductors, and ice growth at hardware connections. This measured data is acquired using 3D laser scanning and image recognition technology, with a spatial resolution of millimeters and a temporal sampling interval of 30 minutes.

[0066] Ice thickness deviation analysis utilizes a spatial grid comparison method. The three-dimensional data field output by the predicted ice thickness model is aligned with the measured point cloud data, and the relative error percentage is calculated within each grid cell. The error distribution diagram shows that the prediction deviation in the insulator shed groove area is relatively large, reaching an average of 12.3%, while the prediction accuracy in flat surface areas is higher, with an error within 5%. This spatial difference in the system prompts the system to focus on optimizing the ice thickness calculation algorithm for complex geometric structures. The generation of the error correction coefficient takes into account the influence of local microenvironmental factors, and an additional air entrapment factor is introduced in the shed groove area with poor ventilation.

[0067] A dynamic time warping algorithm was used to compare the crack growth data in time series. This method aligns the measured and predicted crack length curves, even when nonlinear deformations exist on the time axis. Analysis revealed that during periods of rapid temperature fluctuations, the prediction model's response to accelerated crack growth experienced a lag of approximately two hours. This lag was compensated for by adjusting the material response time constant within the model. The new parameter settings reduced the phase difference between the predicted curve and the measured data to less than 30 minutes. The length error correction factor was calculated segmentally based on the different rates of temperature change, with a greater correction weight applied during periods of rapid temperature drops.

[0068] The model parameter weights are adjusted using a sensitivity-guided optimization strategy. The system automatically tests the marginal impact of each input parameter on the output, identifying that thermal conductivity and surface roughness parameters contribute most significantly to ice thickness prediction. The weight update algorithm maintains the integrity of core physical relationships and adjusts only the relative importance of minor parameters. While maintaining the underlying thermodynamic framework, the optimized model improves prediction accuracy for insulator sheds by approximately 8 percentage points. This constrained parameter optimization avoids overfitting and ensures the model's generalization capabilities under unknown environmental conditions.

[0069] The optimized model was validated using a cross-validation method. Historical data was divided into a training set and a test set. The training set was used for parameter adjustment, and the test set was used to evaluate the effectiveness of the improvements. The validation process monitored the model's improvements in spatial resolution and temporal response speed. The new model demonstrated greater adaptability to sudden snowfall events and was able to more quickly capture changing trends in snow accumulation patterns. Validation results also showed that the optimized model had a stronger ability to distinguish between different types of ice and snow mixtures, and was able to more accurately identify different types of ice accumulation, such as rime and hoarfrost.

[0070] The real-time scheduling of the environmental simulation control system utilizes an event-driven architecture. The anti-icing solution is converted into an executable sequence of control instructions, including heater start and stop timing, de-icing robot motion trajectory, and chemical spray dosage control. These instructions are distributed to field devices via the Industrial Internet of Things (IIoT) platform, with execution status transmitted back in real time via the 4G / 5G network. The scheduling algorithm dynamically adjusts instruction parameters based on device feedback. If a power anomaly is detected in a heating unit, heating tasks are automatically reallocated to surrounding units to compensate for thermal field distribution.

[0071] Data collection parameters are updated based on the principle of maximizing information entropy. The system analyzes the coverage completeness of the current dataset and prioritizes increasing monitoring density in areas of space and time where information is scarce. For example, high-speed camera monitoring points are added mid-span, where conductors frequently vibrate, and the sampling frequency is doubled during sunrise, when temperature gradients fluctuate dramatically. The new data collection solution is implemented through the dynamic deployment of a wireless sensor network, with some monitoring nodes employing mobile inspections to fill in blind spots left by fixed monitoring.

[0072] The trigger for iterative simulations is based on the rate of change of environmental conditions. When the monitored temperature change rate exceeds 3°C / h or the wind speed increases by 5m / s², the system automatically initiates a new round of multiphysics coupling calculations. The iterative process reuses some intermediate results from previous calculations, employing an incremental update method to improve computational efficiency. The newly generated coupled field characteristics of each iteration are analyzed for differences with the previous results, highlighting areas of significant physical quantity distribution changes.

[0073] Anti-icing solution optimization utilizes a multi-objective Pareto frontier search method. The system systematically explores the effectiveness of various parameter combinations, seeking the optimal balance between de-icing efficiency, energy costs, and equipment safety. The optimization process retains multiple non-inferior solutions for operators to choose from, each accompanied by a comprehensive technical and economic evaluation report. The final output solution is expected to reduce energy consumption by approximately 15% and equipment wear by 20% while maintaining core anti-icing performance.

[0074] Simulation test reports are generated using automated structured document layout technology. The report content is logically organized, including sections such as an overview of environmental conditions, simulation method description, results analysis, discussion of deviations, and improvement suggestions. Key data is presented interactively, allowing readers to dynamically explore interesting data details. A report version management system records the content and reasons for each modification and supports comparison of differences between different versions.

[0075] The system's self-diagnostic capabilities continuously monitor the health of the entire workflow. Metrics such as computing resource utilization, data transmission integrity, and model convergence are analyzed in real time. When localized performance degradation is detected, the system automatically initiates fault isolation and recovery procedures. This self-healing capability significantly improves reliability during extended unattended operation, which is particularly important in extreme weather conditions.

[0076] A knowledge accumulation mechanism transforms the experience gained from each iteration into system intelligence. Solutions for typical scenarios are abstracted into case studies and stored in a knowledge base, enabling rapid recall of historical best practices under similar circumstances. Parameter adjustment records form a decision tree model, helping the system make more reasonable initial choices when faced with new circumstances. This continuous learning capability enables the system's simulation and prediction capabilities to improve over time.

[0077] It should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that includes a list of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus.

[0078] While embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions, and variations may be made to these embodiments without departing from the principles and spirit of the invention, and that the scope of the invention is defined by the appended claims and their equivalents.

Claims

1. A method for simulating ice and snow environment testing based on multi-field coupling, characterized in that: The method comprises: Collect dynamic multi-field data sets in the target ice and snow environment, wherein the dynamic multi-field data sets include periodic temperature change series, relative humidity fluctuation series and real-time monitoring data of surface ice thickness; Performing multi-physics field coupling numerical simulation processing on the dynamic multi-field data set to generate coupled field distribution characteristics of the target ice and snow environment, wherein the coupled field distribution characteristics include an axial heat transfer coefficient gradient, a radial fluid velocity accumulation, and a tangential ice layer stress fluctuation coefficient; Calling a preconfigured ice and snow growth simulation model to perform ice layer evolution analysis and processing on the coupled field distribution characteristics, and generating an ice layer growth rate prediction value and a location identifier of a high-risk icing area; performing an environmental dynamic compensation correction process on the ice layer growth rate prediction value to generate a corrected ice layer growth rate prediction value, wherein the environmental dynamic compensation correction process is implemented by utilizing the correlation characteristics between the relative humidity fluctuation sequence and the heat capacity property of the material; Based on the high-risk icing area location identifier, an anti-icing treatment solution set is generated, and the anti-icing treatment solution set includes a heat source deployment strategy and a surface ice melting technology solution.

2. The ice and snow environment simulation test method based on multi-field coupling according to claim 1 is characterized in that: The performing multi-physics field coupling numerical simulation processing on the dynamic multi-field data set to generate coupling field distribution characteristics of the target ice and snow environment includes: Dividing the periodic temperature variation sequence into a plurality of temperature subsequences according to a preset time window, each temperature subsequence corresponding to a simulation period unit; For each temperature subsequence, perform the following operations: Constructing a three-dimensional ice layer deformation topology model based on the real-time monitoring data of the surface ice layer thickness, wherein the three-dimensional ice layer deformation topology model includes spatial coordinate information of thermal expansion field distribution, fluid velocity field distribution, and stress field distribution; Performing field coupling interaction processing on the three-dimensional ice layer deformation topology model and the temperature subsequence to generate a coupled field simulation output of the current simulation period unit; Performing cumulative integration calculation on the coupled field simulation output of multiple consecutive simulation period units to obtain the axial heat transfer coefficient gradient, radial fluid velocity accumulation and tangential ice layer stress fluctuation coefficient; wherein, The axial heat transfer coefficient gradient is the maximum rate of change of the thermal expansion field along the axial direction. The radial fluid velocity accumulation is the total integral of the fluid velocity field in the normal direction of the ice layer contact surface. The tangential ice layer stress fluctuation coefficient is the coefficient of variation of the stress field distribution over time.

3. The ice and snow environment simulation test method based on multi-field coupling according to claim 2 is characterized in that: The performing field coupling interaction processing on the three-dimensional ice layer deformation topology model and the temperature subsequence to generate a coupled field simulation output of the current simulation period unit includes: Based on the corresponding relationship between the heat conduction law and the temperature subsequence, a thermal field-temperature balance equation is established, and an initial distribution function of the heat conduction component is obtained by solving the thermal field-temperature balance equation; Based on the correlation between fluid dynamics principles and ice deformation data, a fluid-ice coupling calculation framework is constructed, which integrates dynamic adjustment parameters of fluid viscosity coefficient and ice elastic modulus; Combining the surface stress change rate with the friction effect factor, an iterative ice growth prediction process is implemented, which includes a feedback calibration mechanism for the stress increment and ice deformation. A spatial field fusion operation is performed on the output results of the initial distribution function, the fluid-ice layer coupling calculation framework, and the ice layer growth iterative prediction process to generate a three-dimensional coupled field distribution data set containing heat conduction, fluid velocity, and ice layer stress components.

4. The ice and snow environment simulation test method based on multi-field coupling according to claim 1 is characterized in that: The calling of a preconfigured ice and snow growth simulation model to perform ice layer evolution analysis and processing on the coupled field distribution characteristics to generate an ice layer growth rate prediction value and a high-risk icing area location identifier includes: Inputting the axial heat transfer coefficient gradient into the initial characteristic analysis layer of the ice and snow growth simulation model, and determining the distribution coordinates of the heat concentration area and the heat amplitude change trajectory through the heat concentration factor calculation unit; Inputting the radial fluid velocity accumulation into the secondary characteristic analysis layer of the ice and snow growth simulation model, performing fluid scour damage accumulation calculation, and generating the melting probability and crack growth rate prediction value of the ice surface; Inputting the tangential ice layer stress fluctuation coefficient into the final characteristic analysis layer of the ice and snow growth simulation model, and calculating the material wear depth of the friction surface and the ice layer roughness evolution data based on the surface ice layer degradation model; The thermal amplitude variation trajectory, the melting probability, and the material wear depth are integrated to generate a comprehensive evolution index of the target ice and snow environment, and a predicted value of the ice layer growth rate is determined based on a comparison result of the comprehensive evolution index with a preset threshold library; Based on the spatial field superposition results of the distribution coordinates, the crack growth rate prediction value and the ice layer roughness evolution data, the geometric position information of the heat concentration area, the crack path and the high-risk wear area is identified.

5. The ice and snow environment simulation test method based on multi-field coupling according to claim 1 is characterized in that: The performing of environmental dynamic compensation correction processing on the ice layer growth rate prediction value to generate a corrected ice layer growth rate prediction value includes: Extracting the extreme humidity values ​​and humidity change frequencies in the relative humidity fluctuation sequence, and calculating the dynamic correction amount of the material heat capacity property as the humidity fluctuates; performing a thermal stress compensation calculation on the axial heat transfer coefficient gradient according to the dynamic correction amount to generate a compensated heat transfer coefficient gradient; Based on the correlation between the humidity change frequency and the phase change characteristics of the material, performing phase change damage correction processing on the radial fluid velocity accumulation to generate a corrected fluid velocity accumulation; performing surface strength adaptation adjustment processing on the tangential ice layer stress fluctuation coefficient according to the material thermal conductivity variation data under extreme humidity to generate an adjusted ice layer stress fluctuation coefficient; The compensated thermal conductivity gradient, the corrected fluid velocity accumulation and the adjusted ice layer stress fluctuation coefficient are input into the ice and snow growth simulation model for recalculation to generate an ice layer growth rate prediction value including environmental factor compensation.

6. The ice and snow environment simulation test method based on multi-field coupling according to claim 5 is characterized in that: The performing of thermal stress compensation calculation on the axial heat transfer coefficient gradient according to the dynamic correction amount to generate a compensated heat transfer coefficient gradient includes: Obtaining the initial heat capacity property of the target ice and snow environment under reference humidity conditions and the dynamic correction amount, and establishing a heat capacity-humidity correlation function; Calculating an axial thermal stress compensation value according to the heat capacity-humidity correlation function, wherein the axial thermal stress compensation value is the product of the humidity change amplitude and the heat capacity change amplitude; superimposing the axial thermal stress compensation value into the calculation process of the axial thermal conductivity gradient to generate an intermediate thermal conductivity gradient value including thermal stress compensation; A thermal relaxation effect compensation process is performed on the intermediate value of the thermal conductivity gradient, wherein the thermal relaxation effect compensation process is implemented by using a multiplication factor of a material thermal relaxation curve and a humidity retention time.

7. The ice and snow environment simulation test method based on multi-field coupling according to claim 4 is characterized in that: The generating of an anti-icing treatment solution set based on the high-risk icing area location identifier includes: Calculating an optimal heat source deployment path based on the identification information of the heat concentration area, wherein the optimal heat source deployment path is achieved by adjusting the power distribution ratio of adjacent heat sources; Constructing a surface ice melting technical solution based on the identification information of the crack path, wherein the surface ice melting technical solution includes the selection of a microwave heating area and the optimization configuration of energy parameters; generating a lubricant application strategy based on the identification information of the high-risk wear area, wherein the lubricant application strategy dynamically adjusts the application frequency and application thickness according to the wear rate prediction value; The optimal heat source deployment path, the surface ice melting technical solution and the lubricant application strategy are prioritized to generate an anti-icing solution set including execution timing and technical parameters.

8. The ice and snow environment simulation test method based on multi-field coupling according to claim 7 is characterized in that: The surface ice melting technology solution includes: Extracting geometric characteristic data of the crack path and calculating the path curvature radius and expansion direction angle; selecting a coverage density parameter of the microwave heating spot according to the curvature radius, wherein the coverage density parameter is in a nonlinear inverse proportional relationship with the curvature radius; Adjusting the direction of application of microwave heating energy based on the expansion direction angle so that the energy application direction forms a preset angle with the crack path direction; Dynamically adjust the duration of microwave heating pulses based on the test data of the material's surface thermal conductivity to ensure that the heating energy is below the critical value of the material's melting point; Generate an ice melting parameter configuration table including coverage density parameters, application direction and pulse duration.

9. The ice and snow environment simulation test method based on multi-field coupling according to claim 1 is characterized in that: The method further comprises: Collecting actual ice thickness and crack propagation length data of the target ice and snow environment within a preset verification period; performing deviation analysis on the actual ice thickness and the predicted ice thickness to generate a thickness error correction coefficient; Performing a time series comparison process on the crack propagation length and the predicted propagation speed to generate a length error correction coefficient; Adjusting the parameter weights of the ice and snow growth simulation model according to the thickness error correction coefficient and the length error correction coefficient to generate an optimized ice and snow growth simulation model; The optimized ice and snow growth simulation model is applied to ice and snow environment simulation tasks in subsequent test batches.

10. The ice and snow environment simulation test method based on multi-field coupling according to claim 1 is characterized in that: After generating the anti-icing treatment solution set, the method further includes: Inputting the anti-icing treatment plan set into the environmental simulation control system to perform real-time scheduling processing; updating the acquisition parameters of the dynamic multi-field data set according to the scheduling result feedback; Restarting the multi-physics field coupling numerical simulation process using the updated acquisition parameters to generate an iterative coupling field distribution feature; The anti-icing treatment scheme set is optimized based on the iterative coupling field distribution characteristics, and a final simulation test report is output.

Images