# An Integrated Model of Atom Search Optimization-Based Resonance Sparse Signal Decomposition and Cross-Validation SVM for Gearbox Fault Diagnosis

Authors: Fengfeng Bie; Yifan Wu; Ying Zhang; Jian Peng; Hongfei Zhu
Journal: Shock and Vibration
 (2022)

---

## Abstract

In the aspect of gearbox fault diagnosis, the periodic pulse signal containing fault characteristics is often overwhelmed with other strong interference components, which brings a great challenge for gearbox fault detection and status identification. To address these issues, this paper develops a novel resonance-based sparse signal decomposition (RSSD) fault diagnosis method combined with support vector machine (SVM). Based on the key of the decomposition parameters in the resonant sparse decomposition method, the influence of atom search optimization (ASO) on the quality factor in the resonance-based sparse signal decomposition method is primarily studied. The vibration signal of the gearbox was decomposed by resonance sparse decomposition with the optimal quality factor, from which the high and low-resonance components were obtained. Then, the power spectrum entropy, singular spectrum entropy, and time-domain energy entropy of the low-resonance signal were calculated. Finally, the pattern recognition of gearbox fault was completed with a designed cross-validation SVM pattern recognition model. The numerical simulation and gearbox fault experiments demonstrate that the presented method can achieve great recognition accuracies and effect capacities against interference components involved in the gearbox vibration signal.


---

## Body

## 1. Introduction

