CN106706282A - Rotating machine fault diagnosis method based on Fourier decomposition - Google Patents

Rotating machine fault diagnosis method based on Fourier decomposition Download PDF

Info

Publication number
CN106706282A
CN106706282A CN201610964018.2A CN201610964018A CN106706282A CN 106706282 A CN106706282 A CN 106706282A CN 201610964018 A CN201610964018 A CN 201610964018A CN 106706282 A CN106706282 A CN 106706282A
Authority
CN
China
Prior art keywords
signal
fourier
frequency
time
fibfs
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.)
Pending
Application number
CN201610964018.2A
Other languages
Chinese (zh)
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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201610964018.2A priority Critical patent/CN106706282A/en
Publication of CN106706282A publication Critical patent/CN106706282A/en
Pending legal-status Critical Current

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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a rotating machine fault diagnosis method based on Fourier decomposition, comprising the following steps: (1) getting an acoustic emission signal through a rubbing acoustic emission test device; (2) defining a Fourier decomposition method and a Fourier natural frequency band function; (3) adaptively searching for a minimum number of Fourier natural frequency band functions within the Fourier frequency domain of the acoustic emission signal; and (4) analyzing the Fourier natural frequency band functions to get a time-frequency energy spectrum, and judging whether there is a fault according to the change of the time-frequency energy spectrum of the signal. According to the invention, a non-stationary nonlinear signal is decomposed into a number of Fourier natural frequency band functions and a residual component, the residual component is almost zero, and the signal is well divided. Whether a rotating machine has a rubbing fault is judged effectively, and the time and evolution trend of fault can be judged. The time-frequency resolution is high. There is neither endpoint effect nor modal aliasing phenomenon. Based on the Fourier transform theory, the method has a solid theoretical foundation and a complete mathematical model.

Description

A kind of rotary machinery fault diagnosis method based on Fourier decomposition
Technical field
The present invention relates to rotary machinery fault diagnosis field, especially a kind of rotating machinery fault based on Fourier decomposition Diagnostic method.
Background technology
Time frequency analysis are the powerfuls for processing time-domain signal.Conventional Time-Frequency Analysis Method mainly have Hilbert convert, Short Time Fourier Transform, Winger-Ville distributions, wavelet transformation, S-transformation and empirical mode decomposition (Empirical Mode Decomposition, EMD) etc..Because time-frequency characteristics are obtained from the formula of time-frequency distributions, explain comprising many thin The labyrinth of signal is saved, when the analysis of non-stationary signal, it is difficult to obtain better effects.Hilbert conversion can be very The good temporal characteristics for providing signal, there is temporal resolution very high, but the BORDER PROCESSING error of signal is larger;In short-term in Fu Leaf transformation is influenceed by window selection can not be while reach the optimal of temporal resolution and frequency resolution;Winger-Ville points Though cloth has time frequency resolution very high, for multicomponent data processing, due to the presence of cross term interference, its application is received very Big limitation;Wavelet transformation can provide preferable time precision, and its low frequency range has resolution ratio very high, but the resolution of high frequency region Rate is but very weak, and the phase information of part can make it lose physical significance.S-transformation is combined with Short Time Fourier Transform and small echo Conversion, solves the phase Localization Problems of wavelet transformation, but caused due to implying a window when multifrequency signal is differentiated Portray inaccurate in the place of frequency discontinuity, and be difficult to accurately portray amplitude.Norden E Huang etc. were carried in 1998 Gone out EMD, the method can carry out adaptive decomposition to non-stationary signal, but the method still suffers from certain defect, such as it It is a kind of empirical method, still lacks complete theoretical foundation, the orthogonality of resulting modal components still needs after decomposition Demonstration;It passes through the no science of end criterion that envelope is decomposed, and mould occurs in the signal that may result in decomposition The problem of state aliasing;Need just to can obtain complete IMF components by successive ignition, time-consuming, computationally intensive.For EMD's Defect, scholars propose the methods such as set empirical mode decomposition, multiple empirical mode decomposition, compact empirical mode decomposition, this A little methods overcome the problem that the time scale of EMD methods presence is separate to a certain extent, reduce modal overlap, end points effect Should, but the problem of EMD can not be fully solved.
The content of the invention
The technical problems to be solved by the invention are, there is provided a kind of rotary machinery fault diagnosis based on Fourier decomposition Method, can effectively judge the generation of rotating machinery fault, and time frequency resolution is high, in the absence of end effect and modal overlap Phenomenon.
In order to solve the above technical problems, the present invention provides a kind of rotary machinery fault diagnosis side based on Fourier decomposition Method, comprises the following steps:
(1) acoustic emission signal is obtained by touching the Experimental on acoustic emission device that rubs;
(2) Fourier decomposition method and the intrinsic frequency band function of Fourier are defined;
(3) in acoustic emission signal Fourier's frequency domain self adaptation search minimum number the intrinsic frequency band function of Fourier;
(4) time-frequency energy spectrum is obtained by parsing the intrinsic frequency band function of Fourier, by the change of signal time-frequency energy spectrum Carry out the generation of failure judgement.
Preferably, in step (2), it is y to define the intrinsic frequency band function of Fourieri(t)∈C[a, b], it meets following bar Part:
(1) FIBFs is zero-mean function, i.e.,
(2) FIBFs meets orthogonality, i.e.,
(3) instantaneous frequency and instantaneous amplitude of the analytical function AFIBFs of FIBFs be always greater than equal to 0, i.e., Wherein,
Signal decomposition is FIBFs and residual components of minimum number by FDM, is expressed as follows:
Wherein r (t) represents residual error component, yiT () is the M one-component signal of non-stationary, i.e. FIBFs, FIBFs is instantaneous Frequency has the single group sub-signal of physical significance.
Preferably, in step (3), if signal x (t) is interval [t1,t1+T0] the interior real-valued letter for following dirichlet condition Number, x (t) is configured to periodic signalI.e.Work as t1≤t≤t+T0When, ω (t)=1;When t≤ t1,t1+T0During≤t, ω (t)=0, thenFourier expansion formula be:
Wherein frequency is
Formula (3-1) can be write as following expression:
Wherein, ck=ak-jbk,ck *=ak+jbk, obtained by formula (3-2):
Analytical function is For Conjugate complex number,ForReal part, therefore can be byIt is written as:
Following formula can be obtained by formula (3-4) and formula (3-5):
When signal is decomposed, it is thus necessary to determine that maximum NiValue obtains the parsing FIBFs of minimum number;FDM methods Defined in two kinds of ways of search, i.e., forward search AFIBFs (from low to high search for) and reverse search AFIBFs (from high frequency To low-frequency acquisition), thus in the presence of two kinds of time-frequency Energy distributions, according to the difference of signal, two kinds of isolations may disclose identical Or two kinds of different characteristic types, two kinds of algorithms are summarized as follows:
Define to the increased direction search of i to search for AFIBFs forward, i.e., search for (LTH-FS) from low to high;To every I layers, from (Ni-1+ 1) start, and gradually increase NiTo maximum ∞, i.e. (Ni-1+1)≤Ni≤ ∞ andI is since 1:
It can thus be concluded that,
Wherein, N0=0, NM=∞,
Similar, in the reverse search of AFIBFs (from high frequency to low-frequency acquisition, HTL-FS), formula is limited under (13) Limit will become from k=NiTo (Ni-1- 1), corresponding N0=∞, NM=1. at this moment from (Ni-1- 1) start to look for, gradually reduce Ni's Value untill when selecting the FIBFs of minimum number, 1≤Ni≤(Ni- 1), i=1 ..., M;I.e.:
Beneficial effects of the present invention are:The intrinsic frequency band of Fourier of the search minimal amount of self adaptation in Fourier's frequency domain Non-stationary nonlinear properties are decomposed into the intrinsic frequency band function of several Fourier and a residual components, residual components by function Almost 0, signal is divided well;The generation of rotating machinery bump-scrape failure can effectively be judged, can determine whether that failure is sent out Raw time and evolving trend, and time frequency resolution is high, in the absence of end effect and modal overlap phenomenon, is become based on Fourier Theory is changed, there is solid theoretical foundation and complete Mathematical Modeling.
Brief description of the drawings
Fig. 1 (a) is experiment porch simulation drawing of the invention.
Fig. 1 (b) is experiment porch structural representation of the invention.
Fig. 2 is emulation time domain plethysmographic signal figure of the invention.
Fig. 3 decomposes (LTH) schematic diagram for emulation signal FDM of the invention.
Fig. 4 decomposes (HTL) schematic diagram for emulation signal FDM of the invention.
Fig. 5 is emulation signal EMD decomposing schematic representations of the invention.
Fig. 6 is emulation signal FDM time-frequencies energy spectrum (LTH) schematic diagram of the invention.
Fig. 7 is emulation signal FDM time-frequencies energy spectrum (HTL) schematic diagram of the invention.
Fig. 8 is emulation signal EMD time-frequency energy spectrums of the invention.
Fig. 9 touches the AE time domain plethysmographic signal figures that rub for nothing of the invention.
Figure 10 touches the AE signal spectrum figures that rub for nothing of the invention.
Figure 11 touches AE signal FDM time-frequencies energy spectrum (LTH-FS) schematic diagram that rubs for nothing of the invention.
Figure 12 touches AE signal FDM time-frequencies energy spectrum (HTL-FS) schematic diagram that rubs for nothing of the invention.
Figure 13 touches the AE signal EDM time-frequency energy spectrums that rub for nothing of the invention.
Figure 14 touches the AE time domain plethysmographic signal figures that rub for of the invention.
Figure 15 touches the AE signal spectrum figures that rub for of the invention.
Figure 16 touches AE signal FDM time-frequencies energy spectrum (LTH-FS) schematic diagram that rubs for of the invention.
Figure 17 touches AE signal FDM time-frequencies energy spectrum (HTL-FS) schematic diagram that rubs for of the invention.
Figure 18 touches the AE signal EDM time-frequency energy spectrums that rub for of the invention.
Specific embodiment
As shown in figure 1, touching acoustic emission experiment system of rubbing by rotor rubbing experimental bench, sensor, preamplifier, speed regulator Constituted with sound emission acquisition system.Rotor rubbing testing stand is flex rotor testing stand, is used to support having for rotor by three The bearing block of sliding bearing, two touch rub disk and touch mount screw composition.Touching mount screw can be by screw on lid guided wave plate Spindle central is pointed to, and is in contact with disk side.When rotor is rotated with certain rotating speed, regulation is touched mount screw and touches the disk that rubs Touch and rub, the bump-scrape acoustic emission signal of generation is received via guided wave plate by acoustic emission sensor.Mount screw is touched by regulation Screw-in depth is rubbed simulating touching for varying strength, and speed regulator realizes the stepless time adjustment of motor 0-10000r/min scopes.In order to drop The low acoustic emission waveform caused because medium is discontinuous is distorted, and couplant is filled between contact surface.The material of lid guided wave plate Expect to be No. 45 steel.
Acoustic emission signal acquisition system is that, by Polar9300e portable industrial pcs, built-in PCI-2 sound emissions capture card is adopted UT-1000 sensors are used, it is 1MbPS to set AE signal sampling frequencies, and sampled point is 20000, and filtered band is set to 0- 200kHz, gain amplifier is 40db.In order to reduce the aliasing of LAMB ripples and boundary echo, the signal acquisition in experiment is touched and rubbed just The signal of phase.
A kind of rotary machinery fault diagnosis method based on Fourier decomposition, comprises the following steps:
(1) acoustic emission signal is obtained by touching the Experimental on acoustic emission device that rubs;
(2) Fourier decomposition method and the intrinsic frequency band function of Fourier are defined;
(3) in acoustic emission signal Fourier's frequency domain self adaptation search minimum number the intrinsic frequency band function of Fourier;
(4) time-frequency energy spectrum is obtained by parsing the intrinsic frequency band function of Fourier, by the change of signal time-frequency energy spectrum Carry out the generation of failure judgement.
In step (2), it is y to define the intrinsic frequency band function of Fourieri(t)∈C[a, b], it meets following condition:
(1) FIBFs is zero-mean function, i.e.,
(2) FIBFs meets orthogonality, i.e.,
(3) instantaneous frequency and instantaneous amplitude of the analytical function AFIBFs of FIBFs be always greater than equal to 0, i.e., Wherein,
Signal decomposition is FIBFs and residual components of minimum number by FDM, is expressed as follows:
Wherein r (t) represents residual error component, yiT () is the M one-component signal of non-stationary, i.e. FIBFs, FIBFs is instantaneous Frequency has the single group sub-signal of physical significance.
In step (3), if signal x (t) is interval [t1,t1+T0] the interior real-valued signal for following dirichlet condition, by x T () is configured to periodic signalI.e.Work as t1≤t≤t+T0When, ω (t)=1;As t≤t1,t1+ T0During≤t, ω (t)=0, thenFourier expansion formula be:
Wherein frequency is
Formula (3-1) can be write as following expression:
Wherein, ck=ak-jbk,ck *=ak+jbk, obtained by formula (3-2):
Analytical function is For Conjugate complex number,ForReal part, therefore can be byIt is written as:
Following formula can be obtained by formula (3-4) and formula (3-5):
When signal is decomposed, it is thus necessary to determine that maximum NiValue obtains the parsing FIBFs of minimum number;FDM methods Defined in two kinds of ways of search, i.e., forward search AFIBFs (from low to high search for) and reverse search AFIBFs (from high frequency To low-frequency acquisition), thus in the presence of two kinds of time-frequency Energy distributions, according to the difference of signal, two kinds of isolations may disclose identical Or two kinds of different characteristic types, two kinds of algorithms are summarized as follows:
Define to the increased direction search of i to search for AFIBFs forward, i.e., search for (LTH-FS) from low to high;To every I layers, from (Ni-1+ 1) start, and gradually increase NiTo maximum ∞, i.e. (Ni-1+1)≤Ni≤ ∞ andI is since 1:
It can thus be concluded that,
Wherein, N0=0, NM=∞,
Similar, in the reverse search of AFIBFs (from high frequency to low-frequency acquisition, HTL-FS), formula is limited under (13) Limit will become from k=NiTo (Ni-1- 1), corresponding N0=∞, NM=1. at this moment from (Ni-1- 1) start to look for, gradually reduce Ni's Value untill when selecting the FIBFs of minimum number, 1≤Ni≤(Ni- 1), i=1 ..., M;I.e.:
EMD decomposition is carried out to signal, and carries out Hilbert conversion, obtain time-frequency energy diagram.
Emulation signal x (t) is built, the signal is by a linear FM signal x1(t) and FMAM signal x2(t) Composition.
x1(t)=cos (400 π t+200 π t2);
x2(t)=cos (100 π t+0.4 π sin (20 π t))+0.1sin (10 π t) cos (100 π t+0.4 π sin (20 π t));
X (t)=x1(t)+x2(t)
Fig. 2 be emulate signal time domain waveform, Fig. 3 be emulate signal through the LTH algorithms of FDM decomposition result, Fig. 4 be through The decomposition result of the HTL algorithms of FDM, Fig. 5 is to emulate the result that signal is decomposed through EMD.As can be seen that emulation from Fig. 3 and Fig. 4 Signal is effectively decomposed into two FIBFs and residual components, and the two components correspond respectively to linear FM signal And amplitude-modulation frequency-modulation signal, residual components are zero.From figure 5 it can be seen that emulation signal decomposites 7 rank IMF components and one through EMD Individual residual error, the natural mode of vibration component for decompositing and emulation signal are not corresponded and there is serious modal overlap and end effect, It can thus be seen that the effect that EMD is decomposed is undesirable.Fig. 6 is the time-frequency energy spectrum that signal is decomposed by the LTH algorithms of FDM, and Fig. 7 is The time-frequency energy spectrum that signal is decomposed through the HTL algorithms of FDM, Fig. 8 is the time-frequency energy spectrum that signal is decomposed through EMD.Comparison diagram 6~8 As can be seen that Fig. 6 is that, through two kinds of time-frequency energy spectrums of different searching methods decomposition of FDM, its time frequency resolution is higher with Fig. 7, LTH algorithms are higher than HTL-FS algorithms in HFS time frequency resolution, and HTL-FS algorithms are higher than in low frequency part time frequency resolution HTL-FS algorithms.Fig. 8 can be seen that through EMD decompose time-frequency energy spectrum resolution ratio it is relatively low, there is false mode and modal overlap Phenomenon.Therefore, FDM algorithms overcome the problem of EMD presence to a certain extent.
UT-1000 sensors, AE signals is used to be gathered by the industrial computer of built-in PCI-2 sound emissions capture card in experiment, if AE signal sampling frequencies are put for 1MbPS, sampled point is 5120, and filtered band is set to 0-200kHz, and gain amplifier is 40db.For Reduce the aliasing of LAMB ripples and boundary echo, the signal acquisition in experiment touches the signal at the initial stage of rubbing, respectively rubbing without touching and Touch under the state of rubbing and measure acoustic emission signal.Time domain plethysmographic signal and spectrogram are rubbed without touching as shown in Figure 9 and Figure 10.Figure 11, figure 12 are respectively without the AE signals that rub are touched using FDM-LTH algorithms and the time-frequency energy spectrum of FDMM-HTL algorithms, and Figure 13 is that nothing touches the AE that rubs Signal is through the time-frequency power spectrum after EMD.Contrast two methods time-frequency energy spectrum, without touch rub AE signals through FDM LTH algorithms and The time-frequency energy spectrum of HTL algorithms is consistent, and signal is decomposed well, and has time frequency resolution very high.Figure 13 It is evident that time frequency resolution is low, and has modal overlap phenomenon, effect is not so good as Figure 11 and Figure 12.Figure 15 believes to touch the AE that rubs Number spectrogram, as can be seen from the figure the energy of the AE signals be concentrated mainly within 100kHz.Figure 16, Figure 17 have been respectively Touch time-frequency energy spectrum of the AE signals using FDM-LTH algorithms and FDM-HTL algorithms that rub, Figure 18 rubs AE signals after EMD to touch Time-frequency energy spectrum.Time-frequency distributions in Figure 16 and Figure 17 are consistent, and with resolution ratio very high, as can be seen from the figure have There is the obvious frequency distribution of energy feature, thus can determine whether out faulty generation.The frequency component of fault-signal is can be seen that simultaneously Energy concentrate on time leading portion, and decay with the time, it can be seen that this is transient state impact-rub malfunction.Figure 18 time frequency resolutions compared with It is low, it is distributed more disorderly, modal overlap is more, and its frequency distribution can not be corresponding with the spectrogram of Figure 15, hardly results in this and touches and rubs The time-frequency characteristic of AE signals.FDM can be seen that by above-mentioned experimental analysis bright in terms of the sensitivity that impact-rub malfunction occurs is judged It is aobvious to be better than EMD methods.
Although the present invention is illustrated and has been described with regard to preferred embodiment, it is understood by those skilled in the art that Without departing from scope defined by the claims of the present invention, variations and modifications can be carried out to the present invention.

Claims (3)

1. a kind of rotary machinery fault diagnosis method based on Fourier decomposition, it is characterised in that comprise the following steps:
(1) acoustic emission signal is obtained by touching the Experimental on acoustic emission device that rubs;
(2) Fourier decomposition method and the intrinsic frequency band function of Fourier are defined;
(3) in acoustic emission signal Fourier's frequency domain self adaptation search minimum number the intrinsic frequency band function of Fourier;
(4) time-frequency energy spectrum is obtained by parsing the intrinsic frequency band function of Fourier, is sentenced by the change of signal time-frequency energy spectrum The generation of disconnected failure.
2. the rotary machinery fault diagnosis method of Fourier decomposition is based on as claimed in claim 1, it is characterised in that step (2) in, it is y to define the intrinsic frequency band function of Fourieri(t)∈C[a, b], it meets following condition:
(1) FIBFs is zero-mean function, i.e.,
(2) FIBFs meets orthogonality, i.e.,
(3) instantaneous frequency and instantaneous amplitude of the analytical function AFIBFs of FIBFs be always greater than equal to 0, i.e., Wherein,
Signal decomposition is FIBFs and residual components of minimum number by FDM, is expressed as follows:
x ( t ) = Σ i = 1 M y i ( t ) + r ( t ) - - - ( 2 - 1 )
Wherein r (t) represents residual error component, yiT () is the M one-component signal of non-stationary, i.e. FIBFs, FIBFs is instantaneous frequency Single group sub-signal with physical significance.
3. the rotary machinery fault diagnosis method of Fourier decomposition is based on as claimed in claim 1, it is characterised in that step (3) in, if signal x (t) is interval [t1,t1+T0] the interior real-valued signal for following dirichlet condition, x (t) is configured to the cycle SignalI.e.Work as t1≤t≤t+T0When, ω (t)=1;As t≤t1,t1+T0During≤t, ω (t) =0, thenFourier expansion formula be:
x T 0 ( t ) = a 0 + Σ k = 1 ∞ [ a k c o s ( kω 0 t ) + b k s i n ( kω 0 t ) ] - - - ( 3 - 1 )
Wherein frequency is
a k = 2 T 0 ∫ t 1 t 1 + T 0 x T 0 ( t ) c o s ( kω 0 t ) d t
b k = 2 T 0 ∫ t 1 t 1 + T 0 x T 0 ( t ) s i n ( kω 0 t ) d t
Formula (3-1) can be write as following expression:
x T 0 ( t ) = a 0 + 1 2 Σ k = 1 ∞ [ c k exp ( jkω 0 t ) + c k * exp ( - jkω 0 t ) ] - - - ( 3 - 2 )
Wherein, ck=ak-jbk,ck *=ak+jbk, obtained by formula (3-2):
Analytical function is
ForConjugate complex number,ForReal part, therefore can WillIt is written as:
Following formula can be obtained by formula (3-4) and formula (3-5):
Σ k = 1 M a i ( t ) exp ( jφ i ( t ) ) = Σ k = 1 ∞ c k exp ( jkω 0 t ) - - - ( 3 - 6 )
When signal is decomposed, it is thus necessary to determine that maximum NiValue obtains the parsing FIBFs of minimum number;Defined in FDM methods Two kinds of ways of search, i.e., search for forward AFIBFs and reverse search AFIBFs, thus in the presence of two kinds of time-frequency Energy distributions, according to letter Number difference, two kinds of isolations may disclose identical or two kinds of different characteristic types, and two kinds of algorithms are summarized as follows:
Define to the increased direction search of i to search for AFIBFs forward, i.e., search for (LTH-FS) from low to high;To every i layers, From (Ni-1+ 1) start, and gradually increase NiTo maximum ∞, i.e. (Ni-1+1)≤Ni≤ ∞ andI is since 1:
It can thus be concluded that,
Wherein, N0=0, NM=∞,
Similar, in the reverse search of AFIBFs (from high frequency to low-frequency acquisition, HTL-FS), the upper limit is limited under formula (13) will Can become from k=NiTo (Ni-1- 1), corresponding N0=∞, NM=1. at this moment from (Ni-1- 1) start to look for, gradually reduce NiValue it is straight To when selecting the FIBFs of minimum number, 1≤Ni≤(Ni- 1), i=1 ..., M;I.e.:
CN201610964018.2A 2016-11-04 2016-11-04 Rotating machine fault diagnosis method based on Fourier decomposition Pending CN106706282A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610964018.2A CN106706282A (en) 2016-11-04 2016-11-04 Rotating machine fault diagnosis method based on Fourier decomposition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610964018.2A CN106706282A (en) 2016-11-04 2016-11-04 Rotating machine fault diagnosis method based on Fourier decomposition

Publications (1)

Publication Number Publication Date
CN106706282A true CN106706282A (en) 2017-05-24

Family

ID=58940730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610964018.2A Pending CN106706282A (en) 2016-11-04 2016-11-04 Rotating machine fault diagnosis method based on Fourier decomposition

Country Status (1)

Country Link
CN (1) CN106706282A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108731945A (en) * 2018-08-02 2018-11-02 南昌航空大学 A kind of extracting method of aircraft engine rotor system fault-signal characteristic information
CN109034042A (en) * 2018-07-20 2018-12-18 东北大学 Based on the double synchronous nonstationary random response methods for extracting transformation of generalized chirp
CN110611630A (en) * 2019-08-09 2019-12-24 西安电子科技大学 2-FSK awakening receiver of tuning radio frequency architecture and demodulation method thereof
CN113239613A (en) * 2021-04-09 2021-08-10 国网新源控股有限公司 Hydro-turbine set throw alarm judgment method
CN117825520A (en) * 2024-03-05 2024-04-05 中国矿业大学(北京) Detection method, device, medium and electronic equipment for detecting object damage
CN113239613B (en) * 2021-04-09 2024-05-31 国网新源控股有限公司 Method for judging swing degree alarm of water turbine unit

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06258198A (en) * 1993-03-02 1994-09-16 Hitachi Ltd Method and apparatus for fault diagnosis
CN101251445A (en) * 2008-04-16 2008-08-27 邓艾东 Method for analysis of fractal characteristic of rotating machinery bump-scrape acoustic emission signal
CN101435799A (en) * 2008-12-19 2009-05-20 清华大学 Failure diagnosis method and apparatus of hydroturbine based on acoustic emission technology
CN102052998A (en) * 2010-11-09 2011-05-11 东南大学 Rotor rub-impact acoustic emission signal recognition method
WO2015012124A1 (en) * 2013-07-22 2015-01-29 Ntn株式会社 Fault diagnosis device and fault diagnosis method for wind power generation device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06258198A (en) * 1993-03-02 1994-09-16 Hitachi Ltd Method and apparatus for fault diagnosis
CN101251445A (en) * 2008-04-16 2008-08-27 邓艾东 Method for analysis of fractal characteristic of rotating machinery bump-scrape acoustic emission signal
CN101435799A (en) * 2008-12-19 2009-05-20 清华大学 Failure diagnosis method and apparatus of hydroturbine based on acoustic emission technology
CN102052998A (en) * 2010-11-09 2011-05-11 东南大学 Rotor rub-impact acoustic emission signal recognition method
WO2015012124A1 (en) * 2013-07-22 2015-01-29 Ntn株式会社 Fault diagnosis device and fault diagnosis method for wind power generation device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
PUSHPENDRA SINGH等: "The Fourier Decomposition Method for nonlinear and nonstationary time series analysis", 《SIGNAL》 *
张占一等: "基于声振信号EMD分解的轻微碰摩故障诊断方法研究", 《振动与冲击》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109034042A (en) * 2018-07-20 2018-12-18 东北大学 Based on the double synchronous nonstationary random response methods for extracting transformation of generalized chirp
CN109034042B (en) * 2018-07-20 2021-10-01 东北大学 Non-stationary signal processing method based on generalized linear frequency modulation double-synchronous extraction transformation
CN108731945A (en) * 2018-08-02 2018-11-02 南昌航空大学 A kind of extracting method of aircraft engine rotor system fault-signal characteristic information
CN110611630A (en) * 2019-08-09 2019-12-24 西安电子科技大学 2-FSK awakening receiver of tuning radio frequency architecture and demodulation method thereof
CN113239613A (en) * 2021-04-09 2021-08-10 国网新源控股有限公司 Hydro-turbine set throw alarm judgment method
CN113239613B (en) * 2021-04-09 2024-05-31 国网新源控股有限公司 Method for judging swing degree alarm of water turbine unit
CN117825520A (en) * 2024-03-05 2024-04-05 中国矿业大学(北京) Detection method, device, medium and electronic equipment for detecting object damage

