CN106769000B - A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method - Google Patents
A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method Download PDFInfo
- Publication number
- CN106769000B CN106769000B CN201610991785.2A CN201610991785A CN106769000B CN 106769000 B CN106769000 B CN 106769000B CN 201610991785 A CN201610991785 A CN 201610991785A CN 106769000 B CN106769000 B CN 106769000B
- Authority
- CN
- China
- Prior art keywords
- transmission path
- wind turbine
- turbine gearbox
- vibration signal
- power flow
- 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.)
- Expired - Fee Related
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 107
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004458 analytical method Methods 0.000 claims abstract description 28
- 230000011664 signaling Effects 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 238000009434 installation Methods 0.000 claims abstract description 7
- 230000008878 coupling Effects 0.000 claims description 26
- 238000010168 coupling process Methods 0.000 claims description 26
- 238000005859 coupling reaction Methods 0.000 claims description 26
- 238000012544 monitoring process Methods 0.000 claims description 17
- 230000032798 delamination Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 230000005284 excitation Effects 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000012546 transfer Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 4
- 230000007812 deficiency Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method, and the present invention relates to the transmission paths of wind turbine gearbox fault vibration signal to determine method.The present invention is to solve the problems, such as that existing method effectively can not determine main transmission path from a plurality of transmission path of wind turbine gearbox fault vibration signal.Step of the present invention are as follows: one: the characteristics of being impacted according to wind turbine gearbox failure determines the installation site of sensor；Two: determining the transmission path of wind turbine gearbox fault vibration signal；Three: obtaining the calculation formula of the fault vibration signaling path contribution amount based on power flow；Four: establishing the finite element model of wind turbine gearbox, utilize the contribution amount of each transmission path of Finite element arithmetic；Five: sorting to the contribution amount of transmission path, determine the main transmission path of wind turbine gearbox fault vibration signal.The present invention is applied to wind turbine gearbox accident analysis field.
Description
Technical field
The present invention relates to the transmission path determination sides of the wind turbine gearbox fault vibration signal based on power flow finite element method
Method.
Background technique
Critical component of the wind turbine gearbox as current wind power equipment, running environment is extremely severe, once breaking down, ties up
It repairs extremely difficult.Effective condition monitoring and fault diagnosis is carried out to wind turbine gearbox, is carried out before catastrophe failure occurs in it
Maintenance is of great significance to the reliability for improving wind power equipment.But due to the drive mechanism of wind turbine gearbox complexity, and
The quantity of fault vibration signaling path is more and is time-varying, this is the condition monitoring and fault diagnosis band of wind turbine gearbox
Carry out great challenge.On the other hand, since dissipation and disturbing effect, fault vibration signal may substantially decline in transmittance process
Subtract, the feature for the fault-signal being hidden in complex vibration signal also weakens therewith.Therefore, to fault vibration signaling path
Research help to improve the accuracy rate of wind turbine gearbox condition monitoring and fault diagnosis, there is very important engineering significance.
Traditional Transfer Path Analysis Method of Automobile need to usually test transmission function and identification load, this just needs to acquire largely
Data.Conventional method usually requires that each transmission path analyzed is non-time-varying and mutually indepedent each other, simultaneously
It is tested it is required that being easily installed sensor inside equipment under test.And contain Planetary Gear Transmission in the structure of wind turbine gearbox more
There is gear engagement vibration easily to be caused to be superimposed for mechanism, and vibration signal transmission path quantity is more and the features such as be time-varying.Furthermore its
Internal structure is also unsuitable for installing sensor, therefore traditional path analysis method is passed in the fault vibration signal of wind turbine gearbox
Applicability is insufficient in terms of passing path analysis.Transmission path research method based on energy is converted into complicated vector calculus relatively
Simple scalar plus-minus, gives an absolute measure of vibrational energy transmitting, discloses variation and the attenuation law of energy transmission.
But current research only recognizes main transmission path by comparing the power flow of a certain specific interface, lacks each paths to mesh
The vibrational energy of punctuate contributes quantifier elimination, therefore the master of fault vibration signal can not be determined by comparing the height of contribution amount
Want transmission path.
Summary of the invention
The present invention is can not be effectively from a plurality of transmitting road of wind turbine gearbox fault vibration signal in order to solve existing method
The problem of main transmission path is determined in diameter, and a kind of wind turbine gearbox failure vibration based on power flow finite element method proposed
The transmission path of dynamic signal determines method.
A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determine method by with
Lower step is realized:
Step 1: the characteristics of being impacted according to wind turbine gearbox failure determines the installation site of sensor；
Step 2: according to the installation site and fault type of sensor, wind turbine gearbox fault vibration signal is determined
Transmission path；
Step 3: it according to the transmission path for the wind turbine gearbox fault vibration signal that step 2 determines, obtains based on power
The calculation formula of the fault vibration signaling path contribution amount of stream；
Step 4: establishing the finite element model of wind turbine gearbox, utilizes the contribution of each transmission path of Finite element arithmetic
Amount；
Step 5: it sorts to the contribution amount of transmission path, determines the main transmitting road of wind turbine gearbox fault vibration signal
Diameter.
Invention effect:
The biography of the object of the present invention is to provide a kind of wind turbine gearbox fault vibration signal based on power flow finite element method
Pass path analysis method, for conventional transmission path analysis method wind turbine gearbox fault vibration signal transmission path quantity
Disadvantage that is more and being time-varying and gear internal sensor not easy to install etc. applicability deficiency, it is limited based on power flow
First method calculates fault vibration signaling path contribution amount, effectively determines the main transmission path of fault vibration signal.With tradition
Transfer Path Analysis Method of Automobile compare, the invention has the benefit that
1, the present invention is directed to the feature of wind turbine gearbox fault vibration signaling path, makes to transmit road with identical energy
The fault vibration signaling path of diameter keeps the factor by subsystem energy with equal, meets the reality of fault vibration signal
Border transmittance process and reduce the operand in calculating process.
2, the present invention is for conventional transmission path analysis method in wind turbine gearbox fault vibration signaling path quantity
Disadvantage that is more and being sensor not easy to install etc. applicability deficiency inside time-varying and equipment, to each paths pair
The analysis of the vibrational energy contribution amount of target point is set out, and is established each item of wind turbine gearbox fault vibration signal based on power flow and is passed
The calculation formula of path contributions amount is passed, improves Transfer Path Analysis Method of Automobile in the transmission path of wind turbine gearbox fault vibration signal
Analyze the applicability of aspect.
3, the present invention carries out model analysis and harmonic responding analysis to wind turbine gearbox using FInite Element, obtain target point and
The power flow of Coupling point and the time-varying of each transmission path and non-time-varying energy keep the factor, more accurately have compared with conventional method
Effect.
4, the present invention only recognizes main transmission path by comparing the power flow of a certain specific interface for existing method
Deficiency, summation is carried out to the contribution amount of each transmission path under each resonant belt frequency and obtains total contribution amount, and it is arranged
Sequence, precise and high efficiency obtain the main transmission path of wind turbine gearbox fault vibration signal.
Detailed description of the invention
Fig. 1 is the stream of the Transfer Path Analysis Method of Automobile of the wind turbine gearbox fault vibration signal based on power flow finite element method
Journey schematic diagram；
Fig. 2 is 6 transmission path figures of planetary gear fault vibration signal；
Fig. 3 is 2 transmission path figures of planet carrier bearing fault vibration signal；
Fig. 4 is that the time-varying energy of the planetary gear fault vibration signal of different gear positions keeps factor variations curve graph；
Fig. 5 is that the time-varying energy of the planet carrier bearing fault vibration signal of different gear positions keeps factor variations curve
Figure；
Fig. 6 is the change curve for the planetary gear fault vibration signaling path contribution amount that gear teeth number is 9；
Fig. 7 is the change curve for the planetary gear fault vibration signaling path contribution amount that gear teeth number is 27；
Fig. 8 is the change curve for the planet carrier bearing fault vibration signaling path contribution amount that gear teeth number is 1；
Fig. 9 is the change curve for the planet carrier bearing fault vibration signaling path contribution amount that gear teeth number is 10.
Specific embodiment
Specific embodiment 1: as shown in Figure 1, a kind of wind turbine gearbox fault vibration letter based on power flow finite element method
Number transmission path determine method the following steps are included:
Step 1: the characteristics of being impacted according to wind turbine gearbox failure determines the installation site of sensor；
Step 2: according to the installation site and fault type of sensor, wind turbine gearbox fault vibration signal is determined
Transmission path；
Step 3: it according to the transmission path for the wind turbine gearbox fault vibration signal that step 2 determines, obtains based on power
The calculation formula of the fault vibration signaling path contribution amount of stream；
Step 4: establishing the finite element model of wind turbine gearbox, utilizes the contribution of each transmission path of Finite element arithmetic
Amount；
Step 5: it sorts to the contribution amount of transmission path, determines the main transmitting road of wind turbine gearbox fault vibration signal
Diameter.
Specific embodiment 2: the present embodiment is different from the first embodiment in that: wind is determined in the step 2
The detailed process of the transmission path of electrical gearbox fault vibration signal are as follows:
Step 2 one: analysis planetary gear localized delamination failure determines the non-time-varying transmission path of fault vibration signal；
Step 2 two: analysis planetary gear localized delamination failure determines the time-varying transmission path of fault vibration signal；
Step 2 three: the outer ring localized delamination failure of analysis planet carrier bearing determines that the non-time-varying of fault vibration signal passes
Pass path；
Step 2 four: the outer ring localized delamination failure of analysis planet carrier bearing determines the time-varying transmitting of fault vibration signal
Path.
Other steps and parameter are same as the specific embodiment one.
Specific embodiment 3: the present embodiment is different from the first and the second embodiment in that: in the step 3
To the detailed process of the calculation formula of the fault vibration signaling path contribution amount based on power flow are as follows:
Step 3 one: selection each transmission path Coupling point, obtain wind turbine gearbox system capacity transmission figure and by
Subsystem energy transfer profiles；
Step 3 two: the wind turbine gearbox system capacity transmission figure that is obtained according to step 3 one and by subsystem energy
Transmission figure finds the fault vibration signaling path with identical energy transmission path；
Step 3 three: power flow and each item of the excitation at passive subsystems couple point under working condition are obtained and transmits road
The energy of diameter keeps the factor；
Step 3 four: the contribution amount percentage of i-th transmission path is obtained are as follows:
C in formulaiIt is the contribution amount percentage of i-th transmission path,It is by the transmission path of subsystem side i-th
Energy keep the factor, Pi A(ω) is the power flow at i-th transmission path Coupling point.
Other steps and parameter are the same as one or two specific embodiments.
Specific embodiment 4: unlike one of present embodiment and specific embodiment one to three: the step 3
The energy that power flow and each transmission path of the excitation at passive subsystems couple point under working condition are obtained in three is kept
The detailed process of the factor are as follows:
The factor is kept to obtain energy, practical driving source need to be removed, while successively being applied at passive subsystems couple point
The ratio of excitation, while the power flow of measuring target point, target point power flow and Coupling point power flow is that energy keeps the factor.
Step 331: m instruction point is being selected by the position in subsystem distance objective point 5cm, in each coupling
Select 1 monitoring point in position in point 5cm；In the operating condition, the power of each instruction point is obtained by sensor or emulation
Stream, and form m dimension power flow column vector
Step 3 three or two: practical driving source is removed, successively applies excitation at passive subsystems couple point, and pass through sensing
Device or emulation obtain the power flow of each monitoring point and instruction point；When applying excitation at the Coupling point in i-th of path, i-th
The monitoring point power flow in a path is pi(ω), the power flow of j-th of instruction point is p at target pointji(ω), thenFor from
The energy of i-th of monitoring point to j-th of instruction point keeps the factor, and all m × n energy keep factor composition m × n dimension energy to protect
Hold factor matrix δB(ω)；
Step 3 three or three: noteAt corresponding to Coupling point under working condition
N ties up power flow column vector, then:
δB(ω)·PA(ω)=PB(ω) (2)
I.e.
To δB(ω) takes generalized inverse δB(ω)-1It can obtain:
PA(ω)=δB(ω)-1·PB(ω) (4)
I.e.
Other steps and parameter are identical as one of specific embodiment one to three.
Specific embodiment 5: unlike one of present embodiment and specific embodiment one to four: the step 4
In establish finite element model, utilize the detailed process of the contribution amount of each transmission path of Finite element arithmetic are as follows:
Step 4 one: the characteristics of being directed to wind turbine gearbox and analysis need, and remove unwanted feature, establish simplified
Wind turbine gearbox finite element model；
Step 4 two: model analysis and harmonic responding analysis are carried out to wind turbine gearbox；
Step 4 three: being calculated the power flow at each instruction point and monitoring point using FEM Numerical Simulation, and
The non-time-varying and time-varying transmission path energy of fault vibration signal keep the factor；
Step 4 four: keeping the factor according to the energy of power flow and each transmission path at each instruction point and monitoring point,
Calculate the contribution amount of each transmission path of wind turbine gearbox fault vibration signal；
Step 4 five: it sums to each transmission path contribution amount under resonant belt frequency, obtains the total of each paths
Contribution amount.
Other steps and parameter are identical as one of specific embodiment one to four.
Embodiment one:
A kind of transmission path of wind turbine gearbox fault vibration signal based on power flow finite element method of the present embodiment determines
Method is specifically to be prepared according to the following steps:
Wind turbine gearbox generally adds two-stage to determine parallel-axes gears transmission mechanism to form by primary planet pinion, wherein planet tooth
Wheel running part is the core of wind turbine gearbox, and main component has planetary gear, sun gear, planet carrier and ring gear.According to
Impact force caused by the unique texture of wind turbine gearbox, planetary gear localized delamination failure and planet carrier shaft bearing outer-ring failure is equal
It is mostly in sagittal plane, so Sensor is arranged in the top of planet gear transmission section.
For planetary gear localized delamination failure, the source of trouble to sensor shares 6 vibration transfer paths, such as Fig. 2 institute
Show.Transmission path 1,2 and 3 is that failure is formed when being located at planetary gear and sun gear meshing point in figure, and transmission path 4,5 and 6 is
Failure is located at planetary gear and when ring gear meshing point formed.Wherein transmission path 1,3,4 and 6 is non-time-varying path, transmits road
Diameter 2 and 5 is time-varying path.
For the outer ring local fault of planet carrier bearing, the source of trouble to sensor shares 2 vibration transfer paths, such as Fig. 3
It is shown.In figure, transmission path 1 is non-time-varying path, and transmission path 2 is time-varying path.
For planetary gear localized delamination failure, the Coupling point selection of each paths is as follows: path 1 selects Coupling point to exist
On sun gear, path 2 selects Coupling point for the meshing point of planetary gear and ring gear, and path 3 selects Coupling point in planet carrier and row
Between carrier bearing, path 4 selects Coupling point for the meshing point of planetary gear and sun gear, and path 5 selects Coupling point in ring gear
On, path 6 selects Coupling point between planet carrier and planet carrier bearing.
Since the fault vibration signal in path 2 and the fault vibration signal in path 5 can all pass through each wheel of ring gear
Tooth is transmitted to sensor, so path 2 and path 5 are having the same by the subsystem energy holding factor.Similarly, path 1 with
Path 4, path 3 and path 6 are also having the same to keep the factor by subsystem energy.Internal annular gear teeth is numbered, with
The gear teeth of top face sensor are number 1, while considering symmetry of the ring gear about sensor, only to half ring gear
Number.It indicates that the ring gear gear teeth are numbered with variable n, 1 monitoring point is set near each Coupling point, is arranged near target point
3 instruction points, can obtain according to formula (3):
The contribution amount percentage that 6 transmission paths can be obtained according to formula (1), as shown in formula (7):
For planetary gear localized delamination failure, the Coupling point selection of 2 paths is as follows: the coupling of 2 paths is clicked
Select as follows: path 1 selects Coupling point between planet carrier bearing outer ring and cabinet, and path 2 selects Coupling point in planet carrier bearing
Between planet carrier.
Since 3 planetary gears are by 3 equal part of circumference, so need to only be numbered to one third ring gear.In each Coupling point
2 instruction points are arranged near target point, can be obtained according to formula (3) for 1 monitoring point of setting nearby:
In order to which the energy for being different from planetary gear fault-signal transmission path keeps factor delta in formula (8), indicated using symbol η
The energy of planet carrier bearing fault-signal transmission path keeps the factor, then the contribution amount percentage of 2 transmission paths are as follows:
Since bolt, nut and chamfering are unrelated to analyzed area of interest in wind power gear box model, so by it
Removal, establishes simplified wind turbine gearbox finite element model.Model analysis and harmonic responding analysis are carried out to wind turbine gearbox.According to
The power flow of Coupling point and target point is calculated in simulation result and time-varying and non-time-varying transmission path energy keep the factor.
Factor summation is kept to the energy of the different frequency under a certain gear teeth number, the size of value reflects time-varying transmission path
Transmission capacity size, such as Fig. 4 and Fig. 5.The factor is kept for the energy of planetary gear fault-signal transmission path, when the gear teeth are compiled
Number be 9 and 27 when obtain maximum value and minimum value respectively；For planet carrier bearing fault-signal transmission path energy keep because
Son obtains maximum value and minimum value when gear teeth number is 10 and 1 respectively.
Planet wheels fault-signal and planet carrier bearing fault-signal transmission path are calculated separately according to formula (7) and formula (9)
Contribution amount, respectively as shown in Fig. 6-Fig. 9.To three resonant belt frequencies (550Hz~600Hz, 1450Hz~1650Hz and
2200Hz~2850Hz) under the contribution amount of each transmission path sum, obtain total contribution amount.Planetary gear failure is believed
Number, as n=9, total contribution amount in path 1 to path 6 is respectively 2.97,2.33,14.08,0.26,0.23 and 1.14；Work as n=
When 27, total contribution amount in path 1 to path 6 is respectively 2.22,3.38,13.39,0.10,0.16 and 1.75.For planet carrier shaft
Fault-signal is held, as n=1, total contribution amount in path 1 is 12.97, and total contribution amount in path 2 is 8.03；As n=10, road
Total contribution amount of diameter 1 is 14.25, and total contribution amount in path 2 is 6.75.The contribution amount sequence of each transmission path is as shown in table 1,
For planetary gear fault-signal, path 3 is main transmission path；For planet carrier bearing fault-signal, path 1 is mainly to transmit
Path.From the experimental results, the present invention can have the spy of time-varying for wind turbine gearbox fault vibration signaling path
The contribution amount of each transmission path is calculated based on power flow finite element method for point, passes through the sequence accurate and effective to contribution amount
Determine the main transmission path of wind turbine gearbox fault vibration signal.
Table 1 is total contribution amount sequencing table of each transmission path
Claims (4)
1. a kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method, feature
Be, the transmission path of the wind turbine gearbox fault vibration signal determine method the following steps are included:
Step 1: the characteristics of being impacted according to wind turbine gearbox failure determines the installation site of sensor；
Step 2: according to the installation site and fault type of sensor, the transmitting of wind turbine gearbox fault vibration signal is determined
Path；
Step 3: it according to the transmission path for the wind turbine gearbox fault vibration signal that step 2 determines, obtains based on power flow
The calculation formula of fault vibration signaling path contribution amount；
Step 4: establishing the finite element model of wind turbine gearbox, utilizes the contribution amount of each transmission path of Finite element arithmetic；
Step 5: it sorts to the contribution amount of transmission path, determines the main transmission path of wind turbine gearbox fault vibration signal；
The detailed process of the transmission path of wind turbine gearbox fault vibration signal is determined in the step 2 are as follows:
Step 2 one: analysis planetary gear localized delamination failure determines the non-time-varying transmission path of fault vibration signal；
Step 2 two: analysis planetary gear localized delamination failure determines the time-varying transmission path of fault vibration signal；
Step 2 three: the outer ring localized delamination failure of analysis planet carrier bearing determines the non-time-varying transmitting road of fault vibration signal
Diameter；
Step 2 four: the outer ring localized delamination failure of analysis planet carrier bearing determines the time-varying transmission path of fault vibration signal.
2. a kind of transmitting of wind turbine gearbox fault vibration signal based on power flow finite element method according to claim 1
Determining method of path, which is characterized in that the fault vibration signaling path contribution based on power flow is obtained in the step 3
The detailed process of the calculation formula of amount are as follows:
Step 3 one: the Coupling point of each transmission path of selection obtains wind turbine gearbox system capacity transmission figure and by mover
System capacity transmission figure；
Step 3 two: the wind turbine gearbox system capacity transmission figure that is obtained according to step 3 one and by subsystem energy transmission
Figure finds the fault vibration signaling path with identical energy transmission path；
Step 3 three: power flow at the passive subsystems couple point of excitation under working condition and each transmission path are obtained
Energy keeps the factor；
Step 3 four: the contribution amount percentage of i-th transmission path is obtained are as follows:
C in formulaiIt is the contribution amount percentage of i-th transmission path,It is by the energy of i-th transmission path in subsystem side
Amount keeps the factor, PiA(ω) is the power flow at i-th transmission path Coupling point.
3. a kind of transmitting of wind turbine gearbox fault vibration signal based on power flow finite element method according to claim 2
Determining method of path, which is characterized in that obtain excitation under working condition in the step 3 three and be derived from passive subsystems couple point
The energy of the power flow at place and each transmission path keeps the detailed process of the factor are as follows:
Step 331: m instruction point is being selected by the position in subsystem distance objective point 5cm, in each Coupling point 5cm
Select 1 monitoring point in interior position；In the operating condition, the power flow of each instruction point is obtained by sensor or emulation, and
It forms m and ties up power flow column vector
Step 3 three or two: removing practical driving source, successively applies excitation at passive subsystems couple point, and pass through sensor or
Emulation obtains the power flow of each monitoring point and instruction point；When applying excitation at the Coupling point in i-th of path, i-th of tunnel
The monitoring point power flow of diameter is pi (ω), and the power flow of j-th of instruction point is p at target pointji(ω), thenFor from i-th
The energy of a monitoring point to j-th of instruction point keeps the factor, and all m × n energy keep factor composition m × n dimension energy to keep
Factor matrix δB(ω)；
Step 3 three or three: noteTo be tieed up corresponding to n at Coupling point under working condition
Power flow column vector, then:
δB(ω)·PA(ω)=PB(ω) (2)
I.e.
To δB(ω) takes generalized inverse δB(ω)-1It can obtain:
PA(ω)=δB(ω)-1·PB(ω) (4)
I.e.
4. a kind of transmitting of wind turbine gearbox fault vibration signal based on power flow finite element method according to claim 3
Determining method of path, which is characterized in that establish finite element model in the step 4, transmit road using each item of Finite element arithmetic
The detailed process of the contribution amount of diameter are as follows:
Step 4 one: simplified wind turbine gearbox finite element model is established；
Step 4 two: model analysis and harmonic responding analysis are carried out to wind turbine gearbox；
Step 4 three: the power flow and failure at each instruction point and monitoring point is calculated using FEM Numerical Simulation
The non-time-varying and time-varying transmission path energy of vibration signal keep the factor；
Step 4 four: the factor is kept according to the energy of power flow and each transmission path at each instruction point and monitoring point, is calculated
The contribution amount of wind turbine gearbox each transmission path of fault vibration signal；
Step 4 five: it sums to each transmission path contribution amount under resonant belt frequency, obtains total contribution of each paths
Amount.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610991785.2A CN106769000B (en) | 2016-11-10 | 2016-11-10 | A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610991785.2A CN106769000B (en) | 2016-11-10 | 2016-11-10 | A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106769000A CN106769000A (en) | 2017-05-31 |
| CN106769000B true CN106769000B (en) | 2019-03-26 |
Family
ID=58973018
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610991785.2A Expired - Fee Related CN106769000B (en) | 2016-11-10 | 2016-11-10 | A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106769000B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107657094B (en) * | 2017-09-14 | 2020-12-29 | 华北电力大学（保定） | Visual analysis method for energy distribution under collision impact of disc structure |
| CN108073757B (en) * | 2017-09-14 | 2021-07-13 | 华北电力大学（保定） | Beam structure natural frequency analysis method based on power flow |
| ES2754278B2 (en) * | 2018-09-28 | 2022-01-25 | Univ Navarra Publica | DIAGNOSIS OF FAULTS IN PLANETARY GEARS |
| CN113465850B (en) * | 2021-02-07 | 2023-09-08 | 西北工业大学 | Method, test device and test method for identifying mechanical vibration signal path |
| CN112985811B (en) * | 2021-05-12 | 2021-09-07 | 成都飞机工业（集团）有限责任公司 | Structure fault positioning method based on virtual excitation source |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101271022A (en) * | 2008-05-15 | 2008-09-24 | 上海交通大学 | Transmission path detecting system for vehicle system structure vibration and noise |
| CN102419252A (en) * | 2011-08-18 | 2012-04-18 | 黑龙江大学 | Optical fiber online type detection device of gear box of high-speed train |
| CN103076177A (en) * | 2013-01-16 | 2013-05-01 | 昆明理工大学 | Rolling bearing fault detection method based on vibration detection |
| CN104502127A (en) * | 2014-11-19 | 2015-04-08 | 哈尔滨工程大学 | Outfield acoustically-driven ship vibration noise transmission path analysis method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007147405A (en) * | 2005-11-25 | 2007-06-14 | Matsushita Electric Works Ltd | Anomaly monitoring apparatus |
| DE102015207134A1 (en) * | 2015-04-20 | 2016-10-20 | Prüftechnik Dieter Busch AG | Method for detecting vibrations of a device and vibration detection system |
-
2016
- 2016-11-10 CN CN201610991785.2A patent/CN106769000B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101271022A (en) * | 2008-05-15 | 2008-09-24 | 上海交通大学 | Transmission path detecting system for vehicle system structure vibration and noise |
| CN102419252A (en) * | 2011-08-18 | 2012-04-18 | 黑龙江大学 | Optical fiber online type detection device of gear box of high-speed train |
| CN103076177A (en) * | 2013-01-16 | 2013-05-01 | 昆明理工大学 | Rolling bearing fault detection method based on vibration detection |
| CN104502127A (en) * | 2014-11-19 | 2015-04-08 | 哈尔滨工程大学 | Outfield acoustically-driven ship vibration noise transmission path analysis method |
Non-Patent Citations (1)
| Title |
|---|
| 基于传递路径的多级齿轮箱齿轮裂纹故障识别;徐玉秀 等;《仪器仪表学报》;20160531;第37卷(第5期);第1018-1024页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106769000A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106769000B (en) | A kind of transmission path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method | |
| Zhao et al. | Deep convolutional neural network based planet bearing fault classification | |
| He et al. | Multifractal entropy based adaptive multiwavelet construction and its application for mechanical compound-fault diagnosis | |
| Praveenkumar et al. | Pattern recognition based on-line vibration monitoring system for fault diagnosis of automobile gearbox | |
| Chen et al. | Generator bearing fault diagnosis for wind turbine via empirical wavelet transform using measured vibration signals | |
| Bafroui et al. | Application of wavelet energy and Shannon entropy for feature extraction in gearbox fault detection under varying speed conditions | |
| CN101354312B (en) | Bearing failure diagnosis system | |
| CN102913432B (en) | Reciprocating compressor fault diagnosis method based on dimensionless indexes | |
| Wang et al. | Meshing frequency modulation (MFM) index-based kurtogram for planet bearing fault detection | |
| CN105092239B (en) | A kind of initial failure of gear detection method | |
| Wang et al. | A time–frequency-based maximum correlated kurtosis deconvolution approach for detecting bearing faults under variable speed conditions | |
| CN103104509A (en) | Obtaining method of variable frequency water pump full working condition operating state | |
| Feng et al. | Spectral negentropy based sidebands and demodulation analysis for planet bearing fault diagnosis | |
| Li et al. | Multipoint optimal minimum entropy deconvolution adjusted for automatic fault diagnosis of hoist bearing | |
| CN104408322B (en) | Rotating mechanical device fault diagnosis method capable of synthesizing multisource fault probability likelihood credibility | |
| CN106441843B (en) | A kind of rotating machinery fault method for waveform identification | |
| Luo et al. | Vibration mechanism and improved phenomenological model of the planetary gearbox with broken ring gear fault | |
| CN102562565A (en) | Method and device for water pump delivery inspection | |
| Li et al. | Optimal SES selection based on SVD and its application to incipient bearing fault diagnosis | |
| CN102914594A (en) | Method for quantitatively diagnosing transverse cracks of pipeline on basis of frequency contour lines | |
| CN108225762B (en) | A kind of single gear tooth crackle broken teeth fault identification diagnostic method | |
| Xue et al. | Real-time diagnosis of an in-wheel motor of an electric vehicle based on dynamic Bayesian networks | |
| CN109682562B (en) | Multi-axis durability to single-axis durability test method | |
| Yang et al. | Vibration AM-FM sidebands mechanism of planetary gearbox with tooth root cracked planet gear | |
| Zhukovytskyy et al. | Information-measuring test system of diesel locomotive hydraulic transmissions |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190326 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |