Next Article in Journal
Real-Time Object Detection for Autonomous Solar Farm Inspection via UAVs
Next Article in Special Issue
Off-Design Operation and Cavitation Detection in Centrifugal Pumps Using Vibration and Motor Stator Current Analyses
Previous Article in Journal
Remote Heart Rate Estimation Based on Transformer with Multi-Skip Connection Decoder: Method and Evaluation in the Wild
Previous Article in Special Issue
Scale-Fractal Detrended Fluctuation Analysis for Fault Diagnosis of a Centrifugal Pump and a Reciprocating Compressor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Intelligent Ball Bearing Fault Diagnosis System Using Enhanced Rotational Characteristics on Spectrogram

Department of Electronics Engineering, Dong-A University, Busan 49315, Republic of Korea
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(3), 776; https://doi.org/10.3390/s24030776
Submission received: 16 November 2023 / Revised: 4 January 2024 / Accepted: 22 January 2024 / Published: 25 January 2024
(This article belongs to the Special Issue Fault Diagnosis and Vibration Signal Processing in Rotor Systems)

Abstract

:
Faults in the ball bearing are a major cause of failure in rotating machinery where ball bearings are used. Therefore, there is a growing demand for ball bearing fault diagnosis to prevent failures in rotating machinery. Although studies on the fault diagnosis of bearing have been conducted using temperature measurements and sound monitoring, these methods have limitations, because they are affected by external noise. Therefore, many researchers have studied vibration monitoring for bearing fault diagnosis. Among these, mel-frequency cepstral coefficients (MFCCs) and 2D convolutional neural networks (CNNs) have attracted significant attention in vibration monitoring schemes. However, the MFCC in existing studies requires a high sampling rate and an expansive frequency band utilization. In addition, 2D CNNs are highly complex. In this study, a rotational characteristic emphasis (RCE) spectrogram process and an optimized CNN were proposed to solve these problems. The RCE spectrogram process analyzes a narrow frequency band and produces low-resolution images. The optimized CNN was designed with a shallow network structure. The experimental results showed an accuracy of 0.9974 for the proposed system. The optimized CNN model has parameters of 5.81 KB and FLOPs of 1.53 × 10 6 . We demonstrate that the proposed ball bearing fault diagnosis system can achieve high accuracy with low complexity. Thus, we propose a ball bearing fault diagnosis scheme that is applicable to a low sampling rate and changing rotation frequency.

1. Introduction

Ball bearings are responsible for reducing rotational friction and supporting radial and axial loads in rotating machinery. However, they are the components that most commonly develop faults in rotating machinery [1,2]. Therefore, ball bearing faults can cause severe damage to rotating machinery, resulting in economic losses and increased safety risks [3,4,5,6]. Over the past few decades, many studies have been conducted on ball bearing fault diagnosis [7,8,9,10,11,12,13,14,15,16]. Ball bearing fault diagnosis methods can be categorized into temperature measurement, sound monitoring, and vibration monitoring. Temperature measurement utilizes 2D ball bearing images that are acquired using an infrared thermal imaging camera to diagnose ball bearing faults, whereas sound monitoring utilizes ball bearing sound signals that are acquired using a microphone to diagnose the faults. However, sensors used in conventional methods have certain limitations. In the case of infrared thermal imaging cameras, dust and air particles in the field can interfere with the image information and distort readings. In the case of microphones, external noise distorts the signal, and sound attenuates over the measurement distance. However, accelerometers are robust to external factors, because they measure vibration signals when in contact with the rotating machinery. Table 1 summarizes the sensors and the limitations of the measurement methods that are used for bearing fault diagnosis.
Convolutional neural networks (CNNs) have powerful nonlinear feature extraction capabilities, so many researchers have recently applied CNNs to classify imaged bearing vibration signals [17]. Imaged bearing vibration signals include time and frequency information. In particular, the short-time Fourier transform (STFT), wavelet transform (WT), and Hilbert–Huang transform (HHT) are used to analyze vibration signals in terms of time and frequency. Among them, the STFT is a simple method for converting vibration signals into time and frequency. Therefore, the STFT is more efficient than the WT and HHT [18]. Shan et al. [19] proposed Mel-CNN, which combines an MFCC image and a CNN by utilizing sound signals. The Mel spectrogram effectively extracted bearing fault features due to the low correlation between vectors, and the CNN was combined with this for image classification to achieve effective fault classification. However, the sound signal contains not only bearing fault features but also a large amount of white noise. Furthermore, the noise in real industrial plants is critical for ball bearing fault classification. Xia et al. [20] proposed a CNN model for fault diagnosis during vibration monitoring. This method used a method that arranged the collected multichannel vibration signals in a row and preserved the time and space information of the original signals, so that the rich information effectively diagnosed the fault. However, existing CNNs are limited in their application to fault diagnosis because they require numerous training samples to achieve high diagnostic accuracy. In [21], a bearing fault diagnosis method based on Mel frequency cepstral coefficients (MFCCs) and a CNN was proposed. First, an MFCC was introduced to extract the signature features of the fault signal. Subsequently, a CNN with deformable convolutional blocks was used to extract spatial information between categories and improve the extracted spatial frequency features. This shows that the MFCC can be used as a tool for feature extraction in bearing fault diagnosis. The MFCC is effective for the feature extraction of speech signals. However, the MFCC used in existing studies requires a high sampling rate and utilizes a wide frequency band, even when using vibrating signals [22,23,24]. In addition, the CNN increases the complexity of the fault classification process [25,26,27]. Within the same context, [28] extracts the characteristic fault frequencies from the resonant frequency band of the vibration signals. The resonant frequency band is a relatively higher frequency band than the ball bearing’s fault frequency.
To overcome the limitations of existing research, we proposed a ball bearing failure diagnosis system, which consists of two steps: feature extraction and fault classification. The feature extraction method for the proposed system is a rotational characteristic emphasis (RCE) spectrogram, which modifies the MFCC process to emphasize the ball bearing’s fault features. In particular, the Mel scale was modified to improve the resolution of the fault frequency. In addition, the discrete cosine transform (DCT) was not used in the classification, because it reduces feature information and requires additional computation. Next, the fault classification of the proposed system was an optimized CNN, which uses the RCE spectrogram as input image. Furthermore, the optimized CNN was designed to fit the size of the RCE spectrogram and optimized for effective classification. The ball bearing fault diagnosis system is illustrated in Figure 1. The contributions of this study are as follows:
  • The proposed system extracts the fault features of ball bearings that are related to rotational frequency. This method does not require a high sampling rate for the accelerometer when collecting vibration signals.
  • Compared to existing studies that use MFCCs to extract the features of vibration signals, the RCE spectrogram process analyzes a narrow frequency band. This requires fewer resources than existing feature extraction methods. In addition, the optimized CNN architecture for fault classification has lower complexity than existing CNN architectures.
  • The Case Western Reserve University (CWRU) dataset was used to evaluate the ball bearing fault diagnosis. The experimental results indicated that the proposed method classified ball bearing faults with an accuracy of 0.9974. This satisfies the conditions for a ball bearing fault diagnosis system.
