CMES CMES CMES Computer Modeling in Engineering & Sciences 1526-1506 1526-1492 Tech Science Press USA 18123 10.32604/cmes.2022.018123 Article Fault Analysis of Wind Power Rolling Bearing Based on EMD Feature Extraction Fault Analysis of Wind Power Rolling Bearing Based on EMD Feature Extraction Fault Analysis of Wind Power Rolling Bearing Based on EMD Feature Extraction MengDebiao 1 2 3dbmeng@uestc.edu.cn WangHongtao 1 YangShiyuan 1 LvZhiyuan 1 HuZhengguo 1 WangZihao 1 School of Mechanical and Electrical Engineering, University of Electronic Science and Technology of China, Chengdu, 611731, China Institute of Electronic and Information Engineering of UESTC in Guangdong, Dongguan, 523808, China Yangzhou Yangjie Electronic Technology Co., Ltd., Yangzhou, 225008, China Corresponding Author: Debiao Meng. Email: dbmeng@uestc.edu.cn 24 11 2021 130 1 543 558 3062021 0482021 © 2022 Meng et al. 2022 Meng et al. This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In a wind turbine, the rolling bearing is the critical component. However, it has a high failure rate. Therefore, the failure analysis and fault diagnosis of wind power rolling bearings are very important to ensure the high reliability and safety of wind power equipment. In this study, the failure form and the corresponding reason for the failure are discussed firstly. Then, the natural frequency and the characteristic frequency are analyzed. The Empirical Mode Decomposition (EMD) algorithm is used to extract the characteristics of the vibration signal of the rolling bearing. Moreover, the eigenmode function is obtained and then filtered by the kurtosis criterion. Consequently, the relationship between the actual fault frequency spectrum and the theoretical fault frequency can be obtained. Then the fault analysis is performed. To enhance the accuracy of fault diagnosis, based on the previous feature extraction and the time-frequency domain feature extraction of the data after EMD decomposition processing, four different classifiers are added to diagnose and classify the fault status of rolling bearings and compare them with four different classifiers.

Wind turbine rolling bearing fault diagnosis empirical mode decomposition
Introduction

With the development of wind power generation, there will be more and more large-capacity large-scale units. Due to the geographical distribution of wind energy, wind turbines generally work in relatively harsh environments . The gearboxes of wind turbines are often affected by alternating loads, strong gusts, and large environmental temperature changes . Therefore, during the long-term operation of the wind turbine generator set, failures will inevitably occur. The failure of one piece of equipment may cause a series of reactions, resulting in an abnormal state of the entire system [3,4]. If it is possible to detect the operation of wind turbines, analyze and diagnose faults, and find the problems in time, it can avoid the faults in time.

Wind turbines are composed of blades, gearboxes, generators, frequency conversion systems, primary control, etc. Among them, the most prone to faults and the most affected parts of the fault are mostly in critical parts such as gearboxes and generators. However, the gearbox is the part with the highest failure rate. Common failures of gearboxes generally occur in gears and rolling bearings, and the failure rate of rolling bearings is very high.

Even if the rolling bearing is not much different from the factory conditions, some bearings fail before reaching the theoretical service life, but some bearings have exceeded the theoretical service life but are still intact. If the rolling bearings are regularly repaired, this will result in the inability to use the rolling bearings better. Therefore, the detection, diagnosis, and analysis of the bearing state are beneficial to accurately understand the working condition of the rolling bearing and avoid unnecessary losses caused by unknown accidents.