As the core component of mechanical power and motion transmission, gearbox drives are widely used in metallurgical, petrochemical, mining, transportation, and other mechanical systems. The status monitoring and fault diagnosis of gearbox are critically important to ensure the reliability and efficiency of the mechanical system. In recent years, more and more attention has been paid to the research of gear state monitoring and fault diagnosis technology because of its relatively complex, nonlinear, and nonstationary dynamic behavior. Since the failure mode of the gearbox could be mostly figured out as periodic pulse characteristic, while the measured fault signals are bound to be overwhelmed with harmonic signals and background noise signals generated in the gearbox system, the core issue of the diagnosis process is highlighted as periodic pulse feature extraction and evaluation [1].As mentioned in most related vibration-based studies, the fault diagnosis was processed in three main forms [2], i.e., the vibration separation, amplitude-frequency demodulation, and the advanced signal processing-based algorithms. Vibration separation-based algorithms are proposed to deal with the complex structure of the planetary gearbox aiming at eliminating the unrelated vibration from original signal through the time synchronized average (TSA) algorithm and its improved versions [3]. Another way for gearbox fault diagnosis was developed as identifying the amplitude and quantity increment of fault feature sideband around the meshing frequency and the corresponding harmonic in the spectrum. As the research focus, typical signal processing methods are developed in related fields [4], such as, empirical mode decomposition (EMD), ensemble empirical mode decomposition (EEMD), local mean decomposition (LMD), variational mode decomposition (VMD), and intrinsic time-scale decomposition (ITD) [5]. In addition, scholars have also made great progress on methods based on advanced signal processing, such as spectral kurtosis (SK)-based algorithm [6], continuous wavelet transformation (CWT) based on wavelet [7], algorithm adaptive stochastic resonance (ASR) based on stochastic resonance, and so on [8]. Besides, sparse representation-based algorithms are often used, as well as other processing methods, such as dynamic time warping (DTW), adaptive integrated empirical mode decomposition, and a series of algorithms [9]. In most of these algorithms, the waveform characteristics of the components are taken as decomposition gist, which is hard to work when the signal involves strong nonlinear nonstationary characteristics with background noise. To address this issue, resonance-based sparse signal decomposition is proposed [10], in which the vibration signal is decomposed into the high-resonance component through the tunable Q wavelet transform, and the signal components with overlapping frequencies can be separated. This method breaks through the limitation of traditional decomposition method based on waveform and each component characteristic. When a defect emerges in the inner gear of the gearbox, periodic pulses are represented by corresponding fault vibration signals. Due to the coupling between the working environment and other mechanical components, the interference of noise and harmonic is involved in the collected gearbox vibration signal, which is shown as continuous oscillation, while the influence of the fault has lower oscillatory behavior unlike the harmonics and noise. RSSD has a more significant effect on feature extraction of nonlinear nonstationary signals of gearbox. A method based on resonance-based sparse signal decomposition integrated with empirical mode decomposition and demodulation was introduced in the literature [11] that can separate the impulse component from the signal, according to the different Q-factors of impulse component and harmonic component. According to the rotating speed signal, the time-domain impact signal of the gearbox was resampled at constant angle increments. Sun et al. [12] firstly applied the method of RSSD to the gearbox feature extraction process. With the frequency spectrum of the resampled collision signal being analyzed, the step domain analysis is completed, and the final diagnosis result is obtained accordingly. The proposed approach is of good anti-noise ability and is suitable for analyzing the actual vibration signal of a gearbox with rotating speed fluctuation. Tang et al. [13] proposed a fault diagnosis approach for planetary gearboxes using a multi-resonance component fusion-based convolutional neural network; RSSD is used to decompose the vibration signal, and the low-resonance components that may contain the impulse components of bearing faults can be retained. However, the final decomposition results might be greatly influenced by the choice of high-quality and low-quality factors for resonance-based sparse signal decomposition if the high-quality factor and low-quality factor are artificially determined, which is often inaccurate and even brings unreasonable final decomposition result. Thus, the ASO algorithm is adopted to determine the optimal quality factor in the study. The SVM model has been widely applied in mechanical fault pattern recognition, while for the process of gearbox fault recognition, the accuracy of SVM is regularly affected by C (penalty factor) and G (variance in RBF kernel function). Hence, the cross-validation method is employed to determine the final values of C and G to improve the accuracy of the final pattern recognition. Extensive shallow learning models, such as support vector machine (SVM), Naive Bayes classifier, and geometric model, were designed with the extracted inherent fault features from vibration signals [14]. SVM is widely used in the field of fault diagnosis.Based on the above analysis, this paper proposes a new method for gearbox fault identification with improved RSSD and CV-SVM model. The main research contents include the following. The study of ASO made by the RSSD of low kurtosis of resonance components reached a maximum, and high-resonance decomposition is processed for the optimal parameters ofQ1 and Q2 of low-resonance decomposition parameters. Then, power spectrum entropy, singular spectrum entropy, and time-domain energy entropy are analyzed in the feature extraction of low-resonance components for pattern recognition. The C/G parameters of SVM are optimized with cross-validation, and the accuracy of pattern recognition is promoted with this method. The simulated fault signal is obtained through the digital model for the gear fault. Through model simulation, the simulated fault signal is obtained, and the feasibility of the method is verified. The rest of this paper is organized as follows. Section 2 states some theoretical background and presents the fault diagnosis process. Section 3 presents the details of gearbox fault simulation. In Section 4, the feasibility and superiority of this method are verified on three different types of gearbox fault data and compared with the nonoptimized diagnosis technology. Finally, the conclusion and future jobs are summarized in Section 5.


## 2. Algorithm Principle

### 2.1. RSSD Optimization Based on ASO

