Signal de-noising in gear pitting fault identification by an improved singular value decomposition method

In this research a new method of improved singular value decomposition (ISVD) is proposed for the vibration signal de-noising of gear pitting fault identification. In this method, the delay time τ and embedding dimension m of the Hankel matrix for SVD are optimized by autocorrelation function and Cao’s algorithm respectively. Simulation and experiments are conducted to demonstrate the method. In the simulation, the ISVD method is employed to de-noise the artificial vibration signal in a mathematical model of gear pitting fault, the result demonstrates the signal-noise ratio (SNR) value is SNR = 31.3 dB, and the root-mean-square error (RMSE) value is RMSE = 0.34. In the experiment, the ISVD method is adopted to de-noising the vibration signal of gear pitting fault identification, the results demonstrate SNR is SNR >45 dB, and the RMSE value is RMSE <0.4 of the fault characteristic signals at each measuring point position. The results of simulation and experiment show, the ISVD method is efficient to de-noise the vibration signal of gear pitting fault.


Introduction
Gears are critical transmission components in rotating machinery with extensive application field especially in modern industrial areas. The whole machinery's performance and production efficiency are directly affected by the working state of gears [1]. The stable operation is directly affected by the state of gear motion whether the gear run normally or not [2,3]. Therefore, it is very crucial to identify the early faults before the gears develop into serious state. However, in practical applications, it is difficult to predict the formation and development of gear fault in advance because while the vibration around the gearbox is perceived, the gear fault has been a severe condition. Up to now, the research of the gear failure mechanism in early and middle condition is not mature. For the practical applications, it is important to develop the research of gear failure on vibration characteristics of gear box.
Recently, many theoretical types of research are developed on the vibration of gearbox by gear failure, few type are based on practical application. And there is a certain gap between theory and practical application. Generally, fault detection of gear teeth is the collection of fault signal samples site, and the identification of the gear fault type according to the information by the signal collection [4]. Since the measured signal is interfered with strong noise background, de-noising for the extraction of useful characteristics from the original signal is the first important stage. At present, the common methods of signal de-noising are various, such as Empirical mode decomposition method (EMD) [5,6], Local mean decomposition method (LMD) [7,8], Singular value decomposition method (SVD) [9][10][11] wavelet analysis method [12][13][14]. In practical, the methods are limited by some problems. EMD is affected by mode aliasing which is occurred frequently in the decomposition process, and the quality of signal decomposition is reduced in the condition [15,16]. LMD also generates mode aliasing which makes the decomposition and extraction of fault features incomplete, and lead to decomposition failure on the whole time scale [17][18][19]. The effect of Wavelet de-noising method is depended on the calculation and selection of reasonable parameters [20,21].
Beside the de-noising methods above, the singular value decomposition (SVD) is a better method to deal with problems of mode aliasing and endpoint effect. The singular value of effective feature is distinguished from original signal by SVD which is based on the singular value matrix of the characteristic signal. In the method, characteristic information of the noise signal is eliminated, the useful characteristic signal is reconstructed. As it is invariant, stable and effective for the de-noising, SVD is applied to the filtering of nonlinear signals, such as practical application of de-noising in fault identification of gear fault [22]. However, the accuracy of de-noising in fault identification by SVD is depended on the calculations of τ and m, which is highly subjective and unreliable in theory at present. Although suitable parameters of τ and m are calculated by a new method which is called k-SVD, the determination of k is restricted by other conditions, and affects the de-noising effect of SVD [23,24].
It is urgency to reasonably determinate the values of τ and m in SVD respectively. For the value of τ, although is appropriate to any value in theory, the determination in practical applications is directly affected by the vibration signal. In the case that τ is too small, the correlation is concentrated to between different elements of the delay vector, and the reconstructed phase space is meaningless. Contrary, while τ is too large, it leads to loss of noise delay coordinates of mutual information elements, and the reconstructed phase space is unable to reflect the initial dynamic behavior. The parameter m controls the phase space, determines the accuracy of de-noising and identification. While m is large, the phase space is too large, and makes the calculation tedious and complex. Contrary, while m is small, the reconstructed information is unrealizable for the extraction from mutation signals [25,26].
In order to improve the SVD for de-noising, an improved singular value decomposition (ISVD) is developed based on autocorrelation function with Cao's algorithm for the parameter of τ and m. For the determination of τ, the autocorrelation function is closely related to the shape of the attractor in the reconstructed phase space and approximated the principal element direction of the signal system. Empirically, while the decay of autocorrelation function is approximated to 1/e, the optimal delay time is obtained [35]. For the determination of m, Cao's algorithm which is based on the False Nearest Neighbor (FNN) method, makes the calculation of m only related to τ, and overcomes the influence of various thresholds on m [27,28].
In this paper, ISVD is proposed for the vibration signal de-noising of gear pitting fault identification. Firstly, the parameter of τ and m of the Hankel matrix for SVD are optimized by autocorrelation function and Cao's algorithm respectively. Secondly, simulation of ISVD by autocorrelation function and Cao's algorithm is conducted to verify the proposed method. Thirdly, the experiment of tooth surface pitting failure is conducted by ISVD method to verify the effectiveness of the proposed method.   Reconstruct the signal K the phase space of the effective characteristic signals and noise signals are obtained respectively; then the optimal de-nosing order is calculated by the singular value difference spectrum (DSSV) method of the characteristic signal; finally, the inverse operation of Hankel matrix is used to reconstruct the characteristic signal after de-noising. It retains the effective signal value and rid the noise signal value, then reconstructs the effective signal in a suitable order. Essentially, SVD is a method of matrix orthogonalization in the mathematical view, which decomposes the given matrix into two matrixes U m×m and V n×n , as shown: [29][30][31][32].

