An improved local mean decomposition method and its application for fault diagnosis of reciprocating compressor

Gui-juan Chen1 , Long-qing Zou2 , Hai-yang Zhao3 , Yu-qian Li4

1, 2, 3, 4Mechanical Science and Engineering Institute, Northeast Petroleum University, Daqing, China

2Corresponding author

Journal of Vibroengineering, Vol. 18, Issue 3, 2016, p. 1474-1485. https://doi.org/10.21595/jve.2016.16594
Received 28 October 2015; received in revised form 12 February 2016; accepted 3 March 2016; published 15 May 2016

Copyright © 2016 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

Local mean decomposition (LMD) is a new time–frequency analysis method which can decompose a signal adaptively into a set of product function (PF) components, and the construction of local mean function and envelope function plays an important role in the accuracy of its PF components. Aiming at the strong nonstationarity, nonlinearity and multi-component coupling characteristics of reciprocating compressor vibration signals, an improved LMD was proposed by a novel construction method of local mean function and envelope function. By introducing an extreme symmetrical point between two extreme points and using the Monotone Piecewise Cubic Hermite Interpolation (MPCHI) instead of Cubic Spline Interpolation (CSI) to construct the envelopes, a novel construction method of local mean function and envelope function was proposed, and then the improved LMD algorithm was given based on this novel construction method. The improved LMD was applied to decompose the vibration signals of reciprocating compressor fault states, and the comparison of details between different LMD decomposition results verified the superiority of this improved method. The envelope frequency spectrum of PF component gives a more significant peak of fault frequency than that of original signal, which further indicates that this proposed method is competent for the diagnosis of reciprocating compressor oversized bearing clearance fault.

Keywords: LMD, MPCHI, reciprocating compressor, fault diagnosis, bearing clearance.

1. Introduction

Reciprocating compressors are one of the widely used machinery in petroleum and chemical industry, and their operational state has become a key subject of the intensive investigation for the high demand for effective performance, safety and convenient maintenance [1, 2]. Due to the abundant operational state information embodied in vibration signal, vibration analysis technique is highly effective in fault diagnosis of machinery and equipment [3, 4], and some of those techniques have been successfully applied in fault diagnosis of rotating machinery such as gearbox and bearing [5, 6]. However, due to the factors such as nonlinear stiffness of bearings and the unbalanced and time-varying forces of components, the vibration signals of reciprocating compressor present nonlinearity, nonstationarity and multi-component coupling characteristics. The traditional vibration analysis techniques may ineffective to detect fault features from the vibration signals of reciprocating compressor.

In 2005, a new adaptive analysis method, which is the local mean decomposition (LMD), was initially developed by Smith [7]. LMD decomposes signals into a set of product functions (PFs), each of which is the product of an envelope signal and a purely frequency-modulated signal. Moreover, the complete time frequency distribution of the original signal could be acquired by the instantaneous amplitude and instantaneous frequency of PF components. Therefore, LMD is an appropriate method to analyze nonstationary, nonlinear and multi-component coupled vibration signal of reciprocating compressor. Recently, LMD has become a significant subject of the intensive investigation for the vibration signal feature extraction methods. Wang et al. [8] had given a novel LMD demodulation method to extract features of a gas turbine system through vibration signals, and they detected the rub-impact fault accurately. Cheng [9] employed LMD to diagnose bearing and gear faults effectively and accurately. Subsequently, they proposed an ensemble LMD method to solve the mode mixing issue of the LMD and applied it to rotor fault diagnosis.

The construction of local mean function and envelope function is the key step in LMD algorithm, and this plays an important role in the accuracy of the PF components it extracts. In the traditional LMD algorithm, the local mean function and envelope function are calculated by using the moving averaging method according to the local extreme points. While this method may cause phase error of the functions, the error will be significant after several iterations and has an adverse effect on the accuracy of PF components [10]. Inspired by the Empirical mode decomposition (EMD) algorithm, Hu and Ren [11] has presented a novel approach to calculate the local mean function and envelope function using the envelopes of local maximum points and local minimum points which are obtained by CSI. This approach can not only improve the computational efficiency, but also avoid the selection of the step size for moving averaging method. While the envelopes obtained by CSI may occur overshoot and undershoot problems due to its continuous second derivative for vibration signals with strong nonstationary characteristic. Piecewise Cubic Hermite Interpolation (PCHI) is also a widely used interpolation method, and its derivative is only first order continuous and can be set monotonically [12, 13]. Hence, to calculate the local mean function and envelope function by using MPCHI is a suitable approach for strong nonstationary signals for it can avoid the overshoot and undershoot problems.

In addition, the local mean function and envelope function are exclusively determined by extreme, and ignore all other critical signal information, this fact contributes to mislocalization of temporal information and produces spurious phase shifts and distortion in the PF components it extracts. Frei et al. [14] has presented a new method for time-frequency-energy analysis, called the Intrinsic Time-scale Decomposition (ITD). Different from LMD method, ITD extracts a baseline signal from original signal by using extreme and its symmetrical points, and this manner causes the residual to be a proper rotation, for which instantaneous frequency and amplitude are well defined. The symmetrical points of extreme enrich the temporal information of original signal, therefore, improve the accuracy of components it extracts. By introducing the symmetrical point to the construction of envelopes between extremes in LMD algorithm may also have a positive effect on the accuracy of local mean function and envelope function.

This paper presents an improved LMD algorithm by using MPCHI and symmetrical points to construct the local mean function and envelope function, and is organized as follows. And an improved LMD algorithm using MPCHI and symmetrical points is presented in Section 2. Section 3 compares the performance of this method by using the simulation signal. In Section 4, the novel method is applied to the diagnosis of reciprocating compressor bearing clearance fault. Finally, we offer the conclusions in Section 5.

2. An improved LMD method

2.1. Monotone piecewise cubic Hermite interpolated (MPCHI) envelope

PCHI is an effective interpolation method and its performance depends on the first derivatives of the interpolation points. By determining suitable derivative-value at the interpolation points, PCHI can preserve monotonicity of envelopes. Therefore, for signals with strong nonstationary characteristic, to construct envelopes by MPCHI may overcome the overshoot problem of envelopes obtained by CSI.

The MPCHI can be defined as follow [12, 13]. Let (xi, Yi, di), i= 0,..., n be given real data, where a=x0<x1<<xn=b is a partition of the interval [a,b], while yi and di,i= 0,..., n are assigned function values and derivative-values respectively. Further let hi, i, yi, i= 0,..., n, be defined by:

h i = x i + 1 - x i ,         y i = y i + 1 - y i ,         i = y i / h i .

We suppose that the data are monotonic, i.e. yi0i or yi0i and:

d i = d i + 1 = 0 , i = 0 , s g n d i = s g n d i + 1 = s g n i , i 0 ,

where we assume the convention sgn0=sgn(i).

The monotone piecewise cubic function S(x)C1a,b interpolating the given values, i.e. Sxi=yi and S'xi=di, i= 0,…, n, is defined as:

(1)
S ( x ) S i ( x ) = ( d i + d i + 1 - 2 i ) h i 2 ( x - x i ) 3
            + ( - 2 d i - d i + 1 + 3 i ) h i ( x - x i ) 2 + d i x + x i + y i ,

for xxi,xi+1.

A segment of vibration signal with strong nonstationary characteristic is shown in Fig. 1, and the maxima and minimum points were enveloped by CSI and MPCHI respectively. We can observe that the upper and lower envelopes constructed by CSI occur overshoot and undershoot problems, while MPCHI overcome this problem due to its monotone characteristic.

Fig. 1. The comparison of envelopes obtained by CSI and MPCHI

 The comparison of envelopes obtained by CSI and MPCHI

2.2. The symmetrical point of extreme

ITD is an adaptive time-frequency analysis method which decomposes a complex signal into several proper rotation components. For a given signal Xt, ITD extracts a baseline signal from Xt in a manner that causes the residual to be a proper rotation. More specifically, Xt can be decomposed as:

(2)
X t = L t + H t ,

where Lt is the baseline signal and Ht is a proper rotation.

Suppose that Lt and Ht have been defined on [0, τk] and that Xt is available for t[0, τk+2]. We can then define a (piece-wise linear) baseline Lt on the interval (τk, τk+1] between successive extremes as follows:

(3)
L t = L k + L k + 1 - L k X k + 1 - X k ( X t - X k ) ,         t ( τ k , τ k + 1 ] ,

where:

(4)
L k + 1 = α X k + τ k + 1 - τ k τ k + 2 - τ k ( X k + 2 - X k ) + ( 1 - α ) X k + 1 ,

and 0 <a< 1 is typically fixed with a= 1/2.

We defined Ak+1 as:

(5)
A k + 1 = X k + τ k + 1 - τ k τ k + 2 - τ k ( X k + 2 - X k ) ,

and then fixed a with 1/2, the Eq. (4) can be expressed as:

(6)
L k + 1 = 1 2 A k + 1 + X k + 1 .

ITD constructs the baseline signal, Lt, in this manner in order to maintain the monotonicity of Xt between extrema. Fig. 2 shows the construction of baseline signal. Firstly, point Ak+1 is calculated by using two adjacent maxima Xk and Xk+2 and a minimum Xk+1 between them based on Eq. (5). Then, the baseline signal point Lk+1 is obtained by using point Ak+1 and Xk+1 according to Eq. (6). Owing to point Ak+1 and the minimum Xk+1 are symmetrical with respect to the baseline signal point Lk+1, Ak+1 is named as the symmetrical point of extreme. Finally, the baseline signal is created by the piecewise linear interpolation according to Eq. (3) after all baseline signal points are calculated. While for the same signal in Fig. 2, the local mean point Mk+1 in LMD is determined by the minima Xk+1 and the envelope point Ek+1 which derived by using MPCHI and the two adjacent maxima Xk and Xk+2. Compared to the envelope point Ek+1 in LMD which is determined by the interpolation method, the symmetrical point of extreme Ak+1 in ITD enable more information ‘intrinsic’ to the original signal to be passed down to the baseline. Therefore, inserting a symmetrical point between two adjacent extreme points and then constructing the envelopes by using MPCHI may be an effective approach to improve the accuracy of local mean function and envelope function, and hence improve the accuracy of PF components it extracts.

Fig. 2. The comparison of local mean function in LMD and baseline signal in ITD

 The comparison of local mean function in LMD and baseline signal in ITD

2.3. The algorithm and flowchart of the improved LMD

To improve the accuracy of PF components, the novel LMD method constructs the local mean function and envelope function in a new manner compared to the traditional algorithm, in which the upper and lower envelopes are calculated by the extreme points and the symmetrical points using MPCHI. For any signal xt, the local mean function and envelope function can be constructed as following steps:

(1) Set the initial parameters: uit=xt, sijt=uit, i= 0, j= 0;

(2) Determine all the extrema series ni,j,k of the signal sijt;

(3) Calculate all the symmetrical points series Ai,j,k of extrema series ni,j,k based on Eq. (5);

(4) Separate the extrema series ni,j,k into the maximum series nmaxi,j,k and the minimum series nmini,j,k, and then insert the symmetrical points Ai,j,k between two adjacent extrema to form the new maximum series nmaxi,j,k and minimum series nmini,j,k;

(5) Calculate the upper envelope Euijt and lower envelope Elijt by the maximum series nmaxi,j,k and minimum series nmini,j,k using MPCHI based on Eq. (1);

(6) The continuous local mean function mijt and the envelope function aijt are given by the upper envelope Euijt andthe lower envelope Elijt as:

(7)
m i j ( t ) = E u i j t + E l i j t 2 ,
(8)
a i j t = E u i j t - E l i j t 2 .

Fig. 3. Flowchart of improved LMD method

 Flowchart of improved LMD method

Following the rest steps of traditional LMD algorithm [15], the signal xt can be decomposed into a set of PF components. The flowchart of improved LMD algorithm based on above steps is shown in Fig. 3.

3. Simulated signal analysis

To validate the effectiveness of the proposed LMD method, a numerical example is provided here. The simulated signal is displayed as:

(9)
x t = x 1 t + x 2 t , x 1 ( t ) = 0.8 s i n ( 6 π t ) e t / 2 x 2 t = c h i r p t , 5,2 , 60 , ,

where t [0, 2], x1t is an AM sinusoidal signal, and x2t is a chirp signal. Fig. 4 shows the time domain waveform of xt.