In general, vibration signals can be decomposed into three signals by RSSD, i.e., the high-resonance signal series with continuous oscillating behavior, the low-resonance components of transient impact signals without specific shape and duration, and the remaining residual component of the original signal. This is achieved by the tunable Q-factor wavelet transform (TQWT). The greater quality factorQ value is obtained for the better aggregation performance on frequency domain of the high-resonance signal. Low-resonance signal has good time-domain aggregation, and its quality factor Q value obtained is smaller [15]. The multi-channel bandpass decomposition filter banks are employed in the process of the tunable quality factor wavelet transform on the acquired gearbox vibration signals, and the basis function libraries of high Q-transform and low Q-transform are obtained. The principle of the method is shown in Figure 1, where H0ω represents the low-resonance filter, H1ω represents the high pass filter, V0n and V1n are the sub-band signals after filtering, and α and β represent the low-pass and high-pass scale factors, respectively. The quality factor Q, redundancy factor R, and the maximum number of decomposition layers Lmax can be set referring to equations (1)–(3), [16].Figure 1
Bandpass filter banks. (a) Decomposition filter banks. (b) Reconstructed filter.
(a)(b)(1)Q=2−ββ,(2)r=β1−α,(3)Lmax≤logβN/8log1/α.Therefore, the acquired vibration signalx can be expressed as the composition of the high-resonance component x1, the low-resonance component x2, and the participating component n, which can be expressed as [17](4)x=x1+x+2n,where x1 and x2 can be represented by overcomplete basis function libraries S1 and S2, which can be constructed by high-quality factor Q1 and low-quality factor Q2, respectively, and the objective function can be established by constructing MCA equation [18]:(5)JW1,W2=x−S1W1−S2W222+λ1W11+λ2W21,where W1 and W2 are the transformation coefficients under the basis function libraries S1 and S2, respectively, and λ1 and λ2 are the regularization parameters.Equation (5) can be reorganized as(6)JW1,W2=x−S1W1−S2W222+∑j=1J1+1λ1,jW1,j1+∑j=1J1+1λ2,jW2,j1.To achieve the minimum function value, the separation augmented Lagrange contraction algorithm (SALSA) is adopted for iterative calculation.In resonance-based sparse signal decomposition, the choice of quality factorsQ1 and Q2 plays a decisive role in the final decomposition result. The prior knowledge for the selection of quality factor parameters is often needed by the traditional resonance-based sparse signal decomposition method, which greatly increases the cost of operation. To solve the above problems, this paper proposes a quality factor optimization scheme based on ASO. The kurtosis of low-resonance component decomposed by RSSD can be taken as the objective function by this method. The optimal quality factor generates as the kurtosis value reaches the maximum level, since the value of kurtosis could correspondingly reflect the normal distribution of vibration signal [19]. The larger the noise, the smaller the kurtosis value and the less obvious the response. Kurtosis can be used as an important index to measure the number of pulse components in the signal. In the literature [19], kurtosis is effectively applied as an important index to measure the fault of reciprocating compressor. The expression of kurtosis is as follows:(7)K=Ex−μ4σ4,where μ represents the mean value of the vibration signal, σ represents the standard deviation of the vibration signal, and E is the mathematical expectation.Atom search optimization is a novel intelligent algorithm based on molecular dynamics model proposed in 2019 [20]. In a molecular system, there are interaction forces (attraction and repulsion) between neighboring atoms, and the globally optimal atoms have geometric constraints on other atoms. Gravitation urges atoms to explore the whole search space extensively, and repulsion enables them to develop potential regions effectively. It has the characteristics of strong optimization ability and fast convergence.The acceleration of thei−th atom at time t can be expressed as(8)aidt=Fidtmidt+Gidtmidt=−δ1−t−1Tmax3e−20t/Tmax,∑j∈Kbestrandj2×hijt13−hijt7mit,xidt−xjdtxit,xjt+γe−20t/Tmaxxbestdt−xidtmit,where δ represents depth weighting, Tmax represents the maximum number of iterations, hijt represents the distance between the two atoms, γ is the multiplier weight, and mit is the mass of the i−th atom after the t−th iteration. The velocity and position of the atomic movement can be changed by the acceleration of atoms:(9)vidt+1=randidvidt+aidt,xidt+1=xidt+vidt+1,where xidt is the position of the i−th atom after t−th iterations, vidt is the velocity of the i−th atom after t−th iterations, and aidt is the acceleration of the i−th atom after t−th iterations.In order to improve the global search ability of ASO algorithm in the initial iteration, every atom should interact with as many other optimal atoms as possible. To improve its search ability in the later iteration, every atom should interact as few as possible.


### 2.2. Cross-Validation SVM