The remainder of this paper is organized as follows: The background for understanding ball bearing fault characteristics is explained in Section 2. The proposed ball bearing fault diagnosis system is described in Section 3. In Section 4, we describe an experiment using the proposed system and analyze the experimental results. Finally, the conclusions are presented in Section 5.

2. Background Knowledge

Ball bearing fault frequencies occur periodically when a component of a ball bearing is faulty. The ball pass frequency of the outer race (BPFO) is observed when the outer race of a ball bearing is faulty. The ball passes through the fault in the outer race, causing a periodic vibration. The ball pass frequency of the inner race (BPFI) is observed when a fault occurs in the inner race of a ball bearing. The ball spins and passes through the fault in the inner race, causing a periodic vibration. The ball spin frequency (BSF) is observed when a fault occurs in the balls of a ball bearing. As the balls pass through the bearing, they spin around their own axis and generate periodic vibrations. All of these cause periodic vibrations, resulting in a peak being observed in the frequency domain. The fault frequency for different types of faults can change depending on the specifications of the ball bearing and the rotational frequency of the rotating machinery. Equations (1)–(3) show the fault frequencies according to the types of ball bearing faults:
B P F O = n f r 2 1 + d D c o s α
B P F I = n f r 2 1 d D c o s α
B S F = D f r 2 d 1 d D c o s α 2
where D , d , a , n , and f r are the pitch diameter, ball diameter, contact angle between the ball and the cage, number of rolling elements, and the rotating speed of the bearing (Hz), respectively. For the ball bearing fault diagnosis dataset, we used the CWRU dataset, which contains vibration signals from ball bearings that are used when operating rotating machines [29]. The vibration signals were collected using an accelerometer attached to the rotating machine. Table 2 lists the specifications of the CWRU dataset that is used to determine the fault frequency.
The rotational frequency multiples (in Hz) in Table 3 are multiples and have no units, so Table 3 shows that the fault frequency can be detected within six times the rotational frequency. The rotational frequencies in the CWRU dataset are 29.95, 29.53, 29.17, and 28.83 Hz, with a maximum fault frequency of about 162 Hz. Bearing fault frequencies typically occur in the frequency band below 500 Hz [30,31]. This means that a sampling rate of more than 1 kHz is required for the detection of a ball bearing fault. The MFCC has been proven to be effective in extracting ball bearing fault features by existing studies. However, these existing studies require a frequency band of at least 8 kHz, based on the fact that humans can determine bearing faults by sound. Therefore, a sampling rate of at least 16 kHz is required [22,32]. We prove that ball bearing fault diagnosis is possible by effectively extracting the fault frequency features of vibration signals using frequency bands within 500 Hz.

3. Proposed Ball Bearing Fault Diagnosis System

3.1. Feature Extraction Method

