Next Article in Journal
Accurate Global Point Cloud Registration Using GPU-Based Parallel Angular Radon Spectrum
Previous Article in Journal
GaitSG: Gait Recognition with SMPLs in Graph Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive VMD–K-SVD-Based Rolling Bearing Fault Signal Enhancement Study

1
School of Mechanical Engineering and Mechanics, Xiangtan University, Xiangtan 411105, China
2
Engineering Research Center of Complex Trajectory Machining Process and Equipment, Ministry of Education, Xiangtan University, Xiangtan 411105, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(20), 8629; https://doi.org/10.3390/s23208629
Submission received: 9 September 2023 / Revised: 14 October 2023 / Accepted: 20 October 2023 / Published: 22 October 2023
(This article belongs to the Section Fault Diagnosis & Sensors)

Abstract

:
To address the challenges associated with nonlinearity, non-stationarity, susceptibility to redundant noise interference, and the difficulty in extracting fault feature signals from rolling bearing signals, this study introduces a novel combined approach. The proposed method utilizes the variational mode decomposition (VMD) and K-singular value decomposition (K-SVD) algorithms to effectively denoise and enhance the collected rolling bearing signals. Initially, the VMD method is employed to separate the overall noise into intrinsic mode functions (IMFs), reducing the noise content within each IMF. To optimize the mode component, K, and the penalty factor, α, in VMD, an improved arithmetic optimization algorithm (IAOA) is employed. This ensures the selection of optimal parameters and the decomposition of the signal into a set of IMFs, forming the original dictionary. Subsequently, the signals are decomposed into multiple IMFs using VMD, and an original dictionary is constructed based on these IMFs. K-SVD is then applied to the original dictionary to further reduce the noise in each IMF, resulting in a denoised and enhanced signal. To validate the efficacy of the proposed method, rolling bearing signals collected from Case Western Reserve University (CWRU) and thrust bearing test rigs were utilized. The experimental results demonstrate the feasibility and effectiveness of the proposed approach in denoising and enhancing the rolling bearing signals.

1. Introduction

Rolling bearings play a crucial role as key support elements in large-scale industrial production. The health and proper functioning of these bearings directly impact the normal operation of specific organizations [1,2,3]. Consequently, assessing the operational reliability of rolling bearings, monitoring their operating status, and conducting fault diagnosis are essential for ensuring the safe and stable operation of machinery and equipment. By implementing predictive maintenance strategies for rolling bearings, productivity can be improved, and the likelihood of production accidents can be reduced [4,5,6,7,8]. To achieve this, it is necessary to conduct comprehensive assessments of rolling bearing health, monitor their performance, and diagnose faults promptly.
Signal-based time–frequency domain analysis is widely employed as a means of diagnosing faults in rolling bearings. The signals generated during mechanical operation encompass a wealth of information pertaining to equipment performance; through signal analysis, it becomes possible to discern the operational state of rolling bearings [1,9]. Among the various types of signals used in rolling bearing analysis, such as vibration signals, acoustic signals, acceleration signals, and speed signals, vibration and acoustic signals are particularly advantageous due to their ability to capture subtle and higher-scale information. Hence, time–frequency domain analysis based on vibration and acoustic signals offers enhanced accuracy and sensitivity [10,11,12]. However, the vibration and acoustic signals collected often contain a considerable amount of redundant noise. Therefore, it is crucial to minimize the noise present in the raw acoustic signals and amplify the fault characteristic signals to ensure accurate fault diagnosis and analysis.
Variational mode decomposition (VMD), introduced by Dragomiretskiy and Zosso in 2014 [13], is a novel signal decomposition method that offers comparable performance to other methods such as empirical mode decomposition (EMD) [14], complete ensemble empirical mode decomposition (CEEMD) [15], and complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) [16]. VMD stands out due to its enhanced stability, robustness, and the avoidance of mode aliasing. Although the effectiveness of VMD is influenced by the mode component, K, and the penalty factor, α, optimizing these parameters using heuristic intelligent optimization algorithms has been recognized as an effective approach [17]. Numerous researchers have employed various intelligent optimization algorithms to fine-tune the parameters of VMD and analyze the resulting IMF to assess the operational status of bearings. Zhang et al. [18] investigated the mode component, K, and penalty factor, α, of VMD using correlation coefficients and energy ratios. They demonstrated the signal decomposition capability of VMD by decomposing both simulated and actual vibration signals under appropriate parameter settings. Zhang et al. [19] utilized the grasshopper optimization algorithm (GOA) to optimize the parameters in VMD. According to the research conducted by the aforementioned scholars, VMD has been identified as a viable signal processing method for decomposing bearing signals. However, the effectiveness of VMD’s decomposition is influenced by certain parameters. To enhance the decomposition performance to its fullest potential, this paper proposes the utilization of the IAOA to optimize the relevant parameters of VMD. The primary objective is to achieve a more comprehensive decomposition of bearing signals by fine-tuning these parameters.
K-singular value decomposition (K-SVD) is an algorithm originally proposed by M. Aharon [20] for designing an ultra-complete dictionary suitable for sparse representation. In comparison with other algorithms for signal sparse representation, K-SVD employs an iterative process of sparse coding and dictionary atom updating based on the current set of example signals. This iterative process leads to a better fit to the data; the dictionary columns are updated with the sparse representations, resulting in accelerated convergence. One of the advantages of K-SVD is its adaptability, as it can be combined with other tracking methods. Lebrun et al. [21] utilized K-SVD for image denoising and achieved favorable results, demonstrating the algorithm’s stability and noise reduction capabilities. Anaraki et al. [22] focused on dictionary construction in K-SVD and utilized compressed sensing measurements based on signals for sparse representation, showcasing the algorithm’s proficiency in noise reduction. Wu et al. [23] constructed dictionaries using different fabric samples and successfully reconstructed the samples using K-SVD, thus confirming the effectiveness of the algorithm. Wang et al. [24] proposed a method that combines the orthogonal matching pursuit (OMP) algorithm with K-SVD, illustrating the feasibility of K-SVD for noise reduction in signals from rotating machinery, such as rolling bearings. Based on the research conducted by the aforementioned scholars, K-SVD has been recognized as an effective method for noise reduction and signal decomposition. The decomposition effect of K-SVD is influenced by the choice of the original dictionary. To optimize the decomposition effect and enhance its performance, this paper proposes a novel approach. It utilizes the matrix composed of IMFs obtained through VMD as the original dictionary for K-SVD. By employing the IMFs from VMD as the initial dictionary, the aim is to improve the noise reduction capability of K-SVD to its maximum potential.
In conclusion, K-SVD is a powerful algorithm for sparse representation and has demonstrated good stability and noise reduction capabilities in various applications, including signal and image denoising. Its effectiveness is often limited by the constraints of the initially constructed dictionary. On the other hand, VMD is a novel signal processing method that offers advantages in signal decomposition and noise reduction. It overcomes issues such as endpoint effects and mode aliasing, resulting in its wide use in signal processing. However, the decomposition performance of VMD is influenced by the mode component, K, and the penalty factor, α. In this paper, an approach is proposed to optimize the relevant parameters of VMD using the IAOA. The IMFs obtained from VMD are then used to construct the original dictionary in K-SVD. This approach ensures the optimal performance of VMD while leveraging the noise reduction and enhancement capabilities of K-SVD. The proposed method was validated using experimental signals from rolling bearings obtained from Case Western Reserve University and Thrust Bearing Test Bed.

