Sage Journals: Discover world-class research

Abstract

Rotor unbalance faults are one of the high-frequency faults in rotating machinery. As such, their accurate and timely diagnosis is important. In contrast to traditional methods based on static features, the dynamics features and support vector machines (SVM) are combined for the accurate detection and classification of rotor unbalance faults. First, the dynamical trajectories of the rotor system associated with unbalance faults are accurately identified locally based on the deterministic learning theory, which is more sensitive to abnormal changes in the rotor system. Second, entropy dynamics features, including the sample entropy, fuzzy entropy, and permutation entropy, are extracted based on the obtained dynamical trajectory data. Finally, the dynamics features are used to train the fault classifier based on the SVM with a Gaussian kernel function. Experiments on a rotor unbalance fault test rig demonstrate the effectiveness of the proposed method. The accurate detection and classification of rotor unbalance faults were also achieved compared with the results of employing static time or frequency features.

Keywords

Fault detection rotor unbalance dynamics features deterministic learning support vector machine

Introduction

As the key component of rotating machinery, rotor systems are widely used in electric power generation, coal mining, ship, and aerospace equipment. One of the most common faults in rotating machinery is rotor unbalance, which has a significant impact on the efficient and stable operation of rotor systems.¹ Unbalance faults may originate from the rotor manufacturing process, which lead to secondary failures of the rotor system, such as the rub-impact phenomenon and broken rotor bars.² Therefore, the accurate detection and classification of rotor unbalance faults is a crucial step in the safe operation of rotating machinery.

Rotor unbalance fault diagnosis has attracted the attention of researchers in recent years. Fault-specific symptom parameters are key elements in rotor fault diagnosis.¹ Several time-frequency feature extraction methods have been proposed for rotor unbalance fault symptom analysis.^3–5 Zhou et al.⁴ used empirical mode decomposition (EMD) to extract the intrinsic mode to represent the different harmonics. Based on the idea of EMD, Chen et al.⁵ proposed the variational mode decomposition method to diagnose the unbalance faults of wind power gearboxes. However, the diagnosis results cannot be obtained based solely on time-domain or frequency-domain features. The quantitative diagnosis results still need to be evaluated by analysts based on their experience or by other methods.^6,7 To resolve this issue, machine learning methods combined with fault features have been proposed,^8–10 such as support vector machines (SVMs),^11,12 K-nearest neighbor,¹³ and deep learning.¹⁴ The main processes of these methods are feature selection and fault classification. Among these machine learning methods, SVM has the advantage of high classification performance, and has been widely used in rotor fault diagnosis.^15–17 However, the performance of SVM methods is dependent on the fault feature extraction process.¹ The existing time or frequency domain features are static features, while the nonlinear features, for example, sample entropy (SE), fuzzy entropy (FE), and permutation entropy (PE), are obtained based on the status data of rotor systems.¹⁸ Compared with the static features, system dynamics indices are more sensitive to system changes, and can provide a more accurate representation of the original system by including internal dynamics information.¹⁹ The key step in dynamics index extraction is dependent on the locally accurate identification of the system dynamical trajectory based on deterministic learning, which was proposed by Wang and Hill.²⁰ The main advantage of deterministic learning is that the system dynamical trajectory can be approximated by radial basis function (RBF) neural networks, in which the neural weights converge in a stable manner. Therefore, the identified trajectory in a dynamical system can be stored in constant neural networks for further feature extraction. The dynamical system trajectory has been applied in the fault diagnosis and pattern recognition of small abrupt and incipient faults, electrocardiogram signals, and the human gait.^21–23

Inspired by the above methods, the deterministic learning theory and SVM method have been used to detect rotor unbalance faults. First, the dynamical trajectory of rotor systems with different unbalance faults are identified by deterministic learning, while the SE, FE, and PE are used to extract the dynamic features based on the obtained dynamical trajectory. Second, the corresponding dynamic features are input into the SVM with a Gaussian kernel (Gaussian-SVM) to train the classification model. Finally, the unbalance faults are detected and classified based on the training SVM model. The results of the rotor unbalance fault test rig demonstrate the advantages of the proposed method. The main contribution of this study is the improved sensitivity diagnosis based on the deterministic learning theory and Gaussian-SVM method.

The remainder of this paper is organized as follows. Section 2 presents the problem statement. Section 3 describes the experiments and methods employed in this study. The data analysis and diagnosis results are provided in Section 4. Section 5 presents some discussions. Section 6 presents the conclusions and directions for future study.

Problem formulation

A schematic diagram of the rotor unbalance is shown in Figure 1. The black shaft and rotor disk in the figure represent the normal state of the rotor system (the centroid is $O_{1}$ ). The red dotted line shows the unbalance fault states (the centroid is $O$ ). The unbalance parameter $u$ reflects the degree of the unbalance fault in the rotor system.

Figure 1.

The schematic diagram of the rotor unbalance.

The mathematical model of the rotor has been investigated in numerous studies, which can be described in the following form.^7,24

\dot{X} = f (x, y, u) + v (x, y, u)

(1)

where $X = [x, y]$ represents the vibration signal in the x- and y-directions, $u$ represents the unbalance parameter, and different $u$ values will produce different degrees of unbalance faults (i.e. the unbalance parameter $u$ in Figure 1). $f (x, y, u)$ represents the unknown dynamics, and $v (x, y, u)$ represents the system uncertainties and small disturbances.

The unbalance parameter $u$ changes the centroid of the rotor system, and influences the system dynamics behavior. Various unknown system dynamics behaviors can cause the operating state of the rotor system to become uncontrollable. If the unbalance faults cannot be detected accurately, they could produce secondary faults, such as serious rub-impact or shaft bending faults. Therefore, the main goal of this paper is to propose a new fault detection method for the accurate detection and classification of the different unbalance faults.

Assumption 1: The rotor vibration states x and y can be measured by the sensors, and are still bounded after the unbalance force occurs.

Materials and methods

Rotor unbalance experiments

The rotor unbalance test rig is shown in Figure 2. The rotor unbalance device consists of a disk and counterweights. The vibrations of the test rig in the $x$ - and $y$ -directions were sampled by two eddy current sensors. Different counterweights produced different degrees of unbalance faults in the rotor system. Based on the rotor dynamical differential equation function (1), different values of the unbalance parameter $u$ will produced different unbalance faults. Therefore, in the experiment on unbalance faults, the counterweights (0, 2.5, 3.6, and 5.0 g) were taken as the unbalance parameter and installed on the disk to inject different degrees of unbalance faults, respectively. The other main components of the test rig are the rotor motor, signal conditioner, and oscilloscope.

Figure 2.

Rotor unbalance test rig.