SVM is usually used to classify problems with multiple categories [21]. The SVM method was first applied to the fault diagnosis field of gearbox, and its fault classification effect is relatively ideal. In reference [22], it is proved that the cross-validation method can effectively improve the classification accuracy of SVM, and this method is significantly better than the traditional SVM method. In order to obtain a better classification effect, CV-SVM was used for pattern recognition of single fault of gearbox.For a training samplex1,y1,x2,y2,…,xk,yk set with a sample number of k, where xi∈Rd represents the input index and yi∈−1,1 represents the output index (where i = 1,2, …, k), if the optimal hyperplane is set to w∗x+b=0, the weight vector w and bias b must satisfy the following constraints:(10)yiw∗x+b≥1−ζi,where ζi is the relaxation variable and represents the degree of deviation between the SVM classification model and the ideal linear case. The hyperplane decision function of the linear classifier is(11)Mx=sgn∑inci∗yiχ•χ+b∗,where sgn represents the function, ci∗ represents the optimal solution of ci, b∗ is the threshold value of classification, χ is the samples to be classified, and Mx represents the classification function, where the category of the classification sample χ is judged according to the positive or negative value of the classification function.Cross-validation is a process commonly applied in validating the performance of a given classifier, in which the data are firstly grouped, and then the validation set is composed of the training set and the corresponding training parameters [23]. The best parameter C and parameter G can be found with cross-validation for the selection of the best parameter of the training model. In the process of vibration and impact feature extraction of gearbox, in order to extract the most significant fault feature as much, the entropy feature should be extracted from multiple ways. Reference [24] elaborates the advantages of extracting entropy features from multiple aspects over extracting entropy features from a single point. Meanwhile, singular spectral entropy, time-domain energy entropy, and power spectral entropy could reflect signal features in time domain and frequency domain. Therefore, only three entropy features are adopted in this paper to construct feature vectors. The complexity and uncertainty of frequency domain signal can be reflected by power spectrum entropy. The more evenly the energy in the signal is distributed in the frequency domain, the more complex the signal is and the higher the degree of uncertainty is. As a kind of information entropy, singular spectrum entropy reflects the uncertainty degree of each mode of time-domain signal under singular spectrum division [25]. As the gearbox fails, a periodic excitation generates in the gearbox system that could be involved in the gearbox vibration. The periodic excitation in the time-domain signal and the complexity and uncertainty of the time-domain signal energy can be manifested with the time-domain energy entropy. All the symbols and variables are explained in the nomenclature.


### 2.3. Gear Fault Diagnosis Process Based on RSSD Optimized by ASO Algorithm Combined with CV-SVM

The specific process is as follows.Step 1.
The maximum decomposition levelL1 of high-resonance component is set as 30, and the maximum decomposition level L2 of low-resonance component is 12; the redundancy r1 of high-resonance component is 9, and the redundancy r2 of low-resonance component is 3.Step 2.
Given the initial high-resonance quality factorQ1 value 4 and low-resonance quality factor Q2 value 2, the values of the high Q-factor and low Q-factor are set to be 0.8–1.3 and 3–9, respectively.Step 3.
Set the population number of the ASO to be 50, the maximum iteration times are 1000 times, the depth weight is 50, and the multiplier weight is 0.2.Step 4.
Set the target function as the kurtosis of low-resonance component. When the kurtosis reaches the maximum value of iteration, the process stops and theQ1 and Q2 are output. Otherwise, repeat Steps 2 and 3.Step 5.
To obtain the low-resonance component and high-resonance component, the power spectrum entropy, singular spectrum entropy, and time-domain energy entropy of the low-resonance component are calculated. These three kinds of feature entropy are substituted into SVM of cross-validation to recognize patterns.
The specific process is shown in Figure2.Figure 2
Flowchart of fault diagnosis and pattern recognition of gearbox with resonance-based sparse signal decomposition based on ASO.


## 2.1. RSSD Optimization Based on ASO

