Next Article in Journal
Design of a Supraharmonic Monitoring System Based on an FPGA Device
Next Article in Special Issue
Analytical Model of Eccentric Induction Machines Using the Conformal Winding Tensor Approach
Previous Article in Journal
Effective Transfer Learning with Label-Based Discriminative Feature Learning
Previous Article in Special Issue
Detection and Analysis of Corrosion and Contact Resistance Faults of TiN and CrN Coatings on 410 Stainless Steel as Bipolar Plates in PEM Fuel Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intelligent Defect Diagnosis of Rolling Element Bearings under Variable Operating Conditions Using Convolutional Neural Network and Order Maps

by
Syed Muhammad Tayyab
,
Steven Chatterton
* and
Paolo Pennacchi
Department of Mechanical Engineering, Politecnico di Milano, Via G. La Masa 1, 20156 Milan, Italy
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(5), 2026; https://doi.org/10.3390/s22052026
Submission received: 12 January 2022 / Revised: 18 February 2022 / Accepted: 1 March 2022 / Published: 4 March 2022
(This article belongs to the Special Issue Feature Papers in Fault Diagnosis & Sensors Section 2022)

Abstract

:
Vibration analysis is an established method for fault detection and diagnosis of rolling element bearings. However, it is an expert oriented exercise. To relieve the experts, the use of Artificial Intelligence (AI) techniques such as deep neural networks, especially convolutional neural networks (CNN) have gained the attention of researchers because of their image classification and recognition capability. Most researchers convert the vibration signal into representative time frequency vibration images such as spectrograms and scalograms. These images are used as inputs to train the CNN model for fault diagnosis. Commonly, fault diagnosis is performed under same operating conditions, where models are trained and deployed for prediction under the same operating conditions. However, outside the laboratory environment, in real world applications, different operating conditions, such as variable speed, may be encountered. With the change in speed, the characteristic frequencies of the vibration signal will also change, which will result in changing the vibration image. Consequently, the performance of the CNN model may drop significantly for prediction under different operating conditions. Accessing the training data from all potential operating conditions may not be feasible for most real-world applications. Therefore, there is a need to find some signal properties which are invariant to change in operating conditions and only change due to change in health state so that models trained under one set of operating conditions may predict correctly under different operating conditions. This paper proposes a defect diagnosis method for rolling element bearings, under variable operating conditions (speed and load) based on CNN and order maps. These maps exhibit consistent properties under varying speed; therefore, they can be used to train the CNN model for fault diagnosis under variable speed. The effect of load change on these order maps is experimentally studied and it is found that the proposed method can undertake fault diagnosis on rolling element bearings under variable speeds and loads with high accuracy.

1. Introduction