The down-sampled vibration signals were used as inputs for the feature extraction method. The RCE spectrogram, which is the feature extraction method used in our proposed system, enhances the faulty features of the ball bearing and converts it into a 2D image. The RCE spectrogram process is shown in Figure 2.
The down-sampled vibration signal was filtered with a pre-emphasis to balance the frequency spectrum. If the pre-emphasis filtered vibration signal is x 1 [ n ] , the down-sampled vibration signal is x [ n ] , and the pre-emphasis filter coefficient is α . The pre-emphasis filter is then as follows:
x p r e n = x n α x n 1   ,             n 1   ,             n 1 ,   2 ,   3 ,  
Next, we used the STFT to generate an image containing time and frequency information. The STFT result, pre-emphasis filtered vibration signals, and window function are x ( m , w ) , x p r e [ n ] , and w [ n ] , respectively. The STFT result is then as follows:
x ( m , w ) = n x p r e n   w n m exp j w n
where m is the time index, w is the frequency index, e x p ( j w p ) is the complex exponential function, and j is the imaginary unit. The frame length was set to 100 samples, and each frame was set to overlap by 50 samples, resulting in 19 frames. The number of overlapping samples in each frame is related to the features of the STFT. When the STFT spectrum has a large overlapping area in a frame, there is less change or distortion, but more computation. To minimize the amount of computation without losing fault features from windowing, we overlapped 50 vibration signal samples, which is half the number of frame samples. Each frame uses a fast Fourier transform (FFT) for frequency analysis. The results of the FFT are filtered using an RCE filter bank to generate a spectrogram.
The RCE filter bank is a triangular filter bank that is generated based on the RCE scale. For the RCE filter bank to detect ball bearing faults, the RCE scale was designed to have a high resolution at the fault frequency of the ball bearing. The RCE scale used to generate the RCE filter bank is expressed as follows:
R f = 6 a · exp f 10 c + 1 1 exp 3 a 10 c + 1 1 · exp 3 a 10 c + 1 exp f 10 c + 1 + exp 3 a 10 c + 1     ,             0 f 6 a  
where R ( f ) is the RCE scale, f is the input frequency, a is the rotational frequency, and c is a slope factor. The RCE scale decreased as the slope shifted from the center frequency. To control the slope, the slope factor was set to between 0 and 1. When the slope factor was zero, the RCE scale had the highest slope, and when the slope factor was one, the RCE scale approximated a first-order function. This allowed the RCE filter bank to control the resolution in the frequency band. The rotational frequency was set as a parameter in the RCE scale to allow the RCE filter bank to be used on rotational machinery at any speed.
The RCE spectrogram was generated with dimensions of 19 × 19 × 3 . Therefore, an RCE filter bank with a size of 19 was required. One triangular filter has two different lines. To generate the two different lines, we need three points: a start point, a center point, and an end point. If you were to arrange these triangular filters, the start point of the next triangular filter becomes the center point of the first triangular filter, and the center point becomes the end point. Each point is an RCE scale value, so the total number of RCE scale values needed should be added by 2 to the number of triangular filters. If the number of RCE scale values is n , i takes the values 1 ,   2 ,   3 ,   ,   n . We called the set of each RCE scale value an RCE point and defined the RCE points as R [ i ] . The following process is required to generate the RCE filter bank. First, the rotational frequency must be six times that of the 21 equal parts to obtain the values. Each set of values is called a frequency point, which is defined as F [ i ] . Next, F [ i ] is substituted into the RCE scale to obtain the RCE points. Three consecutive values of RCE points generate one triangular filter. This process is repeated to form the RCE filter bank, which can be defined as H R ( x ) . The frequency points, RCE points, and RCE filter bank are expressed as follows:
F i = R ( 6 a ) max i 1 i 1 ,             i 1 ,           i 1 ,   2 ,   3 , ,   n
R i = 6 a · exp F i 10 c + 1 1 exp 3 a 10 c + 1 1 · exp 3 a 10 c + 1 exp F i 10 c + 1 + exp 3 a 10 c + 1      
H R ( x ) =               0 ,                                 x < R i                                                         x     R ( i ) R i   +   1     R ( i ) ,         R i x R i   +   1         R i   +   2     x R i   +   2     R ( i   +   1 ) , R i + 1 x R i + 2             0 ,                                     x > R i + 2                  
Figure 3 shows the RCE scale and RCE filter bank when the rotational frequency is 29.95 Hz (1797 rpm) and the slope factor is 0, 0.3, 0.6, and 0.9.

3.2. Fault Classification

The optimized CNN was designed to achieve accurate image classification with low computational complexity. Table 4 provides detailed information regarding the layers of the CNN model, whereas Figure 4 provides a visual representation of the optimized configuration of the CNN model.
The goal of the optimized CNN is to minimize computational costs, enhance generalization, and reduce the number of features. To achieve this goal, the optimized CNN uses small-sized RCE spectrogram images as input and consists of a cross-channel normalization layer and a global average pooling layer. The small-sized input image simplifies the convolutional operations, and the cross-channel normalization layer aggregates the information of adjacent maps and improves the generalization of the time–frequency spatial information obtained from the STFT [33]. The global average pooling layer improves the correspondence between feature maps and classes and significantly reduces the number of feature maps. Consequently, feature maps can be interpreted as confidence maps, and overfitting is mitigated, because no additional parameters require optimization [34]. Our approach involves a fully connected layer and softmax function to determine the probability distribution for classification, incorporating dropout to enhance generalization. For layers requiring an activation function, the clipped rectified linear unit (ReLU) activation function combines the advantages of the ReLU activation function with reduced computational costs and addresses issues relating to exploding activations [21,24,35].
The optimized CNN consists of an input layer, convolution layer, normalization layer, pooling layer, and fully connected layer. In the optimized CNN, blocks are not reused to achieve low complexity. Initially, the input image of the input layer is an RCE spectrogram image, so the size of the input image is 19 × 19 × 3. In the convolution layer, a 7 × 7 × 32 kernel is used, and the stride is set to 2 and the padding to 3. In the normalization layer, we use cross-channel normalization and generalize between each channel. The window channel size is 5. In the pooling layer, we use global average pooling. After that, we perform classification in the fully connected layer and obtain the classification result.
We chose a strategy for improving generalization and reducing the number of feature maps by following the structure of the existing 2D-CNN while reducing the complexity. To apply this, we used cross-channel normalization in the normalization layer and global average pooling in the pooling layer to reduce the number of feature maps.

4. Experiment

4.1. Experimental Setup and Dataset Description

The proposed diagnosis system was evaluated using the CWRU dataset as input. The dataset consists of ball bearings that are attached to the drive and fan ends and is measured using accelerometers. The classes for the ball bearing fault classification are assigned to the vibration signals in the normal state as the normal class and the vibration signals of the outer ring fault, inner ring fault, and ball fault as fault classes. The vibration signal had a sampling rate of 48 kHz and was provided as a MATLAB extension mat file. The vibration signal was down-sampled to a sampling rate of 1 kHz to analyze the frequency band within 500 Hz. The down-sampled vibration signal was collected from a random starting point of the 1000th sample, and the starting point was not repeated. This collection method is based on collecting ball bearing vibration signals from a plant using a ball bearing fault diagnosis system. Table 5 lists the acquisition conditions and classification details from the CWRU dataset.
The identification of ball bearings according to their faults helps companies that produce and repair ball bearings analyze future ball bearing faults. However, industrial plants that actually use ball bearings need to be able to diagnose quickly to prevent rotor failure when a ball bearing fails, and ball bearings need to be replaced regardless of the type of fault [36,37]. For this reason, in this paper, we categorize them into two classes: normal and fault. In addition, the data included in the fault class include the fault data of the outer race, inner race, and ball. The CWRU dataset has three fault types, inner ring fault, outer ring fault, and ball fault, and each fault type has three fault diameters. In this paper, we use all the provided fault datasets as fault classes.
The collected ball bearing vibration signal was converted into an RCE spectrogram with dimensions of 19 × 19 × 3 , which was then used as input for training the optimized CNN. The total number of RCE spectrogram images was 6840, of which 4653 were used as training images, 1024 as validation images, and 1163 as test images. The parameters used to train the optimized CNN are listed in Table 6.

4.2. Evaluation Details

The accuracy, recall, precision, and F1 score were selected to evaluate the images of the ball bearing fault diagnosis system, and the parameters and floating-point operations (FLOPs) were selected to evaluate the complexity of the CNN architecture. FLOPs measure the total number of operations that are performed within a CNN model, representing the computational workload. The calculation of FLOPs was based on the layer structure of the CNN model. However, the parameters indicate the number of learnable parameters that are involved in the CNN architecture. First, the accuracy, recall, precision, and F1 scores were calculated using the following equations:
A c c u r a c y = T P + T N T P + F N + F P + T N
R e c a l l = T P T P + F N
P r e c i s i o n = T P T P + F P
F 1   s c o r e = 2 · P r e c i s i o n · R e c a l l P r e c i s i o n + R e c a l l
where TP and TN denote the number of positive and negative instances that are correctly classified. Also, FP and FN denote the number of negative and positive instances that are not correctly classified. In detail, TP is the result of predicting a defective image as a fault, and FP is the result of predicting a fault image as normal [38]. Here, positive instances refer to the faulty state of the ball bearing and negative instances refer to the normal state of the ball bearing. To calculate the complexity, the parameters and FLOPs were calculated using the following equations:
C o n v o l u t i o n a l   l a y e r   F L O P s = K 2 × C i n × L o u t 2
F u l l y   c o n n e c t e d   l a y e r   F L O P s = ( 2 N i n 1 ) × N o u t
P a r a m e t e r s = K 2 × C i n × C o u t + B
Equations (14) and (15) are used to calculate the FLOPs, and Equation (16) is used to calculate the number of parameters. In these equations, K is the kernel size, C i n is the number of input channels, C o u t is the number of kernels, L o u t is the output image size, N i n and N o u t are the numbers of input and output nodes, and B is the number of biases.

4.3. Analysis Result

To evaluate whether the image size and slope factor in the proposed model are appropriate for the CWRU dataset, we check the accuracy as a function of the slope factor and image size. Table 7 shows the accuracy and loss of the slope factor, and Table 8 shows the accuracy and loss of the image size. Through the analysis, the accuracy and loss of the ball bearing fault classification were obtained. The results indicated that the highest accuracy and loss accuracy and the loss of the test data were 0.9974 and 0.0182, respectively. Figure 5 shows the training process that were used to obtain the classification accuracy and the loss of training and test data of the proposed model, whereas Table 9 shows the accuracy and loss in the training and test. Table 10 shows the confusion matrix of the test data of the proposed model. Since the training and test data are changed from training to training, the accuracy and loss of the proposed technique shows robustness.
Next, the parameters and FLOPs of the optimized CNN were calculated using Equations (14)–(16). As a result of the calculation, the optimized CNN has a parameter of 5.81 KB and a FLOP of 1.53 × 10 6 .

4.4. Comparison of the Results with Other Methods

