Vibrationbased classification of centrifugal pumps using support vector machine and discrete wavelet transform
Ebrahim Ebrahimi^{1} , Mohammad Javidan^{2}
^{1, 2}Department of Mechanical Engineering, College of Engineering, Islamic Azad University, Kermanshah Branch, Kermanshah, Iran
^{1}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 4, 2017, p. 25862597.
https://doi.org/10.21595/jve.2017.18120
Received 19 December 2016; received in revised form 25 May 2017; accepted 27 May 2017; published 30 June 2017
JVE Conferences
Due to the quick advancement of technology, application of different methods is highly required to maintain the high quality of production and health assessment of production lines. Hence, condition monitoring is widely used in the industry as an efficient approach. The purpose of the present study was to classify faults in centrifugal pumps using the vibration signal analysis and support vector machine (SVM) method. Vibration signals were decomposed in three levels by Daubechies wavelets, and a total of 44 descriptive statistical features were extracted from detail coefficients and approximation coefficients of the wavelets. In order to find the best model for fault classification of centrifugal pumps, parameters such as penalty, degree of polynomial, and width of the Gaussian radial basis function kernel (RBF kernel) were investigated. The classification results using the SVM method indicated that the maximum classification accuracy was 96.67 percent, which was obtained at an RBF kernel width of 0.1 and a penalty parameter value of 1.
Keywords: centrifugal pumps, fault detection, support vector machine (SVM).
1. Introduction
Due to the advancement of technology, industrial equipment is increasingly made more complex, requiring more careful attention to be paid to them as their failure and breakdown can cause significant costs. Hence, factors such as reliability, accessibility, reduction of failure duration, and reparability of equipment are of great importance. To this goal, condition monitoring has been of high interest as an efficient method to enhance the factors of safety, health and optimal performance of machineries.
Condition monitoring has been defined as fault detection and maintenance of equipment without interrupting their activities [1]. Generally, this method is based on regular data acquisition from the dynamic characteristics of the equipment, as well as comparison of the data with those of healthy conditions. In traditional condition monitoring, fault detection was usually based on analysis of either vibration or acoustic data [2]. Various methods have been introduced and implemented for singlesensor condition monitoring based on a single characteristic such as vibration [35] and acoustic [68] using classifiers such as support vector machines (SVMs) [5, 9, 10] and artificial neural networks (ANNs) [11].
Vibration behavior analysis techniques have been widely used in researches on different subjects such as analysis of natural frequencies [12], acquisition of real mechanical behaviors of micro/nano structures [1315] and fault diagnosis in rotational machinery [1618].
Centrifugal pumps play vital role in many critical applications and therefore continuous availability of such mechanical components becomes an absolute essential. Pumps are the key elements in waste water treatment plants, food industries, agriculture, oil and gas industries,
paper and pulp industries, etc. [19]. Hence, fault detection in this type of pump is considerably important in order to prevent further failures and breakdowns.
In recent years numerous researches have been focused on the vibration analysis and intelligent method for fault diagnosis. Kiran et al. [20] have presented a new method for signal processing by using a fast Fourier transform and wavelet transform for fault diagnosis of the gear used in an internal combustion engine. Li et al. [21] have developed a new planetary gearboxes fault diagnosis method based on vibration signals and the support vector machine. Yunlong and Peng [22] have developed a method based on the empirical mode decomposition (EMD) and SVM for fault classification of centrifugal pump. Li et al. [18] presented a novel feature extraction algorithm based on the multiscale permutation entropy (MPE) and SVM for rolling bearing fault diagnosis. Sakthivel et al. [23] have presented an approach based on C4.5 decision tree algorithm for fault detection of monoblock centrifugal pumps. Cui et al. [24] have developed a method based on the information entropy and SVM for compressor valve fault diagnosis. Farokhzad et al. [11] have developed a fault diagnosis technique for fault detection of centrifugal pumps based on their vibration behavior in a centrifugal pump in frequency domain and artifical neural network. Therefore, in this study we presented a method based on the wavelet transform and SVM for detection and classification of faults in centrifugal pumps.
2. Background materials
2.1. Signal processing
The vibration diagnosis is normally carried out in the following main steps: signal measurement, signal analysis, diagnosis and strategic decision, where the signal analysis plays a key role and has the task of extracting useful information, filtering noise from a measured vibration signal and finding the fault feature and its developing trend. Traditional spectral analysis techniques based on the Fourier transform provide a good description of stationary and pseudo stationary signals. Unfortunately, these techniques have several shortcomings. Therefore, in recent years, there has been an increasing interest in the research of signal analysis concerning the timefrequency domain. The timefrequency analyses are able to analyze a non stationary signal and indicate not only which frequencies the signal contains, but also when these frequencies occur [25]. The best approach for analysis of a non stationary vibration signal in a timefrequency domain is a wavelet transform.
2.1.1. Wavelet transform (WT)
The use of a wavelet transform (WT) was very popular for two decades, due to advantageous properties of this transformation and availability of computer software. The WT is used in different fields of science, like medicine, biology and engineering; it is also employed to process signals and images [26]. Generally, conventional data processing is computed in a time or frequency domain. Wavelet processing combines both time and frequency. In simple language, we use the term ‘timefrequency’ analysis. A wavelet is a basis function characterized by two aspects; one is its shape and amplitude that is chosen by the user and the other is its scale (frequency) and time (location) relative to the signal [27]. The wavelet transform of signal $f\left(t\right)\in {L}^{2}\left(R\right)$ is defined by the innerproduct between ${\psi}_{ab}\left(t\right)$ and $f\left(t\right)$ as [28]:
where $a$ and $b$ are the scaling (dilation) and translation (time shift) constants respectively, and $\psi $ is the wavelet function that may not be real as assumed in the above equation for simplicity. The choice of the wavelet function (mother wavelet) is flexible provided that it satisfies the so called admissibility conditions [28]. The discrete wavelet transform (DWT) is derived from the discretization of continuous wavelet transform (CWT). The most common discretization is dyadic. The DWT is given by:
where, the parameters $a$ and $b$ in Eq. (1) are replaced by ${a}_{0}^{m}$ and $k{a}_{0}^{m}$, respectively, with $k$ and $m$ being integer variables [28].
The discrete signal is passed through a high pass filter (H) and a low pass filter (L), resulting in two vectors at the first level; approximation coefficient (A1) and detail coefficient (D1). Application of the same transform on the approximation (A1) causes it to be decomposed further into approximation (A2) and detail (D2) coefficients at the second level. Finally, the signal is decomposed at the expected level. The approximations are the highscale lowfrequency components of the signal, whereas details are its lowscale, highfrequency components. The wavelet decomposition for level 3 is illustrated in Fig. 1 [29].
Each vector ${A}_{j}$ includes approximately $N/{2}^{j}$^{}coefficients, where $N$ is the number of data points in the input signal $s$, and provides information about a frequency band [0, $Fs/{2}^{j+1}$], where $Fs$ is the sampling frequency. In Fig. 1, H and L represent the decomposition filters, and ↓2 denotes a down sampling by a factor of 2. An important property of the DWT is [29]:
Fig. 1. The principle of the DWT decomposition [29]
2.2. Support vector machine (SVM)
This section gives a brief description of SVM. For more details, one can refer to [30], which provides a complete description of the SVM theory.
You can use an SVM when your data has exactly two classes. An SVM classifies data by finding the best hyperplane that separates all data points of one class from those of the other class. The best hyperplane for an SVM means the one with the largest margin between the two classes. Margin means the maximal width of the slab parallel to the hyperplane that has no interior data points. The support vectors are the data points that are closest to the separating hyperplane; these points are on the boundary of the slab. The Fig. 2 illustrates these definitions, with + indicating data points of type 1 and – indicating data points of type –1 [31].
Fig. 2. The schematic model of a linear SVM [32]
2.2.1. Mathematical formulation: primal
The data for training is a set of points (vectors) ${x}_{i}$ along with their categories ${y}_{i}$. For some dimension $d$, ${x}_{i}\u220a{R}^{d}$, and ${y}_{i}=$ ±1. The equation of a hyperplane is $\u2329w,x\u232a+b=$ 0, where $w\u220a{R}^{d}$, $\u2329w,x\u232a$ is the inner product of $w$ and $x$, and $b$ is real. The following problem defines the best separating hyperplane. Find w and b that minimize $\Vert w\Vert $ such that for all data points (${x}_{i},{y}_{i}$), ${y}_{i}\left(\u2329w,{x}_{i}\u232a+b\right)\ge $ 1. The support vectors are ${x}_{i}$ on the boundary, those for which ${y}_{i}\left(\u2329w,{x}_{i}\u232a+b\right)=$ 1. For mathematical convenience, the problem is usually given as the equivalent problem of minimizing $\u2329w,w\u232a/2$. This is a quadratic programming problem. The optimal solution w, b enables classification of a vector z as follows [31]:
Class ($z$) = sign $\left(\u2329w,z\u232a+b\right)$.
2.2.2. Mathematical formulation: dual
It is computationally simpler to solve the dual quadratic programming problem. To obtain the dual, take positive Lagrange multipliers ${\alpha}_{i}$ multiplied by each constraint, and subtract from the objective function [31, 32]:
where you look for a stationary point of ${L}_{P}$ over $w$ and $b$. Setting the gradient of ${L}_{P}$ to 0, you get $w=\sum _{i}{\alpha}_{i}{y}_{i}{x}_{i}$, $\sum _{i}{\alpha}_{i}{y}_{i}=$ 0, and substituting into ${L}_{P}$, you get the dual ${L}_{D}$_{}[31]:
Which you maximize over ${\alpha}_{i}\ge 0$. In general, many ${\alpha}_{i}$ are 0 at the maximum. The nonzero ${\alpha}_{i}$ in the solution to the dual problem defines the hyperplane, as seen in Eq. (4), which gives $w$ as the sum of ${\alpha}_{i}{y}_{i}{x}_{i}$. The data points ${x}_{i}$ corresponding to nonzero ${\alpha}_{i}$ are the support vectors. The derivative of ${L}_{D}$ with respect to a nonzero ${\alpha}_{i}$ is 0 at an optimum. This gives ${y}_{i}\left(\u2329w,{x}_{i}\u232a+b\right)1=$ 0.
In particular, this gives the value of $b$ at the solution, by taking any $i$ with nonzero ${\alpha}_{i}$. The dual is a standard quadratic programming problem.
2.2.3. Nonseparable data
Your data might not allow for a separating hyperplane. In that case, the SVM can use a soft margin, meaning a hyperplane that separates many, but not all data points. There is a standard formulation of soft margins which involves adding slack variables ${s}_{i}$ and a penalty parameter $C$.
The ${L}^{1}$norm problem is:
Such that ${y}_{i}(\u2329w,{x}_{i}\u232a+b)\ge 1{s}_{i},{s}_{i}\ge 0$.
In these formulations, you can see that increasing $C$ places more weight on the slack variables ${s}_{i}$, meaning that the optimization attempts to make a stricter separation between classes. Equivalently, reducing $C$ towards 0 makes misclassification less important.
For easier calculations, consider the ${L}^{1}$ dual problem to this softmargin formulation. Using Lagrange multipliers ${\mu}_{i}$, the function to minimize for the ${L}^{1}$norm problem is [31]:
where you look for a stationary point of ${L}_{P}$ over $w$, $b$ and positive ${s}_{i}$. Setting the gradient of ${L}_{P}$ to 0, you get $b=\sum _{i}{\alpha}_{i}{y}_{i}{x}_{i},\sum _{i}{\alpha}_{i}{y}_{i}=0$, ${\alpha}_{i}=C{\mu}_{i}$. These equations lead directly to the dual formulation [31, 32]:
Subject to the constraints $\sum _{i}{\alpha}_{i}{y}_{i}=\mathrm{0,0}\le {\alpha}_{i}\le C$.
The final set of inequalities, $0\le {\alpha}_{i}\le C$, shows why $C$ is sometimes called a box constraint. $C$ keeps the allowable values of the Lagrange multipliers ${\alpha}_{i}$ in a “box”, a bounded region. The gradient equation for b gives the solution b in terms of the set of nonzero ${\alpha}_{i}$, which corresponds to the support vectors.
2.2.4. Nonlinear transformation with Kernels
Some binary classification problems do not have a simple hyperplane as a useful separating criterion. For those problems, there is a variant of the mathematical approach that retains nearly all the simplicity of an SVM separating hyperplane. This approach uses these results from the theory of reproducing kernels [31]:
There is a class of functions $K\left(x,y\right)$ with the following property. There is a linear space $S$ and a function $\phi $ mapping $x$ to $S$ such that:
The dot product takes place in the space $S$.
This class of functions includes:
Polynomials: For some positive integer $d$:
Radial basis function: For some positive number $\sigma $:
In this study, different parameters including degree of the polynomial ($d$), the RBF kernel and the penalty parameter ($C$) were considered to find the best SVM model.
3. Materials and methods
3.1. Experimental setup
In order to conduct the intended experiments and install the pump and the electrical motor for a better acquisition of data from different faults, the setup of the experiment, as one of the most important parts of data acquisition, was installed according to Fig. 3.
The vibration signals are measured from the pump working under normal condition at a constant rotation speed of 1440 rpm. The accelerometer sensor was attached to the body using magnetic probe (Fig. 4). Vibration data in this study was collected through interruptedsampling technique. To this end, a path was defined using the SpectraPro4 software package, which was then transferred to the data acquisition system, i.e., the Easy_Viber device. This device is equipped with a piezoelectric accelerometer and a tachometer to measure and record the speed.
The considered faults were applied to the impeller of the pump and the mechanical seal, as demonstrated in Fig. 5. In the present study, the following faults were simulated.
3.1.1. Fault in impeller
Two impellers ($d=$ 220 mm) were used. One impeller was assumed to be free from defect. The other impeller was prepared from an outofservice pump. This impeller was rusted and the corrosion existed at the surface and the eye of the impeller as shown in Fig. 5.
Fig. 3. Experimental setup
Fig. 4. Location of the sensor
Fig. 5. a) Faulty seal, b) faulty impeller
a)
b)
3.1.2. Fault in sealing system
Packing’s are one of the oldest sealing devices, named for the way they perform the sealing function. Packing is a ropelike material that is inserted into the cylindrical opening in the rear of the pump casing where the shaft passes through it. This soft material, usually with lubricant in it, consists of a number of rings that wrap around the shaft and compress the inner sides of the stuffing box. This effectively limits leakage out of the pump. The bronze component that separates the two sets of packing is called a lantern ring. The bronze lantern ring serves a number of functions inside the pump packing. The lantern ring is grooved with an oil groove and drilled with holes to allow lubrication to reach the packing material. The lantern ring is also used to distribute cooling water to all packing rings as well as keep the stuffing box clean of containments. Faulty packing was used to reduce the pressure of the pumped liquid and maximize leakage. Sometimes gland packing slightly swells on contact with water. This must be considered and an appropriate clearance should be provided. Otherwise, it would seize the shaft leading to overheating and burning of the packing. In some cases, the shafts have been even broken due to excessively tightened packing. Also, wear in the internal surface of the lantern ring makes friction and heat between the pump shaft and lantern ring, and destroys them. Heat and friction destroys pump shaft and lantern ring. Thus, in this study, we used a burned gland packing and a worn lantern ring, as shown in Fig. 5.
3.1.3. Cavitation
When a pump is under low pressure or high vacuum conditions, suction cavitation occurs. Thus, to simulate the cavitation, the valve at the inlet of the pump is used to make the pressure drop between the suction and the eye of the impeller. The valve control system is used to adjust the flow at the inlet and outlet of the pump. The control of priming and leakage of the pump was done. After closing the delivery valve, the pump started. Then the delivery valve was opened fully and the valve at the suction side was gradually closed. For more information on how to make cavitation, refer to [19].
3.2. Signal processing
As shown in Fig. 6, the collected signals of the first stage are in time domain. Before analyzing the signals, the samples were imported into MATLAB for preprocessing purposes, on each of which then the wavelet transform was applied using the Daubechies mother wavelet at three levels.
Fig. 6. Vibration signals during different pump conditions; a) good condition; b) defect in seal; c) defect in impeller; d) cavitation
a)
b)
c)
d)
3.3. Feature extraction and feature selection
In order to detect the symptoms for each of the faults in the centrifugal pump, the wavelet transform method was used for signal processing.
Vibration signals were processed by the WT at three levels of decomposition. Each signal was decomposed into one approximation signal and three detail signals. These signals are required to be transformed into usable forms; however, due to the large number of data for each fault, their processing procedure using conventional computers is a difficult and timeconsuming task. To overcome the problem, 11 features from each coefficient of the wavelet transform were extracted. Mathematical equations of features are shown in Table 1. These features were extracted from Approximation3 (Ap3), Detail1 (De1), Detail2 (De2) and Detail3 (De3). Since each signal is subdivided into four auxiliary signals, 44 statistical features were extracted from each vibrational signal. However, since not all these features contain beneficial information, superior features should be selected as the input to the classifier. For this purpose, we used the correlationbased feature selection method as one of the most important feature selection techniques available in Weka software.
One hundred samples were collected for each of the pump conditions, 70 percent of which were used for training the SVM, and the remaining 30 percent for testing. The training and testing data were randomly selected, such that the SVM test was conducted using those data never exposed to the classifier.
Table 1. Features and their equations used for feed classifier [11]
$\text{Meanvalue}\text{}\text{\u2192}{F}_{1}=\left(\frac{{\sum}_{n=1}^{N}x\left(n\right)}{N}\right)$

