Adaptive Asymmetric Real Laplace Wavelet Filtering and Its Application on Rolling Bearing Early Fault Diagnosis

.e early fault of rolling bearing is weak andmay not be readily detected. To overcome this issue, the present paper comes up with a rolling bearing fault-diagnosing approach based on adaptive asymmetric real Laplace wavelet (ARLW) filtering, which is on the strength of water cycle optimization algorithm (WCA). Firstly, ARLW is introduced to filter the initial vibration signal since its waveform has the same asymmetric structure as the fault impact. Secondly, the optimum center frequency and bandwidth of ARLW is found out adaptively by applying the WCA through the proposed square envelope fault energy ratio (SEFER). Finally, envelope analysis is conducted to the narrowband signal obtained by the optimum ARLW filtering, and its envelope spectrum presents the rolling bearing fault characteristic frequency apparently. .e proposed approach and two existing approaches are all tested in four signal analysis cases..e results are analyzed, and the conclusion is that the approach proposed by the present paper can detect the early fault of rolling bearing more accurately. .e present research is valuable for diagnosing the early fault of rolling bearing.


Introduction
Rolling bearing is broadly used in rotating equipment, and its fault acts on the safe operation of the whole equipment [1][2][3].At the beginning of the rolling bearing fault, the impact component of vibration signal collected by the sensor is weak and often submerged in strong background noise, bringing challenges to the diagnostic process [4].Consequently, diagnosing the early fault of rolling bearing acts as the focus and difficulty among researchers and scholars.
e periodic impact produced by partial defects on the roiling bearing surface can arouse the resonance between the rolling bearing and its adjacent parts [5].Using resonance demodulation to extract impact response characteristics from vibration signal is a fast and simple approach of rolling bearing fault diagnosis, the pivotal step of which is to precisely identify the resonance frequency band that contains plentiful fault information.e traditional resonance demodulation approach has the limitation of the request of artificially presetting the parameters of the bandpass filter (center frequency and bandwidth).
Antoni et al. creatively put forward the spectral kurtosis theory [6,7] and the fast spectrum kurtosis (FSK) approach [8] that can automatically set the parameters of the bandpass filter, which firstly used short time Fourier transform (STFT) or finite impulse response (FIR) filters to divide frequency bands and then used the kurtosis of signal as an evaluating indicator to identify the frequency band containing the most fault information.However, there are two disadvantages in the FSK approach.
e first one is that kurtosis as an evaluation index cannot distinguish the random impact and periodic impact of the signal, which easily leads to identifying the erroneous resonance band.e second one is that the optimal frequency band selected by FSK may not include the whole resonant frequency band region.For the first disadvantage, several new evaluating indicators are raised to improve the accuracy of selecting the resonant frequency band, such as correlated kurtosis [9], harmonic-to-noise ratio [10], Gini index [11], and spectral L2/L1 [12].
In addition, the noise reduction approach can be used as a preprocessing method of FSK to avoid interference, such as the improved FSK approach with the aid of EEMD [13], ICA [14], and ITD [15].Although the above approach improves the robustness of FSK, more prior knowledge is needed, and the calculation process is complex.For the second shortcoming, many scholars have done useful work.Lei et al. used the wavelet packet transform (WPT) instead of the STFT or FIR filters in the FSK, which segments the frequency band more finely [16][17][18].Comparing with WPT, continuous wavelet transform (CWT) can flexibly and accurately divide the frequency domain.CWT segments the signal through constructing a series of filters with the same property but different center frequency (CF) and bandwidth (BW) by translating and stretching the mother wavelet.
us, designing wavelet parameters properly to obtain the resonance frequency band is key to diagnosing rolling bearing fault by CWT.Qiu et al. [19] combines the Shannon entropy and SVD theory to realize the optimal wavelet transform.Bozchalooi et al. [20] uses the smoothing index to select the CF and BW of Gabor wavelet.e variable step size of the CF is 50 Hz, and the search range is [200,6000].e variable step size of the BW is 0.01, and the search range is [0.01 , 1].
at is to say, we need to calculate 11700 (117 × 110) times to determine the best solution.It is obvious that intelligent algorithms should be used to replace the huge computation process and to reduce the computation time and memory footprint.Su et al. [21] puts forward an approach based on the minimal entropy criterion, using genetic algorithm to obtain the optimum wavelet parameters.Wang et al. [22] obtains the parameters of Morlet wavelet using simulated annealing algorithm that adopted the maximum sparsity as the fitness function.Chen et al. [23] used a particle swarm optimization algorithm on the basis of the correlated kurtosis of squared envelope spectrum as fitness function to obtain the parameters of Morlet wavelet.
At present, there are three main issues in using wavelet filtering to determine the resonance frequency band of rolling bearing: the selection of the wavelet basic function, the selection of evaluation index, and fast determination of wavelet parameters.Accordingly, this paper proposed an adaptive asymmetric real Laplace wavelet (ARLW) filtering approach based on water cycle optimization algorithm (WCA), which has the following improvements: ARLW is used as wavelet basis function; square envelope fault energy ratio (SEFER) is chosen as a new evaluation index; WCA is used to choose the optimum wavelet parameters quickly.Four signal analysis cases are conducted to testify the validity of the proposed approach, the results of which are compared to that of the two existing approach.
e rest of the present paper is structured as follows: Section 2 reviews the theoretical background of wavelet filtering.Section 3 discusses the wavelet parameter optimization process.Section 4 proposes the rolling bearing fault diagnosis approach on the basis of adaptive wavelet filtering.Section 5 verifies the proposed approach through four signal analysis cases.Section 6 summarizes the full text.

