Next Article in Journal
An Evaluation Model for Node Influence Based on Heuristic Spatiotemporal Features
Previous Article in Journal
Furstenberg Family and Chaos for Time-Varying Discrete Dynamical Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intelligent Fault Diagnosis Method for Rotating Machinery Based on Recurrence Binary Plot and DSD-CNN

School of Mechanical Engineering, Xinjiang University, Urumqi 830046, China
*
Author to whom correspondence should be addressed.
Entropy 2024, 26(8), 675; https://doi.org/10.3390/e26080675
Submission received: 26 July 2024 / Accepted: 7 August 2024 / Published: 9 August 2024
(This article belongs to the Section Signal and Data Analysis)

Abstract

:
To tackle the issue of the traditional intelligent diagnostic algorithm’s insufficient utilization of correlation characteristics within the time series of fault signals and to meet the challenges of accuracy and computational complexity in rotating machinery fault diagnosis, a novel approach based on a recurrence binary plot (RBP) and a lightweight, deep, separable, dilated convolutional neural network (DSD-CNN) is proposed. Firstly, a recursive encoding method is used to convert the fault vibration signals of rotating machinery into two-dimensional texture images, extracting feature information from the internal structure of the fault signals as the input for the model. Subsequently, leveraging the excellent feature extraction capabilities of a lightweight convolutional neural network embedded with attention modules, the fault diagnosis of rotating machinery is carried out. The experimental results using different datasets demonstrate that the proposed model achieves excellent diagnostic accuracy and computational efficiency. Additionally, compared with other representative fault diagnosis methods, this model shows better anti-noise performance under different noise test data, and it provides a reliable and efficient reference solution for rotating machinery fault-classification tasks.

1. Introduction

Rolling bearings and gears constitute the crucial components of rotating machinery, and their health conditions significantly impact the performance, stability, and lifespan of mechanical equipment [1,2]. Vibration signals exhibit distinct variations across different fault types in rotating machinery. By implementing the real-time monitoring of vibration signals, it is possible to effectively prevent major accidents [3]. Hence, the development of a practical algorithmic tool to accurately and reliably diagnose the intricate, non-linear correlation between a raw vibration signal and the fault modes of rolling bearings and gears holds immense practical significance.
Traditional intelligent diagnostic algorithms based on expert systems are complex, time-consuming, and lack flexibility, rendering them inadequate for meeting the demands of industrial production in the era of internet big data. Exploring intelligent mechanical fault diagnosis methods can diminish the level of manual intervention and contribute to the safe and dependable operation of mechanical equipment [4]. In recent years, a frontier research direction called deep learning has proposed a novel approach for feature extraction based on neural networks. Compared with the classical shallow machine learning method, deep learning is a hierarchical representation learning method based on data. It learns deep abstract features by establishing neurons with multiple hidden layers, and it comprehensively considers the steps of feature extraction and classification so as to effectively improve the classification or prediction accuracy [5].
For the fault diagnosis tasks of mechanical equipment, a CNN aptly captures the intricate mapping relationship between various signals and corresponding fault states by learning from sensor-collected data [6]. Since the signal collected by the sensor is a one-dimensional time-domain signal, a large number of scholars have directly inputted the signal into a one-dimensional CNN to accomplish a fault diagnosis task. Wu et al. [7] proposed a method based on one-dimensional convolutional neural networks to autonomously learn valuable features from the original vibration signal of a rotor system, facilitating prompt and accurate fault detection. Mo et al. [8] integrated a learnable variational kernel into a one-dimensional CNN framework, with a focus on extracting crucial data features pertaining to faults, thereby achieving a commendable performance on small sample datasets. Junior et al. [9] introduced a multi-head, one-dimensional convolutional neural network that employed accelerometers measuring in two distinct directions to detect six different types of motor faults.
With the vigorous development of modern signal-processing technology, the method of preprocessing vibration signals at the front end has become widespread, yielding promising results. Deng et al. [10] employed an envelope demodulation analysis method for the feature extraction and fault diagnosis of a rolling bearing. Jin et al. [11] introduced a locomotive-bearing fault-diagnosis approach based on variational mode decomposition (VMD) and an improved convolutional neural network. Wang et al. [12] performed swarm decomposition (SWD) to decompose each segment of an original signal into several oscillation components (OC), and then they combined this with the improved MRDE algorithm to accomplish the intelligent identification and classification of various fault signals. Shao et al. [13] proposed an enhanced, modified, stacked autoencoder (MSAE) that employed adaptive Morlet wavelets for the automatic diagnosis of different fault types and severities in rotating machinery. Furthermore, Miao et al. [14] introduced a novel decomposition theory known as feature mode decomposition (FMD) for mechanical fault feature extraction.
Converting one-dimensional signals into feature images creates a link between low-level visual features and high-level visual features, which can emphasize recognizable characteristics. This technique is an efficient means of achieving information recognition and recovering high-dimensional features from low-dimensional information [15]. Wen et al. [16] employed a sliding window to normalize one-dimensional vibration signals into two-dimensional grayscale images which were then combined with a CNN, obviating the process of manual feature extraction. Xu et al. [17] input a time-frequency diagram, which was obtained through the continuous wavelet transform (CWT) of the vibration signal, into a CNN as a fault feature map for fault diagnosis. Shao et al. [18] retrained a CNN pretrained by ImageNet using a time-frequency diagram, effectively shortening the training time while extracting discernible fault features. Zhu et al. [19] transformed the vibration signals from multiple sensors into symmetrized dot pattern images (SDP) and employed a CNN to differentiate the rotor faults. Choudhary et al. [20] utilized a CNN to identify the thermal images of rolling bearings under five fault conditions, suggesting that infrared thermal imaging enables the non-contact early detection of faults without being influenced by speed. Zhi et al. [21] proposed a new entropy-aided meshing-order modulation (EMOM) indicator to capture the most sensitive modulation frequency area embedded in a signal, and they developed a wind-turbine fault-diagnosis method. Tang et al. [22] presented a composite model that combined improved mode decomposition, Gram angular summation field (GASF), and convolutional neural networks for the automated identification of the health statuses of complex mechanical systems. Xiong et al. [23] employed a mutual dimensionless theory and a similar Gram matrix to process bearing-fault vibration signals and subsequently integrated them with a convolutional neural network, which significantly reduced the training time. Bai et al. [24] proposed a frequency-domain Gram angular field (FDGAF) algorithm which could intelligently classify these feature images under the condition of 30 samples through a transfer learning network. Additionally, Bai [25] proposed a spectral Markov transition field (SMTF) algorithm which constructed the first-order Markov transition matrix of the frequency domain signal and represented the spectral characteristics of the vibration signal in image form.
In summary, although the above methods have achieved excellent results, the correlation characteristics between the time-series of fault signals have not been fully explored and utilized. Concurrently, there is still much room for improvement in the accuracy and computational efficiency of such a model. This paper proposes an intelligent fault diagnosis method for rotating machinery called RBP-DSD-CNN, which is based on a recurrence binary plot and a lightweight, deep, separable, dilated convolutional neural network. The main contributions of this paper are summarized as follows:
(a)
A recursive quantization technique is introduced into the field of fault diagnosis. By leveraging the correlation characteristics of time-series data, feature information is extracted from the internal structures of fault signals, and this effectively enhances the expression capability of the features and mitigates the loss of weak information in the fault signals.
(b)
Considering the characteristics of mechanical monitoring signals, a minimum mutual information method and the Cao method are employed to determine the optimal phase space-reconstruction parameters for each category.
(c)
A DSD-CNN is developed for feature extraction and fault classification. This model adopts a lightweight structure, improving computational efficiency without compromising diagnostic accuracy.