$\text{Standarddeviation}\to {F}_{2}=\sqrt{\left(\frac{n{\sum}_{n=1}^{N}x(n{)}^{2}({F}_{1}{)}^{2}}{n\left(n1\right)}\right)}$

$\text{Samplevariance}\to {F}_{3}=\left(\frac{n{\sum}_{n=1}^{N}x(n{)}^{2}({F}_{1}{)}^{2}}{n\left(n1\right)}\right)$

$\text{Rootmeansquare}\left(\text{RMS}\right)\to {F}_{4}=\sqrt{\frac{{\sum}_{n=1}^{N}\left(x\right(n){)}^{2}}{N}}$

$\text{Crestfactor}\to {F}_{5}=\left(\frac{\mathrm{m}\mathrm{a}\mathrm{x}\left(\leftx\left(n\right)\right\right)}{RMS}\right)$

$\text{Skewness}\to {F}_{6}=\left(\frac{n}{n}1\right)\sum _{n=1}^{N}\left(x\right(n){F}_{1}/{F}_{2}{)}^{3}$

$\text{Slippage}\to {F}_{7}=\left(\frac{1}{N}\right)\left({\sum}_{n=1}^{N}\left(x\right(n){F}_{1}{)}^{3}/{\left(\sqrt{\frac{1}{N}{\sum}_{n=1}^{N}\left(x\right(n){F}_{1}{)}^{2}}\right)}^{3}\right)$

$\text{Kurtosis}\to {F}_{8}=\left\{\frac{n(n+1)}{(n1)(n2)(n3)}\sum _{n=1}^{N}{\left(\frac{x\left(n\right){F}_{1}}{{F}_{2}}\right)}^{4}\right\}\frac{3(n1{)}^{2}}{(n2)(n3)}$

$\text{FifthCentralMoment}\to {F}_{10}=\left(\frac{\sum _{n=1}^{N}\left(x\right(n){F}_{1}{)}^{5}}{N}\right)$

$\text{SixthCentralMoment}\to {F}_{11}=\left(\frac{\sum _{n=1}^{N}\left(x\right(n){F}_{1}{)}^{6}}{N}\right)$

$x\left(n\right)$ is a signal series for $n=$ 1, 2, …, $N$. $N$ is the number of data points

4. Classification results using the SVM
The performance of the SVM depends on different factors including type of kernel function, penalty parameter, kernel parameter, and degree of a polynomial. In this study, onevsone classification was used to classify different conditions of the pump. In order to find the most appropriate model for classification of the faults in the centrifugal pump using the SVM, different parameters including penalty parameter ($C$), degree of a polynomial ($d$) (increasing from 1 until the accuracy of the classifier decreases) and RBF kernel ($0.1\le \sigma \le 1$) were considered. Table 2 shows the classification results of the SVM for different values of the parameters. By changing the penalty parameter, the type of the kernel function and its degree, different classification percentages were obtained. The results revealed that the accuracy of the classifier increased as the penalty coefficient ($C$) decreased, such that the best results were obtained by changing $C=$ 10 to $C=$ 1, through which the accuracy of the SVM classifier increased by 10 percent. Moreover, the accuracy of the classifier increased as the degree of the polynomial increased, such that the maximum accuracy was obtained for a degree of 6. It should be noted that this is not always true, and better performances may be achieved for lower degrees. The most accurate results of fault classification using the SVM were obtained as 66.67 percent for a 6degree kernel at a penalty coefficient of 1.
Fault classification using the SVM with the RBF kernel is given in Table 3 for different values of the penalty coefficient and kernel width. The results indicate that by decreasing the kernel width from 1 to 0.1, the accuracy of the SVM in detecting faults significantly increases, such that the best accuracy was obtained for a kernel width of 0.1. Similar to the previous results, the accuracy of the classifier increased with decreases in penalty coefficient. The best accuracy for fault classification using the SVM with the RBF kernel was obtained as 96.67 percent for a penalty coefficient of 1 and a kernel width of 0.1.
Table 2. The accuracy of fault classification for a pump using a SVM with polynomial kernel
The penalty coefficient
The kernel degree of the polynomial ($d$)