Rolling element bearings are a vital part of rotating machines as they support the shaft, take load, and reduce the friction between different components. Their health status has a significant effect on the performance and availability of industrial machinery. As they are one of the most vulnerable components of rotating machines, bearing faults are a major cause of machine defects [1]. Bearing failure may cause serious damage to a machine and may result in the unavailability of important machinery, leading to financial loss and serious safety hazards. Therefore, the early defect diagnosis of rolling element bearings is crucial for the uninterrupted availability and operation of machines. To manage the potential failure of machinery, accurate fault detection and the diagnosis of bearings has gained significant researcher attention. Different methods have been introduced by researchers defectto diagnose the s of rolling element bearings, though vibration signal analysis in the time domain, frequency domain and time frequency domain have been most widely used [2,3]. However, it is an expert oriented task, and human involvement in general is not very effective or efficient in terms of responding quickly to large volumes of data. To minimize human dependency, intelligent condition monitoring techniques are attracting the attention of researchers. Over the past two decades, researchers have proposed many new approaches related to artificial intelligence techniques or conventional machine learning (ML) techniques [4], such as, the use of artificial neural networks (ANN) for rolling element bearing fault diagnosis [5,6], the principal omponent analysis (PCA) -based feature selection scheme for machine defect classification [7], a support vector machine (SVM) -based approach for machine condition monitoring and fault diagnosis [8], a bearing defect diagnosis approach, based on ANN and SVM [9,10], a hierarchical diagnosis network approach for fault pattern recognition in rolling element bearings [11], and the intelligent defect diagnostics of bearings and gears using k-nearest neighbors (KNN) as classifier along with the genetic algorithm (GA) [12].
Furthermore, Zhou et al. presented a study in whichthe K-mean algorithm was used to label the un-labelled signals [13]. A comparison of the performance between ANN and KNN was carried out by Gunerkar et al. [14]. The study was performed on rolling element bearings based upon a feature extracted after wavelet transform. J. P. Patel compared the performance of ANN with SVM for bearing fault diagnosis [15]. Tian et al. presented a bearing defect diagnosis method based on spectral kurtosis and cross correlation [16]. The method was validated by experiments using a machinery fault simulator. R. Zhang et al. undertook a fault diagnosis of a rolling element bearing using a k-mean singular value decomposition dictionary learning algorithm and a back propagation neural network [17]. Amar et al. used vibration spectrum imaging to train the ANN classifier for bearing fault diagnosis [18]. Khan and Kim trained a KNN classifier for bearing fault diagnosis under different speeds using features extracted from grayscale vibration images. The grey scale image textures were encrypted using the local binary pattern operator [19]. Ankush Mehta et al. used a combination of infrared thermography and different machine learning techniques for bearing fault diagnosis. They used two-dimensional discrete wavelet transform for feature extraction from thermal images. PCA was used for dimensionality reduction. KNN, linear discriminant analysis, and SVM were used as classifiers. A comparison of the performance of all three classifiers was also carried out [20]. Wu et al. proposed a bearing fault diagnosis technique through kernel matrix construction-based support vector machine [21].
The process of fault diagnosis through conventional machine learning techniques mainly consists of data acquisition, features extraction, features selection and finally fault classification or health state prediction. Classical machine learning algorithms are considered superficial because of their limited learning and generalization ability; therefore, they cannot extract important information from large, complex, and non-linear data. Although classical machine learning techniques have reduced human intervention in the fault diagnosis process, to apply classical ML techniques for bearing fault diagnosis, the correct features extraction and selection process is crucial, and this is still a labor-oriented exercise and requires expert knowledge to some extent. To further reduce human intervention in the bearing fault diagnosis process, techniques based on deep learning theories have been implemented. These deep learning-based techniques include autoencoder (AE) -based techniques, deep belief network (DBN) -based techniques, and convolutional neural network-based techniques. AE- and DBN-based techniques are unsupervised techniques and are mostly used for features extraction and information fusion. Models based on these techniques can be trained on unlabeled data and for labelled data fine tuning is required. As they are unsupervised techniques, AE-based techniques cannot be used directly for bearing health state identification, and a classification layer needs to be added at the top of the model architecture [22]. Whereas CNN-based techniques are supervised techniques and they are task specific. Due to their capability of joint features and classification learning, they can achieve better classification accuracy as compared to other deep learning models in image classification problems [23]. Deep learning models can automatically extract features from vibration data and manual features extraction is not required. Thus, deep learning techniques can provide end-to-end diagnosis models to further reduce human intervention in defect diagnosis activities. Deep learning models have gained the attention of researchers for bearing defect diagnosis, such as, the deep neural network (DNN) -based scheme for rolling element bearing fault diagnosis [24], the DNN for fault characteristic mining and the intelligent diagnosis of rotating machinery with massive data [25], the short-time Fourier transform-deep learning scheme for rolling bearing fault diagnosis [26], the bearing condition recognition method based on multi-feature extraction and DNN for the intelligent condition monitoring of bearings [27], the multi vibration signals and DBN scheme for bearing fault diagnosis [28], the CNN-based method for bearing and rotor fault detection [29], the hierarchical adaptive deep CNN approach for bearing fault diagnosis [30], the energy-fluctuated multiscale feature mining technique based on a wavelet packet energy image and a deep convolutional network for spindle bearing fault diagnosis [31], and the deep residual learning-based method for bearing fault diagnosis [32]. Moreover, Saucedo-Dorantes et al. [33] presented a fault diagnosis methodology for different bearing technologies (metallic, hybrid and ceramic bearings) based on deep feature learning. Vibration and current data were acquired under different operating conditions. Time domain, Frequency domain and Time Frequency domain features were extracted from vibration and current signals. A staked auto encoder-based feature learning method was introduced. Information from different domains was integrated through feature fusion and softmax layer was used for final classification. Ma et al. [34] presented a multi source information fusion algorithm based on variational autoencoder and random forest for bearing fault diagnosis in case of limited labeling. Hoang and Kang [35] proposed a bearing fault diagnosis method based upon motor current signal, using CNN and the information fusion technique. CNN is used for automatic features extraction from grey images from the current signal. After the application of the information fusion technique, classical machine learning algorithms are used for classification. Three working conditions are considered in this study. However, training and testing is performed under the same operating conditions.
Under the umbrella of deep learning, convolutional neural networks are a special type of DNNs which are known for their image recognition and classification capability. Researchers have used one dimensional (1D) and two dimensional (2D) CNN models for bearing defect diagnosis. In the 1D-CNN model, raw vibration data can be directly used as input to the model, as used by Eren et al. [36,37]. For 2D-CNN models, vibration data cannot be used in their raw form. Rather, they are initially converted into time-frequency image representations such as spectrograms [38,39,40], scalograms [40,41] or other types of vibration images [42]. Then these images which represent the vibration signal in image form are used as input in 2D-CNN model.
For conventional ML models and deep learning models, generally it is assumed that training and testing data sets belong to the same operating conditions. The performance of these models may degrade significantly if the operating conditions under which the models are deployed for prediction differ from the operating conditions under which the models were trained [43]. In real world applications, such as in trains, bearings may undergo different speed and loading conditions. Due to varying speed, the frequency characteristics of vibration signal also change. Therefore, accurate intelligent defect diagnosis of rolling element bearings under variable/inconsistent operating conditions is still a challenge for researchers. Moreover, enough training data may not be available for all potential operating conditions. This situation arises frequently which makes it challenging to use AI methods for fault diagnosis outside the laboratory environment. To solve this issue, transfer learning was introduced by researchers [44]. Zhang et al. proposed a transfer learning approach with neural networks for bearing fault diagnosis in changing working conditions when there is only a small amount of target data available [45]. Shen et al. presented a singular value decomposition-based features extraction and transfer learning-based approach for bearing fault diagnosis under various operating conditions [46]. Chao et al. presented an enhanced least squares support vector machine-based transfer learning approach for bearing fault diagnosis in the case of small target dataset availability [47]. Chunfeng et al. proposed a heterogeneous transfer learning-based approach for the scenarios when the available labelled data in the target domain is less [48]. Cheng et al. presented a deep adversarial transfer learning method based on Wasserstein distance for bearing fault diagnosis with insufficient labelled data [49]. However, most of the studies in transfer learning have normally considered some fixed operating conditions. Whereas, practically, the operating condition may lie between these fixed conditions. Furthermore, in most of the studies there is still a need for a small amount of labelled target domain data.
For bearing fault identification under variable speed, Bruand et al. reconstructed shaft orbits using angle measurements and derived a feature based on shaft orbit shape. It was shown that the signature of the rolling element bearing fault can be observed in shaft orbit. The information was retrieved using a set of angle measurement sensors. In order to remove speed variation, angular reference was used to work in the angular domain [50]. Order tracking-based methods have gained the attention of researchers for rolling element bearing diagnosis under variable conditions. Mishra et al. proposed a technique to use order tracking on the envelop of wavelet denoised estimate of the short duration angle synchronous averaged signal for bearing fault diagnosis under variable conditions. To extract deterministic content from the vibration signal, the Bayesian wavelet denoising approach was adapted and envelop order spectra were utilized to identify the faults [51]. Guo et al. proposed a rolling element bearing fault detection method under variable speed based upon envelop order tracking utilizing envelop analysis, order tracking and spectral kurtosis [52]. However, these methods cannot be used for automatic defect diagnosis as they are not intelligent methods and require human expertise in the field. For the automatic fault diagnosis of roller bearings under variable speed, Yang et al. extracted the features from a vibration signal by combining local mean decomposition and order tracking techniques. The features’ values were used as the input to a variable predictive model-based class discriminate classifier for automatic fault pattern and working condition identification [53]. Farhat et al. updated three frequency domain features, spectrum peak ratio outer, spectrum peak ratio inner, and spectrum peak ratio rolling element, to perform with a nonstationary signal, utilizing the order tracking technique. The updated features were used as an input to a multi-kernel support vector machine classifier for automatic defect classification [54]. However, these order tracking-based methods involve manual features extraction and a modification process.
To avoid the manual features extraction process, researchers have proposed end to end bearing diagnosis techniques based on CNN. Appana et al. proposed a low-speed bearing fault diagnosis method using CNN and envelope spectrums extracted from acoustic emission signal. The maximum speed considered was 500 revolutions per minute (rpm) [55]. Pham et al. presented a fault classification method utilizing the spectrogram images and CNN under four different speeds 1730 rpm, 1750 rpm, 1772 rpm and 1797 rpm [39]. In another study Pham et al. presented a fault classification method under different rotational speeds between 250 rpm and 500 rpm. Spectrogram images and CNN were utilized [38].
Researchers have mainly undertaken the defect diagnosis of bearings under variable operating conditions considering only one operating condition i.e., speed. However, bearings may experience variable loading conditions in some applications. Furthermore, speed variation is too low, for example between 1730 and 1797 rpm, or the proposed model is tested only on low speeds. Moreover, the considered speeds are steady; ramp or varying conditions are not considered. However, in real world applications, such as trains, continuously varyingied speed may be experienced which includes acceleration and deceleration. Furthermore, seeded defects are mostly researched and these may be an easier task for the machine learning or deep learning model as compared to diagnosesing the bearing defects which are encountered during real field operation. Therefore, there is a need to propose a robust method for bearing defect diagnosis under variable operating conditions, which can perform the task of fluctuations of speed and load. Moreover, to access the training data from all potential operating conditions may not be possible for most real-world applications. Therefore, it is necessary to find some signal properties which are invariant to changing operating conditions and only change due to changing health state so that AI models trained under one set of operating conditions may correctly predict, under different operating conditions, whose fault data for training is not available.
In this research a novel method for the fault diagnosis of rolling element bearings is proposed utilizing a combination of the order analysis and image classification and recognition capabilities of CNN. Order maps are computed from the vibration signal of rolling element bearings. These order maps show consistent properties over varying operating conditions. They vary with respect to the type of defect but show consistent configurations under changing operating conditions, therefore, they can be considered as operating condition (speed and load) invariant vibration images. The order maps computed under one set of operating condition are used to train the CNN model. The trained CNN model is deployed for defect diagnosis under different operating conditions (speed and load). Performance of the proposed method is compared with five techniques: (1) a technique based on spectrograms and CNN as utilized in [38,39,40], (2) KNN in combination with GA as utilized by [12], (3) ANN/multilayer perceptron, (4) SVM, (and5) KNN + +order maps.

2. Proposed Defect Diagnosis Methodology

The intelligent fault diagnosis of rolling element bearings using a vibration signal requires features extraction if we use conventional machine learning techniques. For accurate diagnosis, correct features extraction and selection is of paramount importance [12]. However, deep learning models such as convolutional neural networks (2D-CNN) have the capability to learn discriminative features directly from vibration characteristic images such as spectrograms, scalograms etc. Therefore, tedious features extraction and selection can be avoided by using CNN models for fault diagnosis. However, due to a change in operating conditions, the vibration images of the signal may also change, resulting in degradation of the performance of the CNN model for fault diagnosis under different operating conditions. Therefore, in this study a fault diagnosis method for rolling element bearings is proposed based upon the operating conditions of invariant vibration images and a deep convolutional neural network. Variation in two operating conditions: speed and load, are considered in this study.
Vibration signals of locomotive rolling element bearings are acquired under different speeds (steady and variable) and loads. Large fluctuations of speed including transient conditions are considered. Tachometer signal is also simultaneously acquired. Different types of real field and seeded defects of locomotive rolling element bearings, including combined defects, are considered in this study. Vibration signals are segmented into 2 segments with an overlap of 1.6 s. Order maps are computed from each segment. These order maps are used as input for CNN model which automatically learns the features from these images/order maps to classify the different defect types of rolling element bearings.
Initially, the CNN model is trained and tested under the same operating conditions. Data are divided into the training, validation, and testing data sets in the ratio of 70%, 15% and 15%, respectively. After completing the training and testing process under the same operating conditions, the trained model is deployed for defect diagnosis under all of the other operating conditions considered in this study. The same process is repeated for the other combinations of speed and load as well. The proposed methodology is depicted in Figure 1. Details of order maps computation, the CNN model and its implementation are given in following subsections:

2.1. Order Maps

For rolling element bearings’ fault diagnosis, characteristic frequencies play an important role for identifying a specific type of defect. As frequency indicates the repetition times per second, in the case of varying shaft speed, spectral line smearing may cause difficulty in fault diagnosis using spectrograms [56]. Consequently, traditional Fourier Transform-based methods are no longer effective for bearing fault diagnosis under varying rotational speed. Thus, the time-frequency vibration images generated from the vibration signal, based on short-time Fourier Transform such as spectrograms will also undergo change if the shaft speed is changed. Consequently, a CNN model trained on such vibration images for fault diagnosis may not be able to predict accurately under varying speed.
The problem can be solved by converting a non-stationary signal into a stationary signal in the angular domain by resampling at constant angular intervals. This technique is called order tracking. An order indicates a frequency which is a specific multiple of the rotational speed. Order tracking has been identified as a trustworthy and practical approach to mitigate the effects of spectral line smearing caused by varying rotational speed [57,58]. In this study, order maps are computed from the vibration signals utilizing the tachometer pulse according to the procedure shown in Figure 2. These order maps show consistent patterns under varying speed and are used as input to the CNN model.
Order maps computation mainly involves the following three sections:

2.1.1. Tachometer Signal Processing and rpm Extraction

Tachometer signal processing is very important in order tracking. If a good stable tachometer signal is not available, then accurate order maps cannot be computed [59]. Tachometer pulse is converted into an rpm signal as given by Equation (1) [60].
r p m ( t ) = 60 N p ( t 2 t 1 )
where, N p is number of pulses per revolution and ( t 2 t 1 ) is the time instance between two pulses. For this purpose, initially low and high states of tachometer signal are determined. Time for each pulse is determined by averaging the start and end time readings of pulse. Time interval between pulse centres is determined to find the rpm at the interval midpoint by rpm = 60/Δt. Afterwards, instantaneous rpm values are interpolated linearly onto the time axis of original signal.

2.1.2. Synchronous Resampling in the Order Domain

Phase angle ( A ( t ) ) is determined as the time integral of rotational speed as shown in Equation (2).
A ( t ) = 0 t r p m ( t ) 60 d t
The vibration signal is resampled onto the new time axis instead of the original constant time axis. Resampling of the vibration signal is done at non sampled time points by utilizing an upsampled vibration signal and interpolating it linearly. Mostly, 10–20 times oversampling gives adequately good accuracy [60]. In this study signal is upsampled by a factor of 15 and then linearly interpolated onto a consistent grid in the phase domain. Thus, a constant number of samples per cycle is achieved to produce a stationary sinusoid for each order. After synchronous resampling, the vibration signal comes in the order domain. The signal frequencies which are constant multiples of the rotational speed are now converted into constant tones. The smearing of spectral components which occurs due to rapid frequency change with time, reduces because of this transformation. The relationship between the highest accessible order ( O m a x ), sampling frequency ( f s ) and the maximum rotational speed of the shaft is given by Equation (3).
O m a x = f s 2 ( max ( r p m ) 60 )
To accurately capture the maximum order, the angular domain sampling rate ( f a ) must be at least two times that of O m a x . i.e., f a 2 O m a x .

2.1.3. Short-Time Fourier Transform of Resampled Signal in the Order Domain

Short-time Fourier transform of the interpolated signal is computed to generate a spectral map of order versus rpm. Since each order is a fixed multiple of the reference rotational speed, the order map has a straight track as a function of rpm for each order, as shown in Figure 3. Therefore, these maps show consistent patterns under varying speed. The maximum order which can be accurately captured is dependent on sampling frequency and maximum rotational speed, as shown by Equation (3). Therefore, to use the order map images as input to a CNN, it is essential to keep the same maximum order limit for the whole diagnosis process lower than the O m a x related to maximum speed under consideration and the sampling frequency.
To ensure the normal operation of rolling element bearings, sufficient friction force is required by the race ways. If the bearing is not properly loaded, slip page may appear in the contact zone and relative skid may occur between the rolling elements and the inner/outer race ways [61]. This phenomenon generally occurs in low loaded roller bearings [62]. Because of this phenomenon, frequencies may be slightly affected [63], which in return may affect the order maps. Consequently, performance of the deep CNN model for fault diagnosis under different loading conditions may degrade. Han et al. studied the skidding behavior of cylindrical roller bearings under variable loads. It was reported that an increase in radial load and bending moment reduces the skidding. By increasing the load, the friction force of the rollers increases, and the maximum slipping velocity of the rollers reduces, leading to the pure rolling rotation of the cage. Therefore, after certain radial load, no overall skidding phenomenon was present in the bearing [64]. Feng et al. reported that an increase in radial load can reduce the skidding of ball bearings [65]. Deep groove ball bearings are less sensitive to skidding because of an extra degree of freedom of the rolling elements. Therefore, if the bearing is properly loaded and diagnosis is not performed under a load or very low loading conditions, the order maps may not be significantly affected, and the CNN model will perform defect diagnosis correctly under the changing loading conditions as well. The effect of load change on the order maps for deep groove ball bearings is observed experimentally. Order maps of a healthy bearing, a bearing with combined defects and a bearing with outer race defect are shown in Figure 4 for three different speeds (1000 rpm, 2000 rpm and 3000 rpm) and two different loads (5 KN and 15 KN) which are the minimum and maximum loads considered in this study. It can be observed that the patterns of order maps do not change with changing speed and load. However, they show different patterns for different fault conditions, making themselves insensitive to operating conditions (speed and load) but fault discriminative. Therefore, they will be used for fault diagnosis under variable speeds and loads.

2.2. Convolutional Neural Network

Convolutional neural networks are known for their high capability in the field of image identification and classification. They consist of an input layer, hidden layers, and an output layer. Generally hidden layers in convolutional neural networks consist of convolutional layers, Rectified Linear Linear Unit (ReLU) layers, pooling or subsampling layers and fully connected layers [22]. The 2D-convolutional layer uses its kernels to convolve the input by moving the kernels vertically and horizontally and getting the dot product of kernels and input and then addings a bias term. Its input is the output of the previous layer. Kernels extract the local features of the input region. An activation function, such as ReLU, is used to obtain output from the results of the convolutional operation which are called features. Sometimes ReLU is shown as a separate layer from the convolutional layer. The ReLU layer performs a threshold operation in which all the values less than zero are set to zero. The mathematical model of the convolutional layer is given in Equation (4).
X j m = f ( i C j X i m 1 k i j m + b j m )
C j represents the input map selection, m is the mth layer in the network. X i m 1 is the input of the convolutional channel. k is the kernel matrix, and b is the bias matrix. f is a nonlinear activation function such as ReLU. After the convolutional layer, a pooling or subsampling layer is applied which reduces the size of the input features and network parameters. The pooling layer can be represented by Equation (5).
X j m = f ( β j m   d o w n ( X j m 1 + b j m ) )
where, d o w n (.) represents the pooling function. b j m is additive bias and β j m is multiplicative bias. Based upon the pooling function, the pooling operation can be maximum pooling or average pooling. After staking the convolutional and pooling layers multiple times, the output is fed to a fully connected layer at the final stage. It is a feed forward neural network (multi-layer perceptron) which uses softmax as an activation function in the output. The softmax activation function can be described by Equation (6).
σ ( y ) i = e y i j = 1 N e y j
where, σ is the softmax, y is the input, e y i is the exponential function for the input, e y j is the exponential function for the output and N is the number of classes. The purpose of the fully connected layer is to collect all the features learned from the previous layers to identify patterns. Therefore, all the neurons in the fully connected layer are connected to all neurons in the previous layer. Moreover, in order to increase the training speed of the convolutional neural network and to reduce the sensitivity to network initialization, a batch normalization layer is used between the convolutional layer and nonlinearities. This layer independently normalizes the minibatch data across all channels. Furthermore, to avoid the possibility of the network memorizeing some specific features, a dropout layer can be added prior to the fully connected layer which sets the input to zero as per set probability.
The convolutional neural network architecture proposed in this study is shown in Figure 5. The order maps were saved in Bitmap Image file format in original size (656 × 875 × 3) without compression to avoid loss of information. Batch normalization layers are added between 2D-convolutional layers and ReLU layers. A combination of the 2D convolutional layer, the batch normalization layer and the ReLU layer makes one convolutional block. Overlapping max-pooling layers are added after 1st, 2nd, 3rd, the first second third and sixth convolutional blocks. A dropout layer with 50% drop out probability is added before the fully connected layer in order to avoid over fitting. Adaptive moment estimation optimizer is used for optimization of hyperparameters because of lesser memory and tuning requirement and faster optimization capability as compared to other optimization algorithms [40]. The initial learn rate was set as 0.0003 and it was dropped by a factor of 0.1 after every 10 epochs. The minibatch size was set as 15.

3. Experimental Setup

The proposed method for defect diagnosis of rolling element bearings under different operating conditions was validated using the vibration data from locomotive rolling element bearings. A total of two case studies were conducted utilizing two different test rigs developed to test bearings for railway electrical traction. Details of the same are appended below.

3.1. Case Study-1

The first test rig (shown in Figure 6) was designed to test locomotive motor bearings under different operating conditions (speed and load). Operating speed and load can be varied in order to get data under different operating conditions. Vibration data were acquired at 25,600 Hz sampling frequency by two accelerometers installed at different angular positions in radial direction. A tachometer was installed at the bearing shaft to give one pulse per revolution signal. For each operating condition, data were acquired for 50 s. SKF 6318, deep groove ball bearings were tested under nine different operating conditions which are described in Table 1. A total of three health states including one normal and two faulty conditions with real field defects are considered in this study, as shown in Figure 7 and Table 2. Raw vibration signals for all three types of bearings considered in case study-1, at 1000 rpm and 5 KN load, are shown in Figure 8.

3.2. Case Study-2

An actual traction system of a high-speed train with an ability to reach 250 km/h speed was installed on the test bench, as shown in Figure 9 [66]. The traction motor is 265 kW converter driven 4-poles asynchronous motor. The rotor of the traction motor is mounted on two types of bearings: a single row, groove ball bearing (BB), SKF-6214, and a cylindrical roller bearing (RB), SKF-N214, on the driven end, as shown in Figure 9. Vibration data at 20,000 Hz sampling frequency are acquired for both support bearings of the traction motor at different speeds (constant speed, acceleration, and deceleration), as presented in Table 3. Accelerometers were installed in the proximity of the bearings, on the inner surface of motor flanges holding the respective bearings, as shown in Figure 9. Tachometer signal was also acquired simultaneously. For both types of the bearings, one healthy and four defective bearings with seeded defects of different types were considered in this study, as shown in Figure 10 and Figure 11 and Table 4.

4. Results and Discussion

During the first case study on test rig 1, vibration data were acquired using two accelerometers under nine different operating conditions (different combinations of speed and load). One healthy and two defective deep groove ball bearings with outer race defects and combined defects were considered. Vibration signals were divided into segments of 2 s with an overlap of 1.6 s. Order maps were computed from the vibration signals of both sensors separately and were combined to train the proposed CNN model for fault diagnosis under one set of operating conditions. The trained model was tested for defect diagnosis under all of the other combinations of speed and load considered during the first case study. The same procedure was repeated for all sets of operating conditions considered in this study. The performance of the proposed methodology for defect diagnosis under different operating conditions in terms of prediction accuracy is given in Table 5, and the same is compared with five methods: (1) the methodology used in [38,39,40], in which time frequency images and spectrograms were used as input to the CNN, (2) KNN as a classifier in combination with GA for features selection using the time domain and spectral kurtosis-based features as utilized in [12], (3) ANN/multilayer perceptron using the same features as utilized in [12], (4) SVM using the same features as utilized in [12], and (5) KNN + order maps using Histogram of Oriented Gradient (HOG) features. Initially, the models were trained at 1000 rpm and 5 KN load. The trained models were deployed for predictions under all other combinations of speed and load. The proposed CNN model predicted with 100% accuracy under all loading conditions at 1000 rpm for both cases i.e., with order maps and spectrograms. Afterwards, for spectrograms, performance degraded drastically when the model trained at 1000 rpm and 5 KN load was deployed for predictions at 2000 rpm (5, 10 and 15 KN load) and 3000 (5, 10 and 15 KN load). However, performance of the proposed methodology in terms of fault diagnosis/prediction accuracy remained excellent under all the loading conditions at 2000 rpm and 3000 rpm with a minimum fault diagnosis accuracy of 96% at 3000 rpm and 10 KN load. However, the minimum fault diagnosis accuracy was 46% when spectrograms were used as input for the CNN model. The fault diagnosis accuracy ofthe other four methods considered for comparison i.e., KNN + GA, ANN, SVM and KNN + order maps was 98.3%, 99.2%, 99.6% and 98.3%, respectively, under thesame operating conditions (when the models were trained and tested at 1000 rpm and 5 KN load). However, when these trained models were deployed for diagnosis under different operating conditions as compared to those under which they were trained, their performance degraded drastically except for KNN + order maps. By changing the load under same speed, the performance of KNN + GA, ANN and SVM was degraded to 87.9%, 89.1% and 88.2%, respectively. However, performance degradation was more drastic by changing the speed which dropped to 46.9%, 35.4% and 58.3% for KNN, ANN and SVM, respectively. For KNN + order maps the performance degradation was not significant by changing the load, however, fault diagnosis accuracy dropped to a minimum of 89.7% by changing the speed which is much better than KNN + GA, ANN, SVM and CNN + spectrograms. However, the proposed method performed better as compared to KNN + order maps. A performance comparison of all models when trained at 1000 rpm and 5 KN load is given in Figure 12.
The same trend was observed when the procedure of training under one set of operating conditions and testing under all other combinations of operating conditions was repeated for all operating conditions considered in this study. It is observed that change in the loading condition at steady speed did not affect the performance of the CNN model in both cases i.e., with order maps and with spectrograms. Similarly, the performance of the KNN model with order maps was also not affected greatly by changing the load. This depicts that spectrograms and order maps are not affected much by changing the load. Although, the performance of KNN + GA, ANN and SVM models was degraded at an average of about 13%. When the speed under which the models were trained differed from the speed under which the models were deployed to make predictions, the performance of the CNN model was adversely affected in the case of spectrograms, which shows that by changing the operating speed spectrograms undergo change and thus the performance of the CNN model for fault diagnosis is affected. Similarly, drastic performance degradation was observed for KNN + GA, ANN and SVM models. The performance of KNN + order maps also degraded but not drastically. However, the proposed methodology showed very good performance under changing speed as well, because of the invariant behaviour of order maps under varying speeds the good image classification and recognition ability of the deep CNN model. The proposed method did diagnosis the faults of a locomotive rolling element bearing under different operating conditions with an overall average accuracy of 98.4% which is much higher compared to the average accuracy of 73.7%, 67.4%, 64.7% and 65.3% for CNN + spectrograms, KNN + GA, ANN and SVM, respectively. For KNN + order maps the overall average prediction accuracy was 92.3% which is much higher compared to other methods used for comparison, but it is less than the average fault detection accuracy of the proposed methodology. The explanation of the better performance of CNN and KNN in the case of order maps is that the order maps exhibit consistent patterns under variable speed, as described in Section 2.1, and for deep groove ball bearings under the load range considered in this study they did not exhibit much change. Therefore, the CNN and KNN models were able to correctly classify them. However, because of the better image classification abilities of the proposed deep CNN model, the proposed method performed better as compared to KNN + order maps. In order to use the KNN model for this fault diagnosis task we had to extract HOG features from order maps, whereas the deep CNN model can automatically learn discriminative features from the order maps for defect diagnosis under variable speeds and loads with high accuracy.
In the second case study, vibration data were acquired for two types of rolling element bearings installed in the locomotive traction motor, under five different variable speed conditions (constant speed, acceleration, and deceleration). The load was constant which corresponds to the motor shaft weight. Order maps were computed from vibration signal segments of 2 s to use as input for the proposed CNN model. An overlap of 1.6 sec was used for each segment. Initially, the CNN model was trained at a steady speed of 3170 rpm and was tested under all other operating conditions considered during this study. Later, the model was trained under the transient condition of variable speed between 1000 and 2000 rpm, and was deployed for defect diagnosis under all other operating conditions considered in this study. The performance of the proposed methodology in terms of prediction accuracy is given in Table 6. The proposed defect diagnosis method undertook a fault diagnosis of ball bearings and cylindrical roller element bearings installed in the locomotive traction motor under variable speeds with an average accuracy of 99.2% and 98%, respectively. In this case study, the overall accuracy of the proposed methodology for defect diagnosis of rolling element bearings under variable speed remained 98.6%, which proves that the proposed method is capable for the defect diagnosis of rolling element bearings with good accuracy at a high speed with large variations. A confusion matrix for ball bearings, when the model was trained at 3170 rpm and tested at 4955 rpm is shown in Figure 13, and when the model was trained at a variable speed between 1000 and 2000 rpm and tested at variable speed between 3525 and 4125 rpm is shown in Figure 14.
For the implementation of the proposed methodology, in addition to the vibration signal, a tachometer signal is also required. However, most of the other methods, including thefour methods used for comparison in first case study i.e., CNN + spectrograms, KNN + GA, ANN and SVM, do not require a tachometer signal. Therefore, the necessity of a very accurate tachometer signal is the limitation of the proposed methodology. If there is an error in the tachometer, then resampling may fail and correct order maps cannot be computed, which will adversely affect the diagnosis performance of the proposed methodology.