2. VMD–K-SVD-Based Signal Enhancement Method for Fault Characterization

2.1. VMD Basic Principle

The operation of the VMD algorithm can be understood as a variational process aimed at decomposing the original signal into IMFs. Each IMF represents a bandwidth-limited AM (amplitude modulation) and FM (frequency modulation) function. The expression for an IMF can be written as follows:
μ k ( t ) = A k ( t ) cos ( φ k ( t ) )
where A k ( t ) is the instantaneous amplitude of U k ( t ) and ϕ k ( t ) is the phase.
The VMD algorithm involves the construction and solution of variational problems. The construction of variational problems in VMD can be described as follows:
(1)
The IMF components are subjected to a Hilbert transform, which introduces a unit pulse signal, δ ( t ) . This process allows us to obtain the expression of the analyzed signal with K IMF components. The expression for the analyzed signal with K IMF components can be represented using Equation (2):
F k = ( δ t + j π t ) μ k ( t )
(2)
The resolved signal of the K IMF components is shifted to the baseband, and the center frequency of each IMF is adjusted using an exponential term, e j w k t . This process results in an expression for the shifted spectrum, which can be represented as shown in Equation (3):
H k = [ ( δ t + j π t ) μ k ( t ) ] e j ω k t
where ω k ( t ) is the center frequency of each IMF component.
(3)
The bandwidth of the signal is estimated using Gaussian smoothness, which involves calculating the square of the two-parameter gradient. This estimation helps in obtaining the expression for the constrained variational model. The constrained variational model expression can be represented as shown in Equation (4):
{ min { μ k } , { ω k } { k = 1 K t [ ( δ ( t ) + j π t ) × μ k ( t ) ] e j ω k t 2 2 } s . t . k = 1 K μ k ( t ) = f ( t )
In the construction of variational problems in VMD, a penalty factor is introduced, and an augmented generalized Lagrange function is constructed to solve for the optimal solution of each component. The expression of the Lagrange function can be represented as follows:
L ( { μ k } , { ω k } , λ ) = α k t [ δ ( t ) + j / π t × u k ( t ) ] e j w k t 2 2 + f ( t ) K = 1 K μ ( t ) 2 2 + λ ( t ) , f ( t ) k u k ( t )
The computation of the variational problem, as described by the Lagrange function, involves using the alternating direction multiplier method. This method produces different values for the mode components, { μ k } , and the center frequency, { ω k } , during the updating process. The iterative formulas for updating the mode components, { μ k } , and the center frequency, { ω k } , are shown in Equation (6) and Equation (7), respectively. After updating the values of these two parameters, the multiplier value λ ^ is updated according to Equation (8). Subsequently, Equation (9) is verified to determine if it holds. When Equation (9) is satisfied, update of the multiplier operator, λ ^ , is completed.
μ ^ k n + 1 ( ω ) = f ^ ( ω ) i k μ ^ i ( ω ) + λ ^ ( ω ) 2 1 + 2 α ( ω ω k ) 2
ω k n + 1 = 0 ω | μ ^ k ( ω ) | 2 d ω 0 | μ ^ k ( ω ) | 2 d ω , k { 1 , 2 , 3 K }
λ ^ n + 1 ( ω ) = λ ^ n ( ω ) + τ ( f ^ ( ω ) k = 1 K μ ^ k n + 1 ( ω ) )
k = 1 K μ ^ k n + 1 μ ^ k n 2 2 / k = 1 K μ ^ 2 2 < ε
where μ ^ k n + 1 ( ω ) denotes the Wiener filtering of the current signal; the real part is { μ k ( t ) } after Fourier transform.
Indeed, finding the optimal parameter combination that matches the signal being analyzed is crucial for the VMD method. The effectiveness and performance of VMD heavily rely on selecting appropriate parameter values for the mode component, K, and penalty factor, α.

2.2. K-SVD Basic Principle

The K-SVD algorithm is an iterative process that alternates between sparse coding and dictionary updating. Its core idea is to update the atoms of the dictionary one by one. The objective optimization function of the K-SVD algorithm can be represented as follows:
arg min D , X Y D X F 2 , s . t . i , x j T 0
where Y is the sample matrix formed by the original signal; D is the learning dictionary matrix; X is the sparse coefficients matrix; x j is the jth row of X ; and T 0 is the sparsity.
After fixing the dictionary, D , and updating the sparse coefficients, x j , using the OMP algorithm, the next step in the K-SVD algorithm is to update the dictionary column by column. This involves calculating the error matrix of the updated atoms, E k . The error matrix can be calculated as follows:
Y D X F 2 = Y j = 1 k d j x j F 2 = ( Y j k d j x j ) d k x k F 2 = E k d k x k F 2
E k = Y j k a j d T j
where d j is the jth column of the dictionary, D ; x j is the j th row of the sparse matrix, X ; d k is the k th column of the dictionary, D ; and x k is the k th row of the sparse matrix, X . Determine the set of positions, Ω , of all non-zero entries, x j , in the sparse matrix, X , and according to Ω , select the corresponding column of E in E k , and perform the singular value decomposition of E , i.e.:
E = U Δ V T
The atoms in the dictionary matrix are updated using the first column of the matrix, U . The value of the first column of the matrix, V , multiplied by the corresponding sparse coefficient, x j , is used to update the atoms in the dictionary. This process is repeated for each column of the matrix, completing the update of the sparse representation dictionary.

2.3. Improved Arithmetic Optimization Algorithm (IAOA)

In the paper by Abualigah [25], the arithmetic optimization algorithm (AOA) is introduced as a new heuristic optimization algorithm based on physical methods. The AOA incorporates the concept of arithmetic operators as an exploration strategy, which enhances the algorithm’s global exploration capability and helps avoid trapping in local optima to some extent. It has been found to be more effective compared with other algorithms such as grey wolf optimization (GWO), particle swarm optimization (PSO), cuckoo search (CS), etc. In this study, AOA was used as the base algorithm, and four arithmetic operators, namely, multiplication (M), division (D), subtraction (S), and addition (A), were employed to further enhance the exploration ability of the algorithm. Similarly to most heuristic optimization algorithms, the AOA operates by generating a set of random solutions (X) and iteratively updating these solutions through the algorithm’s exploration process until an optimal solution is achieved. The stochastic expression for the solution is represented in Equation (14):
X = [ x 1 , 1 x 1 , 2 x 1 , n x 2 , 1 x 2 , 1 x 2 , n x N , 1 x N , 2 x N , 3 x N , n ]
where N is the number of important parameters in the problem to be solved by the algorithm and n is the number of predefined populations.
The exploration capability in the arithmetic optimization algorithm (AOA) is regulated by the mathematical optimizer acceleration (MOA) function. The MOA function is determined by a linear time, t, a predefined minimum value, MOAmin, a predefined maximum value, MOAmax, and the maximum number of iterations, T. The formula for calculating the MOA function is shown in Equation (15):
M O A ( t ) = M O A min + t × ( M O A max M O A min T )
In addition to the MOA function, another key parameter in AOA is the mathematical optimizer probability (MOP). The MOP parameter is utilized to balance the exploration and exploitation aspects of the algorithm. It is defined by the equation shown in Equation (16):
M O P ( t ) = 1 ( t T ) 1 α
where α is the sensitivity parameter.
In the AOA, the exploration behavior differs from other optimization algorithms. The AOA uses the M operator or the D operator to jointly determine its exploration behavior. On the other hand, the S operator and the A operator are employed to generate high-density results. The outcomes obtained from different operators exhibit distinct characteristics, ensuring the effectiveness of the AOA’s exploration capability and exploration behavior. To determine the selection of operators during the exploration process, random values r 1 , r 2 , and r 3 , which range from 0 to 1, are introduced. These random values are used to guide the decision-making process and determine which operator to apply at a particular stage of exploration. By incorporating randomization, the AOA introduces a stochastic element in the selection of operators, enhancing the algorithm’s exploration capability, and promoting diverse exploration paths.
If M O A < r 1 , M and D operators are used in the AOA, the position update rule during this search phase is described by Equation (17).
x i , j ( t + 1 ) = { x b e s t ( t ) × M O P × ( ( u b j l b j ) × μ + l b j ) ,             else x b e s t ( t ) ÷ ( M O P + ε ) × ( ( u b j l b j ) × μ + l b j ) ,   r 2 < 0.5
where x i , j ( t + 1 ) is the calculated new position, x b e s t ( t ) is the best position obtained so far, ε is the epsilon number used to avoid division by zero, l b j and u b j denote the lower and upper bounds of the jth dimension of the variable, respectively, and μ is a parameter empirically set to 0.5.
If M O A > r 1 is used, it is calculated using the S operator and the A operator. When r 3 < 0.5 , the S operator is used; vice versa, the A operator is used. The specific position update rule during this search phase is represented by Equation (18):
x i , j ( t + 1 ) = { x b e s t ( t ) + M O P × ( ( u b j l b j ) × μ + l b j ) ,     e l s e x b e s t ( t ) M O P × ( ( u b j l b j ) × μ + l b j ) ,   r 3 < 0.5
The initialization of the population is a critical step in the AOA that can significantly impact the algorithm’s performance. Random distribution in the initialization process can lead to an uneven distribution of individuals in the initial population, which may increase the risk of the algorithm becoming trapped in local optima and limit its global search ability. To address this issue, chaotic variables are introduced in the population initialization phase. Chaotic variables exhibit characteristics such as randomness, ergodicity, and regularity, which can enhance the exploration capability of the algorithm [26]. In this study, the combination of three types of chaotic mapping functions, namely, logistic, sine, and cosine, was employed. This combination, known as logistic–sine–cosine (LSC) mapping, offers more complex chaotic performance [27]. By utilizing LSC mapping to initialize the population in the AOA, the algorithm benefits from the enhanced exploration capabilities provided by chaotic variables. This approach aims to improve the algorithm’s ability to escape local optima and enhance its global search ability, ultimately leading to better optimization results.
The expression for logistic mapping is described by Equation (19):
x i + 1 = 4 r x i ( 1 x i ) r [ 0 , 1 ]
The expression for sine mapping is described by Equation (20):
x i + 1 = r sin ( π x i ) r [ 0 , 1 ]
The expression for tent mapping is described by Equation (21):
x i + 1 = { 2 r x i 2 r ( 1 x i ) x i < 0.5 x i 0.5
The expression for logistic–sine–cosine mapping is described by Equation (22):
x i + 1 = cos ( π ( 4 r 4 x i ( 1 x i ) + ( 1 r 4 ) sin ( π x i ) 0.5 ) )   r 4 [ 0 , 1 ]
In the context of opposition-based learning (OBL), a learning strategy that directly affects population initialization, there is a technique known as the inverse approximation of the number of opposites. This technique aims to address the issue of uneven distribution of population positions more effectively while reducing computational complexity.
We define a hypothesis, x [ a , b ] , whose inverse proximal opposites, x q r , are defined as follows:
x q r = r a n d ( c , x )
where c = ( a + b ) / 2 is the center of [ a , b ] [28].
To expedite the convergence towards the optimal solution and achieve a more uniform distribution of solution seekers in the population initialization of the AOA, a multi-step process is employed. Initially, a set of initial values is generated using logistic–sine–cosine (LSC) chaotic mapping. Then, contrastive learning is applied to generate a new set of samples using the initial values as references. These samples are substituted into a specific function, and those with the minimum fitness value are continuously recorded. This process yields a set of positions that are closer to the optimal solution compared with the previous methods. These positions are encoded into variable X and utilized as the initial solution seekers for subsequent computations in the algorithm.
The exploration behavior of different optimization algorithms can vary, leading to different optimization processes for various functions. To further enhance the optimization capability of AOA, an adaptive Gaussian distribution is introduced to improve its exploration process [29]. When the global optimal solution of the AOA remains unchanged for a prolonged period, indicating that further computation and improvement are unlikely to yield better results, Gaussian distribution is employed to generate a new solution. This helps to introduce diversity and exploration into the search process. The specific equation depicting the usage of Gaussian distribution for creating a new solution is presented in Equation (24).
x i n e w = { G a u s s i a n ( x i , j , μ , σ ) + ( x i , j r 5 × x i ) , r 6 < 0.5 G a u s s i a n ( x i , j , μ , σ ) , else }
where r 5 and r 6 are the random values in [0, 1], G a u s s i a n ( x i , j , μ , σ ) is the Gaussian distribution, x i is the current position, μ is the mean, and σ is the standard deviation. In order to improve the exploration ability of the algorithm, the Gaussian distribution is no longer calculated with the global optimal position as the sample, but with all the global positions as the reference and redistributed, so as to obtain a new set of global variables again, which, in turn, improves the exploration ability of the algorithm.
The flow of the improved artificial organism algorithm (IAOA) is illustrated in Figure 1.

2.4. Flow Chart for the Signal Enhancement of Rolling Bearing Fault Characteristics

In this study, a fault feature signal enhancement process is proposed, which involves several key steps. Firstly, the acquired raw signals are subjected to decomposition using VMD. To optimize the VMD parameters, specifically the mode component, K, and the penalty factor, α, the IAOA is employed. The IAOA ensures effective parameter optimization for VMD.
The resulting IMFs obtained from the VMD serve as the original dictionary, capturing the underlying fault features present in the signal. Subsequently, the K-SVD algorithm, a widely used dictionary learning technique, is applied for noise reduction on the original dictionary. Through iterative updates of the dictionary and sparse coefficients, K-SVD achieves noise reduction and signal enhancement.
The noise-reduced dictionary, obtained from the K-SVD process, is then utilized to reconstruct the signal by multiplying it with the sparse coefficient matrix. This multiplication effectively combines the dictionary elements, weighted by the sparse coefficients, resulting in the reconstruction of a noise-reduced and enhanced signal.
Overall, the proposed fault feature signal enhancement process, as depicted in Figure 2, involves the decomposition of raw signals using VMD with parameter optimization through IAOA. The resulting IMFs form the original dictionary, which undergoes noise reduction using K-SVD. The reconstructed signal, obtained through multiplication of the noise-reduced dictionary with the sparse coefficients, represents an enhanced version of the fault feature signal.

3. Experimentation and Analysis

3.1. Based on Case Western Reserve University Data Analysis

The proposed method in this paper was validated using experimental vibration signal data from the Bearing Signal Library Center of Case Western Reserve University (CWRU) [30]. The data pertained to the SKF 6205-2RSJEM SKF bearing, with a motor speed of 1797 r/min. The analyzed signal represented the inner ring failure signal, with a loss diameter of 0.1778 mm. The vibration signals were sampled at a frequency of 12 kHz for a duration of 1s. The inner ring failure frequency, f i , was determined to be 162 Hz, while the bearing rotation frequency, f o , was measured at 30 Hz.
Figure 3 displays the time–domain representation of the fault signal. Due to a significant amount of noise present in the original signal, an initial step involved decomposing the original signal using the IAOA–VMD method. The iterative convergence of the fitness function for the IAOA–VMD algorithm is depicted in Figure 4. It is evident from the graph that the IAOA algorithm, enhanced by LSC mapping and adversarial learning, achieves the optimal value during the second iteration. The optimized VMD parameters consist of a mode component, K, value 12, and a penalty factor, α, value 7250.
Under the premise of guaranteeing parameter optimization, the VMD is decomposed to obtain 12 IMF components. The feature matrix is constructed with this IMF component to obtain a 12 × 12,000 feature matrix. The obtained feature matrix is used as the original dictionary for K-SVD. By multiplying the obtained dictionary and sparse coefficient matrix to obtain the signal components after noise reduction, the method is used to realize the noise reduction and enhancement of the bearing fault feature signal. The time–domain diagram of the enhanced signal is shown in Figure 5, from which it can be seen that the noise has been significantly reduced.
To assess the effectiveness of noise reduction and signal enhancement, the squared-envelope spectrum analysis of the faulty signal processed by VMD–K-SVD is performed and compared with the original signal and the signal processed by VMD. The results are depicted in Figure 6a–c. Additionally, to evaluate the impact of different signal decomposition methods combined with K-SVD, the EMD and CEEMDAN algorithm are employed in conjunction with the K-SVD algorithm to perform noise reduction of the vibration signals. The squared-envelope spectrum analysis is then applied; the outcomes are presented in Figure 6d,e. Figure 6a illustrates that, aside from the fault eigenfrequency, the signal contains a significant amount of irrelevant noise that can potentially interfere with fault diagnosis outcomes. In Figure 6b, the squared-envelope spectrum following VMD processing exhibits more prominent bearing rotation frequency, rotation multiplier frequency, and inner ring fault frequency when compared with the original signal. However, the amplitude of the inner ring fault remains less conspicuous in comparison to the rotation frequency. Figure 6c demonstrates that, in contrast to the original signal and the squared-envelope spectrum after VMD processing, the squared-envelope spectrum subsequent to VMD–K-SVD processing exhibits more pronounced bearing rotation frequency, multiplicative frequency, and fault frequency. A comparison of the amplitude of the processed fault characteristic frequency and the rotation frequency allows for the determination of the bearing’s operational status and fault type.

3.2. Signal Analysis of Rolling Bearings in a Thrust Bearing Test Rig

We considered the limitations of the bearing dataset available in Case Western Reserve University (CWRU) and, in order to demonstrate the applicability of the proposed method to signals with high sampling frequency and high noise, this study focused on acoustic signals generated during the operation of a thrust bearing test bench. The VMD–K-SVD method was employed to perform noise reduction and enhancement on these signals.
The thrust bearing test rig, depicted in Figure 7a, was primarily designed to measure changes in oil film pressure and oil film temperature of large plain bearings during operation. The test rig incorporated 29412E type roller bearings in its bearing box. Acoustic signals from this rolling bearing were collected by placing an acoustic signal sensor on the outer side of the box, as illustrated in Figure 7b. The specific model and style of the sensor are presented in Figure 7c. During the experiment, the motor speed was set to 2970 r/min, and the acoustic signals were sampled at a frequency of 1000 kHz for a duration of 1s. This resulted in a dataset containing 1,000,000 data points, as shown in Figure 7d. To ensure the experiment’s authenticity, a naturally damaged rolling bearing was used, reflecting wear and tear accumulated during long-term usage. Table 1 provides the relevant parameters associated with the rolling bearing. Based on these parameters, the rotation frequency, f o , and the outer ring failure frequency, f i , of the rolling bearing were calculated to be 49.5 Hz and 323 Hz, respectively.
Figure 8 presents the time–domain diagram of the acquired acoustic signals from the rolling bearing. Comparison between this diagram and the time–domain diagram of the vibration signals shown in Figure 3 reveals that the acoustic signals in Figure 8 contain a significant amount of impact noise and periodic noise. The presence of such noise components can significantly impact the accuracy and reliability of subsequent analysis and fault diagnosis processes.
To reduce noise and enhance the fault signal in the acoustic signal of the rolling bearing, an adaptive decomposition and reduction process was performed using IAOA–VMD. The iterative diagram of the fitness function for the IAOA–VMD algorithm is depicted in Figure 9. It can be observed that the IAOA algorithm achieved the optimal value during the second iteration. The optimized parameters obtained from the IAOA–VMD optimization search, namely, the mode component, K (12), and the penalty factor, α (9910), were then used in the VMD of the original acoustic signal. This decomposition step aims to effectively separate the signal components for further analysis and processing.
After substituting the optimal parameters into the VMD, 12 IMFs components were obtained. These IMF components were then used to construct a feature matrix with dimensions of 12 × 12,000. This feature matrix served as the original dictionary for the subsequent K-SVD.
By applying the K-SVD using the constructed feature matrix, the signal after enhancing the fault signal was obtained. Figure 10 illustrates the time–domain diagram of the enhanced signal. Comparing it with the original signal, it is evident that the enhanced signal successfully removed a substantial amount of shock noise while significantly improving the smoothness of the periodic signal.
To further investigate the amplitude relationship between the corresponding fault frequency and the rotation frequency in the enhanced signal, squared-envelope analysis was performed on the original signal, the VMD-processed signal, and the VMD–K-SVD-processed signal.
Figure 11a–c display the squared-envelope spectra of the original signal, the VMD-processed signal, and the VMD–K-SVD-processed signal, respectively.
From Figure 11a, only the rotation frequency of the rolling bearing and its octave frequency are visible, while the fault frequency of the bearing cannot be distinguished clearly. Figure 11b demonstrates that the amplitude of the fault signal was significantly improved compared with the original signal. However, the fault frequency is still not clearly discernible. In contrast, Figure 11c shows that the amplitude of the fault characteristic frequency is significantly enhanced compared with the previous spectra. It is now possible to clearly distinguish the rotation frequency of the rolling bearing and accurately determine its operational status. This result confirms the effectiveness of the proposed method for enhancing weak fault signals.
The squared-envelope spectrum of the original signal, VMD-processed signal, VMD–K-SVD-processed signal, EMD–K-SVD-processed signal, and CEEMDAN–K-SVD-processed signal are presented in Figure 11a, b, c, d, and e, respectively.
In Figure 11a, only the rotation frequency of the rolling bearing and its octave frequency are visible, while the fault frequency of the bearing cannot be distinguished. Figure 11b demonstrates a noticeable improvement in the amplitude of the fault signal; however, it is still not clearly evident. Conversely, Figure 11c shows a significant improvement in the amplitude of the fault characteristic frequency compared with the previous cases. The rotation frequency of the rolling bearing is clearly distinguishable, enabling the determination of its operational status. This substantiates the effectiveness of the proposed method for weak fault signals. Based on the comparison of Figure 11c–e, it can be observed that the amplitude of the fault frequency of the rolling bearing in the signals processed by EMD–K-SVD and CEEMDAN–K-SVD is enhanced compared with that of VMD. However, the effect is still weaker than that achieved by the VMD–K-SVD method proposed in this paper.

4. Conclusions

This paper addresses the challenge of separating noise in rolling bearing signals by proposing an effective method that combines VMD with K-SVD for the signal enhancement of bearing fault characteristics.
(1)
The optimized VMD algorithm based on the IAOA is capable of effectively separating the overall noise into individual IMF components, thereby improving the clarity of each IMF.
(2)
By subjecting the original dictionary, constructed using the IMF components, to K-SVD processing, the embedded noise is significantly reduced. This noise reduction step contributes to the enhancement of the fault characteristic signals, resulting in improved signal quality.
(3)
The effectiveness of the proposed method is demonstrated through two sets of experiments involving low-sampling-frequency vibration signals and high-sampling-frequency acoustic signals. These experiments validate the method’s ability to reduce noise and enhance fault characteristic signals, consequently enhancing the accuracy of fault diagnosis. However, it is important to note that the datasets used in this paper are experimental data and may not fully represent the high noise conditions encountered in complex machinery. Therefore, further testing and validation under real working conditions with highly noisy data are necessary to confirm the method’s effectiveness.
It is crucial to conduct additional research and experimentation to evaluate the proposed method’s performance in real-world scenarios with complex machinery and high-noise environments.

Author Contributions

Conceptualization, M.M., K.Z. and Z.T.; Formal analysis, Z.Z.; Funding acquisition, M.M., C.Q. and X.C.; Investigation, X.C.; Methodology, Z.T.; Software, Z.H.; Writing—original draft, K.Z. and M.M.; Writing—review and editing, Z.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of Hunan, grant number 2021JJ50125, 2022JJ50133 and 2022JJ50134.

Data Availability Statement

The CWRU data can be found at https://engineering.case.edu/bearingdatacenter/12k-drive-end-bearing-fault-data (accessed on 20 December 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yan, R.; Gao, R.X. Hilbert-Huang Transform-Based Vibration Signal Analysis for Machine Health Monitoring. IEEE Trans. Instrum. Meas. 2006, 55, 2320–2329. [Google Scholar] [CrossRef]
  2. Sohn, H.; Farrar, C.R. Damage Diagnosis Using Time Series Analysis of Vibration Signals. Smart Mater. Struct. 2001, 10, 446. [Google Scholar] [CrossRef]
  3. Salunkhe, V.G.; Desavale, R.G.; Khot, S.M.; Yelve, N.P. A Novel Incipient Fault Detection Technique for Roller Bearing Using Deep Independent Component Analysis and Variational Modal Decomposition. J. Tribol. 2023, 145, 074301. [Google Scholar] [CrossRef]
  4. Salunkhe, V.G.; Desavale, R.G.; Jagadeesha, T. A Numerical Model for Fault Diagnosis in Deep Groove Ball Bearing Using Dimension Theory. Mater. Today Proc. 2021, 47, 3077–3081. [Google Scholar] [CrossRef]
  5. Jakubek, B.; Grochalski, K.; Rukat, W.; Sokol, H. Thermovision Measurements of Rolling Bearings. Measurement 2022, 189, 110512. [Google Scholar] [CrossRef]
  6. Salunkhe, V.G.; Khot, S.M.; Desavale, R.; Yelve, N. Unbalance Bearing Fault Identification Using Highly Accurate Hilbert-Huang Transform Approach. J. Nondestruct. Eval. Diagn. Prognost. Eng. Syst. 2023, 1, 1–28. [Google Scholar] [CrossRef]
  7. Chen, X.; Zhang, B.; Gao, D. Bearing Fault Diagnosis Based on Multi-Scale CNN and LSTM Model. J. Intell. Manuf. 2021, 2, 971–987. [Google Scholar] [CrossRef]
  8. Peng, B.; Bi, Y.; Xue, B.; Zhang, M.; Wan, S. A Survey on Fault Diagnosis of Rolling Bearings. Algorithms 2022, 15, 347. [Google Scholar] [CrossRef]
  9. Peng, Z.; Chu, F.; He, Y. Vibration Signal Analysis and Feature Extraction Based on Reassigned Wavelet Scalogram. J. Sound Vib. 2002, 253, 1087–1100. [Google Scholar] [CrossRef]
  10. König, F.; Jacobs, G.; Stratmann, A.; Cornel, D. Fault Detection for Sliding Bearings Using Acoustic Emission Signals and Machine Learning Methods. In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2021; Volume 1097, p. 012013. [Google Scholar]
  11. Elforjani, M.; Shanbr, S. Prognosis of Bearing Acoustic Emission Signals Using Supervised Machine Learning. IEEE Trans. Ind. Electron. 2017, 65, 5864–5871. [Google Scholar] [CrossRef]
  12. Ganesan, S.; David, P.W.; Balachandran, P.K.; Samithas, D. Intelligent Starting Current-Based Fault Identification of an Induction Motor Operating under Various Power Quality Issues. Energies 2021, 14, 304. [Google Scholar] [CrossRef]
  13. Dragomiretskiy, K.; Zosso, D. Variational Mode Decomposition. IEEE Trans. Signal Process. 2013, 62, 531–544. [Google Scholar] [CrossRef]
  14. Boudraa, A.O.; Cexus, J.C. EMD-Based Signal Filtering. IEEE Trans. Instrum. Meas. 2007, 56, 2196–2202. [Google Scholar] [CrossRef]
  15. Yeh, J.R.; Shieh, J.S.; Huang, N.E. Complementary Ensemble Empirical Mode Decomposition: A Novel Noise Enhanced Data Analysis Method. Adv. Adapt. Data Anal. 2010, 2, 135–156. [Google Scholar] [CrossRef]
  16. Torres, M.E.; Colominas, M.A.; Schlotthauer, G.; Flandrin, P. A Complete Ensemble Empirical Mode Decomposition with Adaptive Noise. In Proceedings of the 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Prague, Czech Republic, 22–27 May 2011; IEEE: Piscataway Township, NJ, USA, 2011; pp. 4144–4147. [Google Scholar]
  17. Jin, Z.; Chen, G.; Yang, Z. Rolling Bearing Fault Diagnosis Based on WOA-VMD-MPE and MPSO-LSSVM. Entropy 2022, 24, 927. [Google Scholar] [CrossRef]
  18. Zhang, M.; Jiang, Z.; Feng, K. Research on Variational Mode Decomposition in Rolling Bearings Fault Diagnosis of the Multistage Centrifugal Pump. Mech. Syst. Signal Process. 2017, 93, 460–493. [Google Scholar] [CrossRef]
  19. Zhang, X.; Miao, Q.; Zhang, H.; Wang, L. A Parameter-Adaptive VMD Method Based on Grasshopper Optimization Algorithm to Analyze Vibration Signals from Rotating Machinery. Mech. Syst. Signal Process. 2018, 108, 58–72. [Google Scholar] [CrossRef]
  20. Aharon, M.; Elad, M.; Bruckstein, A. K-SVD: An Algorithm for Designing Overcomplete Dictionaries for Sparse Representation. IEEE Trans. Signal Process. 2006, 54, 4311–4322. [Google Scholar] [CrossRef]
  21. Lebrun, M.; Leclaire, A. An Implementation and Detailed Analysis of the K-SVD Image Denoising Algorithm. Image Process. OnLine 2012, 2, 96–133. [Google Scholar] [CrossRef]
  22. Anaraki, F.P.; Hughes, S.M. Compressive K-SVD. In Proceedings of the 2013 IEEE International Conference on Acoustics, Speech and Signal Processing, Vancouver, BC, Canada, 26–31 May 2013; IEEE: Piscataway Township, NJ, USA, 2013; pp. 5469–5473. [Google Scholar]
  23. Wu, Y.; Lou, L.; Wang, J. Cotton Fabric Defect Detection Based on K-SVD Dictionary Learning. J. Nat. Fibers 2022, 19, 10764–10779. [Google Scholar] [CrossRef]
  24. Wang, L.; Li, X.; Xu, D.; Ai, S.; Wang, C.; Chen, C. Bearing Fault Feature Extraction Based on Adaptive OMP and Improved K-SVD. Processes 2022, 10, 675. [Google Scholar] [CrossRef]
  25. Abualigah, L.; Diabat, A.; Mirjalili, S.; Elaziz, M.A.; Gandomi, A.H. The Arithmetic Optimization Algorithm. Comput. Methods Appl. Mech. Eng. 2021, 376, 113609. [Google Scholar] [CrossRef]
  26. Chai, X.; Gan, Z.; Chen, Y.; Zhang, Y. A Visually Secure Image Encryption Scheme Based on Compressive Sensing. Signal Process. 2017, 134, 35–51. [Google Scholar] [CrossRef]
  27. Hua, Z.; Zhou, Y.; Huang, H. Cosine-Transform-Based Chaotic System for Image Encryption. Inf. Sci. 2019, 480, 403–419. [Google Scholar] [CrossRef]
  28. Xu, Y.; Yang, Z.; Li, X.; Kang, H.; Yang, X. Dynamic Opposite Learning Enhanced Teaching-Learning-Based Optimization. Knowl.-Based Syst. 2020, 188, 104966. [Google Scholar] [CrossRef]
  29. Çelik, E. IEGQO-AOA: Information-Exchanged Gaussian Arithmetic Optimization Algorithm with Quasi-Opposition Learning. Knowl.-Based Syst. 2023, 260, 110169. [Google Scholar] [CrossRef]
  30. Smith, W.A.; Randall, R.B. Rolling Element Bearing Diagnostics Using the Case Western Reserve University Data: A Benchmark Study. Mech. Syst. Signal Process. 2015, 64, 100–131. [Google Scholar] [CrossRef]
Figure 1. The IAOA optimization flowchart.
Figure 1. The IAOA optimization flowchart.
Sensors 23 08629 g001
Figure 2. Flowchart depicting the fault characteristic signal enhancement process.
Figure 2. Flowchart depicting the fault characteristic signal enhancement process.
Sensors 23 08629 g002
Figure 3. Time–domain diagram of the original fault signal.
Figure 3. Time–domain diagram of the original fault signal.
Sensors 23 08629 g003
Figure 4. Iteration of the fitness function.
Figure 4. Iteration of the fitness function.
Sensors 23 08629 g004
Figure 5. Enhanced signal time–domain diagram.
Figure 5. Enhanced signal time–domain diagram.
Sensors 23 08629 g005
Figure 6. Squared-envelope spectrum analysis.
Figure 6. Squared-envelope spectrum analysis.
Sensors 23 08629 g006aSensors 23 08629 g006b
Figure 7. Experimental setup.
Figure 7. Experimental setup.
Sensors 23 08629 g007
Figure 8. Time–domain diagram of the original signal.
Figure 8. Time–domain diagram of the original signal.
Sensors 23 08629 g008
Figure 9. VMD iteration diagram.
Figure 9. VMD iteration diagram.
Sensors 23 08629 g009
Figure 10. Enhanced signal time–domain diagram.
Figure 10. Enhanced signal time–domain diagram.
Sensors 23 08629 g010
Figure 11. Squared-envelope spectrum analysis.
Figure 11. Squared-envelope spectrum analysis.
Sensors 23 08629 g011aSensors 23 08629 g011b
Table 1. Parameters related to rolling bearings to be tested.
Table 1. Parameters related to rolling bearings to be tested.
ParametersBearing TypeInner Ring Diameter/mmOuter Ring Diameter/mmThickness/mmBearing Contact Angle/°Pitch Diameter/mmNumber of Scrollers
Numerical value29412E6013042309516
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

Mao, M.; Zeng, K.; Tan, Z.; Zeng, Z.; Hu, Z.; Chen, X.; Qin, C. Adaptive VMD–K-SVD-Based Rolling Bearing Fault Signal Enhancement Study. Sensors 2023, 23, 8629. https://doi.org/10.3390/s23208629

AMA Style

Mao M, Zeng K, Tan Z, Zeng Z, Hu Z, Chen X, Qin C. Adaptive VMD–K-SVD-Based Rolling Bearing Fault Signal Enhancement Study. Sensors. 2023; 23(20):8629. https://doi.org/10.3390/s23208629

Chicago/Turabian Style

Mao, Meijiao, Kaixin Zeng, Zhifei Tan, Zhi Zeng, Zihua Hu, Xiaogao Chen, and Changjiang Qin. 2023. "Adaptive VMD–K-SVD-Based Rolling Bearing Fault Signal Enhancement Study" Sensors 23, no. 20: 8629. https://doi.org/10.3390/s23208629

APA Style

Mao, M., Zeng, K., Tan, Z., Zeng, Z., Hu, Z., Chen, X., & Qin, C. (2023). Adaptive VMD–K-SVD-Based Rolling Bearing Fault Signal Enhancement Study. Sensors, 23(20), 8629. https://doi.org/10.3390/s23208629

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