To compare the performance of the improved LMD, the simulated signal xt was decomposed by the traditional LMD, the cubic spline LMD and the improved LMD, and the results are shown in Fig. 5, Fig. 6 and Fig. 7, respectively. Obviously, three methods all succeed in extracting two PF components from signal xt. The first component PF1 corresponds to the chirp signal x2t, and the second component PF2 corresponds to the AM sinusoidal signal x1t in results of all three methods.

Fig. 4. Time-domain waveform of simulated signal

 Time-domain waveform of simulated signal

Fig. 5. Decomposition results of simulated signal with the traditional LMD

 Decomposition results of simulated  signal with the traditional LMD

Fig. 6. Decomposition results of simulated signal with the cubic spline LMD

 Decomposition results of simulated  signal with the cubic spline LMD

Fig. 7. Decomposition results of simulated signal with the improved LMD

 Decomposition results of simulated  signal with the improved LMD

The Mean Squared Error (MSE) between the obtained PF components and its corresponding simulated signals were calculated, which were listed in Table 1 together with the iterations of each PF component and the total calculation time. The MSEs of PF components in the improved LMD are lower than those of the traditional LMD and the cubic spline LMD. This implies the obtained PF components of the improved LMD have the higher accuracy than those of other two methods.

The iterations of PF component reflect the accuracy of local mean function and envelope function directly, and the less iterations one PF component needs, the higher the accuracy of local mean function and envelope function is. So the improved LMD has the least iterations of PF components than other two methods due to the introduction of extreme symmetrical points. Furthermore, the whole decomposition of the improved LMD takes 0.478 s, while the traditional LMD and the cubic spline LMD need 0.832 s and 0.678 s for the same procedure respectively. This indicates the MPCHI used in the improved LMD is a more efficient method.

Table 1. Performance of decomposition results for simulated signal with different LMD methods

Method
PF1
PF2
Time / s
MSE
Iterations
MSE
Iterations
Traditional LMD
0.0061
5
0.0224
4
0.843
CSI LMD LMD
0.0053
4
0.0198
3
0.665
Improved LMD
0.0047
2
0.0191
2
0.478

4. Fault diagnosis of reciprocating compressor

A 2D12 type reciprocating compressor, which is widely used to compress natural gas in chemical and petroleum industry, was selected as research object. The shaft power of reciprocating compressor is 500 kW, and the piston stroke and motor rotation speed are 240 mm and 496 rpm respectively. The pistons were derived by a two-stage crank-connecting rod transmission mechanism, and the connecting rod was assembled to crankshaft and crosshead pin with sliding bearings respectively [16-18]. Due to tolerances and defects in manufacturing process or wearing, bearing clearance will increase after a certain working period [19, 20]. When the bearing clearances are oversized, the contact forces between bearings generate impulsive effect, and this causes vibratory running condition. Thus, the effective and accurate diagnosis of oversized bearing clearance states plays an important role in reducing maintenance costs.

In order to diagnose the bearing clearance states by vibration signal, the oversized clearance fault state of bearing between the crankshaft pin and first stage connecting rod was tested by using the worn bushing, and the normal state of this bearing was also tested for comparison. An ICP acceleration sensor, as the triangle in Fig. 8, was placed on the top of crosshead guide surface to collect vibration signal.

Fig. 8. Distribution of reciprocating compressor measurement points

 Distribution of reciprocating compressor measurement points

The typical vibration acceleration signals of the acceleration sensor in normal bearing clearance state and oversized clearance fault state were shown in Fig. 9 and Fig. 10 for a period of two crank revolutions respectively. Compared to the chaotic vibration acceleration of normal state in Fig. 9, the vibration acceleration of oversized clearance fault in Fig. 10 shows periodical shocks, and the amplitude of shocks increases significantly. In normal bearing clearance state, the impact excitation force is insignificant compared to other excitation forces because the journal and bearing contacts continuously with oil film. While in oversized clearance fault state, the excitation force is enhanced significantly for the strong impacts between journal and bearing. So the vibration acceleration in Fig. 10 shows periodical and strong shocks.

Fig. 9. The vibration acceleration in normal bearing clearance state

 The vibration acceleration  in normal bearing clearance state

