Published: February 28, 2026

Fault diagnosis of rolling bearings based on amplitude modulation of local W transform spectrum

Zongcai Ma1
Yongqi Chen2
Qinge Dai3
Linqiang Wu4
Yitong Qin5
1, 2, 3, 5College of Science and Technology Ningbo University, Ningbo, 315211, China
4Ningbo Donghuang Bearing Co., Ltd., Ningbo, 315211, China
Corresponding Author:
Yongqi Chen
Article in Press
Views 0
Reads 0
Downloads 0

Abstract

As one of the key components of transmission and support in rotating machinery, rolling bearings will directly affect the operating status of the equipment. Therefore, scholars have proposed various methods to process its fault signals. Among them, spectral amplitude modulation (SAM) is a nonlinear filtering method that can effectively extract bearing fault characteristics. However, due to the high-intensity noise interference in the working environment of the equipment, the interference component dominates the monitoring signal, and the fault characteristics are no longer obvious or even undetectable. To solve the above problems, this paper proposes a local W transform spectral amplitude (LWTSAM) modulation method. First, the method divides the original signal into multiple narrowband signals in different frequency bands and selects the narrowband signal with the most abundant fault information among them; then, the selected narrowband signal uses a windowed Fourier transform (WFT) to obtain the amplitude in the time-frequency domain, and uses different weight indices (MO) to correct the amplitude; Finally, the modified amplitude is combined with its original phase to perform inverse window Fourier transform to obtain the modified signal and its square envelope, and the square envelope under the optimal weight is calculated using unbiased autocorrelation and information entropy to complete the local W transform spectrum amplitude modulation. In this paper, this method is verified through fault data sets. The research results show that this method can effectively reduce the interference of noise on fault diagnosis, and the fault characteristic information obtained is clearer. Compared with SAM method, Autogram method and fast spectral kurtosis diagram method, the superiority of this method is proved.

Fault diagnosis of rolling bearings based on amplitude modulation of local W transform spectrum

Highlights

  • The method divides the original signal into multiple narrowband signals in different frequency bands and selects the narrowband signal with the most abundant fault information among them.
  • The selected narrowband signal uses a windowed Fourier transform (WFT) to obtain the amplitude in the time-frequency domain, and uses different weight indices (MO) to correct the amplitude.
  • The modified amplitude and original phase are reconstructed via inverse WFT, and the optimal weight is selected using autocorrelation and entropy to achieve LWTSAM.

1. Introduction

As one of the key components in rotating machinery, rolling bearings need to withstand various dynamic loads and complex working conditions, and internal components are prone to failure [1, 2]. Once a failure occurs, major safety accidents will be caused, resulting in mechanical system failures and economic losses [3, 4], and even casualties. Therefore, bearing failures should be discovered as soon as possible to avoid accidents. However, due to the complex and diverse working environment, the collected vibration signals include a large amount of environmental noise and irrelevant pulses, making faults difficult to identify at the early stage. Therefore, the importance of identifying weak bearings under complex conditions is self-evident [5, 6].

In recent decades, bearing fault diagnosis has gradually developed into an important research field in academia and industry. Among the currently known fault diagnosis technologies, resonance demodulation technology [7] is undoubtedly one of the most classic methods, also known as envelope analysis. The corresponding fault information can be obtained by band-pass filtering the resonance band with high signal-to-noise ratio and then demodulating it using the Hilbert transform [8, 9]. However, the selection of the resonance band is a relatively empirical process and is difficult in practical engineering applications. In order to overcome this problem, Antoni and Randall [10-11] conducted a detailed and meaningful study of spectral kurtosis and proposed a more practical and faster method, Kurtogram [12]. In Kurtogram, a 1/3 binary tree filter bank is used to partition the spectrum and determine the best signal based on kurtosis. However, kurtosis is easily disturbed by noise and random pulses [13]. In order to improve its practical application capabilities, scholars have proposed a large number of improvement methods, such as Protrugram [14], Infogram [15], enhanced Kurtogram [16], improved Kurtogram [17], Autogram [18]. All of the above methods have achieved good results. In addition, other kurtosis methods based on period indicators have been proposed [19-23]. These methods can specify the frequency band in which the periodic component is located, and all suppress noise and irrelevant pulses to a certain extent. However, the biggest drawback of this kind of method is that only one defect can be detected and other defects are ignored. When multiple defects appear in different frequency bands, fault information will be lost and fault information will be lacking, which brings great difficulties to bearing fault diagnosis.

As a common method to separate signal periods and random components, cepstrum pre-whitening is fast and easy to implement and can be used for automatic diagnosis [24, 25]. However, setting the amplitudes of all frequency components to unit amplitude seriously reduces the SNR and greatly limits practical applications. Therefore, Moshrefzadeh et al. introduced spectral amplitude modulation (SAM), a new nonlinear filtering technique [26]. The method performs Fourier transformation on the obtained vibration signal to obtain its amplitude spectrum and phase spectrum and then assigns different weight indices (MO) to its phase spectrum. The resulting edit spectrum is subjected to inverse Fourier transformation to obtain a series of corrected signals. Finally, the square spectrum of the corrected signal is Fourier transformed to obtain a square envelope spectrum to view fault characteristic information. At present, SAM is favored by many scholars for its powerful nonlinear enhancement capabilities. Zhu et al.[27] combined SAM with Savitzky-Golay filtering to suppress noise interference and signal distortion during signal transmission; Moshrefzadeh [28] combined it with support vector machines to expand the application of SAM in online monitoring of variable operating conditions; Liu [29] applied unbiased autocorrelation to corrected signals to optimize SAM results; Wang [30] used SAM to extract fault characteristics under variable speed conditions and test its noise suppression ability; Zhong [31-32] used SAM to highlight and reconstruct the pulse characteristics in the signal; Ma [33] pointed out that at certain amplitude orders, fault information may be concealed by noise, and gave a generalized envelope spectrum to improve fault characteristics, which was verified on a locomotive bearing test stand.