The RCE spectrogram was compared with another spectrogram, the Mel spectrogram, and the MFCC. The results demonstrated that the RCE spectrogram is effective for bearing fault diagnosis. Each image was trained with the optimized CNN and evaluated in terms of accuracy, recall, precision, and F1 score. The evaluation was conducted under identical conditions. The number of training images was 4660, and the size of each image was set to 19 × 19 × 3 . Table 11, Table 12 and Table 13 show the confusion matrices of the spectrogram, Mel spectrogram, and MFCC.
Table 11, Table 12 and Table 13 show that the RCE spectrogram provides the highest classification accuracy when the image size is small. In addition, the Mel spectrogram and MFCC are unable to classify the normal class. We calculated the metric of each confusion matrix, and the results are shown in Table 14.
The CNNs were also compared in terms of the complexity and accuracy of the input images. Table 15 shows the complexities of the proposed CNN and CNNs that are used for ball bearing fault diagnosis in a related work [39].
The existing 2D CNN model has the disadvantages of high memory utilization and long computation time owing to the high number of parameters and FLOPs. The proposed optimized CNN is a 2D CNN model that uses small-sized input images and shallow depth layers. Therefore, the limitations of the high number of parameters and FLOPs of the 2D CNN model were improved.

5. Conclusions

The applications of noncontact bearing diagnosis methods are limited in real industrial plants owing to the influence of external noise. Therefore, this study proposes a ball bearing fault diagnosis system that combines the RCE spectrogram process and an optimized CNN using the contact measurement of vibration signals. The system down-samples the vibration signal to 1 kHz and then applies the RCE spectrogram process to generate an image with a size of 19 × 19 × 3 . The generated image is used as an input to the optimized CNN to classify faults in the ball bearing. The feature extraction of the RCE spectrogram and the complexity of the CNN were evaluated. The evaluation results showed that the accuracy of the RCE spectrogram is 99.74%, which is the highest accuracy among the generated image sets, and the optimized CNN has 5.81 KB of parameters and FLOPs of 1.53 × 10 6 , despite being a 2D CNN model. In conclusion, our evaluation demonstrated that the problems of a high sampling rate, high utilization of existing wideband frequency analysis, and the high complexity of deep neural network models are solved. In particular, the RCE scale in the RCE spectrogram process can control the band of the rotational frequency and the resolution of the RCE filter bank by changing the slope factor. Therefore, the RCE filter bank can be applied even when the rotation frequency of the rotating machine changes. Furthermore, it can be used to analyze ball bearing faults in various rotating machineries.
In our future studies, we plan to design a prototype module with a ball bearing defect system. The prototype module will be used to apply the ball bearing fault system to a practical environment. Furthermore, we will compare and verify the performance of the ball bearing fault diagnosis system through experiments. Based on the experimental results, we will apply an algorithm so that the slope factor of the RCE scale can change adaptively.

Author Contributions