Based on the selected counterweights, the corresponding data samples were divided into four classes of unbalance faults: normal class (0 g counterweight), slight unbalance fault class (2.5 g counterweight), moderate unbalance fault class (3.6 g counterweight), and severe unbalance fault class (5.0 g counterweight). The sample frequency was set as 5.12 kHz, and the sampling time was 2 s. In this study, two unbalance fault datasets (a small dataset #1 and big dataset #2) were constructed as follows. For dataset #1, the sampling speed was set as 1200 rpm. The normal class includes 30 sample data obtained with the benchmark test rig. The slight unbalance fault class also includes 30 sample data associated with the 2.5 g counterweight. The other two classes have 240 sample data (each class contained 120 data associated with the 3.6 and 5.0 g counterweights).

For dataset #2, the unbalance data was collected under three different rotor speeds: 1200, 1600, and 1800 rpm. The total number of data samples is 2100. Detailed descriptions of dataset #1 and dataset #2 are provided in Table 1. To evaluate the proposed method, the training and testing datasets were randomly divided 1:1 for dataset #1 and 6:4 for dataset #2.

Table 1.

Descriptions of dataset #1 and dataset #2.

Fault classes		Normal	Slight unbalance	Moderate unbalance	Severe unbalance	Training dataset	Testing dataset
Counterweights		0 g	2.5 g	3.6 g	5.0 g	Training dataset	Testing dataset
Dataset #1	1200 rpm	30	30	120	120	150	150
Dataset #2	1200 rpm	100	100	200	200	1260	840
	1600 rpm	100	100	200	200
	1800 rpm	100	100	200	200

Deterministic learning

To diagnose the unbalance faults, the nonlinear dynamical ordinary differential equation (1) was first considered to identify the dynamical trajectory. Under various conditions, the system uncertainties and unknown dynamics could not be clearly decoupled; therefore, these two factors were considered as a whole term to be identified.²¹ The following RBF neural network was constructed to estimate the system dynamics.

\dot{\hat{X}} = - A (\hat{X} - X) + {\hat{W}}^{T} S (X, u)

(2)

where $A = diag (a_{x}, a_{y})$ represents the designed constant gain parameter; $\hat{W}$ is the weight of the RBF neural networks; and $S (X, u)$ is the RBF. The weight updating law is set as

\overset{\cdot}{\hat{W}} = \overset{\cdot}{\tilde{W}} = - Γ S (X, u) \tilde{X} - σ Γ \hat{W}

(3)

where $\tilde{X} = X - \hat{X}$ , $\tilde{W} = \hat{W} - W^{*}$ , $W^{*}$ is the ideal weight vector, and $Γ = Γ^{T} > 0$ , and $σ$ is a small value. Based on the dynamics identifier (2) and the weight updating law (3), the unknown rotor dynamics can be accurately identified locally.

Dn = {\hat{W}}^{T} S (X, u) + ξ_{1}

(4)

where $Dn = diag (D n_{x}, D n_{y})$ represents the rotor dynamics, and $ξ_{1}$ is the practical identification error. To obtain and store the dynamical knowledge in constant RBF networks, the neural weights are chosen as the mean value after the stable convergence process.

\bar{Dn} = {\bar{W}}^{T} S (X, u) + ξ_{2}

(5)

where $\bar{W} = mea n_{t \in [t_{a}, t_{b}]} \hat{W} (t)$ , $ξ_{2} = O (ξ_{1})$ , $[t_{a}, t_{b}]$ is the selected time segment between the stable convergence process. Therefore, the system dynamical trajectory can be obtained and stored as $\bar{Dn} = diag ({\bar{Dn}}_{x}, {\bar{Dn}}_{y})$ for the next dynamics features extraction.

Remark 1: A recent study⁷ presented the detailed identification process of the nonlinear system dynamics, which is omitted here for unbalance faults and presented in the Discussion section.

Dynamics features

Based on the obtained dynamical trajectory data $({\bar{Dn}}_{x}, {\bar{Dn}}_{y})$ , the nonlinear features, including the sample entropy (SE), fuzzy entropy (FE), and permutation entropy (PE) as shown in Li et al.,¹⁸ were used to extract the dynamics features of the unbalance fault vibration signals. The SE can measure the regularity of the obtained dynamical trajectory data. The FE and PE can characterize the randomness and complexity of the dynamical trajectory data, respectively. The different SE, FE, and PE values indicate the degree of data irregularity of the corresponding unbalance faults. The greater the data irregularity, the larger the SE, FE, and PE values. The corresponding value of SE, FE, and PE values based on the dynamical trajectory data $({\bar{Dn}}_{x}, {\bar{Dn}}_{y})$ are defined as $S E_{x}, S E_{y}, F E_{x}, F E_{y}, P E_{x}, P E_{y}$ . The subscripts $x$ and $y$ represent the entropy dynamics features in the x- and y-directions, respectively.

Remark 2 : In this study, the dynamics features were extracted based on the dynamical trajectory, which can be locally identified via deterministic learning theory. In contrast, the static features are the features extracted from the system vibration states directly, such as the time or frequency domain features.

Support vector machine

The SVM has been widely used in the field of diagnosis owing to its good optimization and generalization ability, and short training time.^15–17 The basic idea of the SVM is to construct a hyperplane for the optimal separation of different classes of input data. It is suitable for nonlinear and high-dimensional conditions.¹¹ Based on the contents of the SVM, the hyperplane can be described as ${(ω \cdot x) + b = 0 | x ϵ H, ω ϵ H, b ϵ R$ , where $x$ is the input data vector, $ω$ is a vector, and $H$ represents the Hilbert space. To obtain the optimal separation of the input training data with the maximal margin, the main task of the SVM is to find the optimal hyperplane (i.e. $\min \frac{1}{2} | | ω | |^{2} + C \sum_{N} ϑ$ , where $C > 0$ is the penalty parameter, $ϑ \geq 0$ is a non-negative slack variable and N is the total number of training samples) under the different classification conditions.¹²

For the obtained rotor fault dynamics feature training dataset ${S E_{xi}, S E_{yi}, F E_{xi}, F E_{yi}, P E_{xi}, P E_{yi}, z_{i}}$ , $i = 1, \dots, N$ and $z_{i} \in {1, 2, 3, 4}$ (associated with the four classes of unbalance). The construction of a suitable hyperplane is dependent on the linear or nonlinear properties of the input training data. In the rotor fault diagnosis process, because the signals of the rotor system are nonlinear and nonstationary, the SVM should have the ability to deal with the nonlinear data and multiclass classification. To address the nonlinear problems, a kernel function is used to construct the decision function of the optimal hyperplane.

f (x) = sgn (\sum_{i} a_{i} z_{i} k (x_{i}, x) + b)

(6)

where $k (x_{i}, x)$ is the kernel function, $x_{i} = [S E_{xi}, S E_{yi}, F E_{xi}, F E_{yi}, P E_{xi}, P E_{yi}]$ and $a_{i} \geq 0$ is the Lagrange multiplier. Among the kernel functions of the SVM in existing literature (e.g. linear, polynomial, sigmoid, and Gaussian), the Gaussian kernel function is the most commonly used for the wide convergence domain of such functions.¹⁷ $σ$ is the acceptance threshold parameter of the kernel function.

k (x_{i}, x) = e^{- \frac{| | x_{i} - x | |^{2}}{2 σ^{2}}}

(7)

The Gaussian kernel function was input into the SVM (Gaussian-SVM) for the diagnosis of rotor unbalance faults. Based on the Gaussian function, the selection of parameters is highly dependent on the accuracy of classification because the kernel function is directly related to the optimal hyperplane. In this study, the determination of the upper parameter of the Gaussian kernel function is utilized based on the rid-quadtree model selection method, which was proposed by Beltrami and da Silva.¹⁷

To address the multiclass classification problem, the one-against-all (OAA) method was applied to deal with the diagnosis of the four classes of unbalance faults. The reason is because the OAA-SVM is the least time-consuming among the multiclass classification methods (e.g. one-against-one SVM).¹⁵ The detailed OAA-SVM process is as follows.

• For the four classes of rotor unbalance faults, the decision function of the optimal hyperplane is set as follows.

f^{p} (x) = sgn (\sum_{i} a_{i}^{p} z_{i} k (x_{p}, x) + b^{p})

(8)

where $p$ belongs to {1,2,3,4}, the $p$ th class is taken as the positive class, and the other three classes are the negative classes.

• The input fault data $x_{p}$ considered to the $p$ th class, where $p$ is the maximum label among the $\sum_{i} a_{i}^{p} z_{i} k (x_{p}, x) + b^{p}$ .

Based on the selected kernel function (Gaussian-SVM) and the multiclass classifier (OAA-SVM), the accurate diagnosis of the rotor system with four classes unbalance faults was obtained in the following experimental section.

In this study, the deterministic learning and the Gaussian-SVM were combined to detection the rotor unbalance faults more accurately. The proposed method included two phases, which are the dynamics feature extraction phase and the diagnosis phase. The following Figure 3 presents a flowchart of the proposed method, and the algorithm of the proposed method is shown in Table 2. The obtained experimental unbalance data were input into the dynamical RBF neural networks, and the dynamical trajectories of different unbalance fault classes were stored. The nonlinear features (SE, FE, and PE) were extracted based on the dynamical trajectory. Finally, the dynamics features were used to train the classification model based on the Gaussian-SVM.

Figure 3.

Flowchart of the proposed method.

Table 2.

Algorithm of the proposed method.

Dynamics feature extraction phase

① Rotor unbalance data acquisition:

X = [x, y]

② RBF neural network identifiers were constructed as shown in (2).

③ Dynamical trajectories were obtained based on the convergence of neural weights

\bar{Dn} = {\bar{W}}^{T} S (X, u)

④ Dynamics features of

\bar{Dn}

were calculated,

S E_{x}, S E_{y}, F E_{x}, F E_{y}, P E_{x}, P E_{y}

Diagnosis phase

① Classifier was constructed based on the obtained dynamics features of

\bar{Dn}

② The parameters of classification model were determined by using the grid-quadtree selection method Beltrami and da Silva.¹⁷

③ Repeating step 1–4 as shown in the dynamics feature extraction phase for the test data.

④ Obtaining the diagnosis results after inputting the dynamics features of test data to the constructed classifier.

Diagnosis results

Firstly, dataset #1 was selected to represent the detailed diagnosis process. The time-domain waveforms of the four classes of unbalance faults are shown in Figure 4. Figure 4(a) to (d) were related to the normal class, the slight unbalance class, the moderate unbalance class and the severe unbalance class.

Figure 4.

Time-domain wave forms of four classes of unbalance faults: (a) normal, (b) slight unbalance fault, (c) moderate unbalance fault, and (d) severe unbalance fault.

The corresponding frequency and marginal spectrum are shown in Figures 5 and 6, respectively. The time-domain features (e.g. extremum, mean value, variance, and kurtosis) and frequency features (e.g. amplitude-frequency, phase-frequency, and peak difference) which were extracted from the rotor system state are noted as the static features for the further classification comparison. As seen in Figures 5 and 6, the main frequency components of the normal class are concentrated on the fundamental frequency (1X) and second harmonic frequency (2X). The third harmonic frequency (3X) is included in the main frequency components of the slight unbalance class. The third harmonic frequency (3X) and second frequency (2X) are separated more clearly in the moderate unbalance class. The severe unbalance class had a weak quadruple frequency.

Figure 5.

Frequency spectrum of four classes of unbalance faults: (a) normal, (b) slight unbalance fault, (c) moderate unbalance fault, and (d) severe unbalance fault.

Figure 6.

Marginal spectrum of the four classes of unbalance faults: (a) normal, (b) slight unbalance fault, (c) moderate unbalance fault, and (d) severe unbalance fault.

Secondly, the deterministic learning theory was used to identify the rotor system dynamics. The dynamical trajectory was obtained from equation (5). The parameters of RBF neural networks, which are used for deterministic learning, were set as: $a_{x} = a_{y} = 0.02$ , $Γ = 1.98$ , $σ = 2$ , and ${\hat{W}}_{x} (0) = {\hat{W}}_{y} (0) = 0.0$ . All the data were normalized for the further detection and classificat ion. Therefore, the neuron centers were evenly spaced on [−1.09,1.09] × [−1.09,1.09] with the widths of 0.05, and the distribution of RBF neural networks is constructed based on the neuron centers width. The identification results of the different unbalance faults are shown in Figures 7 to 10, where (a) shows the convergence of the neural network weights, and (b) shows the identified dynamical trajectories of the four classes of unbalance faults.

Figure 7.

Dynamical trajectory identification of normal condition. (a) convergence of NN weights and (b) dynamics trajectory.

Figure 8.

Dynamical trajectory identification of slight unbalance fault. (a) convergence of NN weights and (b) dynamics trajectory.

Figure 9.

Dynamical trajectory identification of moderate unbalance fault. (a) convergence of NN weights and (b) dynamics trajectory.

Figure 10.

Dynamical trajectory identification of severe unbalance fault. (a) convergence of NN weights and (b) dynamics trajectory.

Compared with the time-frequency features in Figures 4 to 6, it shows that the dynamical trajectory shapes of the different classes of unbalance faults were clearly distinguished from the identification results. Based on the Figures 7(b) and 10(b), the identified dynamical trajectories can be stored as the constant data owing to the convergence of the neural weights (Figures 7(a) and 10(a)). Finally, based on the obtained dynamical trajectory data, the dynamics features, including $S E_{x}, S E_{y}, F E_{x}, F E_{y}, P E_{x}, P E_{y}$ , were directly extracted. The total number of the extracted features are 150 × 6 = 900 for the training dataset #1. Finally, the extracted dynamics features were input into the Gaussian-SVM to produce the training classification model. To evaluate the effectiveness of the proposed method, four SVM classifiers were constructed, and the diagnosis results were obtained for the SVM classifier 1 (associated with time-domain features), the SVM classifier 2 (associated with frequency-domain features), the SVM classifier 3 (associated with dynamics features), and the SVM classifier 4 (associated with the hybrid features which are combined with the static features and dynamics features), respectively. The parameters (penalty parameter $C$ and acceptance threshold parameter $σ$ ) of the four Gaussian-SVM classifiers were determined based on the grid-quadtree selection method,.¹⁷ The values of the four classifiers are listed in Table 3. The parameter determination process is described in the Discussion section.

Table 3.

The parameters of four Gaussian-SVM classifiers.

	SVM classifier 1	SVM classifier 2	SVM classifier 3	SVM classifier 4
Best $C$	48.36	79.53	64	59.8
Best $σ$	14.14	10.21	17.15	15.08

For test dataset #1, based on the four training classifiers, the detailed classification results and their corresponding confusion matrices are shown in Figures 11 and 12. The figures show that the normal class, slight unbalance fault class, moderate unbalance fault class, and severe class can be accurately classified using the proposed method. It is also shown that the SVM with dynamics features has a higher classification rate (Acc = 97.33%) than the SVM with time-domain features (Acc = 90.67%) or frequency features (Acc = 91.33%). In particular, the detection of rotor unbalance faults reached 100% using SVM classifier 3, with a false alarm rate (FAR) of zero. The detailed comparison of the detection and classification results is shown in Table 4. The static and dynamics features were combined for training SVM classifier 4. The fault detection accuracy was also 100%, while the classification results reached 98%. The similar diagnosis results of SVM classifiers 3 and 4 indicate that the dynamics features have the greatest impact on the classification results and are more sensitive to abnormal unbalance changes in the rotor system.

Remark 3: TP (true positive) represents the number of accurate unbalance fault detections. FP (false positive) represents the number of incorrect detections. TN (true negative) represents the number of accurate normal-state detections, while FN (false negative) represents the opposite. TC (true classification) represents the number of accurate unbalance fault classifications. FC (false classification) represents the number of incorrect classifications. The detection and classification accuracies are $(TF + TN) / (TP + TN + FP + FN)$ and $TC / (TC + FC)$ , respectively.

Figure 11.

The diagnosis results of rotor unbalance faults based on dataset #1. (a) based on the SVM classifier 1, (b) based on the SVM classifier 2, (c) based on the SVM classifier 3, and (d) based on the SVM classifier 4.

Figure 12.

The confusion matrices for the four classifiers based on dataset #1. (a) based on the SVM classifier, (b) based on the SVM classifier 2, (c) based on the SVM classifier 3, and (d) based on the SVM classifier 4.

Table 4.

Comparison results of detection and classification based on the dataset #1.

Classifiers	TP	TN	FP	TC	FC	Accuracy of detection (%)	Accuracy of classification (%)
SVM Classifier 1 (time features)	140	5	5	136	14	96.67	90.67
SVM Classifier 2 (frequency features)	140	5	5	137	13	96.67	91.33
SVM Classifier 3 (dynamics features)	140	10	0	146	4	100	97.33
SVM Classifier 4 (hybrid features)	140	10	0	147	3	100	98

Dataset #2 was also used to verify the effectiveness of the proposed method. The diagnosis confusion matrices are shown in Figure 13. The detailed comparison of the detection and classification results is presented in Table 5. It can be seen that the classification results of dataset #2 are consistent with the results of dataset #1. Most importantly, based on the dynamics or hybrid features, the FAR valued (FAR = 2/105 = 1.9%) are still much lower than the FAR of SVM classifier 1 (FAR = 15/105 = 14.29%) and classifier 2 (FAR = 8/105 = 7.62%). The results of datasets #1 and #2 show that the diagnostic performance of the proposed method has advantages in the case of small samples, which are more suitable for practical applications.

Figure 13.

The confusion matrices for the four classifiers based on dataset #2. (a) based on the SVM classifier, (b) based on the SVM classifier 2, (c) based on the SVM classifier 3, and (d) based on the SVM classifier 4.

Table 5.

The comparison of detection and classification results based on the dataset #2.

Classifiers	TP	TN	FP	FN	TC	FC	Accuracy of detection (%)	Accuracy of classification (%)
SVM classifier1 (time features)	729	90	15	6	783	57	97.5	93.21
SVM classifier2 (frequency features)	732	97	6	3	793	47	98.69	94.4
SVM classifier3 (dynamics features)	735	103	2	0	822	18	99.76	97.86
SVM classifier4 (hybrid features)	735	103	2	0	825	15	99.76	98.21

Discussion

Discussion 1: During the construction process of the training classification model, it should be noted that the penalty parameter $C$ and the acceptance threshold parameter $σ$ of the Gaussian kernel function are the main focus of researchers. The final diagnosis results rely heavily on these two parameters because the penalty parameter and the acceptance threshold can determine the good region (in contrast to the underfitting and overfitting region). Many methods are available to deal with the selection of the parameters ( $C$ , $σ$ ). In this study, the grid-quadtree selection method proposed by Beltrami and da Silva¹⁷ was employed for parameter selection. In the grid-quadtree method, a quadtree is used to identify uniform regions in the hyperparameter space. This process can be divided into three steps. After inputting the training dataset, the first step is to determine the initial solution. The second step is to start the quadtree division. The last step is to optimize the parameters ( $C$ , $σ$ ) based on the first and second steps. The details of this process can be found in Beltrami and da Silva.¹⁷

Discussion 2: Compared with the latest convolutional neural network-based deep learning methods, the main advantages of the proposed method are that the dynamics features are first used to construct the training SVM classifier, and the extracted features have a clear physical meaning. The constructed classification model is more sensitive to the dynamics changes under the unbalance faults condition. Otherwise, the required sample size is not too large, which is suitable for real applications. Under the condition of arranging a large-scale neural network, the main limitation of the proposed method is that it requires a certain amount of computing power to deal with the RBF structure.

Conclusion

In this study, dynamics features including SE, FE, and PE of rotor unbalance faults were firstly extracted based on the dynamical trajectory (in which was obtained by using the deterministic learning theory), and the diagnosis results based on the dynamics features and Gaussian-SVM achieved more accurate classification of different degrees of rotor unbalance faults. Based on the experiment results, it is indicated that the extracted rotor dynamics features were more sensitivity to the changes of rotor unbalance faults. For the further investigation, the prognosis of rotor unbalance faults and the diagnosis of other rotor faults (e.g. rub-impact faults) based on the dynamics features can be considered.

Research Data

sj-rar-1-mac-10.1177_00202940221135917 – Supplemental material for Fault detection and classification of the rotor unbalance based on dynamics features and support vector machine

Supplemental material, sj-rar-1-mac-10.1177_00202940221135917 for Fault detection and classification of the rotor unbalance based on dynamics features and support vector machine by Lan Lan, Xiao Liu and Qian Wang in Measurement and Control

This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/open-access-at-sage).

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was partly supported by the Henan Province Foundation for University Key Teacher [grant numbers 2018GGHS243], partly supported by the Doctoral Research Project of Zhengzhou University of Light Industry [grant number 2020BSJJ078], and partly supported by the Science and Technology Research Project of Henan Province [grant numbers 222102220007].

ORCID iD

Qian Wang

References

Nath

Udmale

Singh

SK.

Role of artificial intelligence in rotor fault diagnosis: a comprehensive review. Artif Intell Rev 2021; 54: 2609–2668.

Ambur

Rinderknecht

Unbalance detection in rotor systems with active bearings using self-sensing piezoelectric actuators. Mech Syst Signal Process 2018; 102: 72–86.

Gao

QL.

Unbalance fault diagnosis of helicopter rotor based on wavelet transform and neutral network. J Nanjing Univ Aeronaut Astronaut 2017; 49(2): 212–218.

Zhou

Tang

YL.

Unbalanced fault feature extraction for wind power gearbox based on improved VMD. J Vib Shock 2020; 39(5): 170–176.

Chen

Wan

Sheng

, et al. Characteristic analysis of blade mass unbalance fault of doubly-fed induction generator based on Hilbert transform. Trans China Electrotech Soc 2021; 36(24): 5225–5236.

Wang

K-nearest neighbors based methods for identification of different gear crack levels under different motor speeds and loads: revisited. Mech Syst Signal Process 2016; 70-71: 201–208.

Wang

Zhang

, et al. Early rub-impact fault detection of rotor systems via deterministic learning. Control Eng Pract 2022; 124: 105190.

Soleymani

Nekoui

Moarefianpour

A novel robust fault detection method for induction motor rotor by using unknown input observer. Syst Sci Control Eng 2019; 7(1): 109–115.

Zheng

Liu

Zhao

, et al. A novel iterative learning control method and control system design for active magnetic bearing rotor imbalance of primary helium circulator in high-temperature gas-cooled reactor. Meas Control 2020; 53(3-4): 474–484.

10.

Liu

Zhao

Liang

, et al. A rotor fault diagnosis method based on BP-adaboost weighted by non-fuzzy solution coefficients. Measurement 2022; 196: 111280.

11.

Hübner

Pinheiro

de Souza

, et al. Detection of mass imbalance in the rotor of wind turbines using support vector machine. Renew Energy 2021; 170: 49–59.

12.

Wang

Zhang

Xing

, et al. Sparse representation theory for support vector machine kernel function selection and its application in high-speed bearing fault diagnosis. ISA Trans 2021; 118: 207–218.

13.

Gohari

Eydi

AM.

Modelling of shaft unbalance: modelling a multi discs rotor using K-nearest neighbor and decision tree algorithms. Measurement 2020; 151: 107253.

14.

Yan

Guo

Rotor unbalance fault diagnosis using DBN based on multi-source heterogeneous information fusion. Procedia Manuf 2019; 35: 1184–1189.

15.

Fei

Bai

Tang

, et al. Quantitative diagnosis of rotor vibration fault using process power spectrum entropy and support vector machine method. Shock Vib 2014; 2014: 957531.

16.

Huang

Fan

Wei

, et al. An intelligent fault identification method of rolling bearings based on SVM optimized by improved GWO. Syst Sci Control Eng 2019; 7(1): 289–303.

17.

Beltrami

da Silva

ACL.

A grid-quadtree model selection method for support vector machines. Expert Syst Appl 2020; 146: 113172.

18.

Wang

, et al. Entropy based fault classification using the case western reserve university data: a benchmark study. IEEE Trans Reliab 2020; 69(2): 754–767.

19.

Wang

Dynamic feature extraction of nonlinear systems with deterministic learning theory and spatio-temporal Lempel-Ziv complexity. Zidonghua Xuebao 2018; 44(10): 1812–1823.

20.

Wang

Hill

DJ.

Deterministic learning theory for identification, recognition and control. Boca Raton, FL: CRC Press, 2009.

21.

Wang

Chen

Rapid detection of small oscillation faults via deterministic learning. IEEE Trans Neural Netw 2011; 22(8): 1284–1296.

22.

Sun

Wang

, et al. The cardiodynamicsgram based early detection of myocardial ischemia using the Lempel-Ziv complexity. IEEE Access 2020; 8: 207894–207904.

23.

Zeng

Yuan

, et al. A novel technique for the detection of myocardial dysfunction using ECG signals based on hybrid signal processing and neural networks. Soft Comput 2021; 25: 4571–4595.

24.

Yang

, et al. Vibration behaviour of a geometrically nonlinear rotor under misalignment-unbalance coupled fault. J Vib Shock 2021; 40(7): 45–52.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

108.45 MB