Classification accuracy (%)


$C=$ 10

$C=$ 5

$C=$ 1


1

50.00

50.00

53.33

2

50.00

53.33

56.67

3

53.33

56.67

56.67

4

53.33

60.00

60.00

5

56.67

63.33

63.33

6

56.67

53.33

66.67

7

50.00

53.33

56.67

8

50.00

53.33

56.67

Table 3. The accuracy of fault classification for a pump using a SVM RBF kernel
The penalty coefficient
RBF kernel Sigma parameter ($\sigma $)

Classification accuracy (%)


$C=$ 10

$C=$ 5

$C=$ 1


0.1

93.33

93.33

96.67

0.2

93.33

93.33

93.33

0.3

90.00

93.33

93.33

0.4

90.00

93.33

93.33

0.5

86.67

86.67

90.00

0.6

83.33

86.67

90.00

0.7

83.33

86.67

86.67

0.8

77.67

83.33

86.67

0.9

77.67

77.67

86.67

1

75.13

76.33

77.67

Table 4. Comparisons between the current work and some published work
References

Fault types

Classifier type

Classified
states

Maximum classification
efficiency (%)

Denoising
technique

Feature selection method

Yunlong and Peng [22]

IA, U, L

SVM

4

92.60

EMD

LS

Sakthivel et al. [23]