G.S. contributed to the methodology, writing, and visualization of this study; performed the implementation and testing; and validated the result; D.K. contributed to the supervision, project administration, and funding acquisition of the research and devised the basic concept of the proposed scheme. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Dong-A University research fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Grebenik, J.; Zhang, Y.; Bingham, C.; Srivastava, S. Roller Element Bearing Acoustic Fault Detection Using Smartphone and Consumer Microphones Comparing with Vibration Techniques. In Proceedings of the 17th International Conference on Mechatronics—Mechatronika (ME), Prague, Czech Republic, 7–9 December 2016; pp. 1–7. [Google Scholar]
  2. Cerrada, M.; Sánchez, R.V.; Li, C.; Pacheco, F.; Cabrera, D.; de Oliveira, J.V.; Vásquez, R.E. A review on data-driven fault severity assessment in rolling bearings. Mech. Syst. Signal Process. 2018, 99, 169–196. [Google Scholar] [CrossRef]
  3. Abdenebi, R.; Fethi, D.; Abdelkrim, N.; Badis, D. Gearbox Fault Diagnosis Using the Short-Time Cepstral Features. In Proceedings of the 2022 2nd International Conference on Advanced Electrical Engineering, ICAEE 2022, Constantine, Algeria, 29–31 October 2022; pp. 1–7. [Google Scholar] [CrossRef]
  4. Wang, Y.; Markert, R.; Xiang, J.; Zheng, W. Research on variational mode decomposition and its application in detecting rub-impact fault of the rotor system. Mech. Syst. Signal Process. 2015, 60–61, 243–251. [Google Scholar] [CrossRef]
  5. Saravanan, N.; Cholairajan, S.; Ramachandran, K.I. Vibration-based fault diagnosis of spur bevel gear box using fuzzy technique. Expert Syst. Appl. 2009, 36, 3119–3135. [Google Scholar] [CrossRef]
  6. Kumar, A.; Vashishtha, G.; Gandhi, C.P.; Zhou, Y.; Glowacz, A.; Xiang, J. Novel Convolutional Neural Network (NCNN) for the Diagnosis of Bearing Defects in Rotary Machinery. IEEE Trans. Instrum. Meas. 2021, 70, 3510710. [Google Scholar] [CrossRef]
  7. 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]
  8. Younus, A.M.D.; Yang, B. Intelligent fault diagnosis of rotating machinery using infrared thermal image. Expert Syst. Appl. 2012, 39, 2082–2091. [Google Scholar] [CrossRef]
  9. Janssens, O.; Schulz, R.; Slavkovikj, V.; Stockman, K.; Loccufier, M.; Van de Walle, R.; Van Hoecke, S. Thermal image based fault diagnosis for rotating machinery. Infrared Phys. Technol. 2015, 73, 78–87. [Google Scholar] [CrossRef]
  10. Liu, Q.; Lu, X.; He, Z.; Zhang, C.; Chen, W.S. Deep convolutional neural networks for thermal infrared object tracking. Knowl. Based Syst. 2017, 134, 189–198. [Google Scholar] [CrossRef]
  11. Glowacz, A. Fault diagnosis of single-phase induction motor based on acoustic signals. Mech. Syst. Signal Process. 2019, 117, 65–80. [Google Scholar] [CrossRef]
  12. Marin, F.B.; Solomon, C.; Marin, M. Bearing Failure Prediction Using Audio Signal Analysis Based on SVM Algorithms. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Galati, Romania, 11–13 October 2018; IOP Publishing: Bristol, UK, 2019; p. 012012. [Google Scholar] [CrossRef]
  13. Nirwan, N.W.; Ramani, H.B. Condition monitoring and fault detection in roller bearing used in rolling mill by acoustic emission and vibration analysis. Mater. Today Proc. 2022, 51, 344–354. [Google Scholar] [CrossRef]
  14. Wang, J.; He, Q.; Kong, F. Adaptive Multiscale Noise Tuning Stochastic Resonance for Health Diagnosis of Rolling Element Bearings. IEEE Trans. Instrum. Meas. 2015, 64, 564–577. [Google Scholar] [CrossRef]
  15. Wu, S.; Gebraeel, N.; Lawley, M.A.; Yih, Y. A Neural Network Integrated Decision Support System for Condition-Based Optimal Predictive Maintenance Policy. IEEE Trans. Syst. Man Cybern. Part A Syst. Hum. 2007, 37, 226–236. [Google Scholar] [CrossRef]
  16. Zhou, Y.; Wang, H.; Wang, G.; Kumar, A.; Sun, W.; Xiang, J. Semi-Supervised Multiscale Permutation Entropy-Enhanced Contrastive Learning for Fault Diagnosis of Rotating Machinery. IEEE Trans. Instrum. Meas. 2023, 72, 3525610. [Google Scholar] [CrossRef]
  17. Chen, Z.; Cen, J.; Xiong, J. Rolling Bearing Fault Diagnosis Using Time-Frequency Analysis and Deep Transfer Convolutional Neural Network. IEEE Access 2020, 8, 150248–150261. [Google Scholar] [CrossRef]
  18. Chen, J.; Jiang, J.; Guo, X.; Tan, L. An Efficient CNN with Tunable Input-Size for Bearing Fault Diagnosis. Int. J. Comput. Intell. Syst. 2021, 14, 625–634. [Google Scholar] [CrossRef]
  19. Shan, S.; Liu, J.; Wu, S.; Shao, Y.; Li, H. A motor bearing fault voiceprint recognition method based on Mel-CNN model. Measurement 2023, 207, 112408. [Google Scholar] [CrossRef]
  20. 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. 2018, 23, 101–110. [Google Scholar] [CrossRef]
  21. Zhao, Y.; Zhang, N.; Zhang, Z.; Xu, X. Bearing Fault Diagnosis Based on Mel Frequency Cepstrum Coefficient and Deformable Space-Frequency Attention Network. IEEE Access 2023, 11, 34407–34420. [Google Scholar] [CrossRef]
  22. Choudakkanavar, G.; Mangai, J.A.; Bansal, M. MFCC based ensemble learning method for multiple fault diagnosis of roller bearing. Int. J. Inf. Technol. 2022, 14, 2741–2751. [Google Scholar] [CrossRef]
  23. Randall, R.B.; Antoni, J. Rolling element bearing diagnostics—A tutorial. Mech. Syst. Signal Process. 2011, 25, 485–520. [Google Scholar] [CrossRef]
  24. Akpudo, U.E.; Hur, J.W. A Cost-Efficient MFCC-Based Fault Detection and Isolation Technology for Electromagnetic Pumps. Electronics 2021, 10, 439. [Google Scholar] [CrossRef]
  25. Pham, M.T.; Kim, J.M.; Kim, C.H. Accurate bearing fault diagnosis under variable shaft speed using convolutional neural networks and vibration spectrogram. Appl. Sci. 2020, 10, 6385. [Google Scholar] [CrossRef]
  26. Zilong, Z.; Wei, Q. Intelligent Fault Diagnosis of Rolling Bearing Using One-Dimensional Multi-Scale Deep Convolutional Neural Network Based Health State Classification. In Proceedings of the 2018 IEEE 15th International Conference on Networking, Sensing and Control (ICNSC), Zhuhai, China, 27–29 March 2018; pp. 1–6. [Google Scholar] [CrossRef]
  27. Liu, R.; Meng, G.; Yang, B.; Sun, C.; Chen, X. Dislocated time series convolutional neural architecture: An intelligent fault diagnosis approach for electric machine. IEEE Trans. Ind. Inform. 2017, 13, 1310–1320. [Google Scholar] [CrossRef]
  28. Zhen, D.; Li, D.; Feng, G.; Zhang, H.; Gu, F. Rolling bearing fault diagnosis based on VMD reconstruction and DCS demodulation. Int. J. Hydromechatron. 2022, 5, 205–225. [Google Scholar] [CrossRef]
  29. Case Western Reserve University Bearing Dataset. Available online: https://engineering.case.edu/bearingdatacenter (accessed on 18 January 2023).
  30. Tandon, N.; Choudhury, A. A review of vibration and acoustic measurement methods for the detection of defects in rolling element bearings. Tribol. Int. 1999, 32, 469–480. [Google Scholar] [CrossRef]
  31. Tandon, N.; Nakra, B.C. Technical Article. Practical articles in shock and vibration technology: Vibration and Acoustic Monitoring Techniques for the Detection of Defects in Rolling Element Bearings—A Review. Shock Vib. Dig. 1992, 24, 3–11. [Google Scholar] [CrossRef]
  32. Jiang, Q.; Chang, F.; Sheng, B. Bearing fault classification based on convolutional neural network in noise environment. IEEE Access 2019, 7, 69795–69807. [Google Scholar] [CrossRef]
  33. Lu, S.; Lu, Z.; Zhang, Y.D. Pathological brain detection based on AlexNet and transfer learning. J. Comput. Sci. 2019, 30, 41–47. [Google Scholar] [CrossRef]
  34. Lin, M.; Chen, Q.; Yan, S. Network in Network. arXiv 2013, arXiv:1312.4400. [Google Scholar]
  35. Wang, S.H.; Muhammad, K.; Hong, J.; Sangaiah, A.K.; Zhang, Y.D. Alcoholism identification via convolutional neural network based on parametric ReLU, dropout, and batch normalization. Neural Comput. Appl. 2018, 32, 665–680. [Google Scholar] [CrossRef]
  36. Shreve, D.H. Introduction to Vibration Technology. In Proceedings of the Predictive Maintenance Technology Conference, Annual Meeting, Columbus, OH, USA, 16 November 1994. [Google Scholar]
  37. Rohani, R.; Jafari, S.M.; Roozban, M. Study of ball bearings failure modes in an eddy current dynamometer. J. Engine Res. 2016, 41, 3–11. [Google Scholar]
  38. Šimundić, A.M. Measures of Diagnostic Accuracy: Basic Definitions. EJIFCC 2009, 19, 203–211. [Google Scholar] [PubMed]
  39. Ma, S.; Liu, W.; Cai, W.; Shang, Z.; Liu, G. Lightweight deep residual CNN for fault diagnosis of rotating machinery based on depthwise separable convolutions. IEEE Access 2019, 7, 57023–57036. [Google Scholar] [CrossRef]
  40. Hakim, M.; Omran, A.A.B.; Inayat-Hussain, J.I.; Ahmed, A.N.; Abdellatef, H.; Abdellatif, A.; Gheni, H.M. Bearing Fault Diagnosis Using Lightweight and Robust One-Dimensional Convolution Neural Network in the Frequency Domain. Sensors 2022, 22, 5793. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic diagram of the proposed ball bearing fault diagnosis system.
