Research on Bearing Fault Identification of Wind Turbines’ Transmission System Based on Wavelet Packet Decomposition and Probabilistic Neural Network

: In order to improve the reliability and life of the wind turbine, this paper takes the rolling bearing in the experimental platform of the wind turbine as the research object. In order to obtain the intrinsic mode function (IMF) of each fault type, the original signals of different fault states of the rolling bearing on the experimental platform are decomposed by using the overall average empirical mode decomposition method (EEMD) and the wavelet packet decomposition method (WPD), respectively. Then the energy ratio of the IMF component of the different types of faults to the total energy value is calculated and the eigenvectors of different types of faults are constructed. The extreme learning machine (ELM) and probabilistic neural network (PNN) are used to learn fault types and eigenvector samples to identify the faults of the rolling bearing. It is found that the bearing fault characteristics obtained by the WPD method are more obvious, and the results obtained by the same recognition method are ideal; and the PNN method is obviously superior to the extreme learning machine method in bearing fault recognition rate.


Introduction
With the continuous progress of the industrial level, the size of the wind turbine continues to increase, coupled with the relatively harsh working environment of the wind turbine, the corresponding faults will also increase.According to the survey, the failure rate of blades is 17%, the failure rate of the generator is 28%, and the failure rate of the gearbox is the highest, which is 57% [1].Once the gearbox fails, the fault detection priority of the bearing is the highest, and the fault of other parts will be analyzed only after it is confirmed that there is no problem with the bearing.The problems of wind turbines' bearing include contact fatigue, wear, defects, dents, and corrosion failures [2].The continuous abrasion leading to the failure of the bearing to operate properly is referred to as wear failure [3].Contact fatigue of bearings is the change in material structure due to repeated stress application.It begins as pitting with a slight shape and size, but with the expansion of pitting comes fatigue spalling [4].Therefore, the accuracy of bearing fault detection and fault type identification of the wind turbine drive system is very important.
Fault diagnosis is divided into two steps.Firstly, signal processing is carried out using time domain, frequency domain and time-frequency domain analysis methods to extract features; then, the support vector machine and other artificial intelligence methods are used for fault diagnosis [5].The time domain analysis method can only reflect the changes of signal waveform and reveal the internal characteristics of the signal.The frequency domain analysis method can reflect the different frequency components of the signal, but it is not suitable for processing non-stationary vibration signals.The time-frequency domain analysis method can analyze the vibration signal from two different aspects of time domain signal, but it is not suitable for processing non-stationary vibration signals.The time-frequency domain analysis method can analyze the vibration signal from two different aspects of time domain and frequency domain, so it can be used to analyze the non-stationary mechanical vibration signal.The time-frequency domain analysis method based on signal decomposition is also widely used in the field of rolling bearing fault diagnosis, such as empirical mode decomposition (EMD), vibrational mode decomposition (VMD), local mean decomposition (LMD) and wavelet packet decomposition (WPD).At present, the fault diagnosis methods of rolling bearing based on artificial intelligence mainly include support vector machine (SVM), extreme learning machine (ELM), convolutional neural network (CNN) and auto encoder (AE) [6].Hoang et al. can accurately diagnose the fault of rolling bearing by converting the original vibration signal into a two-dimensional gray image and combining the strategy of the convolutional neural network [7].Peng et al. proposed a one-dimensional convolutional neural network model containing residual blocks.The model uses wide convolution and dropout technology to enhance the generalization ability of the network [8].Jin proposed a rolling bearing fault diagnosis method based on a two-dimensional image convolutional neural network (CNN) structure [9].Zhao et al. proposed a rolling bearing diagnosis method with one-dimensional CNN [10].Hao et al. proposed an end-to-end pipeline for bearing fault diagnosis based on 1D-CLSTM [11].Chen et al. proposed a novel short-term electric load forecasting method EMD-ELM which is based on empirical mode decomposition (EMD) and extreme learning machine (ELM).And the forecasting results of the EMD-ELM were proved to be better than all the other methods [12].Zhang et al. found that the arc fault diagnostic method with the combination of EMD and ELM identifies arc fault for various loads effectively [13].To sum up, the above-mentioned fault diagnosis method has achieved great success at this stage.But most equipment works under normal operating conditions [14].
In this paper, the frequency of different fault types of rolling bearing in the wind turbine drive system is studied.Through the experimental platform of the wind turbine drive system, the original signal data of bearing in the normal state, cage fault, rolling element fault, inner ring fault and outer ring fault are collected.The eigenvectors of different fault types are constructed by using the EEMD and WPD method, and the recognition accuracy of fault types by using the ELM and PNN method is compared and analyzed.It is found that the WPD method is better than the EEMD method, and PNN is better than the ELM method.

