Intelligent Analysis Method of Gear Faults Based on FRWT and SVM

An intelligent analysis method for gear faults based on fractional wavelet transform (FRWT) and support vector machine (SVM) is proposed. Based on this method, FRWT is used to eliminate noise from the gear vibration signal, and the vibration signal after noise elimination is carried thought wavelet packet decomposition and reconstruction. A sequence corresponding to the signal is constructed consisting of the module with the highest level wavelet coefficients after decomposition and feature vectors corresponding to the energy sequence whichwere obtained by calculation.Then, a particle optimizationmethod is used to optimize SVM parameters, and the feature vectors as training samples are input into SVM for training while the test samples are input for fault recognition. Experimental results show that the gear fault analysis method proposed in this paper is able to effectively extract the weak fault signal. The accuracy rate for identification of the type of gear fault reached 96.7%.


Introduction
As an indispensable component in modern mechanical structure to transmit power and motion, the gear is widely used in various equipment due to its large load capacity, high transmission accuracy, and constant transmission power.As power-transmission equipment develops toward larger scale, greater complexity, automation, and continuous use, gear faults may induce increasingly large losses.Therefore, it is necessary to develop an intelligent gear fault analysis method that can improve the efficiency and accuracy of detection and reduce cost [1][2][3].
Gear faults are always found in the vibration signal produced by meshing gears.This vibration signal also contains a large number of nonstationary signal components.Using classical spectral analysis methods will result in averaged spectrum values across the whole time scale, which inadvertently conceals some details of mutations.In contrast, a timefrequency method can be used effectively for analysis of nonstationary signals, which can fill the gap in the standard FFTbased method and is developing rapidly as a new approach in the fault analysis field.Wavelet transform (WT) [4][5][6] is considered as an ideal time-frequency method capable of providing localized information of time and frequency domains at the same time, which can help identify the mutations in the vibration signal.But wavelet transform has limitations, since WT of various scales is similar to a group of frequency-domain band-pass filters [7] which are limited in only time and frequency domain, which makes WT an imperfect choice for chirp type signals whose energy is not optimally gathered in the frequency domain and which exist widely in nature and in various artificial situations.Therefore, analysis methods like fractional Fourier transform (FRFT) [8][9][10], short-time fractional Fourier transform [11], Radon-Wigner transform [12], and fractional wavelet transform (FRWT) [13,14] are proposed one after another.FRFT has some unique characteristics but cannot present localized features.Short-time fractional Fourier transform has the fault of low resolution.Radon-Wigner transform is based on quadratic transformations, which have cross terms.FRWT was originally proposed by Mendlovic in 1997, who applied FRFT to the signal to obtain a fractional spectrum, which was then transformed by traditional WT.As a generalized form of WT, FRWT merged many advantages of WT and FRFT 2 Shock and Vibration such as no cross terms, linear characteristics, multiresolution analysis, and fractional domain presentation.Thus, FRWT has attracted more attention as a potential application in signal processing field.
Taking the gear vibration signal as the researched object, two key problems are studied in the practical engineering application, including the gear fault feature extraction and the pattern recognition.An intelligent gear fault analysis method based on a combination of FRWT and SVM is proposed.FRWT makes use of the WT's multiresolution analysis theory and FRFT's time-frequency focusing properties to filter the gear vibration signal.So it can ensure the consistency and effectiveness of time-frequency resolution.Then wavelet packet transform was applied to the denoised signal extracting energy values in the frequency band to build feature vectors, which can be used for training and learning by SVM to further realize intelligent identification of the faulty gear.The SVM method has the advantage of identifying any sample according to the limited samples information.It overcomes the shortcoming of poor classification accuracy caused by insufficient samples in faults of detection of gears.So it improves the efficiency and accuracy of faults of detection of gears.
Combining FRWT and SVM, the method can quickly achieve the detection of finished gears.Then qualities and faults of gears are determined according to the detection results.All in all, this method will provide a powerful analysis means for the diagnosis of gears quality.

Fundamental Principles
Here,   (, ) is the transforming kernel function: Here,  is the order of FRFT, which is an arbitrary real number. represents a rotation angle;  = /2. is an integer.(⋅) is an impulse function.The inverse transformation of FRFT is

