Abstract

A new residual life prediction method for complex systems based on Wiener process and evidential reasoning is proposed to predict the residual life of complex systems effectively. Moreover, the better maintenance strategies and decision supports are provided. For the residual life prediction of complex systems, the maximum likelihood method is adopted to estimate the drift coefficient, and the Bayesian method is adopted to update the parameters of Wiener process. The process of parameters estimation and the probability density function (PDF) of the residual life are deduced. To improve the accuracy of the residual life prediction results, the evidential reasoning (ER) is used to integrate the prediction results of Wiener process. Finally, a case study of gyroscope is examined to illustrate the feasibility and effectiveness of the proposed method, compared with fuzzy theory, which provides an important reference for the optimization of the reliability of complex systems and improvement.

1. Introduction

Nowadays, the residual life prediction is becoming an attention and hot topic in engineering practice. Especially in the performance study of pipeline, oil depot, engine, LED lamp control system, storage tanks, gyroscope, and other complex systems, the residual life prediction has been widely applied. If the residual life of complex systems can be estimated at the beginning of the performance degradation of the systems, and if the performance can be evaluated effectively [1], especially when no major damage has been caused, the happening of unnecessary faults will be avoided to some extent and the maintenance costs will be reduced greatly [2]. Therefore, it is essential to study the residual life prediction of complex systems.

In the traditional life experiment, failure data can be used to predict the life of equipment and device. On the one hand, many failure data cannot be obtained in a short time [3]. On the other hand, some of the systems’ characteristics will degrade over time. However, a large amount of information related to reliability and residual life is included in the degradation data. Therefore, the degradation data can be used to predict the residual life of complex systems. Through the acceleration test, both reliability and life information under different environmental conditions can be obtained. At present, the commonly used acceleration equations include Arrhenius model, Eyring model, Power Law Inverse (IPL), and Hallberg-Peck model. Among them, Arrhenius and Eyring models are always used for thermal accelerated stress in the systems, while temperature is the major factor causing aging. IPL is always used for other different kinds of accelerated stresses except thermal accelerated stress. Hallberg-Peck model synthesizes temperature and humidity that can describe the aging test more accurately under conditions of temperature and humidity.

In current study, the prediction methods can be classified as the physics-based prediction method [4] and data-driven prediction method [5]. With the development of related technologies such as signal acquisition and signal processing, abundant data of system operation can be obtained. According to these data, the corresponding mathematical model is utilized. The method is mainly composed of two kinds of techniques, which are named as artificial intelligence and probability statistics. The artificial intelligence methods have a higher degree of the data fitting, but they cannot obtain the uncertainty of the prediction. The probability statistics can predict the future state better.

Firstly, Wiener process is a kind of the most commonly used probability statistical model, which can describe the degradation process accurately [6, 7]. The degradation model based on Wiener process has obvious advantages in mathematics [8]. Wiener process has been widely applied in many fields [911]. Wang proposed a reliability assessment method that fused both prior and on-site degradation data [12]. To update the parameters of Wiener process, the prior data was used to determine the parameters, and the on-site data was integrated by Bayesian method. In [13], the degradation of aeroengine was modeled based on multistage Wiener process. Based on the historical degradation data and the historical failure time data, the prior distribution of the model parameters was estimated by using the expectation maximization (EM) algorithm. In [14], the recursive filtering algorithm and EM algorithm were combined to update the parameters of the Wiener process, and the method was applied to the inertial navigation system to predict the residual life accurately.

Secondly, the evidential reasoning rule (ER rule) proposed by Yang and Xu [15] in 2013 is also a kind of probabilistic statistical model. As a multiattribute decision-making method, the ER rule is widely used to solve the fuzziness and uncertainty of parameters in the evaluation and decision-making process. In [16, 17], the authors have described how to handle the multiattribute decision analysis under uncertainty by the ER rule that is developed based on the decision theory and D-S theory of evidence [15]. Meanwhile, the ER rule is suitable for solving and expressing uncertainty [1820].