Research Object
This paper uses the rolling bearing in the gearbox of the dynamic test platform for the planetary gear transmission system as the research object.The rolling bearing model is ER-16 K.The number of rolling elements of the bearing is 9, the diameter of the rolling element is 7.9375 mm, and the bearing pitch diameter of the bearing is 38.5064 mm.The structural diagram is shown in Figure 1.In Figure 1, D is pitch diameter, d is rolling element diameter, r 1 is inner race radius, r 2 is outer race radius, α is contact angle, Z is number of rolling elements, f i is inner race rotation frequency, f o is outer race rotation frequency, f c is cage rotation frequency, and f b is rolling element rotation frequency.The speed at point i on the inner race is as follows [15]: The speed at point i on the outer race is as follows: The speed at point i on the cage is as follows: Therefore, the characteristic frequency of cage failure is as follows: For a single rolling element, when the outer ring fails, the frequency of the rolling element passing through the outer ring is as follows: When the inner ring fails, the frequency of the rolling element passing through the inner ring is as follows: Because V = θ • r, the conversion frequency is inversely proportional to the radius, and the scale factor is d/2r 1 ; therefore, Thus, the rotation frequency of the rolling element relative to the cage is obtained, that is, the failure frequency of the rolling element is as follows: The inner and outer rings of the bearing have a relative rotation frequency which is f r = f i − f o , because the outer ring is fixed, so f o = 0, f r = f i .When there are Z bearing rolling element, the fault characteristic frequency of each part can be obtained as follows.
Inner ring fault frequency: Outer ring fault frequency: Rolling element failure frequency: Energies 2024, 17, 2581 4 of 16 Cage failure frequency: The frequency corresponding to different frequency spectrum of bearing fault types is different.The more serious the fault is, the greater the corresponding amplitude is.

The Experimental Platform
This paper uses the bearing in the gearbox of the dynamic test platform for the planetary gear transmission system as the research object, and the experimental platform of the transmission system is shown in Figure 2. In Figure 2, the experimental platform is composed of an industrial computer, a data-acquisition system, a controller, a drive motor, a planetary-gear box and a magnetic power brake.The controller controls the rotation frequency of the drive motor.The frequency displayed by the controller is the theoretical rotation frequency of the motor.The torque is controlled by the magnetic power brake.The data-acquisition card of the data-acquisition system is a four-channel data acquisition card.And the acceleration sensors are installed in three mutually perpendicular directions of the gearbox.
Rolling element failure frequency: Cage failure frequency: The frequency corresponding to different frequency spectrum of bearing fault types is different.The more serious the fault is, the greater the corresponding amplitude is.

The Experimental Platform
This paper uses the bearing in the gearbox of the dynamic test platform for the planetary gear transmission system as the research object, and the experimental platform of the transmission system is shown in Figure 2. In Figure 2, the experimental platform is composed of an industrial computer, a data-acquisition system, a controller, a drive motor, a planetary-gear box and a magnetic power brake.The controller controls the rotation frequency of the drive motor.The frequency displayed by the controller is the theoretical rotation frequency of the motor.The torque is controlled by the magnetic power brake.The data-acquisition card of the data-acquisition system is a four-channel data acquisition card.And the acceleration sensors are installed in three mutually perpendicular directions of the gearbox.In the experiment of the bearing fault signal acquisition, bearing with different fault parts are installed in the gearbox for data acquisition.Set the sampling frequency as 5000 Hz, the sampling time as 1 s, and the torque as 22 N•m.Collecting the original signals under five different working conditions: normal state, cage fault, rolling element fault, inner ring fault and outer ring fault of the bearing, and 10 groups of data for each state are collected.