Theory of the SVD method
where S is the singular value matrix of the real matrix H, in S, σi is the singular value of the real matrix H, and For the signals x = .x 1 ,x 2 , ,x m /, the de-noising is described as follows: x 1 +2 x 1 +n : : : : : : : : : : : : where m is the embedded dimension, τ is delay time, N is the length of the signal, and n = N -(m -1) × τ.
While the difference spectrum of singular value is positioned on the mutation point, the spectrum value di is the maximum in the whole sequence Di. It is expressed by: Based on of Eq. 4, the sequence number corresponding to the maximum SVDS dmax = Max (Di) is calculated of the signal. According to the above principle, the new singular value matrix SNew is obtained from the effective characteristic signal, as shown in Eq. 5: where, σ H is the singular value left over.
(vi) Reconstruction of characteristic signals X H i , show as follow:

Characteristic parameters of ISVD method
In this section, the two parameters of τ and m are optimized by the combination of autocorrelation and the advantages of Cao's algorithm.

The parameter of delay time τ
A measured time series is x(t), and the autocorrelation function of the sequence after normalization processing is shown in Eq. 7 [33,34]: where, N is sample size, τ is delay time, N x is average value of the samples, The research shows, the corresponding time is selected while the value of the autocorrelation function exactly decayed to 1/e by the autocorrelation function method, as the optimal value of the τ of the effective phase space feature vector [35].  (i) On the definition, an m-dimension "embedding space", the time series vector E X m .i / at the ith phase point is defined as where, x(i) is the signal corresponding to the time series; i is the time-series sequence number, i = 1, 2, ..., N -(m -1)τ. Then the Euclidean distance Rm(i) from E X m .i / to the nearest neighbor E X NN m .i / is calculated as: where, k k is the 1-norm, and all the normalizations followed are 1-norm.
(ii) While the m-dimensional phase space is extended to m + 1, similarly, the ith phase point is E X m+1 .i /, and the Euclidean distance R m+1 .i / of E X m+1 .i / to the nearest neighbor point E X NN m+1 .i / is calculated: (iii) The relationship between E(m) and m is shown: (iv) Continue increasing the dimension of the phase space m1, and make it larger than m, then repeat the calculation of the first three steps. Based on the results of each calculation, the relationship between E1(mi) and mi is charted. While E1(mi) in the graph reaches a plateau within a certain threshold, the corresponding mi is the optimal m of the E X m .i /. (v) In practical applications, a criterion is defined to judge the smoothness of the finite sequence change of E1(mi) and mi, it is shown as: From the equations, in a range of mi, if E2(m) = 1 the signal is random; if E2(m) within a certain threshold are not all equal to 1, the signal is the deterministic signal.

Evaluation of de-noising effect
The definition of SNR is expressed, shown as in follow: where, z(i) is the de-noising signal.
The definition of RMSE for root-mean-square error evaluation is expressed, shown as in follow: where, N is the length of x(i).

Random noise ξG-noise(t)
The mean is 0 and the variance is 1

