CN113719428B - Method and system for evaluating service life of blade of horizontal axis wind generating set - Google Patents
Method and system for evaluating service life of blade of horizontal axis wind generating set Download PDFInfo
- Publication number
- CN113719428B CN113719428B CN202111007759.9A CN202111007759A CN113719428B CN 113719428 B CN113719428 B CN 113719428B CN 202111007759 A CN202111007759 A CN 202111007759A CN 113719428 B CN113719428 B CN 113719428B
- Authority
- CN
- China
- Prior art keywords
- blade
- thermal insulation
- front edge
- protective coating
- coating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Wind Motors (AREA)
Abstract
The invention discloses a method and a system for evaluating the service life of a blade of a horizontal axis wind generating set, which are used for obtaining a thermal insulation protective coating damage expansion rule curved surface, namely a function relation of the volume reduction rate of a test piece thermal insulation protective coating along with experiment time; obtaining a service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade, namely a function relation of the volume reduction rate of the thermal insulation protection coating on the front edge of the blade along with the spanwise position and the service time of the blade; according to the volume reduction rate of the thermal insulation protection coating at the front edge of the blade, calibrating a curved surface by using the service life of the thermal insulation protection coating at the front edge of the blade, and calculating to obtain the service time of the thermal insulation protection coating coated at the current time on the front edge of the section where the corresponding spanwise position of the blade is located; the service time of a certain spanwise position of the blade is the sum of the service time of each previous coating of the heat-insulating protective coating on the front edge of the section where the blade is located, and the service life of the blade is the average value of the service time of each spanwise position of the blade. The sensor does not need to be installed on the blade, and the service life of the blade can be obtained.
Description
Technical Field
The invention belongs to the field of wind generating sets, and relates to a method and a system for evaluating the service life of a blade of a horizontal axis wind generating set.
Background
The service life evaluation work of the wind generating set part generally comprises the following steps: (1) Determining the potential failure type and the corresponding dangerous position of the component according to the design evaluation report and the type test report; (2) Selecting a main characteristic load causing component failure and a corresponding measurable characteristic variable according to the design evaluation process and the result; (3) Arranging a sensor at a dangerous position of the component to realize the function of acquiring data of the characteristic variables; (4) Establishing a numerical function relation between the main characteristic load and the characteristic variable through a load calibration test; (5) Carrying out real-time measurement on the characteristic variables, and recording and storing data acquisition results; (6) Deducing a main characteristic load time sequence history applied to the component by using a numerical function relation between the main characteristic load and the characteristic variable through data statistics; (7) And comparing with the design evaluation load and process to evaluate the actual service life accumulation condition of the part.
Potential problems with this technical idea include: (1) The sensors arranged at the dangerous positions of the components are easy to generate data deviation or system damage in the long-term operation process and are not easy to find in time, so that the accuracy of data acquisition is influenced, and the error of a service life evaluation result is increased; (2) Because the sensors in different types and different positions are subjected to various and different degrees of working environment interference factors, the generated measurement deviation can also present respective characteristics, and a fixed rule is difficult to summarize, so that a method and a basis for analyzing and rejecting suspicious data by a rule system are difficult to find; (3) If the sensor is reset or replaced, load calibration test work must be carried out again due to the requirement on the compensation setting of the system, so that the normal operation of the unit is influenced.
Especially for wind turbine generator system rotor blade, the problem that faces further includes: (1) In order to not influence the aerodynamic shape of the blade, the sensor can only be installed or embedded inside the cavity of the blade, and once a problem occurs in a place where a worker cannot reach, the monitoring information is lost and the function of monitoring the safety of the local area of the blade is lost if the worker gives up maintenance or replacement, otherwise, windowing operation is required and the structure of the blade is damaged; (2) The disassembly and maintenance work can influence the connection of a series of systems related to the blades, the cost is high, the period is long, certain safety risks exist in related units in the maintenance process, workers are required to carry out high-altitude operation even long-time extravehicular operation by using a hanging basket or a suspension platform without disassembling the maintenance work, the safety risks are greatly increased, and meanwhile, the arrangement precision of the sensor can be greatly influenced by the psychological quality of the workers and the difficulty degree of the working environment; (3) The sensor maintenance work is not easy and can cause problems in the running process, the maintenance or the replacement is needed to be carried out repeatedly at irregular intervals, and the running and maintenance cost is too high; (4) With the capacity maximization of a wind generating set, the length and flexibility of the blades of the wind wheel are increased, the phenomena of flapping vibration bending coupling and bending torsion coupling are more obvious, and a certain main characteristic load is difficult to decouple, so that the precision of the process of deducing the main characteristic load by a characteristic variable is influenced, and the service life evaluation precision is further influenced.
The reasons are relevant, and monitoring and service life evaluation work of wind turbine blades of the wind turbine generator set is very rough for manufacturers of large complete machines, manufacturers of parts, owners and technical service providers at the present stage; correspondingly, after the blade is large in size, the bending moment load in the shimmy direction is obviously increased, and the safety problems of blade trailing edge cracking, connection bolt fracture near the blade root of 0 DEG/180 DEG and the like are more easily caused.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a method and a system for evaluating the service life of a blade of a horizontal-axis wind generating set, which do not need to install a sensor on the blade and can obtain the service life of the blade.
In order to achieve the purpose, the invention adopts the following technical scheme to realize the purpose:
a blade life evaluation method for a horizontal axis wind generating set comprises the following steps:
obtaining a thermal insulation protective coating damage expansion rule curved surface, namely a function relation of the volume reduction rate of the thermal insulation protective coating of the test piece along with the experiment time, through an environmental damage resistance experiment of the thermal insulation protective coating; by the blade leading edge thermal insulation protection coating life calibration method, a blade leading edge thermal insulation protection coating life calibration curved surface is obtained by utilizing a thermal insulation protection coating damage expansion rule curved surface, namely the volume reduction rate of the blade leading edge thermal insulation protection coating is in a function relation with the blade spanwise position and the service time;
according to the volume reduction rate of the thermal insulation protection coating at the front edge of the blade, calibrating a curved surface by using the service life of the thermal insulation protection coating at the front edge of the blade, and calculating to obtain the service time of the thermal insulation protection coating coated at the current time on the front edge of the section where the corresponding spanwise position of the blade is located;
the service time of a certain spanwise position of the blade is the sum of the service time of each previous coating of the heat-insulating protective coating on the front edge of the section where the blade is located, and the service life of the blade is the average value of the service time of each spanwise position of the blade.
Preferably, the calculation formula of the volume reduction rate of the thermal insulation protective coating is as follows:
wherein V% represents a septumThermal protection coating volume reduction rate; v 0 Representing the initial thermal barrier protective coating volume value; v test Indicating the volume value of the thermal insulation protective coating.
Further, the calculation process of the volume value of the heat insulation protective coating is as follows: performing infrared scanning on the measurement region to obtain surface temperature distribution of the measurement region, drawing a surface temperature isotherm distribution diagram of the measurement region, and calculating region area values corresponding to different temperature steps; the surface temperatures of an uncoated heat-insulation protective coating area and a standard coated heat-insulation protective coating area in the surface temperature isotherm distribution diagram are used as references, and heat-insulation protective coating thickness values corresponding to different temperature steps are obtained by inward interpolation; and multiplying and accumulating the area values of the regions corresponding to the different temperature steps and the thickness value of the heat-insulating protective coating to obtain the volume value of the heat-insulating protective coating of the measurement region.
Preferably, the obtaining process of the service life calibration curved surface of the thermal insulation protection coating at the front edge of the blade is as follows: traversing the spanwise positions of the blades to obtain the upper and lower limits of the volume reduction rate of the thermal insulation protective coating at each spanwise position in a certain month after the construction of the thermal insulation protective coating at the front edge of the blades is finished, thereby obtaining the corresponding volume reduction rate of the thermal insulation protective coating; and repeating the process to obtain a function relation curved surface of the volume reduction rate of the thermal insulation protection coating of the front edge of the blade along with the spanwise position and the service time of the blade, namely a service life calibration curved surface of the thermal insulation protection coating of the front edge of the blade.
Further, the process of obtaining the upper and lower limits of the volume reduction rate of the thermal insulation protective coating at a certain spanwise position in a certain month after the construction of the thermal insulation protective coating on the front edge of the blade is as follows: traversing the working condition of the damage coupling of the front edge of the blade at the extended position in a month after the construction of the thermal insulation protective coating on the front edge of the blade is finished, arranging the numerical values of the rate of change of the volume reduction of the thermal insulation protective coating along with the time under the condition that the extended position of the blade is constant in a descending order, and sequentially accumulating the damage conditions generated by hours which occur in the corresponding month in a representative year according to the sequence of the working condition of the damage coupling of the front edge of the blade to obtain the lower limit of the volume reduction rate of the thermal insulation protective coating at the extended position in a month after the construction of the thermal insulation protective coating on the front edge of the blade is finished;
traversing the damage coupling working condition of the front edge of the blade at the spreading position in a month after the construction of the heat insulation protection coating of the front edge of the blade is finished, arranging the numerical values of the change rate of the volume reduction rate of the heat insulation protection coating along with the time under the condition that the spreading position of the blade is constant in an ascending order, arranging the sequence of the damage coupling working conditions of the front edge of the blade, and sequentially accumulating the damage conditions representing the number of hours generated in the month corresponding to the year to obtain the lower limit of the volume reduction rate of the heat insulation protection coating at the spreading position in the month after the construction of the heat insulation protection coating of the front edge of the blade is finished;
and (4) taking the average value of the upper and lower limits of the volume reduction rate of the thermal insulation protective coating at the spanwise position in a month after the construction of the thermal insulation protective coating on the front edge of the blade is finished, and obtaining the corresponding volume reduction rate of the thermal insulation protective coating.
And further, the blade leading edge damage coupling working condition is obtained by coupling an environmental damage factor working condition and a unit operation factor working condition, and the characterization data comprises an environmental damage factor characteristic parameter value, an equivalent inflow angle and an equivalent inflow velocity value corresponding to the section where each spanwise position of the blade is located, and a number of hours representing each month of the year.
Further, the process of acquiring the environmental damage factor working condition and the unit operation factor working condition is as follows: acquiring environmental damage factor composition and characteristic parameters of the front edge of a blade of a wind generating set of a wind power plant, classifying the weather conditions of the wind power plant according to historical data records of a meteorological station, and acquiring the number of hours of each month of a corresponding representative year to form an environmental damage factor working condition;
classifying the unit operation conditions according to the operation data records of the wind generating set, and obtaining the number of hours of each month of the corresponding representative year to form the unit operation factor working conditions;
and calculating to obtain an equivalent inflow angle and an equivalent inflow speed corresponding to the section of each spanwise position of the blade under the operating condition of each unit according to the historical data record of the wind measuring tower of the wind power plant, the operating data record of the wind generating set, the torsional angle distribution of the blade, the initial pitch angle of the blade and the unit control strategy.
A blade life assessment system of a horizontal axis wind generating set comprises:
the function relation obtaining module is used for obtaining a thermal insulation protective coating damage expansion rule curved surface through a thermal insulation protective coating environmental damage resistance experiment, namely the function relation of the volume reduction rate of the thermal insulation protective coating of the test piece along with the experiment time; by the blade leading edge thermal insulation protection coating life calibration method, a blade leading edge thermal insulation protection coating life calibration curved surface is obtained by utilizing a thermal insulation protection coating damage expansion rule curved surface, namely the volume reduction rate of the blade leading edge thermal insulation protection coating is in a function relation with the blade spanwise position and the service time;
the secondary coating service time calculation module is used for calculating the service time of the thermal insulation protection coating of the current coating on the front edge of the section where the corresponding extending position of the blade is located by utilizing the service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade according to the volume reduction rate of the thermal insulation protection coating on the front edge of the blade;
the service life obtaining module of the blade is used for enabling the service time of a certain spanwise position of the blade to be the sum of the service times of the previous thermal insulation protective coatings coated on the front edge of the section where the blade is located, and the service life of the blade is the average value of the service times of the spanwise positions of the blade.
A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor when executing the computer program implementing the steps of the method for blade life assessment of a horizontal axis wind park as claimed in any one of the preceding claims.
A computer-readable storage medium storing a computer program which, when executed by a processor, carries out the steps of the method for blade life assessment of a horizontal-axis wind turbine generator set according to any one of the preceding claims.
Compared with the prior art, the invention has the following beneficial effects:
the service life condition of the blade of the wind generating set is calculated by actually measuring the damage condition of the thermal insulation protection coating of the leading edge region of the blade and combining the service life calibration curved surface of the thermal insulation protection coating of the leading edge of the blade, the damage condition of the blade at different spanwise positions can be measured and evaluated, and according to the discrete condition of the damage result at the different spanwise positions, on one hand, the error of the evaluation life of the blade can be corrected and reduced, on the other hand, the deviation of the actual operation and theoretical design of the blade can be known, and the potential dangerous section and dangerous position in the actual operation can be found; the device and parts can be shared with a blade detection system and a blade ice preventing and removing system, so that the operation and maintenance investment cost of a unit is reduced; the system can share data results with a blade detection system and a monitoring system, so that the working efficiency is improved, and the data availability is improved; the data acquisition does not need a sensor, the technical requirements of arranging and maintaining the sensor on the blade of the wind generating set can be effectively avoided, and the workload, difficulty and risk of operators are reduced.
Furthermore, the volume reduction rate of the thermal insulation protective coating is defined by utilizing the physical property of heat conduction and infrared scanning, the three-dimensional representation of the damage degree of the front edge area of the blade is realized, more abundant feedback information can be obtained compared with the traditional visual observation, and the expansion process of discovering and early warning the damage phenomenon is facilitated; meanwhile, the method can also realize the test and the monitoring of the crack initiation and expansion processes of the structural adhesive at the front edge of the blade, and ensure the structural safety of the blade.
Further, a service life calibration curved surface of the thermal insulation protective coating at the front edge of the blade is constructed, the direct relation between the damage condition of the front edge area of the blade and the service life evaluation of the blade is effectively established, and the effective method for evaluating the service life of the protective coating in the field of wind power is realized.
Further, the environmental damage factors of the front edge of the blade are combined with the running condition of the wind generating set, the equivalent inflow angle and the equivalent inflow speed of the section of the blade under different coupling working conditions are considered, and the damage generation and development process of the front edge area of the blade is simulated in a more detailed and real mode.
Detailed Description
The first part is a service life evaluation system of a blade of a horizontal-axis wind generating set, and the service life evaluation system consists of a heat insulation protective coating, an infrared scanner, a heat collector, an air blower and a flow guide passage. The thermal insulation protective coating covers the front edge area of the blade, the front edge area of the blade is protected from being influenced by environmental damage factors (including erosion, corrosion, aging and the like) in the running process of the wind generating set, and the self damage degree of the thermal insulation protective coating can effectively represent the service life condition of the blade; the infrared scanner is held by a worker or carried by an unmanned aerial vehicle, scans the front edge area of the heating blade, and judges the damage degree of the heat insulation protective coating by sensing the temperature distribution of the surface of the front edge area of the blade; the heat collector is arranged at the root of the blade or in the inner cavity of the hub, and plays a role in heating air in the cavity of the blade; the blower is arranged at the root of the blade or the inner cavity of the hub and plays a role of blowing hot air to the tip of the cavity of the front edge of the blade; the flow guide passage is arranged in the cavity of the front edge of the blade, plays a role in guiding hot air to flow and preserving heat, and realizes the energy consumption saving of the heat collector.
And in the second part, according to the requirements of the blade coating process, a heat-insulating protective coating prepared from the heat-insulating protective coating needs to meet the following requirements: (1) The quality requirements of the protective coating on the front edge of the blade are met in the aspects of thickness deviation, surface roughness deviation and the like; (2) The performance requirements of the protective coating on the front edge of the blade are met in the aspects of hardness, peeling strength, corrosion resistance, aging resistance and the like; (3) In the aspect of heat conduction performance, the heat conduction performance of the composite material is lower than that of the blade composite material.
And in the third part, a thermal insulation protective coating damage degree evaluation method comprises the following steps:
the method comprises the steps of scanning a measurement area by using an infrared scanner, meeting the requirement of obtaining the surface temperature distribution of the measurement area by adjusting the scanning power and resolution of the instrument, drawing a surface temperature isotherm distribution diagram of the measurement area, obtaining area values corresponding to different temperature steps by using an image processing technology, and ensuring that the measurement result meets the requirement of measurement precision.
In the measurement area coated with the heat-insulating protective coating on the surface, under the condition that the back surface is uniformly heated, the difference of the surface temperature is a direct result of different heat-conducting properties shown by different thicknesses of the heat-insulating protective coating. And (3) taking the surface temperatures of the areas which are not coated with the heat-insulating protective coatings and the areas coated with the heat-insulating protective coatings in the standard mode in the surface temperature isotherm distribution diagram as reference, interpolating inwards to obtain the thicknesses of the heat-insulating protective coatings corresponding to different temperature steps, multiplying the thicknesses by the corresponding area values of the areas and accumulating to obtain the volume value of the heat-insulating protective coatings.
Defining the volume reduction rate of the heat insulation protective coating as an evaluation index of the damage degree of the heat insulation protective coating, namely
Wherein V% represents the volume reduction rate of the heat insulation protective coating; v 0 Representing the initial thermal barrier protective coating volume value; v test Indicating the volume value of the thermal insulation protective coating.
And fourthly, performing an environment damage resistance experiment on the heat-insulating protective coating:
(1) The environmental damage factor composition and characteristic parameters of the front edge of the blade of the wind generating set of the wind power plant are investigated, wherein the environmental damage factor composition and the characteristic parameters comprise erosion (sand-blown particle diameter and sand transportation amount), rain erosion (rainfall intensity and rainwater pH value), corrosion (air temperature and humidity and salt mist concentration), aging (ultraviolet intensity) and the like. According to the historical data record of the nearby meteorological station, the weather conditions of the wind power plant are classified, the number of hours of each month of the corresponding representative year is obtained, and the environmental damage factor working condition is formed: for example, erosion conditions, rain erosion conditions, corrosion conditions, aging conditions, combined erosion-aging conditions, combined rain erosion-corrosion conditions, combined corrosion-aging conditions, and the like.
(2) According to the running data record of the wind generating set, the running conditions of the set are classified and the corresponding number of hours of each month in the representative year is obtained, so that the working conditions of the running factors of the set are formed: for example, a shutdown (stationary) condition, a shutdown (idle) condition, a start-stop (brake) condition, a normal operating condition, and the like.
According to the historical data record of the wind measuring tower of the wind power plant, the operation data record of the wind generating set, the torsional angle distribution of the blades, the initial pitch angle of the blades and the unit control strategy (the rotating speed of the wind wheel, the pitch variation logic, the yaw logic and the like), the equivalent inflow angle theta corresponding to the section where the extended position r of the blades under the working condition of each unit operation factor is located is converted in (r) and the equivalent incoming flow velocity v in (r)。
(3) Coupling the environmental damage factor working condition with the unit operation factor working condition, and defining the blade leading edge damage coupling working condition Load i The characterization data comprises the characteristic parameter values of the environmental damage factors and equivalent inflow corresponding to the section where the blade extension position r is locatedAngle theta in (r)| i And equivalent incoming flow velocity v in (r)| i Value representing the number of hours t (m) of occurrence of each month of the year i (where i =1,2,3, \8230;. M =1,2, \8230;. 8230;, 12).
(4) The method comprises the steps of preparing a composite material Test (r) by using blade raw materials according to the geometrical shapes of front edges of sections where different blade spanwise positions r are located, and coating a heat insulation protective coating on the surface of the Test according to the requirements of a blade coating process.
(5) And measuring the volume reduction rate of the heat-insulating protective coating of the test piece by using the method provided in the third part, and recording the volume reduction rate as Vp.
(6) Simulating the damage coupling working condition Load of the blade leading edge defined in the fourth part (3) in the laboratory i And (3) carrying out an environment damage resistance experiment on the heat insulation protective coating at the front edge of the section where the blade spanwise position r is located under the condition to obtain a corresponding heat insulation protective coating damage expansion rule curve, namely a function relation of the volume reduction rate Vp of the heat insulation protective coating of the Test piece Test (r) along with the experiment time T (wherein i =1,2,3, \ 8230; \8230;).
Traversing the unfolding position r of the blade to obtain a corresponding thermal insulation protective coating damage expansion rule curved surface, namely
Vp i ＝f i (r，T) [2]
Wherein, the lower corner mark i represents the damage coupling working condition Load of the front edge of the blade i A serial number; r represents the blade spanwise position; t represents experiment time; vp represents the volume reduction rate of the thermal insulation protective coating of the test piece.
Its inverse function expression is recorded as
T i ＝f i -1 (r，Vp)[3]
And a fifth part, a method for calculating the service life of the thermal insulation protective coating of the blade leading edge comprises the following steps:
(1) During actual operation, the volumetric reduction ratio Vp of the thermal protection coating at the front edge of the blade can also be regarded as a function of the spanwise position r and the service time t of the blade and is recorded as Vp (r, t).
Considering the dispersion of the effect, the volume reduction ratio Vp of the thermal insulation protective coating at the leading edge of the blade is numerically defined as the average value of the upper limit and the lower limit, namely
Wherein, vp + (r, t) represents the volume reduction rate upper limit value of the thermal insulation protective coating at the front edge of the blade; vp - And (r, t) represents the lower limit value of the volume reduction rate of the thermal insulation protective coating at the front edge of the blade.
(2) The construction of the thermal insulation protective coating on the front edge of the blade is used as the initial time t =0, and the volume reduction rate of the thermal insulation protective coating at each spanwise position r is 0, namely
Vp(r，0)＝Vp + (r，0)＝Vp - (r，0)＝0。
(3) The j-1 month after the construction of the thermal insulation protective coating on the front edge of the blade is finished, and the upper limit and the lower limit of the volume reduction rate of the blade are recorded as Vp + (r，j-1)、Vp - (r, j-1) according to the formula [4]The volume reduction ratio Vp (r, j-1) is obtained by calculation.
(4) And the construction of the thermal insulation protective coating on the front edge of the blade is completed in the jth month, which is exactly m months in the year.
Coupling condition Load for damage of blade leading edge i Solving its inverse function expression [3 ]]Partial derivative at coordinate point (r, vp (r, j-1)))
Namely the change rate of the experiment time along with the volume reduction rate under the condition that the spanwise position of the blade is constant, and the larger the value is, the thermal insulation protective coating is in the damage coupling working condition Load of the leading edge of the blade i The less likely damage will occur.
For a certain spanwise position r of the blade * Traversing blade leading edge damage coupling working condition Load i Partial derivative of
The numerical values are arranged in a descending order, and the damage coupling working condition Load of the front edge of the blade is arranged according to the numerical values i Sequentially arranging the counting number t (m) representing m months of the year i Accumulating the generated damage conditions to obtain the spanwise position r of the j month after the construction of the thermal insulation protective coating on the front edge of the blade is finished * Lower limit Vp of volume reduction rate of thermal insulation protective coating - (r * J); load for traversing damage coupling working condition of leading edge of blade i Partial derivative of
The numerical values are arranged in an ascending order, and the damage coupling working condition Load of the front edge of the blade is arranged according to the numerical values i Sequentially arranging the counting number t (m) representing m months of the year i Accumulating the generated damage conditions to obtain the spanwise position r of the j month after the construction of the thermal insulation protective coating on the front edge of the blade is finished * Upper limit Vp of volume reduction rate of heat insulation protective coating + (r * ，j)。
Since damage accumulation is a gradual acceleration process, changing the accumulation order of damage conditions affects the volumetric reduction rate result of the thermal protective coating, thereby defining the upper and lower limit values thereof.
Traversing the unfolding position r of the blade to obtain the upper and lower limits Vp of the volume reduction rate of the thermal insulation protective coating of the blade leading edge in the jth month after the thermal insulation protective coating construction is completed + (r，j)、Vp - (r，j)。
(5) Circulating the steps (3) and (4) to obtain a function relation curved surface of the volume reduction rate Vp of the thermal insulation protection coating at the front edge of the blade along with the unfolding position r and the service time t of the blade, namely a service life calibration curved surface of the thermal insulation protection coating at the front edge of the blade, which is recorded as
Vp＝g(r，t) [5]
Wherein r represents the blade spanwise position; t represents the service time; vp represents the volumetric reduction rate of the thermal barrier coating at the leading edge of the blade.
Its inverse function expression is noted as
t＝g -1 (r，Vp) [6]
(6) Wherein, the damage coupling working condition Load of the front edge of the accumulated blade i Representing the number of hours of occurrence of m months of the year t (m) i The damage condition is generated by adopting the following steps:
i) In the curved surface corresponding to the damage extension rule of the thermal insulation protective coating, an inverse function expression [3 ] is utilized]Find (r, vp) before ) Corresponding experimental time T before I.e. by
T before ＝f i -1 (r，Vp before ) [7]
Wherein, vp before Representing the volume reduction rate of the thermal insulation protective coating before accumulation; t is before The corresponding experimental time before accumulation is indicated.
II) extending the experimental time to t (m) & ltLight i As the experimental time T after accumulation after I.e. by
T after ＝T before +t(m)| i [8]
Wherein, T after The corresponding experimental time after accumulation is indicated.
III) utilizing a function expression [2 ] in a curved surface corresponding to the damage extension rule of the thermal insulation protective coating]Find (r, T) after ) Corresponding thermal insulation protective coating volume reduction rate Vp after I.e. by
Vp after ＝f i (r，T after )[9]
Wherein, vp after Showing the volumetric reduction rate of the thermal protective coating after accumulation.
The sixth part of the invention relates to a method for evaluating the service life of a blade of a horizontal axis wind generating set, which comprises the following steps:
(1) Selecting a heat insulation protective coating, and preparing a heat insulation protective coating environmental damage resistance experiment test piece according to the blade coating process requirement.
(2) And defining a service life calibration curved surface of the thermal insulation protection coating at the front edge of the blade.
(3) And configuring hardware of a blade life calculating device of the horizontal-axis wind generating set. The method comprises the step of preparing a standard heat-insulation protective coating in the front edge area of the blade according to the coating process requirement of the blade.
(4) The construction of the thermal insulation protective coating on the front edge of the blade is finished as the initial time t =0, and the volume reduction rate Vp (r, 0) =0 of the thermal insulation protective coating corresponding to each stretching position r.
(5) The method comprises the following steps of (1) normally operating the wind generating set, periodically evaluating the damage degree of the thermal insulation protective coating on the front edge of the blade, and adopting the following steps:
i) Starting the heat collector to heat air in the blade cavity; and starting the air blower, and enabling the hot air to be fully distributed in the cavity of the front edge of the blade by utilizing the flow guide passage. The temperature of the hot air at each spanwise position in the cavity of the leading edge of the blade is ensured to be obviously higher than the ambient temperature.
And II) operating the wind generating set to stop feathering, so that the wind wheel is in an idle state.
And III) starting a mechanical brake system of the wind generating set, and locking the wind wheel and the pitch angle limiting mechanism.
IV) the infrared scanner is held by the staff or carried by the unmanned aerial vehicle, scans the blade leading edge area and records the measured data.
V) closing the heat collector and the blower.
VI) releasing the wind wheel and the pitch angle limiting mechanism, and closing the mechanical brake system of the wind generating set.
VII) operating the wind generating set to start normally.
VIII) processing the measurement data obtained in the step (5) IV) according to a thermal insulation protective coating damage degree evaluation method, and calculating the volume reduction rate Vp (r, t) of the thermal insulation protective coating of the front edge of the blade so as to display the damage degree of the thermal insulation protective coating of the front edge of the blade.
(6) Calibrating the curved surface according to the service life of the thermal insulation protection coating at the front edge of the blade determined in the step (2), and utilizing an inverse function expression [6 ] of the curved surface]And (5) calculating the volume reduction rate Vp (r, t) result of the thermal insulation protective coating at the front edge of the blade obtained in the step (5) VIII), and calculating the service time t (r) = g of the thermal insulation protective coating coated at the extended position r of the blade -1 (r，Vp)。
(7) According to the technical requirements of wind power plants, the damage degree of the thermal insulation protective coating on the front edge of the blade exceeds a corresponding threshold value, and in order to avoid influencing the aerodynamic performance of the wind turbine blade of a unit and protecting the structural safety of the wind turbine blade, the whole or partial spread position of the front edge area of the blade needs to be recoated with a standard thermal insulation protective coating. And (5) returning to the step (4) for the extended position needing to be recoated, and starting the service time evaluation work of the new-time coating thermal insulation protective coating.
(8) Service time t at blade deployment position r all (r) is the sum of the service time t (r) of the previous thermal-insulation protective coating coated on the front edge of the section where the thermal-insulation protective coating is positioned, namely
t all (r)＝∑t(r) [10]
Wherein, t all (r) represents the time of service at the blade deployment position r; t (r) represents the service time of the thermal insulation protective coating coated at a certain spreading position r of the blade.
Service life t of blade of wind generating set blade Theoretically corresponding to the service time t at each deployment position r all (r) are equal, but because of the influence of uncertain factors such as calculation model deviation, measurement deviation and the like, the service time t is defined as the service time t at each spanwise position r all Average value of (r), i.e.
Wherein, t blade Representing the service life of the blade; l represents the blade length.
The seventh part of the invention relates to a blade life evaluation system of a horizontal axis wind generating set, which comprises:
the function relation obtaining module is used for obtaining a thermal insulation protective coating damage expansion rule curved surface through a thermal insulation protective coating environmental damage resistance experiment, namely the function relation of the volume reduction rate of the thermal insulation protective coating of the test piece along with the experiment time; by the method for calibrating the service life of the thermal insulation protection coating on the front edge of the blade, the service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade, namely the functional relation of the volume reduction rate of the thermal insulation protection coating on the front edge of the blade along with the spanwise position and service time of the blade, is obtained by utilizing the damage expansion rule curved surface of the thermal insulation protection coating.
And the secondary coating service time calculation module is used for calculating the service time of the thermal insulation protection coating of the current coating on the front edge of the section where the corresponding extended position of the blade is located by utilizing the service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade according to the measured volume reduction rate of the thermal insulation protection coating on the front edge of the blade.
The service life obtaining module of the blade is used for enabling the service time of a certain spanwise position of the blade to be the sum of the service times of the previous thermal insulation protective coatings coated on the front edge of the section where the blade is located, and the service life of the blade is the average value of the service times of the spanwise positions of the blade.
In an eighth part, the computer device of the present invention includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor, when executing the computer program, implements the steps of the method for estimating the blade life of the horizontal axis wind turbine generator system.
Ninth, the invention provides a computer readable storage medium, which stores a computer program, wherein the computer program when executed by a processor implements the steps of the method for estimating the service life of the blade of the horizontal axis wind turbine generator system.
The above-mentioned contents are only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited thereby, and any modification made on the basis of the technical idea of the present invention falls within the protection scope of the claims of the present invention.
Claims (9)
1. A service life evaluation method for a blade of a horizontal axis wind generating set is characterized by comprising the following steps:
obtaining a thermal insulation protective coating damage expansion rule curved surface, namely a function relation of the volume reduction rate of the thermal insulation protective coating of the test piece along with the experiment time, through an environmental damage resistance experiment of the thermal insulation protective coating; by the method for calibrating the service life of the thermal insulation protection coating on the front edge of the blade, a service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade is obtained by utilizing a damage extension rule curved surface of the thermal insulation protection coating, namely a function relation of the volume reduction rate of the thermal insulation protection coating on the front edge of the blade along with the spanwise position and the service time of the blade;
the obtaining process of the service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade comprises the following steps: traversing the span-wise position of the blade to obtain the volume reduction rate of the thermal insulation protective coating at each span-wise position in a certain month after the construction of the thermal insulation protective coating at the front edge of the blade is finished; repeating the process to obtain a function relation curve of the volume reduction rate of the thermal insulation protection coating of the front edge of the blade along with the spanwise position and the service time of the blade by using the duration variable, namely a service life calibration curve of the thermal insulation protection coating of the front edge of the blade;
according to the volume reduction rate of the thermal insulation protection coating at the front edge of the blade, calculating the service time of the thermal insulation protection coating coated on the front edge of the section where the corresponding unfolding position of the blade is located at the current time by using the service calibration curved surface of the thermal insulation protection coating at the front edge of the blade;
the service time of a certain spanwise position of the blade is the sum of the service time of each previous coating of the heat-insulating protective coating on the front edge of the section where the blade is located, and the service life of the blade is the average value of the service time of each spanwise position of the blade.
2. The method for evaluating the life of a blade of a horizontal-axis wind turbine generator system according to claim 1, wherein the volume reduction rate of the thermal insulation protective coating is calculated by the formula:
wherein V% represents the volume reduction rate of the heat insulation protective coating; v 0 Representing the initial thermal barrier protective coating volume value; v test Indicating the volume value of the thermal insulation protective coating.
3. The method for evaluating the service life of the blade of the horizontal-axis wind turbine generator set according to claim 2, wherein the volume value of the thermal insulation protective coating is calculated by the following steps: performing infrared scanning on the measurement region to obtain surface temperature distribution of the measurement region, drawing a surface temperature isotherm distribution diagram of the measurement region, and calculating region area values corresponding to different temperature steps; the surface temperatures of an uncoated heat-insulation protective coating region and a standard coated heat-insulation protective coating region in the surface temperature isotherm distribution diagram are used as references, and the thickness values of the heat-insulation protective coating corresponding to different temperature steps are obtained by inward interpolation; and multiplying and accumulating the area values of the regions corresponding to the different temperature steps and the thickness value of the heat-insulation protective coating to obtain the volume value of the heat-insulation protective coating of the measurement region.
4. The method for evaluating the service life of the blade of the horizontal axis wind turbine generator set according to claim 1, wherein the step of obtaining the volume reduction rate of the thermal insulation protective coating at a certain spanwise position in a month after the construction of the thermal insulation protective coating at the front edge of the blade is as follows: traversing the damage coupling working condition of the front edge of the blade at the spreading position in a month after the construction of the heat insulation protection coating of the front edge of the blade is finished, arranging the numerical values of the volume reduction rate of the heat insulation protection coating along with the time under the condition that the spreading position of the blade is constant in a descending order, and sequentially accumulating the damage conditions of the front edge of the blade generated by hours corresponding to the month in the representative year according to the sequence of the damage coupling working condition of the front edge of the blade to obtain the lower limit of the volume reduction rate of the heat insulation protection coating at the spreading position in the month after the construction of the heat insulation protection coating of the front edge of the blade is finished;
traversing the damage coupling working condition of the front edge of the blade at the spreading position in a month after the construction of the heat insulation protection coating of the front edge of the blade is finished, arranging the numerical values of the change rate of the volume reduction rate of the heat insulation protection coating along with the time under the condition that the spreading position of the blade is constant in an ascending order, and sequentially accumulating the damage conditions of the front edge of the blade generated by hours corresponding to the month in the representative order according to the sequence of the damage coupling working condition of the front edge of the blade to obtain the upper limit of the volume reduction rate of the heat insulation protection coating at the spreading position in the month after the construction of the heat insulation protection coating of the front edge of the blade is finished;
and (3) taking the average value of the upper and lower limits of the volume reduction rate of the thermal insulation protective coating at the spreading position in a certain month after the construction of the thermal insulation protective coating on the front edge of the blade is finished, and obtaining the corresponding volume reduction rate of the thermal insulation protective coating.
5. The method for evaluating the service life of the blade of the horizontal-axis wind turbine generator system according to claim 4, wherein the damage coupling working condition of the leading edge of the blade is obtained by coupling an environmental damage factor working condition and a unit operation factor working condition, and the characterization data of the damage coupling working condition comprises an environmental damage factor characteristic parameter value, an equivalent inflow angle and an equivalent inflow velocity value corresponding to a section where each spanwise position of the blade is located, and a number of hours of occurrence in each month of the year.
6. The method for evaluating the service life of the blade of the horizontal-axis wind turbine generator system according to claim 5, wherein the obtaining process of the environmental damage factor working condition and the unit operation factor working condition comprises the following steps: acquiring environmental damage factor composition and characteristic parameters of the front edge of a blade of a wind generating set of a wind power plant, classifying the weather conditions of the wind power plant according to historical data records of a meteorological station, and acquiring the number of hours of each month of a corresponding representative year to form an environmental damage factor working condition;
classifying the unit operation conditions according to the operation data records of the wind generating set, and obtaining the number of hours of each month of the corresponding representative year to form the unit operation factor working condition;
and calculating to obtain an equivalent inflow angle and an equivalent inflow speed corresponding to the section of each spanwise position of the blade under the operating condition of each unit according to the historical data record of the wind measuring tower of the wind power plant, the operating data record of the wind generating set, the torsional angle distribution of the blade, the initial pitch angle of the blade and the unit control strategy.
7. A system for evaluating service life of a blade of a horizontal-axis wind turbine generator system is characterized by comprising:
the function relation acquisition module is used for obtaining a heat insulation protection coating damage expansion rule curved surface through a heat insulation protection coating environmental damage resistance experiment, namely, a function relation of the volume reduction rate of the heat insulation protection coating of the test piece along with the experiment time; by the method for calibrating the service life of the thermal insulation protection coating on the front edge of the blade, a service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade is obtained by utilizing a damage extension rule curved surface of the thermal insulation protection coating, namely a function relation of the volume reduction rate of the thermal insulation protection coating on the front edge of the blade along with the spanwise position and the service time of the blade;
the obtaining process of the service life calibration curved surface of the thermal insulation protection coating on the front edge of the blade comprises the following steps: traversing the span-wise position of the blade to obtain the volume reduction rate of the thermal insulation protective coating at each span-wise position in a certain month after the construction of the thermal insulation protective coating at the front edge of the blade is finished; repeating the process to obtain a function relation curved surface of the volume reduction rate of the thermal insulation protection coating at the front edge of the blade along with the spanwise position and the service time of the blade by using a duration variable, namely a service life calibration curved surface of the thermal insulation protection coating at the front edge of the blade;
the secondary coating service time calculation module is used for calculating the service time of the secondary coating of the front edge of the section where the corresponding unfolding position of the blade is located by using the service life calibration curved surface of the thermal insulation protection coating at the front edge of the blade according to the volume reduction rate of the thermal insulation protection coating at the front edge of the blade;
the service life obtaining module of the blade is used for enabling the service time of a certain spanwise position of the blade to be the sum of the service times of all previous thermal insulation protective coatings coated on the front edge of the section where the blade is located, and the service life of the blade is the average value of the service times of all spanwise positions of the blade.
8. Computer arrangement comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor when executing the computer program performs the steps of the method for blade life assessment of a horizontal axis wind park according to any of claims 1 to 6.
9. A computer-readable storage medium, in which a computer program is stored, which computer program, when being executed by a processor, carries out the steps of a method for assessing the life of a blade of a horizontal-axis wind park according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111007759.9A CN113719428B (en) | 2021-08-30 | 2021-08-30 | Method and system for evaluating service life of blade of horizontal axis wind generating set |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111007759.9A CN113719428B (en) | 2021-08-30 | 2021-08-30 | Method and system for evaluating service life of blade of horizontal axis wind generating set |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113719428A CN113719428A (en) | 2021-11-30 |
| CN113719428B true CN113719428B (en) | 2022-11-04 |
Family
ID=78679345
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111007759.9A Active CN113719428B (en) | 2021-08-30 | 2021-08-30 | Method and system for evaluating service life of blade of horizontal axis wind generating set |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113719428B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001032724A (en) * | 1999-07-19 | 2001-02-06 | Toshiba Corp | On-line life diagnostic system |
| CN101526424A (en) * | 2008-03-06 | 2009-09-09 | 通用电气公司 | A system and a method for detecting gas turbine blade or aircraft engine blade problems in real time |
| CN101776645A (en) * | 2010-01-06 | 2010-07-14 | 湘潭大学 | Simulation test method of thermal fatigue failure of blade with thermal barrier coating |
| KR102099614B1 (en) * | 2018-11-08 | 2020-04-10 | 군산대학교산학협력단 | Method and apparatus for fatigue life evaluation of wind turbine composite blade |
| CN111779638A (en) * | 2020-07-10 | 2020-10-16 | 天津智惠未来科技有限责任公司 | Wind power blade defect danger level assessment method based on fatigue stress strain |
| CN112373072A (en) * | 2020-10-23 | 2021-02-19 | 沈阳航空航天大学 | Repair of composite material wind power blade damage and health monitoring method of repaired structure |
| CN113029405A (en) * | 2021-02-04 | 2021-06-25 | 南京航空航天大学 | Blade bending moment decoupling and calibrating method based on optical fiber strain combined bridge circuit |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6730918B2 (en) * | 2001-12-20 | 2004-05-04 | General Electric Company | Apparatus for determining past-service conditions and remaining life of thermal barrier coatings and components having such coatings |
| CN113011109B (en) * | 2021-01-15 | 2022-05-17 | 浙江大学 | Fatigue analysis method for wind driven generator blade coating considering raindrop erosion |
-
2021
- 2021-08-30 CN CN202111007759.9A patent/CN113719428B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001032724A (en) * | 1999-07-19 | 2001-02-06 | Toshiba Corp | On-line life diagnostic system |
| CN101526424A (en) * | 2008-03-06 | 2009-09-09 | 通用电气公司 | A system and a method for detecting gas turbine blade or aircraft engine blade problems in real time |
| CN101776645A (en) * | 2010-01-06 | 2010-07-14 | 湘潭大学 | Simulation test method of thermal fatigue failure of blade with thermal barrier coating |
| KR102099614B1 (en) * | 2018-11-08 | 2020-04-10 | 군산대학교산학협력단 | Method and apparatus for fatigue life evaluation of wind turbine composite blade |
| CN111779638A (en) * | 2020-07-10 | 2020-10-16 | 天津智惠未来科技有限责任公司 | Wind power blade defect danger level assessment method based on fatigue stress strain |
| CN112373072A (en) * | 2020-10-23 | 2021-02-19 | 沈阳航空航天大学 | Repair of composite material wind power blade damage and health monitoring method of repaired structure |
| CN113029405A (en) * | 2021-02-04 | 2021-06-25 | 南京航空航天大学 | Blade bending moment decoupling and calibrating method based on optical fiber strain combined bridge circuit |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113719428A (en) | 2021-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3139038B1 (en) | A method for estimating the surface condition of a rotating blade | |
| CN101718255B (en) | Method for controlling a wind turbine using a wind flow model | |
| Shu et al. | Study of ice accretion feature and power characteristics of wind turbines at natural icing environment | |
| JP5984791B2 (en) | Wind power generator monitoring system and monitoring method | |
| EP2494192B1 (en) | Wind sensor system using blade signals | |
| JP5033033B2 (en) | Method for measuring turbulence intensity of horizontal axis wind turbine | |
| KR101706508B1 (en) | Data Analysis System of Mechanical Load Measurement Data for Wind Turbine | |
| Ehrmann et al. | Effect of surface roughness on wind turbine performance | |
| CN103745024B (en) | Wind turbines tail portion wind speed power characteristic Evaluation Method is corrected based on three-dimensional wake flow model | |
| EP3181897B1 (en) | Operating a wind turbine | |
| CN108431404B (en) | Method and system for controlling a plurality of wind turbines | |
| JP2020502424A (en) | Wind turbine power station level load management control strategy | |
| CN105403588A (en) | Test system for assessing freeze resistance of wind turbine blade coating, and method thereof | |
| Bai et al. | Aerodynamic design and analysis of a 10 kW horizontal-axis wind turbine for Tainan, Taiwan | |
| Hulsman et al. | Turbine power loss during yaw-misaligned free field tests at different atmospheric conditions | |
| CN113719428B (en) | Method and system for evaluating service life of blade of horizontal axis wind generating set | |
| Machielse et al. | ECN test farm measurements for validation of wake models | |
| Ehrmann | Effect of surface roughness on wind turbine performance | |
| CN116484652B (en) | Wake flow interference detection method in wind power plant based on blade root load | |
| CN205317705U (en) | Anti test system that freezes performance of aassessment wind -powered electricity generation blade coating | |
| Ehrmann et al. | Effect of Blade Roughness on Transition and Wind Turbine Performance. | |
| Li et al. | Research on the Change of Airfoil Geometric Parameters of Horizontal Axis Wind Turbine Blades Caused by Atmospheric Icing | |
| Visbech et al. | Aerodynamic effects of leading edge erosion in wind farm flow modeling | |
| Noppe et al. | High frequent scada-based thrust load modeling of wind turbines | |
| CN114399060B (en) | Wind driven generator blade maintenance method based on wind and rain co-invasion index |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |