# Vibration performance prediction and reliability analysis for rolling bearing

## Fannian Meng1, Xiaoyun Gong2, Wenliao Du3, Liangwen Wang4, Feng Zhao5, Liwei Li6

1, 2, 3, 4, 5, 6Henan Key Laboratory of Intelligent Manufacturing of Mechanical Equipment, Zhengzhou University of Light Industry, Zhengzhou, 450002, China

1Corresponding author

Journal of Vibroengineering, Vol. 23, Issue 2, 2021, p. 327-346. https://doi.org/10.21595/jve.2020.21463
Received 9 May 2020; received in revised form 8 October 2020; accepted 7 November 2020; published 9 January 2021

Copyright © 2021 Fannian Meng, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 111
Abstract.

The bearing vibration signal is a rich dynamic symptom of bearing wear, and the vibration signal of rolling bearing presents chaotic characteristics. Input and output variables of vibration signal can be constructed through phase space reconstruction, the Input and output variables can be imported into the prediction model for prediction. The prediction accuracy of the extreme learning machine (ELM) model, Kriging model and RBF model are compared, the results show that ELM has higher accuracy, so ELM chaos model is used to predict the future vibration time series data, and the forecasting error can be obtained by comparing the prediction value with the actual values so as to verity the feasibility of the ELM model. The prediction results of the future state of the bearing are processed as the grey-bootstrap method, and the performance reliability prediction of the bearing is realized by the Poisson counting process. The experimental data show that with the deepening of the fault degree, the reliability performance decreases gradually. The reliability performance of the bearing without fault is 100 %, and the reliability performance is 47.56 % when the inner ring faulty size is 0.72 mm.

Highlights
• Rolling bearing prediction sample Y generated by ELM model
• Large sample YB in statistical sense is generated from small sample Y based on grey bootstrap method
• Reliability R(θ) of bearing for each state is achieved according to the Poisson formula
• The experimental data show that with the deepening of the fault degree, the reliability performance decreases gradually

Keywords: rolling bearing, chaotic theory, grey-bootstrap method, reliability analysis, ELM.

#### 1. Introduction

Rolling bearings are an important support for shafts and other rotating components, and their performance is vital for the normal equipment operation. As the core component of mechanical equipment, the vibration performance of rolling bearings is especially applied in the status assessment of high-precision equipment such as satellites and space shuttles. Therefore, the prediction of the rolling bearing vibration provides an important guarantee for the safe equipment operation.

Various factors are intertwined in the rolling bearing to cause their vibration performance nonlinear, these factors include: morphological characteristics of the roller contact surface and the viscosity-temperature effect of the lubricant . The ununiformity of complex factors ultimately leads to the vibration performance extreme sensitivity. The method based on the chaotic phase space reconstruction theory allows accurate prediction of results of future rolling bearing vibration.

The phase space reconstruction made as the dynamic system method is the basis for analyzing chaotic time series. Since the 1980s, many domestic and foreign scholars have studied phase space reconstruction techniques. Among which the most widely used method is a delayed coordinate state space reconstruction theory proposed by Packard and Stewart . Embedding-dimension and Delay-time parameters are used in phase-space reconstruction. The determination of these two parameters will affect the quality of phase space. Therefore, determining its value has very important theoretical and practical significance . The commonly used methods to calculate the embedding dimension are saturation correlation dimension (GP), pseudo nearest neighbor method, Cao method, etc. [4-10]. Mutual information method, autocorrelation function method, and C-C method are used to calculate the delay time [11, 12].

Chaotic dynamics can be used to analyze time series with equal time intervals, the dynamic characteristics of chaotic systems can be indirectly obtained through the study of these time series, and these time series contain rich dynamic information. Extracting and using this information to study the characteristics of the system is one of the important aspects of chaos research. Many scholars have studied the chaotic dynamics of bearings, and the chaos performance of rolling bearing is studied in reference [13-16]. The authors  carried out a variable prediction research of rolling bearings. Professor Xia Xintao  used the Chaos-prediction method to predict bearing vibration time series, and the prediction error can be obtained through comparison and analysis between prediction value and true value.

The bearing reliability theory mainly involves the fatigue failure, and the establishment of models depends on the failure data greatly, which often ignores a large amount of evolutionary information of time series signals (such as vibration, temperature, friction torque, etc.) that promote the bearing degradation during its service life. By digging and extracting useful information of this type of time series, bearing performance prediction and reliability detection can be realized.

At present, bearing performance prediction and reliability research models are mostly based on the failure data or classical statistics, and there are relatively few research based on time series. In references [19, 20], it is recommended to apply a reliability prediction method based on the state information, the prediction model is used to accurately calculate the degradation index during the bearing service. The model establishment breaks through the limitations of traditional reliability large sample failure data. The authors  established an empirical probability density function for the friction torque parameters of aerospace bearings, the theory achieves accurate prediction for friction torque time series with the help of fuzzy sets. The authors  used nonlinear state estimation methods to effectively predict the temperature performance of gearbox bearings and accurately detect the operating state of gearboxes.

The prediction of bearing state is relatively single in the above research, and the prediction of bearing performance and reliability has not been carried out at the same time. In the process of bearing operation, an unexpected large vibration signal may appear, which may cause the peak value of time series, so monitoring and predicting the bearing performance is difficulty, the combination of performance prediction and reliability analysis has certain advantages, which can dig out the hidden information and discover the actual hidden danger in time. In the article, a chaotic prediction and reliability analysis of rolling bearing vibration time series is made. A chaos prediction model taken from the extreme learning machine (ELM) is constructed by solving parameters embedding dimension, delay time and using phase space reconstruction theory, and the chaotic prediction of rolling bearing vibration sequence is carried out. Secondly, a large number vibration variation sample is produced based on the grey-bootstrap-method. The Poisson counting theory is applied under a given threshold, and the corresponding variation intensity is obtained. According to the Poisson process, the predicted value of bearing reliability is obtained. The chaotic prediction and reliability analysis can be effectively applied to engineering practice.