The Fault Feature Extraction of Bearing Based on the EEMD Method
Since EMD uses cubic spline interpolation to fit the curve, frequency aliasing will occur, and the error will become larger.Therefore, the EEMD method was adopted to extract fault features of the bearing which is a noise aided data analysis method [16].In the experiment of the bearing fault signal acquisition, bearing with different fault parts are installed in the gearbox for data acquisition.Set the sampling frequency as 5000 Hz, the sampling time as 1 s, and the torque as 22 N•m.Collecting the original signals under five different working conditions: normal state, cage fault, rolling element fault, inner ring fault and outer ring fault of the bearing, and 10 groups of data for each state are collected.

The Fault Feature Extraction of Bearing Based on the EEMD Method
Since EMD uses cubic spline interpolation to fit the curve, frequency aliasing will occur, and the error will become larger.Therefore, the EEMD method was adopted to extract fault features of the bearing which is a noise aided data analysis method [16].
The EEMD method adds a uniformly distributed Gaussian white noise background to the original signal, and the spectrum of Gaussian white noise is uniformly distributed.Therefore, the original signal can be decomposed into continuous IMF of different scales to achieve the effect of suppressing mode aliasing.
Firstly, the EEMD method needs to set a reasonable overall average decomposition number N and the ratio of Gaussian white noise to the standard deviation of the original signal ε.The larger the value of N is, the smaller the noise error is, and the better the decomposition effect is, but the calculation time is doubled.If the value of the ratio of Gaussian white noise to the standard deviation of the original signal ε is too large, the components with low frequency will overlap.If the value of ε is too small, the components with high frequency will overlap.Therefore, it is necessary to make reasonable requirements for the determination of N and ε.
When the expected error is fixed to e = 0.01, the EEMD method is used to decompose the original signal.After decomposition, the amplitude standard deviation of the first signal is σ 1 , and the ratio of this value to the amplitude standard deviation σ 0 of the original signal is ζ [17].
Determine parameters ε according to white noise addition criteria.And then N obtained by using the viewpoint given by Norden E. Huang [17].
Taking the rolling element fault signal as an example, the original signal is decomposed by the EEMD method, and the amplitude standard deviation of the original signal and the first signal after decomposition is calculated.The results are shown in Table 1.Finally, ε is 0.12, N is 200.Assuming that the original signal is x(t), the combined signal x'(T) can be obtained by adding Gaussian white noise g(T) to the original signal [18].
Then, the IMF can be obtained through decomposition of the combined signal with the EEMD method.Decompose the cycle N times to obtain N IMF, and take its mean value as the final IMF: where, i is the number of cycles.
The EEMD method is used to decompose the original signal of the rolling bearing under the rolling element fault, and the decomposition results are shown in Figure 3.
From Figure 3, it can be seen that there are a number of IMF after decomposition by the EEMD method, but there are still meaningless components in IMF components.The last components are low-frequency components and residual components.If these components are also used as feature elements, it will cause interference and reduce the accuracy of the fault identification algorithm.
Therefore, the correlation coefficient method is used to filter out unimportant IMF components.According to the threshold value of the set coefficient, if the coefficient of the IMF component is smaller than the threshold value, it will be eliminated as an unimportant component.The calculation formula is as follows:  From Figure 3, it can be seen that there are a number of IMF after decomposition by the EEMD method, but there are still meaningless components in IMF components.The last components are low-frequency components and residual components.If these components are also used as feature elements, it will cause interference and reduce the accuracy of the fault identification algorithm.
Therefore, the correlation coefficient method is used to filter out unimportant IMF components.According to the threshold value of the set coefficient, if the coefficient of the The correlation coefficient of the IMF component of the rolling bearing under the rolling element fault after the calculation is shown in Table 2.The IMF components with a correlation coefficient less than 0.02 in Table 2 are eliminated, and the first eight IMF components are retained as useful components.Different fault states of rolling bearings have different frequency ranges, and different frequency ranges correspond to different frequency band energies.Therefore, the frequency band energies under different fault states are regarded as the characteristics of fault diagnosis.
The calculation formula of IMF band energy is as follows [19]: where, N is the number of sampling points; n is the sequence of sampling points, n = 1, 2, 3 . . .n; and IMF n,i is the i-th IMF component.The total band energy of the signal is as follows: where, M = 8.Taking the energy ratio of IMF component to total energy as the feature element of the feature vector to construct the feature vector, the ratio of each IMF component to total energy is as follows: The fault feature vectors of the rolling bearing under the rolling element fault are shown in Table 3.In the same way, the fault feature vectors of the rolling bearing in normal state, cage fault state, inner ring fault state and outer ring fault state can be obtained.Each state contains 10 groups of feature vectors.Some results are shown in Table 4. Wavelet packet decomposition (WPD) is used to decompose the fault signal into three layers, and the ratio of the signal energy in the decomposition frequency band to the total signal energy is also taken as the fault feature vector element to construct the fault feature vector of each fault type of the bearing.The schematic diagram of signal decomposition frequency band division by the WPD method is shown in Figure 4.
In the same way, the fault feature vectors of the rolling bearing in normal state, cage fault state, inner ring fault state and outer ring fault state can be obtained.Each state contains 10 groups of feature vectors.Some results are shown in Table 4.