To solve above problems, the prediction framework of the residual life of complex systems is constructed, and the corresponding model of complex systems based on Wiener process and evidential reasoning is established in this paper. The other contents of this paper are mainly organized as follows. In Section 2, the structure of the residual life prediction model of complex systems is constructed, and the problems to be solved are described in detail. The residual life prediction model of complex systems based on Wiener process and evidential reasoning is established in Section 3, and the details of modeling steps are organized and summarized systematically. In the Section 4, a case study about the residual life prediction of gyroscope is done, based on fuzzy theory and evidential reasoning, respectively, and the simulation results are analyzed. Finally, the practicability of the prediction methods proposed and future work are summarized in the Section 5, as well as the reference of other similar products.

2. New Residual Life Prediction Model and Problem Formulation

2.1. Structure of the New Residual Life Prediction Model

The residual life of complex systems is influenced by many factors. They can be divided into two categories, which are named as the internal factors and the external factors, respectively. The internal factors are composed of materials, structure, and so on, while the external factors include debugging, operation, and maintenance. Thus, it is hard and complex to establish an accurate and complete mathematical prediction model of residual life of complex systems [21]. However, the residual life of the systems can be reflected by some so-called health indicators [22]. Suppose denotes the th health indicator for the monitoring, . By integrating all the health indicators, a new prediction model of residual life is established, whose structure is shown as Figure 1.

Figure 1 shows the structure of the new model that mainly includes two parts. In the first part, the residual life prediction model is established based on the selected health indicators and Wiener process. In the second part, the prediction results are integrated by the ER rule to obtain more accurate prediction results.

2.2. Problem Description for Residual Life Prediction

Based on Figure 1, the following three problems should be studied in order to predict the residual life of complex systems.

Problem 1. In Figure 1, the prediction of residual life is necessary. Thus, it is essential to establish a prediction model to get the indicators’ prediction value with the obtained indicators’ information. Therefore, Problem 1 concentrates on how to establish the following model:where . denotes the prediction value of the residual life of indicator , which is a time series in nature. represents a nonlinear function.

Problem 2. In practice, in order to represent the relative importance of the health indicators, it is of much necessity to allocate an adaptive weight to each indicator. Thus, Problem 2 concentrates on how to establish the following model to adjust the weight and update the weight coefficient according to the whole available observation data of all indicators:where . represents a nonlinear function to determine the weighting coefficient.

Problem 3. To obtain the comprehensive prediction results, it is necessary to integrate various indicators and reprocess the prediction value of the residual life. Therefore, Problem 3 concentrates on how to establish the following integration model to obtain the residual life level of complex systems:where represents a nonlinear function.

Based on the above three problems, the concrete models are shown as follows.

3. Residual Life Prediction Model

In this paper, the following two models are used to solve the above three problems. Firstly, based on Wiener process, a degradation model of complex systems is established to analyze and describe the performance degradation process. In order to eliminate the tendency items, the incremental degradation is adopted. Secondly, the parameters of the degradation model are estimated and updated, with the maximum likelihood method to estimate the parameters, and the Bayesian method to update the parameters. Thirdly, the residual life prediction model of complex systems is established based on evidential reasoning, compared with fuzzy theory, and the prediction results of Wiener process are integrated to obtain more accurate prediction results.

3.1. Prediction Model of Residual Life of Complex Systems Based on Wiener Process
3.1.1. Degradation Model of Complex Systems

Wiener process has been widely used in the modeling of degradation process because of the excellent capability of analyzing and describing the process of degradation [23, 24]. So, it is used to model the performance of complex systems in this paper.

If is the degradation measure of the performance of the systems at time , the degradation model of Wiener process based on the drift parameter iswhere is drift coefficient, is diffusion coefficient, is a short time interval, and is the Standard Brownian Motion which has two features as follows.

Feature 1. .

Feature 2. For any two different times and , , , , and are independent.

According to the characteristic of , it is assumed that obeys the normal distribution, that is, . So, its degradation amount mean is , and degradation amount variance is . The mean reflects the trend of the degradation data affected by the accelerated stress. So, it is considered that the drift coefficient is related to the acceleration stress. The variance reflects the degree of deviation between the degenerate data and the mean, which is usually unrelated to the acceleration stress. So it is considered that the diffusion coefficient has nothing to do with the acceleration stress. In this paper, the random parameters and are used to describe the individual difference of complex systems. Since there is no sufficient prior information of the parameters to be estimated, the prior distribution without information is adopted.

The failure threshold of the performance of some complex systems is often a fixed value, which is set as a constant . When the current degradation data reaches the threshold value for the first time, the system is failed. So, its residual useful life can be defined as

It is proved that its first failure threshold time is subject to inverse Gaussian distribution, and its probability density function is

And the cumulative distribution function iswhere is probability density function of the standard normal distribution.

In the accelerated degradation test, temperature and humidity are the most common accelerated stresses that can aggravate the reaction to degenerate the product. Both temperature and humidity are important factors for reliability of some complex systems. High temperature and humid environment may make the service life decrease, so the entire life of the system is seriously affected. In this way, the Hallberg-Peck acceleration model is used to construct the relationship between the drift coefficient with the temperature and humidity stresses. Considering the influence of temperature and humidity comprehensively, the Hallberg-Peck model can describe the accelerated degradation test of the product under the condition of temperature and humidity accurately. Its expression iswhere is a constant, , is relative humidity, and is activation energy in the unit of electron volts (). is Boltzmann constant whose value is 8.617 × 10−5 eV/°C. So, the unit of is °C. In addition, is relative temperature. The degradation model can be obtained through substituting (8) into (4) as

3.1.2. The Estimation of Parameters

The accelerated stresses in the accelerated degradation test are set to . There are six samples for each acceleration stress, and the measured time gradient of each sample is . The degradation amount of the th sample under the stress is at time , and the performance degradation is 0 at the initial time; is the increment of degradation. By the nature of , , we have

In the period of , the probability density function (PDF) of the performance degradation data is

The likelihood function of and can be obtained in

The maximum likelihood estimation can be used to estimate the values of and in (13). But the solving process is very complicated. Therefore, the prior distribution of the parameters is determined firstly, and the parameters are updated by the Bayesian method secondly; then the distribution of the parameters can be obtained. Based on the Bayesian theory, the posterior distribution can be represented as where is the probability density function of posterior distribution, is likelihood function, and is the probability density function of prior distribution.

If the distribution of the measured data of the first set sample in stress is a normal distribution, where the parameters of and are noninformative prior distribution, the prior distribution probability density function can be expressed as

Let , , then , the joint posterior distribution of is:

First of all, to update the parameters of the estimated parameters , can be regarded as redundant parameters. To integrate the posterior distribution of the integral, the excess parameters are removed, shown in

The equation is finished as

is transformed as and then put it into (18): where ,  . From (19), it can be seen that the marginal posterior distribution of follows the mean value of , and the scale parameter fits the normal distribution of .

Update the estimated parameter in

From (20), it can be seen that the posterior distribution of is proportional to the probability density function of the inverse Gamma distribution. The shape parameter is and the scale parameter is . Based on (19) and (20), the mean value of the posterior distribution of is , and the scale parameter is . The shape parameter of the posterior distribution of is , and the scale parameter is . The estimated values of the parameters and under the stress of are

In the same way, the estimated parameters can be obtained under .

To , evaluate the logarithmic of (10):

So, and are combined:

Substitute the estimated value of (21) into (23):

The estimated value of can be calculated as

Based on (8), the estimated value of is

Since the diffusion coefficient is independent of the acceleration stress, its estimated value is

After the values of the parameters in the degradation model have been obtained, the mean and variance of the degradation data under normal stress level can be obtained according to the accelerated degradation model. And the corresponding failure distribution can be obtained by substituting them into (6).

3.2. Residual Life Prediction Model Based on Evidential Reasoning

As mentioned in Section 1, considering the advantages of evidential reasoning in theoretical research and practical application, the prediction model of residual life of complex systems based on evidential reasoning is established.

