CN113719428A - 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
- CN113719428A CN113719428A CN202111007759.9A CN202111007759A CN113719428A CN 113719428 A CN113719428 A CN 113719428A CN 202111007759 A CN202111007759 A CN 202111007759A CN 113719428 A CN113719428 A CN 113719428A
- 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.)
- Granted
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 implementing data acquisition on the characteristic variable; (4) establishing a numerical function relationship 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 the actual service life of the part 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 always easy, problems can continuously occur in the operation process, the maintenance or the replacement is needed to be carried out repeatedly at irregular intervals, and the operation and maintenance cost is too high; (4) with the large capacity of the wind generating set, the length and flexibility of the blades of the wind wheel are increased, the phenomena of flap-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 from the characteristic variable is influenced, and the service life evaluation precision is further influenced.
The reasons are related, and monitoring and service life evaluation work of wind turbine blades of the wind turbine generator set is very rough for various manufacturers of complete machines, manufacturers of components, 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 service life evaluation method for a blade of 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 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;
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 the volume reduction rate of the heat insulation protective coating; v0Representing the initial thermal barrier protective coating volume value; vtestIndicating 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 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-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 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.
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 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, 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 (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.
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 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 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 process 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 protection coating at the front edge of the blade is constructed, a direct relation between the damage condition of the front edge area of the blade and service life evaluation of the blade is effectively established, and an 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 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; v0Representing the initial thermal barrier protective coating volume value; vtestIndicating 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: such as a shutdown (rest) condition, a shutdown (idle) condition, a start-stop (brake) condition, a normal operating condition, and so forth.
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 convertedin(r) and the equivalent incoming flow velocity vin(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 LoadiThe characterization data comprises characteristic parameter values of the environmental damage factors and equivalent inflow angles theta corresponding to the cross sections of the spanwise positions r of the bladesin(r)|iAnd equivalent incoming flow velocity vin(r)|iValue, representing the number of hours t (m) of occurrence of each month of the yeari(wherein i is 1, 2, 3, … …; m is 1, 2, … …, 12).
(4) Preparing a composite material test (r) by using blade raw materials according to the geometrical shapes of the front edges of the sections of different blade spanwise positions r, and coating a heat-insulating protective coating on the surface of the test according to the requirements of the 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 laboratoryiAnd (3) carrying out an environment damage resistance experiment on the heat insulation protection coating at the front edge of the section where the blade extending position r is located to obtain a corresponding heat insulation protection coating damage extension rule curve, namely a function relation of the volume reduction rate Vp of the heat insulation protection coating of the test piece test (r) along with the experiment time T (wherein i is 1, 2, 3 and … …).
Traversing the unfolding position r of the blade to obtain a corresponding thermal insulation protective coating damage expansion rule curved surface, namely
Vpi＝fi(r，T) [2]
Wherein, the lower partThe angle mark i represents the damage coupling condition Load of the leading edge of the bladeiA 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
Ti＝fi -1(r，Vp)[3]
And a fifth part, a blade leading edge heat insulation protection coating life calculation method:
(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) Taking the construction completion of the thermal insulation protective coating on the front edge of the blade as the initial time t equal to 0, and setting the volume reduction rate of the thermal insulation protective coating at each spanwise position r to 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 rate 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 edgeiSolving 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 bladeiThe less likely damage will occur.
For a certain spanwise position r of the blade*Traversing blade leading edge damage coupling working condition LoadiPartial 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 valuesiSequentially arranging the air holes representing the number of generating hours of m months per year t (m)iAccumulating the generated damage conditions to obtain the spanwise position r of the j month after the construction of the thermal insulation protective coating of the front edge of the blade is completed*Lower limit Vp of volume reduction rate of heat insulation protective coating-(r*J); load for traversing damage coupling working condition of leading edge of bladeiPartial 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 valuesiSequentially arranging the air holes representing the number of generating hours of m months per year t (m)iAccumulating the generated damage conditions to obtain the spanwise position r of the j month after the construction of the thermal insulation protective coating of the front edge of the blade is completed*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 curve of the volume reduction rate Vp of the thermal insulation protection coating of the front edge of the blade along with the spanwise position r and the service time t of the blade, namely a service life calibration curve of the thermal insulation protection coating of 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 protective coating at the leading edge of the blade.
Its inverse function expression is recorded as
t＝g-1(r，Vp) [6]
(6) Wherein, accumulating the damage coupling working condition Load of the leading edge of the bladeiRepresenting the number of hours t (m) of occurrence of year m monthsiThe 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 experiment time TbeforeI.e. by
Tbefore＝fi -1(r，Vpbefore) [7]
Wherein, VpbeforeRepresenting the volume reduction rate of the thermal insulation protective coating before accumulation; t isbeforeThe corresponding experimental time before accumulation is indicated.
II) prolonging the experimental time by t (m) hours of occurrenceiAs the experimental time T after accumulationafterI.e. by
Tafter＝Tbefore+t(m)|i [8]
Wherein, TafterThe 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 VpafterI.e. by
Vpafter＝fi(r，Tafter)[9]
Wherein, VpafterShowing 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) and selecting a heat-insulating protective coating, and preparing a heat-insulating protective coating environmental damage resistance experiment test piece according to the coating process requirement of the blade.
(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-insulating protective coating on the front edge area of the blade according to the coating process requirement of the blade.
(4) And taking the construction completion of the thermal insulation protective coating at the front edge of the blade as the initial time t as 0, and setting the volume reduction rate Vp (r, 0) of the thermal insulation protective coating at each spanwise position r as 0.
(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) an infrared scanner is held by a worker or carried by an unmanned aerial vehicle, the leading edge area of the blade is scanned and measured, and measurement data is recorded.
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 at the front edge of the blade so as to display the damage degree of the thermal insulation protective coating at the front edge of the blade.
(6) According to the service life calibration curved surface of the thermal insulation protection coating at the front edge of the blade determined in the step (2), utilizing the inverse function expression of the service life calibration curved surface[6]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 (VIII), and accounting the service time t (r) g of the thermal insulation protective coating coated at the blade spreading position r-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 (4) returning to the step (4) for the spanwise position needing to be recoated, and starting the service time evaluation work of the new coating thermal insulation protective coating.
(8) Service time t at blade deployment position rall(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
tall(r)＝∑t(r) [10]
Wherein, tall(r) represents the time in service at the blade deployment position r; t (r) represents the service time of the thermal insulation protective coating coated at the position r of the blade in the unfolding direction.
Service life t of blade of wind generating setbladeTheoretically corresponding to the service time t at each deployment position rall(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 rallAverage value of (r), i.e.
Wherein, tbladeRepresenting 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 is obtained by utilizing the damage extension rule curved surface of the thermal insulation protection coating, namely the volume reduction rate of the thermal insulation protection coating on the front edge of the blade is in a function relation with the spanwise position and the service time of the blade.
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 present invention provides a computer readable storage medium storing a computer program, which when executed by a processor, implements the steps of the method for estimating the blade life of a horizontal axis wind turbine generator system as described above.
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 (10)
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;
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.
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; v0Representing the initial thermal barrier protective coating volume value; vtestIndicating 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 calculation process of the volume value of the thermal 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 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-insulating protective coating to obtain the volume value of the heat-insulating protective coating of the measurement region.
4. The method for evaluating the service life of the blade of the horizontal axis wind generating set according to claim 1, wherein the obtaining process of the service life calibration curved surface of the thermal insulation protection coating at the leading 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 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.
5. The method for evaluating the service life of the blade of the horizontal axis wind turbine generator set according to claim 4, wherein the step 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 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, 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 (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.
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 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.
7. The method for evaluating the service life of the blade of the horizontal axis wind generating set according to claim 6, wherein the obtaining process of the environmental damage factor condition and the unit operation factor 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.
8. A life-span evaluation system of a blade of a horizontal axis wind generating set is characterized by comprising:
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, 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 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.
9. 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 7.
10. 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 the method for blade life assessment of a horizontal-axis wind park according to any one of claims 1 to 7.
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 true CN113719428A (en) | 2021-11-30 |
| CN113719428B 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 (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001032724A (en) * | 1999-07-19 | 2001-02-06 | Toshiba Corp | On-line life diagnostic system |
| US20030115941A1 (en) * | 2001-12-20 | 2003-06-26 | General Electric Crd | Thermal barrier coatings, components, method and apparatus for determining past-service conditions and remaining life thereof |
| 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 |
| CN113011109A (en) * | 2021-01-15 | 2021-06-22 | 浙江大学 | Fatigue analysis method for wind driven generator blade coating considering raindrop erosion |
| CN113029405A (en) * | 2021-02-04 | 2021-06-25 | 南京航空航天大学 | Blade bending moment decoupling and calibrating method based on optical fiber strain combined bridge circuit |
-
2021
- 2021-08-30 CN CN202111007759.9A patent/CN113719428B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001032724A (en) * | 1999-07-19 | 2001-02-06 | Toshiba Corp | On-line life diagnostic system |
| US20030115941A1 (en) * | 2001-12-20 | 2003-06-26 | General Electric Crd | Thermal barrier coatings, components, method and apparatus for determining past-service conditions and remaining life thereof |
| 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 |
| CN113011109A (en) * | 2021-01-15 | 2021-06-22 | 浙江大学 | Fatigue analysis method for wind driven generator blade coating considering raindrop erosion |
| 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 |
|---|---|
| CN113719428B (en) | 2022-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10202965B2 (en) | Method and system for estimating the surface condition of a rotating blade | |
| Shu et al. | Study of ice accretion feature and power characteristics of wind turbines at natural icing environment | |
| CN101718255B (en) | Method for controlling a wind turbine using a wind flow model | |
| JP5984791B2 (en) | Wind power generator monitoring system and monitoring method | |
| CN103745024B (en) | Wind turbines tail portion wind speed power characteristic Evaluation Method is corrected based on three-dimensional wake flow model | |
| JP5033033B2 (en) | Method for measuring turbulence intensity of horizontal axis wind turbine | |
| EP2494192B1 (en) | Wind sensor system using blade signals | |
| CN107194097B (en) | Analysis method based on wind power plant aerodynamic simulation and wind speed and direction data | |
| Ehrmann et al. | Effect of surface roughness on wind turbine performance | |
| CN108431404B (en) | Method and system for controlling a plurality of wind turbines | |
| EP3181897B1 (en) | Operating a wind turbine | |
| CN105403588A (en) | Test system for assessing freeze resistance of wind turbine blade coating, and method thereof | |
| 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 | |
| Ronsten et al. | State-of-the-art of Wind Energy in Cold Climates | |
| CN205317705U (en) | Anti test system that freezes performance of aassessment wind -powered electricity generation blade coating | |
| Sant et al. | Estimating the angle of attack from blade pressure measurements on the National Renewable Energy Laboratory phase VI rotor using a free wake vortex model: yawed conditions | |
| Pedersen et al. | Using wind speed from a blade-mounted flow sensor for power and load assessment on modern wind turbines | |
| Dai et al. | Study on Obtaining Real Power Curve of Wind Turbines Using SCADA Data | |
| CN112036663A (en) | Wind power plant wind energy resource assessment method in alpine mountain area | |
| Li et al. | Research on the Change of Airfoil Geometric Parameters of Horizontal Axis Wind Turbine Blades Caused by Atmospheric Icing | |
| Noppe et al. | High frequent scada-based thrust load modeling of wind turbines |
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 |