#### 2.1. Forecasting theory

Suppose the rolling bearing vibration is given as $\mathbf{X}=\left\{e\left(i\right),i=1,2,\cdots ,T\right\}$, $T$ is the number of time series. The matrices $\mathbf{X}$ and $\mathbf{Y}$ can be obtained using the delay coordinate method and are shown as follows:

(1)
$\mathbf{X}=\left[\begin{array}{cccc}e\left(1\right)& e\left(1+\tau \right)& \cdots & e\left(1+\left(m-1\right)\tau \right)\\ e\left(2\right)& e\left(2+\tau \right)& \cdots & e\left(2+\left(m-1\right)\tau \right)\\ ⋮& ⋮& ⋮& ⋮\\ e\left(N\right)& e\left(N+\tau \right)& \cdots & e\left(N+\left(m-1\right)\tau \right)\end{array}\right],$
(2)
$\mathbf{Y}=\left[\begin{array}{c}e\left(2+\left(m-1\right)\tau \right)\\ e\left(3+\left(m-1\right)\tau \right)\\ ⋮\\ e\left(N+1+\left(m-1\right)\tau \right)\end{array}\right],$

where $N=T-1-\left(m-1\right)\tau$, $\tau$ is delay time which is obtained by the autocorrelation function method, $m$ is embedding-dimension which is calculated by the Cao-method.

According to the Takens theory, the reconstructed phase trajectory is dynamically equivalent to the original system in the homeomorphic sense. The rolling bearing vibration signal $\mathbf{Y}\left(i\right)$ can be predicted based on the point $\mathbf{X}\left(i\right)$ in the phase space, and the mapping function $f$ is shown as follows, from which the next data point of time series can be gotten:

(3)
$\mathbf{Y}\left(i\right)=f\left(e\left(i\right)\right).$

The own chaotic characteristics of the above prediction model, and the mapping function $f$ has a non-linear structure. Traditional statistical methods such as auto-regression, moving average, and ARIMA are not suitable for solving it. Machine learning model ELM can establish non-linear mapping, which can act as mapping function to do such a prediction work, the specific prediction principle is shown in Fig. 1.

Fig. 1. Chaos prediction theory #### 2.2. Delay time calculation theory

The autocorrelation function method is used to calculate a delay time parameter, which is constructed based on the linear correlation degree of two motion trajectory parameters at time $t$ and $t+\tau$. The autocorrelation function $C\left(\tau \right)$ of vibration series is shown as follows:

(4)
$C\left(\tau \right)=\frac{\sum _{i=1}^{n-1}\left({x}_{i}-\overline{x}\right)\left({x}_{i+\tau }-\overline{x}\right)}{{\sum _{i=1}^{n}\left({x}_{i}-\overline{x}\right)}^{2}},$

where $\overline{x}$ is the mean value of the sample.

If the autocorrelation function shows significant attenuation trend with the growth of delay-time, the optimal delay time is equal to the time when the autocorrelation function reaches to $1-1/e$ times the initial value for the first time.

#### 2.3. Cao method for embedding dimension

Cao Liangyue et al proposed the Cao method, which has the advantages of not-relying on the subjectivity and high calculation efficiency, let:

(5)
$a\left(i,m\right)=\frac{‖{x}_{i}\left(m+1\right)-{x}_{n\left(i,m\right)}\left(m+1\right)‖}{‖{x}_{i}\left(m\right)-{x}_{n\left(i,m\right)}\left(m\right)‖},$

where, $‖•‖$ is phasor norm, ${x}_{i}\left(m+1\right)$ is the $i$th quantity of reconstructed phase space, where the embedding dimension is $m+1$, ${x}_{n\left(i,m\right)}$ is an integer greater than 1 and less than or equal to $N-m$, and the sequence with the shortest distance from ${x}_{i}\left(m+1\right)$ is ${x}_{n\left(i,m\right)}\left(m+1\right)$.

The mean value of all $a\left(i,m\right)$ is calculated using the above Eq. (5):

(6)
$E\left(m\right)=\frac{1}{N-m}\sum _{i=1}^{N-m}a\left(i,m\right).$

The embedding dimension is determined by observing the change of mean $E\left(m\right)$, along with $m$, the change ratio is shown as follows:

(7)
${E}_{1}\left(m\right)=\frac{E\left(m+1\right)}{E\left(m\right)}.$

When the change is stable, and the value of $m$ plus 1 at this time is the required embedding dimension.

The criteria for determining the stable change of ${E}_{1}\left(m\right)$ are given below, the specific assessment process is as follows:

(1) ${\mathrm{\Delta }}_{i}$ is calculated as:

(8)

(2) An initial threshold $e$ is chosen according to the fluctuation of ${E}_{1}$, here, $e=\overline{{\mathrm{\Delta }}_{i}}$, $\overline{{\mathrm{\Delta }}_{i}}$ represents the average value of ${\mathrm{\Delta }}_{i}$, and the subscript $u$ of the first ${\mathrm{\Delta }}_{i} is found, the value of ${\mathrm{\Delta }}_{j}$ is calculated:

(9)

(3) The values of $e$, $i$ are reset as: $e=\overline{{\mathrm{\Delta }}_{i}}$, $j\le i\le N-1$, where $\overline{{\mathrm{\Delta }}_{i}}$ represents the average value of ${\mathrm{\Delta }}_{i}$

(4) $j\le i\le N-2$ is taken when ${\mathrm{\Delta }}_{i}>{\mathrm{\Delta }}_{i+1}$ and ${\mathrm{\Delta }}_{i+1}>{\mathrm{\Delta }}_{i+2}$ and ${\mathrm{\Delta }}_{i+1}, embedding dimension is equal to $m=i+1$.

Through the above criteria, a relatively objective evaluation can be made when ${E}_{1}\left(m\right)$ tends to be stable, the embedding dimension calculation process is more scientific and rigorous.