Continuous Wavelet Transform.
e CWT of onedimensional signal x(t) is defined as where W(a, b) represents the wavelet coefficient, a represents the scale parameter, b represents the shift parameter, φ(•) represents the wavelet basis function, and * represents the conjugate.In the frequency domain, Equation ( 1) is expressed as where X(f) refers to the Fourier transform of x(t), ψ(f ) refers to the Fourier transform of φ(t), and IFT refers to the inverse Fourier transform.

Asymmetric Real Laplace Wavelet.
In the fault diagnosis of the mechanical system, more fault features can be extracted when the wavelet function matches the impact response in the dynamic signal.Among the existing wavelet basic functions, the time domain waveform (TW) of ARLW and Morlet wavelet is similar to that of impact response caused by local fault of rolling bearing.Figures 1(a)-1(c) depict the TW of ARLW, Morlet wavelet, and the actual rolling bearing signal, respectively.Compared with the Morlet wavelet with symmetrical structure, the waveform of ARLW with asymmetrical structure is more similar to fault waveform, so ARLW is more suitable for extracting rolling bearing fault features.e ARLW φ(t) can be expressed as the following equation [24,25]: Fourier transform of φ(t) is displayed as where c represents the center frequency (CF) and σ represents the bandwidth (BW).
e ARLW transform of one-dimensional vibration signal x(t) is expressed as where X(f) represents the Fourier transform of x(t), W x (c, σ) represents the result of ARLW filtering, and IFT represents the inverse Fourier transform.According to Equation (4), ARLW transform is the product of the analyzed signal and the ARLW in the frequency domain.Figure 1(d) illustrates that the frequency spectrum (FS) of ARLW has the window form.erefore, the ARLW transform has the same function as the bandpass filter.Based on the above analysis, designing wavelet parameters properly to obtain the resonance frequency band is the key to diagnosing rolling bearing fault by ARLW filtering.In addition, it should be 2 Shock and Vibration noted that c and σ are not completely independent in the ARLW filtering process.If the two parameters are optimized separately, the result is not the optimal solution.

Wavelet Parameter Adaptive Optimization
Water cycle optimization algorithm (WCA) proposed by Eskandar is a new embedded optimization approach [26,27], which can optimize multiple parameters in parallel.
Compared with genetic algorithm, simulated annealing algorithm, and other traditional optimization algorithms, WCA is more efficient and unlikely to fall into local optimal solution.Based on this, the optimal ARLW filter can be constructed adaptively with the help of WCA optimization algorithm.In the iterative process of WCA, the wavelet obtained by each iteration should be evaluated and determined whether it is the optimal solution.erefore, it is particularly necessary to select an appropriate evaluation index.

