Next Article in Journal
Constructing a Measurement Method of Differences in Group Preferences Based on Relative Entropy
Previous Article in Journal
Perturbative Treatment of the Non-Linear q-Schrödinger and q-Klein–Gordon Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Misalignment Fault Diagnosis of DFWT Based on IEMD Energy Entropy and PSO-SVM

School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Entropy 2017, 19(1), 6; https://doi.org/10.3390/e19010006
Submission received: 27 September 2016 / Revised: 17 November 2016 / Accepted: 30 November 2016 / Published: 1 January 2017

Abstract

:
Misalignment is an important cause for the early failure of large doubly-fed wind turbines (DFWT). For the non-stationary characteristics of the signals in the transmission system of DFWT and the reality that it is difficult to obtain a large number of fault samples, Solidworks and Adams are used to simulate the different operating conditions of the transmission system of the DFWT to obtain the corresponding characteristic signals. Improved empirical mode decomposition (IEMD), which improves the end effects of empirical mode decomposition (EMD) is used to decompose the signals to get intrinsic mode function (IMF), and the IEMD energy entropy reflecting the working state are extracted as the inputs of the support vector machine (SVM). Particle swarm optimization (PSO) is used to optimize the parameters of SVM to improve the classification performance. The results show that the proposed method can effectively and accurately identify the types of misalignment of the DFWT.

1. Introduction

The development of the wind energy is very rapid with growing concerns about the energy and environmental problems nowadays. The global cumulative installed capacity reached 432,419 MW at a compound growth rate of 17% per year [1]. Wind farms mainly use the large capacity variable speed constant frequency (VSCF) doubly-fed wind turbines (DFWT) as the main models [2]. Its structure mainly includes the rotor, shaft, gearbox, the doubly-fed induction generator (DFIG), the converter, etc. Because of the long transfer chain of DFWT and the installation in the cabin at high altitude of tens of meters, or even hundreds of meters, precise alignment is rather difficult; on the other hand, due to the fluctuation of wind speed, start-stop of the units frequently happened. As time goes on, there will be a shift or deformation in some components, causing the misalignment between the generator and the gearbox [3]. When misalignment happens, gearbox high-speed shaft and generator bearings will produce heavy dynamic load, which will not only increase axial and radial vibration, but also cause the bearing oil leakage, high temperature of bearing and attachment bolts and fastening bolts loosening. The accumulation of eccentric error even causes damage to bearings of the high-speed end and generator. Because misalignment is an important cause of early failure of large DFWT and may cause damage to two core parts including the gearbox and generator, it is necessary to study the misalignment fault diagnosis methods to ensure longstanding and stable running of the DFWT.
At present, there are many methods used in the wind turbines fault diagnosis. For example, the back propagation (BP) neural network based on the fruit flies algorithm optimization was used by Qi Liwan in fault diagnosis for the wind turbine gearbox [4]; a fault diagnosis method based on learning vector quantization (LVQ) neural network was proposed by Ding at al. [5]; the fuzzy clustering analysis method was applied to the wind motor gearbox fault by Li [6]; the fault classification method based on the ART2 neural network and C-average clustering algorithm of wind turbine gearbox was used by Li et al. [7]; an expert system model for the fault diagnosis of a wind turbine cabin was established by Zhang [8]; a classification method based on support vector machine (SVM) was proposed by Liu et al. to classify the operating mode of wind turbines [9]; a wind turbine main shaft bearing fault diagnosis based on SVM was studied by Huang [10]; machine learning and data mining techniques organized in the framework of a general scheme that achieves fault diagnosis of wind turbines were studied by Precup et al. [11]; a hybrid dynamic classifier for drift-like fault diagnosis was applied to wind turbine converters by Toubakh [12]; a mixed Bayesian/Set-membership approach was used by Rosa et al. to solve the problem of fault detection and isolation of wind turbines [13]; the interval nonlinear parameter-varying (NLPV) parity equations were used by Blesa et al. to solve the problem of fault diagnosis of a wind farm [14]; and the interval observer approach and virtual actuators/sensors were applied by Blesa et al. in literature [15], where the problem of Fault Detection and Isolation (FDI) and Fault Tolerant Control (FTC) of wind turbines were addressed. Compared with other intelligent algorithms, SVM has many advantages in learning ability, generalization performance, small sample processing capacity and the nonlinear processing capacity. Therefore, in this paper, SVM is used in the misalignment fault diagnosis.
There are three types of misalignment: parallel misalignment, angle misalignment and the comprehensive misalignment [16]. Due to the short running time of large DFWT in China, the data of misalignment available for diagnostic studies are lacking. Thus, a 3D model was established by Solidworks (Solidworks 2014, version SP0) for a 1.5 MW wind turbine in this paper, and the model is then imported into Adams (Adams 2013, version 2013.1), where the normal working and the misalignment of the system are simulated. The vibration signals (angular acceleration) of the gearbox high-speed shaft are obtained for analysis. The methods for extracting fault features from vibration signals mainly include traditional and modern signal processing methods [17]. The traditional signal processing methods include time domain analysis and frequency domain analysis. Signal stability hypothesis is the basis of the traditional time domain or frequency domain analysis. However, the vibration signals of DFWT usually have the characteristics of being nonlinear and non-stationary. Therefore, the modern signal processing methods should be used to analyze them. Many scholars have done a lot of work in this area. For example, wavelet decomposition was used by Zhang o process the vibration signals of gearbox in DFWT [18]; while wavelet packet was used by Wang [19], a method of multi resolution empirical mode decomposition (EMD) combined with the average in frequency domain was proposed by Ju et al. to extract gear fault frequency [20]; EMD and the envelope spectrum of the intrinsic mode function (IMF) was used by Zheng et al. to identify the fault types of rolling bearing [21] and so on. The advantages and shortcomings of modern signal processing methods now commonly used are summarized in Table 1.
Among these commonly used processing techniques summarized in Table 1, EMD proposed by Hibert–Huang is a signal processing technique to decompose a data set into several IMFs by a sifting process, reducing the coupling between the signal characteristic information. Using the EMD method, the feature information of the original data can be grasped accurately and effectively, which is conducive to mining deep feature [35]. However, there will be divergent phenomenon at both ends of the date when using the EMD method, which may result in the entire data sequence gradually being contaminated and the decomposition effect becoming severely distorted [36,37,38,39]. To solve this, a mirror extension algorithm is used to improve the EMD in this paper. After the vibration signals of DFWT processed by the IEMD, the energy entropy of the IMF are calculated and taken as the inputs of SVM. Meanwhile, PSO is used to optimize the parameters of SVM to get better classification performance. The results show that the proposed method can effectively identify the different misalignment types.