Although the above studies have achieved good research results, SAM still faces huge challenges in practical applications. First of all, even if SAM can reduce the impact of interference components on fault components, if it is in a strong interference environment, the interference component is still the dominant component, which will cause fault information to be submerged in a large amount of interference; Secondly, because SAM processes signals, Scanning the entire frequency band amplifies fault information while also enhancing interference information, which means that the final result may be erased due to interference. Based on the above problems, a new local W transform spectral amplitude modulation (LWTSAM) method is proposed. In this paper, first, the maximum discrete overlapped wavelet transform (MODWT) is used to decompose the original signal into multiple narrowband signals in different frequency bands to remove some interference components; then the windowed Fourier transform (WFT) is used to calculate the phase and amplitude of the narrowband signal with the richest fault information and select the square envelope under the optimal weight, which will get more accurate results, because the amplitude is not calculated over the entire frequency band. Experimental results show that this method has good results in bearing fault diagnosis and analysis.

To sum up, a review table of the comparative literature is given, as shown in Table 1.

Table 1Comparison between existing bearing fault diagnosis methods and the proposed LWTSAM

Method
Core idea
Noise robustness
Multi-fault capability
Automatic band selection
Limitation
Kurtogram [10-12]
Kurtosis-based band selection
Medium
No
No
Sensitive to impulsive noise
Autogram [18]
Thresholded kurtosis
Medium
No
Yes
Noise-dominated envelope
SAM [26]
Global spectral amplitude modulation
Medium
No
No
Global enhancement amplifies noise
LWTSAM
MODWT + local WFT-SAM
High
Yes
Yes
Slightly increased complexity

2. SAM theory

The SAM method was proposed by Moshrefzadeh [26]. As an advanced nonlinear filtering technique, it extracts frequency components of different energies in the signal by setting different MOs without requiring the input of any parameters in advance, thereby realizing bearing fault diagnosis., made great contributions in the field of bearing fault diagnosis.

First, FT is used to convert the original signal from the time domain to the frequency domain:

1
Xf=FTx(t)=AfejΦf,

where j is the imaginary unit, FT is the Fourier transform, and Af and Φf represent the amplitude and phase of the spectrum respectively.

Subsequently, the phase remains unchanged, and different weight indices (MO) are assigned to the amplitude spectrum to adjust the amplitude in different frequency bands. Then the original phase is combined with the modulated amplitude to obtain the edited spectrum. Finally, inverse Fourier transform is performed to obtain the correction signal. This nonlinear filtering process is as follows:

2
xmt,MO=IFTA(f)MOejΦ(f),

where xm is the correction signal under different weight indices, and IFT is the inverse Fourier transform. The recommended value for MO is -0.5MO1.5, with a step size of 0.1. The advantage of this is that it can adjust the size of MO to diagnose faults under different situations and effectively shield interference components. Among them, there are three main effects of MO value on the signal: when MO<0, if the noise component does not dominate, the fault component with very low energy level can be extracted; when 0<MO<1, the dominant position of high amplitude frequency gradually decreases, but the high-energy components and low-energy components will be retained to a certain extent; when MO>1, the larger amplitude component is further amplified to conceal the smaller amplitude component; when MO=1, it is the original signal. The signal spectra under different MO are shown in Fig. 1.

Then, the SES for each modified signal is calculated as follows:

3
SESxmt,MO=FTxmt,MO+j*Hilbertxmt,MO2.

Next, the SES are normalized between [0, 1] for explicit comparisons:

4
NSESxmt,MO=SESxmt,MOmax(SESxmt,MO).

Finally, use a three-dimensional diagram to display the results. Where the x-axis represents the modulation frequency, the y-axis represents MO, and the z-axis represents the normalized amplitude. In addition, other visual view demonstration results will also be given. One view is a top view of the three-dimensional map, which is a two-dimensional color map along the normalized amplitude axis, and the color intensity corresponds to the normalized amplitude magnitude. This view is clear for analyzing the relationship between the intensity of each modulation frequency and the corresponding MO value. The other view is a front view of the three-dimensional graph, which is the logarithmic maximum square envelope spectrum plotted along the MO axis, which provides the modulation frequency and corresponding amplitude for different MO values. As shown in Fig. 2.

Fig. 1Signal spectrum under different MO

Signal spectrum under different MO
Signal spectrum under different MO
Signal spectrum under different MO
Signal spectrum under different MO

SAM restores the original signal by assigning different amplitudes and different weight indices, and combines it with the original phase, and determines fault components according to different angles. This has proven to be a convenient, efficient and automated method. However, in complex industrial environments, the application of SAM will face some challenges. First, MO is assigned to the entire amplitude spectrum, but if the fault information only appears in a certain narrow band at this time, this will mean that the weight of the fault information will increase and the noise information will also be enhanced. Then after SAM processing, the fault information will be covered up by a large amount of enhanced noise; In addition, since the Fourier transform is performed over the entire time domain, which means that the amplitude of each frequency in the amplitude spectrum is based on the contribution of the entire time process, inaccuracies will undoubtedly be amplified in subsequent calculations.

To address the above problems, it is necessary to explore a more ideal method.

Fig. 2SAM example

SAM example
SAM example
SAM example

3. Recommended method: local W-transform spectral amplitude modulation

A method called LWTSAM was proposed to solve the above problems. First, the entire spectrum is segmented using the Maximum Discrete Overlapped Wavelet Transform (MODWT), separated into narrowband signals in different frequency bands, and selected the narrowband with the richest fault information; secondly, the narrowband signal is used to obtain its phase and amplitude and reconstruct it. The reconstructed signal can be used to extract the characteristic frequency more obviously. Based on the above statement, the LWTSAM flow chart is shown in Fig. 3, and the following section details the implementation and explanation of each step.

3.1. Signal segmentation via maximum discrete overlapping wavelet transform

As mentioned earlier, SAM can enhance the fault frequency component by modifying the weighting index. However, for a strong interference environment, the interference component and the fault component may have the same amplitude. In this case, the fault component and the noise component will be enhanced at the same time, and the SAM diagnosis will be invalid. Bearing fault components usually exist in the fault resonance frequency band, so the original signal can be divided into several narrowband signals to maximize the separation of the fault components from most of the noise components.

Fig. 3LWTSAM method flowchart

LWTSAM method flowchart

The Maximum Discrete Overlap Wavelet Transform [34] can decompose a signal into multiple frequency bands, and it is convolved through maximum overlap, which means that at each scale, the convolution window of the wavelet function and the scale function covers the entire signal, thereby retaining more information. This section uses the maximum discrete overlapping wavelet transform to segment the original signal into multiple narrowband signals. First, the characteristics of the wavelet transform are defined by selecting appropriate wavelet basis functions ψ(t) and scale functions ϕ(t); secondly, at each scale j, the signal Xf is decomposed into approximation coefficients cj,k and detail coefficients dj,k:

5
cj,k=fXf*ϕj,kf,dj,k=fXf*ψj,kf,

where k is the translation parameter, ϕj,kf, ψj,kf are defined as follows:

6
ϕj,kf=2-j2ϕ2-jf-k,ψj,kf=2-j2ψ2-jf-k.

Finally, the original signal is reconstructed from the MODWT coefficients to obtain the reconstructed signal Xrf:

7
Xrf=kcJ,k*ϕJ,kf+j=1Jkdj,k*ψj,kf.

The reconstructed signal can be decomposed into multiple narrowband signals through the filter bank included in MODWT:

8
Xrn=jkdjkg~j,kn+jkωjkh~j,kn,

where n is the time index, djk is the scale coefficient of the j layer, ωjk is the wavelet coefficient of the j layer, g~j,kn and h~j,kn are the inverse scale filter and the inverse wavelet filter respectively.

3.2. Signal segmentation via maximum discrete overlapping wavelet transform

Although the optimal narrowband signal obtained by MODWT segmentation has filtered out most irrelevant components, there are inevitably some interference components in a single narrowband. Therefore, in order to observe and extract more comprehensive information, a windowed Fourier transform is used for the optimal narrowband signal. For the selected optimal narrowband signal x(t) and window function ω(t), it is defined as:

9
Xωω,τ=-xtωt-τe-jωtdt=At,ωejϕt,ω,

where Xωω,τ is the result of the windowed Fourier transform, which is a function of time shifted τ and frequency ω; ωt-τ is a window function shifted by time; e-jωt is a complex exponential function, and j is an imaginary unit. Then, the inverse window Fourier transform is calculated by extracting the IFT signal of each modified FT vector in a certain window of the WFT and overlapping and adding them. However, in practical applications, this process can be achieved by providing a different MO value for each amplitude value in the time-frequency matrix. It can be expressed as:

10
xmLWTτ,MO=12π--At,ωMOejωt+ϕt,ωωt-τdtdω,

where xmLWTτ,MO represents the correction signal obtained by the proposed method.

Window Fourier transform can analyze local frequency components of a signal by sliding a window on the signal, which is particularly useful for non-stationary signals; Fourier transform provides global frequency information of the entire signal and cannot distinguish between frequency changes in different time periods, and the amplitude obtained by Fourier transform is the average of the amplitudes obtained over the entire time, which can lead to significant differences in results. A comparison of these two methods is shown in Fig. 4. For the same MO value, it can be observed that LWTSAM reveals the resonance band well, while the resonance band in SAM is submerged by noise, which also leads to completely different diagnostic results for the envelope spectrum.

Fig. 4a) SAM modified signal; b) LWTSAM modified signal; c) SAM modified signal spectrum; d) LWTSAM modified signal spectrum; e) SAM modified signal envelope spectrum; f) LWTSAM modified signal envelope spectrum

a) SAM modified signal; b) LWTSAM modified signal; c) SAM modified signal spectrum;  d) LWTSAM modified signal spectrum; e) SAM modified signal envelope spectrum;  f) LWTSAM modified signal envelope spectrum

a)

a) SAM modified signal; b) LWTSAM modified signal; c) SAM modified signal spectrum;  d) LWTSAM modified signal spectrum; e) SAM modified signal envelope spectrum;  f) LWTSAM modified signal envelope spectrum

b)

a) SAM modified signal; b) LWTSAM modified signal; c) SAM modified signal spectrum;  d) LWTSAM modified signal spectrum; e) SAM modified signal envelope spectrum;  f) LWTSAM modified signal envelope spectrum

c)

a) SAM modified signal; b) LWTSAM modified signal; c) SAM modified signal spectrum;  d) LWTSAM modified signal spectrum; e) SAM modified signal envelope spectrum;  f) LWTSAM modified signal envelope spectrum

d)

a) SAM modified signal; b) LWTSAM modified signal; c) SAM modified signal spectrum;  d) LWTSAM modified signal spectrum; e) SAM modified signal envelope spectrum;  f) LWTSAM modified signal envelope spectrum

e)

a) SAM modified signal; b) LWTSAM modified signal; c) SAM modified signal spectrum;  d) LWTSAM modified signal spectrum; e) SAM modified signal envelope spectrum;  f) LWTSAM modified signal envelope spectrum

f)

Therefore, in this study, WFT was used to improve spectral amplitude modulation. The amplitude calculated by WFT no longer corresponds to the entire time domain, so a more accurate corrected amplitude can be obtained in subsequent calculations. As shown in Fig. 5, it can be observed that by assigning different weight indices to the time-frequency domain, fault characteristics can also be highlighted for certain specific MO values, which verifies the feasibility of the method in this paper.

3.3. Select the best weights through unbiased autocorrelation and information entropy

As mentioned earlier, when a signal is subjected to high-intensity interference, a large amount of noise will flood into the envelope spectrum, which can easily drown out the fault characteristic frequency. At this point, it is undoubtedly a huge challenge to manually analyze the harmonic bands and correctly select the best MO. Therefore, an indicator that can automatically select the best MO is needed to quantify the fault information in the envelope spectrum. Although harmonic spectrum kurtosis [35], spectral negative entropy [15] and sparsity [36] are excellent indicators for detecting pulse components in bearing vibration signals, in order to highlight signal characteristics and strengthen feature information, this paper adopts an indicator based on unbiased autocorrelation (AC) and information entropy adaptively selecting the optimal MO [37].

Fig. 5a) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= –0.5; b) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 0.5; c) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1; d) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1.5

a) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= –0.5;  b) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 0.5;  c) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1;  d) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1.5

a)

a) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= –0.5;  b) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 0.5;  c) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1;  d) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1.5

b)

a) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= –0.5;  b) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 0.5;  c) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1;  d) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1.5

c)

a) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= –0.5;  b) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 0.5;  c) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1;  d) Time-frequency domain spectrum, modified signal and envelope spectrum at MO= 1.5

d)

First, the square envelope of the modified signal is calculated through the Hilbert transform:

11
Ste,MO=xmLWTte,MO+HilbertxmLWTte,MO2,

where S is the square envelope of the modified signal, e is the signal length, and the interval is [0,l].

Unbiased autocorrelation is then used to enhance the periodic component of the signal and weaken the random pulse component:

12
Rxτ,MO=1l-ek=1l-eStk,MOStk+τ,MO,

where τ is the delay coefficient, calculated as follows:

13
τ=efs,

where, fs represents the sampling frequency. In Eq. (12), the number of data samples used to calculate AC will decrease as τ increases, which may lead to insufficient estimation variance in the result part. Therefore, only the first half of the AC is selected, and the uncertainty of the result is represented by information entropy, which is defined as follows:

14
IMO=-i=1l/2pRxi,MOlog2pRxi,MO,

where, IMO corresponds to the result of each suggested indicator, and p is the probability distribution:

15
p=hAl,

where, h represents the amplitude under a certain range of sample numbers A. Since random pulse components in the vibration signal will cause energy fluctuations, and periodic fault components will cause entropy reduction, MO with the smallest index value is used as the best weight. The results are shown in Fig. 6. It can be seen that the fault component dominates and the interference component can be ignored.

Fig. 6a) Indicator analysis results; b) 2D diagram of LWTSAM; c) envelope spectrum at MO= 0.5

a) Indicator analysis results; b) 2D diagram of LWTSAM; c) envelope spectrum at MO= 0.5

a)

a) Indicator analysis results; b) 2D diagram of LWTSAM; c) envelope spectrum at MO= 0.5

b)

a) Indicator analysis results; b) 2D diagram of LWTSAM; c) envelope spectrum at MO= 0.5

c)

4. Experimental verification

This section uses real data from the Southeast University Bearing Failure Test Bench [38] to verify the feasibility of LWTSAM. It will be compared with SAM method, Autogram method and fast spectral kurtosis diagram method.

4.1. Experimental device

The vibration signal was collected at the motor speed of 1800 r/min and the sampling frequency of 12 kHz. The experiment will use outer ring, inner ring and composite fault data sets of inner and outer rings to verify this research. The fault types are shown in Table 2, and samples under a load of 30 Hz-2 V are used.

Table 2Description of SEU bearing faults

Speed-load
Fault type
Description
Label
30 Hz-2 V
Outer fault
Outer crack
1
Inner fault
Inner crack
2
Comb fault
Comb crack
3

4.2. Fault diagnosis of bearing outer ring

This section will provide experimental verification of bearing with outer ring failure. The actual running speed of the motor is 1768.6 r/min, and the actual rotating frequency is 29.91 Hz. After calculation, the characteristic frequency fo = 93.75 Hz can be obtained. Fig. 7 shows the time domain diagram and envelope spectrum of the original signal. Due to the interference of strong noise, there is no particularly obvious periodic effect in the time domain; the fault frequency and its doubling cannot be observed in the envelope spectrum. It is not easy to find more useful information.

First, the SAM method is used to process the signal, as shown in Fig. 8. It can be observed that no obvious fault components are observed in either the 3D image or the 2D image. fo in log-MSES is concealed by a large amount of noise, and only weak 2fo is observed. The interference of invalid components makes fault feature extraction difficult.

The LWTSAM method was then used to process the signal. First, the signal is divided using the filter bank that comes with MODWT. The number of filters is set to 8 and the bandwidth is set to 600 Hz. Random impact interference can be successfully separated. According to the parameters, narrowband signal #2 contains fault information, and all narrowband signals are shown in Fig. 9.

Fig. 7Time domain diagram and envelope spectrum of outer ring original signal

Time domain diagram and envelope spectrum of outer ring original signal
Time domain diagram and envelope spectrum of outer ring original signal

Fig. 8SAM outer ring fault diagnosis

SAM outer ring fault diagnosis
SAM outer ring fault diagnosis
SAM outer ring fault diagnosis

Narrow-band signal #2 and its envelope spectrum are shown in Fig. 10, and fault-related information is hidden in harmonic interference.

Although the fault frequency and its doubling can be observed in the segmented narrowband signal, it still appears weak compared with the interference component. So, we proceed to the next step. WFT is used to reconstruct narrowband signal #2. The reconstructed signal and its envelope spectrum are shown in Fig. 11. It can be seen that the noise is greatly reduced, and three harmonics of the characteristic frequency can be found in log-MSES, which shows that the LWTSAM method has obvious advantages over the SAM method. The recommended indicators are then used to select the best weights for extracting individual SES, and the minimum value of the recommended indicators appears at MO= 0.4, which is considered the best weights. As shown in Fig. 12, the characteristic frequency and its doubling are further enhanced, enough to determine that it is an outer ring fault.

Fig. 9Outer band signal

Outer band signal
Outer band signal
Outer band signal
Outer band signal
Outer band signal
Outer band signal
Outer band signal
Outer band signal
Outer band signal

Fig. 10Outer band signal #2 and its envelope spectrum

Outer band signal #2 and its envelope spectrum
Outer band signal #2 and its envelope spectrum

Fig. 11Outer ring reconstructed signal and its envelope spectrum

Outer ring reconstructed signal and its envelope spectrum
Outer ring reconstructed signal and its envelope spectrum

Fig. 12Outer ring reconstructed signal and its envelope spectrum

Outer ring reconstructed signal and its envelope spectrum
Outer ring reconstructed signal and its envelope spectrum
Outer ring reconstructed signal and its envelope spectrum
Outer ring reconstructed signal and its envelope spectrum

For further comparison, the fast spectral kurtosis method and the Autogram method continue to be used to process signals. The diagnosis results of the fast spectral kurtosis method are shown in Fig. 13. Where Kmax is the maximum kurtosis value, Bw is the bandwidth, and fc is the center frequency. From the kurtosis diagram, we can know that Kmax is 8.9, Bw is 500 Hz, and fc is 5750 Hz. The best frequency band is found in Level 3.5, and no useful fault information is observed in the envelope spectrum.

Fig. 13Fast spectral kurtosis method outer ring diagnostic results

Fast spectral kurtosis method outer ring diagnostic results
Fast spectral kurtosis method outer ring diagnostic results

Fig. 14Autogram method outer ring diagnostic results

Autogram method outer ring diagnostic results
Autogram method outer ring diagnostic results

The diagnosis results of the Autogram method are shown in Fig. 14. Where Kmax is the maximum kurtosis value, Bw is the bandwidth, and fc is the center frequency. From the kurtosis diagram, we can know that Kmax is 31.8, Bw is 750 Hz, and fc is 5625 Hz. The best frequency band is found at Level 3. The envelope spectrum results under the best frequency band are given. For further verification, the envelope spectra after the rising and falling thresholds are given respectively. It can be observed that although the fault frequency can be observed after processing, a large number of noise components dominate.

To sum up, the LWTSAM method is significantly better than the SAM method in processing outer ring fault signals. The fast spectral kurtosis diagram method and Autogram method verify the feasibility of the method in this paper.

4.3. Fault diagnosis of bearing inner ring

This section will conduct experimental verification of bearing with inner ring failure. The actual running speed of the motor is 1779.2 r/min, and the actual rotating frequency is 29.93 Hz. After calculation, the characteristic frequency of bearing inner ring fault fi = 107.61 Hz can be obtained. Fig. 15 shows the time domain diagram and envelope spectrum of the original signal. Due to the interference of a large amount of noise, there is no particularly obvious fault information in the time domain; the fault frequency and its doubling cannot be observed in the envelope spectrum, so it is difficult to extract more useful information.

Fig. 15Time domain diagram and envelope spectrum of inner ring original signal

Time domain diagram and envelope spectrum of inner ring original signal
Time domain diagram and envelope spectrum of inner ring original signal

Fig. 16SAM iuter ring fault diagnosis

SAM iuter ring fault diagnosis
SAM iuter ring fault diagnosis
SAM iuter ring fault diagnosis

First, the signal is processed using the SAM method, as shown in Fig. 16. Fault information can be clearly observed in both the 3D and 2D graphs. fi and 2fi can also be observed in log-MSES. However, there are some interference components with amplitude equivalent to the fault frequency, so it is very necessary to eliminate these interference.

The LWTSAM method was then used to process the signal. First, the signal is divided using the filter bank that comes with MODWT. The number of filters is set to 9 and the bandwidth is set to 600Hz. Random impact interference can be successfully separated. According to the parameters, narrowband signal #5 contains fault information, and all narrowband signals are shown in Fig. 17.

Fig. 17Inner band signal

Inner band signal
Inner band signal
Inner band signal
Inner band signal
Inner band signal
Inner band signal
Inner band signal
Inner band signal
Inner band signal

Narrow-band signal #5 and its envelope spectrum are shown in Fig. 18, and fault-related information is hidden in harmonic interference.

Obvious fi and 2fi can be observed from the obtained narrowband signal #5, but it can also be seen that there are many irrelevant components with amplitudes similar to the fault component, which can easily lead to misdiagnosis. So, we proceed to the next step. WFT is used to reconstruct narrowband signal #5. The reconstructed signal and its envelope spectrum are shown in Fig. 19. It can be seen that the noise is greatly weakened, but the fault component is also weakened. The characteristic frequency can be found in log-MSES. 2 harmonics, but the amplitude intensity is far less than that of the irrelevant component. The recommended indicators are then used to select the best weights for extracting individual SES, and the minimum value of the recommended indicators appears at MO= 0.1, which is considered the best weights. As shown in Fig. 20, the irrelevant components are eliminated and the fault components are enhanced, which is enough to determine that it is an inner ring fault.

Fig. 18Inner circle narrowband signal #5 and its envelope spectrum

Inner circle narrowband signal #5 and its envelope spectrum
Inner circle narrowband signal #5 and its envelope spectrum

Fig. 19Inner ring reconstructed signal and its envelope spectrum

Inner ring reconstructed signal and its envelope spectrum
Inner ring reconstructed signal and its envelope spectrum

Next, we continue to use the fast spectral kurtosis method and the Autogram method to process the signal. Fig. 21 shows the fast spectral kurtosis method. Where Kmax is the maximum kurtosis value, Bw is the bandwidth, and fc is the center frequency. From the kurtosis diagram, we can know that Kmax is 8.9, Bw is 500 Hz, and fc is 5750 Hz. The best frequency band is found in Level 3.5, and no useful fault information is observed in the envelope spectrum.

Fig. 20Final diagnostic results of LWTSAM method inner circle

Final diagnostic results of LWTSAM method inner circle
Final diagnostic results of LWTSAM method inner circle
Final diagnostic results of LWTSAM method inner circle
Final diagnostic results of LWTSAM method inner circle

The diagnosis results of Autogram method are shown in Fig. 22. Where Kmax is the maximum kurtosis value, Bw is the bandwidth, and fc is the center frequency. From the kurtosis diagram, we can know that Kmax is 11.2, Bw is 375 Hz, and fc is 5812.5 Hz. The best frequency band is found at Level 4. The envelope spectrum results under the best frequency band are given. For further verification, the envelope spectra after the rising and falling thresholds are given respectively. It can be observed that fault information is submerged by a large number of noise components.

To sum up, the LWTSAM method has also achieved good results in processing inner ring fault signals, which is significantly superior to other classic methods, which verifies the feasibility and superiority of the method in this paper.

Fig. 21Inner circle diagnosis results of fast spectral kurtosis method

Inner circle diagnosis results of fast spectral kurtosis method
Inner circle diagnosis results of fast spectral kurtosis method

Fig. 22Autogram method inner circle diagnosis results

Autogram method inner circle diagnosis results
Autogram method inner circle diagnosis results

Fig. 23Composite original signal time domain diagram and envelope spectrum

Composite original signal time domain diagram and envelope spectrum
Composite original signal time domain diagram and envelope spectrum

4.4. Bearing composite fault diagnosis

This section will conduct experimental verification of composite fault bearings. The actual running speed of the motor is 1758.8 r/min, and the actual rotating frequency is 29.86 Hz. After calculation, the characteristic frequency of bearing outer ring fault fo = 70.2 Hz; the characteristic frequency of inner ring fault fi = 188.4 Hz can be obtained. Fig. 23 shows the time domain diagram and envelope spectrum of the original signal. Due to the interference of a large amount of noise, although part of the impact signal can be observed in the time domain, the fault frequency and its doubling cannot be observed in the envelope spectrum, and no more useful information can be extracted.

First, the signal is processed using the SAM method, as shown in Fig. 24. The outer ring fault information can be observed in both the 3D diagram and the 2D diagram, but the inner ring fault information cannot be observed. Only fo can be observed in log-MSES, and all fault information cannot be accurately diagnosed.

Fig. 24Composite diagnosis results of SAM method

Composite diagnosis results of SAM method
Composite diagnosis results of SAM method
Composite diagnosis results of SAM method

The LWTSAM method was then used to process the signal. First, the signal is divided using the filter bank that comes with MODWT. The number of filters is set to 9 and the bandwidth is set to 600Hz. Random impact interference can be successfully separated. According to the parameters, narrowband signal #7 contains fault information, and all narrowband signals are shown in Fig. 25.

Narrow-band signal #7 and its envelope spectrum are shown in Fig. 26, and fault-related information is hidden in harmonic interference.

From the obtained narrowband signal #7, fo and fi can be observed, but the amplitude is relatively low and lower than other irrelevant components. So, we proceed to the next step. WFT is used to reconstruct narrowband signal #7. The reconstructed signal and its envelope spectrum are shown in Fig. 27. It can be seen that the noise is greatly weakened, but the fault component is also weakened. The harmonics of the inner and outer rings can be found respectively in log-MSES., but the amplitude intensity is lower than that of the original narrowband signal, because no suitable MO was selected. The recommended indicators are then used to select the best weights for extracting individual SES, and the minimum value of the recommended indicators appears at MO= 1.1, which is considered the best weights. As shown in Fig. 28, the irrelevant components are weakened and the fault components are enhanced, which is enough to determine that they are inner and outer rings.

Fig. 25Composite fault narrowband signal

Composite fault narrowband signal
Composite fault narrowband signal
Composite fault narrowband signal
Composite fault narrowband signal
Composite fault narrowband signal
Composite fault narrowband signal
Composite fault narrowband signal
Composite fault narrowband signal
Composite fault narrowband signal

Fig. 26Narrowband signal #7 and its envelope spectrum

Narrowband signal #7 and its envelope spectrum
Narrowband signal #7 and its envelope spectrum

Fig. 27Composite reconstructed signal and its envelope spectrum

Composite reconstructed signal and its envelope spectrum
Composite reconstructed signal and its envelope spectrum

Fig. 28Final diagnosis result of LWTSAM composite fault

Final diagnosis result of LWTSAM composite fault
Final diagnosis result of LWTSAM composite fault
Final diagnosis result of LWTSAM composite fault
Final diagnosis result of LWTSAM composite fault

The fast spectral kurtosis method and Autogram method are used to process the signal. Fig. 29 shows the fast spectral kurtosis method. Where Kmax is the maximum kurtosis value, Bw is the bandwidth, and fc is the center frequency. It can be known from the kurtosis diagram that Kmax is 3.9, Bw is 6000 Hz, and fc is 3000 Hz. The best frequency band is found at Level 0. The envelope spectrum is all irrelevant information, and the effect is very poor.

Fig. 29Composite diagnosis results of fast spectral kurtosis method

Composite diagnosis results of fast spectral kurtosis method
Composite diagnosis results of fast spectral kurtosis method

Fig. 30Autogram method composite diagnosis results

Autogram method composite diagnosis results
Autogram method composite diagnosis results

The diagnosis results of the Autogram method are shown in Fig. 30. Where Kmax is the maximum kurtosis value, Bw is the bandwidth, and fc is the center frequency. From the kurtosis diagram, we can know that Kmax is 2434.9, Bw is 1500 Hz, and fc is 2250 Hz. The best frequency band is found at Level 2. The envelope spectrum results under the best frequency band are given. For further verification, the envelope spectra after the rising and falling thresholds are given respectively. It can be observed that only the frequency of fault information is observed in the envelope spectrum after lowering the threshold, and the fault information of the rest of the envelope spectrum is concealed by noise.

To sum up, the LWTSAM method can also show good performance when processing composite fault signals, which verifies the feasibility and superiority of the proposed method.

4.5. Quantitative comparison

Through the above experiments and calculations, a comparative analysis of quantitative indicators is given, mainly including: SNR improvement (dB); Fault-to-noise ratio (FNR); Entropy of envelope spectrum; Detected fault frequencies.

The analysis is shown in Table 3.

Table 3Quantitative indicator comparison

Method
SNR (dB)
FNR
Envelope entropy
Detected fault frequencies
Original
–6.2
0.31
4.82
No
Kurtogram
–1.4
0.58
3.91
Partial
Autogram
–0.8
0.63
3.74
Partial
SAM
0.9
0.71
3.26
Single
LWTSAM
3.8
0.89
2.11
All

It can be seen from the table that the index data of this method are better than the comparison method, which is enough to prove the feasibility of this method.

4.6. Ablation study and statistical validation

Based on the above experiments and results, the ablation experimental analysis is given, as shown in Table 4.

Statistical verification is shown in Table 5, the indicators are as follows: Mean peak amplitude; Standard deviation (Std); Coefficient of Variation (CV).

Table 4Analysis of ablation experiments

Case
MODWT
WFT
MO selection
Result
Case 1
Failed
Case 2
+
Weak
Case 3
+
+
Partial
Case 4
+
+
+
Best

Table 5Statistical verification

Method
Mean peak amplitude
Std
CV
Kurtogram
1.23
0.55
0.44
Autogram
1.37
0.46
0.34
SAM
1.68
0.41
0.33
LWTSAM
2.98
0.29
0.10

5. Conclusions

In this study, a method called LWTSAM was proposed to diagnose bearing faults under complex conditions. First, the original signal is divided into several narrowband signals by using the filter bank included in the MODWT to achieve the effect of separating some interference signals; then, the narrowband signal with the richest fault information is selected and subjected to window Fourier transform to obtain the narrowband signal in the time domain. The amplitude value, and more comprehensive and detailed amplitude information is obtained; then, unbiased autocorrelation and information entropy are used to select the best weight and extract the single square envelope to further reduce the interference of irrelevant components; Finally, the proposed method is applied to outer ring fault, inner ring fault and composite fault signals respectively. The results prove that the method is fast, effective and adaptable, can automatically and accurately extract fault information, and successfully diagnose bearing faults. It also compares it with some classic methods, highlighting the advantages of the research method in this paper. In future research, we will continue to improve the LWTSAM method to expand its application in practical projects. Future work will focus on extending LWTSAM to variable-speed conditions, reducing computational cost for online monitoring, and integrating the proposed framework with data-driven learning models for intelligent fault diagnosis.

References

  • C. He, T. Wu, R. Gu, Z. Jin, R. Ma, and H. Qu, “Rolling bearing fault diagnosis based on composite multiscale permutation entropy and reverse cognitive fruit fly optimization algorithm – Extreme learning machine,” Measurement, Vol. 173, p. 108636, Mar. 2021, https://doi.org/10.1016/j.measurement.2020.108636
  • J. Li, Q. Yu, X. Wang, and Y. Zhang, “An enhanced rolling bearing fault detection method combining sparse code shrinkage denoising with fast spectral correlation,” ISA Transactions, Vol. 102, pp. 335–346, Jul. 2020, https://doi.org/10.1016/j.isatra.2020.02.031
  • G. He, J. Li, K. Ding, and Z. Zhang, “Feature extraction of gear and bearing compound faults based on vibration signal sparse decomposition,” Applied Acoustics, Vol. 189, p. 108604, Feb. 2022, https://doi.org/10.1016/j.apacoust.2021.108604
  • K. Zhang, Y. Xu, and P. Chen, “Feature extraction by enhanced analytical mode decomposition based on order statistics filter,” Measurement, Vol. 173, p. 108620, Mar. 2021, https://doi.org/10.1016/j.measurement.2020.108620
  • L. Xiang, X. Yang, A. Hu, H. Su, and P. Wang, “Condition monitoring and anomaly detection of wind turbine based on cascaded and bidirectional deep learning networks,” Applied Energy, Vol. 305, p. 117925, Jan. 2022, https://doi.org/10.1016/j.apenergy.2021.117925
  • H. Su, L. Xiang, A. Hu, Y. Xu, and X. Yang, “A novel method based on meta-learning for bearing fault diagnosis with small sample learning under different working conditions,” Mechanical Systems and Signal Processing, Vol. 169, p. 108765, Apr. 2022, https://doi.org/10.1016/j.ymssp.2021.108765
  • Y. Liu, S. Yang, and Y. Liao, “An adaptive resonance demodulation method and its application in early fault diagnosis of rolling bearings,” (in Chinese), Journal of Vibration Engineering, Vol. 29, No. 2, pp. 366–370, 2016.
  • X. Wang, J. Zheng, Q. Ni, H. Pan, and J. Zhang, “Traversal index enhanced-gram (TIEgram): A novel optimal demodulation frequency band selection method for rolling bearing fault diagnosis under non-stationary operating conditions,” Mechanical Systems and Signal Processing, Vol. 172, p. 109017, Jun. 2022, https://doi.org/10.1016/j.ymssp.2022.109017
  • B. Cai, L. Zhang, and G. Tang, “Encogram: An autonomous weak transient fault enhancement strategy and its application in bearing fault diagnosis,” Measurement, Vol. 206, p. 112333, Jan. 2023, https://doi.org/10.1016/j.measurement.2022.112333
  • J. Antoni and R. B. Randall, “The spectral kurtosis: application to the vibratory surveillance and diagnostics of rotating machines,” Mechanical Systems and Signal Processing, Vol. 20, No. 2, pp. 308–331, Feb. 2006, https://doi.org/10.1016/j.ymssp.2004.09.002
  • J. Antoni, “The spectral kurtosis: a useful tool for characterising non-stationary signals,” Mechanical Systems and Signal Processing, Vol. 20, No. 2, pp. 282–307, Feb. 2006, https://doi.org/10.1016/j.ymssp.2004.09.001
  • J. Antoni, “The spectral kurtosis of nonstationary signals: formalisation, some properties, and application,” in 12th European Signal Processing Conference, 2004.
  • C. Ma, Z. Sheng, and Y. Xu, “Fault diagnosis of rolling bearings based on adaptive frequency slice wavelet transform,” (in Chinese), Journal of Agricultural Engineering, Vol. 35, No. 10, pp. 34–41, Oct. 2019.
  • T. Barszcz and A. Jabłoński, “A novel method for the optimal band selection for vibration signal demodulation and comparison with the kurtogram,” Mechanical Systems and Signal Processing, Vol. 25, No. 1, pp. 431–451, Jan. 2011, https://doi.org/10.1016/j.ymssp.2010.05.018
  • J. Antoni, “The infogram: entropic evidence of the signature of repetitive transients,” Mechanical Systems and Signal Processing, Vol. 74, pp. 73–94, Jun. 2016, https://doi.org/10.1016/j.ymssp.2015.04.034
  • D. Wang, P. W. Tse, and K. L. Tsui, “An enhanced kurtogram method for fault diagnosis of rolling element bearings,” Mechanical Systems and Signal Processing, Vol. 35, No. 1-2, pp. 176–199, Feb. 2013, https://doi.org/10.1016/j.ymssp.2012.10.003
  • Y. Lei, J. Lin, Z. He, and Y. Zi, “Application of an improved kurtogram method for fault diagnosis of rolling element bearings,” Mechanical Systems and Signal Processing, Vol. 25, No. 5, pp. 1738–1749, Jul. 2011, https://doi.org/10.1016/j.ymssp.2010.12.011
  • A. Moshrefzadeh and A. Fasana, “The Autogram: An effective approach for selecting the optimal demodulation band in rolling element bearings diagnosis,” Mechanical Systems and Signal Processing, Vol. 105, pp. 294–318, May 2018, https://doi.org/10.1016/j.ymssp.2017.12.009
  • X. Gu, S. Yang, Y. Liu, and R. Hao, “Rolling element bearing faults diagnosis based on kurtogram and frequency domain correlated kurtosis,” Measurement Science and Technology, Vol. 27, No. 12, p. 125019, Dec. 2016, https://doi.org/10.1088/0957-0233/27/12/125019
  • Y. Liao, P. Sun, B. Wang, and L. Qu, “Extraction of repetitive transients with frequency domain multipoint kurtosis for bearing fault diagnosis,” Measurement Science and Technology, Vol. 29, No. 5, p. 055012, May 2018, https://doi.org/10.1088/1361-6501/aaae99
  • W. Liu, Y. Liu, and S. Yang, “Fault diagnosis method for rolling bearings based on typical spectral correlation kurtosis diagram,” (in Chinese), Vibration and Shock, Vol. 37, No. 8, pp. 87–92, 2018.
  • P. Borghesani, P. Pennacchi, and S. Chatterton, “The relationship between kurtosis – and envelope-based indexes for the diagnostic of rolling element bearings,” Mechanical Systems and Signal Processing, Vol. 43, No. 1-2, pp. 25–43, Feb. 2014, https://doi.org/10.1016/j.ymssp.2013.10.007
  • X. Zhang, J. Kang, L. Xiao, J. Zhao, and H. Teng, “A new improved kurtogram and its application to bearing fault diagnosis,” Shock and Vibration, Vol. 2015, pp. 1–22, Jan. 2015, https://doi.org/10.1155/2015/385412
  • D. Ibarra-Zarate, O. Tamayo-Pazos, and A. Vallejo-Guevara, “Bearing fault diagnosis in rotating machinery based on cepstrum pre-whitening of vibration and acoustic emission,” The International Journal of Advanced Manufacturing Technology, Vol. 104, No. 9-12, pp. 4155–4168, Aug. 2019, https://doi.org/10.1007/s00170-019-04171-6
  • P. Borghesani, P. Pennacchi, R. B. Randall, N. Sawalhi, and R. Ricci, “Application of cepstrum pre-whitening for the diagnosis of bearing faults under variable speed conditions,” Mechanical Systems and Signal Processing, Vol. 36, No. 2, pp. 370–384, Apr. 2013, https://doi.org/10.1016/j.ymssp.2012.11.001
  • A. Moshrefzadeh, A. Fasana, and J. Antoni, “The spectral amplitude modulation: A nonlinear filtering process for diagnosis of rolling element bearings,” Mechanical Systems and Signal Processing, Vol. 132, pp. 253–276, Oct. 2019, https://doi.org/10.1016/j.ymssp.2019.06.030
  • D. Zhu, B. Yin, and C. Teng, “An improved spectral amplitude modulation method for rolling element bearing fault diagnosis,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, Vol. 45, No. 5, Apr. 2023, https://doi.org/10.1007/s40430-023-04184-z
  • A. Moshrefzadeh, “Condition monitoring and intelligent diagnosis of rolling element bearings under constant/variable load and speed conditions,” Mechanical Systems and Signal Processing, Vol. 149, p. 107153, Feb. 2021, https://doi.org/10.1016/j.ymssp.2020.107153
  • W. Liu, S. Yang, Q. Li, Y. Liu, and X. Gu, “An enhanced spectral amplitude modulation method and its application to rolling element bearings fault diagnosis under complex interference,” (in Chinese), Journal of Vibration Engineering, Vol. 34, pp. 1064–1075, 2021.
  • Z. Wang, J. Yang, and Y. Guo, “Unknown fault feature extraction of rolling bearings under variable speed conditions based on statistical complexity measures,” Mechanical Systems and Signal Processing, Vol. 172, p. 108964, Jun. 2022, https://doi.org/10.1016/j.ymssp.2022.108964
  • Z. Du, X. Chen, H. Zhang, and Y. Yang, “Low-rank enhanced convolutional sparse feature detection for accurate diagnosis of gearbox faults,” Mechanical Systems and Signal Processing, Vol. 150, p. 107215, Mar. 2021, https://doi.org/10.1016/j.ymssp.2020.107215
  • X. Zhong, Q. Mei, X. Gao, and T. Huang, “Fault diagnosis of rolling bearings based on improved direct fast iterative filtering and spectral amplitude modulation,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, Vol. 236, No. 9, pp. 5111–5123, Jan. 2022, https://doi.org//10.1177/09544062211055657
  • Z. Ma, M. Zhao, X. Dai, B. Ma, and H. Bi, “Optimized spectral amplitude modulation based on generalized envelope for bearing compound fault diagnosis,” in 5th International Conference on Intelligent Control, Measurement and Signal Processing (ICMSP), pp. 1–6, May 2023, https://doi.org/10.1109/icmsp58539.2023.10171055
  • A. G. Hafez and E. Ghamry, “Geomagnetic sudden commencement automatic detection via MODWT,” IEEE Transactions on Geoscience and Remote Sensing, Vol. 51, No. 3, pp. 1547–1554, Mar. 2013, https://doi.org/10.1109/tgrs.2012.2207962
  • K. Zhang, P. Chen, M. Yang, L. Song, and Y. Xu, “The Harmogram: A periodic impulses detection method and its application in bearing fault diagnosis,” Mechanical Systems and Signal Processing, Vol. 165, p. 108374, Feb. 2022, https://doi.org/10.1016/j.ymssp.2021.108374
  • D. Wang, Y. Zhao, C. Yi, K.-L. Tsui, and J. Lin, “Sparsity guided empirical wavelet transform for fault diagnosis of rolling element bearings,” Mechanical Systems and Signal Processing, Vol. 101, pp. 292–308, Feb. 2018, https://doi.org/10.1016/j.ymssp.2017.08.038
  • Z. Jiang, K. Zhang, L. Xiang, G. Yu, and Y. Xu, “A time-frequency spectral amplitude modulation method and its applications in rolling bearing fault diagnosis,” Mechanical Systems and Signal Processing, Vol. 185, p. 109832, Feb. 2023, https://doi.org/10.1016/j.ymssp.2022.109832
  • S. Shao, S. Mcaleer, R. Yan, and P. Baldi, “Highly accurate machine fault diagnosis using deep transfer learning,” IEEE Transactions on Industrial Informatics, Vol. 15, No. 4, pp. 2446–2455, Apr. 2019, https://doi.org/10.1109/tii.2018.2864759

About this article

Received
December 20, 2025
Accepted
February 11, 2026
Published
February 28, 2026
SUBJECTS
Fault diagnosis based on vibration signal analysis
Keywords
spectral amplitude modulation
time-frequency analysis
residual signal
rolling bearing
fault diagnosis
Acknowledgements

This research work was supported by the National Natural Science Foundation of China (611175054); Natural Science Foundation of Zhejiang Province (LY16E050001); Cixi City Science and Technology Plan Project:(2023CXS080009); Cixi City Science and Technology Project (CZ2023006); Cixi City Science and Technology Project (CZ2025003).

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author Contributions

Zongcai Ma: writing-original draft preparation, writing-review and editing. Yongqi Chen: conceptualization. Qinge Dai: data curation. Linqiang Wu: software. Yitong Qin: formal analysis.

Conflict of interest

The authors declare that they have no conflict of interest.