To accurately extract the fault characteristics of vibration signals of rotating machinery is of great significance to the unit online monitoring and evaluation. However, because the current feature extraction methods are mainly for single channel, the results of feature extraction are often inaccurate. To this end, a coupling fault feature extraction method based on bivariate empirical mode decomposition (BEMD) and full spectrum is proposed for rotating machinery. Firstly, the two-dimensional orthogonal signal obtained by orthogonal sampling technique is decomposed by bivariate empirical mode decomposition to obtain the intrinsic mode function with phase information. In order to obtain the sensitive modal components, the sensitivity coefficients are constructed on the basis of mutual information. Then, the sensitivity coefficient of each intrinsic mode function is calculated, and the intrinsic mode function with the larger sensitive coefficient is selected as the sensitive component. Finally, the full spectrum of the sensitive component is obtained using the full vector envelope technique, so as to get a comprehensive and accurate characteristic component. The results of simulations experiment and an application example show that this method can extract the fault characteristic component of the rotating machinery comprehensively and accurately. It is of great significance to realize the accurate diagnosis of coupling faults of rotating machinery.
National Natural Science Foundation of China5177920651279161Water Resources Department of Shaanxi Province2015slkj-041. Introduction
Rotating machinery is one of the key parts of large mechanical equipment, such as hydraulic turbine, steam turbine, and wind turbine, which have found wide application in various industrial fields [1]. Due to the influence of load, damping, friction, and other factors, the rotating machinery usually exhibits complex dynamic behavior, which is related to nonstationarity and nonlinearity [2]. And the vibration signals of rotating machinery often show coupling. In order to ensure the safe and reliable operation of the equipment, the equipment condition monitoring and fault diagnosis technology needs to have higher precision. Whether it can accurately extract the characteristic signal of coupling fault is the key to the condition monitoring and fault diagnosis of rotating machinery [3, 4].
At present, many scholars have done some research on the feature extraction of rotating machinery fault signal, especially in the matter of signal processing methods. The main methods are wavelet transform (WT) [5, 6], Hilbert-Huang transform, empirical mode decomposition (EMD) [7, 8], variational mode decomposition (VMD) [9, 10], and so on, and the characteristics of different fault feature extraction methods are shown in Table 1. Wavelet transform is the improvement of the Fourier transform and has a certain ability to deal with nonlinear signals. However, its essence is the inner product principle characteristic waveform signal decomposition based on basis function, but there exists the problem of basis functions selection, and the different vibration will show the characteristics of different waveforms [11, 12]. Therefore, the adaptive ability of wavelet transform is poor. EMD has strong adaptability and good local analysis ability and has been widely used in the field of fault diagnosis. However, the EMD method lacks rigorous theoretical basis, and there is a serious phenomenon of modal aliasing in the decomposition process [13, 14]. Variational mode decomposition (VMD) is a new adaptive signal decomposition method, which has a solid theoretical foundation and can solve the problem of mode mixing in EMD, but there are parameter optimization problems in VMD [15].
The characteristics of different fault feature extraction method.
Method
Characteristic
WT
The choice problem for basis functions exists and the different vibration will show the characteristics of different waveforms.
EMD
It lacks rigorous theoretical basis and there is a serious phenomenon of modal aliasing in the decomposition process.
VMD
It solves the problem of mode mixing in EMD, but there are parameter optimization problems in VMD.
BEMD
It can detect synchronous features contained in the two-dimensional signal and comprehensively extract the information of fault signal.
However, due to the complex mechanism of the fault of rotating machinery, as well as the vibration signal being often nonlinear, the vibration signals in different directions may represent different fault characteristic information [16]. Through single-channel signal feature extraction and fault diagnosis, it often leads to misjudgment and omission [17, 18]. Therefore, it is more accurate to collect the signal of multiple channels, through the effective information fusion. Therefore, taking into account the rotating machinery, which generally has two sensors in horizontal and vertical directions, this paper utilizes bivariate empirical mode decomposition (BEMD) to decompose the signal collected from two directions.
In 2007, bivariate empirical mode decomposition (BEMD) was firstly introduced by Rilling et al. [19]. It is not limited to single direction processing of real valued signals, and it can analyze complex signals consisting of two orthogonal directions containing phase information and synchronization [20, 21].
Molla et al. [22] applied bivariate empirical mode decomposition to decompose complex climate signals into simple components and then combined BEMD with the classical FFT spectrum to clearly observe the climate change. Then Yang et al. [23] used bivariate empirical mode decomposition on condition monitoring for wind turbine; this method can overcome the problem of inability to carry out an information fusion function in EMD, which is of importance to reach a reliable condition monitoring conclusion.
The joint information between different sensors such as phase information and synchronization is of significance for rotating machinery condition to evaluate [24]. The signal obtained from multiple sensors is evaluated to contribute to the effective extraction of fault information, so it has received increasing attention in the past years. Okatan et al. [25] use the normalized matrix of the spectral norm to propose the new Kalman filter and propose a sensor information fusion Kalman filter test method, which improves the fusion accuracy. Liu et al. [26] proposed a fuzzy fault measurement method and a fuzzy integral information fusion method for rotating machinery. By using multiple classifiers, fuzzy measures, and fuzzy integral theory, the diagnosis results are fused and the final diagnosis results are obtained. The full spectrum is a fault diagnosis method based on homologous information fusion technology, which can acquire more comprehensive fault characteristics and reveal the direction of the vibration relative to the direction of the vibration of rotating machinery [24, 27].
In this paper, we proposed a coupling fault feature extraction method based on bivariate empirical mode decomposition and full spectrum for rotating machinery. In order to suppress noise interference in the signal component, the sensitivity coefficients are constructed on the basis of mutual information and through the sensitive coefficient to extract unique IMFs of the signal. The specific arrangement of this paper is organized as follows. A review on bivariate empirical mode decomposition (BEMD) is illustrated in Section 2. Section 3 describes the sensitive feature extraction based on mutual information. Section 4 gives brief introductions of full spectrum. Section 5 describes the coupling fault feature extraction method based on BEMD and full spectrum. Section 6 applies the proposed method to both simulated data and experimental data obtained from the turbine guide bearing. Conclusions come in Section 7.
2. Bivariate Empirical Mode Decomposition
The essence of bivariate empirical mode decomposition (BEMD) is to decompose a two-dimensional signal adaptively into intrinsic mode components with physical meaning. The modal component obtained from the decomposition is a series of single component signals from high frequency to low frequency in two directions. This paper takes Algorithm II in [19] to perform BEMD. For a two-dimensional signal x(t), the basic decomposition process is as follows.
Step 1.
Determine the direction of projection φk=2kπ/N, where 1≤k≤N.
Step 2.
The two-dimensional signal x(t) is projected onto the ϕk.(1)pφkt=Ree-φkxt.
Step 3.
Extract the corresponding moment for the local maximum of Xtjk; then the set tjk,eiφkpjk is interpolated. Get the maximum envelope eφk′(t) in the direction φk.
Step 4.
Calculate the mean of the maximum envelope eφk′(t) in each direction.(2)mt=2N∑k=1Neφk′t.
Step 5.
Similar to the EMD decomposition process, the residual component is calculated:(3)St=xt-mt.
Determine whether the S(t) meets the requirements of IMF. If satisfied, proceed to Step 6, and if not, repeat Steps 2–6 until S(t) satisfies the conditions of the intrinsic mode function IMF.
Step 6.
Record the resulting IMF and remove it from the original signal. And get IMF1 as c1(t)=h(t), residual component as(4)r1t=xt-c1t.
Step 7.
Repeat the above steps until you get all the IMFs. The original signal can be expressed as(5)xt=∑k=1Kckt+rKt,where K represents the total number of IMFs.
3. Sensitive Feature Extraction Based on Mutual Information3.1. Mutual Information
Mutual information can measure the degree of interdependence between the two variables, which means that the information content is shared between the two variables [28, 29]. Given two random variables X and Y, calculate their respective marginal probability distribution and joint probability distribution px,p(y), and p(x,y). The mutual information between them is(6)MIx,y=∑x∑ypx,ylogpx,ypxpy.
3.2. Sensitivity Coefficient Based on Mutual Information
The IMFs can be obtained by bivariate empirical mode decomposition (BEMD) for the fault vibration signals. However, not all the IMFs are important to the evaluation condition of machinery rotating, so it is significant to select the sensitive IMFs [30]. Firstly, orthogonal sampling technique is used to collect the vibration signals in the normal and fault state, which are recorded as xz(t) and x(t). Then the fault vibration signal x(t) is decomposed by BEMD. Get modal component IMFs, c1(t),…,cn(t).
(1) According to formula (6), the mutual information MIi between the intrinsic mode components IMF ci(t)(i=1,2,…,n) and the fault vibration signals x(t) is calculated. And do the normalization as follows:(7)ai=MIimaxMIi.
(2) Calculate the mutual information MIi′ between the intrinsic mode components IMF ci(t)(i=1,2,…,n) and the normal vibration signals xz(t). The same normalization is done as follows:(8)bi=MIi′maxMIi′.
(3) Calculate the sensitivity coefficient of each modal component IMFs as follows:(9)λi=aibi,where ai represents the mutual information between the nth IMFs and their original signals and bi represents the mutual information between the nth IMFs and the normal signals. In order to make the component better reflect the characteristics of the original signal, ai should be as large as possible. In order to better monitor the characteristics of the fault signals, bi should be as small as possible. Considering the characteristics of the early fault signal of rotating machinery, the amplitude of the signal is smaller, and the mutual information between the IMFs and the original signal is relatively small. Therefore, this paper reflects the change of signal characteristic components by increasing the multiple relations. That is to say, the higher the λi value is, the more sensitive the IMF is.
4. Full Spectrum Theory
There are usually two vertical sensors in the same section of the rotor for rotating machinery to extract vibration information. The full spectrum can fuse the information of vibration signals in two orthogonal directions and comprehensively express the intensity and spectrum structure of the rotor vibration [31, 32].
The core idea of the full spectrum is that rotating machinery shows whirling phenomena under the combination of harmonic frequencies. The whirling trajectory is a series of ellipses. The length of the long axis of the ellipse is defined as the principal vibration vector to evaluate vibration intensity. The length of the short axis of the ellipse is defined as the auxiliary vector to evaluate vibration intensity as auxiliary. And through the whirling intensity under the harmonic frequency of rotor, the rotating machinery fault is diagnosed and identified [33].
Set X, Y direction of the data sequence as x~ and y~. They can constitute the complex sequence z~=x~+iy~. Carry on the Fourier transform for the complex sequence z~, and get the frequency domain complex signal Z~.(10)Z~=X~+iY~=ReZ+ImZ.
The length of the long axis of the ellipse Rai is defined as the main vibration vector. (11)Rai=12NZi+ZN-i.
The length of the short axis of the ellipse Rbi is defined as vice vibration vector. (12)Rbi=12NZi-ZN-i
and αi is the angle between the main vibration vector and the X axis. (13)tan2αi=ImZiReZN-i-ReZiImZN-iReZiImZN-i-ImZiReZN-i.
ϕi is the initial phase angle of the elliptic trajectory at this frequency, whose calculation formula can be obtained according to the property of Fourier transform (the details of derivation were shown in [24]):(14)tanϕi=ImZi+ImZN-iReZi+ReZN-i,where i=1,2,…,N/2-1.
Therefore, utilize Fourier transform for the vibration signals in two orthogonal directions and obtain the characteristic information of each direction, whose full spectrum is required, which can simplify the calculation between complex Fourier parameters. This method not only greatly reduces the amount of calculation, but also is very accurate and established a contact with the conventional analysis method. When the information source is a single source, the algorithm is still established, fully meeting the requirements of real-time detection and analysis.
5. Coupling Fault Feature Extraction Method Based on BEMD and Full Spectrum
Owing to the complexity, coupling, and uncertainty of rotating machinery faults, the intrinsic dynamic characteristics of the faults are more complicated. The external manifestation is that the vibration signals in different directions may indicate different characteristic information, and the single-channel signal characteristics diagnosis for rotary machinery is prone to misjudgment and leakage judgment. And BEMD and full spectrum have advantages in the processing of two-dimensional signals and information fusion. Therefore, this paper proposes a new coupling fault feature extraction method based on BEMD and full spectrum for rotating machinery. To select the sensitive IMF, use the sensitivity coefficient based on mutual information as the screening criteria. The extraction process of this method is shown in Figure 1. Specific steps are as follows.
Flow chart of the proposed method.
Step 1.
Collect the orthogonal vibration signal of the rotating machine in the normal and fault condition by the sensor.
Step 2.
The orthogonal fault vibration signal is decomposed by BEMD to obtain K intrinsic mode components ckk=1,2,…,K.
Step 3.
Calculate the sensitivity of each IMF.
Step 4.
Select the modal component IMFs with the larger sensitive coefficient as sensitive modal components IMFs.
Step 5.
Full spectrum analysis of the selected modal components IMFs is carried out to detect the characteristic frequency of coupling fault signals.
6. Experiment Results6.1. Simulation Experiment
There are two common phenomena, a change in amplitude at the typical frequency [34] and a new frequency appearance [35], when rotating machinery breaks down. In this paper, the vibration signal of the hydropower unit bearing is simulated. Suppose that the rotating frequency of hydropower unit is f0, It is considered that the vibration fault signal of the actual hydroelectric generating set often contains f0, 2f0, 3f0, 50 Hz, 100 Hz, and some common characteristic signals and noise. In order to make the simulation result more fitting with the actual situation and verify the effectiveness of the proposed method, this paper, in the processing of two-dimensional signals, simulates the normal and fault signals of the hydropower unit. The specific simulation vibration signal is as follows:(15)xz=8sinωt+5sin5ωt+3sin100πt+2rand1,N,yz=8cosωt+5cos5ωt+3cos100πt+2rand1,N,x=8sinωt+5sin5ωt+5sin100πt+2rand1,N+5sin2ωt,y=8cosωt+5cos5ωt+5cos100πt+2rand1,N+5cos2ωt,where the xz and yz are used to simulate the normal signals, which are composed of a rotating component at a frequency of f0 with amplitude of 8, a rotating component at a frequency of 5f0 with amplitude of 5, a rotating component at a frequency of 50 Hz with amplitude of 3, and the noise generated randomly with amplitude of 2. Each component of the normal vibration signal is shown in Figure 2(a). Figure 2(a) separately represents the rotating component at a frequency of f0, 5f0, 50 Hz, and the noise generated randomly in turn.
Components of normal and fault vibration signal.
Components of normal vibration signal
Components of fault vibration signal
The x and y are used to simulate the fault signals, which not only are composed components of the normal signal, but also have the rotating component at a frequency of 2f0 with amplitude of 5 and the rotating component at a frequency of 50 Hz with increased amplitude of 2. Each component of the fault vibration signal is shown in Figure 2(b). And Figure 2(b) separately represents the rotating component at a frequency of f0, 5f0, 50 Hz, the noise generated randomly, and the rotating component at a frequency of 2f0 in turn, where ω represents the angular velocity, in actual hydroelectric generating set, the rated speed is 107.1 r/min, the sampling frequency is 227 Hz, and the sampling point is 1024. Considering the consistency with the actual operation of hydropower units and the simplicity of simulation, it can be assumed that the unit speed is 107 r/min, sampling frequency is 250 Hz, and sampling point N is 1000. The normal and the fault raw signal are shown in Figure 3, which separately represent the normal signal of the X direction, the normal signal of the Y direction, the fault signal of the X direction, and the fault signal of the Y direction in turn.
Observation signal of normal and faulty.
The proposed method is used to extract the feature of the simulation signal in Figure 3. Firstly, the orthogonal fault raw signal is decomposed by BEMD, and the results are shown in Figure 4. Where the blue curve represents the real part and the black one represents the imaginary part of the decomposition results by BEMD. Then calculate the sensitivity coefficient of each IMF, as shown in Table 2.
Sensitivity coefficient of each characteristic component.
IMF
IMF1
IMF2
IMF3
IMF4
IMF5
Sensitivity coefficient
2.881
0.693
1.078
0.767
0.654
Results of BEMD decomposition for fault signal.
The IMF1 and IMF3, whose sensitivity coefficient is greater than 1, are selected as sensitive IMFs. Utilize the full vector envelope technique to fuse the sensitive IMFs selected and obtain full spectrum of the sensitive IMFs as shown in Figure 5.
Full spectrum of sensitive IMFs.
Based on the full spectrum of the sensitive IMFs in Figure 5, the method proposed in the paper can fuse the orthogonal signals in two directions. And it can accurately detect the new frequency appearance at a frequency of 5f0 and the change in amplitude at the typical frequency of 50 Hz. Obviously, the method has strong compatibility and effectiveness.
6.2. Coupling Fault Feature Extraction of Turbine Guide Bearing
To verify the effectiveness of the proposed method in coupling fault feature extraction for rotating machinery, the experimental studies on a hydroelectric turbine in upper reaches of Yellow River are conducted. The experimental signals are acquired from the prototype of hydroelectric turbine with 5 blades and 16 guide vanes, the maximum water head is 25.7 m, the rated head is 16 m, the rated power of the turbine is about 49 MW, and the rated speed is 107.1 r/min (1.79 Hz). Figure 6 describes specific layout of measuring points of pressure fluctuation signals in the turbine [36].
Specific layout of measuring point in hydroelectric turbine.
The measured data are collected turbine guide bearing by the vibration and throw monitoring sensors turbine guide bearing (As shown in Figure 6), which are from horizontal and vertical direction. Randomly select a segment of signal from normal and fault condition to verify the effectiveness of the proposed method, as shown in Figure 7: from top to bottom, the normal signal of the X direction, the normal signal of the Y direction, the fault signal of the X direction, and the fault signal of the Y direction, where the sampling frequency is 227 Hz and the sampling point M is 1024. Under the same load, the fault raw data is collected when the unit is abnormal.
Raw vibration data of normal and fault signals.
Using the method proposed in this paper to extract the coupling fault feature of turbine guide bearing, firstly, the fault signal is decomposed by BEMD, and the results of the decomposition are shown in Figure 8 (where the blue represents the real part and the black represents the imaginary part). As shown in Figure 8, the BEMD decomposition can effectively preserve the phase information of the orthogonal two-dimensional signal, and the decomposition results have good synchronization. In order to illustrate the advantages of BEMD, the real and imaginary parts of the fault signal are separately decomposed by EMD, and the results are shown in Figure 9. Obviously, when utilizing EMD to decompose the orthogonal two-dimensional signal, the phase information contained in the orthogonal signal is lost. In addition, compared to the decomposition results of BEMD, the decomposition component of EMD is too much, that is, the phenomenon of overdecomposition.
BEMD decomposition results of measured fault signal.
EMD decomposition results of real and imaginary part of measured fault signal.
The decomposition results of real fault signal
The decomposition results of imaginary part fault signal
On the basis of BEMD decomposition, the sensitivity coefficients of each IMFs are calculated, as shown in Table 3. And select the sensitive coefficient of large IMF2, IMF3, and IMF5 as the sensitive component IMFs. Then utilize the full vector envelope technique to fuse the sensitive IMFs selected, and obtain the full spectrum of the sensitive IMFs as shown in Figure 10. Obtain the decomposition results of EMD to compare them with the decomposition results of BEMD as shown in Figure 11. Calculate the sensitivity coefficients of the IMFs obtained by EMD in two directions, as shown in Table 4. In the IMFs of the real part, IMF2, IMF3, and IMF6 with the larger sensitivity coefficient are selected as the sensitive IMFs, and In the IMFs of the imaginary part, IMF2, IMF4, and IMF5 with the larger sensitivity coefficient are selected as the sensitive IMFs. The frequency spectrum analysis of the sensitivity coefficient in each direction is carried out. The frequency spectrum of the EMD sensitive IMFs is shown in Figure 11 (where the upper part represents the real part and the lower part represents the imaginary part of the decomposition results by BEMD).
Sensitivity coefficients of BEMD modal components.
IMFs
IMF1
IMF2
IMF3
IMF4
IMF5
Sensitivity coefficient
0.661
1.043
1.172
0.564
0.984
Sensitivity coefficients of EMD modal components.
IMFs
IMF1
IMF2
IMF3
IMF4
IMF5
IMF6
IMF7
IMF8
IMF9
Real part
0.831
1.634
1.875
0.73
0.455
1.344
0.495
0.393
0.543
Imaginary part
0.76
2.292
0.877
1.12
1.853
0.742
0.436
0.874
0.432
Full vector envelope spectrum of BEMD.
Spectrum of EMD decomposition sensitive component.
It can be seen from Figure 11 that the vibration signal (real part) in the horizontal direction can be divided into 2f0, a small number of the 3-frequency 3f0 and 15 f0. However, the vibration signal in the vertical direction (imaginary part) contains 1f0, 15 f0, and 50 Hz vibration signal. Therefore, the extraction results for vibration signals in the two orthogonal directions are different due to the certain interference, which often brings trouble for condition monitoring of rotating machinery. Compared with this, the method proposed in this paper can solve the problem of the inconsistency between the two channels. As shown in Figure 10, this method can better fuse the information in two orthogonal directions of vibration signal. This method has good compatibility, which can be more comprehensive and accurate to detect the vibration signal in the rotating machinery, and the feature extraction results are more reliable.
7. Conclusions
This paper proposed a method to extract a feature from orthogonal signals in two-direction sensors for the condition monitoring of rotating machinery based on bivariate empirical mode decomposition (BEMD) and full spectrum.
BEMD is employed to decompose signals from two orthogonal sources together; thus a complicated rotation can be represented by a set of simpler rotation components. BEMD is proved to outperform standard EMD for two orthogonal signals, because the intrinsic mode functions (IMFs) derived by BEMD preserve the signal phase information, detect synchronous features contained in the two-dimensional signal, and effectively solve the problem of misjudgment and leakage judgment in extract feature method by single-channel signal for rotary machinery. Thus the comparison of IMFs for different health conditions is made reasonable and easy.
A criterion based on mutual information is proposed for selecting the most sensitive IMF. The definition of this standard is to better represent the IMF of the original signal and penalize IMFs that cannot distinguish between different health conditions. Therefore, the IMF selected always ensures that it retains unique information about a particular health condition.
The performance of the proposed method has been evaluated for the hydroelectric turbine in upper reaches of Yellow River. Experiments results verified the effectiveness of the proposed method in two orthogonal directions, fault feature extraction and fault diagnosis for rotating machinery.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (nos. 51779206, 51279161) and the science and technology projects of Water Resources Department of Shaanxi Province (no. 2015slkj-04). Furthermore, the authors are grateful to the staff of China Electric Power Research Institute and State Grid Gansu Province Electric Power Research Institute.
HeX.-H.WangD.LiY.-F.ZhouC.-H.A novel bearing fault diagnosis method based on gaussian restricted boltzmann machine201620168295708310.1155/2016/2957083LiX.ZhengA.ZhangX.LiC.ZhangL.Rolling element bearing fault detection using support vector machine with improved ant colony optimization20134682726273410.1016/j.measurement.2013.04.0812-s2.0-84879089831LangG.LinJ.LiaoY.ZhaoM.Impact of system anisotropy on vibration reduction of rotating machinery and its evaluation method2017932993112-s2.0-8501532702410.1016/j.ymssp.2017.02.017LiangL.LiuF.LiM.XuG.Manifold learning with self-organizing mapping for feature extraction of nonlinear faults in rotating machinery201520151187390510.1155/2015/8739052-s2.0-84945299276CohenA.2003ElsevierFengY.CaoJ.-R.LuB.-C.ZhangD.-F.Fault feature extraction method for rolling bearing based on wavelet transform optimized by continuous kurtosis2015341427322-s2.0-8493920568610.13465/j.cnki.jvs.2015.14.006HuangN. E.ShenZ.LongS. R.WuM. C.ShihH. H.ZhengQ.YenN.TungC. C.LiuH. H.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis199845490399510.1098/rspa.1998.0193MR1631591Zbl0945.620932-s2.0-5444236478JunshengC.DejieY.YuY.A fault diagnosis approach for roller bearings based on EMD method and AR model20062023503622-s2.0-2774458064010.1016/j.ymssp.2004.11.002DragomiretskiyK.ZossoD.Variational mode decomposition2014623531544MR31602932-s2.0-8489341875610.1109/TSP.2013.2288675LiuS.TangG.WangX.HeY.Time-frequency analysis based on improved variational mode decomposition and teager energy operator for rotor system fault diagnosis201620169171304610.1155/2016/17130462-s2.0-85003874725KeskesH.BrahamA.LachiriZ.Broken rotor bar diagnosis in induction machines through stationary wavelet packet transform and multiclass wavelet SVM2013971511572-s2.0-8487385718110.1016/j.epsr.2012.12.013PercivalD. B.WaldenA. T.2006Cambridge, UKCambridge University PressSaidiL.AliJ. B.FnaiechF.Bi-spectrum based-EMD applied to the non-stationary vibration signals for bearing faults diagnosis20145351650166010.1016/j.isatra.2014.06.002ChengJ.YuD.YangY.A fault diagnosis approach for gears based on IMF AR model and SVM2008200864713510.1155/2008/6471352-s2.0-45749124824LiZ.ChenJ.ZiY.PanJ.Independence-oriented VMD to identify fault feature for wheel set bearing fault diagnosis of high speed locomotive2017855125292-s2.0-8499542706810.1016/j.ymssp.2016.08.042HouL.ChenY.CaoQ.LuZ.Nonlinear vibration analysis of a cracked rotor-ball bearing system during flight maneuvers20161055155282-s2.0-8498035053110.1016/j.mechmachtheory.2016.07.024MandicD. P.GohS. L.AiharaK.Sequential data fusion via vector spaces: fusion of heterogeneous data in the complex domain2006481-29910810.1007/s11265-006-0025-62-s2.0-34547327934GoldmanP.MuszynskaA.Application of full spectrum to rotating machinery diagnostics, Orbit19991721RillingG.FlandrinP.GoncalvesP.LillyJ. M.Bivariate empirical mode decomposition200714129369392-s2.0-3674905225010.1109/LSP.2007.904710LooneyD.MandicD. P.Multiscale image fusion using complex extensions of EMD20095741626163010.1109/TSP.2008.2011836MR2590469TsoumasI.GeorgoulasG.SafacasA.VachtsevanosG.Empirical mode decomposition of the stator start-up current for rotor fault diagnosis in asynchronous machinesProceedings of the 2008 International Conference on Electrical Machines, ICEM'08September 2008Vilamoura, Portugal1610.1109/ICELMACH.2008.47999872-s2.0-65249124691MollaM. K. I.GhoshP. R.HiroseK.Bivariate EMD-based data adaptive approach to the analysis of climate variability201120112193503410.1155/2011/9350342-s2.0-80052654825YangW. X.CourtR.TavnerP. J.CrabtreeC. J.Bivariate empirical mode decomposition and its contribution to wind turbine condition monitoring2011330153766378210.1016/j.jsv.2011.02.0272-s2.0-79955522082PatelT. H.DarpeA. K.Use of full spectrum cascade for rotor rub identification2009821391512-s2.0-67849107260OkatanA.HajiyevC.HajiyevaU.Fault detection in sensor information fusion Kalman filter20096397627682-s2.0-6774908467710.1016/j.aeue.2008.06.003LiuX.MaL.MathewJ.Machinery fault diagnosis based on fuzzy measure and fuzzy integral data fusion techniques20092336907002-s2.0-5804913640610.1016/j.ymssp.2008.07.012ElbhbahK.SinhaJ. K.Vibration-based condition monitoring of rotating machines using a machine composite spectrum201333211283128452-s2.0-8487521010110.1016/j.jsv.2012.12.024WangX. J.CaiY. P.LinX. Z.ICE fault diagnosis method based on mutual information and WVD time-frequency analysis20145257417452-s2.0-8489685970610.4028/www.scientific.net/AMM.525.741LiF.WangJ.ChyuM. K.TangB.Weak fault diagnosis of rotating machinery based on feature reduction with Supervised Orthogonal Local Fisher Discriminant Analysis20151685055192-s2.0-8493781444610.1016/j.neucom.2015.05.076ShaoH.JiangH.WangF.ZhaoH.An enhancement deep feature fusion method for rotating machinery fault diagnosis20171192002202-s2.0-8500821965010.1016/j.knosys.2016.12.012LuN.XiaoZ.MalikO. P.Feature extraction using adaptive multiwavelets and synthetic detection index for rotor fault diagnosis of rotating machinery201552-5313934152-s2.0-8492289704110.1016/j.ymssp.2014.07.024ZhaoX.PatelT. H.ZuoM. J.Multivariate EMD and full spectrum based condition monitoring for rotating machinery201227171272810.1016/j.ymssp.2011.08.0012-s2.0-82255179549GoldmanP.A Muszynska . Application of full spectrum to rotating machinery diagnostics19991721WangH.-Q.HouW.TangG.YuanH.-F.ZhaoQ.-L.CaoX.Fault detection enhancement in rolling element bearings via peak-based multiscale decomposition and envelope demodulation201420141132945810.1155/2014/329458ChuF.LuW.Experimental observation of nonlinear vibrations in a rub-impact rotor system20052833–562164310.1016/j.jsv.2004.05.0122-s2.0-15944386375DangJ.JiaR.LuoX.WuH.ChenD.Partly duffing oscillator stochastic resonance method and its application on mechanical fault diagnosis2016201614310938510.1155/2016/31093852-s2.0-85006049665