2. The Related Theories

2.1. The Improved EMD

The basic idea of EMD is to decompose a complex non-stationary signal into a finite number of IMF, and the IMF must satisfy the following:
  • The amount of extreme values and passing zeros of the function must be equal to or differ by one at most;
  • The average value of the upper envelope formed by the local maximum points and the lower envelope formed by the local minimum points of the function is zero.
The specific steps of decomposing the original signal x ( t ) by EMD are as follows:
(1)
Find out all the maximum points and the minimum points of the original signal firstly. Then, fit the upper envelope e u p p ( t ) and the lower envelope e l o w ( t ) by the cubic spline interpolation. The average value m 1 ( t ) = e u p p ( t ) + e l o w ( t ) 2 can be obtained. Assume:
w 1 ( t ) = x ( t ) m 1 ( t ) .
(2)
Judge whether w 1 ( t ) meets the two conditions of IMF. If it is satisfied, then w 1 ( t ) is the first IMF component; if not, then take w 1 ( t ) as a new x ( t ) . Repeat step (1) until w 1 k ( t ) = w 1 ( k 1 ) ( t ) m 1 k ( t ) meets the two requirements.
(3)
After the first IMF component c 1 ( t ) is decomposed from the original signal x ( t ) , separate c 1 ( t ) from x ( t ) to get the rest r 1 ( t ) :
r 1 ( t ) = x ( t ) c 1 ( t ) .
Take r 1 ( t ) as the original data and repeat steps (1), (2) and (3). The second IMF component of the original signal can be then obtained. Repeat n times, and r 1 ( t ) , r 2 ( t ) ,…, r N ( t ) can be obtained. The process will last until the remainder r N ( t ) becomes a monotonic function.
Finally, x ( t ) can be expressed as follows:
x ( t ) = 1 N c i ( t ) + r N ( t ) .
The IMF components contain different frequency band elements of the original signal x ( t ) , while the remaining part r N ( t ) represents the central tendency of x ( t ) .
According to the maximum values and minimum values provided by the original signal, cubic spline interpolation is performed to obtain the upper and lower envelope in the EMD processing. However, the two end points of the signal may not be extreme value points, which will produce the fitting error. The error will gradually spread to the inside signal with the decomposition process, causing the waveform of IMF component to become seriously distorted. This is called the end effects of the EMD. In order to maintain the authenticity of the data, the mirror extending method is used to improve the EMD in this paper.
The mirror extending method is based on the characteristics of mirror symmetry mapping. To minimize the side effect, the mirror is placed in the position that ensures that the signal has symmetrical extreme values. The specific steps are:
(1)
To draw the curves at both ends and mark the corresponding extreme values (the maximum and minimum);
(2)
To decide the position of the mirror according to the distribution characteristics of the curves, and make the signal in the mirror and the original signal connected end-to-end to form a closed loop;
(3)
To get the new signal with the length of two times the original data.
By extension, the upper and lower envelope obtained by fitting are completely determined by the internal data, and the end effects are solved in essence.

2.2. The Principles of SVM

SVM was proposed by Vapnik in 1995, which takes the training error as the constraint condition of the optimization problem and minimizes the confidence range as the objective of optimization. SVM is a kind of learning method based on structuring risk minimization principle, with better promotion ability than other traditional methods [40]. Since the solving process of SVM is finally transformed into a quadratic programming resolving, its solution is the only optimal one. SVM shows many unique advantages in solving nonlinear and high-dimensional pattern recognition problems [41,42,43]. The principle of SVM can be briefly described as giving training samples ( x 1 , y 1 ) , ( x 2 , y 2 ) , …, ( x l , y l ) , where x i R n , y i { 1 , 1 } , and l is the number of training samples. They can be separated by a hyper plane, and the hyper plane can be expressed as follows:
ω x + b = 0 ,
where ω and b , respectively, represent the normal vector and the constant term. When ω and b are the best, it means the optimal hyper plane has been found, which makes the distance between the positive and negative samples the largest.
For linearly inseparable problems, the samples can be mapped to the hyper plane by the kernel function K ( x i , y j ) to realize linear separability. Its objective function is as follows:
W ( a ) = i = 1 l a i 1 2 i = 1 l j = 1 l a i a j y i y j K ( x i , x j ) , s .   t . i = 1 l a i y i = 0 ,   0 a i C
where a i is the Lagrange operator; C is the punishment factor, and its basic function is to control the penalty of the wrong samples; K ( x i , x j ) is the kernel function, and its basic function is to transform the vectors of low-dimensional to the inner product in the high-dimensional.
The kernel functions used in SVM are mainly linear kernel function, polynomial kernel function, radial basis kernel function and sigmoid kernel function. For the selection of kernel function, there is no definite criterion. Many studies and experiments show that the radial basis function (RBF) is a better choice when there is not enough prior knowledge [44]. In this paper, the RBF is used as the kernel function. The function is as follows:
K ( x i , x j ) = exp ( g x i x j 2 ) ,   g > 0 ,
where g is the kernel parameter.
After determining the type of kernel function, the penalty parameter C and the RBF parameter g should be decided. Parameter C determines the training error and the generalization ability of the classifier; parameter g affects the distribution form of the samples in the feature space. The selection of them is very important for the performance of the classifier. The traditional grid search is a way to find the optimal parameters, but the search result may not be necessarily good. Thus, in this paper, PSO is used to select the parameters of C and g to improve the classification performance.