BF, IF, SF, IBF, Ca

C4.5 decision tree algorithm

6

100

–

–

Farokhzad et al. [11]

IF, SF, Ca

ANN

4

100

FFT

–

Zouari et al. [33]

IA, Ca, PF, AI

ANN, FNN

5

99.77

–

PCA

Muralidharan and Sugumaran [19]

BF, Ca, IF, IBF

J48 algorithm

5

99.76

WT

–

Ebrahimi (Present work)

SF, IF, Ca

SVM

4

96.67

WT

CFS

In order to illustrate the potential application of the proposed methodology in centrifugal pump fault diagnosis, a comparative study between the present work and published literature is presented in Table 4 [11, 19, 22, 23, 33]. The comparing items include the machine elements used, fault type, classifier type, classified states, maximum classification efficiencies, denoising technique and feature selection method. Therefore, the comparison results demonstrate that the proposed approach is effective and reliable in detecting centrifugal pump faults.
Also, the results of this study were in agreement with those of Konar [27], Yin and Hu [34] and, Bordoloi and Tiwari [35]. Comparison of the polynomial kernel method and the RBF kernel indicates that the RBF kernel has more accuracy for faults classification.
5. Conclusions
The aim of the present study was to use the wavelet transform and SVM to detect faults in centrifugal pumps. Five classical states viz., good condition, defect in seal, defect in impeller and cavitation, were simulated on the centrifugal pump. The vibration signals were decomposed by the Daubechies wavelets, and statistical features were extracted from the detail coefficients and approximation coefficients of the wavelet. Superior features should be selected as the input to the classifier. For this purpose, we used the correlationbased feature selection method as one of the most important feature selection techniques. Various parameters were used to increase the accuracy of the classifier, the effect of each of which is given in the following:
1) The accuracy of the classifier increased as the degree of kernel increased. However, this is not always true and lower degrees may result in better performances.
2) The accuracy of the classifier increased as the penalty coefficient ($C$) decreased, such that the best results were obtained by changing $C=$ 10 to $C=$ 1, through which the accuracy of the SVM classifier increased by 10 percent.
3) Reducing the width of the RBF kernel results in further accuracy of the classifier.
The best accuracy for fault classification using the SVM and RBF kernel was obtained as 96.67 percent for a penalty coefficient of 1 and a kernel width of 0.1. The results demonstrate that the proposed method can achieve good performance for fault detection of the centrifugal pump.
Acknowledgements
This research is supported by the Islamic Azad University, Kermanshah Branch, Iran. This assistance is gratefully appreciated.
References
 Yang M., Makis V. ARX modelbased gearbox fault detection and localization under varying load conditions. Sound and Vibration, Vol. 329, Issue 24, 2010, p. 52095221. [Publisher]
 Miao Q., Makis V. Condition monitoring and classification of rotating machinery using wavelets and hidden Markov models. Mechanical Systems and Signal Processing, Vol. 21, Issue 2, 2007, p. 840855. [Publisher]
 Niu G., Han T., Yang B. S., Tan A. C. C. Multiagent decision fusion for motor fault diagnosis. Mechanical Systems and Signal Processing, Vol. 21, Issue 3, 2007, p. 12851299. [Publisher]
 Wang X., Makis V., Yang M. A wavelet approach to fault diagnosis of a gearbox under varying load conditions. Sound and Vibration, Vol. 329, Issue 9, 2010, p. 15701580. [Publisher]
 Sun H. C., Huang U. C. Support vector machine for vibration fault classification of steam turbinegenerator sets. Procedia Engineering, Vol. 24, 2011, p. 3842. [Publisher]
 Bagheri B., Ahmadi H., Labbafi R. Implementing discrete wavelet transform and artificial neural networks for acoustic condition monitoring of gearbox. Elixir Mechanical Engineering, Vol. 35, 2011, p. 29092911. [Search CrossRef]
 Van Hecke B., Yoon J., He D. Low speed bearing fault diagnosis using acoustic emission sensors. Applied Acoustics, Vol. 105, 2016, p. 3544. [Publisher]
 Wu W., Lin T. R., Tan A. C. C. Normalization and source separation of acoustic emission signals for condition monitoring and fault detection of multicylinder diesel engines. Mechanical Systems and Signal Processing, Vols. 6465, 2015, p. 479497. [Publisher]
 Yuan S., Chu F. Fault diagnosis based on support vector machines with parameter optimisation by artificial immunisation algorithm. Mechanical Systems and Signal Processing, Vol. 21, Issue 3, 2007, p. 13181330. [Publisher]
 Wang J. D., Gao Y. Q., Zhao H. Y., Cong R. Sectional normalization and recognization on the PV Diagram of reciprocating compressor. Journal of Measurements in Engineering, Vol. 3, Issue 2, 2015, pp. 35 41. [Search CrossRef]
 Farokhzad S., Ahmadi H., Jaefari A., Asadi Asad Abad M. R., Ranjbar Kohan M. Artificial neural network based classification of faults in centrifugal water pump. Journal of Vibroengineering, Vol. 14, Issue 4, 2012, p. 17341744. [Search CrossRef]
 Yayli M. O. Free vibration behavior of a gradient elastic beam with varying cross section. Shock and Vibration, Vol. 2014, 2014, p. 801696. [Publisher]
 Yayli M. O. On the axial vibration of carbon nanotubes with different boundary conditions. Micro and Nano Letters, Vol. 9, Issue 11, 2014, p. 807811. [Publisher]
 Hosseini A. H., Rahmani O., Nikmehr M., Golpayegani I. F. Axial vibration of cracked nanorods embedded in elastic foundation based on a nonlocal elasticity Model. Sensor Letters, Vol. 14, Issue 10, 2016, p. 10191025. [Publisher]
 Yayli M. O. An analytical solution for free vibrations of a cantilever nanobeam with a spring mass system. International Journal of Engineering and Applied Sciences, Vol. 7, Issue 4, 2015, p. 1018. [Search CrossRef]
 Li Y., Xu M., Wang R., Huang W. A fault diagnosis scheme for rolling bearing based on local mean decomposition and improved multiscale fuzzy entropy. Journal of Sound and Vibration, Vol. 360, 2016, p. 277299. [Publisher]
 Li Y., Xu M., Wei Y. An improvement EMD method based on the optimized rational Hermite interpolation approach and its application to gear fault Diagnosis. Measurement, Vol. 63, 2015, p. 330345. [Publisher]
 Li Y., Xu M., Wei Y., Huang W. A new rolling bearing fault diagnosis method based on multiscale permutation entropy and improved support vector machine based binary tree. Measurement, Vol. 77, 2016, p. 8094. [Publisher]
 Muralidharan V., Sugumaran V. Feature extraction using wavelets and classiﬁcation through decision tree algorithm for fault diagnosis of monoblock centrifugal pump. Measurement, Vol. 46, Issue 1, 2013, p. 353359. [Publisher]
 Kiran V., Hemantha K., Gangadharan K. V. Gear fault detection using vibration analysis and continuous wavelet transform. Procedia Materials Science, Vol. 5, 2014, p. 18461852. [Publisher]
 Li Y., Yang Y., Li G., Xu M., Huang W. A fault diagnosis scheme for planetary gearboxes using modified multiscale symbolic dynamic entropy and mRMR feature selection. Mechanical Systems and Signal Processing, Vol. 91, 2017, p. 295312. [Publisher]
 Yunlong Z., Peng Z. Vibration fault diagnosis method of centrifugal pump based on EMD complexity feature and least square support vector machine. Energy Procedia, Vol. 17, 2012, p. 939495. [Publisher]
 Sakthivel N. R., Sugumaran V., Babudevasenapati S. Vibration based fault diagnosis of monoblock centrifugal pump using decision tree. Expert Systems with Applications, Vol. 37, Issue 6, 2010, p. 40404049. [Publisher]
 Cui H., Zhang L., Kang R., Lan X. Research on fault diagnosis for reciprocating compressor valve using information entropy and SVM method. Journal of Loss Prevention in the Process Industries, Vol. 22, 2009, p. 864867. [Publisher]
 Zhao G., Jiang D., Diao J. Application of wavelet timefrequency analysis on fault diagnosis for steam turbine. 5th International Conference Acoustical and Vibratory Surveillance Methods and Diagnostic Techniques, Senlis, France, 2004. [Search CrossRef]
 Jedlinski L., Jonak J. Early fault detection in gearboxes based on support vector machines and multilayer perceptron with a continuous wavelet transform. Applied Soft Computing, Vol. 30, 2015, p. 636641. [Publisher]
 Konar P., Chattopadhyay P. Bearing fault detection of induction motor using wavelet and support vector machines (SVMs). Applied Soft Computing, Vol. 11, Issue 6, 2011, p. 42034211. [Publisher]
 Choudhury M., Ganguly A. Transmission line fault classification using discrete wavelet transform. International Conference on Energy, Power and Environment: Towards Sustainable Growth (ICEPE), 2015. [Publisher]
 Bendjama H., Bouhouche S., Boucherit M. S. Application of wavelet transform for fault diagnosis in rotating machinery. International Journal of Machine Learning and Computing, Vol. 2, Issue 1, 2012, p. 8287. [Publisher]
 Cristianini N., Taylor J. S. An Introduction to Support Vector Machines and Other KernelBased Learning Methods. Cambridge University Press, 2000. [Publisher]
 Bioinformatics Toolbox, Getting Started, Support Vector Machine (SVM). Matlab Help, 2011. [Search CrossRef]
 Shen L., Chen H., Yu Z., Kang W., Zhang B., Li H., Yang B., Liu D. Evolving support vector machines using fruit fly optimization for medical data classification. KnowledgeBased Systems, Vol. 96, Issue 15, 2016, p. 6175. [Publisher]
 Zouari R., SiegZieba S., Sidahmed M. Fault detection system for centrifugal pumps using neural networks and neurofuzzy techniques. Surveillance 5 CETIM, Senlis, 2004. [Search CrossRef]
 Yin Z., Hou J. Recent advances on SVM based fault diagnosis and process monitoring in complicated industrial processes. Neurocomputing, Vol. 174, 2016, p. 643650. [Publisher]
 Bordoloi D. J., Tiwari R. Optimum multifault classification of gears with integration of evolutionary and SVM algorithms. Mechanism and Machine Theory, Vol. 73, 2014, p. 4960. [Publisher]
Articles Citing this One
IEEE Access
Zahoor Ahmad, Akhand Rai, Andrei S. Maliuk, JongMyon Kim

2020

IEEE Access
Zahoor Ahmad, Alexander E. Prosvirin, Jaeyoung Kim, JongMyon Kim

2020