The Fault Feature Extraction of Bearing Based on the WPD Method
Wavelet packet decomposition (WPD) is used to decompose the fault signal into three layers, and the ratio of the signal energy in the decomposition frequency band to the total signal energy is also taken as the fault feature vector element to construct the fault feature vector of each fault type of the bearing.The schematic diagram of signal decomposition frequency band division by the WPD method is shown in Figure 4. Let sequence { } Z n n h ∈ satisfy the following [20]:  Let sequence {h n } n∈Z satisfy the following [20]: Define a set of recursive functions ω n ∈ L 2 (R) (n = 1, 2, 3,. ..), which are generated by scaling function φ(t) and wavelet function ψ(t).The relationship is as follows: where, g k = (−1) k h 1−k .When n = 0, ω 0 (t) and ω 1 (t), respectively, correspond to φ(t) and ψ(t).The sequence {ω n (t)} n∈Z generated by the formula is called the wavelet packet determined by the basis function ω 0 (t) = ϕ(t).
The fast orthogonal wavelet packet transform in multiscale analysis is introduced into the wavelet packet algorithm.The detail signal d in wavelet space is further decomposed, Energies 2024, 17, 2581 9 of 16 and the wavelet packet decomposition formula and reconstruction formula of signal x(t) are obtained.They are as follows: According to the multiresolution analysis relation L 2 (R) = ⊕ j∈Z W j , j ∈ Z, the decomposition relation in wavelet packet subspace W n j is obtained as follows: The wavelet packet decomposes wavelet subspace W j (j = 1, 2, • • •) step by step, and the decomposition expression is as follows: It can be seen from the above formula that the wavelet packet realizes the further decomposition of the high-frequency sequence.If j = 0, it indicates the original signal x(t) itself at the resolution j level, and is recorded as x 1 ; if x 1 is decomposed once, the decomposed signal signals x 2 and x 3 will be obtained at the first layer of wavelet packet decomposition; and if x 1 is decomposed twice, the decomposed signals x 4 , x 5 , x 6 and x 7 will be obtained at the second layer of wavelet packet decomposition, and so on.
Taking the rolling element fault signal as an example, the rolling element fault signal is decomposed by the WPD method.The decomposition diagrams of the third level are shown in Figure 5.
After the rolling element fault signal is decomposed by the WPD method, the fault feature vectors are obtained by taking the ratio of each frequency band energy to the total energy as an element.As shown in Table 5.Similarly, the fault feature vectors of the inner ring, outer ring, cage fault signal and normal signal are obtained, and some calculation results are shown in Table 6.
After obtaining the fault feature vectors of the five states of the bearing, these fault feature vectors can be used as samples to identify the fault state.After the rolling element fault signal is decomposed by the WPD method, the fault feature vectors are obtained by taking the ratio of each frequency band energy to the total energy as an element.As shown in Table 5.