e Selection of Evaluation
Index.Aiming at improving the performance of envelope analysis in resonance demodulation, Randall [28] proposed a method, namely, square envelope analysis, which can suppress noise and highlight fault characteristics.Based on this advantage, the square envelope fault energy ratio (SEFER) is proposed, and the specific calculation process is as follows.
Step 1. Performing Hilbert transform on the signal x(t).
Step 2. Constructing the analytic signal z(t) and its conjugate complex analytic signal z ′ (t).
Step 3. e square envelope signal y(t) obtained by multiplying the analytic signal with the conjugate analytic signal.
Step 4. Fourier transform is conducted on the square envelope signal y(t).
where SES refers to squared envelope spectrum sequence and FT refers to the Fourier transform.
Step 5. e SEFER value of different fault type is calculated by the follow equation: where FE and SE denote the fault energy and signal total energy, respectively.SES (f) refers to the amplitude of each frequency in squared envelope spectrum.SES (nf ′ ) refers to the amplitude of the fault feature frequency and its frequency doubling in squared envelope spectrum.f ′ represents the feature frequency of different fault types.e periodic impact signal overlaying the Gaussian noise signal can be served to simulate the rolling bearing fault

Shock and Vibration
signal.e mathematical model of the periodic impact signal is defined as where the amplitude is A 1 � 0.5, the damping coefficient is g � 0.05, the resonance frequency is f n � 2500 Hz, the single impact sampling time is t 0 � 0.025 s, and the number of impact is K � 20. e sampling frequency and sampling points is 8096 Hz and 4096, respectively.Figure 2 displays the simulation signal with different SNRs, where the periodic impact characteristics (red waveform line) of signals become more prominent with the increase of SNR. e variation tendency of SEFER values for different SNR signals is presented in Figure 3, where the SEFER value increases monotonically with SNR increase.To sum up, we can say that the larger the SEFER value is, the more obvious impact characteristic the signal has.

Adaptive Optimization Process.
According to the analysis of Section 3.1, the larger the TFR value is, the more fault information the signal has.By comparing the SEFER value, we can judge the matching degree between ARLW and impulse components of the signal and determine the optimal ARLW filter according to the maximal value criterion.It should be noted that WCA is an algorithm for solving the minimization problem, and the process of finding the maximum SEFER should be converted into the minimization problem.e detailed steps of the adaptive optimization process of ARLW parameters based on WCA algorithm are as follows: Step 1. Setting the reciprocal of SEFER of the filtered signal obtained by ARLW filtering as the objective function of WCA as follows: where f ′ represents the feature frequency of different fault type.
Setting the constraint condition of the c and σ of ARLW.e literature shows that the optimal bandwidth of ARLW filter should be 3 times greater than the rolling bearing fault feature frequency.In addition, according to the sampling theory, the high-pass cutoff frequency of the ARLW filter needs to be less than half the time of the sampling frequency.
us, the constraint conditions of c and σ are as follows: where f ′ refers to the feature frequency of different fault type, f s refers to the sampling frequency, and σ min refers to the minimal value in the range of σ.
Step 2. Initializing the parameters of WCA algorithm.e number of design variables is N var � 2, the population size is N pop � 10, the number of rivers and sea is N sr � 2, the evaporation condition constant is d max � 1e − 5, and the maximal number of iterations is Max _it � 50.
Step 3. Generating the initial raindrops and determining the number of river and stream.e initial raindrops can be expressed as follows.
where x � LB + rand × (UB − LB), LB and UB represent the upper and lower limitations of the variables, and the rand represents a random number evenly distributed between 0 and 1.
Use Equation ( 15) to calculate the fitness function (objective function) of each raindrop: e best individual (the minimal fitness value) is chosen as the sea, and some better individuals are selected as the river.
e number of sea (optimum solution) is 1, the number of rivers is N r � N sr − 1 � 1, and the number of streams flowing to rivers or sea is determined by.
Step 4. In the process of water confluence, that is, raindrops form streams, some streams flow into the river and the others flow into the sea.e algorithm assumes that all rivers and streams become merged in the sea eventually.If the fitness value of a stream is greater than that of the adjacent river, exchange the location of the river and stream (i.e., stream becomes river and river becomes stream); if the fitness value of a river is greater than that of the adjacent sea, exchange the location of the sea and river (i.e., river becomes sea and sea becomes river).e location updating of stream and river can be expressed as follows: where the rand definition is the same as above, C refers to numbers between 1 and 2. X i stream , X i river , and X i river refer to the position of the stream, river, and sea at the i times' iteration, respectively.
e process of evaporation and precipitation.Evaporation condition can be expressed as follows: where d max can be reduced adaptively, and its calculation formula is as follows: 4

Shock and Vibration
When fulfilling the evaporation criterion, the new raindrops will form and then enter the process of water confluence again.
Step 6. Checking the terminating condition.When living up to the maximal iteration, stop the algorithm and output the optimal solution.Or else, go back to the Step 4 and continue to solve.

