CN103900815A - Rolling bearing fault diagnosis method based on EEMD and distribution fitting testing - Google Patents
Rolling bearing fault diagnosis method based on EEMD and distribution fitting testing Download PDFInfo
- Publication number
- CN103900815A CN103900815A CN201410130757.2A CN201410130757A CN103900815A CN 103900815 A CN103900815 A CN 103900815A CN 201410130757 A CN201410130757 A CN 201410130757A CN 103900815 A CN103900815 A CN 103900815A
- Authority
- CN
- China
- Prior art keywords
- eemd
- fault
- signal
- fault diagnosis
- rolling bearing
- 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
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a rolling bearing fault diagnosis method based on EEMD and distribution fitting testing, and relates to the field of rotary machine fault diagnosis. The method includes the steps that original signals are decomposed into a series of IMF components by means of EEMD; then each IMF component is sampled to obtain a sampling point; testing is carried out through a normal probability plot and a Jarque-Bera test, and whether data accord with normal distribution or not is judged; the IMF components according with white noise characteristics are removed, and signals after other components are added and noise reduction is carried out are reserved; finally, fault diagnosis is carried out on rolling bearings by means of an envelope spectrum. According to the method, the aim of signal noise reduction can be well achieved, signal fault features can be more obvious, ultimately fault frequency is identified by means of the envelope spectrum, and fault types of the rolling bearings can be well diagnosed.
Description
Technical field
The present invention relates to rotary machinery fault diagnosis field, in particular a kind of Fault Diagnosis of Roller Bearings based on EEMD and fitting of distribution check.
Background technology
Rolling bearing is widely used in rotating machinery, and it is also the critical component of large rotating machinery.Due to its architectural characteristic, rolling bearing is to hold flimsy part, and in the polytype fault occurring at rotating machinery, many all have relation with the damage of rolling bearing.The fault of plant equipment 70% is vibration fault, in these vibration faults, approximately has 30% to be caused by bearing.This is because the working environment very severe of rolling bearing, once rolling bearing breaks down, will cause a series of chain reactions of whole plant equipment, brings huge potential safety hazard to plant equipment.Therefore,, in mechanical fault diagnosis technology, the fault diagnosis technology of rolling bearing is important ingredient.
The method for diagnosing faults of rolling bearing has a lot, and rolling bearing fault diagnosis technology based on vibration signal is the most frequently used, is also effective method.Vibration signal has comprised abundant plant equipment extremely or failure message, and it can reflect the status flag that plant equipment is moved.Traditional most widely used Digital Signal Processing is exactly Fourier transform.Invent FFT(Fast Fourier Transform (FFT) from nineteen sixty-five Cooley-Tu Ji, FastFourier Transformation) afterwards, Spectral Analysis Method is all widely used in every field.Wavelet transformation is the novel Digital Signal Processing growing up in the later stage eighties, and its essence is a Fourier transform that window is adjustable, and its time window and frequency window can change, and is a kind of analytical approach of time frequency localization.FFT can not provide frequency situation over time, can only provide separately the frequency distribution situation of signal, lacks the positioning function with the time.Wavelet transformation, due to wavelet basis function limited length, therefore, concerning the distribution of signal energy relative time and frequency, carries out accurate analysis to it and will produce larger difficulty.Once wavelet base is chosen, just must analyze all data to be analyzed with this wavelet basis function, therefore wavelet transformation does not possess adaptive feature.
EMD(empirical mode decomposition, empirical mode decomposition) can be limited intrinsic mode function component a complicated signal decomposition according to local time's characteristic of signal.But EMD also has shortcoming, topmost problem is exactly mode Aliasing Problem.EEMD(gathers empirical mode decomposition, ensemble empirical mode decomposition) be a kind of novel digital signal processing method, can be by a series of being decomposed into of complex fault signal adaptive intrinsic mode function components, be applicable to analysis and processing non-linear, non-stationary signal.Fitting of distribution check is whether sample value for testing one group of data is from the method for inspection of certain given distribution.
Summary of the invention
Technical matters to be solved by this invention is for the deficiencies in the prior art, and a kind of Fault Diagnosis of Roller Bearings based on EEMD and fitting of distribution check is provided.
Technical scheme of the present invention is as follows:
Based on a Fault Diagnosis of Roller Bearings for EEMD and fitting of distribution check, its step is as follows:
(1) utilize EEMD that original signal is decomposed into a series of IMF components;
(2) each IMF component is sampled, obtain sampled point; Then utilize normal probability paper figure and Jarque-Bera test to test, judge whether data meet normal distribution;
(3) remove the IMF component that meets white noise feature, retain after all the other components are added and obtain the signal after noise reduction;
(4) utilize envelope spectrum to carry out fault diagnosis to rolling bearing.
The present invention adopts EEMD to carry out the decomposition of sophisticated signal, can solve classic method and not possess adaptive shortcoming, also can alleviate the problem of the mode aliasing that EMD occurs, is more applicable for analysis pre-service non-linear, non-stationary signal; Adopt fitting of distribution check method of inspection, screen useful IMF component from the probability density characteristics of white noise IMF component, remove useless IMF component simultaneously, have good theoretical foundation; The present invention can reach the object of good signal de-noising, makes the fault signature of signal more obvious, finally utilizes envelope spectrum to identify failure-frequency, can well diagnose out the fault type of bearing.
Accompanying drawing explanation
Fig. 1 is that s (t) utilizes EMD decomposition result figure.
Fig. 2 is that s (t) utilizes EEMD decomposition result figure.
Fig. 3 is the decomposition result figure of inner ring fault-signal through EEMD.
Fig. 4 is the normal probability paper figure of the IMF component of inner ring fault-signal.
Fig. 5 is the result after inner ring fault-signal noise reduction.
Fig. 6 is inner ring fault-signal envelope spectrum.
Fig. 7 is the decomposition result figure of outer ring fault-signal through EEMD.
Fig. 8 is the normal probability paper figure of the IMF component of outer ring fault-signal.
Fig. 9 is the result after the fault-signal noise reduction of outer ring.
Figure 10 is outer ring fault-signal envelope spectrum.
Embodiment
Below in conjunction with specific embodiment, the present invention is described in detail.
1, suppose a simulate signal s (t), its formation is the stack of 7Hz sinusoidal component and Gauspuls pulse component, and Fig. 1 is the result after s (t) decomposes through EMD.
Obviously can find out from Fig. 1, clearly there is mode aliasing in IMF1 and IMF2, has been difficult to find out IMF1 and the represented physical meaning of IMF2 from figure.
2, the simulate signal s (t) in 1 is carried out to EEMD decomposition, result as shown in Figure 2.As can be seen from Figure 2, s (t) has obtained 4 IMF components and 1 surplus after decomposing through EEMD, wherein c1 and c2 represent to decompose Gauspuls pulse component out, c4 expresses 7Hz sinusoidal component well, therefore, after decomposing through EEMD, can well decomposite the composition that original signal comprises, each IMF component has the physical significance of oneself.Therefore, EEMD can solve the problem of mode aliasing well, can well reflect the essential characteristic of s (t).
Rolling bearing inner ring fault-signal comes from bearing data center of CWRU of the U.S., and wherein: sample frequency is 12000Hz, the rotating speed of motor is 1750r/min.First rolling bearing inner ring fault-signal is carried out to EEMD decomposition, obtain some intrinsic mode function components.Then according to the probability density characteristics of IMF component, utilize normal probability paper figure and Jarque-Bera method of testing to carry out the removal of choosing of useful IMF component and useless IMF component.Useful IMF component is reconstructed to the inner ring fault-signal obtaining after noise reduction, finally utilizes envelope spectrum to identify Rolling Bearing Fault Character frequency.
Fig. 3 is the result of rolling bearing inner ring fault-signal after EEMD decomposes.Fig. 4 is the normal probability paper figure of the inner ring fault IMF component that decomposes out.
For normal probability paper figure, if sampled point ("+" point in figure) approaches the dotted line in figure, can think that this group sampled point is similar to Normal Distribution, therefore, for c4, c5, c6, its characteristic distributions meets the distribution character of the IMF of white noise, therefore can set it as useless IMF component and remove.In order to improve the accuracy of check, also adopt Jarque-Bera test herein, test result is as shown in table 1.
The Jarque-Bera test result of the IMF component of table 1 inner ring fault-signal
For Jarque-Bera test, there is a H value, in the time of H=0, show the hypothesis that this group data receiver sample is normal distribution, otherwise, in the time of H=1, the hypothesis that this group data rejection sample is normal distribution.Therefore, can remove equally c4, c5, c6 component, retains all the other components, consistent with the result that normal probability paper figure obtains.Therefore, reconstruct c1, c2, c3 and c7, as shown in Figure 5, the signal after noise reduction has not only well retained the primitive character of signal, and makes fault signature more obvious.Fig. 6 is the envelope spectrum of inner ring fault-signal, and this has many failure-frequencies, the computing formula according to Bearing Analysis Theory failure-frequency: inner ring failure-frequency f
i=0.5Zf (1+dcos α D), wherein, gyro frequency is f, and roller diameter is d, and the pitch diameter of bearing is D, and roller number is Z, bearing pressure angle is α.Rolling bearing turn frequently f=29.17Hz, inner ring failure-frequency is f
i=157.9Hz, as can be seen from Figure 6, f
1=29.3Hz, f
2=58.6Hz, f
3=158.2Hz, f
4=316.5Hz.Therefore, f
1and f
2be approximately equal to respectively rolling bearing turn frequently and 2 times turn frequency, f
3and f
4be approximately equal to respectively inner ring failure-frequency and 2 times of inner ring failure-frequencies of rolling bearing.There are some errors in the calculated results and results of calculation, is because the workmanship precision of rolling bearing is not accurate enough affected.
Housing washer fault-signal still comes from bearing data center of CWRU of the U.S., and wherein: sample frequency is 12000Hz, the rotating speed of motor is 1750r/min.First housing washer fault-signal is carried out to EEMD decomposition, obtain some intrinsic mode function components.Then according to the probability density characteristics of IMF component, utilize normal probability paper figure and Jarque-Bera method of testing to carry out the removal of choosing of useful IMF component and useless IMF component.Useful IMF component is reconstructed to the outer ring fault-signal obtaining after noise reduction, finally utilizes envelope spectrum to identify Rolling Bearing Fault Character frequency.
Fig. 7 is the result of housing washer fault-signal after EEMD decomposes.Fig. 8 is the normal probability paper figure of the outer ring fault IMF component that decomposes out.Table 2 is the Jarque-Bera test result of outer ring fault IMF component.
The Jarque-Bera test result of the IMF component of table 2 outer ring fault-signal
From Fig. 8 and table 2, can find out, c1, c4, c6 and c7 do not meet the probability density characteristics of white noise IMF, therefore, retain these components, remove remaining component simultaneously.Fig. 9 is the result after the failure reconfiguration of outer ring, and Figure 10 is outer ring fault-signal envelope spectrum.Computing formula according to Bearing Analysis Theory failure-frequency: outer ring failure-frequency f
o=0.5Zf (1-dcos α D), the represented implication of correlation parameter as shown in Example 2.Rolling bearing turn frequently f=29.17Hz, outer ring failure-frequency is f
o=104.57Hz, as can be seen from Figure 10, f
1=11.7Hz, f
2=29.3Hz, f
3=105.5Hz, f
4=211Hz.Therefore, f
1substantially equal the failure-frequency (11.6Hz) of retainer, f
2for turning frequently of rolling bearing, f
3and f
4substantially equal outer ring failure-frequency and 2 times of outer ring failure-frequencies.
Should be understood that, for those of ordinary skills, can be improved according to the above description or convert, and all these improvement and conversion all should belong to the protection domain of claims of the present invention.
Claims (1)
1. the Fault Diagnosis of Roller Bearings based on EEMD and fitting of distribution check, is characterized in that, its step is as follows:
(1) utilize EEMD that original signal is decomposed into a series of IMF components;
(2) each IMF component is sampled, obtain sampled point; Then utilize normal probability paper figure and Jarque-Bera test to test, judge whether data meet normal distribution;
(3) remove the IMF component that meets white noise feature, retain after all the other components are added and obtain the signal after noise reduction;
(4) utilize envelope spectrum to carry out fault diagnosis to rolling bearing.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410130757.2A CN103900815A (en) | 2014-04-02 | 2014-04-02 | Rolling bearing fault diagnosis method based on EEMD and distribution fitting testing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410130757.2A CN103900815A (en) | 2014-04-02 | 2014-04-02 | Rolling bearing fault diagnosis method based on EEMD and distribution fitting testing |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103900815A true CN103900815A (en) | 2014-07-02 |
Family
ID=50992279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410130757.2A Pending CN103900815A (en) | 2014-04-02 | 2014-04-02 | Rolling bearing fault diagnosis method based on EEMD and distribution fitting testing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103900815A (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155108A (en) * | 2014-07-21 | 2014-11-19 | 天津大学 | Rolling bearing failure diagnosis method base on vibration temporal frequency analysis |
CN104361242A (en) * | 2014-11-20 | 2015-02-18 | 青岛理工大学 | Bearing fault diagnosis method based on data driving and random intuition fuzzy strategy |
CN105092023A (en) * | 2015-08-11 | 2015-11-25 | 西安科技大学 | Bolt vibration signal correction method based on white noise statistical characteristics |
CN105973603A (en) * | 2016-06-29 | 2016-09-28 | 潍坊学院 | EEMD and rational spline smooth envelope analysis method for rotating machine |
CN106053061A (en) * | 2016-06-29 | 2016-10-26 | 潍坊学院 | Envelope analysis method based on non-linear mode decomposition and spectrum kurtosis |
CN106404396A (en) * | 2016-08-30 | 2017-02-15 | 四川中烟工业有限责任公司 | Rolling bearing fault diagnosis method |
CN106644484A (en) * | 2016-09-14 | 2017-05-10 | 西安工业大学 | Turboprop Engine rotor system fault diagnosis method through combination of EEMD and neighborhood rough set |
CN108446629A (en) * | 2018-03-19 | 2018-08-24 | 河北工业大学 | Rolling Bearing Fault Character extracting method based on set empirical mode decomposition and modulation double-spectrum analysis |
CN111855178A (en) * | 2020-07-23 | 2020-10-30 | 贵州永红航空机械有限责任公司 | Diagnosis method for running state of rotary product |
CN113390631A (en) * | 2021-06-15 | 2021-09-14 | 大连理工大学 | Fault diagnosis method for gearbox of diesel engine |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW201015345A (en) * | 2008-10-10 | 2010-04-16 | Univ Nat Central | Data decomposition method and computer system therefrom |
CN102840907A (en) * | 2012-09-18 | 2012-12-26 | 河南省电力公司电力科学研究院 | Rolling bearing vibration signal characteristic extracting and analyzing method under early fault state |
CN103091096A (en) * | 2013-01-23 | 2013-05-08 | 北京信息科技大学 | Extraction method for early failure sensitive characteristics based on ensemble empirical mode decomposition (EEMD) and wavelet packet transform |
-
2014
- 2014-04-02 CN CN201410130757.2A patent/CN103900815A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW201015345A (en) * | 2008-10-10 | 2010-04-16 | Univ Nat Central | Data decomposition method and computer system therefrom |
CN102840907A (en) * | 2012-09-18 | 2012-12-26 | 河南省电力公司电力科学研究院 | Rolling bearing vibration signal characteristic extracting and analyzing method under early fault state |
CN103091096A (en) * | 2013-01-23 | 2013-05-08 | 北京信息科技大学 | Extraction method for early failure sensitive characteristics based on ensemble empirical mode decomposition (EEMD) and wavelet packet transform |
Non-Patent Citations (2)
Title |
---|
LEI YAGUO, ET AL: "Application of the EEMD method to rotor fault diagnosis of rotating machinery", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》, no. 23, 24 November 2008 (2008-11-24), pages 1327 - 1338 * |
ZHANG YOUPENG, ET AL: "Improved Ensemble Empirical Mode Decomposition for Rolling Bearing Fault Diagnosis", 《TELKOMNIKA INDONESIAN JOURNAL OF ELECTRICAL ENGINEERING》, vol. 12, no. 1, 31 January 2014 (2014-01-31), pages 2 - 4 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155108A (en) * | 2014-07-21 | 2014-11-19 | 天津大学 | Rolling bearing failure diagnosis method base on vibration temporal frequency analysis |
CN104155108B (en) * | 2014-07-21 | 2017-07-07 | 天津大学 | A kind of Fault Diagnosis of Roller Bearings based on vibration time frequency analysis |
CN104361242B (en) * | 2014-11-20 | 2017-04-05 | 青岛理工大学 | Bearing fault diagnosis method based on data driving and random intuition fuzzy strategy |
CN104361242A (en) * | 2014-11-20 | 2015-02-18 | 青岛理工大学 | Bearing fault diagnosis method based on data driving and random intuition fuzzy strategy |
CN105092023A (en) * | 2015-08-11 | 2015-11-25 | 西安科技大学 | Bolt vibration signal correction method based on white noise statistical characteristics |
CN105973603A (en) * | 2016-06-29 | 2016-09-28 | 潍坊学院 | EEMD and rational spline smooth envelope analysis method for rotating machine |
CN106053061A (en) * | 2016-06-29 | 2016-10-26 | 潍坊学院 | Envelope analysis method based on non-linear mode decomposition and spectrum kurtosis |
CN105973603B (en) * | 2016-06-29 | 2018-03-13 | 潍坊学院 | The EEMD and rational spline smoothed envelope analysis method of a kind of rotating machinery |
CN106053061B (en) * | 2016-06-29 | 2018-03-13 | 潍坊学院 | A kind of envelope Analysis Method for decomposing and composing kurtosis based on nonlinear model |
CN106404396A (en) * | 2016-08-30 | 2017-02-15 | 四川中烟工业有限责任公司 | Rolling bearing fault diagnosis method |
CN106644484A (en) * | 2016-09-14 | 2017-05-10 | 西安工业大学 | Turboprop Engine rotor system fault diagnosis method through combination of EEMD and neighborhood rough set |
CN108446629A (en) * | 2018-03-19 | 2018-08-24 | 河北工业大学 | Rolling Bearing Fault Character extracting method based on set empirical mode decomposition and modulation double-spectrum analysis |
CN111855178A (en) * | 2020-07-23 | 2020-10-30 | 贵州永红航空机械有限责任公司 | Diagnosis method for running state of rotary product |
CN111855178B (en) * | 2020-07-23 | 2022-04-19 | 贵州永红航空机械有限责任公司 | Diagnosis method for running state of rotary product |
CN113390631A (en) * | 2021-06-15 | 2021-09-14 | 大连理工大学 | Fault diagnosis method for gearbox of diesel engine |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103900815A (en) | Rolling bearing fault diagnosis method based on EEMD and distribution fitting testing | |
CN105510023B (en) | Variable working condition wind power planetary gear box fault diagnosis method based on divergence index | |
Han et al. | Fault feature extraction of low speed roller bearing based on Teager energy operator and CEEMD | |
CN205067090U (en) | Antifriction bearing fault detection and diagnostic system | |
CN104596766B (en) | Early fault determining method and device for bearing | |
Wang et al. | Multi-scale enveloping order spectrogram for rotating machine health diagnosis | |
Wang et al. | A time–frequency-based maximum correlated kurtosis deconvolution approach for detecting bearing faults under variable speed conditions | |
CN105928702B (en) | Variable working condition box bearing method for diagnosing faults based on form PCA | |
CN104359674A (en) | High-speed rolling bearing fault diagnosing method based on time domain and frequency domain state monitoring | |
CN109029999B (en) | Rolling bearing fault diagnosis method based on enhanced modulation bispectrum analysis | |
Guo et al. | Enhancing the ability of ensemble empirical mode decomposition in machine fault diagnosis | |
CN104374575A (en) | Wind turbine main bearing fault diagnosis method based on blind source separation | |
Wang et al. | Weak fault diagnosis of rolling bearing under variable speed condition using IEWT-based enhanced envelope order spectrum | |
Lin et al. | A review and strategy for the diagnosis of speed-varying machinery | |
CN109655266A (en) | A kind of Wind turbines Method for Bearing Fault Diagnosis based on AVMD and spectral coherence analysis | |
CN105547627B (en) | Rotating machinery feature extracting method based on WPT-CEEMD | |
CN111582248B (en) | SVD-based gearbox signal noise reduction method | |
CN109900447A (en) | A kind of circulation impact method for detecting vibration and system based on harmonic signal decomposition | |
Makowski et al. | Parametric time-frequency map and its processing for local damage detection in rotating machinery | |
Mauricio et al. | Cyclostationary-based tools for bearing diagnostics | |
CN112345248B (en) | Fault diagnosis method and device for rolling bearing | |
Pang et al. | A Fault Feature Extraction Method for a Gearbox with a Composite Gear Train Based on EEMD and Translation-Invariant Multiwavelet Neighbouring Coefficients. | |
Chen et al. | Time-frequency demodulation analysis for gearbox fault diagnosis under nonstationary conditions | |
Zhang et al. | Fault diagnosis technology of rolling bearing based on LMD and BP neural network | |
Lu et al. | Bearing fault diagnosis based on redundant second generation wavelet denoising and EEMD |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140702 |