The Fault Pattern Recognition Based on ELM 4.1. The Extreme Learning Machine (ELM)
The ELM can solve the problems of slow training speed, the difficulty reaching the global minimum, and the large impact of learning rate on the network.The structure of a typical single hidden layer feedforward neural network is shown in Figure 6.

The Extreme Learning Machine (ELM)
The ELM can solve the problems of slow training speed, the difficulty reaching the global minimum, and the large impact of learning rate on the network.The structure of a typical single hidden layer feedforward neural network is shown in Figure 6. Figure 6 shows that the network is composed of the input layer, hidden layer and output layer.Suppose there are n neurons, there are n input characteristic variables, and there are m target characteristic variables in the output layer, and there are l neurons in the middle hidden layer.When there is a training set with Q samples, the matrix is n * Q, and the rows represent features, that is, suppose each sample has n features, and the columns represent samples.For the specific calculation principle, refer to [21].H is the output matrix of the hidden layer.
( ) Figure 6 shows that the network is composed of the input layer, hidden layer and output layer.Suppose there are n neurons, there are n input characteristic variables, and there are m target characteristic variables in the output layer, and there are l neurons in the middle hidden layer.When there is a training set with Q samples, the matrix is n * Q, and the rows represent features, that is, suppose each sample has n features, and the columns represent samples.For the specific calculation principle, refer to [21].H is the output matrix of the hidden layer.
where T ′ is the transpose of matrix T. H is the matrix of Q * l, that is, Q samples multiplied by l neurons.Each row represents the output of a sample on l neurons, and each column represents the output of each neuron to Q samples.

The Experimental Analysis
The input part of the fault feature vectors of the bearing under different fault conditions into the ELM randomly select six groups of data for each state as training samples for learning, and the remaining samples as test samples for identification.Since the training samples and test samples are randomly selected each time, the results of each test are not necessarily the same.The results of 10 tests of the fault feature vectors extracted by EEMD are shown in Table 7.It can be seen that the average recognition accuracy of the ELM for the fault feature vectors extracted by EEMD is 93.5% from Table 7.This is obtained in the case of a small number of experimental data sets.In practical applications, because the collected data is more complex, it is likely to be lower than this accuracy.Several fault type identification errors are shown in Figure 7.
feature vectors extracted by the ELM is concentrated between 3 and 5, that is, between cage fault and outer ring fault, while the rolling element fault, inner ring fault and normal state can be identified more accurately.However, this method did not achieve the expected effect.In 10 tests, there are only three times when the samples of the test set are accurately classified, which is obviously not suitable for application to actual conditions which are more complex.From Figure 7, it can be found that the fault type identification errors of the fault feature vectors extracted by the ELM is concentrated between 3 and 5, that is, between cage fault and outer ring fault, while the rolling element fault, inner ring fault and normal state can be identified more accurately.However, this method did not achieve the expected effect.In 10 tests, there are only three times when the samples of the test set are accurately classified, which is obviously not suitable for application to actual conditions which are more complex.
The ELM method is used to test the fault feature vectors obtained by WPD.The results of 10 tests are shown in Table 8.It can be seen that the average recognition accuracy of the ELM for the fault feature vectors extracted by WPD is 98.5%.The recognition accuracy of these data is 5% higher than that of the EEMD extracted by the fault feature vectors, and the lowest accuracy is 95% in these 10 tests.Therefore, it can be explained that the effect of the fault feature vectors obtained by the WPD method is better than that obtained by the EEMD method.Among Table 8, there are three fault-type identification errors, as shown in Figure 8.It can be seen that the three misclassification situations are 4 pairs of 3, 5 and 2 misclassifications, that is, the inner ring fault is incorrectly identified as the cage fault, outer ring fault and the rolling element fault, respectively, in the three processes.It shows that this method is not ideal for fault identification.
Among Table 8, there are three fault-type identification errors, as shown in Figure 8.It can be seen that the three misclassification situations are 4 pairs of 3, 5 and 2 misclassifications, that is, the inner ring fault is incorrectly identified as the cage fault, outer ring fault and the rolling element fault, respectively, in the three processes.It shows that this method is not ideal for fault identification.