The Diagnosis Approach on the Basis of Adaptive Wavelet Filtering
Rolling bearing fault can be diagnosed quickly and effectively through the demodulation analysis of the resonance frequency band.Based on this idea, this paper came up with a rolling bearing fault diagnosis approach on the basis of adaptive ARLW filtering through WCA. Figure 4 describes the diagnosis process, and the detailed steps are as follows: Step 1. Loading the original signal and calculating the feature frequency of different fault type.e mathematical formulas are listed below [29].
where d and D represent the diameter of the balls and the pitch, respectively; α represents the contact angle between the ball and the raceway; Z represents the number of rolling element; and n represents the rotating speed.f i refers to the fault frequency of the inner ring (IRFF), f o refers to the fault frequency of the outer ring (ORFF), f e refers to the fault frequency of the rolling element (REFF), and f c refers to the fault frequency of the cage (CFF).
Step 2. Constructing the optimal bandpass filter according to the result of WCA optimization.e corresponding optimizing process can be seen in Section 3.2.
Step 3. Filtering the original signal with the constructed bandpass filter.
Step 4. Performing envelope demodulation on the filtered signal obtained by the optimal bandpass and comparing these frequencies corresponding to the spectral lines with larger amplitude in envelope spectrum with the fault characteristic frequency to judge the fault kind.

Validations for the Proposed Approach
In the following content, the proposed approach, the improved fast spectral kurtogram based on EEMD (EEMD-FSK) [13], and the improved kurtogram based on WPT (WPT-FSK) [16] are used to analyze a simulation signal, an artificial single fault signal, an artificial compound fault signal, and a life cycle fault signal, respectively.

Simulation Signal Case.
e rolling bearing fault model [30,31] is conducted for simulating the impact produced by the inner ring defect, and the strong white noise is overlaid for simulating the early fault.e mathematical model of simulation signal is as follows: where s(t) refers to the periodic impact component, n(t) refers to the Gaussian white noise, the amplitude refers to A 0 � 0.5, τ i � 0 refers to the slight fluctuation of the i-th time impact relative to periodic T, the attenuation coefficient refers to C � 1000, the resonance frequency refers to f n � 3 kHz, the shaft rotation frequency refers to f r � 20 Hz, the fault feature frequency refers to f i � 1/T �130 Hz, the random fluctuation obeys the zero-mean normal distribution, the standard deviation is 0.5% of the shaft rotation frequency, the SNR of the simulation signal is −12 dB, the sampling frequency refers to f s � 8192 Hz, and the sampling point number refers to N � 4096.Figure 5(a) presents the TW of the simulation signal, where the periodic impact components (red waveform lines) are masked by strong noise and thus difficult to observe fault information.Envelope demodulation is conducted on the simulation signal.Figure 5(b) illustrates the ES, where the IRFF f i is detected, but its amplitude is not prominent comparing with noise.e proposed approach is conducted to analyze the simulation signal.e optimization range of BW is set as σ ∈ [0, 390], and the optimization range of CF is set as c ∈ [180, 3916].WCA algorithm is initialized and 50 iterations are performed.Figure 6(a) shows the optimization curve of function value with evolution algebra.It can be seen that the minimum value is 29.75, obtained after the 8 times' iteration calculation.
e parameter combination consisting of CF and BW corresponding to the minimal value is chosen to construct the optimal ARLW filter.Figure 6(b) illustrates the FS of the optimal ARLW filter, where the CF � 3046 Hz and BW � 1349 Hz. e TW of the narrowband signal acquired from the optimal ARLW filter is displayed as Figure 6(c).Figures 6(d) and 6(e) depict the FS and ES of the filtered signal, respectively.Figure 6(d) shows that the CF of the optimal resonance frequency band obtained by adaptive optimization is basically consistent with the resonance frequency of the simulation signal, and the BW is 3 times significantly greater than fault feature frequency.As shown in Figure 6(e), the IRFF f i and its frequency doubling 2f i −3f i are detected productively.6 Shock and Vibration e EEMD-FSK approach is conducted to analyze the simulation signal.According to the diagnostic process in [13], the simulation signal is firstly decomposed into a group of IMFs by EEMD.e two IMFs with obvious impact characteristics are chosen by the measure factor to reconstruct signal.Figure 7  Shock and Vibration e WPT-FSK approach is conducted to analyze the simulation signal.According to the diagnostic process in [16], a number of distinct frequency band signals are obtained after the initial vibration signal is processed by WPT algorithm.e kurtosis value of different frequency band signals are calculated and then presented in the kurtgram, where the lateral and vertical axes represent the frequency and the decomposition level of WPT, respectively.Each node of kurtgram represents the narrowband signal obtained by WPT decomposition, and these narrowband signals are in the different frequency bands.
e signal corresponding to the node with the maximal kurtosis is chosen for envelope demodulation.Figure 8(a) presents that node (4, 13) has the maximal kurtosis.Figures 8(b) and 8(c) depict the TW and ES of the frequency band signal corresponding to node (4,13).As seen, although the impact feature is seen in TW, there is no IRFF f i in ES. e impact feature shown in Figure 8(b) is not caused by rolling bearing fault.
In the simulation signal case above, the proposed approach can efficiently detect the fault information and accurately determine that the rolling bearing is under inner ring defect, while the EEMD-FSK approach and the WPT-FSK approach cannot realize such functionality.