In general, vibration signals can be decomposed into three signals by RSSD, i.e., the high-resonance signal series with continuous oscillating behavior, the low-resonance components of transient impact signals without specific shape and duration, and the remaining residual component of the original signal. This is achieved by the tunable Q-factor wavelet transform (TQWT). The greater quality factorQ value is obtained for the better aggregation performance on frequency domain of the high-resonance signal. Low-resonance signal has good time-domain aggregation, and its quality factor Q value obtained is smaller [15]. The multi-channel bandpass decomposition filter banks are employed in the process of the tunable quality factor wavelet transform on the acquired gearbox vibration signals, and the basis function libraries of high Q-transform and low Q-transform are obtained. The principle of the method is shown in Figure 1, where H0ω represents the low-resonance filter, H1ω represents the high pass filter, V0n and V1n are the sub-band signals after filtering, and α and β represent the low-pass and high-pass scale factors, respectively. The quality factor Q, redundancy factor R, and the maximum number of decomposition layers Lmax can be set referring to equations (1)–(3), [16].Figure 1
Bandpass filter banks. (a) Decomposition filter banks. (b) Reconstructed filter.
(a)(b)(1)Q=2−ββ,(2)r=β1−α,(3)Lmax≤logβN/8log1/α.Therefore, the acquired vibration signalx can be expressed as the composition of the high-resonance component x1, the low-resonance component x2, and the participating component n, which can be expressed as [17](4)x=x1+x+2n,where x1 and x2 can be represented by overcomplete basis function libraries S1 and S2, which can be constructed by high-quality factor Q1 and low-quality factor Q2, respectively, and the objective function can be established by constructing MCA equation [18]:(5)JW1,W2=x−S1W1−S2W222+λ1W11+λ2W21,where W1 and W2 are the transformation coefficients under the basis function libraries S1 and S2, respectively, and λ1 and λ2 are the regularization parameters.Equation (5) can be reorganized as(6)JW1,W2=x−S1W1−S2W222+∑j=1J1+1λ1,jW1,j1+∑j=1J1+1λ2,jW2,j1.To achieve the minimum function value, the separation augmented Lagrange contraction algorithm (SALSA) is adopted for iterative calculation.In resonance-based sparse signal decomposition, the choice of quality factorsQ1 and Q2 plays a decisive role in the final decomposition result. The prior knowledge for the selection of quality factor parameters is often needed by the traditional resonance-based sparse signal decomposition method, which greatly increases the cost of operation. To solve the above problems, this paper proposes a quality factor optimization scheme based on ASO. The kurtosis of low-resonance component decomposed by RSSD can be taken as the objective function by this method. The optimal quality factor generates as the kurtosis value reaches the maximum level, since the value of kurtosis could correspondingly reflect the normal distribution of vibration signal [19]. The larger the noise, the smaller the kurtosis value and the less obvious the response. Kurtosis can be used as an important index to measure the number of pulse components in the signal. In the literature [19], kurtosis is effectively applied as an important index to measure the fault of reciprocating compressor. The expression of kurtosis is as follows:(7)K=Ex−μ4σ4,where μ represents the mean value of the vibration signal, σ represents the standard deviation of the vibration signal, and E is the mathematical expectation.Atom search optimization is a novel intelligent algorithm based on molecular dynamics model proposed in 2019 [20]. In a molecular system, there are interaction forces (attraction and repulsion) between neighboring atoms, and the globally optimal atoms have geometric constraints on other atoms. Gravitation urges atoms to explore the whole search space extensively, and repulsion enables them to develop potential regions effectively. It has the characteristics of strong optimization ability and fast convergence.The acceleration of thei−th atom at time t can be expressed as(8)aidt=Fidtmidt+Gidtmidt=−δ1−t−1Tmax3e−20t/Tmax,∑j∈Kbestrandj2×hijt13−hijt7mit,xidt−xjdtxit,xjt+γe−20t/Tmaxxbestdt−xidtmit,where δ represents depth weighting, Tmax represents the maximum number of iterations, hijt represents the distance between the two atoms, γ is the multiplier weight, and mit is the mass of the i−th atom after the t−th iteration. The velocity and position of the atomic movement can be changed by the acceleration of atoms:(9)vidt+1=randidvidt+aidt,xidt+1=xidt+vidt+1,where xidt is the position of the i−th atom after t−th iterations, vidt is the velocity of the i−th atom after t−th iterations, and aidt is the acceleration of the i−th atom after t−th iterations.In order to improve the global search ability of ASO algorithm in the initial iteration, every atom should interact with as many other optimal atoms as possible. To improve its search ability in the later iteration, every atom should interact as few as possible.


## 2.2. Cross-Validation SVM

SVM is usually used to classify problems with multiple categories [21]. The SVM method was first applied to the fault diagnosis field of gearbox, and its fault classification effect is relatively ideal. In reference [22], it is proved that the cross-validation method can effectively improve the classification accuracy of SVM, and this method is significantly better than the traditional SVM method. In order to obtain a better classification effect, CV-SVM was used for pattern recognition of single fault of gearbox.For a training samplex1,y1,x2,y2,…,xk,yk set with a sample number of k, where xi∈Rd represents the input index and yi∈−1,1 represents the output index (where i = 1,2, …, k), if the optimal hyperplane is set to w∗x+b=0, the weight vector w and bias b must satisfy the following constraints:(10)yiw∗x+b≥1−ζi,where ζi is the relaxation variable and represents the degree of deviation between the SVM classification model and the ideal linear case. The hyperplane decision function of the linear classifier is(11)Mx=sgn∑inci∗yiχ•χ+b∗,where sgn represents the function, ci∗ represents the optimal solution of ci, b∗ is the threshold value of classification, χ is the samples to be classified, and Mx represents the classification function, where the category of the classification sample χ is judged according to the positive or negative value of the classification function.Cross-validation is a process commonly applied in validating the performance of a given classifier, in which the data are firstly grouped, and then the validation set is composed of the training set and the corresponding training parameters [23]. The best parameter C and parameter G can be found with cross-validation for the selection of the best parameter of the training model. In the process of vibration and impact feature extraction of gearbox, in order to extract the most significant fault feature as much, the entropy feature should be extracted from multiple ways. Reference [24] elaborates the advantages of extracting entropy features from multiple aspects over extracting entropy features from a single point. Meanwhile, singular spectral entropy, time-domain energy entropy, and power spectral entropy could reflect signal features in time domain and frequency domain. Therefore, only three entropy features are adopted in this paper to construct feature vectors. The complexity and uncertainty of frequency domain signal can be reflected by power spectrum entropy. The more evenly the energy in the signal is distributed in the frequency domain, the more complex the signal is and the higher the degree of uncertainty is. As a kind of information entropy, singular spectrum entropy reflects the uncertainty degree of each mode of time-domain signal under singular spectrum division [25]. As the gearbox fails, a periodic excitation generates in the gearbox system that could be involved in the gearbox vibration. The periodic excitation in the time-domain signal and the complexity and uncertainty of the time-domain signal energy can be manifested with the time-domain energy entropy. All the symbols and variables are explained in the nomenclature.


## 2.3. Gear Fault Diagnosis Process Based on RSSD Optimized by ASO Algorithm Combined with CV-SVM

The specific process is as follows.Step 1.
The maximum decomposition levelL1 of high-resonance component is set as 30, and the maximum decomposition level L2 of low-resonance component is 12; the redundancy r1 of high-resonance component is 9, and the redundancy r2 of low-resonance component is 3.Step 2.
Given the initial high-resonance quality factorQ1 value 4 and low-resonance quality factor Q2 value 2, the values of the high Q-factor and low Q-factor are set to be 0.8–1.3 and 3–9, respectively.Step 3.
Set the population number of the ASO to be 50, the maximum iteration times are 1000 times, the depth weight is 50, and the multiplier weight is 0.2.Step 4.
Set the target function as the kurtosis of low-resonance component. When the kurtosis reaches the maximum value of iteration, the process stops and theQ1 and Q2 are output. Otherwise, repeat Steps 2 and 3.Step 5.
To obtain the low-resonance component and high-resonance component, the power spectrum entropy, singular spectrum entropy, and time-domain energy entropy of the low-resonance component are calculated. These three kinds of feature entropy are substituted into SVM of cross-validation to recognize patterns.
The specific process is shown in Figure2.Figure 2
Flowchart of fault diagnosis and pattern recognition of gearbox with resonance-based sparse signal decomposition based on ASO.


## 3. Analogue Simulation

In order to verify the rationality of the method proposed above, a single fault model of gearbox is established. The research object of this paper is a straight bevel gear, whose structural parameters are shown in Table1: the face width coefficient ΦR=0.3, the coefficient of addendum height ha∗=1.0, and the coefficient of clearance c∗=0.2. The failure of driving gear is shown in Figure 3. According to the parameter set, the system of 3D model of gear, gear shaft, bearing, and gearbox is established and assembled by SolidWorks. The structural parameters are the same as those of the experimental gearbox shown in Figure 4.Table 1
Parameters of spur bevel gears.
ParameterDriving gearDriven gearNumber of teeth1727Diametral pitch1010Reference circle1.72.7Pressure angle2020Reference cone angle32.257.8Face width0.50.5Addendum0.130.07Dedendum0.090.15Tip angle37.4461.07Root angle28.9352.56Figure 3
Tooth fracture of driving gear.Figure 4
Structure of the gearbox.The model is imported into ADAMS for simulation. The sensor is set above the meshing position of the gear. The data of radial vibration of gearbox were measured. In order to make the simulation signal more realistic, the case of the gearbox box is set to be flexible before the simulation starts, the simulation time is 1.6 s, and the step is set to 2000. The time-domain diagram of the final simulation signal is shown in Figure5.Figure 5
Time-domain simulation signal.The ASO is employed to find the optimal high Q-factorQ1 and low Q-factor Q2. After the iterative search of atom search algorithm, the high Q-factor Q1=7.1767 and the low Q-factor Q2=1.2184 are finally determined. In the subsequent resonance-based sparse signal decomposition, Q1 and Q2 are applied to process the simulation signal. The high-resonance component, low-resonance component, and residual component are obtained as shown in Figure 6. Conventional parameters using the time domain of the vibration signal (kurtosis and skewness) were used for detecting and diagnosing the faults by applying them to filtered signals [26]. In order to compare the results of kurtosis and skewness as optimization objectives, RSSD with skewness as optimization objective is applied to decompose the simulation signal, and the decomposition result is shown in Figure 7.Figure 6
Resonance-based sparse signal decomposition with kurtosis as optimization objective.Figure 7
Resonance-based sparse signal decomposition with skewness as optimization objective.Figure8 illustrates that the obvious peaks could be found at the meshing frequency and higher harmonics. For example, the amplitude reaches 0.008288 m/s2 at fc=1146Hz, which obviously indicates the local abnormality of the gearbox. The spectral lines of 1, 2, and 3 times the conversion frequency show that the modulation frequency of the impact pulse reaches the maximum value, which is consistent with the rotation frequency of the gear shaft; therefore, tooth fracture mode of the gearbox could be identified. Envelope spectrum of RSSD with skewness as optimization objective is shown in Figure 9. From the comparison of envelope spectra of these two typical indicators, it can be obviously observed that the RSSD with kurtosis performs better in the signal decomposition process.Figure 8
Envelope spectrum of the low-resonance component with kurtosis as optimization objective.Figure 9
Envelope spectrum of the low-resonance component with skewness as optimization objective.


## 4. Experimental Validation