5. Conclusions

In this study a method for the defect diagnosis of rolling element bearings under variable operating conditions (speed and load) using order maps and convolutional neural networks is proposed. A non stationery signal was resampled synchronously at con-stant angular intervals to convert it into a stationery signal in the order domain. Short-time Fourier Transform of the resampled signal was computed to generate the order maps. These maps show consistent patterns under variable speed but different patterns for different types of defects. The sensitivity of order maps for changing load was studied experimentally for deep groove ball bearings and it was found that they remain consistent under varying loads if the bearings are properly loaded. Therefore, order maps can be termed as operating condition (speed and load) invariant but fault discriminative vibration images. In addition, due to this property, they can be used for fault diagnosis under varying speeds and loads. A deep CNN model was proposed which can automatically extract fault discriminating features from the order maps for defect classification. Order maps were used as input to the CNN model for fault diagnosis under varying speeds and loads. The proposed method conducted the fault diagnosis of different types of locomotive rolling element bearings under a huge fluctuation of operating conditions (speed and load) with an average prediction accuracy of 98.4% and 98.6% in two separate case studies. The proposed method outperformed when it was compared with other CNN, KNN, ANN and SVM-based methods. The limitation of the proposed method is the requirement of a very accurate tachometer signal. The proposed method can be implemented for rolling element bearings’ fault diagnosis under variable speeds and loads, such as in the transportation industry. In future studies, the efficacy of the proposed methodology needs to be investigated for such applications where a tachometer signal is not available along with the vibration signal, by extracting the rpm information from the vibration signal using signal processing techniques.

Author Contributions

