CN106769000A - A kind of bang path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method - Google Patents

A kind of bang path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method Download PDF

Info

Publication number
CN106769000A
CN106769000A CN201610991785.2A CN201610991785A CN106769000A CN 106769000 A CN106769000 A CN 106769000A CN 201610991785 A CN201610991785 A CN 201610991785A CN 106769000 A CN106769000 A CN 106769000A
Authority
CN
China
Prior art keywords
omega
bang path
wind turbine
path
turbine gearbox
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
Application number
CN201610991785.2A
Other languages
Chinese (zh)
Other versions
CN106769000B (en
Inventor
黄文涛
孙宏健
窦宏印
王伟杰
赵学增
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201610991785.2A priority Critical patent/CN106769000B/en
Publication of CN106769000A publication Critical patent/CN106769000A/en
Application granted granted Critical
Publication of CN106769000B publication Critical patent/CN106769000B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/028Acoustic or vibration analysis

Abstract

A kind of bang path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method, and the bang path the present invention relates to wind turbine gearbox fault vibration signal determines method.The present invention is to solve the problems, such as that existing method effectively cannot determine main bang path from a plurality of bang path of wind turbine gearbox fault vibration signal.Step of the present invention is:One:The characteristics of being impacted according to wind turbine gearbox failure, determines the installation site of sensor;Two:Determine the bang path of wind turbine gearbox fault vibration signal;Three:Obtain the computing formula of the fault vibration signaling path contribution amount based on power flow;Four:The FEM model of wind turbine gearbox is set up, using the contribution amount of each bar bang path of Finite element arithmetic;Five:Contribution amount sequence to bang path, determines the main bang path of wind turbine gearbox fault vibration signal.The present invention is applied to wind turbine gearbox accident analysis field.

Description

A kind of transmission of the wind turbine gearbox fault vibration signal based on power flow finite element method Determining method of path
Technical field
Bang path determination side the present invention relates to be based on the wind turbine gearbox fault vibration signal of power flow finite element method Method.
Background technology
Wind turbine gearbox as current wind power equipment critical component, running environment is extremely severe, once break down, tie up Repair extremely difficult.Effective condition monitoring and fault diagnosis are carried out to wind turbine gearbox, was carried out before it occurs catastrophe failure Safeguard, the reliability to improving wind power equipment is significant.But because the drive mechanism of wind turbine gearbox is complicated, and The quantity of fault vibration signaling path is more and be time-varying, and this is the condition monitoring and fault diagnosis band of wind turbine gearbox Carry out great challenge.On the other hand, due to dissipation and disturbing effect, fault vibration signal may significantly decline in transmittance process Subtract, the feature of the fault-signal being hidden in complex vibration signal also weakens therewith.Therefore, to fault vibration signaling path Research be favorably improved the accuracy rate of wind turbine gearbox condition monitoring and fault diagnosis, with very important engineering significance.
Traditional Transfer Path Analysis Method of Automobile need to generally test transmission function and identification load, and it is substantial amounts of that this is accomplished by collection Data.Conventional method usually requires that each bar bang path analyzed is non-time-varying and separate each other, while It is required that sensor is easily installed inside equipment under test being tested.And contain Planetary Gear Transmission in the structure of wind turbine gearbox more Mechanism, easily causes vibration to be superimposed with gear engagement, and vibration signal bang path quantity is more and the features such as be time-varying.In addition its Internal structure is also unsuitable for install sensor, therefore traditional path analysis method is passed in the fault vibration signal of wind turbine gearbox Pass path analysis aspect applicability not enough.Bang path research method based on energy is converted into complicated vector calculus relatively Simple scalar plus-minus, gives an absolute measure of vibrational energy transmission, discloses change and the attenuation law of energy transmission. But current research only recognizes main bang path by the power flow of relatively more a certain specific interface, lacks each paths to mesh The vibrational energy contribution quantifier elimination of punctuate, therefore the master of fault vibration signal cannot be determined by contrasting the height of contribution amount Want bang path.
The content of the invention
The present invention is cannot effectively from a plurality of transmission road of wind turbine gearbox fault vibration signal in order to solve existing method The problem of main bang path is determined in footpath, and a kind of wind turbine gearbox failure based on power flow finite element method for proposing is shaken The bang path of dynamic signal determines method.
A kind of bang 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 one: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 Bang path;
Step 3:The bang path of the wind turbine gearbox fault vibration signal determined according to step 2, obtains based on power The computing formula of the fault vibration signaling path contribution amount of stream;
Step 4:The FEM model of wind turbine gearbox is set up, using the contribution of each bar bang path of Finite element arithmetic Amount;
Step 5:Contribution amount sequence to bang path, determines the main transmission road of wind turbine gearbox fault vibration signal Footpath.
Invention effect:
It is an object of the invention to provide a kind of biography of the 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 bang path quantity It is many and be time-varying, and gear internal is difficult the shortcoming of the aspect such as install sensor applicability deficiency, it is limited based on power flow First method calculates fault vibration signaling path contribution amount, effectively determines the main bang path of fault vibration signal.With tradition Transfer Path Analysis Method of Automobile compare, beneficial effects of the present invention are:
1st, the present invention makes to transmit road with identical energy for the feature of wind turbine gearbox fault vibration signaling path The fault vibration signaling path in footpath have it is equal the factor is kept by subsystem energy, meet the reality of fault vibration signal Border transmittance process and reduce the operand in calculating process.
2nd, the present invention is directed to conventional transmission path analysis method in wind turbine gearbox fault vibration signaling path quantity It is many and be time-varying, and device interior is difficult the shortcoming of the aspect such as install sensor applicability deficiency, to each paths pair The analysis of the vibrational energy contribution amount of impact point is set out, and is set up each bar of wind turbine gearbox fault vibration signal based on power flow and is passed The computing formula of path contributions amount is passed, bang path of the Transfer Path Analysis Method of Automobile in wind turbine gearbox fault vibration signal is improved Applicability in terms of analysis.
3rd, the present invention carries out model analysis and harmonic responding analysis using FInite Element to wind turbine gearbox, obtain impact point and The time-varying of the power flow of Coupling point and each bang path keeps the factor with non-time-varying energy, more accurately has compared with conventional method Effect.
4th, the present invention only recognizes main bang path for existing method by the power flow of relatively more a certain specific interface Deficiency, to each resonant belt frequency under the contribution amount of each bar bang path carry out summation and obtain total contribution amount, and it is arranged Sequence, the main bang path for obtaining wind turbine gearbox fault vibration signal of precise and high efficiency.
Brief description of the drawings
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 bang path figures of planetary gear fault vibration signal;
Fig. 3 is 2 bang 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 map;
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 of the planetary gear fault vibration signaling path contribution amount that gear teeth numbering is 9;
Fig. 7 is the change curve of the planetary gear fault vibration signaling path contribution amount that gear teeth numbering is 27;
Fig. 8 is the change curve of the planet carrier bearing fault vibration signaling path contribution amount that gear teeth numbering is 1;
Fig. 9 is the change curve of the planet carrier bearing fault vibration signaling path contribution amount that gear teeth numbering is 10.
Specific embodiment
Specific embodiment one:As shown in figure 1, a kind of wind turbine gearbox fault vibration letter based on power flow finite element method Number bang path determine that method is comprised the following steps:
Step one: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 Bang path;
Step 3:The bang path of the wind turbine gearbox fault vibration signal determined according to step 2, obtains based on power The computing formula of the fault vibration signaling path contribution amount of stream;
Step 4:The FEM model of wind turbine gearbox is set up, using the contribution of each bar bang path of Finite element arithmetic Amount;
Step 5:Contribution amount sequence to bang path, determines the main transmission road of wind turbine gearbox fault vibration signal Footpath.
Specific embodiment two:Present embodiment from unlike specific embodiment one:Determine wind in the step 2 The detailed process of the bang path of electrical gearbox fault vibration signal is:
Step 2 one:Planetary localized delamination failure is analyzed, the non-time-varying bang path of fault vibration signal is determined;
Step 2 two:Planetary localized delamination failure is analyzed, the time-varying bang path of fault vibration signal is determined;
Step 2 three:The outer ring localized delamination failure of planet carrier bearing is analyzed, determines that the non-time-varying of fault vibration signal is passed Pass path;
Step 2 four:The outer ring localized delamination failure of planet carrier bearing is analyzed, the time-varying transmission of fault vibration signal is determined Path.
Other steps and parameter are identical with specific embodiment one.
Specific embodiment three:Present embodiment from unlike specific embodiment one or two:In the step 3 The detailed process of computing formula to the fault vibration signaling path contribution amount based on power flow is:
Step 3 one:The Coupling point of each bar bang path is selected, wind turbine gearbox system capacity transmission figure and quilt is obtained 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 bang path;
Step 3 three:Excitation comes from power flow and each bar the transmission road at passive subsystems couple point under acquisition working condition The energy in footpath keeps the factor;
Step 3 four:The contribution amount percentage for obtaining i-th bang path is:
C in formulaiIt is the i-th contribution amount percentage of bang path,It is by the bang path of subsystem side i-th Energy keep the factor, Pi A(ω) is the power flow at i-th bang path Coupling point.
Other steps and parameter are identical with specific embodiment one or two.
Specific embodiment four:Unlike one of present embodiment and specific embodiment one to three:The step 3 Excitation comes from the energy holding of the power flow and each bar bang path at passive subsystems couple point under acquisition working condition in three The detailed process of the factor is:
The factor is kept to obtain energy, actual driving source need to be removed, while applying at passive subsystems couple point successively Excitation, while the ratio of the power flow of measuring target point, impact point power flow and Coupling point power flow is energy keeps the factor.
Step 331:The selected m instruction point in position in by subsystem distance objective point 5cm, in each coupling Select 1 monitoring point in position in point 5cm;Under in working order, each power for indicating point is obtained by sensor or emulation Stream, and constitute m dimension power flow column vectors
Step 3 three or two:Actual driving source is removed, applies excitation successively at passive subsystems couple point, and by sensing Device or emulation obtain each monitoring point and indicate the power flow of point;When excitation is applied at i-th Coupling point in path, i-th The monitoring point power flow in individual path is pi(ω), the power flow that j-th indicates point at impact point is pji(ω), thenBe from I-th monitoring point indicates the energy of point to keep the factor to j-th, and all m × n energy keep factor composition m × n dimension energy to protect Hold factor matrix δB(ω);
Step 3 three or three:NoteIt is corresponding at 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(ω)-1Can obtain:
PA(ω)=δB(ω)-1·PB(ω) (4)
I.e.
Other steps and parameter are identical with one of specific embodiment one to three.
Specific embodiment five:Unlike one of present embodiment and specific embodiment one to four:The step 4 In set up FEM model, the detailed process using the contribution amount of each bar bang path of Finite element arithmetic is:
Step 4 one:The characteristics of for wind turbine gearbox and analysis need, and remove unwanted feature, set up simplification Wind turbine gearbox FEM model;
Step 4 two:Model analysis and harmonic responding analysis are carried out to wind turbine gearbox;
Step 4 three:The power flow at each instruction point and monitoring point is calculated using FEM Numerical Simulation, and The non-time-varying and time-varying bang path energy of fault vibration signal keep the factor;
Step 4 four:Energy according to each power flow indicated at point and monitoring point and each bang path keeps the factor, Calculate the contribution amount of each bar bang path of wind turbine gearbox fault vibration signal;
Step 4 five:Each bar bang path contribution amount under to resonant belt frequency is sued for peace, and obtains the total of each paths Contribution amount.
Other steps and parameter are identical with one of specific embodiment one to four.
Embodiment one:
A kind of bang 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 following steps:
Wind turbine gearbox typically adds two-stage to determine parallel-axes gears transmission mechanism to constitute by primary planet pinion, wherein planet tooth Wheel running part is the core of wind turbine gearbox, and critical piece has planetary gear, sun gear, planet carrier and ring gear.According to Impulsive force produced 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 setting Sensor in the top of planet geared portions.
For planetary localized delamination failure, the source of trouble to sensor has 6 vibration transfer paths, such as Fig. 2 institutes Show.Bang path 1,2 and 3 is that failure is formed when being located at planetary gear and sun gear meshing point in figure, and bang path 4,5 and 6 is What failure was formed when being located at planetary gear with ring gear meshing point.Wherein bang path 1,3,4 and 6 is non-time-varying path, transmits road Footpath 2 and 5 is time-varying path.
For the outer ring local fault of planet carrier bearing, the source of trouble to sensor has 2 vibration transfer paths, such as Fig. 3 It is shown.In figure, bang path 1 is non-time-varying path, and bang path 2 is time-varying path.
For planetary localized delamination failure, the coupling point selection of each paths is as follows:The selection Coupling point of path 1 exists On sun gear, the selection Coupling point of path 2 is the meshing point of planetary gear and ring gear, the selection of path 3 Coupling point planet carrier and row Between carrier bearing, the selection Coupling point of path 4 is the meshing point of planetary gear and sun gear, and path 5 selects Coupling point in ring gear On, between the selection Coupling point planet carrier of path 6 and planet carrier bearing.
Because the fault vibration signal in path 2 and the fault vibration signal in path 5 can all be taken turns by each of ring gear Tooth is delivered to sensor, so path 2 with path 5 there is identical to keep the factor by subsystem energy.Similarly, path 1 with Also there is identical to keep the factor by subsystem energy with path 6 for path 4, path 3.Internal annular gear teeth is numbered, with The upper upright gear teeth to sensor are numbering 1, while in view of symmetry of the ring gear on sensor, only to half ring gear Numbering.Represent that the ring gear gear teeth are numbered with variable n, 1 monitoring point is set near each Coupling point, set near impact point 3 instruction points, can obtain according to formula (3):
6 contribution amount percentages of bang path can be obtained according to formula (1), as shown in formula (7):
For planetary 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 is selected between Coupling point planet carrier bearing outer ring and casing, the selection of path 2 Coupling point planet carrier bearing Between planet carrier.
Because 3 planetary gears are by the decile of circumference 3, so only need to be to 1/3rd ring gear numberings.In each Coupling point 1 monitoring point is nearby set, 2 instruction points are set near impact point, can be obtained according to formula (3):
Factor delta is kept in order to be different from the energy of planetary gear fault-signal bang path in formula (8), is represented using symbol η The energy of planet carrier bearing fault-signal bang path keeps the factor, then the contribution amount percentage of 2 bang paths is:
Due to bolt in wind power gear box model, nut and chamfering are unrelated to analyzed area of interest, thus by its Removal, sets up simplified wind turbine gearbox FEM model.Model analysis and harmonic responding analysis are carried out to wind turbine gearbox.According to Simulation result is calculated the power flow of Coupling point and impact point, and time-varying keeps the factor with non-time-varying bang path energy.
The energy of the different frequency under being numbered to a certain gear teeth keeps factor summation, the size reflection time-varying bang path of value Transmission capacity size, such as Fig. 4 and Fig. 5.Energy for planetary gear fault-signal bang path keeps the factor, when the gear teeth are compiled Number for 9 and 27 when obtain maximum and minimum value respectively;For planet carrier bearing fault-signal bang path energy keep because Son, maximum and minimum value are obtained when gear teeth numbering is 10 and 1 respectively.
Calculate planet wheels fault-signal and planet carrier bearing fault-signal bang path respectively 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 bar bang path sued for peace, obtain total contribution amount.Believe for planetary gear failure 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 in footpath 1 is 14.25, and total contribution amount in path 2 is 6.75.The contribution amount sequence of each bar bang path is as shown in table 1, For planetary gear fault-signal, path 3 is main bang path;For planet carrier bearing fault-signal, path 1 is main transmission Path.Knowable to experimental result, the present invention can be directed to wind turbine gearbox fault vibration signaling path has the spy of time-varying Point, the contribution amount of each bar bang path is calculated based on power flow finite element method, by the sequence accurate and effective to contribution amount Determine the main bang path of wind turbine gearbox fault vibration signal.
Table 1 is total contribution amount sequencing table of each bang path

Claims (5)

1. a kind of bang path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method, its feature It is that the bang path of the wind turbine gearbox fault vibration signal determines that method is comprised the following steps:
Step one: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 transmission of wind turbine gearbox fault vibration signal is determined Path;
Step 3:The bang path of the wind turbine gearbox fault vibration signal determined according to step 2, obtains based on power flow The computing formula of fault vibration signaling path contribution amount;
Step 4:The FEM model of wind turbine gearbox is set up, using the contribution amount of each bar bang path of Finite element arithmetic;
Step 5:Contribution amount sequence to bang path, determines the main bang path of wind turbine gearbox fault vibration signal.
2. the transmission of a kind of wind turbine gearbox fault vibration signal based on power flow finite element method according to claim 1 Determining method of path, it is characterised in that the tool of the bang path of wind turbine gearbox fault vibration signal is determined in the step 2 Body process is:
Step 2 one:Planetary localized delamination failure is analyzed, the non-time-varying bang path of fault vibration signal is determined;
Step 2 two:Planetary localized delamination failure is analyzed, the time-varying bang path of fault vibration signal is determined;
Step 2 three:The outer ring localized delamination failure of planet carrier bearing is analyzed, the non-time-varying transmission road of fault vibration signal is determined Footpath;
Step 2 four:The outer ring localized delamination failure of planet carrier bearing is analyzed, the time-varying bang path of fault vibration signal is determined.
3. a kind of wind turbine gearbox fault vibration signal based on power flow finite element method according to claim 1 and 2 Bang path determines method, it is characterised in that the fault vibration signaling path based on power flow is obtained in the step 3 The detailed process of the computing formula of contribution amount is:
Step 3 one:The Coupling point of each bar bang path is selected, wind turbine gearbox system capacity transmission figure is obtained 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 bang path;
Step 3 three:Obtain excitation under working condition and come from power flow at passive subsystems couple point and each bar bang path Energy keeps the factor;
Step 3 four:The contribution amount percentage for obtaining i-th bang path is:
C i = δ i B ( ω ) · P i A ( ω ) Σ j = 1 n δ j B ( ω ) · P j A ( ω ) × 100 % - - - ( 1 )
C in formulaiIt is the i-th contribution amount percentage of bang path,It is by the energy of bang path of subsystem side i-th Keep the factor, Pi A(ω) is the power flow at i-th bang path Coupling point.
4. the transmission of a kind of wind turbine gearbox fault vibration signal based on power flow finite element method according to claim 3 Determining method of path, it is characterised in that excitation comes from passive subsystems couple point under working condition is obtained in the step 3 three The energy of the power flow at place and each bar bang path keeps the detailed process of the factor to be:
Step 331:The selected m instruction point in position in by subsystem distance objective point 5cm, in each Coupling point 5cm Select 1 monitoring point in interior position;Under in working order, each power flow for indicating point is obtained by sensor or emulation, and Composition m dimension power flow column vectors
Step 3 three or two:Remove actual driving source, apply to encourage successively at passive subsystems couple point, and by sensor or Emulation obtains each monitoring point and indicates the power flow of point;When excitation is applied at i-th Coupling point in path, i-th tunnel The monitoring point power flow in footpath is pi(ω), the power flow that j-th indicates point at impact point is pji(ω), thenIt is from i-th Individual monitoring point indicates the energy of point to keep the factor to j-th, and all m × n energy keep factor composition m × n dimension energy to keep Factor matrix δB(ω);
Step 3 three or three:NoteIt is to be tieed up corresponding to n at Coupling point under working condition Power flow column vector, then:
δB(ω)·PA(ω)=PB(ω) (2)
I.e.
δ 11 B ( ω ) ... δ 1 i B ( ω ) ... δ 1 n B ( ω ) . . . . . . . . . . . . . . . δ j 1 B ( ω ) ... δ j i B ( ω ) ... δ j n B ( ω ) . . . . . . . . . . . . . . . δ m 1 B ( ω ) ... δ m i B ( ω ) ... δ m n B ( ω ) P 1 A ( ω ) . . . P i A ( ω ) . . . P n A ( ω ) = P 1 B ( ω ) . . . P j B ( ω ) . . . P m B ( ω ) - - - ( 3 )
To δB(ω) takes generalized inverse δB(ω)-1Can obtain:
PA(ω)=δB(ω)-1·PB(ω) (4)
I.e.
P 1 A ( ω ) . . . P i A ( ω ) . . . P n A ( ω ) = δ 11 B ( ω ) ... δ 1 i B ( ω ) ... δ 1 n B ( ω ) . . . . . . . . . . . . . . . δ j 1 B ( ω ) ... δ j i B ( ω ) ... δ j n B ( ω ) . . . . . . . . . . . . . . . δ m 1 B ( ω ) ... δ m i B ( ω ) ... δ m n B ( ω ) + P 1 B ( ω ) . . . P j B ( ω ) . . . P m B ( ω ) - - - ( 5 ) .
5. a kind of wind turbine gearbox fault vibration signal based on power flow finite element method according to claim 1,2 or 4 Bang path determine method, it is characterised in that FEM model is set up in the step 4, using each bar of Finite element arithmetic The detailed process of the contribution amount of bang path is:
Step 4 one:Set up simplified wind turbine gearbox FEM model;
Step 4 two:Model analysis and harmonic responding analysis are carried out to wind turbine gearbox;
Step 4 three:The power flow at each instruction point and monitoring point, and failure are calculated using FEM Numerical Simulation The non-time-varying and time-varying bang path energy of vibration signal keep the factor;
Step 4 four:Energy according to each power flow indicated at point and monitoring point and each bang path keeps the factor, calculates The contribution amount of each bar bang path of wind turbine gearbox fault vibration signal;
Step 4 five:Each bar bang path contribution amount under to resonant belt frequency is sued for peace, and obtains total contribution of each paths Amount.
CN201610991785.2A 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 Active CN106769000B (en)

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 true CN106769000A (en) 2017-05-31
CN106769000B CN106769000B (en) 2019-03-26

Family

ID=58973018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610991785.2A Active 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)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107657094A (en) * 2017-09-14 2018-02-02 华北电力大学(保定) The visual analysis method of Energy distribution under disc structure collision impact
CN108073757A (en) * 2017-09-14 2018-05-25 华北电力大学(保定) A kind of girder construction natural frequencies analysis method based on power flow
CN112985811A (en) * 2021-05-12 2021-06-18 成都飞机工业(集团)有限责任公司 Structure fault positioning method based on virtual excitation source
ES2754278R1 (en) * 2018-09-28 2021-07-01 Univ Navarra Publica FAULT DIAGNOSIS IN PLANETARY GEARS

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007147405A (en) * 2005-11-25 2007-06-14 Matsushita Electric Works Ltd Anomaly monitoring apparatus
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
WO2016169554A1 (en) * 2015-04-20 2016-10-27 Prüftechnik Dieter Busch AG Method for detecting vibrations of a device and vibration detection system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007147405A (en) * 2005-11-25 2007-06-14 Matsushita Electric Works Ltd Anomaly monitoring apparatus
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
WO2016169554A1 (en) * 2015-04-20 2016-10-27 Prüftechnik Dieter Busch AG Method for detecting vibrations of a device and vibration detection system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐玉秀 等: "基于传递路径的多级齿轮箱齿轮裂纹故障识别", 《仪器仪表学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107657094A (en) * 2017-09-14 2018-02-02 华北电力大学(保定) The visual analysis method of Energy distribution under disc structure collision impact
CN108073757A (en) * 2017-09-14 2018-05-25 华北电力大学(保定) A kind of girder construction natural frequencies analysis method based on power flow
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
ES2754278R1 (en) * 2018-09-28 2021-07-01 Univ Navarra Publica FAULT DIAGNOSIS IN PLANETARY GEARS
CN112985811A (en) * 2021-05-12 2021-06-18 成都飞机工业(集团)有限责任公司 Structure fault positioning method based on virtual excitation source
CN112985811B (en) * 2021-05-12 2021-09-07 成都飞机工业(集团)有限责任公司 Structure fault positioning method based on virtual excitation source

Also Published As

Publication number Publication date
CN106769000B (en) 2019-03-26

Similar Documents

Publication Publication Date Title
CN106769000A (en) A kind of bang path of the wind turbine gearbox fault vibration signal based on power flow finite element method determines method
CN103234748B (en) Klingelnberg bevel gear fault diagnosis method based on sensitive IMF (instinct mode function) components
Wang et al. A new SKRgram based demodulation technique for planet bearing fault detection
CN101354312B (en) Bearing failure diagnosis system
CN103499443B (en) A kind of gear distress is without key phase angular domain average computation order analysis method
CN105092239B (en) A kind of initial failure of gear detection method
Wang et al. Meshing frequency modulation (MFM) index-based kurtogram for planet bearing fault detection
CN103162959B (en) Multifunctional gear sample plate
CN102736562B (en) Knowledge base construction method oriented to fault diagnosis and fault prediction of numerical control machine tool
CN107228772A (en) A kind of shield tunnel construction method for estimating damage
CN102562565B (en) Method and device for water pump delivery inspection
CN106052743A (en) Method for evaluating effects on frequency response functions by transducer quality
James III et al. Modal parameter extraction from large operating structures using ambient excitation
CN109596349A (en) A kind of decelerator trouble diagnostic method based on VMD and PCT
Gui et al. Detection and localization of tooth breakage fault on wind turbine planetary gear system considering gear manufacturing errors
Yu et al. Analytical vibration signal model and signature analysis in resonance region for planetary gearbox fault diagnosis
Yoon et al. Planetary gearbox fault diagnosis using a single piezoelectric strain sensor
CN105115916A (en) Method for monitoring equipment wear during CRH train running
CN109682562B (en) Multi-axis durability to single-axis durability test method
El Morsy et al. Vehicle gearbox fault diagnosis based on cepstrum analysis
CN109470481A (en) A kind of engine rig test method
Saidi et al. Application of higher order spectral analysis for rotor broken bar detection in induction machines
Zhukovytskyy et al. Information-measuring test system of diesel locomotive hydraulic transmissions
CN203238609U (en) Drilling platform lifting transmission device capable of monitoring climb gear shaft torque in real time
CN104699905B (en) Governing system gear drive identification modeling method based on frequency domain response

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