Numerical simulation
The other parameters are shown in Table 1. Base on the signal, the operation of de-nosing by the simulation is listed as:  Meanwhile, in the singular value difference spectrum sequence, the first three large difference spectra are d1 = 10.39, d2 = 13.74, and d3 = 0.12 respectively, the maximum difference spectrum is d2. (iv) Determine de-noising order Fig. 7 is the simulation results of the relationship between σ and DSSV. In the case that σ1 and σ2 are selected for the reconstructed signal, many information is lost; in the case of selecting σ1, σ2 and σ3, it contains a little noise information; in the case of  selecting σ1, σ2, σ3, and σ4, the reconstructed signal is distorted; In the case of selecting from σ1 to σ5, a large amount of noise information is mixed in the reconstructed signals. As a result, the selection of σ1 σ2 and σ3 is suitable because the effective information of the obtained characteristic signal remains relatively completely and has a good effect of de-noising.
The signals are decomposed of f0(t) by the principle of the SVD method, and the σ1 to σ4 are obtained. Combined with Cao's algorithm, the optimum parameters are m = 5 and τ = 8. Finally, as shown in Fig. 8, the inverse solution of singular value is used to reconstruct the effective characteristic signal after de-noising. From the figure the reconstructed effective signal is proximity to the pure signal, which indicates that ISVD has the potential for de-noising.
As shown in Table 2, after the de-noising by ISVD. And the SNR = 31.31 dB, RMSE = 0.34 are acquired by Eqs. 15 and 16. It reminds ISVD has a good de-noising potential. Fig. 9a is shown the experiment layout. Fig. 9b is the photograph of experiment. On the outer side of the gearbox, four measurement points are set as Points I, II, III and IV with four uniaxial acceleration sensors for the collection of vibration signal. The speed of input is 1500 r/min; the input current of the magnetic powder is P = 38 W; the modulus of gear is mo = 2 mm. The gear material is 18CrNiMo7 steel. Fig. 10 is shown the waveforms of vibration signal which were collected at four measurement points. The signals are affected by a large numbers of interference signals. Table 3 is the calculation results. The minimum value of τ is τ = 8, which is calculated by the autocorrelation function method as while τ = 8 at each measurement point, the autocorrelation function C(τ) 1/e.

The ISVD of experimental data
Based on the calculation above, the parameters are set as: τ = 8 and m = 6. Fig. 12 is the singular value of the experimental data which is decomposed at the measurement Points I, II, III, and IV by ISVD method. Fig. 12a is the singular value distribution graph. While m < 2, the singular values of the four measurement points are σ > 800, while m > 2, the singular value in each point is smooth stability to a certain value as σI = 900, σII = 900, σIII = 700, σIV = 900 the mutated point is m = 2 or m = 3. Fig. 12b is the Singular value difference spectrum to supplement the determination of reconstruction order. In the figure, while m = 3, 5 and 6, the difference spectrum values are approximately close to 50. With the combination of results by Fig. 12a,b, the optimum reconstruction order is DSSV = 3.  Fig. 13 is the reconstructed signal with fault characteristics of four measuring points in time-domain by Eq. 6. In the figures, the green wave is original signal, the brown wave is the de-noising signal. From the figures, the impact vibration generated by backlash and transmission error is dominant and effectively reduced after reconstruction. From the figures, the acceleration value of de-noising signal is much smaller than original one in each time. In addition, the trend of de-noising signal and original signal are same. The effect of signal de-noising is evaluated by SNR method and RMSE, the results are showed in Table 4. From the results, SNR ≥45 dB, and RMSE <0.4 in all measurement points. It reveals a better effect of de-noising by ISVD. In addition, it reveals that the de-noising effect of signals at Points I and II is better than Points III and IV.

Experimental results and analysis
Fig. 14 is the frequency response diagram of de-noising signals and original signals in each measurement point from Fig. 13 by Fourier transform. In Fig. 14, Fig. 14a,c and d are similar, because Point 1 Point III and Point IV are meshed. In the figures, the effective signals are retained in the frequency range of 300 Hz < f < 800 Hz and 1.2 kHz < f < 1.7 kHz as the different rotation speeds of each axis, signals of other frequencies are noise which is denoised by ISVD. Fig. 14b are the signals in Point II which near the input source and affected by more vibrations. In the figure, the signal of noise is de-noised in ineffective range of frequency, and the effective signal is retained.

Conclusions
In this research, a new method of ISVD is proposed for the vibration signal de-noising of gear pitting fault identification. Simulation and experiments have been conducted to prove the efficiency of ISVD. The conclusions are listed as follows: (i) The theory of ISVD is expressed, the method for the determination of parameters τ and m for SVD are optimized by autocorrelation function and Cao's algorithm respectively. The signal after de-noising is reconstructed by the inverse of the Hankel matrix. (ii) In the simulation, the ISVD method is employed to de-noise the artificial vibration signal in a mathematical model of gear pitting fault, the results are SNR = 31.3 dB and RMSE = 0.34. The simulation results reveal that ISVD has a good effect on de-noising of the fault vibration signal. (iii) In the experiment, the ISVD method is adopted to de-noising the vibration signal of gear pitting fault identification, the results demonstrate SNR >45 dB, and RMSE <0.4. The experiment results reveal that ISVD has a good effect on de-noising of the fault vibration signal. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4. 0/.