2.3. Particle Swarm Optimization

PSO is a kind of swarm intelligence optimization algorithm besides the ant colony algorithm and the fish swarm algorithm, which was first proposed by Kennedy and Eberhart in 1995. PSO is used to initialize a group of particles in the solution space, and each particle represents a potential optimal solution to the optimization problem. The characteristics of the particles are represented by three indexes of position, velocity and fitness value. Fitness value can be calculated by the fitness function, and the size of the value represents the goodness and badness of the particle. The particles’ location are updated by tracking the individual extreme value and the population extreme value when they move in the solution space. The fitness value of each particle is calculated when it updates its location. The principle can be briefly described as follows.
Supposing in a D -dimensional search space, there is a population X = ( X 1 , X 2 , , X n ) , consisting of n particles. The i-th particle is represented as a D-dimensional vector of X i = [ x i 1 , x i 2 , , x i D ] , which represents the position of the i-th particle in the D-dimensional search space, and it also represents a potential solution to the problem. According to the objective function, the fitness value of each particle's position X i can be calculated. If the velocity of the i-th particle is V i = [ V i 1 , V i 2 , , V i D ] , the individual extremum is P i = [ P i 1 , P i 2 , , P i D ] , the population extremum is P g = [ P g 1 , P g 2 , , P g D ] , and then the velocity and position of the particle are updated according to the following equations:
V i d k + 1 = V i d k + c 1 r 1 ( P i d k X i d k ) + c 2 r 2 ( P g d k X i d k ) ,
X i d k + 1 = X i d k + V i d k + 1 ,
where d = 1 , 2 , , D ; i = 1 , 2 , , n ; k is the number of iterations; c 1 and c 2 are learning factors, which are non negative constants, making the particles learn from their own or other better particles to achieve the purpose of being close to the better position of itself or the whole group; r 1 and r 2 are random numbers distributed in [0, 1] to maintain the diversity of group. In order to prevent the blind search of particles, it is generally recommended to limit the position and velocity to a certain range of [ X max , X max ] and [ V max , V max ] .

3. Acquisition of Samples for Diagnosis

A 1.5 MW wind turbine system in a certain factory is studied in this paper. Gearbox, coupling, generator, etc., are modeling and assembling by SolidWorks firstly. Secondly, the models are saved in SolidWorks as Parasolid(.x_t) format. Then, they are imported into Adams. According to the actual status of the wind turbine system, material definition, motion constraints, drive and contact force for components and parts are added. In this paper, four states of normal working condition, parallel misalignment, angle misalignment and the comprehensive misalignment are studied, respectively. The models have been verified before the vibration signals are extracted from the gearbox high-speed shaft.

3.1. Verification of the Models

3.1.1. Normal Working State

At normal working state, the model is verified by the speed of the gearbox, generator shaft and the meshing frequency. The comparison of theoretical results and simulation results are listed in Table 2.
It can be seen the relative error is very small by comparing the theoretical results with the simulation results in Table 2. This demonstrates that the model of normal working state is correct.

3.1.2. At Fault Conditions

The shaft radial force and axial force are used to test the validity of the failure models. In theory [45], the characteristics of misalignment can be summarized in Table 3.
The actual outputs of the failure models would be described below:
  • Parallel Misalignment
    In Adams, the local coordinate system (Marker) can be created in the left half coupling mass center. In addition, the rotation axis of the revolute joint relative to the ground is arranged on a z-axis of the Marker global coordinate system. Then, the Marker global coordinate system is translated a certain distance along the y-axis to make the mechanism rotate. The deviation between the center of mass and the center of rotation can lead to the eccentric mass excitation to the mechanism; thus, the parallel misalignment of the rotor shaft of the generator in this distance can be simulated. Some results of the parallel misalignment are shown in Table 4.
    From Table 4, it can be found that the fundamental frequency and the frequency-doubled of the radial force and axial force are obvious when parallel misalignment exists. With the increment of parallel misalignment, the amplitude increment of the fundamental frequency and frequency-doubled of the radial force is far bigger than axial force, that is, the radial vibration is the main. This is in accordance with the theoretical analysis shown in Table 3.
  • Angle Misalignment
    The local coordinate system (marker) rotates a certain angle around the y-axis, and the z-axis of the local coordinate system is set up as the orientation of revolution joint, and then the corresponding angle misalignment can be simulated. Some results of the angle misalignment are shown in Table 5.
    From Table 5, it can be found the fundamental frequency of radial force and axial force are obvious when angle misalignment exists. With the increment of angle misalignment, the fundamental frequency amplitude is also increasing, but priority is given to the axial vibration. This accords with the theoretical analysis shown in Table 3.
  • Comprehensive Misalignment
    Adding the parallel misalignment and angle misalignment in the local coordinate system (maker) of the left half coupling at the same time, comprehensive misalignment can be simulated. Some results of the comprehensive misalignment are shown in Table 6.
    From Table 6, it can be seen that the radial vibration and axial vibration are obvious at the same time, as well as the fundamental frequency and the frequency doubling when the comprehensive misalignment happens. The simulation results verify the theoretical conclusion very well.
    Through validation, it can be concluded the simulation results accord with theoretical conclusions. This is effectively guaranteeing the validity of the extraction of vibration signals.