Figure 1. Schematic diagram of the proposed ball bearing fault diagnosis system.
Sensors 24 00776 g001
Figure 2. Schematic diagram of the RCE spectrogram process.
Figure 2. Schematic diagram of the RCE spectrogram process.
Sensors 24 00776 g002
Figure 3. RCE scale and RCE filter bank depending on the slope factor: (a) c = 0 ; (b) c = 0.3 ; (c) c = 0.6 ; and (d) c = 0.9 .
Figure 3. RCE scale and RCE filter bank depending on the slope factor: (a) c = 0 ; (b) c = 0.3 ; (c) c = 0.6 ; and (d) c = 0.9 .
Sensors 24 00776 g003
Figure 4. Schematic diagram of the optimized CNN.
Figure 4. Schematic diagram of the optimized CNN.
Sensors 24 00776 g004
Figure 5. Classification training process: (a) accuracy; (b) loss.
Figure 5. Classification training process: (a) accuracy; (b) loss.
Sensors 24 00776 g005
Table 1. Limitations of sensor types for bearing fault diagnosis.
Table 1. Limitations of sensor types for bearing fault diagnosis.
MethodSensorLimitation
Temperature measurement [7,8,9,10]Infrared thermal
imaging camera
Distortion in readings
from dust and air particles
Sound monitoring [11,12,13]MicrophonesExternal noise
and signal attenuation
Table 2. Specifications of the ball bearings from the CWRU dataset.
Table 2. Specifications of the ball bearings from the CWRU dataset.
Specification6205-2RS JEM SKF
(Drive-End Bearing)
6203-2RS JEM SKF
(Fan-End Bearing)
Inside diameter [mm]25.017.0
Outside diameter [mm]52.040.0
Pitch diameter [mm]39.028.5
Ball diameter [mm]7.946.75
Thickness [mm]15.012.0
Number of balls98
Contact angle [ ° ] 3.1343.126
Table 3. Fault frequencies according to the ball bearing fault type from the CWRU dataset.
Table 3. Fault frequencies according to the ball bearing fault type from the CWRU dataset.
Types of Ball Bearing Fault6205-2RS JEM SKF
(Drive-End Bearing)
6203-2RS JEM SKF
(Fan-End Bearing)
Multiple of Rotational Frequency (in Hz)
Inner ring5.41524.9469
Outer ring3.58483.0530
Rolling element (ball)4.71353.9874
Table 4. Details of the optimized CNN architecture.
Table 4. Details of the optimized CNN architecture.
LayerData SizeKernel SizeActivation Function
Input data 19 × 19 × 3 Not neededNot needed
Convolutional 10 × 10 × 32 7 × 7 × 32 Clipped rectified linear unit
Cross-channel normalizationWindow channel size = 5Not needed
Global average pooling Not needed
Fully connected32 32 × 32 Clipped rectified linear unit
Dropout0.5Not needed
Fully connected2 2 × 32 Not needed
Output data2Not neededSoftmax
Table 5. Vibration signal acquisition conditions and classification details from the CWRU dataset.
Table 5. Vibration signal acquisition conditions and classification details from the CWRU dataset.
Fault Diameter (mm)Motor Load
(HP)
Motor Speed
(Hz)
Sampling Frequency (kHz)Class
0.1778029.9548Fault
129.53
229.17
328.83
0.3556029.9548Fault
129.53
229.17
328.83
0.5334029.9548Fault
129.53
229.17
328.83
0.000029.9548Normal
129.53
229.17
328.83
Table 6. The optimized CNN model’s parameters.
Table 6. The optimized CNN model’s parameters.
ClassValue
Learning rate0.01
Batch size50
Max epoch40
Table 7. The accuracy and loss of the slope factor of the training data and test data.
Table 7. The accuracy and loss of the slope factor of the training data and test data.
Slope Factor 0.10.20.30.40.50.60.70.80.91.0Mean Std.
TrainingAccuracy84.38%82.81%84.38%93.75%99.38%95.31%90.63%82.81%71.11%67.36%85.19%9.671
Loss0.37740.41230.40220.12330.04550.12070.22370.47440.75820.79390.37310.244
TestAccuracy75.06%68.31%83.19%91.13%99.01%90.03%85.17%81.21%68.54%64.75%80.64%10.70
Loss0.49950.57360.39800.18650.06910.21320.26260.50810.75470.91120.06910.252
Table 8. The accuracy and loss of the image size of the training data and test data.
Table 8. The accuracy and loss of the image size of the training data and test data.
Image Size 5 ×   5 × 3 7 ×   7 × 3 11 ×   11 × 3 19 ×   19 × 3 35 ×   35 × 3 67 ×   67 × 3 131 ×   131 × 3 259 ×   259 × 3Mean Std.
TrainingAccuracy81.81%90.63%92.19%99.38%99.22%95.70%98.44%77.34%91.96%7.606
Loss0.39950.22370.19800.04550.04650.12020.09950.53640.20870.165
TestAccuracy75.06%83.58%93.64%99.01%99.14%95.31%97.59%75.06%89.80%9.688
Loss0.56010.26260.19460.06910.07090.13500.08330.56160.24220.194
Table 9. The accuracy and loss of the training data and test data of proposed model.
Table 9. The accuracy and loss of the training data and test data of proposed model.
Count 12345678910Mean Std.
TrainingAccuracy100%98.44%99.22%100%99.22%99.22%100%99.22%99.22%99.22%99.38%0.469
Loss0.01020.22680.03210.01240.04910.03180.01610.01610.02580.03490.04550.061
TestAccuracy99.74%95.96%99.40%99.91%98.19%98.62%99.74%99.94%99.40%99.05%99.01%1.158
Loss0.01820.12500.05790.05450.05770.04650.08630.08630.06840.09060.06910.028
Table 10. Confusion matrix of the experimental results using the RCE spectrogram.
Table 10. Confusion matrix of the experimental results using the RCE spectrogram.
Actual
FaultNormal
PredictedFault8730
Normal3287
Table 11. Confusion matrix of the experimental results using spectrogram.
Table 11. Confusion matrix of the experimental results using spectrogram.
Actual
FaultNormal
PredictedFault660286
Normal44264
Table 12. Confusion matrix of the experimental results using Mel spectrogram.
Table 12. Confusion matrix of the experimental results using Mel spectrogram.
Actual
FaultNormal
PredictedFault9460
Normal3080
Table 13. Confusion matrix of the experimental results using MFCC.
Table 13. Confusion matrix of the experimental results using MFCC.
Actual
FaultNormal
PredictedFault9460
Normal3080
Table 14. Performance metric for each input image.
Table 14. Performance metric for each input image.
Image RCE Spectrogram
(Proposed)
SpectrogramMel SpectrogramMFCC
Metric
Accuracy0.99740.73680.75440.7544
Recall0.99660.93750.75440.7544
Precision1.00000.69771.00001.0000
F1 score0.99830.80000.86000.8600
Table 15. Complexity comparison of CNNs.
Table 15. Complexity comparison of CNNs.
Classification TypesModel# of Parameters [KB]FLOPsAccuracy [%]
2D CNNOptimized CNN (Proposed)5.81 1.53 × 10 6 99.01
[20]206.32 1.01 × 10 7 99.89
[25]258.52 4.69 × 10 6 98.90
[27]367.81 2.00 × 10 8 98.77
[26]565.16 8.08 × 10 7 99.27
[21]1028.00 4.20 × 10 6 100
1D CNN[40]9.22 2.6 × 10 6 100
[17]50371.07 1.45 × 10 8 99.97
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

Seong, G.; Kim, D. An Intelligent Ball Bearing Fault Diagnosis System Using Enhanced Rotational Characteristics on Spectrogram. Sensors 2024, 24, 776. https://doi.org/10.3390/s24030776

AMA Style

Seong G, Kim D. An Intelligent Ball Bearing Fault Diagnosis System Using Enhanced Rotational Characteristics on Spectrogram. Sensors. 2024; 24(3):776. https://doi.org/10.3390/s24030776

Chicago/Turabian Style

Seong, Gyujin, and Dongwan Kim. 2024. "An Intelligent Ball Bearing Fault Diagnosis System Using Enhanced Rotational Characteristics on Spectrogram" Sensors 24, no. 3: 776. https://doi.org/10.3390/s24030776

APA Style

Seong, G., & Kim, D. (2024). An Intelligent Ball Bearing Fault Diagnosis System Using Enhanced Rotational Characteristics on Spectrogram. Sensors, 24(3), 776. https://doi.org/10.3390/s24030776

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