Fig. 10. The vibration acceleration in oversized bearing clearance state

 The vibration acceleration  in oversized bearing clearance state

According to the theory of structural dynamics, the vibration acceleration of a structure can be considered as the superposition of different natural frequencies of this structure modulated by the excitation forces. The envelope of vibration acceleration can reflect the feature of the excitation forces, so it is an effective approach to study the envelope of vibration acceleration for the diagnosis of bearing clearance state. The journal and bearing usually occurs two impacts for a period of one crank revolution. The envelope frequency spectrum of vibration acceleration in oversized clearance fault state was shown in Fig. 11, and we can see a peak value in double rotation frequency. While many excitation sources exist in the reciprocating compressor, some strong peaks values also appear in other frequencies, which may disturb the accuracy of fault diagnosis. So the fault information should be further extracted from the collected vibration acceleration signal for an accurate diagnosis of bearing clearance state.

As mentioned above, the collected vibration signal can be considered as a multi-component signal which is superimposed by different modulated natural frequencies signals, and each modulated natural frequencies signal is a mono-component signal. While LMD method can decompose a multi-component signal into a series of mono-component signals, that is, PF components, each of which is the product of an envelope signal and a purely frequency-modulated signal. Therefore, LMD is an appropriate method to extract fault information from the vibration signal of reciprocating compressor.

Fig. 11. The envelope frequency spectrum of vibration acceleration in oversized bearing clearance fault state

 The envelope frequency spectrum of vibration acceleration  in oversized bearing clearance fault state

However, the vibration signal of reciprocating compressor has characteristics of strong nonstationarity, and using the traditional LMD method to decompose the vibration signal may have an adverse effect on the accuracy of PF components. This new LMD method is proposed to improve the accuracy of local mean function and envelope function for strong nonstationary signals by introducing the symmetrical points and using the MPCHI. The improved LMD was applied to decompose the typical vibration acceleration of oversized clearance fault state in Fig. 10, and 7 PF components were derived. Owing to the first several PF components contains the main fault information, the first three PF components of this decomposition results were shown in Fig. 12. To compare the performance of the improved LMD, the vibration signal was also decomposed by the traditional LMD and the CSI LMD, and 9 and 11 PF components were derived respectively. The first three PF components of two methods were illustrated in Fig. 13 and Fig. 14 respectively.

Fig. 12. Decomposition results of vibration signal with the improved LMD

Decomposition results of vibration  signal with the improved LMD

Fig. 13. Decomposition results of vibration signal with the traditional LMD

 Decomposition results of vibration  signal with the traditional LMD

Fig. 14. Decomposition results of vibration signal with the CSI LMD

 Decomposition results of vibration signal with the CSI LMD

To clearly show the details of decomposition results for different LMD methods, 100 vibration signal samples are selected from oversized clearance fault state, and they were decomposed by different LMD methods respectively. The average value and 3δ (three times standard deviation) of iterations of each PF component, the Average Index of Orthogonality (IOav) and the Index of Energy Conservation (IEC) were listed in Table 2.

As mentioned above, the iterations of PF component reflect the accuracy of local mean function and envelope function directly. By introducing the symmetrical points and using the MPCHI, all the three PF components of the improved LMD have the least iterations. While the CSI LMD needs the most iteration to finish the decomposition, and we can see that the third PF component has a mutation. This is because the envelope of extreme interpolated by cubic spline curve occur overshoot and undershoot problems due to its continuous second derivative for this strong nonstationary vibration signal, and therefore decrease the accuracy of local mean function and envelope function.

Table 2. Performance of decomposition results for vibration acceleration signal with different LMD methods

Method
Iterations
I O a v e
I E C
PF1
PF2
PF3
CSI LMD
14.15±0.52
10.15±0.23
18.28±0.48
0.344±0.053
0.924±0.063
Traditional LMD
7.32±0.43
7.08±0.54
8.51±0.71
0.213±0.032
0.977±0.045
Improved LMD
5.54±0.23
3.47±0.31
5.27±0.027
0.189±0.021
0.984±0.023