3.2. Acquisition of the Vibration Signals of the Gearbox High-Speed Shaft

In this paper, the vibration signals (angular acceleration) of the gearbox high-speed shaft are obtained for analysis.
According to the requirements of the selected unit, the allowance range of the parallel misalignment is 0–0.15 mm. The allowance range of the angle misalignment is 0°–0.05°. In these ranges, the wind turbine system runs normally. In this paper, the range of parallel misalignment simulated is 1–10 mm. The range of angle misalignment that simulated is 0.10°–10°. The samples used are simulated by the method described above from four cases of normal working conditions, parallel misalignment, angle misalignment and the comprehensive misalignment. Each kind of sample is 68, with a total of 272 samples.

4. Extracting IEMD Energy Entropy from the Vibration Signals

Entropy is a measure of the uncertainty degree of the information, including energy entropy, singular entropy, and the approximate entropy (they are generally referred to as the information entropy). Energy entropy has a very good effect on the extraction of non-stationary and nonlinear complex signal statistical characteristics [46]. Its definition is as follows:
Assuming
E i = | c i ( t ) | 2 d t = k = 1 n | x i k 2 | ,   E = i = 1 n E i ,   p i = E i E ,
where i = 1 , 2 , 3 , , n , x i k is the amplitude of each discrete point.
Then, the energy entropy is the following:
P i = i = 1 n p i lg p i .
When misalignment occurs in DFWT, the frequency of the vibration signal will change, and so will the energy distribution. Therefore, the energy entropy can be used to reflect the characteristics of vibration signals. In this paper, vibration signals are analyzed by the IEMD first, and then the energy entropy of the IMF components are calculated. Different numbers of the IMF may be obtained from different vibration signals, and the correlation coefficient method is used to determine the number of IMF components to be put into the SVM. The correlation coefficient method is to determine the correlation coefficient between the IMF information and the original signal [47]. Assuming that a non-stationary signal x ( t ) can be decomposed into a finite mutually uncorrelated component x i ( t ) , that is:
x ( t ) = x 1 ( t ) + x 2 ( t ) + + x N ( t ) ,
the correlation coefficient of each component and the original signal is calculated as:
ρ x i x = t = 0 x i ( t ) x ( t ) [ t = 0 x i 2 ( t ) t = 0 x 2 ( t ) ] 2 ,   i = 1 , 2 , , N ,
and the larger the value of correlation coefficient, the greater the relevance, and vice versa.
For example, after processing the vibration signal, which is shown in Figure 1 by IEMD, the correlation coefficients of all IMF components are calculated as shown in Table 7.
The contribution rate of the first eight IMF components in 13 IMFs reaches 92.23%. Thus, it can be concluded that the first eight IMFs contain the main information of the signals. Through a lot of experiments, it is known that the correlation between the first several IMF components and the original signal are relatively big. The contribution rates of the first eight IMFs are over 90%. Therefore, they are used to calculate the energy entropy to characterize the different states in this paper. The partial data of the first eight IEMD energy entropies of the vibration signals are shown in Table 8.

5. Fault Diagnosis Based on IEMD Energy Entropy and PSO-SVM

After extracting the IEMD energy entropy of the vibration signals, the samples were divided into two equal groups, one was the training set, containing 45 × 4 samples, and the other was the testing set, containing 23 × 4 samples. The training samples were used to obtain the optimized SVM classifier (3-fold cross validation). The PSO algorithm was used to optimize the penalty parameter C and kernel parameter g of SVM, and the whole fault diagnosis processing can be referred to as IEMD-PSO-SVM; the specific implementation is shown in Figure 2.
Learning factors c 1 and c 2 control the degree of interaction between the individual experience and the population experience, which reflects the information exchange between the groups. The appropriate choice of c 1 and c 2 can not only speed up the search but also avoid the whole population plunging into the local optimum. Learning factors generally take a fixed constant value [48], according to experience. In this paper, let c 1 = c 2 = 2 . The classification results of DFWT using IEMD-PSO-SVM are shown in Figure 3 and Figure 4.
From Figure 3 and Figure 4, it can be seen that the testing accuracy of the SVM optimized by PSO is 95.5556%, and the optimal parameters are: C = 59.966 and g = 1.75228 . Only one sample is misclassified. The training accuracy at this time reaches 100% and the testing accuracy at this time also reaches 98.913%, with the correct rate of fault classification being very high. In order to illustrate the superiority of the proposed algorithm, the same energy entropy of the same vibration signals are classified by GridSearch-SVM (the parameters of SVM are optimized by GridSearch) and genetic algorithm (GA)-SVM (the parameters of SVM are optimized by Genetic Algorithm), and the recognition results are shown in Figure 5.
It can be seen that four samples are misclassified by the IEMD-GA-SVM model, while three samples are misclassified by IEMD-GridSearch-SVM. The testing accuracies are, respectively, 95.6522% and 96.7391%.
In a similar way, the energy entropy calculated after the same vibration signals being decomposed by EMD can be put into the PSO-SVM, GridSearch-SVM and GA-SVM, and the classification results are shown in Figure 6 and Table 9.
It can be seen that the promotion ability of EMD-GA-SVM, EMD-GridSearch-SVM and EMD-PSO-SVM was not high. Furthermore, the corresponding recognition accuracy decomposed by EMD was much lower than that by IEMD. The promotion ability of IEMD-GA-SVM and IEMD-GridSearch-SVM was less than that of IEMD-PSO-SVM. Thus, on the whole, the IEMD-PSO-SVM algorithm is better in diagnostic performance.