Artificial Single Fault Signal Case. Case Western Reserve
University discloses the rolling bearing vibration signal under different operation conditions to the outside world [32].Figure 9 depicts the overall appearance of the test platform.e experiment objective is SKF 6203-2RS rolling bearing, and Table 1 illustrates its structural parameters.Aiming at simulating the early fault signal of the rolling element, a single-point dent is machined on the surface of the rolling element by adopting electron discharge machining.
e damaged diameters are divided into 0.1778 mm, 0.3556 mm, and 0.5334 mm.To show the superiority of the proposed approach, the vibration data corresponding to the least degree (0.1778 mm) of fault are selected for analysis.
e shaft rotary speed refers to n � 1478 r/min, acceleration sensors collect the experimental data, and the sampling frequency refers to f s � 12000 Hz.Shock and Vibration e parameters shown in Table 1 and the shaft rotary speed are introduced into Equation ( 20), where the REFF f e ≈ 118 Hz can be obtained.
e TW of the single fault experimental signal is displayed as Figure 10(a), where no obvious periodic impact characteristics can be seen.A further envelope demodulation is conducted on of the analyzed signal.ere are no frequency components related to rolling element fault in the ES (Figure 10(b)).
e proposed approach is conducted to analyze the single fault experimental signal.e optimization range of BW is set as σ ∈ [357, 12000], and the optimization range of CF is set as c ∈ [178.5, 11821.5].WCA is initialized and 50 iterations are performed.Figure 11(a) shows the optimization curve of the function value with evolution algebra.It can be seen that the minimal value is 118.2 obtained after the 8 times' iteration calculation.e parameter combination consisting of CF and BW corresponding to the minimal value is chosen to construct the optimal ARLW filter.Figure 11(b) illustrates the FS of the optimal ARLW filter, where the CF � 4390 Hz and BW � 911 Hz. e TW of the filtered signal acquired from the optimum ARLW filter is displayed as Figure 11(c).As seen, the filtered signal has and obvious periodic impact feature.e FS and the ES of the filtered are presented in Figures 11(d) and 11(e), respectively.As shown in Figure 11(e), the REFF f e and its frequency doubling 2f e are detected productively.
Aiming at validating the superiority of the proposed approach, the EEMD-FSK approach and the WPT-FSK approach are conducted to analyze the single fault experimental signal.e reconstruct signal obtained after EEMD processing is displayed as Figure 12(a), and its ES is displayed as Figure 12(b).As seen, no useful information is detected.Figure 12(c) presents the kurtogram of the reconstruct signal.As seen, the maximum kurtosis is 0.1, corresponding to the bandpass filter (CF � 2062.5 Hz, level � 7) that is circled by the red dotted line.Figure 12(d) depicts the filtered signal with the maximal kurtosis, where the impact component in the filtered signal has no obvious periodicity.Envelope demodulation is conducted on the filtered signal.Figure 12(e) illustrates the ES of the filtered signal, where the REFF f e cannot be detected.e impact    In the artificial single fault signal case above, the proposed approach can efficiently detect the fault information and accurately determine that the rolling bearing is under the rolling element defect, while the EEMD-FSK approach and the WPT-FSK approach cannot realize such functionality.
e analysis results of Section 5.1 and Section 5.2 illustrate the following conclusions: (1) SEFER index can overcome the influence of random impulse noise, thus it can accurately identify the periodic impact in the signal.10