The core of the ER rule is a kind of distributed type with a belief structure. In this way, any input information, whether it is accurate data or inaccurate date with uncertainty, subjective or objective judgment, can be consistently modeled based on the power set. Moreover, the ER rule also handles with the probabilistic uncertainty by integrating the reliability of all indicators in (24). The details are as follows.

(1) Transformation of All Kinds of Input Indicators. On the one hand, different indicators can be described in different ways. On the other hand, there is always ignorance and vagueness in the prediction of residual life of complex systems. In addition, both qualitative and quantitative information are considered with uncertainties [25]. So, it is of much importance to transform various types of health indicators into belief distribution.

Define distinctive residual life grades as represented bywhere is the th residual life grade. It is assumed that is inferior to . In addition, could be determined by circumstances, including special conditions and personal experience. For example, three grades “Long” (), “Average” (), and “Short” () can be chosen in the prediction of residual life of complex systems.

Thus, the distribution can be described aswhere , , and denotes a degree of belief. The above distributed assessment reads that the indicator is assessed to the grade with the degree of belief , which can be calculated aswhere , denotes a quantitative value that can be judged as a referential value .

(2) Integration of Multi-Indicators. For complex systems, the residual life at time relies on all health indicators, which means the residual life can be obtained through multiple health indicators. The information in (28) can be integrated by the ER rule. The results of the residual life prediction can be expressed as follows:where the distribution represents the assessment of the residual life of complex systems at time , from which the detailed life grades and belief degrees can be seen clearly. can be calculated as follows [15]:where denotes the weight coefficient of the th indicator at time and denotes the reliability of the th indicator.

(3) Utility Evaluation. In order to visually represent the distributed assessment results, the expected utility is introduced. is the utility of evaluation grade , and the expected utility of the prediction result is defined as follows:

(4) Ways to Calculate the Adaptive Weighting Coefficient. The weighting coefficient represents the importance of the indicators for prediction results [15, 26]. In this paper, the adaptive weight is calculated with CVBW method [27] shown as follows.

Suppose the observed data of indicator within is denoted by , , whose mean value is . So, can be calculated by

denotes the mean square error of the observed data of indicator within time , which can be calculated by

Also, there is

The weighting coefficient of within can be calculated by

3.3. A New Algorithm for Residual Life Prediction of Complex Systems Based on Wiener Process and Evidential Reasoning

Based on the above analysis, a new algorithm for the residual life prediction of complex systems based on Wiener process and evidential reasoning is proposed in this section. The details are shown as follows.

Step 1. Construct the framework of the residual life prediction model of complex systems, and comb through the problems to be solved.

Step 2. Based on the working mechanism of complex systems, establish the degradation model of the systems, as is shown in (9).

Step 3. Estimate and update the parameters of the degradation model, using the maximum likelihood method and Bayesian method, respectively, as is shown in (21) and (27).

Step 4. Select the life indicators of complex systems, and obtain the residual life of complex systems based on Wiener process.

Step 5. Establish the residual life prediction model of complex systems based on evidential reasoning, as is shown in the Section 3.2.

Step 6. Calculate the adaptive weight, as is shown in (37). Based on the ER rule, make a precise prediction for the residual life prediction results of Wiener process.

4. Case Study

To prove the practicability of the proposed prediction models of complex systems, the residual life of gyroscope is studied in this paper. Considering the fact that the residual life of gyroscope is mainly determined by its drift coefficient which changes as time, the drift coefficient for zero and once term are chosen to measure the residual life of gyroscope.

4.1. Residual Life Prediction of Gyroscope Based on Wiener Process

Under the condition of constant temperature 20°C and humidity 70%, the drift coefficient of the gyroscope is tested once a month, and the degradation values are obtained. To get more accurate results, the test is repeated. Figure 2 shows the results.

In Figure 2, denotes the drift coefficient for zero item and denotes the drift coefficient for once item. Obviously, the two groups of degradation data are roughly linear, random processes, whose increments are not strict. Thus, the degradation process can be determined in Wiener process. In order to fit the data, the increment is introduced based on (4). According to the characteristics of Wiener process, , it can be described by the probability distribution . The fitting results are shown as Figures 3 and 4.

