Simulation methods, devices, terminals and media for icing morphology of wind turbine blades
By using the Lagrange method and distortion point removal technology in the simulation of icing on wind turbine blades, combined with moving average and linear fitting, the problem of low accuracy in icing simulation was solved, achieving more accurate simulation of icing morphology and supporting more effective anti-icing design.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- ELECTRIC POWER RES INST CHINA SOUTHERN POWER GRID CO LTD
- Filing Date
- 2023-07-21
- Publication Date
- 2026-06-30
AI Technical Summary
Existing simulation technologies for icing on wind turbine blades have low accuracy, resulting in significant discrepancies between icing data and actual conditions, which affects the effectiveness of anti-icing designs.
The Lagrange method was used to simulate the structure of wind turbine blades. The coordinates of local water droplet collision points were calculated, and distortion points were removed. Distortion points were identified and removed by changing the slope angle. The icing morphology curve was optimized by combining moving average processing and linear fitting to improve the simulation accuracy.
The accuracy of wind turbine blade icing simulation has been improved, and simulation results of icing morphology that are closer to the actual situation have been obtained, supporting more effective anti-icing design.
Smart Images
Figure CN116738760B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of wind power generation simulation technology, and in particular to a method, device, terminal and medium for simulating the icing morphology of wind turbine blades. Background Technology
[0002] As a clean and renewable energy source, wind energy is gradually being accepted and widely applied by more and more countries. With the increasing maturity of wind power generation technology, many developed countries have increased their investment in the construction of the wind power industry, and my country's wind power generation industry is also growing rapidly.
[0003] In recent years, the wind power industry in Sichuan, Chongqing, Yunnan and other regions of my country has been increasing year by year. These regions are located in the southwest of my country, with a humid climate and prominent plateau and mountainous terrain. Winter temperatures are around 0°C, which provides environmental conditions for icing. The problem of wind turbine blade icing has gradually become prominent. The average annual power generation loss caused by wind turbine blade icing can be as high as 30%. Secondly, wind turbine blade icing can change the structure of the wind turbine, making the technical components of the wind turbine brittle and reducing the service life of the wind turbine. Furthermore, when the temperature rises, the ice on the wind turbine blades will fall off and be thrown out under the action of centrifugal force, which seriously endangers the personal safety of wind farm workers.
[0004] Currently, the most common method for preventing and controlling icing on wind turbine blades is to simulate the growth pattern and rate of icing based on environmental conditions. This provides a technical reference for the anti-icing design of wind turbine blades and has important guiding significance for the anti-icing and de-icing work of wind turbine blades. However, in practical applications, there are significant errors between the icing data calculated by simulation and the actual icing situation, which restricts the further development of simulation technology in the prevention and control of icing on wind turbine blades. Summary of the Invention
[0005] This application provides a method, device, terminal, and medium for simulating the icing morphology of wind turbine blades, which solves the technical problem of low accuracy in existing wind turbine blade icing simulation technologies.
[0006] To address the aforementioned technical problems, the first aspect of this application provides a method for simulating the icing morphology of wind turbine blades, comprising:
[0007] Based on the structural information of the wind turbine blade, simulation calculations were performed using the Lagrange method to obtain the set of coordinates of local water droplet collision points on the surface of the wind turbine blade.
[0008] Based on the set of local water droplet collision coordinates, calculate the local water droplet collision coefficient of each adjacent water droplet collision point;
[0009] Based on the coordinates of the water droplet collision points, the local water droplet collision coefficients of adjacent collision points are connected to obtain a collision coefficient curve. By comparing the slope angle changes of adjacent connections, if the slope angle changes exceed a preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connections are determined as distortion points, and the distortion points are removed.
[0010] Based on the collision coefficient curve with distortion points removed, the first icing morphology curve and the first icing thickness corresponding to the first icing morphology curve are calculated using the normal local water droplet collision coefficient contained in the collision coefficient curve, so as to obtain the icing morphology simulation result according to the first icing morphology curve and the first icing thickness.
[0011] Preferably, after obtaining the set of coordinates of local water droplet collision points on the surface of the wind turbine blade, the method further includes:
[0012] If the difference between the vertical coordinates of adjacent collision points does not meet the preset coordinate difference threshold, the water droplet collision point coordinates are interpolated using linear interpolation to obtain the interpolated set of local water droplet collision point coordinates.
[0013] Preferably, it further includes:
[0014] Based on the first icing morphology curve, and combined with the moving average processing method, the collision coefficients of several local water droplets at the front and rear ends of the first icing morphology curve are linearly fitted to obtain two zero points of water droplet collision coefficients through fitting.
[0015] Based on the zero point of the water droplet collision coefficient and the normal local water droplet collision coefficient, the second icing morphology curve and the second icing thickness corresponding to the second icing morphology curve are calculated, so as to obtain the icing morphology simulation result based on the second icing morphology curve and the second icing thickness.
[0016] Preferably, the step of linearly fitting several local water droplet collision coefficients at the front and rear ends of the first icing morphology curve based on the first icing morphology curve and using a moving average processing method to obtain the zero points of the two water droplet collision coefficients specifically includes:
[0017] Based on the collision coefficient curve with distortion points removed, the normal local water droplet collision coefficient of the collision coefficient curve is processed by moving average through a preset moving average window.
[0018] Based on the collision coefficient curve after moving average processing, linear fitting is performed on several local water droplet collision coefficients at the front and back ends of the collision coefficient curve to obtain two zero points of water droplet collision coefficient through fitting.
[0019] Preferably, the formula for calculating the number of local droplet collision coefficients for linear fitting is:
[0020] m = 20% N;
[0021] In the formula, m is the number of local water droplet collision coefficients that are linearly fitted at the front or back end of the collision coefficient curve, and N is the number of water droplet collision points.
[0022] Preferably, after calculating the second icing morphology curve and the second icing thickness corresponding to the second icing morphology curve based on the zero point of the water droplet collision coefficient and the normal local water droplet collision coefficient, the method further includes:
[0023] Based on the error comparison results between the first icing thickness and the second icing thickness, the simulation results of the icing morphology are optimized or output according to the error comparison results.
[0024] Preferably, the step of optimizing or outputting the icing morphology simulation result based on the error comparison result between the first icing thickness and the second icing thickness specifically includes:
[0025] Based on the error comparison results of the first ice thickness and the second ice thickness, if the average error of the first ice thickness and the second ice thickness exceeds a preset error threshold, the parameters of the moving average window are modified, and then the normal local water droplet collision coefficient of the collision coefficient curve is returned to be processed by moving average. If it does not exceed the threshold, the ice morphology simulation result is obtained based on the second ice morphology curve and the second ice thickness.
[0026] Meanwhile, the second aspect of this application also provides a simulation device for the icing morphology of wind turbine blades, including:
[0027] The collision point coordinate acquisition unit is used to perform simulation calculations based on the wind turbine blade structure information and the Lagrange method to obtain the set of local water droplet collision point coordinates on the surface of the wind turbine blade.
[0028] The water droplet collision coefficient calculation unit is used to calculate the local water droplet collision coefficient of each adjacent water droplet collision point based on the local water droplet collision coordinate set.
[0029] The distortion point processing unit is used to connect the local water droplet collision coefficients of adjacent collision points according to the coordinates of the water droplet collision points to obtain a collision coefficient curve. By comparing the slope angle changes of adjacent connections, if the slope angle changes exceed a preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connections are determined as distortion points, and the distortion points are removed.
[0030] The icing data calculation unit is used to calculate the first icing morphology curve and the first icing thickness corresponding to the first icing morphology curve based on the collision coefficient curve after removing distortion points and the normal local water droplet collision coefficient contained in the collision coefficient curve, so as to obtain the icing morphology simulation result according to the first icing morphology curve and the first icing thickness.
[0031] A third aspect of this application provides a simulation terminal for icing morphology of wind turbine blades, including: a memory and a processor;
[0032] The memory stores computer programs that can run on the processor;
[0033] When the processor executes the computer program, it implements the wind turbine blade icing morphology simulation method as provided in the first aspect of this application.
[0034] The fourth aspect of this application provides a computer-readable storage medium storing computer-executable instructions, which, when executed by a control processor, implement the wind turbine blade icing morphology simulation method provided in the first aspect of this application.
[0035] As can be seen from the above technical solutions, the embodiments of this application have the following advantages:
[0036] The solution provided in this application uses the Lagrange method to perform preliminary simulation based on the structure of wind turbine blades to obtain the set of coordinates of local water droplet collision points on the blade surface. Then, along the blade surface, the local water droplet collision coefficients of the collision points are calculated. Based on the coordinates of the water droplet collision points, the local water droplet collision coefficients of adjacent collision points are connected to obtain the collision coefficient curve. By comparing the slope angle changes of adjacent connections, if the slope angle change exceeds a preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connections are identified as distortion points. These distortion points are removed to avoid the impact of the existence of distortion points on the accuracy of the icing simulation results, thereby solving the technical problem of low accuracy in existing wind turbine blade icing simulation technology. Attached Figure Description
[0037] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0038] Figure 1 This is a flowchart illustrating an embodiment of a method for simulating the icing morphology of wind turbine blades provided in this application.
[0039] Figure 2 This is a flowchart illustrating another embodiment of the wind turbine blade icing morphology simulation method provided in this application.
[0040] Figure 3 The image shows the trajectory tracking of water droplets on the outside of the wind turbine blade in a simulation method for icing morphology of wind turbine blades provided in this application.
[0041] Figure 4 This is a schematic diagram of the distortion points of the local collision coefficient of water droplets on the surface of a wind turbine blade in a simulation method for icing morphology of wind turbine blades provided in this application.
[0042] Figure 5 This is a schematic diagram illustrating the determination of the zero-value point of the local collision coefficient of water droplets on the surface of a wind turbine blade in a simulation method for icing morphology of wind turbine blades provided in this application.
[0043] Figure 6 Comparison images before and after the reconstruction of the icing morphology simulation of the wind turbine blade surface, which is provided in this application for the simulation method of icing morphology simulation of wind turbine blade.
[0044] Figure 7 This is a schematic diagram of an embodiment of a wind turbine blade icing morphology simulation device provided in this application.
[0045] Figure 8 A schematic diagram of the structure of a simulation terminal for icing morphology of wind turbine blades provided in this application. Detailed Implementation
[0046] Wind turbine blades have airfoil structures of varying sizes from hub to tip. Water droplets at different locations on the blade freeze at different rates, resulting in icing structures of varying thicknesses. When simulating icing on blade cross-sections, the icing thickness needs to be calculated based on local droplet collision and freezing coefficients. The Lagrange method is a commonly used method for simulating droplet motion and collisions. However, due to the dispersed nature of droplet collision positions, the calculated local collision coefficients often exhibit distortion points, causing significant difficulties for subsequent numerical calculations of icing morphology and mesh updates. This leads to deviations in the calculation results, ultimately resulting in the low accuracy of existing wind turbine blade icing simulation technologies.
[0047] In view of this, in order to solve the above-mentioned technical problems, this application provides a method, device, terminal and medium for simulating the icing morphology of wind turbine blades, which is used to solve the technical problem of low accuracy of existing wind turbine blade icing simulation technology.
[0048] To make the inventive objectives, features, and advantages of this application more apparent and understandable, the technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the embodiments described below are only some embodiments of this application, and not all embodiments. Based on the embodiments in this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.
[0049] First, a detailed description of the basic and further embodiments of the wind turbine blade icing morphology simulation method provided in this application is as follows:
[0050] Please see Figure 1 This embodiment provides a method for simulating the icing morphology of wind turbine blades, including:
[0051] Step 101: Based on the structural information of the wind turbine blade, perform simulation calculations using the Lagrange method to obtain the set of coordinates of local water droplet collision points on the surface of the wind turbine blade.
[0052] It should be noted that, firstly, based on the blade structure information of the wind turbine generator to be simulated, a conventional Lagrange multiplication algorithm is used for simulation calculations to obtain the local water droplet collision points on the surface of the wind turbine blades. Based on the coordinates of these local water droplet collision points, a corresponding coordinate set A is formed, as shown in the following example. Figure 3 As shown.
[0053] Step 102: Calculate the local droplet collision coefficient at each adjacent droplet collision point based on the local droplet collision coordinate set;
[0054] Based on the set of local water droplet collision coordinates obtained in the previous step, the local water droplet collision coefficients of each adjacent water droplet collision point are calculated respectively. In this embodiment, the local water droplet collision coefficient can be calculated by using the ratio of the initial position interval of the water droplet to the collision position interval of the water droplet to calculate the local water droplet collision coefficient on the blade surface. The local water droplet collision coefficient β1 can be calculated according to equation (1).
[0055]
[0056] Where dY is the initial distance between adjacent water droplets, in meters; and dL is the relative distance between adjacent water droplets after they collide with the fan blades, in meters.
[0057] Step 103: Based on the coordinates of the water droplet collision points, connect the local water droplet collision coefficients of adjacent collision points to obtain the collision coefficient curve. By comparing the slope angle changes of adjacent connections, if the slope angle changes exceed the preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connections are determined as distortion points, and the distortion points are removed.
[0058] Step 104: Based on the collision coefficient curve with distortion points removed, and based on the normal local water droplet collision coefficient contained in the collision coefficient curve, calculate the first icing morphology curve and the first icing thickness corresponding to the first icing morphology curve, so as to obtain the icing morphology simulation result according to the first icing morphology curve and the first icing thickness.
[0059] like Figure 4 As shown, based on the coordinates of the water droplet collision points, the local water droplet collision coefficients of adjacent collision points are connected to obtain the collision coefficient curve. By comparing the slope angle changes of adjacent connections, if the slope angle change exceeds a preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connections are identified as distortion points. More specifically, the water droplet collision coefficient β1 along the surface of the wind turbine blade is iterated one by one, and the slope angle of the line connecting two adjacent collision points β1 is compared with the slope angle of the line connecting the previous two water droplet collision points β1. If the slope angle change of β1 exceeds 90°, the position is identified as an abnormal change in the local water droplet collision coefficient and marked. For example, three adjacent collision points a(x1, y1)~β 11 b(x2, y2)~β 21 c(x3, y3)~β 31 The slope change value of β1 at point b is dθ. If dθ ≥ 90°, then point b is a distorted point; otherwise, point b is a normal point. This step is repeated to find the remaining distorted points and remove them to obtain a collision coefficient curve without distorted points. The formula for calculating the slope change value of β1 mentioned in this embodiment can be referred to as follows:
[0060]
[0061] Finally, using the processed collision coefficient curve, the first icing morphology curve and the corresponding first icing thickness are calculated according to conventional methods, so as to obtain the icing morphology simulation results based on the first icing morphology curve and the first icing thickness.
[0062] The above is a detailed description of a basic embodiment of a wind turbine blade icing morphology simulation method provided in this application. The scheme provided in this embodiment uses the Lagrange method to perform preliminary simulation based on the wind turbine blade structure to obtain the set of coordinates of local water droplet collision points on the blade surface. Then, along the blade surface, the local water droplet collision coefficient of each water droplet collision point is calculated. Based on the coordinates of the water droplet collision points, the local water droplet collision coefficients of adjacent collision points are connected to obtain a collision coefficient curve. By comparing the slope angle changes of adjacent connections, if the slope angle change exceeds a preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connections are determined as distortion points. These distortion points are removed to avoid the influence of the existence of distortion points on the accuracy of the icing simulation results. This allows for the acquisition of continuous, smooth, and more realistic wind turbine blade icing morphology simulation results, solving the technical problem of low accuracy in existing wind turbine blade icing simulation technologies.
[0063] Based on the detailed description of the above basic embodiments, the following is a detailed description of a further embodiment of the wind turbine blade icing morphology simulation method provided in this application, as follows:
[0064] Please see Figure 2 and Figure 3 Furthermore, after obtaining the set of coordinates of local water droplet collision points on the surface of the wind turbine blade in step 101, the process may further include:
[0065] Step 1001: Based on the water droplet collision point coordinates in the local water droplet collision point coordinate set, if the difference in the vertical axis coordinates of adjacent collision points does not meet the preset coordinate difference threshold, then the water droplet collision point coordinates are interpolated using the linear interpolation method to obtain the interpolated local water droplet collision point coordinate set.
[0066] It should be noted that, based on the set of local water droplet collision point coordinates A obtained in step 101, the water droplet collision point coordinates are interpolated using the linear interpolation method to obtain the set of interpolated water droplet collision point coordinates B. After interpolation, the difference between the vertical axis coordinates of two adjacent collision points is not less than 1 / 30 of the blade thickness h.
[0067] Furthermore, after step 104, the following is also included:
[0068] Step 105: Based on the first icing morphology curve, and combined with the moving average processing method, perform linear fitting on several local water droplet collision coefficients at the front and rear ends of the first icing morphology curve, so as to obtain two zero points of water droplet collision coefficient through fitting.
[0069] Step 106: Calculate the second icing morphology curve and the second icing thickness corresponding to the second icing morphology curve based on the zero point of the water droplet collision coefficient and the normal local water droplet collision coefficient, so as to obtain the icing morphology simulation results based on the second icing morphology curve and the second icing thickness.
[0070] More specifically, step 105 described above may include the following steps:
[0071] Based on the first icing morphology curve, the normal local water droplet collision coefficient of the collision coefficient curve is processed by moving average through a preset moving average window.
[0072] Based on the collision coefficient curve after moving average processing, linear fitting is performed on the collision coefficients of several local water droplets at the front and back ends of the collision coefficient curve to obtain the zero points of the two water droplet collision coefficients through fitting.
[0073] It should be noted that, as Figure 5 As shown, in steps 105 to 106 of this embodiment, based on the first icing morphology curve obtained in step 104, the moving average subset window size is set to w, and the local water droplet collision coefficient value corresponding to the first icing morphology curve is subjected to reasonable moving average processing. The first and last m β1 value points of the first and last m positions of the local water droplet collision coefficient value after the averaging process are extracted and linearly fitted to determine the positions P1 and P2 where the local collision coefficients at the first and last ends are zero. Based on the obtained P1 and P2, combined with the remaining normal local water droplet collision coefficients, the second icing morphology curve and the second icing thickness corresponding to the second icing morphology curve are calculated. According to the second icing morphology curve and the second icing thickness, the value of m is determined according to the total number of collision points, and can be taken as m = 20%N.
[0074] Furthermore, step 106 includes the following:
[0075] Step 107: Based on the error comparison results of the first icing thickness and the second icing thickness, optimize or output the icing morphology simulation results.
[0076] Furthermore, step 107 specifically includes:
[0077] Based on the error comparison results of the first ice thickness and the second ice thickness, if the average error between the first ice thickness and the second ice thickness exceeds the preset error threshold, the parameters of the moving average window are modified, and then the normal local water droplet collision coefficient of the collision coefficient curve is returned to be processed by moving average. If it does not exceed the threshold, the ice morphology simulation result is obtained based on the second ice morphology curve and the second ice thickness.
[0078] It should be noted that the thickness data of the first icing shape Shape1 and the second icing shape Shape2, calculated using the local droplet collision coefficient values obtained in steps 104 and 106, are compared. Figure 6 As shown, the average error Error corresponding to the ice thickness of the two ice cladding patterns is calculated;
[0079] Determine whether the comparison result meets the minimum error ε requirement. If not, adjust the size parameter w of the moving average window, and then return to step 1051. Based on the adjusted moving average window, re-execute steps 1051 to 107 until the obtained comparison result meets the requirement.
[0080] The above is a detailed description of a further embodiment of the wind turbine blade icing morphology simulation method provided in this application. The following is a detailed description of an embodiment of the wind turbine blade icing morphology simulation device provided in this application.
[0081] Please see Figure 7 The third embodiment of this application provides a simulation device for the icing morphology of wind turbine blades, including:
[0082] The collision point coordinate acquisition unit 201 is used to perform simulation calculations based on the wind turbine blade structure information and the Lagrange method to obtain the set of local water droplet collision point coordinates on the surface of the wind turbine blade.
[0083] The water droplet collision coefficient calculation unit 202 is used to calculate the local water droplet collision coefficient of each adjacent water droplet collision point based on the local water droplet collision coordinate set.
[0084] The distortion point processing unit 203 is used to connect the local water droplet collision coefficients of adjacent collision points according to the coordinates of the water droplet collision points to obtain the collision coefficient curve. By comparing the slope angle changes of adjacent connections, if the slope angle changes exceed the preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connections are determined as distortion points and the distortion points are removed.
[0085] The icing data calculation unit 204 is used to calculate the first icing morphology curve and the first icing thickness corresponding to the first icing morphology curve based on the collision coefficient curve after removing the distortion points and the normal local water droplet collision coefficient contained in the collision coefficient curve, so as to obtain the icing morphology simulation result according to the first icing morphology curve and the first icing thickness.
[0086] The above is a detailed description of an embodiment of a wind turbine blade icing morphology simulation device provided in this application. The following is a detailed description of an embodiment of a wind turbine blade icing morphology simulation terminal and a computer-readable storage medium provided in this application.
[0087] Please see Figure 8 The fourth embodiment of this application provides a simulation terminal embodiment for the icing morphology of wind turbine blades. The types of terminals include, but are not limited to, personal computers, industrial computers, servers, and embedded intelligent devices. The simulation terminal for the icing morphology of wind turbine blades provided in this embodiment includes: a memory 33 and a processor 31. The memory and the processor can be connected via a bus 34.
[0088] The memory 33 stores a computer program that can run on the processor 31;
[0089] When processor 31 executes a computer program, it implements the wind turbine blade icing morphology simulation method provided in the first or second embodiment of this application.
[0090] The fifth embodiment of this application provides a computer-readable storage medium storing computer-executable instructions. When the computer-executable instructions are executed by a control processor, they implement the wind turbine blade icing morphology simulation method provided in the first or second embodiment of this application.
[0091] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the terminals, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.
[0092] In the several embodiments provided in this application, it should be understood that the disclosed terminals, devices, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative; for instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection between devices or units through some interfaces, and may be electrical, mechanical, or other forms.
[0093] The terms “first,” “second,” “third,” “fourth,” etc. (if present) in the specification and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that embodiments of the application described herein can be implemented, for example, in orders other than those illustrated or described herein. Furthermore, the terms “comprising” and “having,” and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.
[0094] It should be understood that in this application, "at least one (item)" means one or more, and "more than" means two or more. "And / or" is used to describe the relationship between related objects, indicating that three relationships can exist. For example, "A and / or B" can represent three cases: only A exists, only B exists, and both A and B exist simultaneously, where A and B can be singular or plural. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. "At least one (item) of the following" or similar expressions refer to any combination of these items, including any combination of single or plural items. For example, at least one (item) of a, b, or c can represent: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, and c can be single or multiple.
[0095] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.
[0096] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.
[0097] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.
[0098] The above-described embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application.
Claims
1. A method for simulating the icing morphology of wind turbine blades, characterized in that, include: Based on the structural information of the wind turbine blade, simulation calculations were performed using the Lagrange method to obtain the set of coordinates of local water droplet collision points on the surface of the wind turbine blade. Based on the water droplet collision point coordinates in the local water droplet collision point coordinate set, if the difference in the vertical axis coordinates of adjacent collision points does not meet the preset coordinate difference threshold, then the water droplet collision point coordinates are interpolated using linear interpolation to obtain the interpolated local water droplet collision point coordinate set. Based on the set of local water droplet collision point coordinates, calculate the local water droplet collision coefficient for each adjacent water droplet collision point; Based on the coordinates of the water droplet collision point, the local water droplet collision coefficients of adjacent collision points are connected to obtain the collision coefficient curve. By comparing the slope angle changes of adjacent connecting lines, if the slope angle changes exceed a preset angle threshold, the coordinates of the water droplet collision point corresponding to the adjacent connecting line are determined as distortion points, and the distortion points are removed. Based on the collision coefficient curve with distortion points removed, the first icing morphology curve and the corresponding first icing thickness are calculated using the normal local droplet collision coefficients contained in the collision coefficient curve. Based on the first icing morphology curve and the first icing thickness, an icing morphology simulation result is obtained. Alternatively, based on the first icing morphology curve, and using a moving average processing method, linear fitting is performed on several local droplet collision coefficients at the front and rear ends of the first icing morphology curve to obtain two zero points for droplet collision coefficients. Based on the zero points for droplet collision coefficients and the normal local droplet collision coefficients, a second icing morphology curve and the corresponding second icing thickness are calculated. Based on the second icing morphology curve and the second icing thickness, an icing morphology simulation result is obtained.
2. The method for simulating the icing morphology of wind turbine blades according to claim 1, characterized in that, The step of linearly fitting several local water droplet collision coefficients at the front and rear ends of the first icing morphology curve, based on the first icing morphology curve and using a moving average processing method, to obtain the zero points of the two water droplet collision coefficients through fitting, specifically includes: Based on the collision coefficient curve with distortion points removed, the normal local water droplet collision coefficient of the collision coefficient curve is processed by moving average through a preset moving average window. Based on the collision coefficient curve after moving average processing, linear fitting is performed on several local water droplet collision coefficients at the front and back ends of the collision coefficient curve to obtain two zero points of water droplet collision coefficient through fitting.
3. The method for simulating the icing morphology of wind turbine blades according to claim 2, characterized in that, The formula for calculating the number of local droplet collision coefficients for linear fitting is: m = 20%N; In the formula, m is the number of local water droplet collision coefficients that are linearly fitted at the front or back end of the collision coefficient curve, and N is the number of water droplet collision points.
4. The method for simulating the icing morphology of wind turbine blades according to claim 2, characterized in that, Based on the zero point of the water droplet collision coefficient and the normal local water droplet collision coefficient, the calculation of the second icing morphology curve and the second icing thickness corresponding to the second icing morphology curve further includes: Based on the error comparison results between the first icing thickness and the second icing thickness, the simulation results of the icing morphology are optimized or output according to the error comparison results.
5. The method for simulating the icing morphology of wind turbine blades according to claim 4, characterized in that, The step of optimizing or outputting the icing morphology simulation result based on the error comparison result between the first icing thickness and the second icing thickness specifically includes: Based on the error comparison results of the first ice thickness and the second ice thickness, if the average error of the first ice thickness and the second ice thickness exceeds a preset error threshold, the parameters of the moving average window are modified, and then the normal local water droplet collision coefficient of the collision coefficient curve is returned to be processed by moving average. If it does not exceed the threshold, the ice morphology simulation result is obtained based on the second ice morphology curve and the second ice thickness.
6. A simulation device for icing morphology of wind turbine blades, characterized in that, include: The collision point coordinate acquisition unit is used to perform simulation calculations based on the wind turbine blade structure information and the Lagrange method to obtain the set of local water droplet collision point coordinates on the surface of the wind turbine blade. The interpolation processing unit is used to interpolate the water droplet collision point coordinates in the local water droplet collision point coordinate set according to the water droplet collision point coordinates in the local water droplet collision point coordinate set. If the difference between the vertical axis coordinates of adjacent collision points does not meet the preset coordinate difference threshold, the unit will interpolate the water droplet collision point coordinates using a linear interpolation method to obtain the interpolated local water droplet collision point coordinate set. The water droplet collision coefficient calculation unit is used to calculate the local water droplet collision coefficient of each adjacent water droplet collision point based on the set of local water droplet collision point coordinates. The distortion point processing unit is used to connect the local water droplet collision coefficients of adjacent collision points according to the coordinates of the water droplet collision points to obtain a collision coefficient curve. By comparing the slope angle changes of adjacent connecting lines, if the slope angle changes exceed a preset angle threshold, the coordinates of the water droplet collision points corresponding to the adjacent connecting lines are determined as distortion points, and the distortion points are removed. The icing data calculation unit is used to calculate a first icing morphology curve and a first icing thickness corresponding to the first icing morphology curve based on a collision coefficient curve with distortion points removed, using the normal local water droplet collision coefficients included in the collision coefficient curve, so as to obtain an icing morphology simulation result based on the first icing morphology curve and the first icing thickness; and / or, based on the first icing morphology curve, combined with a moving average processing method, to perform linear fitting on several local water droplet collision coefficients at the front and rear ends of the first icing morphology curve respectively, so as to obtain two water droplet collision coefficient zero points through fitting; and to calculate a second icing morphology curve and a second icing thickness corresponding to the second icing morphology curve based on the water droplet collision coefficient zero points and the normal local water droplet collision coefficients, so as to obtain an icing morphology simulation result based on the second icing morphology curve and the second icing thickness.
7. A simulation terminal for icing morphology of wind turbine blades, characterized in that, include: Memory and processor; The memory stores computer programs that can run on the processor; When the processor executes the computer program, it implements the wind turbine blade icing morphology simulation method as described in any one of claims 1 to 5.
8. A computer-readable storage medium storing computer-executable instructions, characterized in that, When the computer-executable instructions are executed by the control processor, the method for simulating the icing morphology of wind turbine blades as described in any one of claims 1 to 5 is implemented.