2. Methodology

2.1. Recurrence Binary Plot

A recurrence plot (RP) is a recursive visualization method proposed by Eckmann [26] et al. on the basis of Poincare’s theorem, which is applied for the qualitative analysis of nonlinear dynamical systems. An RP portrays the recursive phenomenon of phase points and the evolution with time of the one-dimensional signals after the phase-space reconstruction in a two-dimensional image, which can intuitively characterize the dynamic characteristics implied by a system.
Time-series data can be classified using unique repetitive behaviors (e.g., periodicity and irregular periodicity). The essence of an RP is that it is a signal-processing method which is used to encode the cyclic behavior in a time-series into a two-dimensional image for representation. By applying the Takens embedding theorem, the phase-space reconstruction of the original time-series   { x i   , i = 1 , 2 , , N } with a length of N is performed, and the phase-space vector obtained is as follows:
Y   i = ( x i , x i + τ , , x i + ( m 1 ) τ , ) i = 1 , 2 , , N ( m 1 ) τ ,
where τ denotes the delay time and m represents the embedding dimension. The determination of τ and m will be described in detail in Section 2.2. On the basis of the reconstructed phase-space vector, the recursive matrix can be defined as follows:
R i j = Θ ( r X i Y j ) ,
where Θ ( ) is the Heavside function (as shown in Equation (3)), represents the Euclidean norm, X i Y j is the distance between any two points in the reconstructed phase-space vector, and r is the distance threshold for determining whether there is a recursive phenomenon between two points in the phase-space. Obviously, the recursive matrix is a binary matrix consisting of 0 and 1.
Θ ( x ) = { 1 , x 0 0 , x < 0 .
When traversing any two points in the phase-space and after calculating all instances of R i j , a recurrence binary plot (RBP) with black and white points symmetrically distributed around the main diagonal can be obtained.
The flowchart illustrating the conversion of the raw signal to an RBP is shown in Figure 1. Figure 2 illustrates an RBP obtained through the recursive coding of diverse signals. With the increase in the non-stationarity, the detail texture becomes more significant.

2.2. Determination of the Phase-Space Reconstruction Parameters

The phase-space reconstruction serves as the foundation for a recursive quantitative analysis and recurrence plot construction. Therefore, determining the delay time τ and embedding dimension m is crucial. Takens posits that τ and m are theoretically independent. Thus, τ can be initially determined, upon which m is computed.
The minimum mutual information method was employed in this study to determine the delay time for two discrete information systems ( S = { s i } i = 1 m and Q = { q j } j = 1 n ), and their mutual information is given by the following:
I ( S , Q ) = H ( S ) + H ( Q ) H ( S , Q ) ,
where H ( S ) and H ( Q ) represent the information entropies of the two systems defined by Shannon’s information theory and H ( S , Q ) denotes the joint probability distribution. When computing an actual time-series, S corresponds to the original time series x ( n ) and Q corresponds to the time-series x ( n + τ ) after being delayed by τ . When the mutual information value I ( S , Q ) reaches its first local minimum, it can be inferred that the corresponding τ value at this point is the optimal delay time, indicating that the time-series x ( n ) and x ( n + τ ) are maximally uncorrelated.
Subsequently, the Cao method [27] was utilized to estimate the embedding dimension, and this was an improvement to the False Nearest Neighbor (FNN) method, which was designed to reduce the influence of the manual threshold parameter selection. The definition of the Cao method is as follows:
E ( m ) = 1 N m τ i = 1 N m τ a ( i , m ) ; E * ( m ) = 1 N m τ i = 1 N m τ | x ( i + m τ ) x N N ( i + m τ ) | .
By calculating E 1 ( m ) = E ( m + 1 ) / E ( m ) under different values of m and introducing E 2 ( m ) = E * ( m + 1 ) / E * ( m ) as an auxiliary criterion, we could determine the optimal embedding dimension. If E 1 ( m ) and E 2 ( m ) are stabilized when m exceeds a certain value, the corresponding m value at this point can be considered the optimal embedding dimension.
In Figure 3a, it is evident that the mutual information reaches its first minimum value when the delay time is 10, indicating the optimal delay time to be 10. Additionally, in Figure 3b, the stabilization of the changes in E1 and E2 occurs when the embedding dimension equals 3, suggesting the optimal embedding dimension to be 3.

2.3. Deep, Separable, Dilated Convolutional Neural Network (DSD-CNN)

Recently, deep learning methods have been widely used to unveil the hidden mapping between monitoring data and mechanical health status. However, these models are typically constructed for large-scale image classification tasks, and as the network depth increases, the computational and parameter complexity increase significantly. Considering the characteristics of mechanical monitoring signals, deep and separable convolutions (DSC) and dilated convolutions (DC) have been introduced into construct models.
A DSC [28] decomposes convolution operations into the following two steps: depth-wise convolution and point-wise convolution. Firstly, independent convolution kernels are applied to each channel of an input feature map to extract channel-specific features. Then, 1 × 1 convolution kernels are used for point-wise convolution, linearly combining and merging features from different channels to obtain the final output feature map, as follows:
C o n v ( W D , W P , x ) ( i , j ) = k K W k x ( i , j ) [ W P , m , n M , N W ( m , n ) x ( i + m , j + n ) ] ,
where W represents the parameters of the kernel weight matrix. The input is denoted as x , while ( i , j ) signifies the coordination of output features. m , n , k , respectively, denote the width, height, and number of channels of the convolution kernel. The subscripts D and P represent depth-wise convolution and point-wise convolution operations, respectively.
A DC [29] expands a receptive field effectively and enhances a model’s performance by introducing gaps and holes within a convolutional kernel to augment the distance between the kernel and the input. The computational process of a DC is as follows;
F i = F i 1 + ( k i 1 ) r i n = 1 i 1 s n ,
where F i represents the receptive field of the convolution kernel in the i-th convolutional layer, k i denotes the size of the convolution kernel, r i is the dilation factor, and s n signifies the stride.
A coordinate attention (CA) [30] is often regarded as a computing unit used to enhance the feature representation ability, and it can take any intermediate tensor as an input and obtain an output of the same size with enhanced representation ability. The structure of a CA is illustrated in Figure 4. In the first step, the coordinate information is embedded into the channel attention. The second step is to generate a coordinate attention.
In this study, a lightweight convolutional neural network embedded with a CA was developed. The proposed model was designed to meet the requirements for both model performance and computational efficiency. The attention module enhanced the sensitivity of the feature location information, enabled the network to concentrate on learning essential features, and improved the training aggregation rate. The structure of the model is illustrated in Figure 5.

3. Overall Framework of the Proposed Approach

Based on the above, a fault diagnosis framework was constructed as shown in Figure 6. The main steps were as follows.
Step 1: Data acquisition. The raw vibration signals of the different health states of the rolling bearings and gearboxes were collected from the test bench.
Step 2: Signal converted to an RBP. The raw vibration signal was divided into a series of one-dimensional samples through overlapping sampling, and then it was recursively encoded into an RBP as the input of the model.
Step 3: Build model learning system. The RBP, obtained in the previous step, was randomly divided into a training set and a test set at a ratio of 8:2, and afterwards, the model was trained with the training set.
Step 4: Fault diagnosis. The test set was input into the trained model to obtain the diagnosis results.

4. Experimental Study

In this study, experimental verification was conducted utilizing the bearing dataset of the Case Western Reserve University (CWRU) [31] and the gear dataset of the Wind Turbine Drivetrain Diagnostics Simulator (WTDS), which was collected in the authors’ laboratory.

4.1. Dataset Information

(1) CWRU: The data utilized in this study were sampled at a frequency of 12 KHz under the rotational speed of 1797 rpm from a drive-end bearing, and they contained normal (N), inner race fault (IF), outer race fault (OF), and ball fault (BF) types. The fault types could be further subdivided based on the fault sizes (14 mils and 21 mils). Consequently, seven different healthy states are considered. The specific information can be found in Table 1.
(2) WTDS: The data utilized in this study were collected via a three-axis acceleration sensor at a frequency of 20,480 Hz, which covered four fault types and contained a total of 2800 samples. The fault types could be further subdivided into seven health conditions, and the rotational speed was 2000 rpm. The detailed information is presented in Table 2.
The mechanical vibration signals collected above were reconstructed in phase space, as described of Section 2.1, and encoded into a recurrence binary plot. As depicted in Figure 7 and Figure 8, the corresponding health status is shown above each image.

4.2. Experimental Parameters

The model was built based on Pytorch-1.12.0 and trained by GPU. The operating system was Windows 11, with an Intel Core i5-13400F CPU running at 2.50 GHz. The GPU employed was an NVIDIA GeForce RTX 3060. For the selection of the initial learning rate, according to Leslie’s theory [32], we first set a smaller learning rate for the optimizer, and then we increased the learning rate in the data training of each batch until the loss began to become larger. The loss calculated by each batch was counted, and so the curves of the loss and learning rate are shown in Figure 9a. The optimal learning rate range could be obtained, as shown in the figure, and so the initial learning rate in this paper was set to 0.001 and the learning rate reduction strategy adopted “Reduce LR On Plateau”. In order to evaluate the influence of the different batch sizes on the classification performance of the model, overall accuracy (OA), average accuracy (AA), and a Kappa coefficient were introduced as quantitative evaluation indexes. The model was tested with batch sizes of 8, 16, 20, 24, and 32, respectively, and the results are shown in Figure 9b. When the batch size was 20, the OA, AA, and Kappa coefficient values were the highest. Therefore, the batch size was set to 20. While the model’s weight-update back-propagation algorithm was small batch gradient descent, the loss function was the cross-entropy loss function, and the optimizer was Adam.

4.3. Model Rationality Analysis

The ablation experiment was carried out by gradually removing some components in the DSD-CNN to observe their effects on the performance of the model so as to evaluate the importance of these components and explain the rationality of the model’s settings. The dataset shown in Table 1 was selected as the model’s input, and the experimental results are shown in Figure 10 and Table 3, respectively.
It can be seen from Table 3 that the SeparableConv block had the greatest influence on the performance and computational efficiency of the model. In Figure 10a, after removing this component, the model could only clearly classify the normal and 21 mils outer race faults, and the remaining health state boundaries were blurred and difficult to distinguish. Compared with the SeparableConv block, the DepthwiseConv block and the CA module had less impact on the performance of the model, and especially after removing the CA module, the model could still achieve a classification accuracy of over 96%. However, considering the influence of computational efficiency, the existence of the above two components was necessary. In general, the proposed DSD-CNN was reasonable and could balance computational efficiency and classification performance.

4.4. Comparison Method

To verify the effectiveness of the proposed method, the diagnostic performance and computational efficiency of the proposed model were compared with other algorithms, including classical deep convolutional neural networks and representative lightweight models.
(1)
ResNet [33]: In this study, the widely used ResNet50 was selected as a comparison method, and it consisted of 4 conv blocks, 12 identity blocks, a total of 49 convolutional layers, and 1 Softmax classifier.
(2)
ShuffleNetV2 [34]: ShufflenetV2 with a 2.0× output channel was constructed for comparison, and it consisted of 17 convolution layers and 2 pooling layers with the original convolution kernel size of 112 × 112, and it was halved to 7 × 7 layer by layer.
(3)
EfficientNet [35]: Considering the model’s depth and computational efficiency, EfficientNetB3 was selected as a comparison method, and it was mainly composed of two convolutional layers and seven stacked MBConv layers.
(4)
DenseNet [36]: The DenseNet 161 used in this study mainly consisted of four dense blocks and three transition layers.
(5)
Light-weight CNN (LW-CNN) [37]: The LW-CNN comprised a convolution layer consisting of eight filters of 3 × 3 size, followed by a maxpool layer and a fully connected layer with a softmax activation.
Moreover, all of the above methods shared the same experimental parameters as the proposed methods.

4.5. Case Study I

In this section, the gear-fault data collected from the WTDS test rig were employed to validate the model. The maximum number of iterations was set to 100.

4.5.1. Model Performance Analysis

The model loss and accuracy curves for the WTDS dataset are shown in Figure 11. It could be concluded that the DSD-CNN proposed in this paper demonstrated the characteristics of rapid model aggregation and high diagnostic accuracy. For the training set, after nine iterations, the training accuracy reached 99.56%. Subsequently, the training accuracy reached 100% after 32 iterations, indicating that the model had converged on the training set. For the test set, after approximately 26 iterations, the model had achieved a 99% test accuracy. Furthermore, after 38 iterations, the test accuracy reached 99.82%, indicating that the proposed model had been trained to the optimal level and had good classification ability for gear-fault signals.
For the purpose of evaluating the performance of the proposed model, indicators such as precision, recall, and F1-score were introduced in the study. The classification results of the model on the WTDS dataset are presented in Table 4.
The corresponding confusion matrix is depicted in Figure 12a, and the results demonstrated that the model could accurately classify fault states. T-distributed Stochastic Neighbor Embedding (T-SNE) [38] is a nonlinear, unsupervised dimensionality reduction method which could reduce the vector dimension and capture the complex manifold structure of the original data. By means of the T-SNE, the output features at the dense layer of the model were visualized in two-dimensional space, as shown in Figure 12b. Despite the fact that there were some errors in reducing the high-dimensional feature information to two-dimensional space, the purpose of the two-dimensional visualization was to provide an intuitive means for judging whether the features extracted by the model were discriminative. As presented in Figure 10b, the model misclassified one sample in the 1.2 Nm gear-tooth crack as 1.2 Nm tooth-surface wear (as shown in the red ellipse annotation in Figure 12b). A possible reason is that the two fault categories were distributed on the gear surface. When the gear was engaged, there was a gap in the meshing line, and the collected vibration signals were similar; however, these were different from the vibration impact generated when the broken tooth was engaged, resulting in a misclassification by the model.
An ROC curve refers to a receiver operating characteristic curve, which comprehensively reflects the sensitivity and specificity of continuous variables. A curve was drawn with 1-specificity as the abscissa and sensitivity as the ordinate. The larger the area under the curve, namely, the larger the AUC value, the higher the accuracy of the experimental results. The results demonstrated that the AUC value of the model for each type of fault on the WTDS test rig was 1, as shown in Figure 13.
Through a comprehensive analysis, it could be concluded that the network proposed in this paper possessed a strong feature extraction ability and could focus on the essential information contained in an input feature map. Simultaneously, the inter-class separation and intra-class compactness were achieved, which indicated the superiority and feasibility of the proposed model for the fault diagnosis of rotating machinery.

4.5.2. Comparison between Different Methods

Each model experiment was repeated 10 times to reduce the effect of randomness. The comparison experiment results are depicted in Figure 14, and the average test accuracy of the experiment is represented in Table 5.
In a comparative experiment conducted on the WTDS test rig, all six models exhibited good diagnostic accuracy and achieved test accuracies of over 95%. Compared with the deep models ResNet and DenseNet, the proposed method achieved better diagnostic accuracy than the above models while saving on computational efficiency. In addition, compared with the lightweight models ShuffleNet and LW-CNN, the proposed method achieved better results at the expense of some computational efficiency. In general, the proposed method achieved a higher classification accuracy and faster training speed through using a lightweight convolutional neural network embedded with an attention module compared with the above classic deep algorithm and representative lightweight network.

4.5.3. Anti-Noise Analysis

In practical industrial production, mechanical equipment usually runs under variable working conditions, and the data collected from the equipment contain certain background noise interference. In order to further study the performance of the algorithm in a noisy environment, different levels of Gaussian white noise were added to the original test dataset, and the signal-to-noise ratios (SNRs) were 8 dB and 4 dB, respectively. The noisy test data were input into the above well-trained model for classification. Table 6 shows the classification results of the model under different noise levels.
It can be seen from the results that the diagnostic accuracy of all methods was greatly reduced compared with that before the noise addition. This was because the discriminant information in the original vibration signal would be disturbed and covered by noise, resulting in a large difference in the distribution of the training data and test data, thus further increasing the difficulty of the model in learning the discriminant features. When the signal-to-noise ratio was high (SNR = 8 dB), all methods obtained an average accuracy higher than 87%, showing a certain anti-noise ability and generalization performance. The proposed method achieved an average diagnostic accuracy of more than 95%, indicating that the proposed model could maintain good feature learning ability under small noise interference.

4.6. Case Study II

In recent years, the CWRU dataset has been widely adopted by numerous researchers in fault diagnosis experiments to validate the performance of the studied models. The experimental conditions remained the same as those of the WTDS experiment described above.

4.6.1. Model Performance Analysis

The loss and accuracy curves of the model based on the CWRU dataset are shown in Figure 15. For the training set, after approximately 18 iterations, the model had achieved a 99.5% training accuracy. Furthermore, after 33 iterations, the training accuracy reached 100%. For the test set, after 26 iterations, the test accuracy reached 99%. This further improved to 100% after approximately 33 iterations, marking the optimal level of model training. Precision, recall and F1-score were still employed to evaluate the performance of the model. The results are represented in Table 7.
The classification results of the model are depicted in Figure 16. The results exhibited that the model could accurately classify each fault state, with the value of the three indicators for each fault type being 100%.

4.6.2. Comparison between the Different Methods

Ten comparative experiments were carried out. The comparison test results are shown in Figure 17, and the average test accuracy is presented in Table 8.
In the above comparative experiments, the diagnostic accuracy was higher than 91%. Compared with the other approaches, the model proposed in this paper demonstrated a superior performance, achieving a diagnostic accuracy of 100%. Concurrently, the computational efficiency was improved under the premise of ensuring the performance of the model.

4.6.3. Anti-Noise Analysis

After the model was well-trained on the training dataset, two different Gaussian white noises were also added to the test data, and the anti-noise performance of the different methods was studied. The diagnostic results are shown in Table 9.
The classification accuracy of all algorithms decreased with the increase in noise level. The anti-noise performance of the model trained by the CWRU dataset was significantly weaker than that of the WTDS dataset. The reason may have been that the gear to be diagnosed in the WTDS test rig was connected to the motor through another pair of gears, and additional modulation components were introduced in the transmission path. Therefore, the data distribution was quite different, which was equivalent to improving the anti-noise performance of the model.

5. Conclusions

In this paper, an intelligent diagnosis method based on a recurrence binary plot and a DSD-CNN is proposed for the condition-monitoring of mechanical equipment. In order to better extract the data structure of the mechanical monitoring signal, a recursive quantization technique was used to encode the vibration signal, and optimal phase-space reconstruction parameters were determined. A lightweight convolutional neural network was developed for feature extraction and fault classification. Through the comparative analysis of two experimental cases, the superiority and robustness of the proposed method were verified. Overall, the experimental results show that the proposed method significantly improves the accuracy and computational efficiency of mechanical fault classification. In addition, the proposed method shows better anti-noise performance under different noise test data compared with other approaches.
In practical engineering scenarios, acquiring sufficient labeled data for model training often proves infeasible. Therefore, the research work of domain adaptation algorithms based on transfer-learning for variable working conditions and cross-machines has been started.

Author Contributions

Funding acquisition, H.W. and W.S.; methodology, Y.S.; software, Y.S.; supervision, W.S.; visualization, R.B.; writing—original draft, Y.S.; writing—review and editing, H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Natural Science Foundation of the Xinjiang Uygur Autonomous Region under grant 2022D01C390 and in part by the Key Research and Development Program of the Xinjiang Uygur Autonomous Region under grant 2022B02016-1.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The public dataset from CWRU can be found at https://engineering.case.edu/bearingdatacenter (accessed on 1 March 2022). The signals collected from the WTDS test rig in the authors’ laboratory are owned by the authors’ institution, and they are not allowed to be used in a public way.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Liang, G.; Gumabay, M.V.N.; Zhang, Q.; Zhu, G. Smart Fault Diagnosis of Rotating Machinery with Analytics Using Deep Learning Algorithm. In Proceedings of the 2022 5th International Conference on Advanced Electronic Materials, Computers and Software Engineering (AEMCSE), Wuhan, China, 22–24 April 2022; pp. 660–665. [Google Scholar]
  2. Kong, X.; Li, X.; Zhou, Q.; Hu, Z.; Shi, C. Attention Recurrent Autoencoder Hybrid Model for Early Fault Diagnosis of Rotating Machinery. IEEE Trans. Instrum. Meas. 2021, 70, 2505110. [Google Scholar] [CrossRef]
  3. Ye, M.; Yan, X.; Jia, M. Rolling Bearing Fault Diagnosis Based on VMD-MPE and PSO-SVM. Entropy 2021, 23, 762. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, H.; Sun, W.; He, L.; Zhou, J. Rolling Bearing Fault Diagnosis Using Multi-Sensor Data Fusion Based on 1D-CNN Model. Entropy 2022, 24, 573. [Google Scholar] [CrossRef] [PubMed]
  5. Li, Y.; Song, Y.; Jia, L.; Gao, S.; Li, Q.; Qiu, M. Intelligent Fault Diagnosis by Fusing Domain Adversarial Training and Maximum Mean Discrepancy via Ensemble Learning. IEEE Trans. Ind. Inform. 2020, 17, 2833–2841. [Google Scholar] [CrossRef]
  6. Xia, M.; Li, T.; Xu, L.; Liu, L.; de Silva, C.W. Fault Diagnosis for Rotating Machinery Using Multiple Sensors and Convolutional Neural Networks. IEEE/ASME Trans. Mechatron. 2017, 23, 101–110. [Google Scholar] [CrossRef]
  7. Wu, X.; Penga, Z.; Renc, J.; Chenga, C.; Zhanga, W.; Wang, D. Rub-Impact Fault Diagnosis of Rotating Machinery Based on 1-D Convolutional Neural Networks. IEEE Sensors J. 2019, 20, 8349–8363. [Google Scholar] [CrossRef]
  8. Mo, Z.; Zhang, Z.; Tsui, K.-L. The Variational Kernel-Based 1-D Convolutional Neural Network for Machinery Fault Diagnosis. IEEE Trans. Instrum. Meas. 2021, 70, 3523210. [Google Scholar] [CrossRef]
  9. Junior, R.F.R.; Areias, I.A.d.S.; Campos, M.M.; Teixeira, C.E.; da Silva, L.E.B.; Gomes, G.F. Fault detection and diagnosis in electric motors using 1d convolutional neural networks with multi-channel vibration signals. Measurement 2022, 190, 110759. [Google Scholar] [CrossRef]
  10. Deng, W.; Li, Z.; Li, X.; Chen, H.; Zhao, H. Compound Fault Diagnosis Using Optimized MCKD and Sparse Representation for Rolling Bearings. IEEE Trans. Instrum. Meas. 2022, 71, 3508509. [Google Scholar] [CrossRef]
  11. Jin, Z.; Chen, D.; He, D.; Sun, Y.; Yin, X. Bearing Fault Diagnosis Based on VMD and Improved CNN. J. Fail. Anal. Prev. 2022, 23, 165–175. [Google Scholar] [CrossRef]
  12. Wang, H.; Sun, W.; He, L.; Zhou, J. Intelligent Fault Diagnosis Method for Gear Transmission Systems Based on Improved Multi-Scale Reverse Dispersion Entropy and Swarm Decomposition. IEEE Trans. Instrum. Meas. 2022, 71, 3502813. [Google Scholar] [CrossRef]
  13. Shao, H.; Xia, M.; Wan, J.; de Silva, C.W. Modified Stacked Autoencoder Using Adaptive Morlet Wavelet for Intelligent Fault Diagnosis of Rotating Machinery. IEEE/ASME Trans. Mechatron. 2021, 27, 24–33. [Google Scholar] [CrossRef]
  14. Miao, Y.; Zhang, B.; Li, C.; Lin, J.; Zhang, D. Feature Mode Decomposition: New Decomposition Theory for Rotating Machinery Fault Diagnosis. IEEE Trans. Ind. Electron. 2022, 70, 1949–1960. [Google Scholar] [CrossRef]
  15. Liu, Z.; Pei, M.; He, Q.; Wu, Q.; Jackson, L.; Mao, L. A novel method for polymer electrolyte membrane fuel cell fault diagnosis using 2D data. J. Power Sources 2020, 482, 228894. [Google Scholar] [CrossRef]
  16. Wen, L.; Li, X.; Gao, L.; Zhang, Y. A New Convolutional Neural Network-Based Data-Driven Fault Diagnosis Method. IEEE Trans. Ind. Electron. 2017, 65, 5990–5998. [Google Scholar] [CrossRef]
  17. Xu, G.; Liu, M.; Jiang, Z.; Söffker, D.; Shen, W. Bearing Fault Diagnosis Method Based on Deep Convolutional Neural Network and Random Forest Ensemble Learning. Sensors 2019, 19, 1088. [Google Scholar] [CrossRef]
  18. Shao, S.; McAleer, S.; Yan, R.; Baldi, P. Highly Accurate Machine Fault Diagnosis Using Deep Transfer Learning. IEEE Trans. Ind. Inform. 2018, 15, 2446–2455. [Google Scholar] [CrossRef]
  19. Zhu, X.; Hou, D.; Zhou, P.; Han, Z.; Yuan, Y.; Zhou, W.; Yin, Q. Rotor fault diagnosis using a convolutional neural network with symmetrized dot pattern images. Measurement 2019, 138, 526–535. [Google Scholar] [CrossRef]
  20. Choudhary, A.; Mian, T.; Fatima, S. Convolutional neural network based bearing fault diagnosis of rotating machine using thermal images. Measurement 2021, 176, 109196. [Google Scholar] [CrossRef]
  21. Zhi, S.; Wu, H.; Shen, H.; Wang, T.; Fu, H. Entropy-Aided Meshing-Order Modulation Analysis for Wind Turbine Planetary Gear Weak Fault Detection under Variable Rotational Speed. Entropy 2024, 26, 409. [Google Scholar] [CrossRef]
  22. Tang, H.; Liao, Z.; Chen, P.; Zuo, D.; Yi, S. A Novel Convolutional Neural Network for Low-Speed Structural Fault Diagnosis Under Different Operating Condition and Its Understanding via Visualization. IEEE Trans. Instrum. Meas. 2020, 70, 3501611. [Google Scholar] [CrossRef]
  23. Xiong, J.; Li, C.; Wang, C.-D.; Cen, J.; Wang, Q.; Wang, S. Application of Convolutional Neural Network and Data Preprocessing by Mutual Dimensionless and Similar Gram Matrix in Fault Diagnosis. IEEE Trans. Ind. Inform. 2021, 18, 1061–1071. [Google Scholar] [CrossRef]
  24. Bai, Y.; Yang, J.; Wang, J.; Li, Q. Intelligent Diagnosis for Railway Wheel Flat Using Frequency-Domain Gramian Angular Field and Transfer Learning Network. IEEE Access 2020, 8, 105118–105126. [Google Scholar] [CrossRef]
  25. Bai, Y.; Yang, J.; Wang, J.; Zhao, Y.; Li, Q. Image representation of vibration signals and its application in intelligent compound fault diagnosis in railway vehicle wheelset-axlebox assemblies. Mech. Syst. Signal Process. 2020, 152, 107421. [Google Scholar] [CrossRef]
  26. Eckmann, J.-P.; Kamphorst, S.O.; Ruelle, D. Recurrence Plots of Dynamical Systems. EPL Europhys. Lett. 1987, 4, 973–977. [Google Scholar] [CrossRef]
  27. Kennel, M.B.; Brown, R.; Abarbanel, H.D.I. Determining embedding dimension for phase-space reconstruction using a geometrical construction. Phys. Rev. A At. Mol. Opt. Phys. 1992, 45, 3403–3411. [Google Scholar] [CrossRef] [PubMed]
  28. Sandler, M.; Howard, A.; Zhu, M.; Zhmoginov, A.; Chen, L.C. MobileNetV2: Inverted Residuals and Linear Bottlenecks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Salt Lake City, UT, USA, 18–23 June 2018. [Google Scholar]
  29. Yu, F.; Koltun, V. Multi-Scale Context Aggregation by Dilated Convolutions. arXiv 2016, arXiv:1511.07122. Available online: http://arxiv.org/abs/1511.07122 (accessed on 25 July 2024).
  30. Hou, Q.; Zhou, D.; Feng, J. Coordinate Attention for Efficient Mobile Network Design. In Proceedings of the 2021 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR), Nashville, TN, USA, 20–25 June 2021. [Google Scholar]
  31. Bearings Data Center. Seeded Fault Test Data, Case Western Reserve University. Available online: http://csegroups.case.edu/bearingdatacenter/pages/download-data-file (accessed on 25 July 2024).
  32. Smith, L.N. Cyclical learning rates for training neural networks. In Proceedings of the 2017 IEEE Winter Conference on Applications of Computer Vision (WACV), Santa Rosa, CA, USA, 24–31 March 2017; pp. 464–472. [Google Scholar] [CrossRef]
  33. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep Residual Learning for Image Recognition. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2015. [Google Scholar]
  34. Ma, N.; Zhang, X.; Zheng, H.-T.; Sun, J. ShuffleNet V2: Practical Guidelines for Efficient CNN Architecture Design. In Proceedings of the European Conference on Computer Vision (ECCV), Munich, Germany, 8–14 September 2018; pp. 116–131. [Google Scholar] [CrossRef]
  35. Tan, M.; Le, Q. EfficientNet: Rethinking Model Scaling for Convolutional Neural Networks. In Proceedings of the Proceedings of the 36th International Conference on Machine Learning, (PMLR), Long Beach, CA, USA, 9–15 June 2019. [Google Scholar]
  36. Huang, G.; Liu, Z.; Maaten, L.; Weinberger, K. Densely Connected Convolutional Networks. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR), Honolulu, HI, USA, 21–26 July 2017. [Google Scholar]
  37. Mukherjee, I.; Tallur, S. Light-Weight CNN Enabled Edge-Based Framework for Machine Health Diagnosis. IEEE Access 2021, 9, 84375–84386. [Google Scholar] [CrossRef]
  38. Cieslak, M.C.; Castelfranco, A.M.; Roncalli, V.; Lenz, P.H.; Hartline, D.K. t-Distributed Stochastic Neighbor Embedding (t-SNE): A tool for eco-physiological transcriptomic analysis. Mar. Genom. 2019, 51, 100723. [Google Scholar] [CrossRef]
Figure 1. The flow chart of RBP.
Figure 1. The flow chart of RBP.
Entropy 26 00675 g001
Figure 2. An RBP of diverse signals. (a) Sinusoidal signal. (b) X component of the Lorentz curve. (c) Gaussian white noise.
Figure 2. An RBP of diverse signals. (a) Sinusoidal signal. (b) X component of the Lorentz curve. (c) Gaussian white noise.
Entropy 26 00675 g002
Figure 3. Determination of the delay time and embedding dimension. (a) Delay time. (b) Embedding dimension.
Figure 3. Determination of the delay time and embedding dimension. (a) Delay time. (b) Embedding dimension.
Entropy 26 00675 g003
Figure 4. The structure of a CA.
Figure 4. The structure of a CA.
Entropy 26 00675 g004
Figure 5. The structure of the DSD-CNN.
Figure 5. The structure of the DSD-CNN.
Entropy 26 00675 g005
Figure 6. The proposed fault diagnosis framework.
Figure 6. The proposed fault diagnosis framework.
Entropy 26 00675 g006
Figure 7. The RBP of the CWRU dataset.
Figure 7. The RBP of the CWRU dataset.
Entropy 26 00675 g007
Figure 8. The RBP of the WTDS dataset.
Figure 8. The RBP of the WTDS dataset.
Entropy 26 00675 g008
Figure 9. The determination of the model’s parameters. (a) Initial learning rate. (b) Batch size.
Figure 9. The determination of the model’s parameters. (a) Initial learning rate. (b) Batch size.
Entropy 26 00675 g009
Figure 10. A visualization of the results of the ablation experiment. (a) Removal of the SeparableConv block. (b) Removal of the DepthwiseConv block. (c) Removal of the CA module. (d) The DSD-CNN.
Figure 10. A visualization of the results of the ablation experiment. (a) Removal of the SeparableConv block. (b) Removal of the DepthwiseConv block. (c) Removal of the CA module. (d) The DSD-CNN.
Entropy 26 00675 g010
Figure 11. Loss and accuracy curves of the DSD-CNN based on the WTDS dataset. (a) Training set. (b) Test set.
Figure 11. Loss and accuracy curves of the DSD-CNN based on the WTDS dataset. (a) Training set. (b) Test set.
Entropy 26 00675 g011
Figure 12. Classification results of the WTDS dataset. (a) Confusion matrix. (b) Dimensionality reduction visualization by the T-SNE.
Figure 12. Classification results of the WTDS dataset. (a) Confusion matrix. (b) Dimensionality reduction visualization by the T-SNE.
Entropy 26 00675 g012
Figure 13. The ROC curves of the DSD-CNN based on the WTDS dataset.
Figure 13. The ROC curves of the DSD-CNN based on the WTDS dataset.
Entropy 26 00675 g013
Figure 14. Comparison of the different methods on the WTDS test rig. (a) Test accuracy curve. (b) Boxplots.
Figure 14. Comparison of the different methods on the WTDS test rig. (a) Test accuracy curve. (b) Boxplots.
Entropy 26 00675 g014
Figure 15. Loss and accuracy curves of the DSD-CNN based on the CWRU dataset. (a) Training set. (b) Test set.
Figure 15. Loss and accuracy curves of the DSD-CNN based on the CWRU dataset. (a) Training set. (b) Test set.
Entropy 26 00675 g015
Figure 16. Classification results of the CWRU dataset. (a) Confusion matrix. (b) Dimensionality reduction visualization by the T-SNE.
Figure 16. Classification results of the CWRU dataset. (a) Confusion matrix. (b) Dimensionality reduction visualization by the T-SNE.
Entropy 26 00675 g016
Figure 17. Comparison of the different methods on the CWRU test rig. (a) Test accuracy curve. (b) Boxplots.
Figure 17. Comparison of the different methods on the CWRU test rig. (a) Test accuracy curve. (b) Boxplots.
Entropy 26 00675 g017
Table 1. Details of the experimental dataset selected from the CWRU data.
Table 1. Details of the experimental dataset selected from the CWRU data.
LabelFault InformationSpeedLength of SampleSize of SampleIndex
Ball1B014_01797 rpm10002400
Ball2B021_010002401
IR1IR014_010002402
IR2IR021_010002403
NN_010002404
OR1OR014@6_010002405
OR2OR021@6_010002406
Table 2. Details of the experimental dataset selected from the WTDS data.
Table 2. Details of the experimental dataset selected from the WTDS data.
LabelFault InformationLoadLength of SampleSize of SampleIndex
BT1Broken tooth 0.4 Nm12004000
BT2Broken tooth 1.2 Nm12004001
GTC1Gear tooth crack 0.4 Nm12004002
GTC2Gear tooth crack 1.2 Nm12004003
NNormal 0 Nm12004004
TSW1Tooth surface wear 0.4 Nm12004005
TSW2Tooth surface wear 1.2 Nm12004006
Table 3. The results of the ablation experiment.
Table 3. The results of the ablation experiment.
TypeAccuracyPrecisionRecallF1-ScoreTotal
Parameters
Total FLOPs
Removal of the SeparableConv block73.25%73.31%73.25%73.25%5,716,980631.29 MFLOAPs
Removal of the DepthwiseConv block87.84%89.65%87.84%87.84%6,517,728719.71 MFLOAPs
Removal of the CA module96.35%96.61%96.35%96.35%6,467,472714.15 MFLOAPs
DSD-CNN100%100%100%100%6,558,192724.18 MFLOAPs
Table 4. The classification results of the DSD-CNN on the WTDS test rig (%).
Table 4. The classification results of the DSD-CNN on the WTDS test rig (%).
LabelBT1BT2GTC1GTC2NTSW1TSW2Average
Precision100.00100.00100.00100.00100.00100.0098.7599.82
Recall100.00100.00100.0098.73100.00100.00100.0099.82
F1-score100.00100.00100.0099.36100.00100.00100.0099.82
Support79797979797979553
Table 5. Accuracy comparison of the different classification methods on the WTDS test rig.
Table 5. Accuracy comparison of the different classification methods on the WTDS test rig.
MethodAccuracy (%)Total ParametersTotal Flops
ResNet98.22 ± 0.2423,528,5224.12 GFLOAPs
EfficientNet97.81 ± 0.3710,706,991987.18 MFLOAPs
ShuffleNetV296.04 ± 0.415,359,339591.8 MFLOAPs
DenseNet98.93 ± 0.2928,681,0007.82 GFLOAPs
LW-CNN95.63 ± 0.332,397,087363.59 MFLOAPs
DSD-CNN99.78 ± 0.076,558,192724.18 MFLOAPs
Table 6. Anti-noise analysis of the different methods using the WTDS dataset (%).
Table 6. Anti-noise analysis of the different methods using the WTDS dataset (%).
MethodAverage Accuracy (4 dB)Average Accuracy (8 dB)
ResNet79.44 ± 3.1292.79 ± 1.89
EfficientNet72.29 ± 5.2691.71 ± 3.22
ShuffleNetV277.89 ± 5.0890.87 ± 4.23
DenseNet85.14 ± 3.8794.31 ± 1.83
LW-CNN75.67 ± 3.7991.29 ± 2.61
DSD-CNN87.82 ± 2.8495.17 ± 1.14
Table 7. The classification results of the DSD-CNN on the CWRU test rig (%).
Table 7. The classification results of the DSD-CNN on the CWRU test rig (%).
LabelBall1Ball2IR1IR2NOR1OR2Average
Precision100.00100.00100.00100.00100.00100.00100.00100.00
Recall100.00100.00100.00100.00100.00100.00100.00100.00
F1-score100.00100.00100.00100.00100.00100.00100.00100.00
Support47474747474747329
Table 8. Accuracy comparison of the different classification methods on the CWRU test rig.
Table 8. Accuracy comparison of the different classification methods on the CWRU test rig.
MethodAccuracy (%)Total ParametersTotal Flops
ResNet99.39 ± 0.2923,528,5224.12 GFLOAPs
EfficientNet99.57 ± 0.1610,706,991987.18 MFLOAPs
ShuffleNetV296.17 ± 0.465,359,339591.8 MFLOAPs
DenseNet99.61 ± 0.3228,681,0007.82 GFLOAPs
LW-CNN96.54 ± 0.452,397,087363.59 MFLOAPs
DSD-CNN100.00 ± 0.006,558,192724.18 MFLOAPs
Table 9. Anti-noise analysis of the different methods on the WTDS dataset (%).
Table 9. Anti-noise analysis of the different methods on the WTDS dataset (%).
MethodAverage Accuracy (4 dB)Average Accuracy (8 dB)
ResNet73.43 ± 2.4591.89 ± 2.82
EfficientNet70.92 ± 4.7792.25 ± 1.93
ShuffleNetV274.27 ± 3.3990.57 ± 4.21
DenseNet86.58 ± 3.1595.01 ± 2.18
LW-CNN74.43 ± 4.2591.12 ± 2.64
DSD-CNN86.97 ± 2.6495.45 ± 1.72
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shi, Y.; Wang, H.; Sun, W.; Bai, R. Intelligent Fault Diagnosis Method for Rotating Machinery Based on Recurrence Binary Plot and DSD-CNN. Entropy 2024, 26, 675. https://doi.org/10.3390/e26080675

AMA Style

Shi Y, Wang H, Sun W, Bai R. Intelligent Fault Diagnosis Method for Rotating Machinery Based on Recurrence Binary Plot and DSD-CNN. Entropy. 2024; 26(8):675. https://doi.org/10.3390/e26080675

Chicago/Turabian Style

Shi, Yuxin, Hongwei Wang, Wenlei Sun, and Ruoyang Bai. 2024. "Intelligent Fault Diagnosis Method for Rotating Machinery Based on Recurrence Binary Plot and DSD-CNN" Entropy 26, no. 8: 675. https://doi.org/10.3390/e26080675

APA Style

Shi, Y., Wang, H., Sun, W., & Bai, R. (2024). Intelligent Fault Diagnosis Method for Rotating Machinery Based on Recurrence Binary Plot and DSD-CNN. Entropy, 26(8), 675. https://doi.org/10.3390/e26080675

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

Article Metrics

Back to TopTop