Conceptualization, P.P. and S.M.T.; methodology, S.M.T.; software, S.M.T. and S.C.; validation, S.M.T. and S.C.; formal analysis, S.M.T.; investigation, S.M.T.; resources, P.P. and S.C.; data curation, S.M.T. and S.C.; writing—original draft preparation, S.M.T.; writing—review and editing, P.P. and S.C.; visualization, S.M.T.; supervision, P.P.; project administration, S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Luo, J.Q.; Zhu, L.C.; Li, Q.F.; Liu, D.P.; Chen, M.Y. Imbalanced Fault Diagnosis of Rotating Machinery Based on Deep Generative Adversarial Networks with Gradient Penalty. Processes 2021, 9, 1751. [Google Scholar] [CrossRef]
  2. Kiral, Z.; Hira, K. Simulation and analysis of vibration signals generated by rolling element bearing with defects. Tribol. Int. 2003, 36, 667–678. [Google Scholar] [CrossRef]
  3. Orhan, S.; Nizami, A.; Celik, V. Vibration monitoring for defect diagnosis of rolling element bearings as a predictive maintenance tool: Comprehensive case studies. Ndt E Int. 2006, 39, 293–298. [Google Scholar] [CrossRef]
  4. Chang, C.W.; Lee, H.W.; Liu, C.H. A review of artificial intelligence algorithms used for smart machine tools. Inventions 2018, 3, 41. [Google Scholar] [CrossRef] [Green Version]
  5. Li, B.; Chow, M.Y.; Yodyium, T.; James, C.H. Neural-network-based motor rolling bearing fault diagnosis. IEEE Trans. Ind. Electron. 2000, 47, 1060–1069. [Google Scholar] [CrossRef] [Green Version]
  6. Samanta, B.; Al-Balushi, K.R. Artificial neural network based fault diagnostics of rolling element bearings using time-domain features. Mech. Syst. Signal Processing 2003, 17, 317–328. [Google Scholar] [CrossRef]
  7. Malhi, A.; Gao, R.X. PCA-based feature selection scheme for machine defect classification. IEEE Trans. Instrum. Meas. 2004, 53, 1517–1525. [Google Scholar] [CrossRef]
  8. Su, H.; Chong, K.T. Induction machine condition monitoring using neural network modeling. IEEE Trans. Ind. Electron. 2007, 54, 241–249. [Google Scholar] [CrossRef]
  9. Seryasat, O.R.; Haddadnia, J.; Arabnia, Y.; Zeinali, M.; Abooalizadeh, Z.; Taherkhani, A.; Tabrizy, S.; Maleki, F. Intelligent fault detection of ball-bearings using artificial neural networks and support-vector machine. Life Sci. 2012, 9, 4186–4189. [Google Scholar]
  10. Mao, W.; He, J.; Yuan, L.; Yan, Y.G. Bearing fault diagnosis with auto-encoder extreme learning machine: A comparative study. Proc. Inst. Mech. Eng. C J. Mech. Eng. Sci. 2017, 231, 1560–1578. [Google Scholar] [CrossRef]
  11. Gan, M.; Cong, W. Construction of hierarchical diagnosis network based on deep learning and its application in the fault pattern recognition of rolling element bearings. Mech. Syst. Signal Processing 2016, 72, 92–104. [Google Scholar] [CrossRef]
  12. Tayyab, S.; Muhammad, E.; Asghar, P.P.; Steven, C. Intelligent fault diagnosis of rotating machine elements using machine learning through optimal features extraction and selection. Procedia Manuf. 2020, 72, 266–273. [Google Scholar] [CrossRef]
  13. Zhou, Q.C.; Shen, H.H.; Zhao, J.; Liu, X.C.; Xiong, X.L. Degradation state recognition of rolling bearing based on K-means and CNN algorithm. Shock. Vib. 2019, 2019, 8471732. [Google Scholar] [CrossRef]
  14. Gunerkar, R.S.; Jalan, A.K.; Belgamwar, S.U. Fault diagnosis of rolling element bearing based on artificial neural network. J. Mech. Sci. Technol. 2019, 33, 505–511. [Google Scholar] [CrossRef]
  15. Patel, J.P.; Upadhyay, S.H. Comparison between artificial neural network and support vector method for a fault diagnostics in rolling element bearings. Procedia Eng. 2016, 144, 390–397. [Google Scholar] [CrossRef] [Green Version]
  16. Jing, T.; Morillo, C.; Michael, H.A.; Pecht, M. Motor bearing fault detection using spectral kurtosis-based feature extraction coupled with K-nearest neighbor distance analysis. IEEE Trans. Ind. Electron. 2015, 63, 1793–1803. [Google Scholar]
  17. Zhang, R.X.; Fang, Y.; Zhou, Z.F. Fault diagnosis of rolling bearing based on k-svd dictionary learning algorithm and BP Neural Network. In Proceedings of the 2019 International Conference on Robotics, Intelligent Control and Artificial Intelligence, Shanghai, China, 20 September 2019; pp. 341–345. [Google Scholar] [CrossRef]
  18. Muhammad, M.; Gondal, I.; Wilson, C. Vibration spectrum imaging: A novel bearing fault classification approach. IEEE Trans. Ind. Electron. 2014, 62, 494–502. [Google Scholar] [CrossRef]
  19. Khan, S.A.; Kim, J.M. Automated bearing fault diagnosis using 2D analysis of vibration acceleration signals under variable speed conditions. Shock. Vib. 2016, 2016, 8729572. [Google Scholar] [CrossRef] [Green Version]
  20. Ankush, M.; Goyal, D.; Choudhary, A.; Pabla, B.S.; Belghith, S. Machine Learning-Based Fault Diagnosis of Self-Aligning Bearings for Rotating Machinery Using Infrared Thermography. Math. Probl. Eng. 2021, 2021, 9947300. [Google Scholar] [CrossRef]
  21. Wu, C.X.; Chen, T.F.; Jiang, R. Bearing fault diagnosis via kernel matrix construction based support vector machine. J. Vibroengineering 2017, 19, 3445–3461. [Google Scholar] [CrossRef] [Green Version]
  22. Lei, Y.G.; Yang, B.; Jiang, X.W.; Feng, J.; Li, N.P.; Asoke, K.N. Applications of machine learning to machine fault diagnosis: A review and roadmap. Mech. Syst. Signal Processing 2019, 138, 106587. [Google Scholar] [CrossRef]
  23. Gu, J.X.; Wang, Z.H.; Kuen, J.; Ma, L.Y.; Shahroudy, A.; Bing, S.; Liu, T.; Wang, X.; Wang, G.; Cai, J.; et al. Recent advances in convolutional neural networks. Pattern Recognit. 2018, 77, 354–377. [Google Scholar] [CrossRef] [Green Version]
  24. Lu, W.N.; Wang, X.Q.; Yang, C.C.; Zhang, T. A novel feature extraction method using deep neural network for rolling bearing fault diagnosis. In Proceedings of the 27th Chinese Control and Decision Conference, Qingdao, China, 23–25 May 2015; pp. 2427–2431. [Google Scholar]
  25. Jia, F.; Lei, Y.; Lin, J.; Zhou, X.; Lu, N. Deep neural networks: A promising tool for fault characteristic mining and intelligent diagnosis of rotating machinery with massive data. Mech. Syst. Signal Processing 2016, 72, 303–315. [Google Scholar] [CrossRef]
  26. Liu, H.M.; Li, L.F.; Ma, J. Rolling bearing fault diagnosis based on STFT-deep learning and sound signals. Shock. Vib. 2016, 2016, 6127479. [Google Scholar] [CrossRef] [Green Version]
  27. Guo, L.; Gao, H.; Huang, H.; He, X.; Li, S. Multifeatures fusion and nonlinear dimension reduction for intelligent bearing condition monitoring. Shock. Vib. 2016, 2016, 4632562. [Google Scholar] [CrossRef] [Green Version]
  28. Tao, J.; Liu, Y.; Yang, D. Bearing fault diagnosis based on deep belief network and multisensor information fusion. Shock. Vib. 2016, 2016, 9306205. [Google Scholar] [CrossRef] [Green Version]
  29. Janssens, O.; Slavkovikj, V.; Vervisch, B.; Stockman, K.; Loccufier, M.; Verstockt, S.; Van de Walle, R.; Van Hoecke, S. Convolutional neural network based fault detection for rotating machinery. J. Sound Vib. 2016, 377, 331–345. [Google Scholar] [CrossRef]
  30. Guo, X.; Chen, L.; Shen, C. Hierarchical adaptive deep convolution neural network and its application to bearing fault diagnosis. Measurement 2016, 93, 490–502. [Google Scholar] [CrossRef]
  31. Ding, X.; He, Q. Energy-fluctuated multiscale feature learning with deep convnet for intelligent spindle bearing fault diagnosis. IEEE Trans. Instrum. Meas. 2017, 66, 1926–1935. [Google Scholar] [CrossRef]
  32. Zhang, W.; Li, X.; Ding, Q. Deep residual learning-based fault diagnosis method for rotating machinery. ISA Trans. 2019, 95, 295–305. [Google Scholar] [CrossRef]
  33. Saucedo-Dorantes, J.J.; Arellano-Espitia, F.; Delgado-Prieto, M.; Osornio-Rios, R.A. Diagnosis Methodology Based on Deep Feature Learning for Fault Identification in Metallic, Hybrid and Ceramic Bearings. Sensors 2021, 21, 5832. [Google Scholar] [CrossRef] [PubMed]
  34. Ma, J.; Li, C.; Zhang, G. Rolling Bearing Fault Diagnosis Based on Deep Learning and Autoencoder Information Fusion. Symmetry 2022, 14, 13. [Google Scholar] [CrossRef]
  35. Hoang, D.T.; Kang, H.J. A motor current signal-based bearing fault diagnosis using deep learning and information fusion. IEEE Trans. Instrum. Meas. 2019, 69, 3325–3333. [Google Scholar] [CrossRef]
  36. Eren, L. Bearing fault detection by one-dimensional convolutional neural networks. Math. Probl. Eng. 2017, 2017, 8617315. [Google Scholar] [CrossRef] [Green Version]
  37. Eren, L.; Ince, T.; Kiranyaz, S. A generic intelligent bearing fault diagnosis system using compact adaptive 1D CNN classifier. J. Signal Processing Syst. 2019, 91, 179–189. [Google Scholar] [CrossRef]
  38. Pham, M.-T.; Kim, J.-M.; Kim, C.-H. 2D CNN-Based Multi-Output Diagnosis for Compound Bearing Faults under Variable Rotational Speeds. Machines 2021, 9, 199. [Google Scholar] [CrossRef]
  39. 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]
  40. Verstraete, D.; Ferrada, A.; Droguett, E.L.; Meruane, V.; Modarres, M. Deep learning enabled fault diagnosis using time-frequency image analysis of rolling element bearings. Shock. Vib. 2017, 2017, 5067651. [Google Scholar] [CrossRef]
  41. Hoang, D.T.; Tran, X.T.; Van, M.; Kang, H.J. A Deep Neural Network-Based Feature Fusion for Bearing Fault Diagnosis. Sensors 2021, 21, 244. [Google Scholar] [CrossRef]
  42. Hoang, D.-T.; Kang, H.-J. Rolling element bearing fault diagnosis using convolutional neural network and vibration image. Cogn. Syst. Res. 2019, 53, 42–50. [Google Scholar] [CrossRef]
  43. Tayyab, S.M.; Chatterton, S.; Pennacchi, P. Fault detection and severity level identification of spiral bevel gears under different operating conditions using artificial intelligence techniques. Machines 2021, 9, 173. [Google Scholar] [CrossRef]
  44. Wei, D.; Wang, K.; Heyns, S.; Zuo, M.J. Convolutional neural networks for fault diagnosis using rotating speed normalized vibration. In Proceedings of the International Conference on Condition Monitoring of Machinery in Non-Stationary Operation, Santander, Spain, 20–22 June 2018; Springer: Cham, Switzerland; Berlin/Heidelberg, Germany, 2018. [Google Scholar] [CrossRef]
  45. Zhang, R.; Tao, H.; Wu, L.; Guan, Y. Transfer learning with neural networks for bearing fault diagnosis in changing working conditions. IEEE Access 2017, 5, 14347–14357. [Google Scholar] [CrossRef]
  46. Shen, F.; Chen, C.; Yan, R.; Gao, R.X. Bearing fault diagnosis based on SVD feature extraction and transfer learning classification. In Proceedings of the 2015 Prognostics and System Health Management Conference (PHM), Beijing, China, 21–23 October 2015; pp. 1–6. [Google Scholar] [CrossRef]
  47. Chen, C.; Shen, F.; Yan, R.Q. Enhanced least squares support vector machine-based transfer learning strategy for bearing fault diagnosis. Chin. J. Sci. Instrum. 2017, 38, 33–40. [Google Scholar]
  48. Chunfeng, W.; Zheng, L.; Jun, Z.; Wei, W. Heterogeneous transfer learning based on stack sparse auto-encoders for fault diagnosis. In Proceedings of the 2018 Chinese Automation Congress (CAC), Xi’an, China, 30 November–2 December 2018; pp. 4277–4281. [Google Scholar] [CrossRef]
  49. Cheng, C.; Zhou, B.; Ma, G.; Wu, D.; Yuan, Y. Wasserstein distance based deep adversarial transfer learning for intelligent fault diagnosis with unlabeled or insufficient labeled data. Neurocomputing 2020, 409, 35–45. [Google Scholar] [CrossRef]
  50. Bruand, G.; Chatelain, F.; Granjon, P.; Martin, N.; Duret, C. Reconstructing shaft orbit using angle measurement to detect bearing faults. Mech. Syst. Signal Processing 2020, 139, 106561. [Google Scholar] [CrossRef] [Green Version]
  51. Mishra, C.; Samantaray, A.K.; Chakraborty, G. Rolling element bearing defect diagnosis under variable speed operation through angle synchronous averaging of wavelet de-noised estimate. Mech. Syst. Signal Processing 2016, 72, 206–222. [Google Scholar] [CrossRef]
  52. Guo, Y.; Liu, T.-W.; Na, J.; Fung, R.-F. Envelope order tracking for fault detection in rolling element bearings. J. Sound Vib. 2012, 331, 5644–5654. [Google Scholar] [CrossRef]
  53. Yang, Y.; Wang, H.; Cheng, J.; Zhang, K. A fault diagnosis approach for roller bearing based on VPMCD under variable speed condition. Measurement 2013, 46, 2306–2312. [Google Scholar] [CrossRef]
  54. Farhat, M.; Chiementin, X.; Chaari, F.; Bolaers, F.; Haddar, M. Order-Based Identification of Bearing Defects under Variable Speed Condition. Appl. Sci. 2021, 11, 3962. [Google Scholar] [CrossRef]
  55. Appana, D.K.; Prosvirin, A.; Kim, J.-M. Reliable fault diagnosis of bearings with varying rotational speeds using envelope spectrum and convolution neural networks. Soft Comput. 2018, 22, 6719–6729. [Google Scholar] [CrossRef]
  56. Xu, L.; Chatterton, S.; Pennacchi, P.; Liu, C. A Tacholess Order Tracking Method Based on Inverse Short Time Fourier Transform and Singular Value Decomposition for Bearing Fault Diagnosis. Sensors 2020, 20, 6924. [Google Scholar] [CrossRef]
  57. Luo, J.; Zhang, S.; Zhong, M.; Lin, Z. Order spectrum analysis for bearing fault detection via joint application of synchrosqueezing transform and multiscale chirplet path pursuit. Shock. Vib. 2016, 2016, 2976389. [Google Scholar] [CrossRef] [Green Version]
  58. Gao, Y.; Guo, Y.; Chi, Y.L.; Qin, S.R. Order tracking based on robust peak search instantaneous frequency estimation. In Proceedings of the International Symposium on Instrumentation Science and Technology, Harbin, China, 8–12 August 2006; Journal of Physics: Conference Series. IOP Publishing: Bristol, UK, 2006; Volume 48, p. 091. [Google Scholar] [CrossRef]
  59. Blough, J.R. Improving the Analysis of Operating Data on Rotating Automotive Components. Ph.D. Thesis, University of Cincinnati, Cincinnati, OH, USA, 1998. [Google Scholar]
  60. Brandt, A. Rotating Machinery Analysis. In Noise and Vibration Analysis: Signal Analysis and Experimental Procedures; John Wiley & Sons: Chichester, UK, 2011. [Google Scholar]
  61. Gao, S.; Chatterton, S.; Pennacchi, P.; Chu, F. Behaviour of an angular contact ball bearing with three-dimensional cubic-like defect: A comprehensive non-linear dynamic model for predicting vibration response. Mech. Mach. Theory 2021, 163, 104376. [Google Scholar] [CrossRef]
  62. Gao, S.; Chatterton, S.; Naldi, L.; Pennacchi, P. Ball bearing skidding and over-skidding in large-scale angular contact ball bearings: Nonlinear dynamic model with thermal effects and experimental results. Mech. Syst. Signal Processing 2021, 147, 107120. [Google Scholar] [CrossRef]
  63. Prashad, H. Solving Tribology Problems in Rotating Machines; Woodhead Publishing: Sawston, UK, 2006. [Google Scholar]
  64. Han, Q.; Li, X.; Chu, F. Skidding behavior of cylindrical roller bearings under time-variable load conditions. Int. J. Mech. Sci. 2018, 135, 203–214. [Google Scholar] [CrossRef]
  65. Li, F.; Deng, S.; Zhang, W. Skidding Characteristics of Ball Bearings Withstanding a Frequent Oscillating Condition. Chin. J. Mech. Eng. 2021, 57, 168–178. [Google Scholar]
  66. Pennacchi, P.; Bruni, S.; Chatterton, S.; Borghesani, P.; Ricci, R.; Marinis, D.; Didonato, A.; Unger-Weber, F. A test rig for the condition-based maintenance application on the traction chain of very high speed trains. In Proceedings of the 9th World Congress on Railway Research, Lille, France, 22–26 May 2011; Rail Safety and Standards Board Ltd.: London, UK, 2011; pp. 1–11. [Google Scholar]