6. Conclusions

In this paper, aiming at the fact that misalignment fault samples are not enough to diagnosis, Solidworks and Adams are used to simulate the working conditions of DFWT to obtain the characteristic signals. After eliminating the end effects of EMD by using IEMD, the IEMD energy entropy of the vibration signals are extracted as the fault features. SVM, which is good at small sample learning, is used as the classifier to identify the misalignment. The parameters of SVM are optimized by the PSO algorithm. The results show that the fault diagnosis accuracy of misalignment of DFWT is higher. The effectiveness of the model is proved.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (51577008). The authors are grateful to the anonymous reviewers for their constructive comments on the previous version of the paper.

Author Contributions

Yancai Xiao and Na Kang contributed to paper writing and the whole revision process. Yi Hong and Guangjian Zhang helped collect some samples. All the authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, T. Statistics on global wind power installed capacity in 2015. Wind Energy Ind. 2016, 2016, 51–56. [Google Scholar]
  2. Bak-Jensen, B.; Kawady, T.A.; Abde1-Rahman, M.H. Coordination between fault-ride-through capability and over-current protection of DFIG generators for wind farms. J. Energy Power Eng. 2010, 4, 20–29. [Google Scholar]
  3. Liao, M.; Liang, Y.; Wang, S.; Wang, Y. Failure analysis of Misalignment of wind turbine failure. Mech. Sci. Technol. 2011, 30, 173–190. [Google Scholar]
  4. Wan, Q.; Geng, L.; Guo, T. The BP neural network based on fruit flies for gearbox fault diagnosis. Power Syst. Clean Energy 2014, 9, 006. [Google Scholar]
  5. Ding, S.; Chang, X.; Wu, Q. Wind turbine gearbox fault diagnosis based on LVQ neural network. J. Mod. Electron. Tech. 2014, 10, 150–152. (In Chinese) [Google Scholar]
  6. Li, H. Large Wind Turbine Gearbox Fault Fuzzy Diagnosis Technology Research. Ph.D. Thesis, Shenyang University of Technology, Shenyang, China, 2014. [Google Scholar]
  7. Li, Z.; Ma, Z.; Jiang, R.; Liu, Y. Fault classification method research for wind turbine gearbox. Mach. Des. Manuf. 2015, 2, 177–180. [Google Scholar]
  8. Zhang, Y.; Zhang, L. The construction of a wind turbine engine fault diagnosis expert system. J. Wind Power 2008, 12, 97–99. [Google Scholar]
  9. Liu, Y.; Wang, F.; Shi, W. Wind turbines operating conditions based on support vector machine (SVM) classification method. Acta Energiae Sol. Sin. 2010, 9, 1191–1197. [Google Scholar]
  10. Huang, Y. The wind turbines spindle bearing fault diagnosis based on support vector machine. J. Electron. Instrum. Cust. 2016, 11, 88–92. [Google Scholar]
  11. Precupa, R.E.; Angelovb, P.; Costad, B.S.J.; Sayed-Mouchawehe, M. An Overview on Fault Diagnosis and Nature-Inspired Optimal Control of Industrial Process Applications. Elsevier J. Comput. Ind. 2015, 74, 75–94. [Google Scholar] [CrossRef]
  12. Toubakh, H.; Mouchaweh, M.S. Hybrid Dynamic Classifier for Drift-Like Fault Diagnosis in a Class of Hybrid Dynamic Systems: Application to Wind Turbine Converters. Neurocomputing 2016, 171, 1496–1516. [Google Scholar] [CrossRef]
  13. Rosa, M.; Blesa, J.; Tornil-Sin, S. Fault detection and isolation for a wind turbine benchmark using a mixed Bayesian/Set-membership approach. Annu. Rev. Control 2015, 40, 59–69. [Google Scholar] [Green Version]
  14. Joaquim, B.; Pedro, J.; Damiano, R. An Interval NLPV Parity Equations Approach for Fault Detection and Isolation of a Wind Farm. IEEE Trans. Ind. Electron. 2015, 62, 3794–3805. [Google Scholar]
  15. Blesa, J.; Rotondo, D.; Puig, V.; Nejjari, F. FDI and FTC of wind turbines using the interval observer approach and virtual actuators/sensors. Control Eng. Pract. 2014, 24, 138–155. [Google Scholar] [CrossRef]
  16. Liu, Z.; Zhao, G.; Long, X. Summary of the research on the misalignment of rotor system coupling. Turbine Technol. 2007, 5, 321–325. [Google Scholar]
  17. Li, M. Research on Feature Extraction and Fault Diagnosis of Diesel Engine Vibration Signal. Ph.D. Thesis, Northwest A&F University, Yangling, China, 2012. [Google Scholar]
  18. Zhang, J.; Zhai, Y.; Wang, S. Fault diagnosis of wind turbine gearbox based on wavelet decomposition and least square support vector machine. Trans. Microsyst. Technol. 2011, 1, 41–43. [Google Scholar]
  19. Wang, X. On Line Fault Diagnosis of Wind Turbine Based on Wavelet Transform and Data Mining. Ph.D. Thesis, Lanzhou University of Technology, Lanzhou, China, 2010. [Google Scholar]
  20. Ju, P.; Qin, S.; Qin, Y. Research on multi resolution EMD method and frequency domain average in gear incipient fault diagnosis. J. Vib. Shock 2009, 28, 97–101. [Google Scholar]
  21. Zheng, X.; Xu, H.; Fu, Y. Research on early fault diagnosis of wind turbine bearing based on empirical mode decomposition. Power Gener. Technol. 2013, 2, 471–474. [Google Scholar]
  22. Lu, W.; Du, R. Mechanical Engineering Test, Information and Signal Analysis; Press of Central China University of Science & Technology: Wuhan, China, 1999. [Google Scholar]
  23. He, Z.; Sun, Z.; Zhang, X. Modern Signal Processing and Engineering Application; Xi’an Jiao Tong University Press: Xi’an, China, 2007. [Google Scholar]
  24. Wang, Z.; Li, W. Wigner-Ville distribution and its application in signal analysis. Sichuan Ordnance J. 2008, 29, 1517. [Google Scholar]
  25. Zou, H. Non existence of time frequency analysis without cross interference and WVD aggregation. Chin. J. Sci. E 2001, 31, 348–354. [Google Scholar]
  26. Cui, J. An Introduction to Wavelet Analysis M. Cheng Zhengxing, Translation; Xi’an Electronic and Science University Press: Xi’an, China, 1994. [Google Scholar]
  27. Zhang, R. Fault Feature Extraction Technique and System of Rotating Machinery. Ph.D. Thesis, Zhejiang University, Hangzhou, China, 2004. [Google Scholar]
  28. Zhao, Z. Research on Mechanical Fault Feature Extraction and Diagnosis Based on Vibration Signal. Ph.D. Thesis, Beijing Jiaotong University, Beijing, China, 2012. [Google Scholar]
  29. Zhu, L. Research on Vibration Fault Diagnosis of Induction Motor Based on Wavelet Neural Network. Ph.D. Thesis, Taiyuan University of Technology, Taiyuan, China, 2008. [Google Scholar]
  30. Li, S.; Guo, H.; Li, D. Summary of methods of vibration signal processing. J. Sci. Instrum. 2013, 8, 1907–1915. [Google Scholar]
  31. Yu, J.; Xie, S.; Zang, G.; Li, Q. Overview of signal processing technology in fault diagnosis of rotating machinery. Mach. Tool Hydraul. 2011, 12, 107–110. [Google Scholar]
  32. Cheng, J. Research on Fault Diagnosis Method of Rotating Machinery Based on Hilbert-Huang Transform D. Ph.D. Thesis, Hunan University, Changsha, China, 2005. [Google Scholar]
  33. Wang, H. HHT Method and Its Application Research D. Ph.D. Thesis, Hefei University of Technology, Hefei, China, 2009. [Google Scholar]
  34. Chen, J.; Cui, Y.; Wang, L.; Zhu, Q. Research on fault feature extraction method of gear box based on improved EMD. Mech. Transm. 2014, 5, 102–106. [Google Scholar]
  35. Huang, N.E.; Wu, M.C.; Long, S.R.; Shen, S.S.P.; Qu, W.; Gloersen, P.; Fan, K.L. A confidence limit for the empirical mode decomposition and the Hilbert spectral analysis. Proc. R. Soc. Lond. Ser. A 2003, 459, 2317–2345. [Google Scholar] [CrossRef]
  36. Zhao, J. Study on the effects of abnormal events to empirical mode decomposition method and the removal method for abnormal signal. J. Ocean Univ. Qingdao 2001, 31, 805–814. [Google Scholar]
  37. Zhang, Y.; Liang, J.; Hu, Y. Application of AR mode to improve end effect of EMD. Prog. Natl. Sci. 2003, 13, 1054–1059. [Google Scholar]
  38. Zeng, K.; He, M. A simple boundary process technique for empirical decomposition. In Proceedings of the IGARSS 2004, Anchorage, AK, USA, 20–24 September 2004; pp. 4258–4261.
  39. Gai, Q.; Ma, X.; Zhang, H. New method for processing end effect in local wave method. J. Dalian Univ. Technol. 2002, 41, 115–117. [Google Scholar]
  40. Ma, L.; Xiao, B.; Wang, C. Sparse Representation Based on K-Nearest Neighbor Classifier for Degraded Chinese Character Recognition; Springer: Berlin/Heidelberg, Germany, 2010; pp. 506–514. [Google Scholar]
  41. Wu, D. Gear box fault diagnosis method based on support vector machine. J. Vib. Meas. Diagn. 2008, 28, 339–342. [Google Scholar]
  42. Vladimir, V.N.; Zhang, X. The Essence of Statistical Learning Theory; Tsinghua University Press: Beijing, China, 1999. [Google Scholar]
  43. Azriel, R.; Wechser, H. Pattern recognition: Historical perspective and future directions. Int. J. Imaging Syst. Technol. 2000, 11, 101–116. [Google Scholar]
  44. Chen, Y.; Ma, H. Research on fault diagnosis of shearer rocker arm gear box based on PSO-SVM. J. Coal Mine Mach. 2015, 36, 303–306. [Google Scholar]
  45. Liao, M.; Liang, Y.; Wang, S.; Wang, Y. Misalignment in Drive Train of Wind Turbines. Mech. Sci. Technol. 2011, 30, 173–180. [Google Scholar]
  46. Zhang, C.; Chen, J.; Guo, X. Gear fault diagnosis method based on EEMD energy entropy and support vector machine. J. Cent. South Univ. 2012, 3, 932–939. [Google Scholar]
  47. Su, W.; Wang, F.; Zhang, Z.; Guo, Z.; Li, H. EMD denoising and spectral kurtosis in the application of fault diagnosis of rolling bearing. J. Vib. Shock 2010, 23, 34–46. [Google Scholar]
  48. Ratnaweeraa, H.; Watson, H.C. Self-organizing Hierarchical Particle Swarm Optimizer with Time-Varying Acceleration Coefficients. IEEE Trans. Evol. Comput. 2004, 8, 240–255. [Google Scholar] [CrossRef]