As shown in Figures10(a) and 10(b), the experimental system is established for the validation of the proposed method based on SQI-MFS mechanical fault simulation platform. The sensor is set above the meshing position of the gear, the motor frequency is set to 15 Hz, and the sampling frequency is set to 5120 Hz. The radial vibration signal of gearbox is measured. Three typical failure modes of the gearbox are simulated, i.e., as abrasive, notched gear, and tooth fracture. The fault setting of the driving gear is shown in Figure 11.Figure 10
1, motor drive; 2, gearbox; 3, acceleration sensor; 4, data acquisition instrument; 5, laptop and signal analysis software. (a) Test rig. (b) Data acquisition system sketch.
(a)(b)Figure 11
Failure mode setting.The collected vibration signals of normal gearbox, abrasion, notched gear, and tooth fracture are shown in Figure12. In the normal mode, the vibration amplitude is small and stable without obvious pulse feature, while in the other three fault cases, the amplitude increases significantly, and the notched gear’s signal is particularly obvious and its impact component in the signal is also distinctive.Figure 12
The collected vibration signals. (a) Time-domain signal of normal gear. (b) Time-domain signal of abrasion. (c) Time-domain signal of notched gear. (d) Time-domain signal of tooth fracture.
(a)(b)(c)(d)The resonance-based sparse signal decomposition parameters of the three kinds of failure signals were optimized by ASO. In the case of abrasion, the original high Q-factorQ1 was set to 4, and the low Q-factor Q2 was set to 2 initially. With the ASO optimized, when the high Q-factor Q1=7.2767 and the low Q-factor Q2=1.2184, the kurtosis of the low-resonance component can reach the maximum. Similarly, Q1=6.8516 and Q2=1.0297 for notched gear, and Q1=6.9860 and Q2=1.0004 for tooth failure. The optimal quality factor is used for resonance sparse decomposition, and the results are shown in Figure 13.Figure 13
Resonance-based sparse signal decomposition of fault signals. (a) Resonance-based sparse signal decomposition result of vibration signal of abrasion. (b) Resonance-based sparse signal decomposition result of vibration signal of notched gear. (c) Resonance-based sparse signal decomposition result of vibration signal of tooth fracture.
(a)(b)(c)The envelope spectrum of the low-resonance component obtained by the resonance-based sparse signal decomposition is shown in Figure14. The existence of faults can be observed. However, the specific fault type cannot be well distinguished. In view of this, the entropy feature of each group of fault signals is extracted for CV-SVM.Figure 14
Envelope spectrum of low-resonance component. (a) Envelope spectrum of low-resonance component for abrasion. (b) Envelope spectrum of low-resonance component for notched gear. (c) Envelope spectrum of low-resonance component for tooth fracture.
(a)(b)(c)30 groups of vibration signals of normal, abrasion, notched gear, and tooth fracture are extracted, respectively, and the resonance-based sparse signal decomposition based on the optimization of ASO is processed. As the resonance-based sparse signal decomposition of various fault signals is optimized, the low-resonance components are obtained, and their power spectrum entropy, singular spectrum entropy, and time-domain energy entropy are extracted as the eigenvalues, respectively. 36 groups of data were randomly selected as the training set from 120 groups, and the remaining 84 were selected as the testing set. The normal state was marked as +1, the wear state was marked as +2, the missing tooth state was marked as +3, and the broken tooth state was marked as +4. The results of CV-SVM fault recognition are shown in Figure15, and the final accuracy rate of the proposed method is 98.8095%, while the accuracy rate of CV-SVM which uses the original vibration data for feature extraction is only 85.7143%. As shown in Figure 15, compared with the optimization method, the effectiveness of the original method is not as ideal as that of the proposed method.Figure 15
Identification result of CV-SVM model. (a) Identification result of CV-SVM model of processed signal. (b) Identification result of CV-SVM of model of original signals.
(a)(b)


## 5. Conclusion

In this paper, an improved gearbox fault diagnosis method based on optimized resonant sparse signal decomposition combined with CV-SVM is developed. Through simulation and experimental verification, the following conclusions are drawn:(1)
The atom search optimization algorithm is employed in the optimization of the resonance-based sparse signal decomposition. The kurtosis of the low-resonance component is taken as the optimization objective parameter, with which the decomposition ability of the resonant sparse signal decomposition method for fault signals was improved. The effectiveness of this method is verified with simulation and experiments.(2)
The power spectrum entropy, singular spectrum entropy, and time-domain energy entropy are used as eigenvalues to apply to CV-SVM, which distinctively improved the accuracy of gearbox fault pattern recognition.(3)
The integrated model of ASO-based resonance sparse signal decomposition and CV- SVM is effective and feasible for gearbox fault diagnosis.In the future, we will further explore the case of gearbox under compound fault mode, and the effectiveness of the proposed method could be further evaluated or optimized.



---

*Source: 4069180-2022-02-27.xml*