Rolling bearing detection and diagnosis technology have experienced more than sixty years of development . After Cooley and Tukey proposed the Fast Fourier Transform (FFT) algorithm theory in 1965 , the bearing diagnosis technology based on spectrum analysis has been developed rapidly . Developed by Svenska Kullager-Fabriken (SKF) in 1969, a bearing diagnostic instrument called shock pulse meter . Dr. Harting applied for a patent for a “resonance demodulation system” in 1974 . Later McFadden et al.  comprehensively summarized the application of this technology in bearing diagnosis. Many Chinese scholars are also committed to the fault diagnosis of wind power bearings. Literature  uses the Short-time Fourier Transform (STFT) and generative neural networks methods to diagnose rolling bearing faults. Literature  uses the Complementary Ensemble Empirical Mode Decomposition (CEEMD) and Linearly Decreasing Particle Swarm Optimization Probabilistic Neural Network (LDWPSO-PNN) methods to analyze and compare the vibration signals of rotating machinery. Literature  uses the Variational Modal Decomposition Fractional Fourier Transform (VMD-FRFT) method to diagnose rolling bearing faults.

This paper takes the rolling bearings as the research object and analyzes the states of rolling bearings. Make full use of signal processing methods such as EMD to research feature extraction and failure analysis. After the signal is obtained for preprocessing, the kurtosis criterion is used to select Intrinsic Mode Functions (IMF), and then envelope spectrum analysis and feature extraction are performed to classify faults. Finally, a rolling bearing fault analysis method based on EMD K-Nearest Neighbor (EMD-KNN) is given.

Rolling Bearing Vibration Mechanism and Characteristic Signal Frequency

The factors affecting the vibration of rolling bearings are internal and external. The vibration of rolling bearings can be divided into three categories: (1) Natural vibration; (2) Forced vibration caused by errors in the processing and assembly of parts of each part; (3) The outer and inner ring grooves or the surface of the ball have impact vibration caused by damages such as wear, scratches, pitting and spalling. Generally, the vibration of a rolling bearing is a superposition of the above three types of vibration, so its actual structure is very complicated . Therefore, to conduct an effective fault diagnosis, it is necessary to decompose the vibration signal, strip out irrelevant signals caused by noise, etc., and extract the vibration signal caused by the fault.

Here, the SKF6205 deep groove ball bearing is utilized. Acceleration sensors measure the vibration signals of rolling bearings in different states. This paper uses a drive end bearing with a diameter of 0.17781 mm and a speed of 1797 r/min; the frequency of bearing rotation under this condition is 29.95 Hz; the sampling frequency is 12 kHz.

In normal working mode, the natural frequency of rolling bearing is only related to the material and structure of the bearing itself. So the natural frequency of each part will be calculated by the following formula :

The natural frequencies of the outer and the inner rings are: fn=n(n21)2πn2+1×4D2EIgρA where: fn is natural frequency; n is vibration order; g is gravity acceleration, g=9800 mm/s2; ρ is material density; D is the diameter of the neutral shaft in the cross-section of the rolling bearing ring; I is the moment of inertia; E is modulus of elasticity, steel is 210 GPa; A is the cross-sectional area of the rolling bearing ring.

The natural frequency of the rolling element is: fbn=0.424rE2ρ where: fbn is natural frequency of rolling element; r is radius of rolling elements.

It should be noted that using Eqs. (1) and (2) can only obtain the natural frequency in the free state. Generally, the natural frequency value of rolling bearings is in the range of thousands to tens of thousands of hertz, and the size is relatively large. When the rolling bearing is running after a failure, different vibration frequencies will be generated when the fault point contacts and collides with other non-faulty components. This frequency can directly reflect the operating conditions of the rolling bearing. When the fault point is different, its value is also different.

The rotation speed of the inner ring is: Vi=2πr1fi=πfi(Ddcosα) where: Vi is rotation speed of inner ring; D is the diameter of the circle where the center of the rolling element is located; d is diameter of rolling element; r1 is average radius of inner ring raceway (2r1= Ddcosα); fi is inner ring fault characteristic frequency; α is the contact angle.

The rotation speed of the outer ring is: Vo=2πr2fo=πfo(D+dcosα) where: Vo is rotation speed of outer ring; fo is outer ring fault characteristic frequency; r2 is the flat radius of the outer ring raceway (2r2=D+dcosα).