Theoretically, each two PF components of LMD decomposition results should be orthogonal, and the Index of Orthogonality (IO) between two PF components is equal to zero, but in fact there are errors of the local mean function and envelope function, the orthogonality between two PF components is relative. So the more accurate the PF components are, the closer to zero the IO is. To compare all PF components of a LMD method, the IOav is used. From another point of view, if the PF components of LMD decomposition results are orthogonal, the energy of a signal before and after decomposition should be conserved, and the IEC is equal to one. Usually the more accurate the PF components are, the closer to one the IEC is [21, 22]. The IOav and IEC are defined as:

I O a v e = A v e r a g e p f j , p f k p f j 2 p f k 2 ,       j k ,
I E C = t = 0 T j = 1 n p f j t 2 t = 0 T x t - u n t 2 ,

where xt is the original signal, pfjt is PF component and unt is the residual component. Based on the results as summarized in Table 2, we can see that the improved LMD indeed give very better results compare to other two methods. This further validates the superiority of improved LMD for strong nonstationary vibration signal.

In addition, the 3δ of PF component iterations, IOav and IEC in the improved LMD are all smaller than those for other two methods. The results indicate that the proposed method is stable and reliable to deal with signals contaminated with noise.

As discussed above, LMD method can decompose a multi-component signal into a series of PF component, each of which is the product of an envelope signal and a purely frequency-modulated signal, and the envelope signal can reflect the feature of the excitation forces. So it is an effective approach to study the envelope of PF component for a more accurate diagnosis of bearing clearance state. The first PF component of the improved LMD decomposition results in Fig. 12 contains more bearing clearance fault information for it was dominated with the periodical and strong shocks compare to other PF components. The envelope frequency spectrum of the first PF component was shown in Fig. 15, and we can also see a peak value in double rotation frequency. Though the envelope frequency spectrum of original vibration signal appears the peak value of double rotation frequency in Fig. 11, this peak value in Fig. 15 is prominent compared to other peaks, in which is easy to find the fault frequency of bearing clearance state.

Fig. 15. The envelope frequency spectrum of the first PF component in oversized bearing clearance fault state

 The envelope frequency spectrum of the first PF component  in oversized bearing clearance fault state

5. Conclusions

According to the strong nonstationarity, nonlinearity and multi-component coupling characteristics of vibration signals, an improved LMD was proposed through the novel construction method of local mean function and envelope function, and this method diagnosed the bearing clearance fault of reciprocating compressor effectively.

By introducing the extreme symmetrical point and using the MPCHI instead of cubic spline to construct the envelopes, a novel construction method of local mean function and envelope function was proposed, and an improved LMD algorithm was give based on the novel construction method.

The improved LMD was applied for the diagnosis of reciprocating compressor bearing clearance fault. Compared to the traditional LMD and the cubic spline LMD, the improved LMD has advantages on the iterations of each PF component, IOav and IEC of decomposition results, and the envelope frequency spectrum of PF component given a more significant peak value of fault frequency.

Acknowledgements

This work was partly supported by General Financial Grant from The China Postdoctoral Science Foundation (2015M581423), School Cultivate Fund of Northeast Petroleum University in China (XN2014105) and Natural Science Foundation of Heilongjiang Province in China (E2015037).