Similar Documents

Publication Publication Date Title
He et al. Time-reassigned synchrosqueezing transform: The algorithm and its applications in mechanical signal processing
Han et al. Fault feature extraction of low speed roller bearing based on Teager energy operator and CEEMD
Zhang et al. Weak feature enhancement in machinery fault diagnosis using empirical wavelet transform and an improved adaptive bistable stochastic resonance
CN106568607A (en) Rub-impact sound emission fault diagnosis method based on empirical wavelet transformation
Cui et al. Matching pursuit of an adaptive impulse dictionary for bearing fault diagnosis
CN106706282A (en) Rotating machine fault diagnosis method based on Fourier decomposition
Li et al. Bearing fault detection and diagnosis based on order tracking and Teager-Huang transform
CN102323480B (en) Electric energy quality analyzing method based on Hilbert-Huang transform
CN107631877A (en) A kind of rolling bearing fault collaborative diagnosis method for casing vibration signal
Cao et al. Zoom synchrosqueezing transform and iterative demodulation: Methods with application
Zhou et al. Study on intra-wave frequency modulation phenomenon in detection of rub-impact fault
Ming et al. Fault feature extraction and enhancement of rolling element bearing in varying speed condition
CN106568589A (en) Rubbing acoustic emission denoise method based on empirical wavelet transform
CN103018044A (en) Bearing combined failure diagnosis method based on improved impact dictionary matching pursuit
Wei et al. Time-varying envelope filtering for exhibiting space bearing cage fault features
CN109520738A (en) Rotating machinery Fault Diagnosis of Roller Bearings based on order spectrum and envelope spectrum
Cheng et al. Adaptive sparsest narrow-band decomposition method and its applications to rolling element bearing fault diagnosis
CN111256993A (en) Method and system for diagnosing fault type of main bearing of wind turbine generator
Fang et al. Comparison of EMD and EEMD in rolling bearing fault signal analysis
Xu et al. An adaptive spectrum segmentation method to optimize empirical wavelet transform for rolling bearings fault diagnosis
Han et al. Second-order multisynchrosqueezing wavelet transform for bearing fault detection
Li et al. Revision of bearing fault characteristic spectrum using LMD and interpolation correction algorithm
Wu et al. Application of multi-base fusion generalized chirplet basis transform in vibration signal analysis of multiple rotor rotating machinery
CN105303033A (en) Rolling bearing fault diagnosis method based on integral inherent time scale decomposition algorithm
Li et al. Gear fault detection based on Teager-Huang transform

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Deng Aidong

Inventor after: Zhang Rui

Inventor after: Si Xiaodong

Inventor after: Wang Sheng

Inventor after: Liu Dongying

Inventor after: Li Jing

Inventor after: Zhu Jing

Inventor after: Di Yimeng

Inventor before: Deng Aidong

Inventor before: Zhang Rui

Inventor before: Si Xiaodong

Inventor before: Liu Dongying

Inventor before: Li Jing

CB03 Change of inventor or designer information
RJ01 Rejection of invention patent application after publication

Application publication date: 20170524

RJ01 Rejection of invention patent application after publication