Shock and Vibration
Combining ARLW filter with WCA, the optimal resonance frequency band of rolling bearing can be determined fast and accurately.(2) Kurtosis index is easily disturbed by random noise.Although the noise reduction method is the preprocessing of FSK method, the resonance frequency band cannot be accurately found because the random noise is unavoidable.

Artificial Compound Fault Signal
Case. e rolling bearing compound fault test is conducted on the QPZZ-II rotating machinery fault test-bed.Figure 14(a) depicts the overall appearance of the test-bed.e rotating shaft of the test-bed is connected with the driving motor through a coupling and a pulley.A pressure-loading device and two disks are fixed on the rotating shaft.Normal bearing is fixed on the middle bearing pedestal, and the defective bearing is fixed on the right bearing pedestal.e experiment objective is SKF 6203 rolling bearing, and Table 2 illustrates its structural parameters.e grooves with a width of 0.2 mm and a depth of 1.5 mm are, respectively, machined on the inner and outer ring by wire cutting technology to simulate the compound fault of rolling bearing.e PCB piezoelectric acceleration sensor is fixed on the right bearing pedestal to the collected vibration signal.Figures 14(b) and 14(c) present the damage rolling bearing and the sensor installation location, respectively.During experiment, the driving motor speed refers to n � 1466 r/min and the sampling frequency refers to f s � 12800 Hz.
e parameters shown in Table 2 and   Shock and Vibration

Shock and Vibration 13
the 16 times' and 9 times' iteration calculation, respectively.Figure 16(c) depicts the FS of the optimal ARLW filter for detecting the inner race defect information, where the CF � 5114 Hz and BW � 619 Hz.Aiming at validating the superiority of the proposed approach, the EEMD-FSK approach and WPT-FSK approach are conducted to analyze the compound fault experimental signal.
e reconstruct signal obtained after EEMD processing is displayed as Figure 17(a), and its ES is displayed as Figure 17(b).As seen, the ORFF f o and its According to the above, the filtered signal obtained by the EEMD-FSK or WPT-FSK approach can only increase the SNR, while it still cannot express the periodic impact feature of the weaker fault.
In the artificial compound fault signal case above, the proposed approach can efficiently detect the fault information and accurately determine that the rolling bearing is under inner ring and outer ring defect, while the EEMD-FSK approach and the WPT-FSK approach can only determine that the rolling bearing is under the outer ring defect.
e analysis result of this section illustrates that the SEFER index with directivity can identify resonance frequency bands caused by different fault types according to different fault characteristic frequencies.However, in the absence of random impact interference, the kurtosis index only can identify the resonance frequency band caused by the periodic impact corresponding to the strong fault in the signal.

Life Cycle Fault Signal Case.
e experimental data are obtained from the life cycle acceleration test of rolling bearings in the center of the NSFI/UCR intelligent maintenance system [19].Figure 19 presents the overall appearance of the test-bed.Four ZA2115 rolling bearings are installed on the rotating shaft of the test bench, and Table 3 illustrates the structural parameters.During experiment, the radial load is about 2671N, the shaft rotary speed refers to n � 1478 r/min, the experimental data are collected by ICP  e RMS variation tendency of the bearing life cycle (0∼9790 min) is shown in Figure 20.As seen, after about 5100 minutes, the RMS began to increase, but the amplitude of fluctuation was not large.is stage is generally called the early stage of bearing fault.After about 7020 minutes, the RMS began to change significantly and reached the extreme value at 9790 min.In this stage, the bearing fault gradually aggravated until the limit of life.

