materials Article
Bearing Fault Detection Based on Empirical Wavelet Transform and Correlated Kurtosis by Acoustic Emission Zheyu Gao 1 , Jing Lin 1,2, *, Xiufeng Wang 1 and Xiaoqiang Xu 1 1
2
*
Shanxi Key Laboratory of Mechanical Product Quality Assurance and Diagnostics, Xi’an Jiaotong University, Xi’an 710049, China;
[email protected] (Z.G.);
[email protected] (X.W.);
[email protected] (X.X.) State Key Laboratory of Manufacturing System Engineering, Xi’an Jiaotong University, Xi’an 710049, China Correspondence:
[email protected]
Academic Editor: Victor Giurgiutiu Received: 15 April 2017; Accepted: 17 May 2017; Published: 24 May 2017
Abstract: Rolling bearings are widely used in rotating equipment. Detection of bearing faults is of great importance to guarantee safe operation of mechanical systems. Acoustic emission (AE), as one of the bearing monitoring technologies, is sensitive to weak signals and performs well in detecting incipient faults. Therefore, AE is widely used in monitoring the operating status of rolling bearing. This paper utilizes Empirical Wavelet Transform (EWT) to decompose AE signals into mono-components adaptively followed by calculation of the correlated kurtosis (CK) at certain time intervals of these components. By comparing these CK values, the resonant frequency of the rolling bearing can be determined. Then the fault characteristic frequencies are found by spectrum envelope. Both simulation signal and rolling bearing AE signals are used to verify the effectiveness of the proposed method. The results show that the new method performs well in identifying bearing fault frequency under strong background noise. Keywords: acoustic emission; correlated kurtosis; Empirical Wavelet Transform; bearing fault detection
1. Introduction Rolling bearings are widely used in industry equipment and are key components of rotating machinery. Bearings are inevitably damaged because of overload and long-time operation. Fault bearings make vibration of equipment exceeding operation standard. Badly damaged bearings even lock the whole system and destroy the rotating axle, which initiates downtime and financial loss [1]. Therefore monitoring bearing fault is of great importance to ensure a safe operation of a machine. Various kinds of monitoring technologies are employed to detect bearing faults, such as encoder, oil, ferro-graphic, vibration and acoustic emission (AE) [2]. Among them, the vibration signal is the most widely used. However, since AE signal can detect early faults of machinery components, the bearing fault monitoring based on AE has been of interest to researchers in recent years [3]. AE is transient stress wave rapidly released by localized energy [4]. During the operation of bearing, main sources of AE include impacting and friction of bearing rollers passing over outer and inner race and fatigue damage of bearing components [5]. Both traditional AE characteristic parameters and newer signal processing methods are employed to analyze AE signals. Traditional AE characteristic parameters are studied thoroughly to detect rolling bearing faults [6–8]. The research results indicate that AE technology has better performance over vibration analysis for monitoring early deterioration of rolling bearing [3]. Morhain and Mba studied the relationship of RMS, energy and count values with rotational speed, load and the size of the faults [8].
Materials 2017, 10, 571; doi:10.3390/ma10060571
www.mdpi.com/journal/materials
Materials 2017, 10, 571
2 of 11
Abdullah first used duration time of AE pulses to indicate the defect size of the outer race [9]. Since AE signal is always disturbed by noise, Khamis first proposed the Energy Index (EI) technique to enhance the AE burst buried in random noise. The EI index is utilized to detect incipient faults of bearings [5]. A variety of signal processing methods based on AE technology are used to monitoring bearing condition. Kurtogram and spectral kurtosis are applied to enhance AE events [10]. Kilundu first investigated cyclo-stationarity characteristic of AE bearing signals and used it to find outer race defects [11]. Pandya utilized Hilbert-Huang Transform (HHT) to obtain the time-frequency distribution of AE bearing signals. Then the K-nearest neighbor algorithm based on asymmetric proximity function was employed to classify the signal features, by which the reliability of bearing fault diagnosis increased [12]. Zvokelj proposed a new method, Ensemble Empirical Mode Decomposition-based multiscale Independent Component Analysis, to monitoring bearing fault. The algorithm can not only detect bearing faults but also de-noise AE signals [13]. Farzad Hemmati further developed the bearing fault detection technology by optimizing the ratio of Kurtosis and calculating Shannon entropy [14]. Mourad Kedadouche first used EWT to AE bearing signals to find frequency band which is excited more than other bands. The proposed method successfully detected defects as low as 40 microns on bearing outer race [15]. Although a lot of investigations are devoted to detect bearing faults, there are still two problems to be solved. Firstly, most of the methods are only effective in outer race defect, but ineffective in inner race fault due to signal attenuation. Secondly, AE signals are easily disturbed by background noise, so method with good anti-noise performance is needed. Therefore, in this paper a new method based on Empirical Wavelet Transform (EWT) and correlated kurtosis (CK) is proposed. Analysis on simulated signal proves that the resonant frequency band and fault characteristic frequency can be effectively obtained with this method. The rolling bearing experiment is carried out and AE signals of defect bearing are used to prove the efficiency of the proposed method. The rest of the paper is organized as follows. The EWT algorithm and CK are introduced in Section 2. In Section 3, the proposed method is applied to simulation signal to validate its effectiveness. Then, in Section 4, rolling bearing experiment is carried out and the method is employed to detect the inner race fault of bearing. Section 5 summarizes the conclusion. 2. Introduction to the Theory 2.1. Empirical Wavelet Transform EWT divides signals into frequency sub-bands adaptively according to the information located in the spectrum of the original signal [16,17]. The signal can be decomposed into sub-bands by a set of wavelet construction. The detail of the EWT algorithm is as follows. Suppose that the signal has a Fourier support within [0, π ]. The spectrum can be divided into Nth contiguous partitions. Each component has the frequency bandwidth as: Λn = [ωn−1 , ωn ], where ω0 = 0, ω N = π, ∪nN=1 Λn = [0, π]. The empirical scaling function and the empirical wavelets are shown in Equations (1) and (2) respectively. 1 h if |ω | ≤ ωn − τn i ∧ π 1 Φn (ω ) = cos 2 β 2τn (|ω | − ωn + τn ) if ωn − τn ≤ |ω | ≤ ωn + τn , 0 otherwise 1 h i if ωn + τn ≤ |ω | ≤ ωn−1 − τn+1 π 1 cos 2 β 2τ (|ω | − ωn+1 + τn+1 ) if ωn+1 − τn+1 ≤ |ω | ≤ ωn+1 + τn+1 ∧ i h n −1 Ψn (ω ) = 1 π sin 2 β 2τn (|ω | − ωn + τn ) if ωn − τn ≤ |ω | ≤ ωn + τn 0 otherwise
(1)
,
(2)
Materials 2017, 10, 571
3 of 11
Here, Tn is the transition phase corresponding to the ωn , which has the width of 2τn . The β(χ) is an arbitrary function. In this paper, it is defined as follows: β(χ) = χ4 35 − 84χ + 70χ2 − 20χ3 ,
(3)
Similar to the continuous wavelet transformation (CWT), the detail coefficients and approximation coefficients can be obtained by the empirical wavelets and scaling function, the functions are shown as follows: ! Wχ (n, t) =
Wχ (0, t) =
Z
Z
∧
∧
χ(τ )ψn (τ − t)dτ = F −1 X (ω )ψn (ω ) , ∧
∧
(4)
!
χ(τ )φ1 (τ − t)dτ = F −1 X (ω )φ1 (ω ) ,
(5)
The extracted modes can be obtained as below: f 0 (t) = Wχ (0, t) × φ1 (t),
(6)
f k (t) = Wχ (k, t) × ψk (t),
(7)
In order to ensure the decomposition has a tight frame, Equation 8 must be satisfied. γ < min( where γ =
τn ωn .
ω n +1 − ω n ), ω n +1 + ω n
(8)
More detailed information can be found in reference [16].
2.2. Correlated Kurtosis According to reference [18], the definition of CK is shown as follows: CK M ( T ) =
M 2 ∑nN=1 (∏m =0 yn−mT )
(∑nN=1 y2n )
M +1
,
(9)
Here, M behalves the order of shift. N is the length of the input signal. T is the sampling point which can be calculated as following: fs T= , (10) fc where f s is the sampling frequency; f c is the fault character frequency. In practice, the T varies in a certain range due to the fluctuation of rotational speed. From the definition, it is known that CK reveals the periodicity of signal. If T is the period of fault, the CK can indicate the occurrence of fault by its magnitude. The frequency band with maximum value of CK is the frequency band which is excited most by the fault characteristic frequency corresponding to T. For rolling bearing, the fault frequencies are always modulated to the resonant frequency and the higher order resonant frequencies. Since the fault frequencies of rolling bearing have obvious periodicity, the CK at time interval T that equals to the fault period will achieve higher value at resonant frequency bands. Therefore, the frequency bands with the highest CK value can be considered as the resonant frequency or the higher order resonant frequencies. According to the analysis above, a new method is proposed to find resonant frequency and fault frequency by EWT and CK. The whole procedure of the proposed algorithm is shown as a flowchart (Figure 1).
periodicity, the CK at time interval T that equals to the fault period will achieve higher value at resonant frequency bands. Therefore, the frequency bands with the highest CK value can be considered as the resonant frequency or the higher order resonant frequencies. According to the analysis above, a new method is proposed to find resonant frequency and fault frequency by EWT Materials 2017, 10, 571and CK. The whole procedure of the proposed algorithm is shown as a flowchart 4 of 11 (Figure 1).
Figure Figure1.1.Flowchart Flowchart of of algorithm. algorithm.
3. Simulation Verification 3. Simulation Verification To To verify thethe proposed is employed. employed.When Whena alocal localfault fault a rolling verify proposedmethod, method,aasimulation simulation signal signal is inin a rolling bearing emerges, harmonic components and periodical impulses will be generated. These frequency bearing emerges, harmonic components and periodical impulses will be generated. These frequency components areare modulated of bearing bearingand anddecay decayexponentially exponentially [19]. components modulatedby byresonance resonancefrequency frequency of [19]. TheThe simulation signal of defect bearing is expressed as follows: simulation signal of defect bearing is expressed as follows: MM
t ) =∑ t −kT kT) )++nn((tt)),, x (tx) (= A0As0(st(−
(11) (11)
k=k1=1
s(t) = e−ζt cos(2π f t + ϕ ), (12) s ( t ) = e −ζ t cos(2πnf n t + ϕ00 ) , (12) where A0 is the impulse response amplitude; ζ is the damping ratio of the impulse; M is the number of A0 T is the period of the impulses which isζ related to the fault characteristic frequency fM impulses; n isthe is the impulse response amplitude; is the damping ratio of the impulse; c ; f is where the resonance frequency; ϕ is the primary phase; n ( t ) is white noise. The value of these parameters 0 number of impulses; T is the period of the impulses which is related to the fault characteristic are displayed as following:
f
f
ϕ
n(t)
is white noise. The frequency c ; n is the resonance frequency; 0 is the primary phase; ζ = 160; Mare = 3; T = 0.05 as s; following: f c = 20 Hz; f n = 1300 Hz; ϕ0 = 0; SNR = −20 dB. 0 = 1;parameters value ofAthese displayed First, is utilized simulated adaptively. Then=the A0 =the 1; EWT ζ = 160; M = 3;toTdecompose = 0.05s; fcthe = 20 Hz; fn =signal 1300Hz ; ϕ0 = 0; SNR −20CK dB.of these modes are calculated at the given value of T. The waveforms of the simulated signal with and without noise are shown in Figure 2. Figure 3 displays the spectrum and envelope spectrum of the simulation signal with noise. It shows that the frequency 20 Hz is buried in strong background noise. The spectrum segmentation by EWT is shown in Figure 4. It is hard to find resonant frequency only by the amplitude of the spectrum. Figure 5 is the CK value of each component at time interval T. The peak frequency appears at 1275 Hz. It approximates to the true value of resonant frequency, which is 1300 Hz. The error between them can be ascribed to the strong noise. Figure 6 is the envelope spectrum of frequency band
modes are calculated at the given value of T. The waveforms of the simulated signal with and without noise are shown in Figure 2. Figure 3 displays the spectrum and envelope spectrum of the simulation signal with noise. shows that the frequency Hz is buried in strong background noise. The noise are shown in It Figure 2. Figure 3 displays the20 spectrum and envelope spectrum of the simulation signal with noise. It shows that the frequency 20 Hz is buried in strong background noise. The spectrum segmentation by EWT is shown in Figure 4. It is hard to find resonant frequency onlyThe by signal with noise. It shows that the frequency 20 Hz is buried in strong background noise. spectrum segmentation by EWT is shown in Figure 4. It is hard to find resonant frequency only by the amplitude of the spectrum. 5 is in theFigure CK value of hard each to component at time intervalonly T. The spectrum segmentation by EWTFigure is shown 4. It is find resonant frequency by the amplitude of the spectrum. Figure 5 is the CK value of each component at time interval T. The peak frequency of appears at 1275 Hz. It approximates to theof true value of resonant frequency, is the amplitude the spectrum. Figure 5 is the CK value each component at time intervalwhich T. The Materials 2017, 10,appears 571 5 of 11 is peak frequency at 1275 Hz. It approximates to the true value of resonant frequency, which 1300 The error between them can be ascribedtotothe the strong Figurefrequency, 6 is the envelope peakHz. frequency appears at 1275 Hz. It approximates true valuenoise. of resonant which is 1300 Hz. The error between them can be ascribed to the strong noise. Figure 6 is the envelope spectrum frequency band centered 1275 Hz. In strong Figure noise. 6, the Figure frequency and its 1300 Hz. of The error between them canaround be ascribed to the 6 is 20 theHz envelope spectrum of frequency band centered around 1275 Hz. In Figure 6, the frequency 20 Hz and its centeredof around 1275 Hz. Figure 6, the Hzwell andFigure its identifying harmonics are obvious.frequency is clearly harmonics are obvious. Itband isIn clearly that thefrequency CK performs in resonant spectrum frequency centered around 127520Hz. In 6, the frequency 20It Hz andand its harmonics are obvious. It is clearly that the CK performs well in identifying resonant frequency and that the CK performs well in identifying resonant frequency and periodic frequency under low SNR. periodic frequency underIt low SNR. that the CK performs well in identifying resonant frequency and harmonics are obvious. is clearly periodic frequency under low SNR. periodic frequency under low SNR.
Amp [V] Amp [V] Amp [V]
Amp [V] Amp [V] Amp [V]
2 2 2 0 0 0 -2 -20 -20 0 10 10 10 0 0 0 -10 -100 -100 0
0.2 0.2 0.2
0.4 0.6 0.4Time [s]0.6 0.4 Time [s]0.6 Time [s]
0.8 0.8 0.8
1 1 1
Amp [V] Amp [V] Amp [V] Amp [V] Amp [V] Amp [V]
0.4 0.6 0.8 1 0.4 0.6 0.8 1 Time [s]0.6 0.4 0.8 1 Time [s] Time [s] Figure 2. 2. The signal and simulated simulated signalwith with noise. Figure Thesimulated simulatedimpulse impulse signal signal and Figure 2. The simulated impulse simulatedsignal signal withnoise. noise. Figure 2. The simulated impulse signal and simulated signal with noise. 0.2 0.2 0.2
0.1 0.1 0.1
0.1 0.1 0.1
00 0 00 0 0.2 0.2 0.2
00 0 00 0
0.2 0.2 0.2
1000 1000 1000
2000 3000 2000 3000 Freq [Hz] 2000 3000 Freq [Hz] Freq [Hz]
4000 4000 4000
5000 5000 5000
X: X: 20 20 Y: X: 0.07001 20 Y: 0.07001 Y: 0.07001
100 150 200 100 150 200 Freq [Hz] 100 150 200 Freq [Hz] Freq [Hz] Figure spectrum of of simulation simulation signal signal with with noise. noise. Figure 3. 3. The The spectrum spectrum and and envelope envelope spectrum Figure3.3.The Thespectrum spectrumand andenvelope envelope spectrum spectrum of noise. Figure ofsimulation simulationsignal signalwith with noise. 0.12 0.12 0.12 0.1 0.1 0.1
Amp [V] Amp [V] Amp [V]
50 50 50
0.08 0.08 0.08 0.06 0.06 0.06 0.04 0.04 0.02 0.02 0 00 0
1 1
2 3 44 55 2 [kHz]3 4 5 Freq Freq [kHz] Figure 4. The segmentation of spectrum of simulation simulation signal. signal. Figure ofsimulation simulationsignal. signal. Figure4.4.The Thesegmentation segmentation of of spectrum spectrum of
Materials 2017, 10, 571 Materials 2017, 10, 571
6 of 11 6 of 10
Materials 2017, 10, 571
6 of 10 -12
Correlated Kurtosis Correlated Kurtosis
x 10 5 -12 x 10 5 4 4
X: 1275 Y: 4.211e-12 X: 1275 Y: 4.211e-12
3 3 2 2 1 1 0 00 0
1000 1000
2000 3000 2000 3000 Freq [Hz] Freq [Hz]
4000 4000
5000 5000
Figure of components components withtime timeinterval intervalequals equals Figure5.5. 5.Correlated Correlatedkurtosis kurtosis(CK) (CK) value value of components with with time interval equals toto T.T. Figure Correlated kurtosis (CK) value to T.
0.2 0.2
Amp [V] Amp [V]
0.15 0.15 X: 20 20 Y: 0.1057 0.1057 X:40 40 X: Y:0.08086 0.08086 Y:
0.1 0.1
0.05 0.05
000 0
20 20
40 60 40 60 Freq Freq [Hz] [Hz]
80 80
100 100
Figure Envelope spectrum Figure6.6. 6.Envelope Envelopespectrum spectrum of of component component with with maximum CK value. Figure of component withmaximum maximumCK CKvalue. value.
Experimental Verification 4. ExperimentalVerification Verification 4.4. Experimental Our rolling bearing experiment is explained this section. The experimental setup is displayed Ourrolling rollingbearing bearingexperiment experiment is is explained explained in in Our in this this section. section.The Theexperimental experimentalsetup setupisisdisplayed displayed in Figure consisted of motor diver, shaft, bearings, bearing housing, inertia disc and control Figure7.7. 7.ItIt Itconsisted consistedof of motor motor diver, diver, shaft, ininFigure shaft, bearings, bearings, bearing bearinghousing, housing,inertia inertiadisc discand andcontrol control system. The shaft was supported by two bearings. The left one was normal and the right one system. The shaft was supported by two bearings. The left one was normal and the right one defective. system. The shaft was supported by two bearings. The left one was normal and the right one defective. An inner race fault rolling bearings was used in this experiment. defect was caused An innerAn race faultrace rolling bearings used in this experiment. The defectThe was defective. inner fault rollingwas bearings was used in this experiment. Thecaused defect artificially. was caused artificially. The rotation speed of the experiment was 3000 rpm. The rotation of the experiment was 3000 rpm. artificially. Thespeed rotation speed of the experiment was 3000 rpm. The AE signals were picked up in the comparing performance of kinds of The AE signals were picked up in the experiment. experiment. By By comparing performance ofofthree three kinds ofof The AE signals were picked up in the experiment. By comparing performance three kinds AE sensors, Nano 30 was selected. The AE sensor was mounted on the top surface of bearing housing AE sensors, Nano 30 was selected. The AE sensor was mounted on the top surface of bearing housing AE sensors, The Nano 30 was selected. The is AE sensor was mounted on the topsignals surface of bearing housing (Figure (Figure 7). 7). The parameters parameters of of Nano Nano 30 30 is shown shown in in Table Table 1. 1. The The original original AE AE signals were were amplified amplified by by (Figure 7). The parameters of Nano 30 is shown in Table 1. The Acoustics original AE signals were amplified by 40 dB. The AE signal collection system was PCI-II (Physical Corporation, Boston, MA, 40 dB. The AE signal collection system was PCI-II (Physical Acoustics Corporation, Boston, MA, USA). 40USA). dB. The AE signal collection system was PCI-II (Physical Acoustics Corporation, Boston, MA, The sample rate was 500 kHz AE threshold level setdB. to 40 dB. The sample rate was 500 kHz and theand AEthe threshold level was setwas to 40 USA). The sample rate was 500 kHz and the AE threshold level was set to 40 dB. Table Table 1. 1. Parameters Parameters of of Nano30. Nano30. Table 1. Parameters of Nano30.
Parameters Range Parameters Range Parameters Range Operating Frequency Range 125–750 KHz Operating Frequency Range 125–750 KHz Operating Frequency Range 125–750 KHz Resonant Frequency 300 KHz Resonant Frequency 300 KHz Temperature Range −65–177 °C Resonant Frequency 300 KHz Temperature Range −65–177 ◦ C Temperature Range −65–177 °C When defects occur in the rolling bearings, some certain frequency components will be generated to thein rotational speed, geometricsome parameters bearings and the type of faults When according defects occur the rolling bearings, certainoffrequency components will be [20]. The formulae frequencies are as follows: generated accordingoftothe thedefective rotational speed, geometric parameters of bearings and the type of faults [20]. The formulae of the defective frequencies are as follows:
Materials 2017, 10, 571
7 of 11
When defects occur in the rolling bearings, some certain frequency components will be generated according to the rotational speed, geometric parameters of bearings and the type of faults [20]. Materials 2017, 10, 571 7 of 10 The formulae of the defective frequencies are as follows: Ball-pass frequency of outer race: Ball-pass frequency of outer race: n f rnf dd BPFO = 1 − cos φ (13) BPFO =2 r {1−D cosφ} , (13) D 2 , Ball-pass frequency of inner race: Ball-pass frequency of inner race: n fr d nfr1 + d cos φ , BPFI = (14) (14) BPFI =2 {1+D cosφ} ,
2
D
Fundamental train frequency (cage speed): Fundamental train frequency (cage speed): d fr f 1 − d cos φ , FTF = FTF =2 r {1−D cosφ} 2 D , Ball (roller) spin frequency: Ball (roller) spin frequency: ( ) 2 d D f rDf 1 − ( dcos φ) 2 BSF = BSF 2d = r {1 − D ( cos φ ) }
2d
(15) (15)
(16) (16)
D
Here,fr fis theshaft shaftrotation rotationfrequency, frequency, nn is is the contact r is the the number number of ofrolling rollingelements, elements,and anda аisisthe the contact Here, angle of the load. D and d represent the diameter and pitch diameter of rolling element respectively. angle of the load. D and d represent the diameter and pitch diameter of rolling element respectively. Thefault faultfeature featurefrequencies frequencies of of the the rolling rolling bearing The bearingare arelisted listedininTable Table2.2.
Figure 7. Experimental set up. Figure 7. Experimental set up. Table 2. Fault character frequency of rolling bearing. Table 2. Fault character frequency of rolling bearing.
Rotation Speed/rpm Rotation 3000 Speed/rpm 3000
BPFO/Hz BPFO/Hz 152.4 152.4
BPFI/Hz BPFI/Hz 247.5 247.5
BSF/Hz BSF/Hz99.6 99.6
Figure 8 shows waveform and spectrum of rolling bearing with inner race faults. Figure 9 is the envelope spectrum ofwaveform AE bearing signal. It canofbe seen that the with rotating frequency and its harmonics Figure 8 shows and spectrum rolling bearing inner race faults. Figure 9 is the are obvious,spectrum but the BPFI defined in Equation 14 and provided in Tableand 2 cannot be found envelope of AEwhich bearingissignal. It can be seen that the rotating frequency its harmonics inare theobvious, figure. Then the EWT is carried out in to Equation decompose AEinsignal. segmentation but the BPFI which is defined (14)the andbearing provided Table 2The cannot be found inof the is shown in Figure 10.out Figure 11 displays value components at time interval thespectrum figure. Then the EWT is carried to decompose the CK bearing AEofsignal. The segmentation of the T spectrum is shown in Figure 10. Figure 11 displays CKvalue value of components timeItinterval which which equals to inner fault period. The maximum appears at 66katHz. meansTthat this equals toband innerisfault period. maximum value appears at 66kFigure Hz. It 12 means this frequency frequency excited mostThe by the inner race fault frequency. is thethat spectrum envelope of this frequency band. It is clear that the rotating frequency is also modulated to this band. However, in Figure 12, the amplitude of the inner race fault frequency increases and the higher order of inner race fault frequencies are clear. The result shows that the proposed method can improve the inner race fault detection of rolling bearing. It must be pointed out that the time interval T is the reciprocal of fault character frequency and
Materials 2017, 10, 571
8 of 11
band is excited most by the inner race fault frequency. Figure 12 is the spectrum envelope of this frequency band. It is clear that the rotating frequency is also modulated to this band. However, in Figure 12, the amplitude of the inner race fault frequency increases and the higher order of inner race fault frequencies are clear. The result shows that the proposed method can improve the inner race fault detection of rolling bearing. It must be pointed out that the time interval T is the reciprocal of fault character frequency and the fault frequency of rolling bearing is proportional to rotating speed, so the speed variation Materials 2017,character 10, 571 8 of 10 will affect the time interval T. If the speed variation is small, the accurate time interval T changes within Materials 2017, 10, 571 8 of 10 certainwithin range certain centeredrange around the reciprocal the reciprocal fault character frequency. If the speed variation changes centered aroundofthe of the fault character frequency. If is the large, the proposed method will not work. speed variation is large, proposed method work. of the fault character frequency. If the changes within certainthe range centered aroundwill thenot reciprocal speed variation is large, the proposed method will not work.
Amp [V] Amp [V]
10
10
0 0
-10 0 -10 0
0.02 0.02
0.04 0.06 Time [s]0.06 0.04 Time [s]
0.08 0.08
0.5 0.5
1 1.5 1 Freq [Hz]1.5 Freq [Hz]
2 2
Amp [V] Amp [V]
1
0.1 0.1
1
0.5 0.5 0 00 0
2.5 2.5 5 x 105 x 10
Figure 8. 8. The waveform bearing signal. Figure The waveformand andspectrum spectrum of of AE AE bearing bearing signal. signal. Figure 8. The waveform and spectrum of AE
0.04 0.04
Amp [V]
Amp [V]
0.03 0.03
0.02 0.02
X:X:5050 Y:Y:0.02525 0.02525 X:X:250 250 Y:Y:0.01033 0.01033
0.01 0.01
0 00 0
200 200
400 600 400 600 Freq[kHz] [kHz] Freq
800 800
1000 1000
Figure The envelope spectrum of AE Figure 9. Theenvelope envelopespectrum spectrumof of AE AE bearing bearing signal. signal. Figure 9.9. The signal.
0.7 0.7 0.6 0.6
Amp [V] Amp [V]
0.5 0.5 0.4 0.4 0.3 0.3 0.2 0.2
0 0
200
Materials 2017, 10, 571
400 600 Freq [kHz]
800
1000 9 of 11
Figure 9. The envelope spectrum of AE bearing signal.
0.7 0.6
Amp [V]
0.5 0.4 0.3 0.2 0.1 0 0
Materials Materials2017, 2017,10, 10,571 571
50
100 150 Freq [kHz]
200
250
Figure 10. The The segmentation segmentation of of spectrum spectrumof ofbearing bearingsignal. signal.
99ofof10 10
-11 -11
xx10 10 1.2 1.2
X: X:66.08 66.08 Y: Y:1.019e-11 1.019e-11
Correlated Kurtosis Kurtosis Correlated
11 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 00 00
50 50
100 150 100 150 Freq Freq[kHz] [kHz]
200 200
250 250
Figure Figure11. 11.CK CKvalue valueof ofcomponents componentsat attime timeinterval intervalT. T.
11
Amp [V] [V] Amp
0.8 0.8 0.6 0.6 0.4 0.4
X: X:247 247 Y: Y:0.2161 0.2161
0.2 0.2 00 00
200 200
400 600 400 600 Freq Freq[Hz] [Hz]
800 800
1000 1000
Figure Figure12. 12.Envelope Envelopespectrum spectrumof ofcomponent componentwith withmaximum maximumCK CKvalue. value. maximum CK value.
5.5.Conclusions Conclusions In Inthis thispaper, paper,aanew newmethod methodisisproposed proposedto todetect detectdefects defectsin inrolling rollingbearing. bearing.EWT EWTisisutilized utilizedto to decompose AE signals into different components. CK of these components at a certain time interval decompose AE signals into different components. CK of these components at a certain time interval TTare arecalculated. calculated.Components Componentswith withthe thelargest largestCK CKvalue valueisisconsidered consideredto tobe bethe theresonant resonantfrequency. frequency.
Materials 2017, 10, 571
10 of 11
5. Conclusions In this paper, a new method is proposed to detect defects in rolling bearing. EWT is utilized to decompose AE signals into different components. CK of these components at a certain time interval T are calculated. Components with the largest CK value is considered to be the resonant frequency. Then envelope demodulation is carried out to obtain the fault feature characteristics of corresponding frequency bands. Simulated and experimental signals were employed to verify the effectiveness of the method. Two conclusions can be drawn as follows: 1. 2.
With the right time interval T, the resonant frequency of the rolling bearing system can be obtained by calculating the CK value. Weak fault features can be identified by the proposed method under low SNR conditions.
Future research work will be devoted to enhance the fault frequencies modulated by different resonant frequencies and estimation of the accurate time interval T under slight speed fluctuation will also be carried out. Acknowledgments: This work was funded by the National Natural Science Foundation of China [Grant No. 51421004; 51505365]; and the Fundamental Research Funds of the Central Universities of China [Grant No. CXTD2014001]. Author Contributions: Zheyu Gao developed the methodology and wrote the manuscript. Jing Lin provided the professional guidance to complete the research work and manuscript revision. Xiufeng Wang provided great assistance in doing experiments and provided professional guidance. Xiaoqiang Xu provided great advice. Conflicts of Interest: The authors declare no conflict of interest.
References 1. 2. 3.
4. 5. 6. 7. 8. 9.
10. 11. 12. 13.
Hecke, B.E.V. Development of Novel Acoustic Emission Based Methodology and Tools for Bearing Fault Diagnostics. Ph.D. Thesis, University of Illinois, Chicago, IL, USA, 2015. Cao, H.R.; Fan, F.; Zhou, K.; He, Z.J. Wheel-bearing fault diagnosis of trains using empirical wavelet transform. Measurement 2016, 82, 439–449. [CrossRef] Mba, D.; Rao, R. Development of Acoustic Emission technology for the condition monitoring and diagnosis of rotating machines: bearings, pumps, gearboxes, engines and rotating structures. Shock Vib. Dig. 2006, 38, 3–16. [CrossRef] Tandon, N.; Choudhury, A. A review of vibration and acoustic measurement methods for the detection of defects in rolling element bearings. Tribol. Int. 1999, 32, 469–480. [CrossRef] Al-Balushi, K.R.; Addali, A.; Charnley, B.; Mba, D. Energy Index technique for detection of Acoustic Emissions associated with incipient bearing failures. Appl. Acoust. 2010, 71, 812–821. [CrossRef] Balerston, H.L. The detection of incipient failure in bearings. Mater. Eval. 1969, 27, 121–128. Yoshioka, T.; Korenaga, A.; Mano, H.; Yamamoto, T. Diagnosis of rolling bearing by measuring time interval of Acoustic Emission generation. J. Tribol. 1999, 121, 468–472. [CrossRef] Morhain, A.; Mba, D. Bearing Defect Diagnosis and Acoustic Emission. Proc. Inst. Mech. Eng. Part J: J. Eng. Tribol. 2003, 217, 272–275. [CrossRef] Al-Ghamd, A.M.; Mba, D. A comparative experimental study on the use of acoustic emission and vibration analysis for bearing defect identification and estimation of defect size. Mech. Syst. Signal Process. 2006, 20, 1537–1571. [CrossRef] Eftekharnejad, B.; Carrasco, M.R.; Charnley, B.; Mba, D. The application of spectral kurtosis on Acoustic Emission and vibrations from a defective bearing. Mech. Syst. Signal Process. 2011, 25, 266–284. [CrossRef] Kilundu, B.; Chiementin, X.; Duez, J.; Mba, D. Cyclostationarity of Acoustic Emissions (AE) for monitoring bearing defects. Mech. Syst. Signal Process. 2011, 25, 2061–2072. [CrossRef] Pandya, D.H.; Upadhyay, S.H.; Harsha, S.P. Fault diagnosis of rolling element bearing with intrinsic mode function of acoustic emission data using APF-KNN. Expert Syst. Appl. 2013, 40, 4137–4145. [CrossRef] Žvokelj, M.; Zupan, S.; Prebil, I. EEMD-based multiscale ICA method for slewing bearing fault detection and diagnosis. J. Sound Vib. 2016, 370, 394–423. [CrossRef]
Materials 2017, 10, 571
14. 15.
16. 17.
18. 19. 20.
11 of 11
Hemmati, F.; Orfali, W.; Gadala, M.S. Roller bearing acoustic signature extraction by wavelet packet transform, applications in fault detection and size estimation. Appl. Acoust. 2016, 104, 101–118. [CrossRef] Kedadouche, M.; Thomas, M.; Tahan, A. Empirical Mode Decomposition combined with Empirical Wavelets for extracting bearing frequencies in a noisy environment and early detection of defects. In Advances in Condition Monitoring of Machinery in Non-Stationary Operations; Springer: Lyon, France, 2014; pp. 151–165. Gilles, J. Empirical Wavelet Transform. IEEE Trans. Signal Process. 2013, 61, 3999–4002. Pan, J.; Chen, J.L.; Zi, Y.Y.; Li, Y.M.; He, Z.J. Mono-component feature extraction for mechanical fault diagnosis using modified empirical wavelet transform via data-driven adaptive Fourier spectrum segment. Mech. Syst. Signal Process. 2016, 72–73, 160–183. [CrossRef] McDonald, G.L.; Zhao, Q.; Zuo, M.J. Maximum correlated Kurtosis deconvolution and application on gear tooth chip fault detection. Mech. Syst. Signal Process. 2012, 33, 237–255. [CrossRef] McFadden, P.D.; Smith, J.D. Model for the vibration produced by a single point defect in a rolling element bearing. J. Sound Vib. 1984, 96, 69–82. [CrossRef] Randall, R.B.; Antoni, J. Rolling element bearing diagnostics—A tutorial. Mech. Syst. Signal Process. 2011, 25, 485–520. [CrossRef] © 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).