Fractional Wavelet Transform.
Fractional wavelet transform combines wavelet transform and fractional theory, which promotes multiresolution wavelet analysis theory to a time domain and generalized frequency domain.Meanwhile, it merges characteristics of wavelet transform and fractional theory, which makes it more powerful in signal processing.Mendlovic first proposed the definition of fractional wavelet transform [13,15] and applies it to signal (): Here, ℎ *  () is wavelet basis function and   (,   ) is kernel function and can be presented as follows: Here,   is  order Hermite polynomials,  = /2, and  is the order of fractional wavelet transform  () (, ), in its value ranging from 0 to 1.The reconstruction formula of fractional wavelet transform is We can see that the combination and decomposition of fractional wavelet transform are actually the consistence transform in the time-frequency domain, which accurately regulates information within time domain and frequency domain.

Support Vector
Machine.The advantage of support vector machine (SVM) [16][17][18] applied to gear vibration signal analysis is that it fits for a small sample of the decision and can maximally excavate the classification knowledge using limited feature information, which has a great practical significance for gear fault analysis.Therefore, we use the SVM in this study to establish the gear fault model by the identification and analysis of the gear fault type from the filtered signal.
The main idea of SVM is to create a hyperplane as a decision interface that has the ability to correctly classify all the training samples and maximize the distance from the nearest point in the training sample to the classification face.The geometry of the hyperplane in the linearly separable two-dimensional input pattern space is shown in Figure 1.
Here, the pentacle and rectangles represent the two types of samples,  is the classification lines, and straight lines  1 and  2 are the nearest samples to the corresponding classification line and parallel to the classification lines.The distance between  1 and  2 is classification interval.The sample on  1 and  2 is the SVM.Assume that the training sample is {(  ,   )}  =1 , in which input is   ∈   and output is   ∈ {−1, 1}.The hyperplane equation is    +  = 0.  is the adjustable weight vector,  is weight function, and  is the offset.And the constraint condition of the optimization problem is The unknown problem is min: This formula is a mathematical inequality problem solved through constrained optimization.The point that makes the equation balance is called support vector.In order to solve the mathematical optimization problem, Lagrange multipliers are introduced to convert it into the dual problem to solve.
And the unknown problem is max: To solve the optimal coefficients above â, ŵ, b, the optimal classification function is The certain nonlinear mapping can map the sample space into a high-dimensional feature space to convert the liner nonseparable pattern which becomes linearly separable and to construct the optimal classification hyperplane in the high-dimensional feature space.The nonlinear mapping function is called a kernel function.Therefore, the choice of appropriate inner product kernel function can map the low-dimensional to high-dimensional space to achieve linear classification without increasing computational complexity.Compared with linearly separable patterns, the constraints (10) become (parameter  which is used to control the degree of penalty to the wrong sample points) The optimal classification function becomes The Gauss radial primary kernel function (  ,   ) = exp(−‖  −   ‖ 2 ) is used in this paper, where  is a constant.As the particle swarm optimization method is faster and more efficient than the grid search method [19], this particle swarm optimization method is used in this paper to optimize the SVM classification model.
The mathematical description of a particle swarm optimization algorithm is as follows: assuming that in a  dimensional optimization space there are a group of particles , where the speed for the th particle is V  = (V 1 , V 2 , . . ., V  ), whose location is   = ( 1 ,  2 , . . .,   ), and then the optimal position for the th particle is   = ( 1 ,  2 , . . .,   ), and the optimal position of the whole population [20] is (the particle update formula is as follows) where  = 1, 2, . . .,  and  = 1, 2, . . ., .  is the number of iterations,  1 ,  2 are accelerating constant,  1 ,  2 are a random number in [0, 1], and the  is inertia weight.

Intelligent Analysis Method of
Gear Fault Detection Here,   (V  ),   (V  ), and   (V  ) are separately the discrete fractional wavelet transform of (), (), and (), and  is the transform order.The basic steps are as follows: (1) Use a sensor to collect the gear vibration signal () in the acceleration stage, select the range and step order, and apply the fractional Fourier transform to signal acquisition, forming a two-dimensional distribution of the signal energy in the fractional Fourier domain (, ).
(2) In the plane formed in step one, peak points are found through two-dimensional search, where the maximum peaks' corresponding FRFT domain is taken as the best FRFT domain while the optimal order  opt is determined according to (17), and  opt order fractional Fourier transform is applied to signal (): (3) For the signal in the fractional Fourier domain, its wavelet basis function is determined, with proper decomposition layers  selected.After the th wavelet decomposition, the fractional wavelet signal is obtained, including the th low frequency coefficients   and high frequency coefficients from 1st to th layer  1 ,  2 , . . .,   .
(4) The signal was decomposed to 1st to th fractional wavelet domain signal, where adaptive threshold filter was applied to the corresponding high frequency coefficients (i.e. 1 ,  2 , . . .,   ) to eliminate noise.
(5) Wavelet reconstruction was applied to the processed signal, where the basis function and layer number should be ensured the same way as that for wavelet decomposition.Then, wavelet reconstruction is applied to the  layers' frequency coefficients (  ) and high frequency coefficients ( 1 ,  2 , . . .,   ) obtained in the steps above.
(6) − opt order fractional Fourier transform was applied to the wavelet reconstruction signal to get the timedomain signal after noise removal.(2) Solve the energy   on th layer within various frequency bands:

Energy Feature Extraction for
(3) Solve the total energy of feature signals.A group of energy sequences {  ,  = 1, 2, . . ., 2  } corresponding to the signal are obtained after wavelet packet decomposition, where the total energy is (4) Construct the feature vector.As the gear faults will have a great impact on the energy of the signal in each frequency band and   may result in some difficulty in data analysis when the value is large, the feature vector is constructed after frequency band energy is normalized:

Intelligent Recognition for SVM.
The normalized feature vectors in Section 3.2 are used as the sample inputs to support vector machines for training and recognition.In this paper, an open source SVM pattern recognition and regression package LIBSVM is used, which utilizes a particle swarm optimization method to optimize the LIBSVM classification model, and then tests the penalty factor  and kernel function parameter  finally obtaining proper SVM correction parameters.The main steps for particle swarm optimization (PSO) to optimize the parameters of SVM for intelligent fault identification are as follows: (1) Use an initialized particle swarm algorithm, in the search space.Twenty (20) particles are randomly generated, randomly determining each particle's initial position and initial velocity.Initialize error penalty parameters  and Gaussian kernel parameters  of support vector machines, where the optimal ranges of  and  are [0.1 ∼ 100] and [0.1 ∼ 1000], respectively.
(2) The initialized SVM parameter will be used in the SVM algorithm to build a corresponding model, which is used to predict and classify the test sample.The fitness value of each particle will be calculated according to the fitness function.
(3) The initial particle fitness value is used as the model's initial optimal solution, which is compared with the global optimization objective function value.If the initial fitness value of particles is better than the optimal value of the objective function, the former shall be used as the current optimal objective function value, going on to search for the global optimal solution.
Acceleration sensor (4) Update formula of the particle velocity and position is used to update and obtain the best position  best of each particle itself.After comparing all the best positions of the particles, the best position  best of all particles is obtained. (

Experiments
In this paper, the rolling-gear inspection machine is used to carry out the experiment, in which the arrangement of acceleration sensor used to collect gear vibration signal is shown in Figure 2. The gear drive is a one-level drive, whose driving wheel and a driven wheel share a gear tooth number of 17, with a transmission ratio of 1 : 1.The driving wheel rotation speed is 750-1000 rpm.The sampling frequency of the vibration signal sound is 10 kHz.The gear vibration signals were collected in three cases of normal gear, shaft unbalance fault, and pitting fault.When gears are in states of normal meshing, shaft unbalance fault, and pitting fault (Figure 3), respectively, curves of measured vibration signal in the time domain and frequency domain are shown in Figures 4(a), 4(b), and 4(c) separately, where obvious modulation phenomena occur in the time-domain signals for all of three types of gears mentioned above.From the FFT frequency spectrum diagram in Figure 3, background noise has masked the noise generated by gear meshing, which makes the time and frequencydomain curves of the three types of gears so similar that it is difficult to distinguish between them based only on a comparison of the curves in either the time domain or in the frequency domain.
A gear fault analysis method based on a combination of FRWT and SVM proposed in this paper is used to analyze the three kinds of gear vibration signals mentioned in Figure 4.
Step 1. Fractional wavelet transform is applied to the collected gear vibration signal to eliminate background noise, extracting the useful signal.Fractional Fourier transform is applied to the gear vibration signal with a transform order of .According to the peak search and calculation, an optimal order can be determined, by the basis on which a narrowband filter is constructed to mask the optimal FRFT domain signal and further to eliminate the energy of background noise while retaining that of useful vibration signal.Therefore, the gear vibration signal after denoising based on the fractional wavelet transform is obtained.Figure 5 is the 2-dimensional distribution of signal energy before filtering and after narrow band filtering in the fractional Fourier transform domain in three cases, showing that the filtering method proposed in this article can achieve better noise suppression effects.
Step 2. Fractional Fourier transform is applied to the signal filtering to get useful gear vibration signal (), which is decomposed and reconstructed by db6 3-layer wavelet packet and to obtain the decomposition energy spectrum.Frequency of gear vibration signal ranges from 0 to 2000 Hz, and sampling frequency   is 10 kHz.Therefore, when -layer wavelet packet decomposition is applied to the gear vibration signal, the node bandwidth at the th layer is where frequency range for each node is Here,  represents the node number at the th layer in wavelet decomposition.Wavelet packet frequency band and frequency correspondence are shown in Table 1.
A sequence corresponding to the signal which consists of the module of the highest level wavelet coefficients after decomposition and feature vectors corresponding to the energy sequence were obtained according to (20). Figure 6 shows the histogram for the feature vectors of the three cases of normal gear meshing, shaft unbalance fault, and pitting fault.

Shock and Vibration
As is shown in Figure 6, in vibration signal wavelet energy spectrum for different types of gear fault, energy is concentrated mainly within the 1st band because of the vibration due to stiffness stimulation and load variation.When the gear fault exists, the amplitude of the meshing frequency and its frequency band tend to produce large growth.However, it is difficult to directly observe the effects of background noise in the FFT spectrum.
As can be seen clearly in the wavelet energy spectrum diagram, the characteristic value of the 1st frequency band for shaft unbalance gear fault is 2% lower than that of normal gear.At the same time, the characteristic value of the 3rd frequency band is increased by about 90% and the characteristic value of the 6th frequency band is decreased by about 99%.
As for the pitting gear fault, the characteristic value in the 1st band decreases by 29% compared with the normal gear operation.Meanwhile, the characteristic value of the 3rd frequency band is increased by about 97% and the characteristic value of the 6th frequency band is increased by about 40%.The characteristic values of the other bands also appear to increase in varying degrees.
The above analysis shows that there are marked differences in characteristic values of the three cases of normal gear, shaft unbalance fault, and pitting fault, which can be used as feature vectors for gear fault type recognition.

Conclusions
This paper studies a gear fault detection method based on the feature vector of gear vibration signals.For the collected vibration signal, the noise is eliminated by the FRWT method.Using the wavelet packet transform, the feature vectors are constituted by the frequency band energy.On this basis, the feature vectors are trained and classified by support vector machines.Through the experiments, the FRWT can extract the useful signals which reflect the gear fault features, achieving good effects in elimination of background noise.And the automatic recognition of gear fault is realized through the clustering and classification of vibration signals vector by SVM, so that it can solve the problem of faults recognition in time and frequency domain which is caused by the interference of the background noise and the complexity of the side band.Therefore, it is intelligent, efficient, and accurate for the gear vibration analysis method based on FRWT and SVM.
Wavelet Packet Frequency Band.Wavelet packet transform is applied to multiple decompositions and coefficient reconstruction for the filtered signal of the fractional Fourier transform.As the energy projection sequence of gear vibration signal in various frequency bands is a product of the gear transmission mode, analysis of gear vibration signals in different frequency bands can help inconspicuous signal frequency features stand out in a more obvious way as energy change in several subspaces of various resolution.Signal energy after wavelet packet decomposition in each frequency band in various scale spaces is lined up according to their scale as feature vectors.Steps for extracting wavelet packet signal features are as follows:(1) The -layer wavelet packet and coefficient reconstruction are carried out by filtered vibration signal after FRWT, of which the feature signals  , () are decomposition coefficients of filtered vibration signal on the  scale:  = 1, 2, . . . = 1, 2, . . .,   = 0, 1, 2, . . ., 2 −1 .Here  is the time-domain location; data length of original signal is ; data length of feature signal is  = /2  .

Figure 2 :
Figure 2: Layout of acceleration sensors in the experimental stage.

Figure 5 : 2 -
Figure 5: 2-dimensional distribution of signal energy before filtering and after narrow band filtering on fractional Fourier transform domain in three cases.

Figure 7 :Figure 8 :
Figure 7: Optimal results of SVM via particle optimal method.
Check whether the iteration conforms to the ending condition.If a prescribed error requirement or a certain number of iterations are reached, the iteration will be stopped.Otherwise, go back to step (2).
(6)The penalty parameter  and Gauss kernel parameter  obtained through optimization are used as SVM parameters.Input training samples for training and respective category labels are established according to different types of gear faults.Meanwhile, "one against rest" SVC is used for training to ultimately obtain some certain training parameter information (including support Vector, pull-down factor, and network bias.).(7) The test samples are input into support vector machines for recognition, outputting the test results, which are used to determine gear fault types.

Table 1 :
A wavelet packet frequency band-frequency tablet.