The rotation speed of the cage is: Vc=12(Vo+Vi)=πfcD where: fc is the rotation frequency of cage itself.

From Eqs. (3) and (4), the rotation speed of the cage itself can be obtained. The rotation frequency of the cage itself is: fc=Vi+Vo2πD=12(1dDcosα)fi+12(1+dDcosα)fo

The rotation frequency of the cage relative to the inner ring and the rotation frequency relative to the outer ring are, respectively: fic=|fcfi|=12(fifc)(1+dDcosα) foc=|fcfo|=12(fifo)(1dDcosα) where: fic is the rotation frequency of cage relative to inner ring; foc is the rotation frequency of the cage relative to the outer ring. fr=fifo=n60 where: fr is the relative rotation frequency of inner and outer ring; n is rolling bearing speed.

When the inner ring of the rolling bearing fails, the calculation formula of the fault characteristic frequency is: fi=zfic=12z(1+dDcosα)fr where: z is the number of rolling elements.

When the outer ring of the rolling bearing fails, the calculation formula of the fault characteristic frequency is: fo=12z(1dDcosα)fr.

When the rolling element of a rolling bearing fails, the calculation formula of the fault characteristic frequency is: fb=2r1dfic=D2d[1(dD)2cos2α]fr where: fb is the characteristic frequency of rolling element failure.

Failure analysis analyzes the failure mechanism, type and impact, frequency of occurrence, and development and change laws of the researched object. Only on the basis of fault analysis can the appropriate diagnosis method be determined to carry out effective fault diagnosis. Therefore, fault analysis and fault diagnosis are generally developed at the same time.

With the progress of human society and economic development, the equipment used by people has become more sophisticated and complex . So early fault diagnosis methods such as sound and odor have not been able to meet actual needs . The modern development of sensing, measurement technology, and computer-based fault diagnosis technology are developing increasingly prosperous . With the in-depth research of many scholars and experts, many emerging fault diagnosis and analysis techniques have been continuously applied to practical engineering [10,20]. As large-scale and large-capacity mechanical equipment, wind turbines have a relatively complete fault diagnosis system. The fault analysis and diagnosis of wind turbines are mainly to detect the operating status of the gearbox. The rolling bearing is a crucial component in the gearbox, so it is necessary to carry out detection status and failure analysis and diagnosis.

The specific diagnosis process has four steps: (1) Signal acquisition: According to the structure and operating characteristics of the equipment, use suitable sensors and suitable measuring point positions to measure the signal of the working state of the rolling bearing. Finally, display and store it; (2) Signal preprocessing: Because the signal collected by the sensor will have certain errors, such as noise and interference, it is not easy to directly obtain the fault characteristics. Therefore, the signal must be preprocessed first. The preprocessing method in this paper is mainly to perform re-sampling; (3) Signal feature extraction: Extract the fault features through signal analysis methods on the preprocessed signal. In this paper, the EMD method is used for signal feature extraction. This way can further judge the working conditions of rolling bearings and enhance the fault characteristics; (4) State detection and decision-making: make decision analysis by analyzing the signal after feature extraction.

Fault Frequency Analysis Based on EMD

Decomposing the data used in this article into EMD will get several IMF components and then use the kurtosis principle to select the appropriate IMF components. First, the original signal of bearing vibration is decomposed by EMD to obtain several IMFs . Second, the kurtosis value is calculated. The IMF components with a kurtosis value greater than three are retained. Third, the selected IMF components are superimposed. Finally, the envelope spectrum of the signal is drawn.

EMD can decompose complex, nonlinear, and non-stationary signals to obtain multiple IMF through its characteristic time scale . Each IMF corresponds to its characteristic time scale.

As a method of processing signals, Hilbert transform is to convolve a given signal with 1πt. For continuous signal y(t), do Hilbert transform: x(t)=1π+y(τ)tτdτ where: y(τ) is signal y(t).

After Eq. (13), the Hilbert transform produces a new signal. The amplitude and frequency of this new signal remain constant while shifting the phase of the original signal by 90°. According to the new signal, the analytical signal of y(t) can be constructed: z(t)=y(t)+ix(t)=a(t)eiφ(t) φ(t)=arctan(x(t)y(t)),a(t)=x2(t)+y2(t) where: a(t) is the instantaneous amplitude of y(t); φ(t) is the instantaneous phase of y(t).

From this find the instantaneous frequency: f(t)=12πdφ(t)dt

After processing a(t), the envelope spectrum of y(t) can be obtained. The envelope spectrum represents the intensity of each frequency of the signal. For shock vibration, the envelope spectrum is highly sensitive. Therefore, the observation and analysis by the envelope spectrum can obtain more useful characteristic information.

The kurtosis K reflects the numerical statistics of the vibration signal distribution characteristics, expressed as [25,26]: K=E(xμ)4σ4 where: x is analyzed vibration signal; σ is the standard deviation of x; μ is the mean value of x.

In this paper, the kurtosis criterion is used to screen IMF, and IMF components with K > 3 are used for envelope spectrum analysis.

Analyze the bearing in a healthy state. The decomposed IMF components are used to calculate the kurtosis value using MATLAB software. According to the criterion that the kurtosis value is greater than 3, the 9 IMF components after decomposition are screened, so IMF2, IMF3, IMF6, IMF8 are selected. The screened IMF component diagram and envelope spectrum diagram are drawn. The IMF image is shown in Fig. 1.

IMF component image of healthy signal

The envelope spectrum is shown in Fig. 2.

In Fig. 2, the theoretical vibration frequency in a healthy state is 29.95 Hz. Some other peak points can also be seen on the graph, which are noise and interference.

Analyze the signal of the bearing inner ring failure. First, according to the criterion that the kurtosis value is greater than 3, the decomposed 10 IMF components are screened. Therefore, the selected components are IMF1, IMF2, IMF4, IMF5, IMF6, IMF8. Second, draw the filtered IMF component diagram and envelope spectrum diagram. The IMF image is shown in Fig. 3.

Envelope spectrum of healthy IMF components IMF component image of the inner ring fault state signal

The envelope spectrum is shown in Fig. 4.

Inner ring fault envelope spectrum

In Fig. 4, the theoretical vibration frequency in a healthy state is 161.83 Hz. Some other peak points can also be seen on the graph, which are noise and interference.

Analyze the original signal of the bearing rolling element failure. First, according to the criterion that the kurtosis value is greater than 3, the decomposed 10 IMF components are screened. Therefore, the selected components are IMF2, IMF3. Second, draw the filtered IMF component diagram and envelope spectrum diagram. The IMF image is shown in Fig. 5.

IMF component of rolling element fault state signal

The envelope spectrum is shown in Fig. 6.

Rolling element fault envelope spectrum

In Fig. 6, it is obvious that the healthy vibration frequency and its frequency multiplier can be seen, but the vibration frequency of the rolling element failure is not found. Therefore, it cannot be thoroughly analyzed and determined to be the rolling element failure.

Analyze the original signal of the bearing outer ring failure. First, according to the criterion that the kurtosis value is greater than 3, the decomposed 10 IMF components are screened. Therefore, the selected components are IMF1∼IMF9. Second, draw the filtered IMF component diagram and envelope spectrum diagram. The IMF image is shown in Fig. 7.

The kurtosis value of each IMF component

The envelope spectrum is shown in Fig. 8.

Outer ring fault envelope spectrum

In Fig. 8, it can be found that the fault vibration frequency of the envelope spectrum image of this state is the most noticeable image of the fault vibration frequency. The noise and interference also are less.

Fault Analysis of Rolling Bearing Based on EMD-KNN Method

The work of the previous chapter mainly focused on feature extraction, fault analysis, and diagnosis for a single signal. After EMD processing the data, this paper uses the envelope spectrum to analyze the failure of rolling bearings. Although the analysis results of the health status and the failure status of the outer ring are very obvious, the analysis results of the failure status of the rolling elements are disturbed. In the work of this chapter, the input signal is a complex mixed signal. The neural network classifier is added to find the frequency characteristics of each state from the complex signal. In this way, the correctness of the fault analysis and diagnosis can be ensured.

First, preprocess the original data and perform the same preliminary work as the previous chapter to obtain the IMF components selected by the kurtosis criterion. Second, superimpose the IMF components to calculate the time-frequency domain parameters of the superimposed signal. Third, normalize the computed data. Fourth, select four classifiers to train and classify the data. Finally, after comparison, a suitable classifier can be selected.

Signal Feature Extraction in Time-Frequency Domain

Since the previous work is similar to the previous work of Chapter 3, it will not be repeated here. The superimposed signal of the IMF component is obtained directly for calculation.

Signal Time Domain Characteristic Parameters

Peak difference (pk):  The peak difference is the difference between the maximum and the minimum value of the data.

pk=YmaxYmin where: Ymax is maximum signal; Ymin is signal minimum.

Standard deviation (δY(t)): It is also called the mean square error, which is the arithmetic square root of the variance.

δY(t)=limN1Ni=1N[Yi(t)μY(t)]2 where: Yi(t) is given signal; μY(t) is mean of the signal.

Kurtosis value (K):  kurtosis K reflects the numerical statistics of the vibration signal distribution characteristics of rolling bearings. The calculation formula is:

K=E(yμ)4σ4 where: y is given signal; μ mean of the signal; σ is signal variance; E is signal amplitude.

Average (μy(t)): The average is the average of the data. The calculation formula is:

μy(t)=limN1Nyi(t)=E[Y(t)] where: N is the number of data points of the signal.

Signal Frequency Domain Characteristic Parameters

Mean square frequency (msf): The frequency domain characteristic parameter index reflects the relationship between signal frequency and amplitude. The calculation formula is:

msf=14π2Δ20πω2S(ω)dω0πS(ω)dω where: S(ω) is power spectrum of signal y(t), S(ω)=Y(ω)Y¯(ω), Y(ω)=n=0N1y(t)ejnω; ω is the angular frequency of the signal data.

Root mean square frequency (rmsf): The calculation formula is:

rmsf=msf

Frequency variance (vf): To calculate the frequency variance, first calculate the center of gravity frequency (fc). The calculation formula is:

vf=14π2Δ20π(ω2πΔfc)2S(ω)dω0πS(ω)dω=msf(fc)2 where: fc=12πΔ0πωS(ω)dω0πS(ω)dω

Frequency standard deviation (rvf): The calculation formula is:

rvf=vf

Calculation Results of Characteristic Parameters

The calculated data is the superposition of the IMF components selected by the kurtosis criterion in Chapter 3. First, select one hundred thousand of the data points. Second, call the reshape function. Third, turn the hundred thousand data into a 5000 × 20 matrix. The matrix represents each state divided into 20 sets, and each set of data has 5000 data points. When repeating the above steps, four 5000 × 20 matrices of four states are obtained. According to Eqs. (18)(25), calculated in MATLAB software, eight time-frequency domain characteristic parameters will form a 20 × 1 matrix. This matrix is the parameter value of each group of data in each state.

Since each feature parameter has its dimension, it is not convenient to directly join the neural network. Therefore, this article uses the method of data normalization to normalize the data. The data are normalized as shown in Table 1.

Failure Analysis and Simulation of Rolling Bearings Vibration Signal Analysis Method of Rolling Bearing

Naive Bayes (NB) : This method is one of the classic learning algorithms. The idea of this algorithm is based on probability prediction. Its theoretical basis is conditional probability, Set of Word (SOW) model, and Bag of Words (BOW) model . The core of the method is Bayes’ rule, and the key of Bayes’ rule is conditional probability. SOW model is to count whether a certain term appears in the document in a given document. BOW model is used to calculate the frequency of a certain term in a given document. In addition, high-frequency words and stop words with very low importance need to be eliminated. Therefore, The BOW model is more accurate and effective.

Data normalization of time domain parameters
 Bearing status Peak difference Standard deviation Kurtosis mean Health 0.041158 0.062256 0.36854 0.095909 Outer ring failure 0.972764 0.970917 0.980174 0.417068 Inner ring failure 0.400276 0.403035 0.699864 −0.64849 Rolling element failure 0.027969 0.03545 0.490336 0.186732

K-Nearest Neighbor (KNN) : This algorithm has three basic elements: K value, distance measurement, and classification decision rules. If the K value is too small, the accuracy will be reduced. However, if it is too large, irrelevant data will be included, and the classification effect will be reduced. Therefore, a reasonable K value should be selected. The greater the number of variables, the worse the ability to distinguish Euclidean distance; the greater the range of the variable, the variable will play a leading role in the distance calculation. Therefore, the variables must be standardized first. The classification decision rules generally choose the weighted voting method.

Discriminant Analysis Classifier (DAC) : This method is an adaptive segmentation algorithm based on multi-feature discriminant analysis. It is generally used for image segmentation problems. In this method, in order to learn a linear transformation matrix, each input high-dimensional sample is converted into a low-dimensional vector through a linear transformation.

Simulation Results

This paper will divide the raw signal data of each of the four vibration states of rolling bearings into 20 groups with 5000 data in each group, a total of 80 groups of data. Among the 20 sets of data in each state, 15 sets are used to train the classifiers, a total of 60 sets; the remaining 5 sets of data in each state are tested to determine the fault diagnosis accuracy and effects of several classifiers.

Training and testing are implemented in MATLAB software. Use the tic toc function for timing and accuracy calculations. The four classifiers are operated separately. Finally, the classification accuracy of the KNN neural network is 100%, which takes 0.449198 s and the Random Forest classifier (RF) neural network takes 2.161731 s. The characteristic of the K-nearest neighbor classifier is to follow the data to judge completely, and there is no specific fixed mathematical model.

In this way, using the other three classifier operations, the final result is shown in Table 2.

Training and test results
 Neural networks Accuracy (%) Time (s) NB 100 0.630496 KNN 100 0.449198 RF 100 2.161731 DAC 100 0.574576

Since the test data is less and the characteristic parameters are more obvious, the accuracy of the four classifiers is 100%. However, it can be seen that the running time of each classifier is different. Among them, the longest time-consuming is RF, which takes 2.161731 s. The shortest time-consuming is KNN, which takes 0.449198 s.

On the premise of this sample of rolling bearing vibration signal data, the classification effect of the K nearest neighbor classifier is better than the other three classifiers. Therefore, the EMD-KNN method not only improves the accuracy of the typical fault diagnosis of rolling bearings but also only uses one operation to obtain the final result.

Conclusions

This paper takes the rolling bearing in the wind turbine generator as the research object, extracts the characteristics of its vibration signals in four different states, processes the signals based on EMD, and uses different methods for fault analysis. The research results of this paper are as follows: (1) The experimental results show that relative to directly observing the original signal, the signal characteristics after processing are more obvious, and it is easier to perform fault analysis and diagnosis. (2) Based on the EMD processed data, the time-frequency domain feature extraction is performed. Four neural network classifiers are added to obtain a more efficient and fast fault analysis method. It shows that fault analysis and diagnosis can be carried out more quickly and conveniently after joining the neural network.

Funding Statement: The support from the Guangdong Basic and Applied Basic Research Foundation (Grant No. 2021A1515012070), the Sichuan Science and Technology Program (Grant Nos. 2021YFS0336 and 2019YJ0712), the Fundamental Research Funds for the Central Universities (Grant No. ZYGX2019J035), and the Sichuan Science and Technology Innovation Seedling Project Funding Project (Grant No. 2020023) are gratefully acknowledged.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