#### 2.4. ELM theory

ELM model is an algorithm proposed by Huang el al, which has its own advantages for training a single hidden layer feed-forward neural network. Its training model is shown in Fig. 2.

In Fig. 2, $n$ is the number of input-variables, $h$ shows the number of hidden layer neurons, ${x}_{i}$ represents the input variable, $i=1,2,3,\dots ,n$, ${w}_{ij}$ represent the connection weight between the input variable and hidden layer neuron $j$, ${w}_{jo}$ represent the connection weight between output variable and hidden layer neuron $j$, ${b}_{j}$ represent the threshold of hidden layer neuron, $j=1,2,3,\dots ,h$, ${b}_{o}$ represents the threshold of neuron in the output layer, ${f}_{j}\left(•\right)$ represents the activation function of layered neurons, $g\left(•\right)$ represents the activation function of output layer neurons.

Fig. 2. ELM network structure Its mathematical model is:

(10)
$\left\{\begin{array}{l}y=g\left({b}_{o}+\sum _{j=1}^{h}{w}_{jo}{v}_{j}\right),\\ {v}_{j}={f}_{j}\left({b}_{j}+\sum _{j=1}^{h}{w}_{ij}{x}_{j}\right),\end{array}\right\$

where, ${v}_{j}$ represents the output of output layer neurons $j$.

If there are $N$ valid samples, when the output threshold is equal to 0, and the output neurons activation function is a linear activation function, then the Eq. (10) can be written as:

(11)
$y=\left({w}_{o}^{T}v{\right)}^{T},$

where, $y=\left[y\left(1\right),y\left(2\right),\cdots ,y\left(N\right){\right]}^{T}$ is the network output vector, ${w}_{o}=\left[{w}_{1o},{w}_{2o},\cdots ,{w}_{ho}{\right]}^{T}$ represents the output weight vector, $v$ is the output matrix of hidden neurons, input weights and threshold matrix $w$ are randomly generated:

(12)

The expression of output weight vector is:

(13)
${w}_{o}={v}^{*}{y}_{e},$

where, ${v}^{*}$ represents Moore-Penrose generalized inverse of output layer matrix $v\text{,}$${y}_{e}=\left[{y}_{e}\left(1\right),{y}_{e}\left(2\right),\cdots ,{y}_{e}\left(N\right){\right]}^{T}$ is desired output. If $v\in {R}^{N×h}$, $N\ge h$ and $rank\left(v\right)=h$, then the Moore-Penrose generalized inverse ${v}^{*}$ can be expressed as:

(14)
${v}^{*}=\left({v}^{T}v{\right)}^{-1}{v}^{T}.$

Putting Eq. (14) into Eq. (13), one can get:

(15)
${w}_{o}=\left({v}^{T}v{\right)}^{-1}{v}^{T}{y}_{e}.$

The used hidden layer activation functions are Sigmoid function, Sin function, RBF function and Hardlim function.

#### 2.5. Kriging model

The Kriging model can be expressed approximately as a sum of one random distribution function and one polynomial, as shown in Eq. (16):

(16)
$y\left(x\right)=f\left(x\right)\beta +z\left(x\right),$

where $y\left(x\right)$ is an unknown Kriging model, $y\left(x\right)$ is known as the two order regression function of $x$, the global approximation model in the design space is provided, $\beta$ is the undetermined coefficient of regression function, its value can be estimated by the known response value, $z\left(x\right)$ is a stochastic process, it is a local deviation on the basis of global simulation, and the expectation is 0, the variance is ${\sigma }^{2}$, the covariance matrix can be expressed as:

(17)
$cov\left[z\left({x}^{i}\right),z\left({x}^{j}\right)\right]={\sigma }^{2}R\left[R\left({x}^{i},{x}^{j}\right)\right],$

where R is the correlation matrix, $R\left({x}^{i},{x}^{j}\right)$ represents the correlation function of any two sample points $i$, $j=1,2,\dots ,n$, $n$ is the number of data in the sample. $R\left({x}^{i},{x}^{j}\right)$ has a variety of functional forms can be selected, common correlation functions include cubic function, gauss function, linear function, spherical function, spline function and so on.

#### 2.6. Radial basis function

Radial basis function (RBF for short) model is a kind of function which takes the Euclidean distance as the independent variable between the test point and the sample point.

RBF is one of the commonly used surrogate model, the basic form is as follows:

(18)
${f}_{r}\left(x\right)=\sum _{i=1}^{{n}_{s}}\left({\beta }_{r}\right)\varphi \left(‖x-{x}_{i}‖\right)={{\beta }_{r}}^{T}\phi .$

Basic function $\phi ={\left(\varphi \left(‖x-{x}_{1}‖\right),\dots ,\varphi \left(‖x-{x}_{ns}‖\right)\right)}^{T}\text{,}$ weight coefficient ${\beta }_{r}={\left(\left({\beta }_{r}{\right)}_{1},\dots ,\left({\beta }_{r}{\right)}_{ns}\right)}^{T}$, and ${\beta }_{r}$ should meet the interpolation conditions:

(19)
$\left({f}_{r}{\right)}_{i}={y}_{i},\begin{array}{ll}& i=1,2,\dots ,{n}_{s}\end{array},$

where ${y}_{i}$ is the true value, $\left({f}_{r}{\right)}_{i}$ is the prediction value, ${n}_{s}$ is the number of samples.

With:

(20)
${A}_{r}{\beta }_{r}=y,$
(21)
${\beta }_{r}={{A}_{r}}^{-1}y,$
(22)

where $\varphi$ is radial function, common radial functions include Gaussian function, Multi-quadrics function, Reciprocal Multi-quadrics function, Thin-Plate spline function.

#### 2.7. Accuracy evaluation method

Mean Squared Error (MSE), Root Mean Squared Error (RMSE) and Mean Absolute Error (MAE) indicators are used to calculate the prediction accuracy of the model, the expressions are:

(23)
$MSE=\frac{1}{N}\sum _{i=1}^{N}\left({y}_{i}-\stackrel{^}{{y}_{i}}{\right)}^{2},$
(24)
$RMSE=\sqrt{\frac{1}{N}\sum _{i=1}^{N}\left({y}_{i}-\stackrel{^}{{y}_{i}}{\right)}^{2}},$
(25)
$MAE=\frac{1}{N}\sum _{i=1}^{N}\left|{y}_{i}-\stackrel{^}{{y}_{i}}\right|,$

where, ${y}_{i}$ represent the true value, $\stackrel{^}{{y}_{i}}$ represents the prediction value of the machine learning model, $N$ represents the number of test samples.

#### 3.1. Grey-bootstrap method

According to the above prediction model, the vibration state information of bearing at the next $\mathbf{U}$ step can be predicted, with the vector $\mathbf{Y}$ expressed as follows:

(26)
$\mathbf{Y}=\left(y\left(1\right),y\left(2\right),\cdots ,y\left(u\right),\cdots ,y\left(U\right)\right),$

where, $\mathbf{Y}$ is the chaotic prediction data using the above ELM model, $y\left(u\right)$ is the $u$th data of $\mathbf{Y}$ data series, $u=1,2,\dots ,U$.

Using the above method, a number is randomly selected from $\mathbf{Y}$ using equal probability, the total selection number is equal to $q$, so the bootstrap sample ${\mathbf{V}}_{1}$ can be gotten, the above selection method is repeated $\mathbf{B}$ times, and $\mathbf{B}$ samples are obtained, the bootstrap sample is shown as follows:

(27)
${\mathbf{V}}_{Bootstrap}=\left({\mathbf{V}}_{1},{\mathbf{V}}_{2},\cdots ,{\mathbf{V}}_{b},\cdots ,{\mathbf{V}}_{B}\right),$

where, ${\mathbf{V}}_{b}$ represents the $b$th bootstrap-sample, $\mathbf{B}$ represents the total bootstrap-resample number, with:

(28)
${\mathbf{V}}_{b}=\left[{v}_{b}\left(g\right)\right],$

where, $g=1,2,\dots ,q$ and $b=1,2,\dots ,B$.

Based on the grey prediction model GM(1,1), supposing that the first-order accumulated generating operator for ${\mathbf{V}}_{b}$ is as shown below:

(29)
${Y}_{b}=\left[{y}_{b}\left(u\right)\right]=\sum _{j=1}^{g}{v}_{b}\left(j\right).$

The grey generated model can be described as the following differential equation:

(30)
$\frac{d{y}_{b}\left(\mathbf{u}\right)}{d\mathbf{u}}+{c}_{1}{y}_{b}\left(u\right)={c}_{2},$

where $\mathbf{u}$ is time variable, and ${c}_{1}$ and ${c}_{2}$ are the undetermined coefficients.

Use the increment instead of differential, the above formula can be expressed as:

(31)
$\frac{d{y}_{b}\left(\mathbf{u}\right)}{d\mathbf{u}}=\frac{\mathrm{\Delta }{y}_{b}\left(\mathbf{u}\right)}{\mathrm{\Delta }\mathbf{u}}={y}_{b}\left(\mathbf{u}+1\right)-{y}_{b}\left(\mathbf{u}\right)={v}_{b}\left(\mathbf{u}+1\right),$

where, $\mathrm{\Delta }\mathbf{u}$ is equal to the unit time interval. Furthermore, the mean generated sequence vector is set as below:

(32)
${\mathbf{Z}}_{b}=\left[{z}_{b}\left(u\right)\right]=\left[0.5{y}_{b}\left(u\right)+0.5{y}_{b}\left(u-1\right)\right].$

Under the initial condition ${y}_{b}\left(1\right)={v}_{b}\left(1\right)$, the least-square-solution of the grey differential equation is shown as below:

(33)
${\stackrel{^}{y}}_{b}\left(q+1\right)=\left({v}_{b}\left(1\right)-{c}_{2}/{c}_{1}\right){e}^{-{c}_{1}q}+{c}_{2}/{c}_{1},$

where the determined coefficients ${c}_{1}$ and ${c}_{2}$ are shown:

(34)
$\left({c}_{1},{c}_{2}{\right)}^{T}=\left({\mathbf{D}}^{T}\mathbf{D}{\right)}^{-1}{\mathbf{D}}^{T}\left({\mathbf{V}}_{b}{\right)}^{T},$

with:

(35)
$\mathbf{D}=\left(-{\mathbf{Z}}_{b},\mathbf{I}{\right)}^{T},$
(36)
$\mathbf{I}=\left(1,1,\cdots ,1\right).$

According to the inverse-accumulated generating, the $b$th generated data are expressed as follows:

(37)
${\stackrel{^}{v}}_{b}\left(q+1\right)={\stackrel{^}{y}}_{b}\left(q+1\right)-{\stackrel{^}{y}}_{b}\left(q\right).$

Therefore, $\mathbf{B}$ generated data for rolling bearing vibration sample can be expressed as below:

(38)
$\begin{array}{c}{\mathbf{Y}}_{B}=\left({w}_{1},{w}_{2},\cdots ,{w}_{b},\cdots ,{w}_{B}\right)=\left({\stackrel{^}{v}}_{1}\left(q+1\right),{\stackrel{^}{v}}_{2}\left(q+1\right),\cdots ,{\stackrel{^}{v}}_{b}\left(q+1\right),\cdots ,{\stackrel{^}{v}}_{B}\left(q+1\right)\right),\end{array}$

where ${w}_{b}$ is the $b$th produced data.

#### 3.2.1. Counting process

Supposing that the generated-sequence ${\mathbf{Y}}_{\mathbf{B}}$ (in Eq. (38)) for bearing’s future vibration signal has $\mu$ data over the vibration threshold $h$, namely, there are $\mu$ data falling outside the interval $\left[-h,h\right]$ of the best vibration performance, the variation intensity estimated value $\theta$ for ${\mathbf{Y}}_{\mathbf{B}}$ is shown as follows:

(39)
$\theta =\frac{\mu }{B}.$

Variation intensity refers to the frequency of the vibration amplitude for rolling bearing exceeding the optimal vibration performance interval, which is a vital characteristic parameter that affects the vibration variation process of the bearing operation, and it changes along with the different accuracy-threshold.

#### 3.2.2. Rolling bearing reliability dynamic prediction

The Poisson counting process can be expressed as follows:

(40)
$Q=\mathrm{e}\mathrm{x}\mathrm{p}\left(-\theta i\right)\frac{\left(\theta i{\right)}^{e}}{e!},$

where $i$ stands for the time variable with $i=1,2,3,\dots ,i\ge 1$, $\theta$ is the variation-intensity, $e$ is the failure event occurring number with $e=0,1,2,3,\dots$, and $Q$ is the probability of failure events occurring $e$ times. The reliability $R$ for occurring failure events can be obtained using the Poisson Counting process.

When solving the vibration performance reliability $R$, the probability is equal to 0 when the product do not have a vibration failure, that is $e=0$. And $i=1$ is the vibration performance reliability for the current time, namely the possibility of the vibration signal occurrence on the current generated sequence ${\mathbf{Y}}_{\mathbf{B}}$ is within the optimal vibration interval $\left[-h,h\right]$. According to the above Eq. (40), the reliability can be described as:

(41)
$R\left(\theta \right)=\mathrm{e}\mathrm{x}\mathrm{p}\left(-\theta \right),$

where $R\left(\theta \right)$ represents the probability for the rolling bearing to have within the optimal vibration during operation.

#### 3.3. Modeling basic ideas

The theoretical modeling uses a variety of mathematical models such as ELM chaotic-prediction-model, grey-bootstrap method and Poisson counting process. The modeling method is shown in Fig. 3. The concrete steps are as follows:

Step 1. Based on the rolling bearing time series $\mathbf{X}$ of the vibration performance, the embedding- dimension $m$ is obtained by the Cao-method, and the delay-time $\tau$ is obtained by the autocorrelation-method, so the phase-space-reconstruction can be realized based on the chaos theory.

Step 2. 30 step vibration data are predicted using the ELM chaotic prediction model, so the small sample $\mathbf{Y}$ can be constructed with the sample number of $\mathbf{U}$.

Step 3. A statistically large sample ${\mathbf{Y}}_{\mathbf{B}}$is generated from a small sample $\mathbf{Y}$ based on the grey-bootstrap method, so the Poisson count and bearing vibration variation strength can be easily gotten.

Step 4. Under the given vibration threshold $h$, the number $\mu$ is found out beyond the optimal vibration interval $\left[-h,h\right]$ from the large sample ${\mathbf{Y}}_{\mathbf{B}}$, and then the future $\mathbf{U}$-step variation intensity $\theta$ of the bearing is acquired, then, the reliability $R\left(\theta \right)$ of the bearing for each state in the future is obtained according to the Poisson counting formula.

Fig. 3. Vibration prediction and reliability analysis #### 4.1. Vibration test data

The data used in this paper for experimental validation were provided by the Case Western Reserve University (CWRU), Cleveland, Ohio, USA . As shown in Fig. 4, the test stand consists of an electronic motor, torque transducer, dynamometer, and control electronics. The bearing to be tested supports the motor shaft. The driving end bearing is SKF6205, and the fan end bearing is SKF6203. Acceleration sensors are placed above the bearing seat of the fan and the driving end of the motor to collect the vibration acceleration signal of the bearing.

Fig. 4. Vibration experiment platform of rolling bearing Fig. 5. Layout of rolling bearing simulation experiment Vibration signal is collected by a 16-channel data logger, the power and speed are measured by a torque sensor. The test data are obtained by an acceleration sensor on the bearing at the drive end, and the rotation speed is 1796 rpm. Fig. 5 depicts a rolling bearing simulation experiment layout.

The above experimental platform is used to collect the rolling bearing vibration signal. The dynamic time data signals are collected as shown in Fig. 6 to 10. Fig. 6 to Fig. 10 respectively depicts the vibration data under normal bearing and four inner ring fault size.

Fig. 6. Vibration data (normal bearing) Fig. 7. Vibration data (0.18 mm inner ring fault size) Fig. 8. Vibration data (0.36 mm inner ring fault size) It can be seen from Figs. 6-10 that with the deepening of bearing failure, the overall vibration data of rolling bearings show an increasing trend. The vibration data of normal bearings is generally distributed between [–0.2, 0.2], and only a few data exceed this range. The bearing vibration data of the inner ring fault size of 0.18 mm is generally between [–1, 1], and data graph shows a similar saw-tooth shape. The bearing vibration data of the inner ring fault size of 0.36 mm is generally at [–1, 1]. these data near 500 series are significantly larger than other data. The bearing vibration data of the inner ring fault size of 0.54 mm is generally between [–2, 2], and the graph looks like a series of oval-shaped components. The bearing vibration data of the inner ring fault size of 0.72 mm is generally between [–4, 4], and the graph looks like a random distribution.

Fig. 9. Vibration data (0.54 mm inner ring fault size) Fig. 10. Vibration data (0.72 mm inner ring fault size) #### 4.2. Chaotic prediction with time series

The above mentioned bearing vibration time series are identified as ${\mathbf{X}}_{1}$ (Normal bearing vibration data), ${\mathbf{X}}_{2}$(0.18 mm inner ring fault size), ${\mathbf{X}}_{3}$(0.36 mm inner ring fault size), ${\mathbf{X}}_{4}$ (0.54 mm inner ring fault size), ${\mathbf{X}}_{5}$ (0.72 mm inner ring fault size). The chaotic prediction method is used to predict these sequences, and the prediction step is 30, then the 2001-2030 original data are used to verify the accuracy and feasibility of five sequence prediction models.

Phase space parameter calculation: The auto-correlation method and Cao-method respectively have been used to acquire the delay-time and embedding-dimensions, and their results are shown in Table 1.

Table 1. Phase space parameters of 5 series

 Phase space parameter ${\mathbf{X}}_{1}$ ${\mathbf{X}}_{2}$ ${\mathbf{X}}_{3}$ ${\mathbf{X}}_{4}$ ${\mathbf{X}}_{5}$ Delay time 4 1 1 1 1 Embedding dimension 7 13 17 20 20

In order to compare the prediction accuracy of the three methods mentioned above, ELM model, Kriging model and RBF model are used to predict the next 30-steps of bearing vibration data, the comparison result between 30-steps chaotic prediction value and experimental value of above three methods are shown in Fig. 11. Among them, the red line means the experimental data, blue line means the ELM model predicted value, pink line means the Kriging model predicted value, and green line means the RBF model predicted value. The RMSE values for prediction error of three models mentioned above are shown as Fig. 12.

Fig. 11. Prediction results comparison of the three methods Fig. 12. Prediction accuracy evaluation of the three methods (RMSE) From Fig. 11 and Fig. 12, the message can be acquired that ELM model has the highest prediction accuracy compared with Kriging model and RBF model, so the ELM model is used to predict vibration signal for reliability analysis.

The determination of the phase-space parameters of the rolling bearing time series is the basis for its phase-space reconstruction, and it is required for preparing the chaotic prediction model. The comparison results between 30-step chaotic prediction value and experimental value for sequences data ${\mathbf{X}}_{1}$, ${\mathbf{X}}_{2}$, ${\mathbf{X}}_{3}$, ${\mathbf{X}}_{4}$, ${\mathbf{X}}_{5}$ are shown in Figs. 13-17. Among them, the red line means the experimental data, and blue line means the predicted data.

From Fig. 13, the shape of the 30-step prediction result of the sequence ${\mathbf{X}}_{1}$ is similar to the curve of the experimental result, the points that differ greatly from the original data are shown in step 1 and steps 19-21, where the maximum difference is only 0.0568 V.

From Fig. 14, the prediction results of the first 5 steps of the sequence ${\mathbf{X}}_{2}$ are almost consistent with the original data. The steps 6, 13, 20, 22, 26, 28, 29 are significantly different from the original data, and the largest step is step 22 where the value is equal to 0.4152 V, the prediction results of other steps are almost consistent with the original data.

Fig. 13. Vibration prediction data (${\mathbf{X}}_{1}$) Fig. 14. Vibration prediction data (${\mathbf{X}}_{2}$) Fig. 15. Vibration prediction data (${\mathbf{X}}_{3}$) Fig. 16. Vibration prediction data (${\mathbf{X}}_{4}$) From Fig. 15, the prediction result of the sequence ${\mathbf{X}}_{3}$ is relatively different from the original data, the largest step is step 10 equal to 0.2904 V, and the prediction error between steps 2, 13, 19, 26, 28, 29 and original data is relatively small.

From Fig. 16, the prediction results of the sequence ${\mathbf{X}}_{4}$ are relatively similar to the original data. The fitting degree between prediction result line and the experimental result line at steps 1-12 and 19-30 are higher, only the prediction results of steps 12-18 are significantly different from the original data.

From Fig. 17, the prediction results of the sequence ${\mathbf{X}}_{5}$ are similar to the original data, and the prediction result lines of steps 6-15 are almost identical to the original data, the change trend of the prediction lines of steps 15 to 30 is very similar to the original data, with a jagged jump between –2 and 2.

Fig. 17. Vibration prediction data (${\mathbf{X}}_{5}$) In order to scientifically evaluate the prediction accuracy, the RMSE value is used to calculate the prediction accuracy, as shown in Fig. 18.

Fig. 18. Prediction accuracy evaluation (RMSE) From Fig. 18, it is known that the RMSE evaluation value is small in each bearing state. When the inner ring fault size is 0.72 m, the RMSE value is approximately equal to 0.7 at the most.

In the 5 time series prediction of ELM chaotic prediction model, the difference between the predicted value and the experimental result is very small, and the two values maintain good consistency, that indicates that the prediction model is reliable.

#### 4.3. Forecast results generated by grey-bootstrap method

Now, the prediction results of the next 30 steps of each sequence are processed by grey-bootstrap method to simulate the large-scale generated data of the vibration performance for each state of the bearing in the next 30 steps. The generated data is doing Poisson count under the given threshold value, and then find the performance reliability of the bearing in the next 30 steps. In the grey-bootstrap generation, setting the sampling number 30 and repeated executions 10000, A large amount of generated-data ${\mathbf{Y}}_{\mathbf{B}}$ in the next 30 steps for sequence ${\mathbf{X}}_{1}$, ${\mathbf{X}}_{2}$, ${\mathbf{X}}_{3}$, ${\mathbf{X}}_{4}$, ${\mathbf{X}}_{5}$ is shown in Figs. 19-23.

Fig. 19. Generated data of series ${\mathbf{X}}_{1}$ in future 30 steps Fig. 20. Generated data of series ${\mathbf{X}}_{2}$in future 30 steps Fig. 21. Generated data of series ${\mathbf{X}}_{3}$in future 30 steps From the Figs. 19-23, the 10000 generated data simulated by grey-bootstrap-method has obvious randomness, and the corresponding value for sequences ${\mathbf{X}}_{1}$, ${\mathbf{X}}_{2}$, ${\mathbf{X}}_{3}$, ${\mathbf{X}}_{4}$ and ${\mathbf{X}}_{5}$ are concentrated in [–0.02, 0.06], [–0.2, 0.2], [–0.15, 0.3], [–0.4, 0.6] and [–0.5, 1]. The five sets of data show a non-linear increasing trend, which demonstrates that when the degree of failure deepens, the bearing vibration signal value gradually increases, and the inherent change law of the bearing vibration is closely related to the service life of bearing. Then, the evolution mechanism of performance reliability is determined and influenced.

Fig. 22. Generated data of series ${\mathbf{X}}_{4}$in future 30 steps Fig. 23. Generated data of series ${\mathbf{X}}_{5}$ in future 30 steps #### 4.4. Future state reliability assessment

The bearing vibration threshold 0.1 is set according to the generated data ${\mathbf{Y}}_{\mathbf{B}}$ in Figs. 19-23, and the number of generated data is to be found for each sequence exceeding the vibration threshold $h$, that is, the number $\mu$ that falls outside the optimal vibration interval [–0.1, 0.1] for the 10000 generated data of each sequence which can be calculated, so the variation intensity $\theta$ can also be obtained from the Eq. (39), then the reliability dynamic prediction result of the future 30 steps for each sequence can be calculated from Eq. (41). The results are shown in Table 2.

Table 2. Number of $\mu$ exceeding vibration threshold $h$, variation intensity $\theta$ and performance reliability $R$

 Reliability parameter ${\mathbf{X}}_{1}$ ${\mathbf{X}}_{2}$ ${\mathbf{X}}_{3}$ ${\mathbf{X}}_{4}$ ${\mathbf{X}}_{5}$ $M$ 0 1768 3414 5526 7432 $\theta$ 0 0.177 0.341 0.553 0.7432 $R$ / % 100 83.78 71.11 57.52 47.56

As it can be seen from Table 2, the data generated in the next 30 steps fall in the best vibration range for sequence ${\mathbf{X}}_{1}$, the variation intensity is 0, the reliability reaches 100 % that indicates that the bearing running state is relatively stable, no trace of bad behavior variation has occurred, and the state of maintaining the best vibration performance is very good.

In the sequence ${\mathbf{X}}_{2}$, the number of generated data in the next 30 steps that falls within the prescribed optimal vibration interval is 1768, the variation intensity is small and equal to 0.1768, and the reliability reaches 83.78 %, indicating that the bearing has a certain variation.

For the sequence ${\mathbf{X}}_{3}$, the number of generated data in the next 30 steps that falls within the prescribed optimal vibration interval is 3414, the variation intensity is relatively small and equal to 0.341, and the reliability is 71.11 %.

For the sequence ${\mathbf{X}}_{4}$, the number of generated data in the next 30 steps that falls within the specified optimal vibration interval is 5526, the variation intensity is more and equal to 0.553, and the reliability is 57.52 %.

For the sequence ${\mathbf{X}}_{5}$, the number of generated data in the next 30 steps that falls within the prescribed optimal vibration interval is 7432, the variation intensity is large and equal to 0.743, and the reliability is 47.56 %.

In order to describe the change of reliability intuitively, the change curve of the reliability parameter is shown in Fig. 24.

Fig. 24. Change of reliability parameter $R$ As shown in Fig. 24, as the degree of bearing failure deepens, the reliability gradually decreases from 100 % with zero-failure to 47.56 % when the inner ring failure size is 0.7112 mm.

As shown above, in the process of chaos prediction of bearing vibration time series, the prediction method of chaos prediction model based on ELM is accurate and reliable that meets the general prediction requirements of engineering practice. Combined with the grey-bootstrap method, the prediction value of each state sequence data in the next 30 steps is sampled and processed, and a large number of generated signals participated in bearing performance degradation is accurately simulated. The variation intensity obtained by counting process can effectively describe the variation degree of bearing, and can reveal the influence mechanism on the bearing state on the variation process of performance reliability in operation.

#### 5. Conclusions

1) ELM model has the highest prediction accuracy compared with Kriging model and RBF model, so the ELM model is used to predict vibration signal for reliability analysis.

2) The chaotic prediction model based on the ELM is accurate and reliable. It can accurately predict the performance value of the future state of the bearing vibration time series. The predicted value and the true value can maintain a good consistency, which can be better applied to the engineering prediction.

3) The grey-bootstrap principle is integrated into the Poisson process, and a reliability prediction method is proposed based on the vibration time series, and it can ensure the performance reliability prediction of the future state of the bearing, and the variation characteristics of the variation intensity effectively reveal the state and its performance reliability variation process. The experimental data show that with the deepening of the fault degree, the reliability performance decreases gradually.

#### Acknowledgements

This paper is supported by the National Natural Science Foundation of China, (No. U1804141), by the Program for Science and Technology Innovation Talents in Universities of Henan Province (No. 17HASTIT028), by the Key Science and Technology Research Project of the Henan Province (No. 212102210071) and by Educational commission of Henan province (No. 21A460037).

1. Xia Xin Tao, Xu Yong Zhi Modern Statistical Analysis of Variation of Rolling Bearing Performance. Science Press, Beijing, 2016. [Search CrossRef]
2. Prokhorov M. D., Ponomarenko V. I., Karavaev A. S., et al. Reconstruction of time-delayed feedback systems from time series. Physica D, Vol. 203, 2005, p. 209-223. [Publisher]
3. Camplani Massimo, Cannas Barbara Role of embedding dimension and time delay in time series forecasting. IFAC Proceedings Volumes, Vol. 42, Issue 7, 2009, p. 316-320. [Publisher]
4. Qiu Huaxu, Huang Zhangyu, Zhu Hua Application of improved G-P method in time series study. Reconnaissance Science and Technology, Vol. 1, 2014, p. 53-54. [Search CrossRef]
5. Adenan Nur Hamiza, Rani Mohd Salmi Md Noo Predicting time series data at floodplain area using chaos approach. Sains Malaysiana, Vol. 44, Issue 3, 2015, p. 463-471. [Publisher]
6. Small Michael, Tse C. K. Optimal embedding parameters: modeling paradigm. Physica D, Vol. 194, 2004, p. 283-296. [Publisher]
7. Maus A., Sprott J. C. Neural network method for determining embedding dimension of time series. Communications in Nonlinear Science and Numerical Simulation, Vol. 16, Issue 8, 2011, p. 3294-3302. [Publisher]
8. Meng Qing Fang, Zhangqiang, Pan Jin Feng New method for determining minimum embedding dimension based on four-order cumulant. Theory and Practice of System Engineering, Vol. 9, 2005, p. 83-88. [Search CrossRef]
9. Sun Yun, Wang Ying, Meng Xiangfei, et al. Methodology of estimating the embedding dimension in chaos time series based on prediction performance of K-CV_GRNN. International Conference on Information Systems and Computer Aided Education, 2018. [Search CrossRef]
10. Cao L. Practical method for determining the minimum embedding dimension of scalar time series. Physica D: Nonlinear Phenomena, Vol. 110, Issue 1, 1997, p. 43-50. [Publisher]
11. Fraser A. M., Swinney H. L. Independent coordinates for strange attractors from mutual information. Physicalreview A, Vol. 33, Issue 2, 1986, p. 1134. [Search CrossRef]
12. Kim H. S., Eykholt R., Salas J. D. Nonlinear dynamics, delay times, and embedding windows. Physica D, Vol. 127, 1999, p. 48-60. [Publisher]
13. Long Chen, Xin Tao Xia, Hao Tian Zheng, Ming Qiu Chaotic dynamics of cage behavior in high-speed cylindrical roller bearing. Shock and Vibration, Vol. 2016, 2016, p. 9120505. [Search CrossRef]
14. Zhang Shuqing, He Peng, Zuo Yige, Chen Rongfei, Zhang Yun, Liu Wan, Jiang Wanlu Study on characteristics of chaotic singular spectrum and applications in rolling bearing fault diagnosis. China Mechanical Engineering, Vol. 29, Issue 492, 2018, p. 16-22. [Search CrossRef]
15. Li Zhaofei, Renxiaohong, Huang Chencheng Study on nonlinear hyper chaotic characteristics for vibration of rolling bearings. Bearing, Vol. 7, 2016, p. 54-60. [Search CrossRef]
16. Xu Yong Zhi, Xia Xin Tao, Nan Xiang Dynamic analysis of robust test data on rolling bearing vibration based on chaos theory. Journal of Aerospace Power, Vol. 30, Issue 8, 2015, p. 1959-1965. [Search CrossRef]
17. Pan Haiyang, Yang Yu, Zhengjinde, Cheng Junsheng Variable predictive model based RBF class discriminate method. Journal of Aerospace Power, Vol. 32, Issue 2, 2017, p. 500-506. [Search CrossRef]
18. Xiaxintao, Chang Zhen, Li Yunfei Vibration performance prediction and reliability analysis for bearins. Journal of System Simulation, Vol. 30, Issue 4, 2018, p. 1390-1398. [Search CrossRef]
19. Lihongkun, Zhang Zhixin, Li Xiugang, et al. Reliability prediction method based on state space model for rolling element bearing. Journal of Shanghai Jiaotong University (Science), Vol. 20, Issue 3, 2015, p. 317-321. [Publisher]
20. Xiao Wenrong, Cheng Wei, Zi Yangyang Support evidence statistics for operation reliability assessment using running state information and its application to rolling bearing. Mechanical System and Signal Processing, Vol. 60, Issues 61-11, 2015, p. 344-357. [Publisher]
21. Xia Xintao, Jia Chenhui, Wang Zhongyu Fuzzy prediction of rolling bearing friction torque using information poor system theory. Journal of Aerospace Power, Vol. 24, Issue 4, 2009, p. 945-950. [Search CrossRef]
22. Li Dazhong, Chang Cheng, Xu Bingkun Wind turbine gearing temperature prediction based on sample optimization. Journal of System Simulation, Vol. 29, Issue 2, 2017, p. 374-380. [Search CrossRef]
23. Chaozhong Guo, Lin Li, Yuanyuan Hu, Jihong Yan A deep learning based fault diagnosis method with hyperparameter optimization by using parallel computing. IEEE Access, Vol. 8, 2020, p. 131248-131256. [Search CrossRef]

#### Cited By

 IEEE Sensors Journal Xin Li, Fengrong Bi, Xiao Yang, Xiaoyang Bi 2022 M. Sobhana, Gajula Siva Sai Preethi, Guttikonda Harsha Sri, Kancharla Bhanu Sujitha 2022
{"ft_lang_label":"__label__en","ft_lang_prob":0.87650096,"math_prob":0.9987821,"size":36245,"snap":"2022-40-2023-06","text_gpt3_token_len":7795,"char_repetition_ratio":0.18514942,"word_repetition_ratio":0.06417018,"special_character_ratio":0.22091323,"punctuation_ratio":0.14571302,"nsfw_num_words":1,"has_unicode_error":false,"math_prob_llama3":0.9992941,"pos_list":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48],"im_url_duplicate_count":[null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null,2,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2022-10-02T13:46:14Z\",\"WARC-Record-ID\":\"<urn:uuid:83e5e66a-1263-4a6c-a88f-3c0a9e20e2a0>\",\"Content-Length\":\"151825\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:bc64f1cf-02dc-4b5e-a636-d8c9222b6aac>\",\"WARC-Concurrent-To\":\"<urn:uuid:84cb3909-8db9-413a-8304-91e0b067fd75>\",\"WARC-IP-Address\":\"20.50.64.3\",\"WARC-Target-URI\":\"https://www.extrica.com/article/21463\",\"WARC-Payload-Digest\":\"sha1:2O5IAN43GSXMX7SZ2Q7XVXFYIZR34HAD\",\"WARC-Block-Digest\":\"sha1:5ZUKIF7S4C3RLKZI24POJFRP2SULKEJX\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2022/CC-MAIN-2022-40/CC-MAIN-2022-40_segments_1664030337322.29_warc_CC-MAIN-20221002115028-20221002145028-00556.warc.gz\"}"}