As can be seen from Figures 3 and 4, the drift coefficient degradation process of gyroscope is a Wiener process. Using the method described in the Section 3.1.2, four groups of degenerate data are updated to obtain the normal inverse Gamma posterior distribution parameters of and , which are shown in Table 1.

Remark 4. The parameters are set as shown in Table 2.

Finally, the prediction results of the residual life of gyroscope based on and are obtained, which are shown in Tables 3 and 4, respectively.

According to Tables 3 and 4, the prediction results of the residual life of gyroscope based on and are close to the real value, which has proved the accuracy of Wiener process in this paper. Therefore, the method proposed in the Section 3.1 can provide a reference for predicting the residual life of gyroscope in engineering practice.

Remark 5. Based on the actual degradation of gyroscope, the residual life of gyroscope is generally thought to be linear decreasing tendency.

4.2. Residual Life Prediction of Gyroscope Based on Evidential Reasoning

Based on evidential reasoning, the prediction results of the residual life of Wiener process are integrated. The details are as follows.

Firstly, calculate the weighting coefficient. The results are shown in Table 5.

Secondly, calculate the evidence reliability. Based on expert knowledge, the reliability of and are 0.8262 and 1, respectively. In addition, the assessment grades are set as “Very Long” (), “Long” (), “Average” (), “Short” (), and “Very Short” (), shown as follows.

In general, if the total life of gyroscope is 200 hours, then denotes the residual life for 200 hours, for 150, for 100, for 50, and for 0.

The standard values of each grade are 200, 150, 100, 50, and 0, respectively, denoted by , , , , and . Hence, the evidence is shown in Table 6.

Thirdly, integrate the evidence. Based on the method in the Section 3.2, the integration results are shown in Table 7.

It can be seen from Table 7 that when the initial value of the degradation of the gyroscope is 0.10, the probability that the residual life of gyroscope is in the state is 0.0283, and the probability of being in state is 0.9717. So, the belief distribution can be described as follows:

Therefore, we have a probability of 0.0283 that the residual life of gyroscope is very long, 0.9717 that the residual life is long, and 0 for other grades.

Similarly, the other results are shown as follows.

Therefore, the ER rule has well explained the prediction results of Wiener process, which is consistent with the actual observation results. In the engineering practice, when the drift coefficient of gyroscope is reduced to 0.25, it can be considered that the residual life of gyroscope is very short, and a new gyroscope should be replaced for use.

Finally, all the above methods to predict the residual life of gyroscope are shown in Table 8.

From Table 8, the common residual life prediction results of gyroscope can be obtained by Wiener process, while a specific prediction result can be obtained by the ER rule, which can integrate all of the former results.

To clearly demonstrate the superiority of the proposed methods, traditional fuzzy theory has been used to deal with the residual life prediction results of Wiener process [2830]. Assuming that the set , and that the fuzzy subsets are and , respectively, based on Table 6, the membership of each subset is shown in Table 9.

Compared with Table 8, Table 9 shows that the final results obtained by fuzzy theory cannot be explained clearly, which are only some uncertain values. Therefore, the prediction results of the ER rule are proved to be more feasible and effective.

5. Conclusion

In this paper, a new model based on Wiener process and evidential reasoning is proposed to predict the residual life of complex systems. The relative importance of each health indicator is considered, which provides an effective method for predicting the residual life of complex systems. The model proposed in this paper is mainly composed of two parts that are prediction and integration. Through the case study of gyroscope and a comparison with fuzzy theory, it is proved more feasible and effective in engineering practice. As offline prediction proposed in this paper is unable to precisely and timely reflect the systems’ life level, its effectiveness, and validity to deal with online problems and more complex problems should be further tested. At the same time, the reliability of all health indicators is crucial for the residual life prediction of the whole system and it is influenced by various factors, including the noise, sensor’s characteristics, and malfunction. In addition, the correlation between each indicator needs to be considered. Hence, further study on how to obtain the reliability and correlation of indicators by integrating all factors based on online prediction needs to be done in the future.

Conflicts of Interest

The authors declare that they have no conflicts of interest.