Figure 1. Flow chart of proposed methodology.
Figure 1. Flow chart of proposed methodology.
Sensors 22 02026 g001
Figure 2. Order maps’ computation procedure.
Figure 2. Order maps’ computation procedure.
Sensors 22 02026 g002
Figure 3. Order map under variable speed (3525–4125 rpm).
Figure 3. Order map under variable speed (3525–4125 rpm).
Sensors 22 02026 g003
Figure 4. Order maps of bearings with different types of defects at different speeds and loads.
Figure 4. Order maps of bearings with different types of defects at different speeds and loads.
Sensors 22 02026 g004
Figure 5. Architecture of proposed convolutional neural network.
Figure 5. Architecture of proposed convolutional neural network.
Sensors 22 02026 g005
Figure 6. Test rig-1 layout.
Figure 6. Test rig-1 layout.
Sensors 22 02026 g006
Figure 7. Healthy and defected bearings considered in case study-1.
Figure 7. Healthy and defected bearings considered in case study-1.
Sensors 22 02026 g007
Figure 8. Raw vibration signal in the time domain.
Figure 8. Raw vibration signal in the time domain.
Sensors 22 02026 g008
Figure 9. Test rig-2 layout. (a) main components; (b) under study bearings’ position; (c) accelerometer’s position.
Figure 9. Test rig-2 layout. (a) main components; (b) under study bearings’ position; (c) accelerometer’s position.
Sensors 22 02026 g009
Figure 10. Different types of seeded defects of ball bearings (case study-2).
Figure 10. Different types of seeded defects of ball bearings (case study-2).
Sensors 22 02026 g010
Figure 11. Different types of seeded defects of cylindrical roller bearings (case study-2).
Figure 11. Different types of seeded defects of cylindrical roller bearings (case study-2).
Sensors 22 02026 g011
Figure 12. Performance comparison when models were trained at 1000 rpm + 5 KN load and tested under all combinations of speed and load.
Figure 12. Performance comparison when models were trained at 1000 rpm + 5 KN load and tested under all combinations of speed and load.
Sensors 22 02026 g012
Figure 13. Confusion matrix: model trained at 3170 rpm and tested at 4955 rpm.
Figure 13. Confusion matrix: model trained at 3170 rpm and tested at 4955 rpm.
Sensors 22 02026 g013
Figure 14. Confusion matrix: model trained at variable speed between 1000 and 2000 rpm and tested at variable speed between 3525 and 4125 rpm.
Figure 14. Confusion matrix: model trained at variable speed between 1000 and 2000 rpm and tested at variable speed between 3525 and 4125 rpm.
Sensors 22 02026 g014
Table 1. Different operating conditions considered in case study-1.
Table 1. Different operating conditions considered in case study-1.
Operating Condition NumberOperating ConditionsOperating Condition NumberOperating Conditions
Speed (rpm)Load (KN)Speed (rpm)Load (KN)
1100056200015
2100010730005
31000158300010
4200059300015
5200010
Table 2. Bearing defect classes for case study-1.
Table 2. Bearing defect classes for case study-1.
Class123
Health State/Type of FaultNormal/HealthyCombined Defects/Multiple DefectsOuter Race Defect
Table 3. Different operating conditions considered in case study-2.
Table 3. Different operating conditions considered in case study-2.
Operating Condition NumberType of Operating ConditionSpeed Range (rpm)Time Duration
1Constant speed3170100 s
2Constant speed4955100 s
3Variable speed (acceleration)1000–2000100 s
4Variable speed (acceleration)3525–4125100 s
5Variable speed (deceleration)4050–2560100 s
Table 4. Bearings defect classes considered in case study-2.
Table 4. Bearings defect classes considered in case study-2.
ClassDefect Type/Condition
Ball Bearing (BB)Roller Bearing (RB)
1Normal/HealthyNormal/Healthy
2Ball defectCombined defects
3Cage defectInner race defect
4Inner race defectOuter race defect
5Outer race defectRoller defect
Table 5. Performance and comparison of the proposed methodology in terms of prediction accuracy for defect diagnosis under different operating conditions (case study-1).
Table 5. Performance and comparison of the proposed methodology in terms of prediction accuracy for defect diagnosis under different operating conditions (case study-1).
TrainingTestingPrediction Accuracy (%)
Speed (rpm)Load (KN)Speed (rpm)Load (KN)CNN + Order Maps (Proposed)CNN + Spectrograms [38,39,40]KNN + GA [12]ANNSVMKNN + Order Maps
100051000510010098.399.299.698.3
10001010010087.989.188.298
10001510010088.890.490.297.9
2000599.867.271.556.364.895.4
20001097.458.551.54963.591.3
20001598.162.446.943.158.193.1
3000598.463.365.640.664.292.7
300010964665.237.365.289.7
30001596.759.159.835.458.391.4
100010100051001008180.180.297.5
10001010010098.89499.898.8
1000151001008480.882.598.2
20005975862.167.756.392.5
20001098.175.766.766.765.293
20001597.448.658.162.543.190.4
3000595.364.660.268.666.587.3
30001096.272.161.167.156.590.7
30001596.557.352.365.642.788.5
1000151000510010090.2189.891.597.5
10001010099.391.789.491.397
10001510010097.596.999.498.2
2000597.146.163.143.160.691.3
20001098.457.766.0459.664.292
20001598.878.14836.53990.7
3000597.257.564.651.551.589.6
30001096.95157.152.752.784.6
3000159962.243.9637.739.285
200051000598.162.845.84953.888.3
10001096.856.451.535.842.391.8
10001597.148.650.639.449.888.1
2000510010097.198.399.699.6
20001010010084.586.485.199.2
20001599.499.386.8787.187.898.8
3000599.47472.373.67297.6
30001098.17065.665.666.995.7
30001597.67068.366.76896.3
2000101000597.46743.131.432.786
10001098.669.27453931.994.6
10001598.151.742.336.23083.9
2000510098.283.281.280.393.4
20001010010099.6100100100
20001599.698.48688.387.598.1
300059966.747071.374.394.6
30001099.272.668.571.96985.7
3000159870.869.470.168.988.8
2000151000597.373.451.351.245.681.2
10001096.873.447.345.143.882.5
10001599.177.553.349.54383.2
2000599.398.487.9848295.8
20001099.399.390.6285.88196.4
20001510010098.599.299.8100
3000594.746.771.364.277.186.3
3000109749.861.364.271.590.3
30001598.153.973.279.279.885.7
300051000598.463.5741.324.227.189.3
10001095.455.345.840.83189.6
10001597.255.436.337.325.485.4
2000598.97471.274.858.384.2
20001098.168.8566.767.977.397.6
20001595.546.2571.775.475.891.3
30005100100999999.8100
30001099.810087.387.486.299.3
30001510010088.287.887.399.6
3000101000596.157.8326.219.41485.7
10001097.35035.812.310.682.1
1000159528.523.51611.2586.1
2000598.569.160.856.953.387.6
20001099.172.362.767.86696.9
20001597.151.661.566.558.385.5
3000599.698.285.682.983.899.3
30001010010094.692.599.4100
30001510099.684.48985.6100
3000151000597.138.822.527.330.483
10001096.938.935.833.533.181.1
10001599.148.929.431.729.884.3
2000597.673.266.96563.589.2
20001098.972.873.177.378.798.5
20001598.975.87475.675.485.3
3000599.210090.689.582.899.8
30001010010083.684.786.3100
30001510010094.295.899.4100
Overall Average Prediction Accuracy98.473.767.464.765.392.3
Table 6. Performance of the proposed methodology in terms of prediction accuracy for defect diagnosis under variable speed (case study-2).
Table 6. Performance of the proposed methodology in terms of prediction accuracy for defect diagnosis under variable speed (case study-2).
TrainingTestingPrediction Accuracy (%)
Ball BearingCylindrical Roller Bearing
3170 rpm
(Constant Speed)
3170 rpm
(Constant speed)
100100
4955 rpm
(Constant speed)
10099.9
1000–2000 rpm
(acceleration)
99.197.2
3525–4125 rpm
(acceleration)
99.496.8
4050–2560 rpm
(deceleration)
99.196.2
1000–2000 rpm
(acceleration)
3170 rpm
(Constant speed)
98.799.1
4955 rpm
(Constant speed)
98.598.1
1000–2000 rpm
(acceleration)
100100
3525–4125 rpm
(acceleration)
99.497.3
4050–2560 rpm
(deceleration)
9895.1
Average Prediction accuracy for each bearing99.298
Overall average Prediction accuracy98.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tayyab, S.M.; Chatterton, S.; Pennacchi, P. Intelligent Defect Diagnosis of Rolling Element Bearings under Variable Operating Conditions Using Convolutional Neural Network and Order Maps. Sensors 2022, 22, 2026. https://doi.org/10.3390/s22052026

AMA Style

Tayyab SM, Chatterton S, Pennacchi P. Intelligent Defect Diagnosis of Rolling Element Bearings under Variable Operating Conditions Using Convolutional Neural Network and Order Maps. Sensors. 2022; 22(5):2026. https://doi.org/10.3390/s22052026

Chicago/Turabian Style

Tayyab, Syed Muhammad, Steven Chatterton, and Paolo Pennacchi. 2022. "Intelligent Defect Diagnosis of Rolling Element Bearings under Variable Operating Conditions Using Convolutional Neural Network and Order Maps" Sensors 22, no. 5: 2026. https://doi.org/10.3390/s22052026

APA Style

Tayyab, S. M., Chatterton, S., & Pennacchi, P. (2022). Intelligent Defect Diagnosis of Rolling Element Bearings under Variable Operating Conditions Using Convolutional Neural Network and Order Maps. Sensors, 22(5), 2026. https://doi.org/10.3390/s22052026

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