References Liu, Z., Zhang, L. (2020). A review of failure modes, condition monitoring and fault diagnosis methods for large-scale wind turbine bearings. Measurement, 149, 107002. DOI 10.1016/j.measurement.2019.107002. Ilkılıç, C., Aydın, H., Behçet, R. (2011). The current status of wind energy in Turkey and in the world. Energy Policy, 39, 961967. DOI 10.1016/j.enpol.2010.11.021. Liu, S., Jiang, H., Wu, Z., Li, X. (2022). Data synthesis using deep feature enhanced generative adversarial networks for rolling bearing imbalanced fault diagnosis. Mechanical Systems and Signal Processing, 163, 108139. DOI 10.1016/j.ymssp.2021.108139. Tan, J., Fu, W., Wang, K., Xue, X., Hu, W. et al. (2019). Fault diagnosis for rolling bearing based on semi-supervised clustering and support vector data description with adaptive parameter optimization and improved decision strategy. Applied Sciences, 9, 1676. DOI 10.3390/app9081676. Tao, H., Wang, P., Chen, Y., Stojanovic, V., Yang, H. (2020). An unsupervised fault diagnosis method for rolling bearing using STFT and generative neural networks. Journal of the Franklin Institute, 357, 72867307. DOI 10.1016/j.jfranklin.2020.04.024. Huo, Z., Zhang, Y., Shu, L., Gallimore, M. (2019). A new bearing fault diagnosis method based on fine-to-coarse multiscale permutation entropy, laplacian score and SVM. IEEE Access, 7, 1705017066. DOI 10.1109/ACCESS.2019.2893497. Zhang, C., Yao, W., Deng, W. (2020). Fault diagnosis for rolling bearings using optimized variational mode decomposition and resonance demodulation. Entropy, 22, 739. DOI 10.3390/e22070739. McFadden, P. D., Smith, J. D. (1984). Vibration monitoring of rolling element bearings by the high-frequency resonance technique–A review. Tribology International, 17, 310. DOI 10.1016/0301-679X(84)90076-8. Liu, F., Gao, J., Liu, H. (2020). The feature extraction and diagnosis of rolling bearing based on CEEMD and LDWPSO-PNN. IEEE Access, 8, 1981019819. DOI 10.1109/Access.6287639. Li, X., Ma, Z., Kang, D., Li, X. (2020). Fault diagnosis for rolling bearing based on VMD-fRFT. Measurement, 155, 107554. DOI 10.1016/j.measurement.2020.107554. Chen, Y., Zhang, T., Luo, Z., Sun, K. (2019). A novel rolling bearing fault diagnosis and severity analysis method. Applied Sciences, 9, 2356. DOI 10.3390/app9112356. Liu, X., Liu, X., Zhou, Z., Hu, L. (2021). An efficient multi-objective optimization method based on the adaptive approximation model of the radial basis function. Structural and Multidisciplinary Optimization, 63, 13851403. DOI 10.1007/s00158-020-02766-2. Liu, X., Wang, X., Xie, J., Li, B. (2020). Construction of probability box model based on maximum entropy principle and corresponding hybrid reliability analysis approach. Structural and Multidisciplinary Optimization, 61(2), 599617. DOI 10.1007/s00158-019-02382-9. Zhou, Z., Chen, S., Liu, X. (2020). The design of linear magnetic negative stiffness element for engineering application using rectangular permanent magnets. Journal of Magnetics, 25(2), 172180. DOI 10.4283/JMAG.2020.25.2.172. Zhu, S. P., Keshtegar, B., Trung, N. T., Yaseen, Z. M., Bui, D. T. (2021). Reliability-based structural design optimization: Hybridized conjugate mean value approach. Engineering with Computers, 37(1), 381394. DOI 10.1007/s00366-019-00829-7. Zhu, S. P., Keshtegar, B., Bagheri, M., Hao, P., Trung, N. T. (2020). Novel hybrid robust method for uncertain reliability analysis using finite conjugate map. Computer Methods in Applied Mechanics and Engineering, 371, 113309. DOI 10.1016/j.cma.2020.113309. Zhu, S. P., Keshtegar, B., Tian, K., Trung, N. T. (2021). Optimization of load-carrying hierarchical stiffened shells: Comparative survey and applications of six hybrid heuristic models. Archives of Computational Methods in Engineering, 28, 41534166. DOI 10.1007/s11831-021-09528-3. Su, X., Li, L., Qian, H., Mahadevan, S., Deng, Y. (2019). A new rule to combine dependent bodies of evidence. Soft Computing, 23(20), 97939799. DOI 10.1007/s00500-019-03804-y. Su, X., Li, L., Shi, F., Qian, H. (2018). Research on the fusion of dependent evidence based on mutual information. IEEE Access, 6, 7183971845. DOI 10.1109/Access.6287639. Li, Y. H., Sheng, Z., Zhi, P., Li, D. (2019). Multi-objective optimization design of anti-rolling torsion bar based on modified NSGA-III algorithm. International Journal of Structural Integrity. 12(1), 1730. DOI 10.1108/IJSI-03-2019-0018. Nasir, N. N. M., Singh, S., Abdullah, S., Haris, S. M. (2019). Accelerating the fatigue analysis based on strain signal using hilbert–Huang transform. International Journal of Structural Integrity. 10(1), 118132. DOI 10.1108/IJSI-06-2018-0032. Theotokoglou, E. E., Balokas, G., Savvaki, E. K. (2019). Linear and nonlinear buckling analysis for the material design optimization of wind turbine blades. International Journal of Structural Integrity. 10(6), 749765. DOI 10.1108/IJSI-02-2018-0011. Yang, J., Huang, D., Zhou, D., Liu, H. (2020). Optimal IMF selection and unknown fault feature extraction for rolling bearings with different defect modes. Measurement, 157, 107660. DOI 10.1016/j.measurement.2020.107660. Huang, N. E., Shen, Z., Long, S. R., Wu, M. C., Shih, H. H. et al. (1998). The empirical mode decomposition and the hilbert spectrum for nonlinear and non-stationary time series analysis. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 454(1971), 903995. DOI 10.1098/rspa.1998.0193. Zhao, J., Zhang, Y., Chen, Q. (2020). Rolling bearing fault feature extraction based on adaptive tunable Q-factor wavelet transform and spectral kurtosis. Shock and Vibration, 2020, 119. DOI 10.1155/2020/8875179. Wan, S., Zhang, X., Dou, L. (2018). Compound fault diagnosis of bearings using an improved spectral kurtosis by MCDK. Mathematical Problems in Engineering, 2018, 112. DOI 10.1155/2018/6513045. Wang, J., Mo, Z., Zhang, H., Miao, Q. (2019). A deep learning method for bearing fault diagnosis based on time-frequency image. IEEE Access, 7, 4237342383. DOI 10.1109/Access.6287639. Yu, J., Ding, B., He, Y. (2018). Rolling bearing fault diagnosis based on mean multigranulation decision-theoretic rough set and non-naive Bayesian classifier. Journal of Mechanical Science and Technology, 32(11), 52015211. DOI 10.1007/s12206-018-1018-7. Zhang, N., Wu, L., Yang, J., Guan, Y. (2018). Naive Bayes bearing fault diagnosis based on enhanced independence of data. Sensors, 18(2), 463. DOI 10.3390/s18020463. Wang, H., Yu, Z., Guo, L. (2020). Real-time online fault diagnosis of rolling bearings based on KNN algorithm. Journal of Physics: Conference Series, 1486(3), 032019. DOI 10.1088/1742-6596/1486/3/032019. Zhou, Y., Yan, S., Ren, Y., Liu, S. (2021). Rolling bearing fault diagnosis using transient-extracting transform and linear discriminant analysis. Measurement, 178, 109298. DOI 10.1016/j.measurement.2021.109298.