Radial load
Aiming at proving verifying the validity of the proposed approach for the early fault of rolling bearing, the 530th data file (measured at 5300 minute) is selected for analysis.e life cycle of the bearing is about 7 days, and the collection time of the selected data file is 3 days earlier than that of the final test shutdown.e data points of the analyzed signal are 4096.Figures 21(a) and 21(b) show the TW and ES of the analyzed signal, respectively, where no periodic impact and no fault feature can be seen.
e proposed approach is conducted to analyze the life cycle fault signal.e optimization range of BW is set as σ ∈ [0, 200], and the optimization range of CF is set as c ∈ [0, 100].WCA algorithm is initialized, and 50 iterations are performed.Figure 22(a) shows the convergent curve of the function value versus the iterative number.It can be seen that the minimum value is 645.5 obtained after the 7 times' iteration calculation.e parameter combination consisting of CF and BW corresponding to the minimal value is chosen to construct the optimal ARLW filter.Figure 22(b) illustrates the FS of the optimal ARLW filter, where the CF � 8370 Hz and BW � 2869 Hz. e TW of the filtered signal acquired from the optimal ARLW filter is displayed as Figure 22(c), where the periodic impact characteristics can be seen.Figures 22(d Aiming at validating the superiority of the proposed approach, the EEMD-FSK approach and the WPT-FSK approach are conducted to analyze the life cycle fault signal.e reconstruct signal obtained after EEMD processing is displayed as Figure 23(a), and its ES is presented as Figure 23(b).As seen, no useful information is detected.e TW and ES of the frequency band signal corresponding to node (4, 5) are depicted as Figures 24(b) and 24(c), respectively.As seen, although the impact feature can be seen in TW, there is no ORFF f o in ES.
In the life cycle fault signal case above, the proposed approach can efficiently detect the fault information and accurately determine that the rolling bearing is under outer ring defective, while the EEMD-FSK approach and the WPT-FSK approach cannot realize such functionality.
Compared with the simulation fault signal and artificial fault signal, the life cycle fault signal can better reflect the evolution process of bearing fault occurrence and development in actual working conditions.erefore, the analysis result of this section can further reflect the practicability of the proposed approach in the early fault detection of rolling bearing.

Conclusions
is paper came up with a rolling bearing diagnosis approach based on ARLW filtering with the help of WCA, which can detect the fault information by demodulation analysis of the determined resonance frequency band.e proposed approach has three improvements: (1) ARLW is used as a bandpass filter to process the initial signal, which can filter out noise efficiently and extract more fault information; (2) WCA is employed to optimize the ARLW parameters in parallel adaptively, which can avoid artificial interference and improve robustness; (3) e proposed SEFER index can evaluate the quality of the bandpass filters constructed for different fault types.
e proposed approach and the traditional resonance demodulation methods, such as EEMD-FSK approach and WPT-FSK approach, are all conducted to analyze the simulation signal, the artificial single fault signal,    Shock and Vibration the artificial compound fault signal, and the life cycle fault signal.e analysis results of the single fault signal show that the proposed approach can efficiently detect the early fault information, while the EEMD-FSK approach and the WPT-FSK approach cannot realize such functionality.e analysis result of the compound fault signal illustrates that the proposed approach can detect both strong and weak fault information, while the EEMD-FSK approach and the WPT-FSK approach can only detect the strong fault information.e analysis result of the life cycle fault signal further reflects the advantages of the proposed approach in the early fault detection of rolling bearing.
In summary, the proposed approach can efficiently detect the fault feature information and accurately determine the rolling bearing fault type, which is valuable in the engineering industry.

Figure 4 :
Figure 4: e diagnosis process of the proposed approach.
(a) presents the reconstruct signal, and Figure7(b) presents its ES.As seen, no useful information is detected.en, the reconstruct signal is analyzed by FSK, and the kurtogram is displayed as Figure7(c).e color block of kurtogram represents the kurtosis of the narrowband signal obtained by different bandpass filters.As seen, the maximal kurtosis is 0.2, corresponding to the bandpass filter (CF � 384 Hz, level � 6), which is circled by the red dotted line.Figure 7(d) depicts the filtered signal acquired by FSK, where the impact component in the filtered signal has no obvious periodicity.Envelope demodulation is conducted on the filtered signal.Figure 7(e) indicates the ES of the filtered signal, where the IRFF f i cannot be detected.e impact feature shown in Figure 7(d) is not caused by rolling bearing fault.

Figure 6 :Figure 5 :
Figure 6: Diagnosis results of the simulation signal by the proposed approach: (a) iteration curves for determining the optimal ARLW parameters; (b) FS of the optimal ARLW for detecting the inner ring defect information; (c) TW of the filtered signal obtained by (b); (d) FS of (c); (e) ES of (d).

Figure 7 :
Figure 7: Diagnosis results of the simulation signal by EEMD-FSK: (a) TW of the reconstruction signal obtained by EEMD; (b) ES of (a); (c) kurtogram; (d) TW of the filtered signal obtained by FSK; (e) ES of (d).

Figure 8 :
Figure 8: Diagnosis results of the simulation signal by WPT-FSK: (a) kurtogram; (b) TW of the filtered signal with the maximal kurtosis; (c) ES of (b).