The Probabilistic Neural Network (PNN)
The PNN was introduced by Specht in the early 1990s which is the probabilistic neural network [22].The PNN can be regarded as a radial basis neural network, which is a  The PNN was introduced by Specht in the early 1990s which is the probabilistic neural network [22].The PNN can be regarded as a radial basis neural network, which is a completely forward calculation process without reverse error calculation.And the PNN has the following advantages: fast operation, good stability, fault tolerance, fast convergence, will not fall into the local optimal situation, can use the linear learning algorithm to achieve the function of the nonlinear learning algorithm, and very suitable for pattern recognition [23].The PNN method is a kind of feedforward neural network.It is generally composed of an input layer, mode layer, summation layer and output layer [24], and its structure is shown in Figure 9.
Energies 2024, 17, x FOR PEER REVIEW 14 of 17 completely forward calculation process without reverse error calculation.And the PNN has the following advantages: fast operation, good stability, fault tolerance, fast convergence, will not fall into the local optimal situation, can use the linear learning algorithm to achieve the function of the nonlinear learning algorithm, and very suitable for pattern recognition [23].The PNN method is a kind of feedforward neural network.It is generally composed of an input layer, mode layer, summation layer and output layer [24], and its structure is shown in Figure 9.The value from the training sample is input to the input layer, which transfers the eigenvector to the network [25].The input layer is connected with the pattern layer through the weight ω ij , and the pattern layer is calculated to obtain the matching relationship between the eigenvector and different patterns.The sum of all the training samples is equal to the number of neurons in the pattern layer.The input samples are weighted and summed in the mode layer, and then transferred to the summation layer after the nonlinear operator g z j operation.
The output of each mode unit in this layer is as follows: where X is the input sample; W i is the connection weight matrix between the input layer and the mode layer; and δ is the smoothing factor [26].
The summation layer accumulates the probabilities belonging to a certain class.
where m is the number of samples.The different states of the rolling bearing correspond to a summing layer unit, respectively, and the outputs of the corresponding mode layer units are added [26].
The output layer is a competitive neuron.Different neurons correspond to different states of bearings.They receive various density function outputs by the summation layer.The number of neurons in the output layer is equal to the number of types of fault types.The neuron with the largest probability density function outputs 1, and the other neurons output 0.

The Experimental Analysis
Part of the fault feature vectors of the bearing under different fault conditions are input into the PNN model.Six groups of data are randomly selected for each state as training samples for learning, and the remaining samples are used as test samples for identification.Since the training samples and test samples are randomly selected each time, the results of each test are not necessarily the same.The results of 10 tests of the fault feature vectors extracted by EEMD are shown in Table 9.The two inaccurate identifications are shown in Figure 10.In the Figure 10, it can be seen that the misclassification of the features extracted by the PNN for EEMD is also concentrated between 3 and 5, that is, between the cage fault and the outer ring fault.The rolling element fault, inner ring fault and normal state can be accurately identified.However, it can be seen from Table 9 that the average recognition accuracy rate of the PNN for the fault feature vectors extracted by EEMD is 99%.Among the 10 times of recognition, there are 8 times that can accurately identify the fault state, and 2 times can also reach 95%.
The PNN model is used to test the fault feature vectors obtained by WPD.The results of 10 tests are shown in Table 10.
From Table 10, it can be found that the test results of the PNN identifying the fault feature vectors obtained by WPD have reached 100% 10 times, that is, when the identification samples are collected in the laboratory, this method can accurately identify the fault type of the rolling bearing.Although the collected data will be more complex in the actual working condition, the effect of this method is expected to achieve the effect of preliminary judgment of the fault type.feature vectors extracted by EEMD are shown in Table 9.The two inaccurate identifications are shown in Figure 10.
In the Figure 10, it can be seen that the misclassification of the features extracted by the PNN for EEMD is also concentrated between 3 and 5, that is, between the cage fault and the outer ring fault.The rolling element fault, inner ring fault and normal state can be accurately identified.However, it can be seen from Table 9 that the average recognition accuracy rate of the PNN for the fault feature vectors extracted by EEMD is 99%.Among the 10 times of recognition, there are 8 times that can accurately identify the fault state, and 2 times can also reach 95%.The PNN model is used to test the fault feature vectors obtained by WPD.The results of 10 tests are shown in Table 10.From Table 10, it can be found that the test results of the PNN identifying the fault feature vectors obtained by WPD have reached 100% 10 times, that is, when the identification samples are collected in the laboratory, this method can accurately identify the fault type of the rolling bearing.Although the collected data will be more complex in the actual working condition, the effect of this method is expected to achieve the effect of preliminary judgment of the fault type.

Conclusions
In this paper, the rolling bearing in the experimental platform of the wind turbine gear transmission system is taken as the research object.The EEMD method and the WPD method are used to decompose the original signals of different fault states of the experimental rolling bearing and construct the fault feature vectors.Then, the ELM and PNN are used to learn fault types and the fault feature vector samples to identify fault types.
(1) Using the same recognition method, it is found that the fault recognition accuracy of the fault feature vectors after the WPD method is higher, that is, the fault characteristics obtained after the WPD method are more obvious.(2) The recognition accuracy of the features extracted by EEMD and WPD using the ELM method were 93.5% and 98.5%, respectively; the recognition accuracy of the PNN method for the bearing fault is 99% and 100%, respectively.The PNN method is obviously superior to the ELM method.(3) Based on the results, the accuracy of the method of combining WPD and PNN to predict bearing fault types is relatively high.

Figure 1 .
Figure 1.Structural diagram of the rolling bearing.Figure 1. Structural diagram of the rolling bearing.

Figure 1 .
Figure 1.Structural diagram of the rolling bearing.Figure 1. Structural diagram of the rolling bearing.

Figure 2 .
Figure 2. Experimental platform of wind turbine drive system.

Figure 2 .
Figure 2. Experimental platform of wind turbine drive system.

Figure 3 .
Figure 3. IMF component of rolling element under rolling element fault.

Figure 3 .
Figure 3. IMF component of rolling element under rolling element fault.

Figure 4 .
Figure 4. Schematic diagram of signal decomposition frequency band division.

Figure 4 .
Figure 4. Schematic diagram of signal decomposition frequency band division.

Figure 5 .
Figure 5.The decomposition diagrams of the third level of the rolling bearing under the rolling element fault.

Figure 5 .
Figure 5.The decomposition diagrams of the third level of the rolling bearing under the rolling element fault.

5 .
The Fault Pattern Recognition Based on the PNN 5.1.The Probabilistic Neural Network (PNN)

Figure 9 .
Figure 9. Structure diagram of PNN.The value from the training sample is input to the input layer, which transfers the eigenvector to the network[25].The input layer is connected with the pattern layer through the weight ij ω , and the pattern layer is calculated to obtain the matching relationship between the eigenvector and different patterns.The sum of all the training samples is equal to the number of neurons in the pattern layer.The input samples are weighted and summed in the mode layer, and then transferred to the summation layer after the nonlinear operator ( ) j z g operation.

Table 1 .
The amplitude standard deviation of the rolling bearing.

Table 2 .
The correlation coefficient of IMF.

Table 3 .
The fault feature vectors of the rolling bearing under the rolling element fault.

Table 4 .
The partial fault feature vectors of the rolling bearing under different fault conditions.
3.2.The Fault Feature Extraction of Bearing Based on the WPD Method

Table 4 .
The partial fault feature vectors of the rolling bearing under different fault conditions.

Table 5 .
The fault feature vectors of the rolling bearing under the rolling element fault.

Table 6 .
The partial fault feature vectors of the rolling bearing under different fault conditions.

Table 5 .
The fault feature vectors of the rolling bearing under the rolling element fault.

Table 7 .
Accuracy of fault type identification.

Table 8 .
Test results of wavelet packet decomposition feature extraction in the ELM recognition.

Table 9 .
Test results of features extracted by the PNN recognition EEMD.

Table 10 .
Test results of features extracted by the PNN recognition WPD.

Table 9 .
Test results of features extracted by the PNN recognition EEMD.

Table 10 .
Test results of features extracted by the PNN recognition WPD.