CN109973325A - Method and apparatus for identifying abnormal vibration - Google Patents
Method and apparatus for identifying abnormal vibration Download PDFInfo
- Publication number
- CN109973325A CN109973325A CN201711384895.3A CN201711384895A CN109973325A CN 109973325 A CN109973325 A CN 109973325A CN 201711384895 A CN201711384895 A CN 201711384895A CN 109973325 A CN109973325 A CN 109973325A
- Authority
- CN
- China
- Prior art keywords
- frequency
- predetermined component
- rotary speed
- speed data
- vibration
- 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
- 230000002159 abnormal effect Effects 0.000 title claims abstract description 231
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000001133 acceleration Effects 0.000 claims abstract description 191
- 238000001228 spectrum Methods 0.000 claims abstract description 99
- 238000006243 chemical reaction Methods 0.000 claims abstract description 21
- 238000005096 rolling process Methods 0.000 claims description 62
- 230000006870 function Effects 0.000 claims description 52
- 210000003781 tooth socket Anatomy 0.000 claims description 42
- 238000012417 linear regression Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 13
- 238000012216 screening Methods 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims 2
- 238000011156 evaluation Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007873 sieving Methods 0.000 description 1
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
-
- 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
- F03D80/00—Details, components or accessories not provided for in groups F03D1/00 - F03D17/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Life Sciences & Earth Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Wind Motors (AREA)
Abstract
A method and apparatus for identifying abnormal vibration are provided, the method including: acquiring operation data of a preset component of a wind generating set in a plurality of time periods, wherein the operation data comprises vibration acceleration data of the preset component and rotating speed data related to the preset component; respectively carrying out frequency domain conversion on the vibration acceleration data of the predetermined component in the plurality of time periods to obtain a plurality of acceleration frequency spectrums corresponding to each other; determining a frequency value for abnormal vibration analysis in each acceleration frequency spectrum; and determining whether the predetermined component has abnormal vibration or not based on the determined frequency value for the abnormal vibration analysis in each acceleration frequency spectrum and the rotation speed data related to the predetermined component in each time period. By adopting the method and the equipment for identifying the abnormal vibration, the preset component with the abnormal vibration in the wind generating set can be timely and accurately positioned, and powerful support is provided for effectively evaluating the vibration state of the preset component.
Description
Technical field
The present invention relates to technical field of wind power generation, more particularly, be related to it is a kind of identify abnormal vibrations method and set
It is standby.
Background technique
Wind power generating set is generally positioned at remote wind power plant, and the main component in wind power generating set is run reliable
Property, stability are most important.To enable a device to safe and stable, long period, oepration at full load, the fortune of equipment need to be understood in time
Row state, trouble saving, prevent accident, extend the equipment cycle of operation, shorten maintenance time, it is maximum excavate equipment productive potentialities with
And it prevents in time, the operation work of, indefinite period safe and stable to the main component in wind power generating set, continuous live load
Condition timely control.
In longtime running, due to the Special geographical position of wind power plant, running environment, the concealment of problem, live in addition
Staff is difficult to timely to find security risk that the main component in wind power generating set occurs in power generation process (if not
Apparent abnormal vibrations).Wherein, the fundamental frequency of wind-driven generator and its multiple-frequency vibration are common vibration modes, can be to wind
Power generator or even entire wind power generating set bring different degrees of damage, will lead to wind-driven generator failure when serious.
Existing fundamental frequency recognition methods general algorithm is complex, needs largely to calculate and signal is analyzed, be used for
Efficiency is lower when high-volume data identify.
Summary of the invention
Exemplary embodiment of the present invention is designed to provide a kind of method and apparatus for identifying abnormal vibrations, to solve
Can not find in time in the prior art the main component in wind power generating set there are abnormal vibrations, to abnormal vibrations recognition efficiency
Low technical problem.
One side according to an exemplary embodiment of the present invention provides a kind of method for identifying abnormal vibrations, the method packet
It includes: obtaining operation data of the predetermined component of wind power generating set within multiple periods, the operation data includes described pre-
Determine component vibration acceleration data and rotary speed data relevant to the predetermined component;To the predetermined component the multiple
Vibration acceleration data in period carry out frequency domain conversion respectively, to obtain corresponding multiple acceleration frequency spectrums;It determines
Frequency values in each acceleration frequency spectrum for abnormal vibrations analysis;It shakes based in determining each acceleration frequency spectrum for abnormal
Frequency values of dynamic analysis and rotary speed data relevant to the predetermined component determines that the predetermined component is in each period
It is no that there are abnormal vibrations.
Optionally, the type of abnormal vibrations may include that the fundamental vibration of the predetermined component is abnormal and the predetermined component
Multiple-frequency vibration is abnormal.
Optionally, the step of frequency values of abnormal vibrations analysis are used in each acceleration frequency spectrum of the determination can include:
Search the Frequency point that frequency range value in the acceleration frequency spectrum is greater than Frequency and Amplitude threshold value;It will be corresponding with the Frequency point found
Frequency values are as the frequency values for being used for abnormal vibrations analysis in the acceleration frequency spectrum.
Optionally, the step of frequency values of abnormal vibrations analysis are used in each acceleration frequency spectrum of the determination can include:
Determine whether Frequency and Amplitude value corresponding with default concern Frequency point is greater than Frequency and Amplitude threshold value in the acceleration frequency spectrum;If with
The corresponding Frequency and Amplitude value of the default concern Frequency point is greater than Frequency and Amplitude threshold value, then will be with the default concern Frequency point pair
The frequency values answered are as the frequency values for being used for abnormal vibrations analysis in the acceleration frequency spectrum.
Optionally, it is described it is default concern Frequency point can for will include in acceleration frequency spectrum all Frequency points according to frequency width
The size descending of angle value arranges, the Frequency point before predetermined quantity.
Optionally, described based on the frequency values for being used for abnormal vibrations analysis in determining each acceleration frequency spectrum and each
Rotary speed data relevant with the predetermined component determines that the step of predetermined component whether there is abnormal vibrations can in period
Comprise determining that whether meet the predetermined linear regularity of distribution between the frequency values and the rotary speed data;When the frequency values with
When meeting the predetermined linear regularity of distribution between the rotary speed data, determine that there are abnormal vibrations for the predetermined component.
Optionally, whether meet the predetermined linear regularity of distribution between the determination frequency values and the rotary speed data
Step can include: draw revolving speed-frequency scatter plot based on the frequency values and the rotary speed data, wherein the revolving speed-frequency
A scatterplot in scatter plot can correspond in the rotary speed data and acceleration frequency spectrum corresponding with the period of a period
A frequency values for abnormal vibrations analysis;Choose the scatterplot in the preset range in preset frequency linear model;It is based on
The scatterplot of selection determines whether meet the predetermined of the preset frequency linear model between the frequency values and the rotary speed data
Linear distribution rule.
Optionally, the scatterplot based on selection determine whether meet between the frequency values and the rotary speed data it is described
The step of predetermined linear regularity of distribution of preset frequency linear model can include: to the frequency values corresponding with the scatterplot of selection
Linear regression is carried out with the rotary speed data, obtains the model parameter of linear regression;The model parameter is calculated to make a reservation for described
The difference of the specified parameter of component;When the difference is not more than the first setting value, determine corresponding with the scatterplot chosen described
Meet the predetermined linear regularity of distribution of the preset frequency linear model between frequency values and the rotary speed data.
Optionally, the scatterplot based on selection determine whether meet between the frequency values and the rotary speed data it is described
The step of predetermined linear regularity of distribution of preset frequency linear model can include: establish objective function, the objective function instruction
Each scatterplot to the preset frequency linear model distance;By by the frequency values corresponding with the scatterplot of selection with it is described
Rotary speed data substitutes into objective function, obtains the value of the objective function;When the value of the objective function is not more than setting value, really
It is fixed to meet the pre- of the preset frequency linear model between the corresponding frequency values of scatterplot and the rotary speed data chosen
The constant linear regularity of distribution.
Optionally, whether meet the predetermined linear regularity of distribution between the determination frequency values and the rotary speed data
Step can include: the revolving speed statistical value of the reflection data characteristics of each period is calculated based on the rotary speed data;Described in determination
Whether meet the predetermined linear regularity of distribution between frequency values and the revolving speed statistical value.
Optionally, the revolving speed statistical value of the reflection data characteristics of each period may include any in following item
: the average value of rotary speed data relevant with the predetermined component in the period, in the period and the predetermined component phase
The median of the rotary speed data of pass, in the period rotary speed data relevant to the predetermined component virtual value.
Optionally, the vibration acceleration data to the predetermined component within the multiple period carry out frequency domain turn
It changes, the step of to obtain corresponding multiple acceleration frequency spectrums can include: determine related to the predetermined component in each period
Rotary speed data whether be within the scope of setting speed;If rotary speed data relevant to the predetermined component in any time period
Within the scope of setting speed, then the vibration acceleration data to the predetermined component in any time period carry out frequency domain
Conversion, to obtain a corresponding acceleration frequency spectrum.
Optionally, whether rotary speed data relevant with the predetermined component is in setting turn in the determination each period
Step in fast range can include: calculate in each period the standard deviation of rotary speed data relevant with the predetermined component or surely
State error;If the standard deviation or steady-state error that calculate for any time period are within the scope of given threshold, it is determined that described
Rotary speed data relevant to the predetermined component is within the scope of setting speed in any time period;If being directed to any time period
The standard deviation or steady-state error of calculating are not within the scope of given threshold, it is determined that in any time period with the reservations
The relevant rotary speed data of part is not within the scope of setting speed.
Optionally, the predetermined component may include any one of following item: wind-driven generator, generator tooth socket, gear
Case, rolling bearing, wherein when the predetermined component is wind-driven generator, the vibration acceleration data of the predetermined component can
For the vibration acceleration data of engine rooms of wind power generators, rotary speed data relevant to the predetermined component can be wind-driven generator
Rotary speed data, the type of abnormal vibrations may include the fundamental frequency abnormal vibration of wind-driven generator and the multiple-frequency vibration of fundamental frequency
It is abnormal;When the predetermined component is generator tooth socket, the vibration acceleration data of the predetermined component can be wind-driven generator
The vibration acceleration data of cabin, rotary speed data relevant to the predetermined component can be the rotary speed data of wind-driven generator, different
The type often vibrated may include the multiple-frequency vibration exception of generator tooth socket frequency abnormal vibration and generator tooth socket frequency;When described
When predetermined component is gear-box, the vibration acceleration data of the predetermined component can be the power gear or driven tooth in gear-box
The vibration acceleration data of wheel, rotary speed data relevant to the predetermined component can be axis where power gear or driven gear
Rotary speed data, the type of abnormal vibrations may include that the multiple-frequency vibration of power gear meshing frequency and power gear meshing frequency is different
Often or the multiple-frequency vibration of driven gear meshing frequency and driven gear meshing frequency is abnormal;When the predetermined component is to roll
When bearing, the vibration acceleration data of the predetermined component can be the vibration acceleration data of the bearing block of rolling bearing, with institute
The rotary speed data that the relevant rotary speed data of predetermined component can be rolling bearing is stated, the type of abnormal vibrations may include rolling bearing event
The multiple-frequency vibration for hindering characteristic frequency abnormal vibration and Rolling Bearing Fault Character frequency is abnormal.
Another aspect according to an exemplary embodiment of the present invention provides a kind of equipment for identifying abnormal vibrations, the equipment
It include: that operation data obtains module, for obtaining operation data of the predetermined component of wind power generating set within multiple periods,
The operation data include the predetermined component vibration acceleration data and rotary speed data relevant to the predetermined component;When
Frequency conversion module turns for carrying out frequency domain respectively to vibration acceleration data of the predetermined component within the multiple period
It changes, to obtain corresponding multiple acceleration frequency spectrums;Frequency values determining module is used for for determining in each acceleration frequency spectrum
The frequency values of abnormal vibrations analysis;Abnormal vibrations analysis module, for abnormal based on being used in determining each acceleration frequency spectrum
Rotary speed data relevant to the predetermined component determines the predetermined component in the frequency values of vibration analysis and each period
With the presence or absence of abnormal vibrations.
Optionally, the type of abnormal vibrations may include that the fundamental vibration of the predetermined component is abnormal and the predetermined component
Multiple-frequency vibration is abnormal.
Optionally, the frequency values determining module can be used for searching frequency range value in the acceleration frequency spectrum and be greater than frequency width
The Frequency point for spending threshold value divides using frequency values corresponding with the Frequency point found as abnormal vibrations are used in the acceleration frequency spectrum
The frequency values of analysis.
Optionally, the frequency values determining module can be used for determining corresponding with default concern Frequency point in the acceleration frequency spectrum
Frequency and Amplitude value whether be greater than Frequency and Amplitude threshold value, if Frequency and Amplitude value corresponding with the default concern Frequency point is greater than
Frequency and Amplitude threshold value, then frequency values determining module will frequency values corresponding with the default concern Frequency point as acceleration frequency
Frequency values in spectrum for abnormal vibrations analysis.
Optionally, it is described it is default concern Frequency point can for will include in acceleration frequency spectrum all Frequency points according to frequency width
The size descending of angle value arranges, the Frequency point before predetermined quantity.
Optionally, the abnormal vibrations analysis module can be used for determining between the frequency values and the rotary speed data whether
Meet the predetermined linear regularity of distribution, when meeting the predetermined linear regularity of distribution between the frequency values and the rotary speed data
When, determine that there are abnormal vibrations for the predetermined component.
Optionally, the abnormal vibrations analysis module can include: scatter plot drafting module, for based on the frequency values and
The rotary speed data draws revolving speed-frequency scatter plot, wherein a scatterplot in the revolving speed-frequency scatter plot can correspond to one
A frequency in the rotary speed data of a period and acceleration frequency spectrum corresponding with the period for abnormal vibrations analysis
Value;Scatterplot screening module, for choosing the scatterplot being in the preset range of preset frequency linear model;Linear distribution determines mould
Block determines that whether to meet the preset frequency between the frequency values and the rotary speed data linear for the scatterplot based on selection
The predetermined linear regularity of distribution of model.
Optionally, the linear distribution determining module can be used for the frequency values corresponding with the scatterplot of selection and institute
It states rotary speed data and carries out linear regression, obtain the model parameter of linear regression, calculate the model parameter and the predetermined component
The difference of specified parameter determine the corresponding frequency of scatterplot with selection when the difference is not more than the first setting value
Meet the predetermined linear regularity of distribution of the preset frequency linear model between value and the rotary speed data.
Optionally, the linear distribution determining module can include: objective function setting up submodule, for establishing target letter
Number, the objective function indicate each scatterplot to the preset frequency linear model distance root-mean-square value;Target function value
Computational submodule, for by the way that the frequency values corresponding with the scatterplot of selection and the rotary speed data are substituted into objective function,
Obtain the value of the objective function;The regularity of distribution determines submodule, for when the value of the objective function be not more than setting value when,
It determines and meets the preset frequency linear model between the corresponding frequency values of scatterplot and the rotary speed data chosen
The predetermined linear regularity of distribution.
Optionally, the abnormal vibrations analysis module can be used for calculating the anti-of each period based on the rotary speed data
The revolving speed statistical value for reflecting data characteristics determines whether meet predetermined linear distribution between the frequency values and the revolving speed statistical value
Rule.
Optionally, the revolving speed statistical value of the reflection data characteristics of each period may include any one of following item: should
The average value of rotary speed data relevant with the predetermined component in period, the period are interior with relevant turn of the predetermined component
The median of fast data, in the period rotary speed data relevant to the predetermined component virtual value.
Optionally, the time-frequency convert module, can be used for determining in each period to relevant turn of the predetermined component
Whether fast data are within the scope of setting speed, if rotary speed data relevant to the predetermined component is in any time period
Within the scope of setting speed, then the vibration acceleration data to the predetermined component in any time period carry out frequency domain turn
It changes, to obtain a corresponding acceleration frequency spectrum.
Optionally, be used for time-frequency convert module, can be used for calculating in each period to relevant turn of the predetermined component
The standard deviation or steady-state error of fast data, if the standard deviation or steady-state error for any time period calculating are in given threshold
In range, then it can determine that rotary speed data relevant to the predetermined component is in setting speed range in any time period
It is interior, if the standard deviation or steady-state error for any time period calculating are not within the scope of given threshold, can determine described
Rotary speed data relevant to the predetermined component is not within the scope of setting speed in any time period.
Optionally, the predetermined component may include any one of following item: wind-driven generator, generator tooth socket, gear
Case, rolling bearing, wherein when the predetermined component is wind-driven generator, the vibration acceleration data of the predetermined component can
For the vibration acceleration data of engine rooms of wind power generators, rotary speed data relevant to the predetermined component can be wind-driven generator
Rotary speed data, the type of abnormal vibrations may include the fundamental frequency abnormal vibration of wind-driven generator and the multiple-frequency vibration of fundamental frequency
It is abnormal;When the predetermined component is generator tooth socket, the vibration acceleration data of the predetermined component can be wind-driven generator
The vibration acceleration data of cabin, rotary speed data relevant to the predetermined component can be the rotary speed data of wind-driven generator, different
The type often vibrated may include the multiple-frequency vibration exception of generator tooth socket frequency abnormal vibration and generator tooth socket frequency;When described
When predetermined component is gear-box, the vibration acceleration data of the predetermined component can be the power gear or driven tooth in gear-box
The vibration acceleration data of wheel, rotary speed data relevant to the predetermined component can be axis where power gear or driven gear
Rotary speed data, the type of abnormal vibrations may include that the multiple-frequency vibration of power gear meshing frequency and power gear meshing frequency is different
Often or the multiple-frequency vibration of driven gear meshing frequency and driven gear meshing frequency is abnormal;When the predetermined component is to roll
When bearing, the vibration acceleration data of the predetermined component can be the vibration acceleration data of the bearing block of rolling bearing, with institute
The rotary speed data that the relevant rotary speed data of predetermined component can be rolling bearing is stated, the type of abnormal vibrations may include rolling bearing event
The multiple-frequency vibration for hindering characteristic frequency abnormal vibration and Rolling Bearing Fault Character frequency is abnormal.
It is according to an exemplary embodiment of the present invention a kind of to be stored with the computer-readable of computer program in another aspect, providing
Storage medium, when the method that the computer program realizes above-mentioned identification abnormal vibrations when being executed by processor.
It is according to an exemplary embodiment of the present invention in another aspect, provide a kind of computing device, the computing device includes: place
Manage device;Memory is stored with computer program, when the computer program is executed by processor, realizes that above-mentioned identification is different
The method often vibrated.
It, being capable of positioning wind promptly and accurately using the method and apparatus of the identification abnormal vibrations of exemplary embodiment of the present
There are the predetermined components of abnormal vibrations in power generator group, provide effectively effectively to assess the vibrational state of the predetermined component
Support.
Detailed description of the invention
Pass through the detailed description carried out below with reference to the attached drawing for being exemplarily illustrated embodiment, exemplary embodiment of the present
Above and other objects, features and advantages will become apparent.
Fig. 1 shows the flow chart of the method for identification abnormal vibrations according to an exemplary embodiment of the present invention;
The step of whether meeting the predetermined linear regularity of distribution Fig. 2 shows determining scatterplot according to an exemplary embodiment of the present invention
Flow chart;
Revolving speed according to an exemplary embodiment of the present invention-frequency scatter plot exemplary diagram is shown respectively in Fig. 3 A to Fig. 3 D;
Fig. 4, which is shown, according to an exemplary embodiment of the present invention determines whether scatterplot meets based on preset frequency linear model
The flow chart of the step of predetermined linear regularity of distribution;
Fig. 5 shows that according to an exemplary embodiment of the present invention to determine whether scatterplot meets based on objective function predetermined linear
The flow chart of the step of regularity of distribution;
Fig. 6 shows the structure chart of the equipment of identification abnormal vibrations according to an exemplary embodiment of the present invention;
Fig. 7 shows the structure chart of abnormal vibrations analysis module according to an exemplary embodiment of the present invention;
Fig. 8 shows the structure chart of linear distribution determining module according to an exemplary embodiment of the present invention.
Specific embodiment
Now, different example embodiments is more fully described with reference to the accompanying drawings, some exemplary embodiments are in the accompanying drawings
It shows.
Fig. 1 shows the flow chart of the method for identification abnormal vibrations according to an exemplary embodiment of the present invention.
Referring to Fig.1, in step slo, operation number of the predetermined component of wind power generating set within multiple periods is obtained
According to.Here, the operation data may include predetermined component vibration acceleration data and rotary speed data relevant to predetermined component.
Preferably, the vibration acceleration data of predetermined component may include the vibration acceleration data and of the first predetermined direction
The vibration acceleration data of two predetermined directions.As an example, the first predetermined direction can refer to from the head of wind power generating set to tail
The direction in portion, the second predetermined direction can refer to the direction vertical with wind direction (for example, field personnel station is on the leeward, towards machine
Head, the left and right directions of field personnel may be defined as the second predetermined direction).
In exemplary embodiments of the present invention, can be made a reservation for by the vibration acceleration data and second to the first predetermined direction
The vibration acceleration data in direction are respectively processed, and whether there is abnormal vibration to judge predetermined component in the first predetermined direction causes
Dynamic or predetermined component whether there is abnormal vibrations on the second predetermined direction.
Here, the type of abnormal vibrations may include that the fundamental vibration of predetermined component is abnormal different with the multiple-frequency vibration of predetermined component
Often.As an example, the predetermined component may include any one of following item: wind-driven generator, generator tooth socket, gear-box,
Rolling bearing.
In the first embodiment, predetermined component can be the wind-driven generator of wind power generating set.In the case, reservations
The vibration acceleration data of part can be the vibration acceleration data of engine rooms of wind power generators, rotary speed data relevant to predetermined component
It can be the rotary speed data of wind-driven generator.Correspondingly, the type of abnormal vibrations may include the fundamental frequency vibration of wind-driven generator
Abnormal and fundamental frequency multiple-frequency vibration is abnormal.
In a second embodiment, predetermined component can be generator tooth socket (such as generator unit stator tooth socket).In the case, in advance
The vibration acceleration data for determining component can be the vibration acceleration data of engine rooms of wind power generators, revolving speed relevant to predetermined component
Data can be the rotary speed data of wind-driven generator.Correspondingly, the type of abnormal vibrations may include that the vibration of generator tooth socket frequency is different
Often and/or the multiple-frequency vibration of generator tooth socket frequency is abnormal.Here, since generator tooth socket frequency is equal to determining for wind-driven generator
Sub- tooth socket number × wind-driven generator revolving speed/60, and the stator tooth slot number of wind-driven generator is greater than the number of magnetic pole pairs of wind-driven generator,
Therefore, generator tooth socket frequency (or frequency multiplication) is greater than the fundamental frequency (or frequency multiplication) of wind-driven generator.That is, can be by dividing
The data of each acceleration frequency spectrum medium-high frequency part (part corresponding with generator tooth socket frequency and/or frequency multiplication) are analysed to determine wind
Power generator tooth socket is with the presence or absence of abnormal vibrations (with the presence or absence of generator tooth socket frequency abnormal vibration and/or generator tooth socket frequency
The multiple-frequency vibration of rate is abnormal), by analyze low frequency part in each acceleration frequency spectrum (with the fundamental frequency of wind-driven generator and/
Or the corresponding part of frequency multiplication) data come determine wind-driven generator with the presence or absence of abnormal vibrations (with the presence or absence of wind-driven generator
Fundamental frequency abnormal vibration and/or the multiple-frequency vibration of fundamental frequency are abnormal).
In the third embodiment, predetermined component can be gear-box, include power gear and driven gear in gear-box.Herein
In the case of, the vibration acceleration data of predetermined component can be the vibration acceleration data or driven tooth of the power gear in gear-box
The vibration acceleration data of wheel, rotary speed data relevant to predetermined component can be the rotary speed data or driven of axis where power gear
The rotary speed data of axis where gear.Correspondingly, the type of abnormal vibrations may include power gear meshing frequency abnormal vibration, power
The multiple-frequency vibration of gear mesh frequency is abnormal, times of driven gear meshing frequency abnormal vibration and/or driven gear meshing frequency
Frequency vibration is abnormal.Here, the meshing frequency of power gear (or driven gear) is equal to the number of teeth (or the driven gear of power gear
The number of teeth) revolving speed/60 of axis where × power gear (or driven gear).
In the fourth embodiment, predetermined component can be rolling bearing, can refer to multiple rollings in wind power generating set here
Any rolling bearing in bearing.In the case, the vibration acceleration data of predetermined component can be the bearing block of rolling bearing
Vibration acceleration data, rotary speed data relevant to predetermined component can be rolling bearing rotary speed data.
Correspondingly, the type of abnormal vibrations may include Rolling Bearing Fault Character frequency abnormal vibration and/or rolling bearing
The multiple-frequency vibration of fault characteristic frequency is abnormal.Here, it should be appreciated that rolling bearing may include bearing inner race, outer ring, rolling element, guarantor
Frame is held, correspondingly, predetermined component can refer to bearing inner race, outer ring, rolling element or retainer.For example, using predetermined component as the axis of rolling
For the rolling element held, can be to identify rolling element based on the method for the identification abnormal vibrations of exemplary embodiment of the present at this time
It is no that there are the multiple-frequency vibration of the fault characteristic frequency abnormal vibrations of rolling element or fault characteristic frequency exceptions.Here, bearing inner race,
Outer ring, rolling element, retainer respectively correspond respective failure coefficient, can get failure spy with the product for turning frequency by failure coefficient
Levy frequency.
In step S20, frequency domain conversion is carried out respectively to vibration acceleration data of the predetermined component within multiple periods,
To obtain corresponding multiple acceleration frequency spectrums.
For example, frequency domain can be carried out to vibration acceleration data of the predetermined component in any time period in multiple periods
Conversion is to obtain acceleration frequency spectrum corresponding with the vibration acceleration data in any time period, that is, in a period
The corresponding acceleration frequency spectrum of vibration acceleration data.As an example, can be existed by Fast Fourier Transform (FFT) to predetermined component
Vibration acceleration data in multiple periods carry out frequency domain conversion, however the invention is not limited thereto, and other modes also can be used
To carry out frequency domain conversion.
It preferably, can be based on rotary speed data relevant to predetermined component in advance to the predetermined component of acquisition in multiple periods
Interior operation data is screened, by vibration acceleration of the predetermined component within multiple periods in the operation data after screening
Data carry out frequency domain conversion respectively.
For example, it may be determined that whether rotary speed data relevant to predetermined component is in setting speed range in each period
It is interior, if rotary speed data relevant to predetermined component is within the scope of setting speed in any time period, exist to predetermined component
Vibration acceleration data in any time period carry out frequency domain conversion, to obtain a corresponding acceleration frequency spectrum.Here,
Since rotary speed data is variations per hour in wind power generating set, the fundamental frequency or frequency multiplication of predetermined component are related to revolving speed, in revolving speed wave
Move to be completed in biggish situation fundamental frequency or frequency multiplication abnormal vibrations identification there are larger difficulty.Therefore in order to improve to fundamental frequency
Or the abnormal vibrations accuracy of identification of frequency multiplication, the fluctuation range of rotary speed data can be limited, that is, so that for carrying out abnormal vibrations knowledge
The fluctuation of speed of other operation data is smaller.
As an example, determining whether rotary speed data relevant to predetermined component is in setting speed range in each period
Interior step can include: calculate the standard deviation or steady-state error of rotary speed data relevant to predetermined component each period Nei, such as
The standard deviation or steady-state error that gynophore calculates any time period are within the scope of given threshold, it is determined that any time period
Interior rotary speed data relevant to the predetermined component is within the scope of setting speed, if the standard calculated for any time period
Difference or steady-state error are not within the scope of given threshold, it is determined that in any time period to relevant turn of the predetermined component
Fast data are not within the scope of setting speed.Here, given threshold range can be set according to the experience of those skilled in the art
It is fixed, the wide precision that can reduce abnormal vibrations identification of given threshold range setting, narrow can subtract of given threshold range setting
The data volume of the rotary speed data calculated less for subsequent analysis, in some instances it may even be possible to lead to that the revolving speed number for analytical calculation is not present
According to therefore, the setting for the given threshold range should can constrain revolving speed within limits, after ensuring that screening
Has the sufficient rotary speed data calculated for subsequent analysis.
In step s 30, the frequency values in each acceleration frequency spectrum for abnormal vibrations analysis are determined.
Here, the abscissa for each acceleration frequency spectrum being converted to can be frequency values, and ordinate can be Frequency and Amplitude value,
It can be chosen compared with Frequency and Amplitude threshold value for abnormal vibrations based on Frequency and Amplitude value in exemplary embodiments of the present invention
The frequency values of analysis.
A kind of situation, determine in each acceleration frequency spectrum for abnormal vibrations analysis frequency values the step of can include: look into
Frequency range value in the acceleration frequency spectrum is looked for be greater than the Frequency point of Frequency and Amplitude threshold value, by frequency corresponding with the Frequency point found
Rate value is as the frequency values for being used for abnormal vibrations analysis in the acceleration frequency spectrum.At this point, can be determined from an acceleration frequency spectrum
Frequency values of the one or more for abnormal vibrations analysis.
Another situation, determine in each acceleration frequency spectrum for abnormal vibrations analysis frequency values the step of can include:
Determine whether Frequency and Amplitude value corresponding with default concern Frequency point is greater than Frequency and Amplitude threshold value in the acceleration frequency spectrum, if with
The default corresponding Frequency and Amplitude value of Frequency point of paying close attention to is greater than Frequency and Amplitude threshold value, then will frequency corresponding with default concern Frequency point
It is worth as the frequency values for being used for abnormal vibrations analysis in the acceleration frequency spectrum.Above-mentioned method of determination is only for default concern Frequency point
Judged, determines that the mode of the frequency values for abnormal vibrations analysis is more accurate and more efficient relative to the first.
As an example, default concern Frequency point can for will include in acceleration frequency spectrum all Frequency points according to Frequency and Amplitude
The size descending of value arranges, the Frequency point before predetermined quantity.For example, being generally believed that frequency range value is most in acceleration frequency spectrum
The corresponding frequency values of big Frequency point are the fundamental frequency of predetermined component, and Frequency and Amplitude is worth the corresponding frequency values of second largest Frequency point and is
2 frequencys multiplication of predetermined component, and so on.At this point, carrying out analysis for default concern Frequency point is the equal of for predetermined component
Fundamental frequency and its frequency multiplication analyzed, can be improved abnormal vibrations identification accuracy.
In step s 40, based on the frequency values for being used for abnormal vibrations analysis in determining each acceleration frequency spectrum and often
Rotary speed data relevant to predetermined component determines predetermined component with the presence or absence of abnormal vibrations in a period.
Particularly, it may be determined that frequency values and each period in each acceleration frequency spectrum for abnormal vibrations analysis
Whether meet the predetermined linear regularity of distribution between interior rotary speed data relevant with predetermined component, when the frequency values and the revolving speed
When meeting the predetermined linear regularity of distribution between data, determine that there are abnormal vibrations for predetermined component.When the frequency values and described turn
When being unsatisfactory for the predetermined linear regularity of distribution between fast data, determining predetermined component, there is no abnormal vibrations.Here, predetermined linear point
Cloth rule can be the regularity of distribution of the linear relationship between the revolving speed and frequency for embodying predetermined component.
Preferably, determine in each acceleration frequency spectrum for abnormal vibrations analysis frequency values and each period in
The step of whether meeting the predetermined linear regularity of distribution between the relevant rotary speed data of predetermined component can include: be based on each period
Interior rotary speed data relevant to predetermined component calculates the revolving speed statistical value of the reflection data characteristics of each period, determines each add
Whether meet the predetermined linear regularity of distribution between frequency values and the revolving speed statistical value in speed spectrum for abnormal vibrations analysis.
As an example, the revolving speed statistical value of the reflection data characteristics of each period may include any one of following item:
The average value of rotary speed data relevant with the predetermined component, the period are interior relevant with the predetermined component in the period
The median of rotary speed data, in the period rotary speed data relevant to the predetermined component virtual value.As an example, revolving speed
The virtual value of data can refer to the maximum value of rotary speed data relevant to the predetermined component period Nei withRatio.
Referring to Fig. 2 by taking rotary speed data relevant to predetermined component is revolving speed statistical value as an example, introduces and determine each add
Whether meet the predetermined linear regularity of distribution between frequency values and revolving speed statistical value in speed spectrum for abnormal vibrations analysis
Step.
The step of whether meeting the predetermined linear regularity of distribution Fig. 2 shows determining scatterplot according to an exemplary embodiment of the present invention
Flow chart.
In step s 201, based on the frequency values and each period for being used for abnormal vibrations analysis in each acceleration frequency spectrum
Interior revolving speed statistical value draws revolving speed-frequency scatter plot.Here, a scatterplot in revolving speed-frequency scatter plot can correspond to one
A frequency in the revolving speed statistical value of period and acceleration frequency spectrum corresponding with the period for abnormal vibrations analysis
Value.
In step S202, the scatterplot in the preset range in preset frequency linear model is chosen.
Preferably, all scatterplots for including in revolving speed-frequency scatter plot can be screened, that is, chooses and is in preset frequency
Scatterplot in the preset range of linear model carries out subsequent abnormal vibrations analysis.Here, those skilled in the art can basis
It is actually needed to define the size of preset range, is carried out using the scatterplot in the preset range in preset frequency linear model
Abnormal vibrations analyze the accuracy that identification can be improved.
For example, preset frequency linear model can be the model for carrying out abnormal vibrations analysis for predetermined component, also
It is to say, preset frequency linear model can be the model that can embody the linear relationship between the revolving speed and frequency of predetermined component.With
Predetermined component is for wind-driven generator, preset frequency linear model can be for for reflecting the fundamental wave frequency of wind-driven generator at this time
The model of linear relationship between rate/frequency multiplication and the revolving speed of wind-driven generator, for example, preset frequency linear model is represented by fn
=n × p × r/60, here, fnFor the fundamental frequency of wind-driven generator or the frequency multiplication of fundamental frequency, p is wind-driven generator magnetic pole
Logarithm, r are wind-driven generator revolving speed (rotor speed), and n is the integer more than or equal to 1.As n=1, f1Indicate wind-driven generator
Fundamental frequency, as n >=2, fnIndicate the frequency multiplication of the fundamental frequency of wind-driven generator.
The case where for wind-driven generator fundamental frequency (that is, n=1), preset range can refer to by scatterplot boundaryOrIn the region for being included.Similarly, for the frequency multiplication of wind-driven generator fundamental frequency
The case where (that is, n >=2), preset range can refer to by scatterplot boundaryOrIncluded
Region in.
As an example, can be based on the corresponding grid-connected range of speeds of wind power generating set, setting speed range (or given threshold
Range) and wind-driven generator number of magnetic pole pairs determine the value range of parameter b and parameter k in above-mentioned scatterplot boundary.
In step S203, the frequency in each acceleration frequency spectrum for abnormal vibrations analysis is determined based on the scatterplot of selection
Whether the predetermined linear regularity of distribution of preset frequency linear model is met between revolving speed statistical value in rate value and each period.
That is, the predetermined linear regularity of distribution can be determined by frequency linearity model, when the frequency values and the revolving speed count
When meeting the predetermined linear regularity of distribution of preset frequency linear model between value, it is linear with preset frequency to determine that predetermined component exists
The corresponding abnormal vibrations of model.
For example, when turn that preset frequency linear model is the fundamental frequency and wind-driven generator for reflecting wind-driven generator
When the model of the linear relationship between speed, if it is linear to meet the preset frequency between the frequency values and the revolving speed statistical value
The predetermined linear regularity of distribution of model, it is determined that there are fundamental frequency abnormal vibrations for wind-driven generator.When the linear mould of preset frequency
Type is the model of the linear relationship for reflecting between the frequency multiplication of the fundamental frequency of wind-driven generator and the revolving speed of wind-driven generator
When, if meeting the predetermined linear distribution rule of the preset frequency linear model between the frequency values and the revolving speed statistical value
Rule, it is determined that there are the abnormal vibrations of the frequency multiplication of fundamental frequency for wind-driven generator.
It should be understood that whether meeting the predetermined linear regularity of distribution between determining frequency values shown in Fig. 2 and revolving speed statistical value
Mode is merely illustrative, and those skilled in the art determine can be used other modes.
Revolving speed according to an exemplary embodiment of the present invention-frequency scatter plot exemplary diagram is shown respectively in Fig. 3 A to Fig. 3 D.
Fig. 3 A and Fig. 3 B are shown respectively when predetermined component is wind-driven generator, and the vibration acceleration data of predetermined component are the
Revolving speed-frequency scatter plot when the vibration acceleration data of one predetermined direction and the second predetermined direction, abscissa are revolving speed statistics
It is worth (such as mean speed), ordinate is frequency values, and curve 1 indicates the preset frequency linear model of wind-driven generator fundamental frequency,
Curve 2 indicates the preset frequency linear model of the frequency multiplication (2 frequency multiplication) of wind-driven generator fundamental frequency.By taking curve 1 as an example, when it is desirable that
When identifying that wind-driven generator whether there is the abnormal vibrations of fundamental frequency, the scatterplot in the preset range around curve 1 can be chosen,
Determine it is linear whether the scatterplot of selection meets preset frequency by the method for determination of linear regression parameters estimation or objective function
The predetermined linear regularity of distribution of model, if the scatterplot chosen meets the predetermined linear regularity of distribution of preset frequency linear model,
Then showing wind-driven generator, there are the abnormal vibrations of fundamental frequency.
Fig. 3 C and Fig. 3 D are shown respectively when predetermined component is wind-driven generator, and the vibration acceleration data of predetermined component are the
Revolving speed-frequency scatter plot when the vibration acceleration data of one predetermined direction or the second predetermined direction, abscissa are revolving speed statistics
It is worth (such as mean speed), ordinate is frequency values, and curve 1 indicates the preset frequency linear model of wind-driven generator fundamental frequency,
Curve 2 indicates that the preset frequency linear model of 2 frequencys multiplication of wind-driven generator fundamental frequency, curve 3 indicate wind-driven generator fundamental wave
The preset frequency linear model of 3 frequencys multiplication of frequency.Identification side shown in the mode and Fig. 3 A and Fig. 3 B of specific identification abnormal vibrations
Formula is identical, and the present invention repeats no more this partial content.
It preferably, can be by for reflecting the preset frequency linear model of the predetermined linear regularity of distribution or by being used to indicate
The objective function of distance of each scatterplot (being formed by scatterplot by frequency values and rotary speed data) to preset frequency linear model Lai
It determines interior relevant to predetermined component for the frequency values of abnormal vibrations analysis and each period in each acceleration frequency spectrum
Whether the predetermined linear regularity of distribution is met between rotary speed data.
It is introduced below with reference to Fig. 4 based on preset frequency linear model and determines whether scatterplot meets predetermined linear distribution rule
The step of rule.
Fig. 4, which is shown, according to an exemplary embodiment of the present invention determines whether scatterplot meets based on preset frequency linear model
The flow chart of the step of predetermined linear regularity of distribution.
As shown in figure 4, in step S401, to the frequency values corresponding with the scatterplot of selection and the rotary speed data into
Row linear regression obtains the model parameter of linear regression.As an example, least square method or maximum likelihood method pair can be used and dissipate
The corresponding frequency values of point and rotary speed data carry out linear regression, obtain linear regression model (LRM), and then determine linear regression model (LRM)
Model parameter.However the invention is not limited thereto, and other modes also can be used and carry out linear regression analysis to scatterplot.
In step S402, the difference of the specified parameter of computation model parameter and predetermined component, and judge that the difference is
It is no to be greater than the first setting value.
As an example, the specified parameter can be the number of magnetic pole pairs of wind-driven generator when predetermined component is wind-driven generator,
When predetermined component is generator tooth socket, which can be the number of magnetic pole pairs of wind-driven generator, when predetermined component is gear
When case, which can be the number of teeth of power gear or the number of teeth of driven gear, and when predetermined component is rolling bearing, this refers to
Determining parameter can be the failure coefficient of rolling bearing (e.g., corresponding to one of the inner ring of rolling bearing, outer ring, rolling element, retainer
Failure coefficient).
For example, by taking predetermined component is wind-driven generator as an example, to frequency values corresponding with the scatterplot of selection and rotary speed data
Linear regression is carried out, to obtain model parameterAgain by model parameterIt is poor to make with the number of magnetic pole pairs p of wind-driven generator.
If the difference of model parameter and the specified parameter of predetermined component is not more than the first setting value, then follow the steps
S403: it determines and meets preset frequency linear model between the corresponding frequency values of scatterplot and the rotary speed data chosen
The predetermined linear regularity of distribution.
If the difference of model parameter and the specified parameter of predetermined component is greater than the first setting value, S404 is thened follow the steps:
It determines and is unsatisfactory for the preset frequency linear model between the corresponding frequency values of scatterplot and the rotary speed data chosen
The predetermined linear regularity of distribution.
The step for determining whether scatterplot meet the predetermined linear regularity of distribution is introduced based on objective function below with reference to Fig. 5
Suddenly.
Fig. 5 shows that according to an exemplary embodiment of the present invention to determine whether scatterplot meets based on objective function predetermined linear
The flow chart of the step of regularity of distribution.
As shown in figure 5, establishing objective function in step S501.Here, objective function can indicate each scatterplot to predetermined
The distance of frequency linearity model.Preferably, objective function can indicate each scatterplot to preset frequency linear model distance it is equal
Root value.
For example, objective function can be indicated by following formula:
In formula (1), y indicates that objective function, N are the quantity of scatterplot, fi' it is i-th of scatterplot in revolving speed-frequency scatter plot
In ordinate value, ri' be and fiAbscissa value in ' corresponding revolving speed-frequency scatter plot, that is, corresponding turn of i-th of scatterplot
Fast statistical value.
In step S502, by the way that frequency values corresponding with the scatterplot of selection and rotary speed data are substituted into objective function, obtain
Obtain the value of objective function.For example, can be by above-mentioned formula (1) come the value of calculating target function y.
In step S503, judge whether the value of objective function is greater than the second setting value.
If the value of objective function is not more than the second setting value, S504 is thened follow the steps: determining corresponding with the scatterplot chosen
The frequency values and the rotary speed data between meet the predetermined linear regularity of distribution of preset frequency linear model.
If the value of objective function is greater than the second setting value, S505 is thened follow the steps: determining corresponding with the scatterplot chosen
The predetermined linear regularity of distribution of preset frequency linear model is unsatisfactory between the frequency values and the rotary speed data.
Preferably, it can determine that above-mentioned first sets based on scatterplot boundary (that is, preset range of preset frequency linear model)
The value range of definite value and the second setting value.
Fig. 6 shows the structure chart of the equipment of identification abnormal vibrations according to an exemplary embodiment of the present invention.
As shown in fig. 6, the equipment of identification abnormal vibrations according to an exemplary embodiment of the present invention includes that operation data obtains
Module 10, time-frequency convert module 20, frequency values determining module 30 and abnormal vibrations analysis module 40.
Particularly, operation data obtains module 10, for obtaining the predetermined component of wind power generating set in multiple times
Operation data in section.Here, the operation data may include predetermined component vibration acceleration data and with predetermined component phase
The rotary speed data of pass.
Preferably, the vibration acceleration data of predetermined component may include the vibration acceleration data and of the first predetermined direction
The vibration acceleration data of two predetermined directions.As an example, the first predetermined direction can refer to from the head of wind power generating set to tail
The direction in portion, the second predetermined direction can refer to the direction vertical with wind direction (for example, field personnel station is on the leeward, towards machine
Head, the left and right directions of field personnel may be defined as the second predetermined direction).
In exemplary embodiments of the present invention, can be made a reservation for by the vibration acceleration data and second to the first predetermined direction
The vibration acceleration data in direction are respectively processed, and whether there is abnormal vibration to judge predetermined component in the first predetermined direction causes
Dynamic or predetermined component whether there is abnormal vibrations on the second predetermined direction.
Here, the type of abnormal vibrations may include that the fundamental vibration of predetermined component is abnormal different with the multiple-frequency vibration of predetermined component
Often.As an example, the predetermined component may include any one of following item: wind-driven generator, generator tooth socket, gear-box,
Rolling bearing.
In the first embodiment, predetermined component can be the wind-driven generator of wind power generating set.In the case, reservations
The vibration acceleration data of part can be the vibration acceleration data of engine rooms of wind power generators, rotary speed data relevant to predetermined component
It can be the rotary speed data of wind-driven generator.Correspondingly, the type of abnormal vibrations may include the fundamental frequency vibration of wind-driven generator
Abnormal and fundamental frequency multiple-frequency vibration is abnormal.
In a second embodiment, predetermined component can be generator tooth socket (such as generator unit stator tooth socket).In the case, in advance
The vibration acceleration data for determining component can be the vibration acceleration data of engine rooms of wind power generators, revolving speed relevant to predetermined component
Data can be the rotary speed data of wind-driven generator.Correspondingly, the type of abnormal vibrations may include that the vibration of generator tooth socket frequency is different
Often and/or the multiple-frequency vibration of generator tooth socket frequency is abnormal.Here, since generator tooth socket frequency is equal to determining for wind-driven generator
Sub- tooth socket number × wind-driven generator revolving speed/60, and the stator tooth slot number of wind-driven generator is greater than the number of magnetic pole pairs of wind-driven generator,
Therefore, generator tooth socket frequency (or frequency multiplication) is greater than the fundamental frequency (or frequency multiplication) of wind-driven generator.That is, can be by dividing
The data of each acceleration frequency spectrum medium-high frequency part (part corresponding with generator tooth socket frequency and/or frequency multiplication) are analysed to determine wind
Power generator tooth socket is with the presence or absence of abnormal vibrations (with the presence or absence of generator tooth socket frequency abnormal vibration and/or generator tooth socket frequency
The multiple-frequency vibration of rate is abnormal), by analyze low frequency part in each acceleration frequency spectrum (with the fundamental frequency of wind-driven generator and/
Or the corresponding part of frequency multiplication) data come determine wind-driven generator with the presence or absence of abnormal vibrations (with the presence or absence of wind-driven generator
Fundamental frequency abnormal vibration and/or the multiple-frequency vibration of fundamental frequency are abnormal).
In the third embodiment, predetermined component can be gear-box, include power gear and driven gear in gear-box.Herein
In the case of, the vibration acceleration data of predetermined component can be the vibration acceleration data or driven tooth of the power gear in gear-box
The vibration acceleration data of wheel, rotary speed data relevant to predetermined component can be the rotary speed data or driven of axis where power gear
The rotary speed data of axis where gear.Correspondingly, the type of abnormal vibrations may include power gear meshing frequency abnormal vibration, power
Multiple-frequency vibration exception, times frequency vibration of driven gear meshing frequency abnormal vibration and driven gear meshing frequency of gear mesh frequency
It is dynamic abnormal.Here, the meshing frequency of power gear (or driven gear) is equal to the number of teeth (or the tooth of driven gear of power gear
Number) revolving speed/60 of axis where × power gear (or driven gear).
In the fourth embodiment, predetermined component can be rolling bearing, can refer to multiple rollings in wind power generating set here
Any rolling bearing in bearing.In the case, the vibration acceleration data of predetermined component can be the bearing block of rolling bearing
Vibration acceleration data, rotary speed data relevant to predetermined component can be rolling bearing rotary speed data.Correspondingly, abnormal vibration
Dynamic type may include the multiple-frequency vibration of Rolling Bearing Fault Character frequency abnormal vibration and/or Rolling Bearing Fault Character frequency
It is abnormal.Here, it should be appreciated that rolling bearing may include bearing inner race, outer ring, rolling element, retainer, and correspondingly, predetermined component can
Refer to bearing inner race, outer ring, rolling element or retainer.For example, can know at this time so that predetermined component is the rolling element of rolling bearing as an example
Other rolling element is abnormal with the presence or absence of the fault characteristic frequency abnormal vibrations of rolling element or the multiple-frequency vibration of fault characteristic frequency.This
In, bearing inner race, outer ring, rolling element, retainer respectively correspond respective failure coefficient, pass through failure coefficient and the product for turning frequency
It can get fault characteristic frequency.
Time-frequency convert module 20, for carrying out frequency respectively to vibration acceleration data of the predetermined component within multiple periods
Domain conversion, to obtain corresponding multiple acceleration frequency spectrums.
For example, time-frequency convert module 20, can be used for predetermined component in any time period in multiple periods
Vibration acceleration data carry out frequency domain conversion to obtain acceleration corresponding with the vibration acceleration data in any time period
Spend frequency spectrum, that is, the corresponding acceleration frequency spectrum of vibration acceleration data in a period.As an example, time-frequency convert mould
Block 20 can carry out frequency domain conversion to vibration acceleration data of the predetermined component within multiple periods by Fast Fourier Transform (FFT),
However the invention is not limited thereto, other modes can be used also to carry out frequency domain conversion.
Preferably, time-frequency convert module 20 can be used for based on rotary speed data relevant to predetermined component in advance to acquisition
Operation data of the predetermined component within multiple periods screened, by the predetermined component in the operation data after screening more
Vibration acceleration data in a period carry out frequency domain conversion.
For example, time-frequency convert module 20, revolving speed relevant to the predetermined component was determined in each period
Whether data are within the scope of setting speed, if rotary speed data relevant with the predetermined component is in setting in any time period
Determine in the range of speeds, then the vibration acceleration data to predetermined component in any time period carry out frequency domain conversion, to obtain
Obtain a corresponding acceleration frequency spectrum.
As an example, time-frequency convert module 20, can be used for calculating revolving speed relevant to predetermined component in each period
The standard deviation or steady-state error of data, if the standard deviation or steady-state error for any time period calculating are in given threshold model
In enclosing, it is determined that rotary speed data relevant to the predetermined component is within the scope of setting speed in any time period, such as
The standard deviation or steady-state error that gynophore calculates any time period are not within the scope of given threshold, it is determined that any time
Rotary speed data relevant to the predetermined component is not within the scope of setting speed in section.Here, given threshold range can basis
The experience of those skilled in the art is set.
Frequency values determining module 30, for determining the frequency values in each acceleration frequency spectrum for abnormal vibrations analysis.
Here, the abscissa for each acceleration frequency spectrum being converted to can be frequency values, and ordinate can be Frequency and Amplitude value,
In exemplary embodiments of the present invention frequency values determining module 30 can based on Frequency and Amplitude value compared with Frequency and Amplitude threshold value come
Choose the frequency values for abnormal vibrations analysis.
A kind of situation, frequency values determining module 30 can be used for searching frequency range value in the acceleration frequency spectrum and be greater than frequency
The Frequency point of rate amplitude threshold, using frequency values corresponding with the Frequency point found as being used for abnormal vibration in the acceleration frequency spectrum
The frequency values of dynamic analysis.At this point, one or more frequencies for abnormal vibrations analysis can be determined from an acceleration frequency spectrum
Value.
Another situation, frequency values determining module 30 are determined in the acceleration frequency spectrum and default concern frequency
Whether the corresponding Frequency and Amplitude value of point is greater than Frequency and Amplitude threshold value, if Frequency and Amplitude value corresponding with default concern Frequency point is big
In Frequency and Amplitude threshold value, then will frequency values corresponding with default concern Frequency point as being used for abnormal vibrations in the acceleration frequency spectrum
The frequency values of analysis.
As an example, default concern Frequency point can for include by acceleration frequency spectrum all Frequency points according to Frequency and Amplitude value
The arrangement of size descending, the Frequency point before predetermined quantity.For example, generally it is believed that frequency range value is maximum in acceleration frequency spectrum
The corresponding frequency values of Frequency point be predetermined component fundamental frequency, it is pre- that Frequency and Amplitude, which is worth the corresponding frequency values of second largest Frequency point,
Determine 2 frequencys multiplication of component, and so on.At this point, carrying out analysis for concern Frequency point is the equal of the fundamental frequency for predetermined component
And its frequency multiplication is analyzed, and the accuracy of abnormal vibrations identification can be improved.
Abnormal vibrations analysis module 40, for based on the frequency for being used for abnormal vibrations analysis in determining each acceleration frequency spectrum
Rotary speed data relevant to predetermined component determines predetermined component with the presence or absence of abnormal vibrations in rate value and each period.
Particularly, abnormal vibrations analysis module 40 is analyzed for determining in each acceleration frequency spectrum for abnormal vibrations
Frequency values and whether meet the predetermined linear regularity of distribution between rotary speed data relevant to predetermined component in each period,
When meeting the predetermined linear regularity of distribution between the frequency values and the rotary speed data, determine that predetermined component has abnormal vibration
It is dynamic.When being unsatisfactory for the predetermined linear regularity of distribution between the frequency values and the rotary speed data, determine that predetermined component is not present
Abnormal vibrations.Here, the predetermined linear regularity of distribution can be the linear relationship between the revolving speed and frequency for embodying predetermined component
The regularity of distribution.
Preferably, abnormal vibrations analysis module 40 can be used for based on relevant to predetermined component turn in each period
Fast data calculate the revolving speed statistical value of the reflection data characteristics of each period, determine in each acceleration frequency spectrum and shake for abnormal
Whether meet the predetermined linear regularity of distribution between the frequency values and the revolving speed statistical value of dynamic analysis.
As an example, the revolving speed statistical value of the reflection data characteristics of each period may include any one of following item:
The average value of rotary speed data relevant with the predetermined component, the period are interior relevant with the predetermined component in the period
The median of rotary speed data, in the period rotary speed data relevant to the predetermined component virtual value.Here, rotary speed data
Virtual value can refer to the maximum value of rotary speed data relevant to the predetermined component period Nei withRatio.
Referring to Fig. 7 by taking rotary speed data relevant to predetermined component is revolving speed statistical value as an example, introduces and determine each add
Whether meet the predetermined linear regularity of distribution between frequency values and revolving speed statistical value in speed spectrum for abnormal vibrations analysis
Process.
Fig. 7 shows the structure chart of abnormal vibrations analysis module according to an exemplary embodiment of the present invention.
As shown in fig. 7, abnormal vibrations analysis module 40 according to an exemplary embodiment of the present invention may include that scatter plot is drawn
Module 41, scatterplot screening module 42 and linear distribution determining module 43.
Particularly, scatter plot drafting module 41, for based in each acceleration frequency spectrum be used for abnormal vibrations analysis
Revolving speed statistical value in frequency values and each period draws revolving speed-frequency scatter plot.Here, in revolving speed-frequency scatter plot
One scatterplot can correspond in the revolving speed statistical value and acceleration frequency spectrum corresponding with the period of a period for abnormal
One frequency values of vibration analysis.
Scatterplot screening module 42, for choosing the scatterplot being in the preset range of preset frequency linear model.
Preferably, scatterplot screening module 42 can be used for sieving all scatterplots for including in revolving speed-frequency scatter plot
Choosing chooses the scatterplot in the preset range in preset frequency linear model to carry out subsequent abnormal vibrations analysis.Here,
Those skilled in the art can define the size of preset range according to actual needs, using pre- in preset frequency linear model
If the accuracy of identification can be improved to carry out abnormal vibrations analysis in the scatterplot in range.
Linear distribution determining module 43 is shaken for being determined in each acceleration frequency spectrum based on the scatterplot of selection for abnormal
Whether the preset lines of preset frequency linear model are met between revolving speed statistical value in the frequency values and each period of dynamic analysis
The property regularity of distribution.That is, the predetermined linear regularity of distribution can be determined by frequency linearity model, when the frequency values and institute
When stating the predetermined linear regularity of distribution for meeting preset frequency linear model between revolving speed statistical value, linear distribution determining module 43,
For determining that predetermined component has abnormal vibrations corresponding with preset frequency linear model.
The process for determining whether scatterplot meets the predetermined linear regularity of distribution is described below based on preset frequency linear model.
Linear distribution determining module 43, for the frequency values corresponding with the scatterplot of selection and the rotary speed data into
Row linear regression obtains the model parameter of linear regression, the difference of computation model parameter and the specified parameter of predetermined component, and sentences
Whether the difference of breaking is greater than the first setting value, if the difference of model parameter and the specified parameter of predetermined component is not more than first
Setting value, it is determined that meet the linear mould of preset frequency between the frequency values corresponding with the scatterplot of selection and the rotary speed data
The predetermined linear regularity of distribution of type.If the difference of model parameter and the specified parameter of predetermined component is greater than the first setting value,
It determines and is unsatisfactory for the preset frequency linear model between the corresponding frequency values of scatterplot and the rotary speed data chosen
The predetermined linear regularity of distribution.
As an example, the specified parameter can be the number of magnetic pole pairs of wind-driven generator when predetermined component is wind-driven generator,
When predetermined component is generator tooth socket, which can be the number of magnetic pole pairs of wind-driven generator, when predetermined component is gear
When case, which can be the number of teeth of power gear or the number of teeth of driven gear, and when predetermined component is rolling bearing, this refers to
Determining parameter can be the failure coefficient of rolling bearing (e.g., corresponding to one of the inner ring of rolling bearing, outer ring, rolling element, retainer
Failure coefficient).
The mistake for determining whether scatterplot meet the predetermined linear regularity of distribution is introduced based on objective function below with reference to Fig. 8
Journey.
Fig. 8 shows the structure chart of linear distribution determining module according to an exemplary embodiment of the present invention.
As shown in figure 8, linear distribution determining module 43 according to an exemplary embodiment of the present invention may include that objective function is built
Vertical submodule 431, target function value computational submodule 432 and the regularity of distribution determine submodule 433.
Particularly, objective function setting up submodule 431, for establishing objective function.Here, objective function can indicate often
A scatterplot to preset frequency linear model distance.Preferably, objective function can indicate each scatterplot to the linear mould of preset frequency
The root-mean-square value of the distance of type.
Target function value computational submodule 432, for by by frequency values corresponding with the scatterplot of selection and rotary speed data
Objective function is substituted into, the value of objective function is obtained.
The regularity of distribution determines submodule 433, for judging whether the value of objective function is greater than the second setting value.If target
The value of function is not more than the second setting value, then the regularity of distribution determines that submodule 433 determines the frequency corresponding with the scatterplot chosen
Meet the predetermined linear regularity of distribution of preset frequency linear model between rate value and the rotary speed data.If the value of objective function
Greater than the second setting value, then the regularity of distribution determines submodule 433, for determine with the corresponding frequency values of scatterplot chosen and
The predetermined linear regularity of distribution of preset frequency linear model is unsatisfactory between the rotary speed data.Preferably, scatterplot side can be based on
Boundary's (that is, preset range of preset frequency linear model) determines the value range of above-mentioned first setting value and the second setting value.
An exemplary embodiment of the present invention also provides a kind of computing device.The computing device includes processor and storage
Device.Memory is for storing computer program.The computer program is executed by processor so that processor executes as described above
Identification abnormal vibrations method computer program.
An exemplary embodiment of the present invention also provides a kind of computer-readable storage medium for being stored with computer program
Matter.The computer-readable recording medium storage has makes processor execute above-mentioned identification abnormal vibrations when being executed by a processor
The computer program of method.The computer readable recording medium is the arbitrary data that can store the data read by computer system
Storage device.The example of computer readable recording medium includes: read-only memory, random access memory, CD-ROM, magnetic
Band, floppy disk, optical data storage devices and carrier wave (such as being transmitted through wired or wireless transmission path by the data of internet).
It, being capable of positioning wind promptly and accurately using the method and apparatus of the identification abnormal vibrations of exemplary embodiment of the present
There are the component of frequency anomaly vibration in power generator group, provided quickly and effectively to assess the vibrational state of predetermined component
The support of power.Further, it is also possible to effectively improve the efficiency of unit fault location, save O&M cost.
Although being particularly shown and describing the present invention, those skilled in the art referring to its exemplary embodiment
It should be understood that in the case where not departing from the spirit and scope of the present invention defined by claim form can be carried out to it
With the various changes in details.
Claims (30)
1. a kind of method for identifying abnormal vibrations, which is characterized in that the described method includes:
Operation data of the predetermined component of wind power generating set within multiple periods is obtained, the operation data includes described pre-
Determine component vibration acceleration data and rotary speed data relevant to the predetermined component;
Frequency domain conversion is carried out respectively to vibration acceleration data of the predetermined component within the multiple period, to be divided
Not corresponding multiple acceleration frequency spectrums;
Determine the frequency values in each acceleration frequency spectrum for abnormal vibrations analysis;
Based in determining each acceleration frequency spectrum be used for abnormal vibrations analysis frequency values and in each period with it is described
The relevant rotary speed data of predetermined component determines the predetermined component with the presence or absence of abnormal vibrations.
2. the method as described in claim 1, which is characterized in that the type of abnormal vibrations includes the fundamental frequency vibration of the predetermined component
Dynamic abnormal and the predetermined component multiple-frequency vibration is abnormal.
3. method according to claim 1 or 2, which is characterized in that shake in each acceleration frequency spectrum of determination for abnormal
The step of frequency values of dynamic analysis includes:
Search the Frequency point that frequency range value in the acceleration frequency spectrum is greater than Frequency and Amplitude threshold value;
Using frequency values corresponding with the Frequency point found as the frequency values for being used for abnormal vibrations analysis in the acceleration frequency spectrum.
4. method according to claim 1 or 2, which is characterized in that shake in each acceleration frequency spectrum of determination for abnormal
The step of frequency values of dynamic analysis includes:
Determine whether Frequency and Amplitude value corresponding with default concern Frequency point is greater than Frequency and Amplitude threshold value in the acceleration frequency spectrum;
If Frequency and Amplitude value corresponding with the default concern Frequency point is greater than Frequency and Amplitude threshold value, will be with the default pass
The corresponding frequency values of Frequency point are infused as the frequency values for being used for abnormal vibrations analysis in the acceleration frequency spectrum.
5. method as claimed in claim 4, which is characterized in that the default concern Frequency point is to include in acceleration frequency spectrum
All Frequency points arranged according to the size descending of Frequency and Amplitude value, the Frequency point before predetermined quantity.
6. the method as described in claim 1, which is characterized in that described abnormal based on being used in determining each acceleration frequency spectrum
Rotary speed data relevant to the predetermined component determines the predetermined component in the frequency values of vibration analysis and each period
Include: with the presence or absence of the step of abnormal vibrations
Determine whether meet the predetermined linear regularity of distribution between the frequency values and the rotary speed data;
When meeting the predetermined linear regularity of distribution between the frequency values and the rotary speed data, the predetermined component is determined
There are abnormal vibrations.
7. method as claimed in claim 6, which is characterized in that be between the determination frequency values and the rotary speed data
No the step of meeting the predetermined linear regularity of distribution includes:
Revolving speed-frequency scatter plot is drawn based on the frequency values and the rotary speed data, wherein the revolving speed-frequency scatter plot
In an a scatterplot corresponding period rotary speed data and acceleration frequency spectrum corresponding with the period in for abnormal
One frequency values of vibration analysis;
Choose the scatterplot in the preset range in preset frequency linear model;
Determine whether meet the linear mould of the preset frequency between the frequency values and the rotary speed data based on the scatterplot of selection
The predetermined linear regularity of distribution of type.
8. the method for claim 7, which is characterized in that the scatterplot based on selection determine the frequency values with it is described
The step of whether meeting the predetermined linear regularity of distribution of the preset frequency linear model between rotary speed data include:
Linear regression is carried out to the frequency values corresponding with the scatterplot of selection and the rotary speed data, obtains the mould of linear regression
Shape parameter;
Calculate the difference of the specified parameter of the model parameter and the predetermined component;
When the difference is not more than the first setting value, the frequency values corresponding with the scatterplot chosen and the revolving speed number are determined
Meet the predetermined linear regularity of distribution of the preset frequency linear model between.
9. the method for claim 7, which is characterized in that the scatterplot based on selection determine the frequency values with it is described
The step of whether meeting the predetermined linear regularity of distribution of the preset frequency linear model between rotary speed data include:
Objective function is established, the objective function indicates the distance of each scatterplot to the preset frequency linear model;
By the way that the frequency values corresponding with the scatterplot of selection and the rotary speed data are substituted into objective function, the target is obtained
The value of function;
When the value of the objective function is not more than setting value, the frequency values corresponding with the scatterplot chosen and described turn are determined
Meet the predetermined linear regularity of distribution of the preset frequency linear model between fast data.
10. the method as described in any one of claim 6-9, which is characterized in that the determination frequency values with it is described
The step of whether meeting the predetermined linear regularity of distribution between rotary speed data include:
The revolving speed statistical value of the reflection data characteristics of each period is calculated based on the rotary speed data;
It determines and whether meets the predetermined linear regularity of distribution between the frequency values and the revolving speed statistical value.
11. method as claimed in claim 10, which is characterized in that the revolving speed of the reflection data characteristics of each period is united
Evaluation includes any one of following item: the average value of rotary speed data relevant to the predetermined component in the period, this when
Between the median of rotary speed data relevant with the predetermined component in section, the period interior revolving speed relevant with the predetermined component
The virtual value of data.
12. the method as described in claim 1, which is characterized in that it is described to the predetermined component within the multiple period
Vibration acceleration data carry out frequency domain conversion, include: the step of corresponding multiple acceleration frequency spectrums to obtain
Determine whether rotary speed data relevant to the predetermined component is within the scope of setting speed in each period;
If rotary speed data relevant to the predetermined component is within the scope of setting speed in any time period, to described pre-
Determine vibration acceleration data of the component in any time period and carry out frequency domain conversion, to obtain corresponding acceleration frequency
Spectrum.
13. method as claimed in claim 12, which is characterized in that in determination each period with the predetermined component phase
The step whether rotary speed data of pass is within the scope of setting speed includes:
Calculate the standard deviation or steady-state error of rotary speed data relevant to the predetermined component in each period;
If the standard deviation or steady-state error that calculate for any time period are within the scope of given threshold, it is determined that described any
Rotary speed data relevant to the predetermined component is within the scope of setting speed in period;
If the standard deviation or steady-state error that calculate for any time period are not within the scope of given threshold, it is determined that described
Rotary speed data relevant to the predetermined component is not within the scope of setting speed in one period.
14. method according to claim 1 or 2, which is characterized in that the predetermined component includes any one of following item:
Wind-driven generator, generator tooth socket, gear-box, rolling bearing,
Wherein, when the predetermined component is wind-driven generator, the vibration acceleration data of the predetermined component are wind-power electricity generation
The vibration acceleration data of machine cabin, rotary speed data relevant to the predetermined component is the rotary speed data of wind-driven generator, different
The type often vibrated includes the fundamental frequency abnormal vibration of wind-driven generator and the multiple-frequency vibration exception of fundamental frequency;
When the predetermined component is generator tooth socket, the vibration acceleration data of the predetermined component are engine rooms of wind power generators
Vibration acceleration data, rotary speed data relevant to the predetermined component is the rotary speed data of wind-driven generator, abnormal vibrations
Type include generator tooth socket frequency abnormal vibration and generator tooth socket frequency multiple-frequency vibration it is abnormal;
When the predetermined component is gear-box, the vibration acceleration data of the predetermined component are the power gear in gear-box
Or the vibration acceleration data of driven gear, rotary speed data relevant to the predetermined component is power gear or driven gear institute
In the rotary speed data of axis, the type of abnormal vibrations includes the multiple-frequency vibration of power gear meshing frequency and power gear meshing frequency
The multiple-frequency vibration of abnormal or driven gear meshing frequency and driven gear meshing frequency is abnormal;
When the predetermined component is rolling bearing, the vibration acceleration data of the predetermined component are the bearing block of rolling bearing
Vibration acceleration data, rotary speed data relevant to the predetermined component is the rotary speed data of rolling bearing, abnormal vibrations
Type includes the multiple-frequency vibration exception of Rolling Bearing Fault Character frequency abnormal vibration and Rolling Bearing Fault Character frequency.
15. a kind of equipment for identifying abnormal vibrations, which is characterized in that the equipment includes:
Operation data obtains module, for obtaining operation data of the predetermined component of wind power generating set within multiple periods,
The operation data include the predetermined component vibration acceleration data and rotary speed data relevant to the predetermined component;
Time-frequency convert module, for being carried out respectively to vibration acceleration data of the predetermined component within the multiple period
Frequency domain conversion, to obtain corresponding multiple acceleration frequency spectrums;
Frequency values determining module, for determining the frequency values in each acceleration frequency spectrum for abnormal vibrations analysis;
Abnormal vibrations analysis module, for based in determining each acceleration frequency spectrum be used for abnormal vibrations analysis frequency values with
And rotary speed data relevant with the predetermined component determines that the predetermined component whether there is abnormal vibrations in each period.
16. equipment as claimed in claim 15, which is characterized in that the type of abnormal vibrations includes the fundamental frequency of the predetermined component
The multiple-frequency vibration of abnormal vibration and the predetermined component is abnormal.
17. the equipment as described in claim 15 or 16, which is characterized in that the frequency values determining module is for searching the acceleration
Spend frequency spectrum in frequency range value be greater than Frequency and Amplitude threshold value Frequency point, using frequency values corresponding with the Frequency point found as
Frequency values in the acceleration frequency spectrum for abnormal vibrations analysis.
18. the equipment as described in claim 15 or 16, which is characterized in that the frequency values determining module is for determining the acceleration
Whether Frequency and Amplitude value corresponding with default concern Frequency point is greater than Frequency and Amplitude threshold value in degree frequency spectrum, if with the default pass
It infuses the corresponding Frequency and Amplitude value of Frequency point and is greater than Frequency and Amplitude threshold value, then frequency values determining module will be with the default concern frequency
The corresponding frequency values of point are as the frequency values for being used for abnormal vibrations analysis in the acceleration frequency spectrum.
19. equipment as claimed in claim 18, which is characterized in that the default concern Frequency point is that will wrap in acceleration frequency spectrum
All Frequency points contained are arranged according to the size descending of Frequency and Amplitude value, the Frequency point before predetermined quantity.
20. equipment as claimed in claim 15, which is characterized in that the abnormal vibrations analysis module, for determining the frequency
Whether the predetermined linear regularity of distribution is met between rate value and the rotary speed data, when between the frequency values and the rotary speed data
When meeting the predetermined linear regularity of distribution, determine that there are abnormal vibrations for the predetermined component.
21. equipment as claimed in claim 20, which is characterized in that the abnormal vibrations analysis module includes:
Scatter plot drafting module, for drawing revolving speed-frequency scatter plot based on the frequency values and the rotary speed data, wherein
The rotary speed data of a scatterplot in a revolving speed-frequency scatter plot corresponding period and it is corresponding with the period plus
A frequency values in speed spectrum for abnormal vibrations analysis;
Scatterplot screening module, for choosing the scatterplot being in the preset range of preset frequency linear model;
Whether linear distribution determining module determines full between the frequency values and the rotary speed data for the scatterplot based on selection
The predetermined linear regularity of distribution of the foot preset frequency linear model.
22. equipment as claimed in claim 21, which is characterized in that the linear distribution determining module, for selection
The corresponding frequency values of scatterplot and the rotary speed data carry out linear regression, obtain the model parameter of linear regression, calculate institute
The difference for stating the specified parameter of model parameter and the predetermined component, when the difference is not more than the first setting value, determine with
Meet the preset lines of the preset frequency linear model between the corresponding frequency values of the scatterplot of selection and the rotary speed data
The property regularity of distribution.
23. equipment as claimed in claim 21, which is characterized in that the linear distribution determining module includes:
Objective function setting up submodule, for establishing objective function, the objective function indicates each scatterplot to the predetermined frequency
The root-mean-square value of the distance of rate linear model;
Target function value computational submodule, for by by the frequency values corresponding with the scatterplot of selection and the rotary speed data
Objective function is substituted into, the value of the objective function is obtained;
The regularity of distribution determines submodule, the scatterplot for when the value of the objective function is not more than setting value, determining with choosing
Meet the predetermined linear regularity of distribution of the preset frequency linear model between the corresponding frequency values and the rotary speed data.
24. the equipment as described in claim 20-23, which is characterized in that the abnormal vibrations analysis module, for based on described
Rotary speed data calculates the revolving speed statistical value of the reflection data characteristics of each period, determines that the frequency values and the revolving speed count
Whether the predetermined linear regularity of distribution is met between value.
25. equipment as claimed in claim 24, which is characterized in that the revolving speed statistical value of the reflection data characteristics of each period
Including any one of following item: the average value of rotary speed data relevant to the predetermined component, the period in the period
The median of interior rotary speed data relevant to the predetermined component, rotary speed data relevant with the predetermined component in the period
Virtual value.
26. equipment as claimed in claim 15, which is characterized in that the time-frequency convert module, for determining each period
Whether interior rotary speed data relevant to the predetermined component is within the scope of setting speed, if in any time period with it is described pre-
Determine the relevant rotary speed data of component to be within the scope of setting speed, then the vibration to the predetermined component in any time period
Dynamic acceleration information carries out frequency domain conversion, to obtain a corresponding acceleration frequency spectrum.
27. equipment as claimed in claim 26, which is characterized in that the time-frequency convert module, for calculating each period
The standard deviation or steady-state error of interior rotary speed data relevant to the predetermined component, if the standard calculated for any time period
Difference or steady-state error be within the scope of given threshold, it is determined that the interior revolving speed relevant to the predetermined component of any time period
Data are within the scope of setting speed,
If the standard deviation or steady-state error that calculate for any time period are not within the scope of given threshold, it is determined that described
Rotary speed data relevant to the predetermined component is not within the scope of setting speed in one period.
28. the equipment as described in claim 15 or 16, which is characterized in that the predetermined component includes any in following item
: wind-driven generator, generator tooth socket, gear-box, rolling bearing,
Wherein, when the predetermined component is wind-driven generator, the vibration acceleration data of the predetermined component are wind-power electricity generation
The vibration acceleration data of machine cabin, rotary speed data relevant to the predetermined component is the rotary speed data of wind-driven generator, different
The type often vibrated includes the fundamental frequency abnormal vibration of wind-driven generator and the multiple-frequency vibration exception of fundamental frequency;
When the predetermined component is generator tooth socket, the vibration acceleration data of the predetermined component are engine rooms of wind power generators
Vibration acceleration data, rotary speed data relevant to the predetermined component is the rotary speed data of wind-driven generator, abnormal vibrations
Type include generator tooth socket frequency abnormal vibration and generator tooth socket frequency multiple-frequency vibration it is abnormal;
When the predetermined component is gear-box, the vibration acceleration data of the predetermined component are the power gear in gear-box
Or the vibration acceleration data of driven gear, rotary speed data relevant to the predetermined component is power gear or driven gear institute
In the rotary speed data of axis, the type of abnormal vibrations includes the multiple-frequency vibration of power gear meshing frequency and power gear meshing frequency
The multiple-frequency vibration of abnormal or driven gear meshing frequency and driven gear meshing frequency is abnormal;
When the predetermined component is rolling bearing, the vibration acceleration data of the predetermined component are the bearing block of rolling bearing
Vibration acceleration data, rotary speed data relevant to the predetermined component is the rotary speed data of rolling bearing, abnormal vibrations
Type includes the multiple-frequency vibration exception of Rolling Bearing Fault Character frequency abnormal vibration and Rolling Bearing Fault Character frequency.
29. a kind of computer readable storage medium for being stored with computer program, when the computer program is held by processor
The method of the identification abnormal vibrations as described in any one in claim 1-14 is realized when row.
30. a kind of computing device, which is characterized in that the computing device includes:
Processor;
Memory is stored with computer program, when the computer program is executed by processor, realizes such as claim 1-14
In any one described in identify abnormal vibrations method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711384895.3A CN109973325B (en) | 2017-12-20 | 2017-12-20 | Method and apparatus for identifying abnormal vibration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711384895.3A CN109973325B (en) | 2017-12-20 | 2017-12-20 | Method and apparatus for identifying abnormal vibration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109973325A true CN109973325A (en) | 2019-07-05 |
CN109973325B CN109973325B (en) | 2020-09-29 |
Family
ID=67072087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711384895.3A Active CN109973325B (en) | 2017-12-20 | 2017-12-20 | Method and apparatus for identifying abnormal vibration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109973325B (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110925135A (en) * | 2019-12-06 | 2020-03-27 | 中国船舶重工集团海装风电股份有限公司 | Wind turbine generator pitch motor selection method, device, equipment and storage medium |
CN110987404A (en) * | 2019-12-23 | 2020-04-10 | 安徽容知日新科技股份有限公司 | Rotating speed identification method and computing equipment |
CN110986264A (en) * | 2019-11-21 | 2020-04-10 | 珠海格力电器股份有限公司 | Air conditioner low-frequency resonance noise identification control method and air conditioner |
CN112985578A (en) * | 2019-12-13 | 2021-06-18 | 中车株洲电力机车研究所有限公司 | Wind turbine generator resonance diagnosis method and device |
WO2021119909A1 (en) * | 2019-12-16 | 2021-06-24 | 远景能源有限公司 | Method and system for monitoring health state of blade root fastener |
CN113432856A (en) * | 2021-06-28 | 2021-09-24 | 西门子电机(中国)有限公司 | Motor testing method, device, electronic equipment and storage medium |
CN114087136A (en) * | 2021-11-22 | 2022-02-25 | 西安热工研究院有限公司 | Wind turbine generator resonance diagnosis method |
CN114184879A (en) * | 2021-11-08 | 2022-03-15 | 陕西千山航空电子有限责任公司 | Abnormal frequency point monitoring method for helicopter vibration data |
CN114216640A (en) * | 2022-02-21 | 2022-03-22 | 蘑菇物联技术(深圳)有限公司 | Method, apparatus and medium for detecting fault status of industrial equipment |
CN114658611A (en) * | 2020-12-23 | 2022-06-24 | 新疆金风科技股份有限公司 | Method and device for detecting abnormality of main bearing of wind power generator |
CN116449688A (en) * | 2023-06-19 | 2023-07-18 | 深圳市大族机器人有限公司 | Robot motor PID parameter determination method, device and computer equipment |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1430063A (en) * | 2001-12-31 | 2003-07-16 | 岩壶卓三 | Diagnosing system for abnormity of rotary machine |
CN101532911A (en) * | 2009-04-24 | 2009-09-16 | 华北电力大学 | Large steam turbine-generator set rotor crack fault real-time diagnosis method |
EP2754889A1 (en) * | 2013-01-15 | 2014-07-16 | ABB Technology Ltd | Method and controller for damping vibrations in a wind power system |
CN105547462A (en) * | 2015-12-16 | 2016-05-04 | 中国航空工业集团公司沈阳发动机设计研究所 | Radial vibration measuring method of engine rotor |
CN105675112A (en) * | 2015-12-31 | 2016-06-15 | 北京金风科创风电设备有限公司 | Method and device for monitoring abnormal vibration of wind turbine generator |
CN107061186A (en) * | 2017-06-09 | 2017-08-18 | 北京金风慧能技术有限公司 | Vibration of wind generating set abnormity early warning method and apparatus |
JP2017181283A (en) * | 2016-03-30 | 2017-10-05 | Ntn株式会社 | Tooth number specification device of single pinion type epicyclic gear and tooth number specification method thereof |
-
2017
- 2017-12-20 CN CN201711384895.3A patent/CN109973325B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1430063A (en) * | 2001-12-31 | 2003-07-16 | 岩壶卓三 | Diagnosing system for abnormity of rotary machine |
CN101532911A (en) * | 2009-04-24 | 2009-09-16 | 华北电力大学 | Large steam turbine-generator set rotor crack fault real-time diagnosis method |
EP2754889A1 (en) * | 2013-01-15 | 2014-07-16 | ABB Technology Ltd | Method and controller for damping vibrations in a wind power system |
CN105547462A (en) * | 2015-12-16 | 2016-05-04 | 中国航空工业集团公司沈阳发动机设计研究所 | Radial vibration measuring method of engine rotor |
CN105675112A (en) * | 2015-12-31 | 2016-06-15 | 北京金风科创风电设备有限公司 | Method and device for monitoring abnormal vibration of wind turbine generator |
JP2017181283A (en) * | 2016-03-30 | 2017-10-05 | Ntn株式会社 | Tooth number specification device of single pinion type epicyclic gear and tooth number specification method thereof |
CN107061186A (en) * | 2017-06-09 | 2017-08-18 | 北京金风慧能技术有限公司 | Vibration of wind generating set abnormity early warning method and apparatus |
Non-Patent Citations (2)
Title |
---|
张学延: "《汽轮发电机组振动诊断》", 30 November 2008, 中国电力出版社 * |
罗良清等: "《统计学 第2版》", 31 January 2014, 湖南大学出版社 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110986264A (en) * | 2019-11-21 | 2020-04-10 | 珠海格力电器股份有限公司 | Air conditioner low-frequency resonance noise identification control method and air conditioner |
CN110925135A (en) * | 2019-12-06 | 2020-03-27 | 中国船舶重工集团海装风电股份有限公司 | Wind turbine generator pitch motor selection method, device, equipment and storage medium |
CN110925135B (en) * | 2019-12-06 | 2021-02-02 | 中国船舶重工集团海装风电股份有限公司 | Wind turbine generator pitch motor selection method, device, equipment and storage medium |
CN112985578B (en) * | 2019-12-13 | 2023-09-26 | 中车株洲电力机车研究所有限公司 | Wind turbine generator set resonance diagnosis method and device |
CN112985578A (en) * | 2019-12-13 | 2021-06-18 | 中车株洲电力机车研究所有限公司 | Wind turbine generator resonance diagnosis method and device |
CN113286944B (en) * | 2019-12-16 | 2022-06-10 | 远景能源有限公司 | Method and system for monitoring health state of blade root fastener |
CN113286944A (en) * | 2019-12-16 | 2021-08-20 | 远景能源有限公司 | Method and system for monitoring health state of blade root fastener |
WO2021119909A1 (en) * | 2019-12-16 | 2021-06-24 | 远景能源有限公司 | Method and system for monitoring health state of blade root fastener |
CN110987404B (en) * | 2019-12-23 | 2022-06-10 | 安徽容知日新科技股份有限公司 | Rotating speed identification method and computing equipment |
CN110987404A (en) * | 2019-12-23 | 2020-04-10 | 安徽容知日新科技股份有限公司 | Rotating speed identification method and computing equipment |
CN114658611A (en) * | 2020-12-23 | 2022-06-24 | 新疆金风科技股份有限公司 | Method and device for detecting abnormality of main bearing of wind power generator |
CN113432856A (en) * | 2021-06-28 | 2021-09-24 | 西门子电机(中国)有限公司 | Motor testing method, device, electronic equipment and storage medium |
CN114184879A (en) * | 2021-11-08 | 2022-03-15 | 陕西千山航空电子有限责任公司 | Abnormal frequency point monitoring method for helicopter vibration data |
CN114184879B (en) * | 2021-11-08 | 2024-01-23 | 陕西千山航空电子有限责任公司 | Abnormal frequency point monitoring method for helicopter vibration data |
CN114087136A (en) * | 2021-11-22 | 2022-02-25 | 西安热工研究院有限公司 | Wind turbine generator resonance diagnosis method |
CN114216640A (en) * | 2022-02-21 | 2022-03-22 | 蘑菇物联技术(深圳)有限公司 | Method, apparatus and medium for detecting fault status of industrial equipment |
CN116449688A (en) * | 2023-06-19 | 2023-07-18 | 深圳市大族机器人有限公司 | Robot motor PID parameter determination method, device and computer equipment |
CN116449688B (en) * | 2023-06-19 | 2023-10-31 | 深圳市大族机器人有限公司 | Robot motor PID parameter determination method, device and computer equipment |
Also Published As
Publication number | Publication date |
---|---|
CN109973325B (en) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109973325A (en) | Method and apparatus for identifying abnormal vibration | |
CN107061186B (en) | Vibration of wind generating set abnormity early warning method and apparatus | |
US20190310129A1 (en) | Detecting faults in rotor driven equipment | |
Bechhoefer et al. | Processing for improved spectral analysis | |
JP6783163B2 (en) | Wing vibration monitoring device and wing vibration monitoring method | |
Yoon et al. | Vibration‐based wind turbine planetary gearbox fault diagnosis using spectral averaging | |
Dong et al. | Vibration source identification of offshore wind turbine structure based on optimized spectral kurtosis and ensemble empirical mode decomposition | |
Haj Mohamad et al. | Rolling element bearing diagnostics using extended phase space topology | |
EP3788328B1 (en) | System and process for pattern matching bearing vibration diagnostics | |
Motte et al. | Operational modal analysis in the presence of harmonic excitations: a review | |
Mo et al. | Conditional empirical wavelet transform with modified ratio of cyclic content for bearing fault diagnosis | |
Antoni et al. | Time-angle periodically correlated processes | |
Teng et al. | Detection and quantization of bearing fault in direct drive wind turbine via comparative analysis | |
CN112345827A (en) | Graphical differentiation of spectral frequency families | |
CN108241761A (en) | Method and device for determining fatigue damage of generator component | |
Yan et al. | Degree of cyclic target protrusion defined on squared envelope spectrum for rotating machinery fault diagnosis | |
CN116558840B (en) | Method, device, equipment and storage medium for monitoring aero-engine blade | |
Samadani et al. | Feature extraction for bearing diagnostics based on the characterization of orbit plots with orthogonal functions | |
Zhang et al. | A multi-mode process monitoring method based on mode-correlation PCA for Marine Current Turbine | |
Lopez-Gutierrez et al. | Induction machine bearing fault detection using empirical wavelet transform | |
JP7169084B2 (en) | Condition monitoring device and condition monitoring system | |
Manzato et al. | Removing the influence of rotor harmonics for improved monitoring of offshore wind turbines | |
US11898994B2 (en) | Waveform acquisition optimization | |
Peeters et al. | Stochastic simulation assessment of an automated vibration-based condition monitoring framework for wind turbine gearbox faults | |
Ruiz de la Hermosa González-Carrato et al. | Methods and tools for the operational reliability optimisation of large-scale industrial 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 |