Figure 9 :
Figure 9: Test platform of Case Western Reserve University.

Table 1 :
Structural parameters of SKF6023-2RS.Ball diameter Pitch diameter Ball number Contact angle 6Figure12(d) may be caused by random noise or background noise.Figure13(a) presents the kurtogram generated after performing WPT-FSK on the original signal, where node (4, 5) has the maximal kurtosis.Figures13(b) and 13(c) depict the TW and ES of the frequency band signal corresponding to node (4, 5), respectively.As seen, although the impact feature can be seen in TW, there is no REFF f e in ES. e impact feature shown in Figure13(b) also may be caused by random noise or background noise.

Figure 11 :Figure 10 :
Figure 11: Diagnosis results of the single fault experimental signal by the proposed approach: (a) iteration curve for determining the optimal ARLW parameters; (b) FS of the optimal ARLW for detecting the rolling element defect information; (c) TW of the filtered signal obtained by (b); (d) FS of (c); (e) ES of (d).

Figure 12 :Figure 13 :Figure 14 :
Figure 12: Diagnosis results of the single fault experimental signal by EEMD-FSK: (a) TW of the reconstruction signal obtained by EEMD; (b) ES of (a); (c) kurtogram; (d) TW of the filtered signal obtained by FSK; (e) ES of (d).
Figure16(i), the IRFF f i and its frequency doubling 2f i are identified productively.As seen from Figure16(j), the ORFF f o and its frequency doubling 2f o −3f o are identified productively.Aiming at validating the superiority of the proposed approach, the EEMD-FSK approach and WPT-FSK approach are conducted to analyze the compound fault experimental signal.ereconstruct signal obtained after EEMD processing is displayed as Figure17(a), and its ES is displayed as Figure17(b).As seen, the ORFF f o and its

Figure 16 :Figure 17 :
Figure 16: Diagnosis results of the compound fault experimental signal using the approach proposed in the paper: (a), (b) iteration curves for determining the optimal ARLW parameters; (c), (d) FS of the optimal ARLW for detecting the inner and outer ring defect information; (e), (f ) TW of the filtered signal obtained by (c) and (d); (g), (h) FS of (c) and (d); (i), (j) ES of (c) and (d).

Figure 18 :
Figure 18: Diagnosis results of the compound fault experimental signal by WPT-FSK: (a) kurtogram; (b) TW of the filtered signal with the maximal kurtosis; (c) ES of (b).
) and 22(e) display the FS and ES of the filtered signal, respectively.As shown in Figure22(e), the ORFF f o and its frequency doubling 2f o −3f o are detected productively.

Figure 23 (
Figure23(c)  illustrates the kurtogram of the reconstruct signal.As seen, the maximal kurtosis is 0.2, corresponding to the bandpass filter (CF � 3125 Hz, level � 4.5), which is circled by the red dotted line.Figure23(d) depicts the filtered signal with the maximal kurtosis, where the impact component in the filtered signal has no obvious periodicity.Envelope demodulation is conducted on the filtered signal.Figure23(e) illustrates the ES of the filtered signal, where the ORFF f o cannot be detected.Figure24(a) displays the kurtogram generated after WPT-FSK performing on the original signal, where node (4, 5) has the maximum kurtosis.eTW and ES of the frequency band signal corresponding to node (4, 5) are depicted as Figures24(b) and 24(c), respectively.As seen, although the impact feature can be seen in TW, there is no ORFF f o in ES.In the life cycle fault signal case above, the proposed approach can efficiently detect the fault information and accurately determine that the rolling bearing is under outer ring defective, while the EEMD-FSK approach and the WPT-FSK approach cannot realize such functionality.Compared with the simulation fault signal and artificial fault signal, the life cycle fault signal can better reflect the evolution process of bearing fault occurrence and

Figure 22 :Figure 23 :
Figure 22: Diagnosis results of the life cycle fault using the approach proposed in the paper: (a) iteration curves for determining the optimal ARLW parameters; (b) FS of the optimal ARLW for detecting the outer ring defect information; (c) TW of the filtered signal obtained by (b); (d) FS of (c); (e) ES of (d).

Figure 24 :
Figure 24: Diagnosis results of the life cycle fault signal by WPT-FSK: (a) kurtogram; (b) TW of the filtered signal with the maximal kurtosis; (c) ES of (b).

Table 3 :
Bearing structure factor of ZA2115.