Dynamic Prediction for Accuracy Maintaining Reliability of Superprecision Rolling Bearing in Service

A dynamic prediction method for accuracy maintaining reliability (AMR) of superprecision rolling bearings (SPRBs) in service is proposed by effectively fusing chaos theory and grey system theory and applying stochastic processes. In this paper, the time series of a vibration signal is used to characterize the state information for SPRB, and four runtime data points can be predicted in the future, which depends on four chaotic forecasting models to preprocess the time series. Using the grey bootstrap method and sampling from the four runtime data, a large amount of generated data (GD) are gained to analyze the changes in information on bearing service accuracy. +en, using a predefined accuracy threshold to match the Poisson count for the GD, the estimated value of variation intensity is obtained. Subsequently, with the help of the Poisson process, the dynamic evolution process is forecast in real time for AMR of the SPRB for each step in the future. Finally, according to a novel concept for maintaining relative reliability in an SPRB, the failure degree of a bearing maintaining an optimum accuracy status (BMOAS) is effectively described. Experimental investigation shows that multiple chaotic forecasting methods are accurate and feasible with all relative errors below 15%; the reliability of each step in the future can truly be described, and the prediction results for AMR over the same subseries show good consistency; dynamic monitoring of the health status of SPRB can be realized by the degree to which a BMOAS fails.


Introduction
A superprecision rolling bearing (SPRB) refers to a kind of accurate rolling bearing with abilities over a wide range of speed, high rotary precision, low fever, high rigidity, low vibration, and low noise, which is a key part of the equipment-manufacturing industries.
e reliability of SPRB is brought into sharp focus with rapid developments in equipment-manufacturing industries.Maintaining superior accuracy for SPRB in service is the basis of accuracy maintaining in manufacturing equipment.Owing to the influence of various factors, the running accuracy of the bearings is nonstationary degradation, and the reduction process possesses nonlinear dynamic characteristics.Information, such as precision recession trajectory, probability distribution, and accuracy reliability function, varies with the movement process.e prediction problem for accuracy reliability involves the interaction between internal factors and the external environment [1][2][3][4][5].According to the existing research, the reliability theory of rolling bearings is primarily concerned with fatigue failure and static reliability problems and assumes that the lifetime data obey the Weibull distribution or lognormal distribution [6][7][8].However, in practical applications, the lifetime evaluation index of an SPRB rests on accuracy rather than fatigue, which means that the bearing may be far from reaching its fatigue lifetime, while the accuracy has failed.e existing lifetime evaluation systems are not suitable for SPRB, and there is no set of proprietary theories for dynamically predicting the accuracy maintaining reliability (AMR) of SPRB at each point in the future.
Domestic and foreign research studies have made some achievements in the field of reliability prediction for product components, which provide theoretical references and some suggestions for this investigation, although there is no ideal model for AMR.Kim et al. [9,10] considered the complex functional structures and changeable environments of research objects, and their performance degradation was often shown to be nonlinear, nonstationary, and non-Gaussian, so the reliability of multiple models was fused into a prediction by making full use of the multisource information.Liu et al. [11] researched a data fusion method for health indicators in degradation modeling, and a narrower confidence interval for the reliability was obtained.Bian and Gebraeel [12,13] viewed the first passage time of a degradation signal for a related device as the remaining life of the equipment and then updated the lifetime signal parameters according to the Bayesian theory.Hong et al. [14] proposed an adaptive method for predicting the health status of bearings and used empirical mode decomposition and self-organizing mapping methods to analyze the bearing vibration information and evaluate the health status at different stages.Tan et al. [15] proposed a fault evolution-test dependency model to improve the testable performance of mechanical systems for tracking fault growth, and this model can provide guidance to engineers searching for failure information in time.Xia [16] proposed a forecasting method for product reliability along with performance data, without any prior information on probability distributions and trends.In the existing literature, there is only a little information available on the reliability of SPRB, so predicting accuracy reliability is the bottleneck that restricts the development of superprecision components and advanced technologies, which is a pressing problem that needs to be solved.In order to effectively solve the problem of accuracy reliability in product components, this paper presents a novel concept for accuracy maintaining reliability and the relative reliability of SPRB in service.
Accuracy maintaining reliability of SPRB is the probability of an SPRB maintaining optimum service accuracy during a test or in the practical application and the probability of it not exceeding its accuracy threshold as well.AMR can be expressed as a function, where the values of this function are called the degree of accuracy maintaining reliability.
e accuracy maintaining relative reliability (AMRR) is an error in AMR in the future for each step relative to the AMR at its optimal time, which is used to characterize the failure degree of an SPRB in maintaining the optimum service accuracy status.Most attention focused on SPRB concerns its accuracy index, so it is generally necessary to carry out vibration tests on bearings or precision equipment in order to evaluate the service accuracy of an SPRB effectively and then to preprocess and extract the features from the test data for subsequent research.If the vibration data are small, it indicates that the SPRB is running smoothly and maintaining good rotary accuracy; if they are large, it shows that the SPRB is running roughly with a greater axial moving or radial moving on the spindle, which will cause the service accuracy to be greatly reduced.erefore, by analyzing the vibration information for the SPRB, the accuracy parameters in the future can be elucidated and predicted, and a dynamic prediction model driven by the vibration data can be established for AMR.
Chaos theory is one of the important methods in nonlinear analysis, due to the vibration performance of time series with nonlinear dynamic characteristics.Analyzing the performance characteristics based on the chaos theory makes the present bearing research at a new height and opens up a new field.Chaotic forecasting models [17,18] can realize dynamic prediction of a time series [19] accurately, extract the evolution information of the system in the future, and explore the inherent deterministic rules of the system.
e grey bootstrap method [20,21] can transform a small sample of data from poor information into large sample data for traditional statistical theory.
e maximum entropy method [22,23] can be used to obtain the probability density of the sample data, whose main idea is that the solution is the most "unbiased" among all feasible solutions.If a large amount of sample data are given, there will be more accurate sample moments and a more authentic probability density.
e Poisson process [24,25] can record the probability of a failure event effectively, but it is difficult to obtain the variation intensity.If one of the above theories is applied only, it will be one-sided and cause many limitations, in which a dynamic forecast of reliability cannot be achieved as well.e innovative nature of this paper is the fusion of the above theories skillfully, where they complement each other.
e predicted values of four chaotic prediction models provide four data, which consist of small sample data with poor information.e grey bootstrap method can be used for the conversion to a large sample of data for statistical theory, and then the true probability density function can be calculated using the maximum entropy theory; the variation intensity across the large sample of data is found according to a predetermined accuracy threshold; finally, a reliability prediction is realized via the Poisson process, which records the probability of a system before failure.
Based on the process above, the vibration time series of an SPRB is used to characterize the service accuracy information, and multiple dynamic forecasting models from the time series are built to research different aspects of future trends for the SPRB via chaos theory.A large amount of data generated for different aspects are simulated by the grey bootstrap method, and the probability density function of the generated data is constructed using the maximum entropy method, and then the estimated true values and intervals are forecasted according to the given significance level.Meanwhile, the counting principle is used to process the large amount of generated data, and the variation intensity can be acquired.en, the concept of AMR for SPRB is proposed using the Poisson process to establish the AMR models, and the dynamic evolution process of AMR is predicted in real time.Finally, the novel concept of AMRR is used to characterize the failure degree of an SPRB maintaining its optimum service accuracy status in the future.e results of the investigation could generate new ideas in new fields of product performance or accuracy reliability, break through the limitations of traditional reliability theory, and enrich the development of basic reliability theories.e proposed models can not only realize dynamic prediction of AMR for SPRB at each point in the future but also effectively forecast the true values and intervals for the future status information.Moreover, the models can monitor the failure degree for future runtime accuracy relative to the optimum service accuracy as well.

Mathematical Models
2.1.Chaotic Forecasting Models.Suppose the time series X for the accuracy information of rolling bearings is given by where n stands for the nth data for the original series X and N is the number of original data.
According to the theory of phase space reconstruction, the phase trajectory of the time series X can be obtained as follows: where t is the tth phase trajectory; x(t + (m − 1)τ) is the delay value; m is the embedding dimension, which can be obtained by the Cao method; τ is the delay time that can be calculated by the mutual information method; and M is the number of phase points.Phase space reconstruction is the basis for predicting the future evolution of the accuracy information for superprecision rolling bearings.
Suppose that X(M) is the center trajectory (viz., the trajectory of prediction started or the phase space trajectory at the end), there are L reference trajectories similar to the center trajectory, and X(M l ) is the lth reference trajectory.
e chaotic dynamic prediction can be obtained as follows.

Adding-Weight Zero-Rank Local-Region Forecasting
Method.According to the theory of phase space reconstruction, the phase trajectory of the time series X for the accuracy information can be obtained as follows: where M stands for the number of points in the reconstructed phase space with M � N − (m − 1)τ, in which N stands for the number of original data.
Based on the adding-weight zero-rank local-region method, the evolution rule of the phase trajectory is with where X(M + 1) is the prediction result; dl is the Euclidean distance between X(M) and X(M l ); d min is the minimum value of dl; k is the prediction parameter, usually k � 1; and L is the number of reference trajectories.e detailed steps are the following: (1) Using the zero-mean method to preprocess the time series X, the sequence x(t) is obtained with t � 1, 2, . . ., N (2) Reconstructing the phase space (3) Finding L reference trajectories expressed as X(M) reference � X(M 1 ), X(M 2 ), . . ., X(M L )   which are the most proximal points of the center trajectory X(M), and X(M) reference can be obtained by using the Euclidean distance (Equation ( 5)) (4) Acquiring the prediction results X(M + 1)

One-Rank Local-Region Forecasting Method.
e most proximal points X(M) reference around the center trajectory X(M) are fitted by the linear model of X(M + 1) � a + bX(M), with where the points of X(M 1 ), X(M 2 ), . . ., X(M L ) are the proximal points of the center trajectory X(M).e coefficients a and b can be solved by the least-squares method, so X(M + 1) is obtained according to the formula X(M + 1) � a + bX(M), and then the prediction value is efficiently separated.

Adding-Weight One-Rank Local-Region Forecasting
Method.Compared with the one-rank local-region forecasting method, the adding-weight one-rank local-region forecasting method considers the influence weight between each proximal point and the center point; namely, the weight term is added, and the weight is as follows: where k is a forecasting parameter, generally k � 1. e linear fitting of the one-rank local-region can be given by with R � [1, 1, . . ., 1] T .Using the least-adding-squares method to solve the coefficients a and b, Shock and Vibration 3 with viz., e coefficients a and b are acquired by solving equations, and the prediction value is obtained by Equation (8).

Improved Adding-Weight One-Rank Local-Region Forecasting Method.
e forecasting method for the improved adding-weight one-rank local-region is proposed based on the adding-weight one-rank local-region method.
e difference between the two is the definition of correlation between the center trajectory X(M) and the proximal points or reference trajectories X(M l ): the correlation of the proximal points of the adding-weight one-rank local-region prediction method is defined by the Euclidean distance, and the correlation of the improved method is defined by the cosine value of the angle, with where cos(l) is the cosine value of the angle between phase points X(M) and X(M l ); thus, the calculation process of the improved adding-weight one-rank local-region method is similar to that of the adding-weight one-rank local-region method; namely, only the Euclidean distance dl is changed into the cos(l).

Grey Bootstrap Method.
e grey bootstrap method can transform a small sample of data from poor information into large sample data for traditional statistical theory by combining the GM(1,1) method and bootstrap method in the grey system theory.
e GM(1,1) method weakens the randomness of sample data obviously and generates into random numbers regularly, which makes the accumulative generation and data fusion by optimizing the parameters of differential equations.Moreover, it does not alter the probability density function of sample data, so the grey bootstrap method by using GM(1,1) can estimate the predicted values of small sample data accurately and effectively.
Using the above four prediction models: adding-weight zero-rank local-region, one-rank local-region, addingweight one-rank local-region, and improved adding-weight one-rank local-region, four service accuracy information points for the ξth step predicted backward can be forecasted for a superprecision rolling bearing, with the vector Y(ξ) expressed as follows: where y ξ (u) is the uth data in Y for the ξth step and ψ is the maximum value of the step number predicted backward for four prediction models.
In order to meet y ξ (u) ≥ 0 according to the grey model GM(1,1) [26], if y ξ (u) < 0 in Equation ( 13), then a constant c should be selected, making y ξ (u) + c ≥ 0. erefore, in the actual analysis, Y is expressed as follows: Using the bootstrap method [26], B simulation samples of size q, namely, the bootstrap resampling samples V bootstrap , can be obtained by an equiprobable sampling with replacement from Equation ( 14), as follows: where V b is the bth bootstrap resampling sample and B is the times of bootstrap resampling and also the number of bootstrap samples, with where g � 1, 2, . . ., q and b � 1, 2, . . ., B.
According to the grey model GM(1,1), suppose the firstorder accumulated generating operator (1-AGO) of V b is given by e grey generated model based on Equation ( 17) can be described as the differential equation as follows: where u is the time variable and c 1 and c 2 are the coefficients to be estimated.Use the increment to replace the differential, viz., where Δu is equal to the unit interval, 1.Furthermore, assume the generated vector of the mean series as follows: e least-squares solution for Equation (18) with the initial condition 1) is given by 4

Shock and Vibration
where the coefficients c 1 and c 2 are as follows: with According to the inverse AGO, the bth generated data are expressed as follows: erefore, B generated data for the service accuracy can be obtained as the following vector: where w b is the bth generated data.

Solving the Probability Density Function.
Transmuting the generated data Y B for the service accuracy in Equation ( 25) into continuous information, the expression of maximum entropy is defined as follows: where p(w) is the probability density function of the data series Y B .By the maximum entropy principle, the optimal estimation of the density function based on sample information can be obtained, and the main idea of maximum entropy is that the solution is the most "unbiased" among all feasible solutions, given as follows: where S represents the integral interval, namely, the feasible region for the random variable w.
It satisfies the constraint condition: where m j stands for the jth-order origin moment, with β for the highest origin moment order.e entropy can reach its maximum by adjusting p(w), and the probability density function p(w) can be obtained by using the Lagrange multiplier method, as follows: where λ 0 , λ 1 , . . ., λ α are Lagrange multipliers and w is a random variable for service accuracy, with Equation ( 29) is the probability density function of generated data Y B , which is constructed by using the maximum entropy principle.With the help of function p(w), the true value and the confidence interval of the generated series can be implemented as follows.

Parameter Estimation.
According to the probability density function p(w) of the random variable w, the estimated true value X 0 of series Y B is given by If the real number α ∈ (0, 1) exists, then the probability of w α is given by where w α is the α quantile of the density function p(w) and α is the significance level.
For the bilateral quantile, the probability is as follows: where X U and X L are the upper and lower boundaries of the generated series Y B , respectively, and [X L , X U ] is the confidence interval under the α level.erefore, four points for accuracy information for superprecision rolling bearings are fused by effectively combining the grey bootstrap method with the maximum entropy method, and then the true value X 0 and the interval [X L , X U ] for its accuracy information can be predicted for each moment in the future.25)) for superprecision rolling bearings has μ data over the accuracy threshold h, namely, there are μ data falling out of the interval [0, h] of optimum service accuracy, the estimated value for variation intensity θ of the generated sequence Y B is expressed as follows:

Poisson Process
Variation intensity refers to the frequency of accuracy amplitude for superprecision rolling bearings exceeding the optimum accuracy interval, which is the important and characteristic parameter for the runtime accuracy affecting Shock and Vibration the variation process, and it varies with the accuracy threshold.

Dynamic Prediction for Accuracy Maintaining
Reliability.Any counting process can be described using the Poisson process as follows: where i stands for the time variable with i � 1, 2, 3, . . .(i ≥ 1); θ for the variation intensity; e for the number of occurring failure events with e � 0, 1, 2, 3, . .., namely, the serious variation in service accuracy that may cause the bearing failure; and Q for the probability of failure events occurring e times.us, the reliability R for failure events can be obtained using the Poisson process.
When solving for the accuracy maintaining reliability R for superprecision rolling bearings in service, let e � 0; viz., the frequency of occurrence of accuracy failure is zero, which indicates that R is the probability before the accuracy failure occurs.Let i � 1, which indicates that R is the accuracy maintaining reliability in real time; namely, it is the probability of the service accuracy being maintained in the interval [0, h] at optimum accuracy for the currently generated sequence Y B .According to Equation (35), reliability can be expressed as follows: where R(θ) stands for the probability of the running status maintaining optimum accuracy for the superprecision rolling bearings in service.en, the reliability R(θ) of sequence Y B is only a function of the variation intensity θ, and θ can be obtained by Equation (34).In a concrete implementation, if the accuracy maintaining reliability (AMR) R(θ) is not less than 90%, then the reliability of bearing service accuracy is deemed to be good and its possibility of maintaining an optimum service accuracy state is great.If the reliability R(θ) is less than 90% and not less than 80%, the reliability of bearing service accuracy is deemed to be general and its possibility of maintaining an optimum service accuracy state is gradually being reduced.If the reliability R(θ) is less than 80% and not less than 50%, the reliability of bearing service accuracy is deemed to be low and the possibility of maintaining an optimum service accuracy state is poor.If the reliability R(θ) is less than 50%, the bearing service accuracy can be regarded as failed.

Service Accuracy Maintaining the Relative Reliability (SAMRR).
e variation intensity θ 1 for the optimum service accuracy is obtained at the optimal time (generally the initial runtime) for a superprecision rolling bearing (SPRB), and the variation intensity of service accuracy at other times is expressed as θ η with η � 2, 3, 4, . ... According to the concept of relative error in measurement theory, service accuracy maintaining the relative reliability d(η) of an SPRB is obtained to characterize the failure degree of maintaining runtime optimum accuracy.
where R(θ 1 ) stands for the SPRB reliability of maintaining optimal service accuracy at the optimal time, R(θ η ) for others' reliability of maintaining optimum service accuracy at different future times, and d(η) for the failure degree of SPRB in maintaining an optimum service status at different future times.35) and (36).erefore, by analyzing the vibration information of the SPRB, the accuracy parameters in the future can be elucidated and predicted, and a dynamic prediction model driven by the vibration data can be established for AMR. e accuracy maintaining relative reliability (AMRR) is an error in AMR in the future for each step relative to the AMR at its optimal time, which is used to characterize the failure degree of an SPRB in maintaining the optimum service accuracy status.e basic principle of service accuracy classification for an SPRB is as follows:

Evaluation of the
(1) According to the principle of significance in hypothesis testing, if SAMRR d(η) of SPRB is not less than 0%, which shows that the AMR at this assessment time is not less than AMR at the optimum time, it cannot deny that the service accuracy of the SPRB has reached its optimal state; otherwise, it can deny that the service accuracy has achieved its optimal state.(2) When SAMRR d(η) of SPRB is less than 0%, if the absolute value of relative error d(η) is in (0%, 10%], this indicates that the error between the evaluation value and the optimum value is very small.If the absolute value of relative error is in (10%, 20%], this indicates that the error between the evaluation value and the optimum value is gradually increasing.If the absolute value of relative error is greater than 20%, this indicates that the error between the evaluation value and the optimum value is very large. Based on that, the service accuracy of an SPRB is divided into S1, S2, S3, and S4 for a total of four levels: S1: if relative reliability d(η) ≥ 0% of the SPRB, it indicates the bearing service accuracy at a future time reaches the optimum, and the accuracy state has almost no failure possibility.According to the four levels of service accuracy for an SPRB, the evolution process of the failure degree in maintaining an optimum accuracy state can be predicted effectively.SAMRR of an SPRB really refers to the attenuation degree in AMR in the future time relative to the optimum service accuracy state.
e negative value of the SAMRR indicates that the bearing state has attenuation; namely, AMR currently is less than AMR at the optimum time, and the positive value indicates no attenuation.e smaller the relative reliability d(η) is, the worse the bearing service accuracy is and the larger the failure probability of the optimum service accuracy status is.
erefore, the period of d(η) � −20% is the critical time where bearing service accuracy becomes poor.Putting appropriate measures in place before the critical time can avoid serious safety accidents that are caused by a state failure in optimum service accuracy.

Basic Modeling Idea.
A variety of mathematical methods including chaotic prediction methods, the grey bootstrap method, the maximum entropy principle, and the Poisson process are used during theoretical modeling.Each model does not exist or function here as an isolated single model but rather function by complementing and interlocking with each other, which breaks through the limitations whereby one method can only solve a certain type of problem.e idea works as follows: Step 1.Based on a time series X of accuracy information, using the Cao method to obtain the embedding dimension m and using the mutual information method to acquire the time delay τ, the phase space is reconstructed.
Step 2. According to the phase space reconstruction to obtain the phase trajectory of time series X, with the help of four chaotic dynamic prediction models including adding-weight zero-rank local-region, one-rank local-region, adding-weight one-rank local-region, and improved adding-weight one-rank local-region, four forecast values for each step in the future are obtained, composing the small sample Y of size four.
Step 3.With the help of the grey bootstrap method, a large sample Y B is generated from the small sample Y to solve the accurate probability density function.
Step 4. Let the large sample Y B be continuous data, and solve the sample moment of each order, according to the maximum entropy principle to obtain the probability density function, and then calculate the true value X 0 and the interval [X L , X U ] under the α level for each step in the future, and the prediction of accurate information is achieved.
Step 5.According to the given accuracy threshold h, find out the number μ of Y B outside the optimum service accuracy interval [0, h] and then obtain the variation intensity θ for the prediction results of each step; then, accuracy maintaining reliability R(θ) of an SPRB for each step in the future is achieved according to the Poisson formula.
Step 6.By the new concept of service accuracy maintaining the relative reliability, the relative reliability d(η) of an SPRB is obtained to characterize the degree of failure of future runtime accuracy in maintaining the optimum status.

Fault Declaration.
e innovative nature of this paper is fusing the chaos theory and grey system theory and applying stochastic processes skillfully, where they complement each other.Certainly, benefits and limits of the fusion approach proposed exist.It can realize the dynamic prediction of a time series with nonlinear dynamic characteristics accurately and effectively, extract the evolution information of the system in the future, and explore the inherent deterministic rules of the system.Furthermore, for small sample data with poor information, the fusion approach weakens the randomness of sample data obviously and generates into random numbers regularly, which does not alter the probability density function of sample data as well.
Besides, for all the potential fault types, the fusion approach proposed is applicable to the performance data sample, without any prior information on probability distributions and trends.
e only drawback is that the sensitivity to the initial values of sample data is great for bearing performance owing to the use of chaos theory, and the selection of the accuracy threshold has an effect on calculating variation intensity.
Many fault types do exist and have an effect on the running accuracy of the bearings, which can be reflected in the performance data collected.e evaluation models realize dynamic prediction of AMR and monitor the failure degree of the optimum service accuracy status based on the performance data collected in the experiment, so they are applicable to all the potential fault types.

Experimental Results and Discussion
is is a strength test of the lifetime accuracy for superprecision rolling bearings, and the test machine model is an ABLT-1A, which mainly includes the test head seat, test head, transmission system, loading system, lubrication system, and computer control system.e test material is a P4 bearing H7008C provided by SKF, and this kind of bearing has a moving grade, which can almost meet the Shock and Vibration accuracy requirements of P2 grade.e research is conducted at a room temperature of 26 °C, a relative humidity of 53%, a motor speed of 4950 r/min, and a radial load of 13.2 kN.Test time and bearing vibration information are automatically recorded by the computer control system: the sampling frequency of bearing vibration amplitude is 5 kHz, and the root mean square values of vibration amplitudes are obtained by each 10 min; then vibration data are collected by the computer with unit in m•s −2 ; namely, the bearing vibration signal is the root mean square value of vibration amplitudes within 10 minutes.From the beginning of the test, if signi cant variation occurs in the bearing ring or the roller or even surface fatigue spalling, the vibration value of the test machine will obviously increase and service accuracy will be reduced.If the vibration value reaches a certain value, the motor will stop running, and the experiment will be over.
e computer collects a total of 8010 signals; that is, the experiment lasts for 8010 × 10 min, and the vibration signal of time series X obtained for the superprecision rolling bearings is shown in Figure 1.It is not di cult to see that, with the increase in runtime, the vibration value becomes larger, which means that the malignant variation in service accuracy becomes more serious, and then accuracy maintaining reliability will gradually decline.erefore, the potential evolution process of service accuracy can be described as analyzing the vibration signals of a superprecision rolling bearing and then predicting the future running conditions and accuracy maintaining reliability.
Solving phase space parameters: the mutual information method and the Cao method were used to, respectively, calculate the time delay τ and the embedding dimension m of the time series of X 1 , X 2 , X 3 , and X 4 , and the results are shown in Table 1.
e calculation of phase space parameters for the time series is the basis of phase space reconstruction, which is helpful for the following chaotic forecasting.Figures 2-5 are the chaotic prediction results and original data for subseries X 1 , X 2 , X 3 , and X 4 for the 10 steps in the future, respectively, using the adding-weight zero-rank local-region method (marked as method 1), one-rank local-region method (method 2), adding-weight one-rank local-region method (method 3), and improved adding-weight one-rank localregion method (method 4).
From Figure 2, the variation trends for the four forecasting results of subseries X 1 are very similar and relatively stable, with a forecasting value oating around 0.70 m•s −2 , and the di erence with the original data is a little large in the 5 th , 6 th , and 8 th steps, but the maximum di erence is only about 0.10 m•s −2 .From Figure 3, the di erence in subseries X 2 between the four forecasting results and the original data is very small, and the di erence is slightly large in the rst three steps of the one-rank local-region method (method 2), but only about 0.10 m•s −2 .From Figure 4, the original  8 Shock and Vibration data from subseries X 3 jump between 1 m•s −2 and 1.1 m•s −2 , the forecasting results of the adding-weight zero-rank local-region method (method 1) and the one-rank localregion method (method 2) oat around 0.94-1.0m•s −2 , and the forecasting results of the adding-weight one-rank localregion method (method 3) and the improved adding-weight one-rank local-region method (method 4) oat around 1.04-1.12m•s −2 , which shows that the di erence between the four prediction results and the actual values is very small.From Figure 5, the prediction results of method 1 and method 2 for subseries X 4 are around 1.40 m•s −2 , and the results of method 3 and method 4 are around 1.30 m•s −2 , which shows that the di erence between the four prediction results and the actual values is also small.us, when the four chaotic forecasting models are used in time-series forecasting, all the di erences between the predicted values and actual values are small, and these models are accurate and reliable.
In order to illustrate the validity of the prediction results, the relative errors are calculated between the predicted values and the original data from subseries X 1 , X 2 , X 3 , and X 4 .e absolute values of the relative errors are shown in Figures 6-9.
From Figure 6, it is not di cult to see that the trends in the relative errors from the four forecasting results from subseries X 1 are very similar, which shows that the prediction results have good consistency.
e minimum error appearing in the rst step of method 2 is 0.13%, showing that the prediction result is very precise, and the maximum error appears in the sixth step of method 2 but is only 14.68%.From Figure 7, in the prediction models of subseries X 2 , the predicted results from method 4 are the best, with the prediction error around 0-4.48%.e prediction results from method 2 are somewhat poor, and the maximum prediction error appearing is 11.43%; in the other 2 methods, the prediction errors are not more than 8.8%.From Figure 8, the uctuation in relative error from the four 0.9 0.95

Shock and Vibration
prediction results from subseries X 3 is violent because the original data uctuation is violent (showing serration), which indicates that the prediction models can re ect the trend change in the original series.e minimum error appearing in method 4 is 0.11%, and the maximum relative error appears in method 1, but is only 12.70%.From Figure 9, in the prediction models of subseries X 4 , the prediction results of method 3 and method 4 are better than those of method 1 and method 2, and the maximum relative error appears in method 1 but is only 13.63%.
In summary, the prediction results for the four subseries show that relative errors from the chaotic forecasting models including the adding-weight zero-rank local-region method, one-rank local-region method, adding-weight one-rank local-region method, and improved adding-weight onerank local-region method are all very small and all less than 15%, so the four forecasting models can be applied to practical engineering.However, it is di cult to nd out which method is the best or the worst of the four forecasting models for the four subseries because the models' advantages or disadvantages are di erent in di erent subseries.A single prediction method can only re ect one facet of service accuracy in the future, and a prediction value for each step is a character of the true value.Only by fusing and eliciting multiple aspects of the information, a forecast of the true value can be realized.
e grey bootstrap method and maximum entropy principle are fused to be applied to the four forecasting results in each step: using the grey bootstrap method for sample processing the four forecast results, a large number of generated data are simulated; by the maximum entropy principle, the probability density of the generated data is calculated, the predicted true value is gained, and the prediction interval is then obtained under a given con dence level.

True Value Evaluation and Interval Prediction.
In the process of grey bootstrap generation, let the sampling number q 4, the times for bootstrap resampling B 20000, and the signi cant level α 0.
e generated data and probability density of the four prediction values X 1,1 [0.7013, 0.6991, 0.6878, 0.6883] of subseries X 1 in the rst step are shown in Figures 10 and 11. us, the estimated true value X 0 0.6937 and the estimated interval [X L , X U ] [0.6677, 0.7121] are obtained in the rst step of the subseries X 1 .en, the estimated true value and the estimated interval can be acquired in the 2 nd through 10 th steps of subseries X 1 .Similarly, the estimated true value and the estimated interval can be obtained for the 10 steps in the future for the other subseries of X 2 , X 3 , and X 4 .e results are shown in Tables 2 and 3.
From Table 2, each subseries is predicted for 10 steps, namely, the runtime situation of a superprecision rolling bearing is forecast into the future 10 × 10 min.e vibration performance of subseries X 1 , X 2 , X 3 , and X 4 is around 0.69 m•s −2 , 0.93 m•s −2 , 1.01 m•s −2 , and 1.33 m•s −2 into the future 10 × 10 min, respectively, showing that the forecast results for each series have good consistency.From Table 3, the di erence between the upper and lower limits of each estimated interval is very small, showing that the prediction interval is accurate and reliable.At the same time, referring to the original data in Figures 2-5 of each subseries, the estimated interval [X L , X U ] envelopes all the original data; namely, the proposed prediction models can describe the uctuation in service accuracy of a superprecision rolling bearing e ectively and can also realize online monitoring of the dynamic information for the bearing components, so the reliability of the interval prediction reaches (1 − 0/15) × 100 100%.In order to verify that the true value of the fusion can thoroughly re ect the future situation and accurately forecast the trend, the relative errors between the true value and the actual value are calculated, and the results are shown in Figure 12. 10

Shock and Vibration
From Figure 12, the true value after fusing the four prediction values in each step can be seen and its prediction accuracy is fairly reliable, and the average error is 4.68% for subseries X 1 .e maximum error for the fusion true value is only 14.02%, which appears in the 6 th step, and the maximum error before fusing is 14.68%.e average error is 4.01% for subseries X 2 , and its maximum error for the fusing true value is only 6.82%, appearing in the first step, and the maximum error before fusing is 11.43%.e average error is 4.36% for subseries X 3 , and its maximum error for the fusing true value is only 7.85%, appearing in the 5 th step, and the maximum error before fusing is as high as 12.70%.e average error is 4.61% for subseries X 4 , and its maximum error for the fusing true value is only 10.58%, appearing in the 9 th step, and the maximum error before fusing is as high as 13.63%.Obviously, through using grey bootstrap and maximum entropy methods, the forecast errors for the fusing true value for each subseries decrease significantly, which indicates that the true value X 0 after fusion can better reflect the trend change in future service accuracy for superprecision rolling bearings, and the forecasting results are more accurate and feasible and can be better applied to practical engineering.

Dynamic Prediction for Accuracy Maintaining Reliability.
e accuracy true value and interval have been accurately predicted for a superprecision rolling bearing for each step in the future, and on this basis, if we want to realize a prediction for accuracy maintaining reliability in each step, it is necessary to use the Poisson counting process: first, according to Section 3.2, there are 20,000 grey bootstrapgenerated data for the forecast values in each step; then setting the accuracy threshold h � 1.0 m•s −2 to the Poisson count (the value of the threshold h depends on the bearing requirements levels of vibration performance of the main axis system, and the following will analyze it in detail), the number of μ outside the interval [0, 1.0] of optimum service accuracy is obtained for 20,000 generated data in each step, and the variation intensity θ is then acquired according to Equation (34), with the results shown in Table 4. Finally, dynamic prediction for accuracy maintaining reliability at each step can be achieved by using Equation (36), and the results are shown in Table 5.
From Table 4, it is not difficult to see that the number of μ exceeding the threshold h is very similar or even equal in a single subseries.For the X 1 and X 4 subseries, the generated data for the former are all in the interval of optimum service accuracy, but those for the latter are all beyond the interval.For the X 2 and X 3 subseries, the number for the former in the interval of optimum service accuracy is around 2550-3925 and that for the latter is around 10065-12091.e difference in variation intensity θ is also small in the single subseries, and the variation intensity is 0 for subseries X 1 in each step, which indicates that the bearing running conditions are very stable without any malignant variation at that time, and the running state is very good for maintaining optimum service accuracy.e variation intensity is 1 for subseries X 4 in each step, which indicates that the bearing running condition is extremely poor, and the bearing may have incurred an accuracy failure.e variation intensity is between 0.1275 and 0.1963 for subseries X 2 in each step, showing that the bearing is beginning to mutate, bearing running condition is becoming poor, and the probability of maintaining the optimum service status is gradually decreasing.In addition, this time point can provide a reference point for online health monitoring.e variation intensity is between 0.5033 and 0.6046 for subseries X 3 in each step, showing that the bearing variation is serious, and the operating condition is very poor, so maintenance measures should be performed in time to avoid serious accidents.
erefore, the prediction results demonstrate good consistency for the single subseries.Meanwhile, subseries X 1 , X 2 , X 3 , and X 4 stand for different running stages of a superprecision rolling bearing in service, showing that with the increase in operation time, the variation intensity will increase gradually, and the degree of malignant variation will be more serious, which results in a decline in the accuracy maintaining reliability.
From Table 5, the accuracy maintaining reliability can reach 100% for subseries X 1 in the next 10 steps, and this is because the time series X 1 is the initial stage of the experiment, the vibration amplitude and spindle runout are small with a high service accuracy, and the operation is very safe and reliable.e accuracy maintaining reliability is between 80 and 90% for subseries X 2 in the next 10 steps, showing that accuracy maintaining reliability is general for the superprecision rolling bearing, and the possibility of maintaining an optimum service state is gradually decreasing.us, the bearing is gradually beginning to mutate, and a potential deterioration in accuracy is occurring in the internal part of the bearing.e accuracy maintaining reliability is between 50 and 80% for subseries X 3 in the next 10 steps, showing that the accuracy reliability is poor for the superprecision rolling bearing, and the malignant variation is serious, so remedial measures should be taken in time.e reliability is less than 50% for subseries X 4 in the next 10 steps, indicating that the superprecision rolling bearing is remarkably unreliable, and the bearing may have incurred an accuracy failure or even significant wear.
By analyzing the variation intensity and accuracy maintaining reliability at different time periods for X 1 , X 2 , X 3 , and X 4 , the dynamic evolution can be recognized Shock and Vibration as illustrating the service accuracy of a superprecision rolling bearing.Compared with the total service length of 80010 × 10 min, the time of prediction for the future 10 steps can be viewed as instantaneous.e instantaneous reliability is the mean value of the prediction values over the next 10 steps, which stands for the accuracy maintaining reliability at the end of each subseries.e results are shown in Figure 13.
From Figure 13, the reliability at the end of each subseries can describe the trend variation in service accuracy well for a superprecision rolling bearing: at the end of each subseries X 1 , the accuracy maintaining reliability is large, showing that the possibility of maintaining an optimum service state is great, and the super-precision rolling bearing has a high running accuracy; at the end of each subseries X 2 , the reliability begins to decline, showing that the possibility of maintaining an optimum service state decreases gradually, and the running accuracy of the bearing is slowing down; at the end of each subseries X 3 , the reliability decreases rapidly, showing that the running accuracy of the bearing drops fast; and at the end of each subseries X 4 , the reliability reaches the minimum, showing that the service accuracy of the bearing has been lost already, and even fatigue failure has set in.
From Tables 4 and 5 and Figure 13, the accuracy threshold h determines the variation intensity in service accuracy for each step in the future, thereby a ecting the reliability of each step.However, the threshold value h depends on the demanding degree of vibration performance for superprecision rolling bearings in machine tools, motors, or other systems.e accuracy maintaining reliability for   12 Shock and Vibration subseries X 1 , X 2 , X 3 , and X 4 is discussed by using different accuracy thresholds in Tables 6-9.
According to the prediction series for accuracy maintaining reliability under different accuracy thresholds h from Tables 6-9, it can be seen that the smaller the accuracy threshold is, the lower the reliability for the next step is during the dynamic prediction process for accuracy maintaining reliability for each subseries; conversely, the higher the reliability is.e value size of the threshold reflects the bearing accuracy sensitivity of the spindle system.erefore, in practice, the corresponding accuracy threshold h should be designed in advance according to the spindle system requirements levels for vibration performance of a superprecision rolling bearing.Real-time monitoring of accuracy information and then obtaining the corresponding reliability can discover failures and hidden dangers in time and thus avoid bad accidents.

Service Accuracy Maintaining the Relative Reliability.
From Tables 6-9, it is not difficult to see that the prediction reliability for time series X 1 is the highest under the same threshold value for the four subseries.us, the best period for maintaining an optimum service accuracy is the time interval X 1 for superprecision rolling bearings, where accuracy maintaining reliability is R(θ 1 ), and the accuracy maintaining reliability of the other subseries X 2 , 3 , and X 4 is R(θ 2 ), R(θ 3 ), and R(θ 4 ) according to Section 2.4.3, respectively.Taking h � 1.0 m•s −2 and h � 1.2 m•s −2 as an example to analyze the relative reliability of maintaining service accuracy, when h � 1.0 m•s −2 , the mean values for prediction reliability are solved for each subseries in the next 10 steps, obtaining R(θ 1 ) � 100%, R(θ 2 ) � 85.07%, R(θ 3 ) � 57.23%, and R(θ 4 ) � 36.79%; when h � 1.2 m•s −2 , obtaining R(θ 1 ) � 100%, R(θ 2 ) � 100%, R(θ 3 ) � 99.06%, and R(θ 4 ) � 38.78%.With Equation (37), the degree of failure d(η) for accuracy maintaining an optimum service state is acquired for superprecision rolling bearings at different service time intervals, and the results are shown in Figure 14.
From Figure 14, when the accuracy threshold h � 1.0 m•s −2 and η � 2, the relative reliability is d(η) � 0% at the end of subseries X 2 , which shows that the accuracy reliability at the end of subseries X 2 is not less than its reliability at its best period of service, and the bearing service accuracy is undeniably at its optimal state at this time.When η � 3, the relative reliability is d(η) � −0.94% ∈ [−10%, 0%) at the end of subseries X 3 and d(η) is close to 0%, which shows that the service accuracy at this time is normal, and the failure probability of maintaining an optimal accuracy state is low; viz., its service accuracy is close to the optimal state.When η � 4, the relative reliability is d(η) � −61.22% < −20% at the end of subseries X 4 , which shows that the service accuracy at this time is very poor, the failure probability of maintaining an optimal accuracy state is quite high, and the service accuracy of the bearing may have failed.When the accuracy threshold h � 1.2 m•s −2 and η � 2, the relative reliability is d (η) � −14.93% ∈ [−20%, −10%) at the end of subseries X 2 , showing that the service accuracy at this time becomes gradually worse, and the failure probability of maintaining an optimal accuracy state increases gradually; when η 3 and η 4, the relative reliability is d(η) −42.77% < −20% and d(η) −63.21% < −20% at the end of subseries X 3 and X 4 , respectively, showing that the service accuracy is very poor, the failure probability of maintaining an optimal accuracy state is quite high, and the service accuracy of the bearing may have failed.e point in time for relative reliability d(η) −20% is the critical point at which the service accuracy begins to worsen.us, as a matter of practical engineering, checking and repairing the bearing in time and then taking appropriate remedial measures are necessary steps before reaching a critical moment for a bearing in service. is can avoid serious safety accidents caused by a failure in accuracy for a superprecision rolling bearing.
In summary, during the process of chaotic forecasting for time series, the adding-weight zero-rank local-region method, one-rank local-region method, adding-weight one-rank localregion method, and improved adding-weight one-rank localregion method are all feasible and accurate, with the maximum prediction error being no more than 15%, thus meeting general forecasting requirements in engineering practice.Combining the grey bootstrap method with the maximum entropy principle to fuse the four prediction values of each step e ectively, the future accuracy true value and estimated interval are obtained for a superprecision rolling bearing.According to the Poisson counting process and grey bootstrap-generated data, a dynamic prediction for accuracy maintaining reliability is achieved for superprecision rolling bearings for each step in the future, and the prediction results can e ectively monitor the variation information on bearing service accuracy.Finally, a new concept is proposed for relative reliability of maintaining accuracy, which can e ectively predict the degree of failure for maintaining an optimal accuracy state for a bearing.

Conclusions
(1) By fusing chaotic prediction models into the grey bootstrap principle and maximum entropy principle, the accuracy true value and uctuation interval can be precisely predicted for superprecision rolling bearings for each step into the future.e prediction values and actual values demonstrate good consistency, and the prediction error is very small with a maximum relative error of only 14.02%.(2) Combining the grey bootstrap principle with the Poisson process, a dynamic prediction method is proposed for accuracy maintaining reliability in superprecision rolling bearings, which can realize the prediction of accuracy maintaining reliability for each step into the future and can reveal the mechanism of in uence between the running time and service accuracy.(3) According to the varying curves of accuracy maintaining reliability at di erent time intervals and the in uence of threshold accuracy on reliability analysis, the evolution information for accuracy maintaining reliability can be predicted in real time for superprecision rolling bearings, and failures or hidden troubles can therefore be discovered in time.(4) Based on the novel concept of accuracy maintaining reliability for superprecision rolling bearings, the possibility and degree of failure in maintaining an optimal accuracy state are predicted e ectively.Some intervention measures for maintenance or replacement can be implemented in a timely manner before the possibility of failure in service accuracy becomes critical.(5) e proposed models can not only realize dynamic prediction of accuracy maintaining reliability for superprecision rolling bearings at each point in time but also realize an estimation of true value and interval for the accuracy attribute.

2. 4 . 1 .
Counting Process.If the accuracy generated sequence Y B (in Equation (

Figure 1 :
Figure 1: Vibration time series X of a superprecision rolling bearing.

Figure 9 :Figure 10 :Figure 11 :
Figure 9: Absolute value of the relative errors from the prediction results from subseries X 4 .

Figure 12 :
Figure 12: Prediction error of the estimated true value X 0 for each subseries.

Figure 14 :
Figure 14: Accuracy maintaining the relative reliability for superprecision rolling bearings.
Failure Degree of the Optimum Service Accuracy State.Accuracy maintaining reliability (AMR) of an SPRB is the probability of it maintaining the optimum service accuracy status and the values of service accuracy not exceeding the accuracy threshold.e concept of AMR is proposed using the Poisson process to establish the AMR models and predict the dynamic evolution process of AMR in real time.e variation intensity of the generated sequence Y B is obtained by recording the number of sample data falling out of the interval [0, h] of optimum service accuracy; then, based on the Poisson process, the estimated value for the AMR of an SPRB can be obtained according to the Equations (

Table 1 :
Phase space parameters of 4 subseries.

Table 2 :
Estimated true value X 0 for each subseries for 10 steps in the future.

Table 3 :
Estimated interval [X L , X U ] for each subseries for 10 steps in the future.

Table 4 :
Number of μ exceeding the threshold h and variation intensity θ.

Table 5 :
Accuracy maintaining reliability for the subseries at each step.

Table 6 :
Accuracy maintaining reliability at different accuracy thresholds for subseries X 1 .

Table 7 :
Accuracy maintaining reliability at different accuracy thresholds for subseries X 2 .

Table 8 :
Accuracy maintaining reliability at different accuracy thresholds for subseries X 3 .

Table 9 :
Accuracy maintaining reliability at different accuracy thresholds for subseries X 4 .