References

  1. Elhaj M., Gub F., Ballb A. D. Numerical simulation and experimental study of a two-stage reciprocating compressor for condition monitoring. Mechanical Systems and Signal Processing, Vol. 22, 2008, p. 374-389. [Publisher]
  2. Almasi A. A new study and model for the mechanism of process reciprocating compressors and pumps. Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering, Vol. 224, 2010, p. 143-148. [Publisher]
  3. Vakharia V., Gupta V. K., Kankar P. K. A multiscale permutation entropy based approach to select wavelet for fault diagnosis of ball bearings. Journal of Vibration and Control, 2014, p. 1-9. [Search CrossRef]
  4. An X. L., Jiang D. X., Chen J. Application of the intrinsic time-scale decomposition method to fault diagnosis. Journal of Vibration and Control, Vol. 18, 2011, p. 240-245. [Publisher]
  5. Yang Y., Pan H. Y., Ma L., Cheng J. S. A fault diagnosis approach for roller bearing based on improved intrinsic timescale decomposition de-noising and kriging-variable predictive model-based class discriminate. Journal of Vibration and Control, 2014, p. 1-16. [Search CrossRef]
  6. Feng Z. P., Liang M. Complex signal analysis for wind turbine planetary gearbox fault diagnosis via iterative atomic decomposition thresholding. Journal of Sound and Vibration, Vol. 333, 2014, p. 5196-5211. [Publisher]
  7. Smith J. S. The local mean decomposition and its application to EEG perception data. Journal of the Royal Society Interface, Vol. 2, 2005, p. 443-454. [Publisher]
  8. Wang Y. X., He Z. J., Zi Y. Y. A demodulation method based on improved local mean decomposition and its application in rub-impact fault diagnosis. Measurement Science and Technology, Vol. 20, 2009, p. 1-10. [Publisher]
  9. Cheng J. S., Yang Y. A rotating machinery fault diagnosis method based on local mean decomposition. Digital Signal Processing, Vol. 22, 2012, p. 356-366. [Publisher]
  10. Zhang K., Cheng J. S., Yang Y. The local mean decomposition method based on rational spline and its application. Journal of Vibration Engineering, Vol. 24, 2009, p. 97-103. [Search CrossRef]
  11. Hu J. S., Yang S. X., Ren D. Q. Spline-based local mean decomposition method for vibration signal. Journal of Data Acquisition and Processing, Vol. 24, 2009, p. 78-83. [Search CrossRef]
  12. Merrien J. L., Sablonnière P. Rational splines for Hermite interpolation with shape constraints. Computer Aided Geometric Design, Vol. 30, 2013, p. 296-309. [Publisher]
  13. Cripps R. J., Hussain M. Z. C1 monotone cubic Hermite interpolant. Applied Mathematics Letters, Vol. 25, 2012, p. 1161-1165. [Publisher]
  14. Frei M., Osorio I. Intrinsic time-scale decomposition: time-frequency-energy analysis and real-time filtering of non-stationary signals. Proceedings of the Royal Society, Vol. 463, 2007, p. 321-342. [Publisher]
  15. Zhao Hai-yang, Wang Jin-dong, Xing Jun-jie, Gao Yi-qi A feature extraction method based on LMD and MSE and its application for fault diagnosis of reciprocating compressor. Journal of Vibroengineering, Vol. 17, Issue 7, 2015, p. 3515-3526. [Search CrossRef]
  16. Gu F., Ball A. D. Automating the diagnosis of valve faults in reciprocating compressors. Maintenance, Vol. 10, 1996, p. 3-8. [Search CrossRef]
  17. Grando F. P., Priest M., Parata A. T. Lubrication in refrigeration system: numerical model for piston dynamics considering oil-refrigerant interaction. Proceedings of the Institution of Mechanical Engineers, Journal of Engineering Tribology, Vol. 220, 2006, p. 245-258. [Search CrossRef]
  18. Prata A. T., Fernandes J. R. S., Fagotti F. Dynamic analysis of piston secondary motion for small reciprocating compressors. Journal of Tribology Transactions of the ASME, Vol. 122, 1995, p. 752-760. [Publisher]
  19. Flores P., Ambrosio J. Revolute joints with clearance in multibody systems. Computers and Structures, Vol. 82, 2006, p. 1359-1369. [Publisher]
  20. Parenti C. V., Venanzi S. Clearance influence analysis on mechanisms. Mechanism and Machine Theory, Vol. 40, 2005, p. 1316-1329. [Publisher]
  21. Huang N. E., Shen Z., Long Wu S. R. M. C., Shih H., Zheng Q., Yen N., Tung C., Liu H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis. Proceedings of Royal Society London A, Vol. 454, 1998, p. 903-995. [Publisher]
  22. Chen Q. H., Huang N. S., Riemenschneider S. A B-spline approach for empirical mode decompositions. Advances in Computational Mathematics, Vol. 24, 2006, p. 171-195. [Publisher]

Cited By

Sensors
Yang Liu, Lixiang Duan, Zhuang Yuan, Ning Wang, Jianping Zhao
2019
Journal of Vibroengineering
Shungen Xiao, Shulin Liu, Shouguo Cheng, Xiaogang Xue, Mengmeng Song, Xin Sun
2017