Figure 1. The collected vibration signals.
Figure 1. The collected vibration signals.
Entropy 19 00006 g001
Figure 2. The IEMD-PSO-SVM fault diagnosis mode.
Figure 2. The IEMD-PSO-SVM fault diagnosis mode.
Entropy 19 00006 g002
Figure 3. Particle fitness curves.
Figure 3. Particle fitness curves.
Entropy 19 00006 g003
Figure 4. IEMD-PSO-SVM testing results.
Figure 4. IEMD-PSO-SVM testing results.
Entropy 19 00006 g004
Figure 5. Testing results of IEMD energy entropy. (a) IEMD-GA-SVM testing results; (b) IEMD-GridSearch-SVM testing results.
Figure 5. Testing results of IEMD energy entropy. (a) IEMD-GA-SVM testing results; (b) IEMD-GridSearch-SVM testing results.
Entropy 19 00006 g005
Figure 6. Testing results of EMD energy entropy. (a) EMD-GA-SVM testing results; (b) EMD-GridSearch-SVM testing results; (c) EMD-PSO-SVM testing results.
Figure 6. Testing results of EMD energy entropy. (a) EMD-GA-SVM testing results; (b) EMD-GridSearch-SVM testing results; (c) EMD-PSO-SVM testing results.
Entropy 19 00006 g006
Table 1. Comparison of various modern signal processing methods.
Table 1. Comparison of various modern signal processing methods.
MethodsAdvantagesShortcomings
Short Time Fourier TransformEasy to implement; avoiding the cross-term interference in the high-order non-stationary analysis [22].High redundancy; limited time-frequency aggregation; lacking of self adaptability [23].
Wiger-VilleHigh time-frequency resolution; concentration of energy; having time-frequency edge characteristics [24].Energy density is not guaranteed to be non-negative in all time and frequency ranges; a serious cross-term interference [25].
Wavelet TransformLocalized in time domain; multi resolution; adjustable analysis precision [26].Leaking of energy; lacking of self adaptability [27].
Wavelet Packet TransformDecomposing the low frequency portion while carrying out a more detailed decomposition to high frequency portion of the signal; no redundancy and no omissions [28].Wavelet packet energy and time-frequency information are not intuitive [29].
Blind Source SeparationIdentifying the original signal directly according to the observation information, solving the problem that the traditional methods can not solve [30,31].Difficulty in identifying the unknown or dynamic change signal source [30,31].
High-order Statistics AnalysisAbility of automatically suppress the influence of Gaussian colored noise and some non-Gaussian colored noise; high-order cyclic statistics can automatically suppress any stationary (Gaussian and non-Gaussian) noise [30].Complicated and no clear physical meaning [30].
EMDSuitable for nonlinear and non-stationary signal analysis; high signal-to-noise ratio; self adaptation; orthogonality and completeness; high efficiency [32,33].End effect problem; mode mixing problem and the limitations of decomposition [34].
Table 2. Comparison of theoretical results and simulation results at a normal working state.
Table 2. Comparison of theoretical results and simulation results at a normal working state.
ItemsSimulation Rotational Speed (°/s)Theoretical Rotational Speed (°/s)Relative ErrorSimulation Fundamental Frequency (Hz)Theoretical Fundamental Frequency (Hz)Relative Error
Low-speed solar gear520.97517.460.68%26.2226.390.64%
Medium-speed solar gear2911.332918.470.24%166.54166.730.11%
High-speed pinion10528.6110527.340.01%818.28818.790.06%
Shaft of generator10520.6810527.340.06%818.89818.790.01%
Table 3. The characteristics of misalignment.
Table 3. The characteristics of misalignment.
Fault TypeSpectral Characteristics
Parallel misalignmentThe fundamental frequency and frequency-doubled of power frequency of rotor system are relatively obvious; the radial vibration is more obvious than the axial vibration; the amplitude increases with the increment of the degree of the misalignment.
Angle misalignmentAngle misalignment produces large axial vibration, the vibration frequency is given priority to fundamental frequency; vibration amplitude increases with the increment of the degree of the misalignment.
Comprehensive misalignmentContains the fault characteristics of parallel misalignment and angle misalignment.
Table 4. Some simulation results of the parallel misalignment.
Table 4. Some simulation results of the parallel misalignment.
The Degree of Parallel Misalignment (mm)Radial Force Times Frequency Amplitude (N)Axial Force Times Frequency Amplitude (N)
Fundermental Frequence AmplitudeFrequency-Doubled AmplitudeFundermental Frequence AmplitudeFrequency-Doubled Amplitude
Normal condition66576871798271
1858620582416265
220,93616,8416364427
336,51518,41111,114405
Table 5. Some simulation results of the angle misalignment.
Table 5. Some simulation results of the angle misalignment.
The Degree of Angle Misalignment (°)Radial Force Times Frequency Amplitude (N)Rxial Force Times Frequency Amplitude (N)
Fundermental Frequence AmplitudeFrequency-Doubled AmplitudeFundermental Frequence AmplitudeFrequency-Doubled Amplitude
Normal condition66576871798271
0.168647163789390
0.274059414314399
0.311,0487644408384
Table 6. Some simulation results of the comprehensive misalignment.
Table 6. Some simulation results of the comprehensive misalignment.
The Degree of Comprehensive MisalignmentRadial Force Times Frequency Amplitude (N)Axial Force Times Frequency Amplitude (N)
Fundermental Frequence AmplitudeFrequency-Doubled AmplitudeFundermental Frequence AmplitudeFrequency-Doubled Amplitude
Normal condition66576871798271
0.1°--1 mm551948051577888
0.2°--2 mm28,29722,46522801114
0.3°--3 mm39,16419,2144834887
Table 7. Correlation coefficient of each IMF component and the original signal.
Table 7. Correlation coefficient of each IMF component and the original signal.
ObjectIMF1IMF2IMF3IMF4IMF5IMF6IMF7
Correlation coefficient0.64030.54400.33910.23740.14030.09030.0913
ObjectIMF8IMF9IMF10IMF11IMF12IMF13
Correlation coefficient0.10860.04480.03320.04900.05790.0006
Table 8. The partial data of the first eight IEMD energy entropies of vibration signals in different working conditions.
Table 8. The partial data of the first eight IEMD energy entropies of vibration signals in different working conditions.
Fault TypeP1P2P3P4P5P6P7P8
Normal condition0.31090.34290.17090.07470.04760.0720.0970.0521
0.31660.3520.15530.07840.04770.06880.05750.0176
0.31860.34980.16570.08140.04780.04520.0690.0172
0.32560.3410.20310.08950.05640.0370.09330.0148
0.33020.33860.16690.07480.05740.05220.07880.0129
Parallel misalignment0.14840.36760.36230.1930.12010.08630.04770.1012
0.14890.36690.36780.22670.09380.04810.0280.1249
0.15130.36780.36620.18390.10770.06650.05170.2216
0.15580.35330.36750.21330.14790.06830.06370.1977
0.15740.36010.36390.1570.08070.0620.20310.1429
Angle misalignment0.23350.30670.23250.14670.08430.07630.13860.3648
0.23610.30240.36790.2890.19340.11570.15940.0478
0.24230.36650.34170.25740.18620.08970.03950.0279
0.24550.34760.36790.30960.10830.08730.03370.0536
0.24870.32980.31340.25650.12540.08040.13090.3447
Comprehensive misalignment0.14560.22850.29120.33220.25210.16720.20010.1162
0.10610.27840.33930.28760.27620.08130.20590.1677
0.10740.25810.33570.35820.16060.11010.13960.1396
0.10820.28530.32710.30220.16850.17650.07510.1556
0.1110.32680.36540.27370.26810.08860.07730.1014
Table 9. Comparison of IEMD-PSO-SVM with other commonly used classifiers.
Table 9. Comparison of IEMD-PSO-SVM with other commonly used classifiers.
Item C g Elapsed Time (s)Accuracy of Training SetAccuracy of Testing SetFalse Alarm Rate
EMD-GA-SVM43.94312.138642.026898.8889% (178/180)69.5652% (64/92)9.7826% (9/92)
EMD-GridSearch-SVM11.31370.514.879086.1111% (155/180)68.4783% (63/92)9.7826% (9/92)
EMD-PSO-SVM68.71241.6815339.235398.8889% (178/180)71.7391% (66/92)8.6956% (8/92)
IEMD-GA-SVM29.81450.9865646.3340100% (180/180)95.6522% (88/92)1.0869% (1/92)
IEMD-GridSearch-SVM16215.3525100% (180/180)96.7391% (89/92)0
IEMD-PSO-SVM59.9661.7522825.6413100% (180/180)98.913% (91/92)0

Share and Cite

MDPI and ACS Style

Xiao, Y.; Kang, N.; Hong, Y.; Zhang, G. Misalignment Fault Diagnosis of DFWT Based on IEMD Energy Entropy and PSO-SVM. Entropy 2017, 19, 6. https://doi.org/10.3390/e19010006

AMA Style

Xiao Y, Kang N, Hong Y, Zhang G. Misalignment Fault Diagnosis of DFWT Based on IEMD Energy Entropy and PSO-SVM. Entropy. 2017; 19(1):6. https://doi.org/10.3390/e19010006

Chicago/Turabian Style

Xiao, Yancai, Na Kang, Yi Hong, and Guangjian Zhang. 2017. "Misalignment Fault Diagnosis of DFWT Based on IEMD Energy Entropy and PSO-SVM" Entropy 19, no. 1: 6. https://doi.org/10.3390/e19010006

APA Style

Xiao, Y., Kang, N., Hong, Y., & Zhang, G. (2017). Misalignment Fault Diagnosis of DFWT Based on IEMD Energy Entropy and PSO-SVM. Entropy, 19(1), 6. https://doi.org/10.3390/e19010006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop