Research Article | Open Access

Volume 2019 |Article ID 9752920 | https://doi.org/10.1155/2019/9752920

Xiaoping Liu, Zhenyu Wu, Dejun Cui, Bin Guo, Lijie Zhang, "A Modeling Method of Stochastic Parameters’ Inverse Gauss Process Considering Measurement Error under Accelerated Degradation Test", Mathematical Problems in Engineering, vol. 2019, Article ID 9752920, 11 pages, 2019. https://doi.org/10.1155/2019/9752920

A Modeling Method of Stochastic Parameters’ Inverse Gauss Process Considering Measurement Error under Accelerated Degradation Test

Accepted30 Apr 2019
Published21 May 2019

Abstract

To solve the problem that the individual differences and the measurement errors affect the accuracy of life estimation in accelerated degradation test, the inverse Gauss process with stochastic parameters is applied in the accelerated degradation test with the consideration of the influence of individual differences, and the analysis of measurement uncertainty is carried out. An inverse Gauss accelerated degradation model considering both individual differences and measurement errors is established. In the maximum likelihood estimation of parameters, Genetic Algorithm and Monte Carlo integral are used to solve the problems caused by complex integral and the unobservable measurement errors in the calculation process. Finally, the proposed method is verified by the Monte Carlo simulation under the constant accelerated stress and step accelerated stress and the illustrative example of electrical connectors under the constant acceleration stress, respectively. The results show that the modeling tool is useful for improving the accuracy of the life prediction and the reliability evaluation.

1. Introduction

Prognostics Health Management (PHM) is a new management method of complex engineering system, which integrates fault diagnosis, life prediction, and health management. And it has been widely applied in different fields, such as communications, electronics, military, and aviation. As a key component of PHM, an appropriate life prediction method is significant for the product reliability evaluation, and it contributes to formulate reasonable health management schemes. It also improves the safety of industrial production; especially it saves the maintenance cost [13].

As the most critical part of the degradation analysis, the degradation models are currently divided into the regression model based on degradation trajectory analysis and the stochastic process model to characterize product performance degradation. Influenced by the changes in the external environment, the manufacturing process and the operation conditions, the degradation process of the products is with some randomness. Therefore, it is reasonable to use the stochastic process to model product degradation, and this is also a widely used modeling method at this stage. Three kinds of commonly used stochastic process models are respectively the Gamma process model, the Wiener process model, and the Inverse Gaussian (IG) process model. However, the degradation trajectories of some products cannot be accurately fitted by the first two models, such as the GaAs laser [4] and the electrical connector [5]. Recently, as a stochastic process with independent increment, the IG process has been increasingly applied in the modeling of degradation. Ye [6] defined the IG process as a composite Poisson process with independent but not necessarily identical distribution in the extreme case, which compensates the physical significance of the IG process. Wang et al. [4] described the degradation characteristics of the products through the IG process and solved the problem of parameter estimation by the maximum likelihood estimation (MLE) method. Peng [7] used the Bayesian method to optimize the parameter estimation in IG process. Then Peng [8] took the individual differences of the similar products into consideration, solved the parameter estimation problem by the EM algorithm, and gave a clear analytical prediction expression of the remaining life. Aiming at handling the impact of catastrophe data on prediction, Xu [9] proposed a self-adaptive IG process which has a great significance on improving the accuracy of prediction.

Many scholars have done a lot of research on the accelerated degradation test using IG process. Wang [10] optimized the generalized IG process under the SSADT. Based on the invariance principle of acceleration coefficient, Wang [11] proposed an accelerated degradation modeling method using IG process with stochastic parameters; Ye [5] represented the individual differences in IG process through the randomization of parameters and optimized the CSADT design; Wu [12] proposed a multiobjective optimization design method and optimized CSADT using IG process. The measurement errors have not been investigated in those studies; however, that affects the fitting degree of the model and the accuracy of life estimation [13]. Therefore, Li [14] integrated the error item into the modeling in SSADT optimization design using IG process; however, the possible influence by the individual differences was neglected in Li’s study. Furthermore, in ADT, the individual differences between similar products and measurement uncertainties are ubiquitous and they considerably affect the accuracy of life prediction [15, 16]. Hence, it is necessary to establish an IG process degradation model with the consideration of both measurement errors and individual differences.

The paper is organized as follows. Section 2 introduces the IG process degradation model considering both measurement errors and individual differences. The acceleration models of the constant stress and the step stress using IG process are introduced in Section 3. In Section 4, a MLE method based on Monte Carlo integral and Genetic Algorithm (GA) is proposed, which solves the problem of complex integrals and unobservable measurement errors effectively in parameter estimation process. In Section 5, the efficiency and reasonability of the proposed approach are validated via Monte Carlo simulation study by two examples: the CSADT and the SSADT. And the proposed approach is validated by an illustrative example of electrical connector in Section 6. The conclusions are finally drawn in Section 7.

In this section, we adopt the IG model with stochastic parameters to characterize individual differences, with the consideration of the effect of measurement errors through incorporating error terms, and establish a degradation model for products.

As a newly introduced stochastic process model, the IG model shows a good fitting characteristic for the degradation process of many products [4, 5]. According to the previous description, the measurement errors of the degraded data are inevitable due to the interference of various noises during the degradation process; therefore, it is necessary to be considered in the degradation modeling. Set to be the observed degradation data and the real degradation data of product performance. In this paper, the increment of degradation data is used as the modeling object. And the degradation model can be expressed as

where and denote the increment of the degradation inspections and the true degradation of the product's performance characteristic at time t, respectively. is the increment of measurement errors which follows a normal distribution as . The follows the IG distribution, and it is expressed as . The probability density function (PDF) of and can be defined as

where is the time scale function and denotes the increment of time scale function. Under the condition that the failure threshold is , it is considered that the first hitting time (FHT) of IG process is product life T:

Therefore, the cumulative distribution function (CDF) of the life for the IG process can be obtained [4]

The corresponding PDF is

To reflect the individual differences between the similar products, Wang and Xu [17] proposed a method to incorporate the random effects in the Inverse Gaussian process by letting , with PDF

where denotes a gamma function, and (3) can be considered as the conditional distribution under the condition of being a constant; the expression of can be calculated by , then

Based on the concept of FHT, the CDF of IG process with random effects can be obtained [4]

where is the student t-distribution with degrees of freedom. The corresponding PDF can be obtained as follows:

From the above exposition, in the case of rated stress, the reliability analysis and the life estimation of the product degradation process described by IG process can be carried out through the analytical solutions.

The acceleration model can reflect the effect of the stress change upon the degradation process by establishing the relationship between accelerated stress and model parameters. Considering that the mean of the IG process is and the variance is , we assume that is a stress dependent function and ignore the effect of stress on the parameter in order to reflect the effect of stress upon the speed and the volatility of degradation simultaneously. This section takes the temperature stress as an example to illustrate the relationship between the parameters and the accelerated stress by the Arrhenius model

where and are undetermined parameters, stress level . is the accelerated stress after standardization, which is defined as

where denotes the normal working stress, is the maximum accelerated stress, .

The CSADT and the SSADT have been widely used in all kinds of product tests, such as LED [18], electrical connector [17], and SLD [19]. We suppose that there are units participating in ADT from the beginning of the experiment to cut-off time .

In the CSADT, let denote m higher stress levels, denote the degradation path of the unit under , and , .

As for SSADT, assume that there are units in test, where each of them is tested under a total number of stress levels , which can be expressed as

Under the stress level , each unit has been monitored for times, . Suppose that the monitoring intervals for each unit is , and , denotes the complete degradation path under SSADT, the relation between and can be given as (14).

4. Statistical Inference

In this paper, the generalized IG process is adopted. For the convenience of parameter estimation, set time function and the model parameters . Since the measurement errors are integrated into model, the solution of the MLE becomes complicated. On the basis of considering individual difference and measurement errors simultaneously, a MLE method based on combination of Monte Carlo integral and GA for the IG process in ADT is proposed. The SSADT is taken as an example to make statistical inference, and the CSADT can draw corresponding calculation results with reference to it.

Let denote the increment of degradation inspections of the unit at time under the stress and denote the increment of the true degradation, and is the increment of measurement errors which is independent and identically distributed. By (1) we can obtain

For the SSADT, denote the complete degradation path under SSADT of units and denote the total measurement times of a single unit, . And the expression of can be obtained as follows:

contains the observation error, which makes and follow different distributions. By (15), let denote and follow IG distribution with random effects. The expression of the maximum likelihood function can be obtained by the full probability formula

During each measurement, the absolute value of each measurement error should be less than the degradation inspection. Thus, the integral bounds of in (17) is . Considering that the direct integral of (17) is complicated, we adopt a simplified method referring to [14] as follows.

Define a two-variable function

can be simplified as

Formula (19) converts into the expectation of the two-variables function to . Then, an approximate calculation is conducted as follows:

In order to calculate the parameter configuration in the case of minimum , the MLE of unknown parameters is proposed based on Monte Carlo integral and Genetic Algorithm:(1)The upper limit of is calculated by formulation (21) and setting the step length to be , a series of can be obtained:(2)For each , simulate times by normal distribution random number. By substituting into (20), the approximate maximum likelihood function expression can be obtained as follows:(3)In the maximum likelihood function, there are no other unknown parameters except the remaining unknown parameters . In this case, the parameter estimation is regarded as an unconstrained minimum optimization problem. Since the maximum likelihood estimation requires to be the maximum, here is regarded as an objective function and as the variable. By calling the gatool in Matlab, the minimum and the corresponding can be obtained by setting parameters ranges, population sizes, and iteration times in the corresponding item columns. A series of and its corresponding can be obtained by performing the operations above under all values of . According to the definition of maximum likelihood estimation, the minimum can be obtained by comparison. The corresponding and are the optimal parameters for the estimation. The flow chart of parameter estimation is shown in Figure 1.

The approximated value of is given by Monte Carlo method, and it is substituted into the expression of maximum likelihood function, which solves the unknown . The optimal values of other parameters are obtained through GA, and the parameter including is selected by comparing the value of a series of . Finally, the unknown parameters estimation of model is realized. Here are three points for explanation.(i)In step (1), for all the units within two adjacent intervals, the variance of is less than variance of because the fluctuation of includes the random effect, individual difference of and the effects of .(ii)In step (3), as the objective function, the maximum likelihood function is searched, and its numerical value is simplified, so there is a certain deviation from the true value, which is only as the basis for the selection of the optimal parameter . The actual value is calculated by substituting the optimal parameter into (17).(iii)In the searching process, since is carried out under the condition that is determined and the variable is only , is denoted as in the description above. However, since is also a variable in the whole process of parameter estimation, we finally express as . There is no difference in numerical values between and .

5. Simulation Study

In order to verify the correctness of the proposed modeling method and the effectiveness of the parameter estimation method, the simulation and verification of the IG process under CSADT and SSADT are carried out, respectively. For the convenience of following description, M0 is used as the IG model with the consideration of both individual differences and measurement errors; M1 is used as the IG model with the consideration of individual differences proposed by Ye in [5], and M2 as the IG model with the consideration of measurement errors proposed by Ye in [14].

In this section a simulation study for the CSADT is conducted. Assume that the observation times of each unit are equal under each accelerated stress, and the time intervals are also equal. 15 groups of degradation data are simulated under three accelerated stresses of 60°C, 80°C, and 100°C, and the working stress is 40°C. The failure threshold and the individual differences are simulated by parameter randomization. The interval is 100h within 12 measurements. The degradation data and the measurement time are shown in Table 1.

 Stress Time/h 100 200 300 400 500 600 700 800 900 1000 1100 1200 60°C 4.241 6.918 8.893 10.505 11.822 13.143 14.238 15.626 16.968 18.094 19.168 20.404 4.527 6.785 8.602 9.979 11.501 12.690 13.956 15.364 16.541 17.671 18.742 19.750 4.491 7.814 9.805 11.338 12.469 13.713 14.502 15.258 16.637 18.604 19.327 20.177 4.218 6.574 8.296 10.041 11.380 12.536 13.866 14.833 16.276 17.569 18.831 19.642 4.329 6.459 8.355 9.767 10.968 12.275 13.572 14.529 15.689 16.943 17.899 19.035 80°C 6.276 9.669 11.771 14.489 16.406 18.583 20.625 22.646 23.863 25.493 26.976 28.476 6.046 9.142 11.719 14.331 16.688 18.725 20.720 22.634 24.066 25.664 27.432 28.888 6.625 10.434 13.049 15.200 17.016 18.922 21.207 23.144 24.970 26.357 27.922 29.478 6.138 9.736 11.984 14.003 16.379 18.483 20.179 21.423 23.084 24.742 26.415 27.524 4.314 9.829 13.640 16.848 19.404 21.811 23.550 25.507 26.828 28.579 29.444 30.861 100°C 7.618 11.411 14.403 17.291 19.591 22.277 24.973 26.862 28.936 31.513 33.204 35.119 7.820 11.702 14.649 17.664 20.021 21.813 23.861 25.671 26.745 29.113 30.272 35.294 7.574 11.425 15.152 19.381 22.349 24.883 27.815 29.779 31.376 33.675 35.321 37.083 6.830 10.088 12.242 15.564 18.013 21.419 23.404 25.524 27.016 29.843 31.435 32.145 7.806 11.775 15.782 18.738 20.897 23.960 26.297 27.949 29.920 31.837 33.655 35.421

As for the three IG models mentioned above, the estimated parameters and real values are separately given in Table 2. Introducing the Akaike Information Criterion (AIC) as the standard to evaluate model fitting degree, and combining with the maximum value of log likelihood function (Log-LF), the three IG models mentioned above can be evaluated. The calculation formula of AIC is given as

 Model Log-LF AIC MTTF/h Real -2.172 1.342 1.415 2.716 0.623 0.050 -185.719 383.438 10362.514 M0 -2.284 1.456 1.225 2.515 0.517 0.042 -195.432 402.864 10312.417 M1 -1.912 1.216 1.514 2.351 0.549 — -200.241 410.482 10192.347 M2 -2.341 1.160 — — 0.721 0.054 -205.329 420.658 10476.214

where denotes the number of undetermined parameters in model and the AIC values of the three IG models are shown in Table 2. To reflect the accuracy of the proposed method in life estimation, the mean time to failure (MTTF) of the models is also listed in Table 2. For proving the correctness of the parameter estimation method proposed above, the total mean square errors (TMSE) of parameter estimated values to real values tested in 500 and 1000 times simulation are respectively shown in Table 3, and the calculation formula of TMSE is given as

 Numbers 500 0.092 0.024 0.035 0.049 0.103 0.063 1000 0.053 0.009 0.012 0.031 0.062 0.026

where denotes the true values of parameters, denotes the estimated values of parameters, and denotes the numbers of simulation.

Table 2 shows that the AIC value of M0 is smaller than the others, which indicates a good fitness. The MTTF is closer to the real value, reflecting its high precision in life estimation. From Table 3, it can be observed that the proposed parameter estimation method has high precision under the results of multiple simulations and can be used to estimate the parameters of proposed model reliably. The reliability curve of each model is shown in Figure 2. It can be observed that the fitting degree of M0 to real model is the highest, and this proves the correctness of proposed method in CSADT.

The simulation study for SSADT is conducted in this section. Assume that units are tested in three accelerated stresses at 60°C, 80°C, and 100°C in sequence and the working stress is 40°C. The failure threshold and the measurement times of each unit are equal at each accelerated stress. The measurement interval is the same, 100h. 8 groups of the degenerate data are simulated in SSADT, and the measurement times are set to be 12 times, 4 times for each stress. The degradation data and the measurement times are shown in Table 4.

 Number Time/h 100 200 300 400 500 600 700 800 900 1000 1100 1200 1 3.810 6.114 8.102 9.399 13.245 15.609 17.547 18.807 22.559 24.889 26.919 28.205 2 4.806 7.028 8.936 10.471 15.218 17.370 19.231 20.684 25.538 27.801 29.664 31.129 3 4.437 6.312 7.985 9.758 14.175 16.139 17.856 19.712 24.078 25.995 27.625 29.443 4 4.223 6.170 8.851 10.568 14.750 16.608 19.131 20.794 24.923 26.851 29.522 31.214 5 4.053 6.824 9.631 10.831 14.880 17.652 20.323 21.475 25.584 28.393 31.043 32.201 6 4.675 6.689 8.557 10.116 14.826 16.888 18.799 20.312 24.962 26.937 28.719 30.228 7 3.938 6.082 8.075 9.642 13.606 15.811 17.792 19.352 23.227 25.462 27.431 29.011 8 4.133 6.349 8.330 10.001 14.131 16.215 18.135 19.838 23.896 26.120 28.052 29.788

By adopting the parameter estimation method proposed above, the results of the parameter estimation of simulation data in SSADT are obtained and shown in Table 5, and compared with the other two models. According to the analysis of Log-LF, AIC, and MTTF, it can be observed that the parameter estimation precision of model M0 is higher and the value of AIC is lower, which reflects a better fitness under the consideration of individual difference and measurement error. For M1, since the measurement errors cannot be considered, the parameter estimation precision of it is reduced and the goodness of fit is lower than that of model M0. Since in the model M2 only the measurement errors are considered and it cannot randomize the parameters, its accuracy of life estimation and fitness are reduced. This is the same as the results from the simulation study of CSADT. It can be observed from the MTTF of the three models that the life estimation of M0 is closer to the one from the real model compared to M1 and M2. The same results can be observed from the reliability curves of three models, as shown in Figure 3.

 Model Log-LF AIC MTTF/h Real -2.172 1.342 1.415 2.716 0.623 0.050 -184.951 381.902 10362.514 M0 -2.041 1.423 1.191 2.842 0.584 0.046 -190.514- 393.028 10410.325 M1 -2.304 1.418 1.251 2.638 0.715 — -198.230 406.460 10223.514 M2 -2.410 1.806 — — 0.493 0.068 -204.244 416.488 10537.241

6. Illustrative Example

Stress relaxation is the phenomenon that the strain on a component remains unchanged while the stress gradually decreases. With the increase of working temperature, the contact of electrical connectors often fails due to stress relaxation. Yang [20] verified that the IG process has a good fitting effect in the degradation test of electrical connectors. In this section, the case of stress relaxation test [20] is used for modeling analysis. Set the time function and the failure threshold . The normal working temperature stress is 40°C, and three accelerated temperature stresses are 65°C, 85°C, and 100°C, respectively. The experimental data and the measurement time are shown in Tables 6 and 7. The degradation trend of stress relaxation is shown in Figure 4.

 Temperature ID Stress Relaxation/% 65°C 1 2.12 2.7 3.52 4.25 5.55 6.12 6.75 7.22 7.68 8.46 9.46 2 2.29 3.24 4.16 4.86 5.74 6.85 7.4 8.14 9.25 10.55 3 2.4 3.61 4.35 5.09 5.5 7.03 8.24 8.81 9.629 10.27 11.11 4 2.31 3.48 5.51 6.2 7.31 7.96 8.57 9.07 10.46 11.48 12.31 5 3.14 4.33 5.92 7.22 8.14 9.07 9.44 10.09 11.2 12.77 13.51 6 3.59 5.55 5.92 7.68 8.61 10.37 11.11 12.22 13.51 14.16 15 85°C 7 2.77 4.62 5.83 6.66 8.05 10.61 11.2 11.98 13.33 15.64 8 3.88 4.37 6.29 7.77 9.16 9.9 10.37 12.77 14.72 16.8 9 3.18 4.53 6.94 8.14 8.79 10.09 11.11 14.72 16.47 18.66 10 3.61 4.37 6.29 7.87 9.35 11.48 12.4 13.7 15.37 18.51 11 3.42 4.25 7.31 8.61 10.18 12.03 13.7 15.27 17.22 19.25 12 5.27 5.92 8.05 9.81 12.4 13.24 15.83 17.59 20.09 23.51 100°C 13 4.25 5.18 8.33 9.53 11.48 13.14 15.15 16.94 18.05 19.44 14 4.81 6.16 7.68 9.25 10.37 12.4 15 16.2 18.24 20.09 15 5.09 7.03 8.33 10.37 12.22 14.35 16.11 18.7 19.72 21.66 16 4.81 7.5 9.16 10.55 13.51 15.55 16.57 19.07 20.27 22.4 17 5.64 6.57 8.61 12.5 14.44 16.57 18.7 21.2 22.59 24.07 18 4.72 8.14 10.18 12.4 15.09 17.22 19.16 21.57 24.35 26.2
 Temperature Measurement time epochs/h 65°C 108 241 534 839 1074 1350 1637 1890 2178 2513 2810 85°C 46 108 212 408 632 764 1011 1333 1517 2586 100°C 46 108 212 344 446 626 729 927 1005 1218

For validating the superiority of the proposed method in the modeling of stress relaxation of electrical connectors, the three proposed models are evaluated by a combination of Log-LF and AIC, as shown in Table 8. And to validate goodness of fit of IG model for the stress relaxation of electrical connectors, the quantile-quantile plot (Q-Q plot) is shown in Figures 57.

 Model Log-LF AIC M0 -2.261 1.715 2.152 3.816 0.449 0.083 -170..373 352.746 M1 -2.291 1.868 2.022 3.761 0.492 — -200.041 410.082 M2 -2.518 2.272 — — 0.459 0.085 -180.725 371.45

From Table 8, it can be observed that comparing with the models M1 and M2, M0 has a higher Log-LF and a smaller AIC in the three stress level CSADT, which indicates that the model considering measurement errors and individual differences have better fitting performance than the other two under three stress levels. Through the comparison of Q-Q plot, it can also be observed that the quantile of model M0 displays a good linearity to test data, which reflects a good fitness compared with the other two models under three stress levels. Thus, it is meaningful to consider individual differences and measurement errors simultaneously.

Through the reliability curve of three IG models, the advantage of the proposed method in reliability analysis can be seen more intuitively, as shown in Figure 8. It can be observed that the model M1 evaluates reliability conservatively compared with M0, and the maintenance time is estimated early, which leads to lower utilization of product. The model M2 has a high reliability but may lead to the lag of maintenance and increase the risk of early failure of product. Therefore, it is reasonable to evaluate the reliability of products via the model M0.

7. Conclusion

In this paper, considering the influence of individual differences and measurement errors of life estimation in ADT, the errors are incorporated into the IG process with stochastic parameters. The IG degradation model considering individual differences and measurement errors simultaneously is proposed in ADT. Considering the characteristics of the unobservable measurement errors and the complex integral problem in the process of MLE, a MLE method is proposed by combining the Monte Carlo integral with GA, and the method has been verified to be accurate and reliable in the simulation study. The simulations study towards CSADT and SSADT have verified the correctness of the proposed model considering individual differences and measurement errors in the fitting degree of the degradation trajectory and the life estimation. And the illustrative example of electrical connector in CSADT has also demonstrated the efficiency and reliability of the proposed model.

Based on the research in this paper, the following points should be further studied: (1) it is worth investigating the analytical expression of life estimation using IG process considering individual differences and measurement errors simultaneously, based on the FHT concept; more details are given in [21]; (2) it may be of interest to reduce the effect of measurement errors by the filtering method; more details can be found in [13]; (3) it is significant to establish an effective maintenance plan based on proposed model in actual engineering; (4) a further study should focus on accelerated degradation tests considering measurement errors when both the mean and scale parameters are stress dependent; (5) considering that there are still some deviations between the estimated parameters and the real parameters of the model, it may be of interest to consider the interval estimation of parameters.

Data Availability

The simulation study and illustrative example data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors are grateful for the financial support from the National Natural Science Foundation of China [No. 51875499] and Innovation Funding for Postgraduates in Hebei Province [No. CXZZSS2018035, CXZZBS2017040].

References

1. L. Hao, L. Bian, N. Gebraeel, and J. Shi, “Residual life prediction of multistage manufacturing processes with interaction between tool wear and product quality degradation,” IEEE Transactions on Automation Science and Engineering, vol. 14, no. 2, pp. 1211–1224, 2017. View at: Publisher Site | Google Scholar
2. X. S. Si, C. H. Hu, X. Y. Kong, and D. H. Zhou, “A residual storage life prediction approach for systems with operation state switches,” IEEE Transactions on Industrial Electronics, vol. 61, no. 11, pp. 6304–6315, 2014. View at: Publisher Site | Google Scholar
3. R. Zhou, N. Serban, and N. Gebraeel, “Degradation-based residual life prediction under different environments,” The Annals of Applied Statistics, vol. 8, no. 3, pp. 1671–1689, 2014. View at: Publisher Site | Google Scholar | MathSciNet
4. X. Wang and D. Xu, “An inverse Gaussian process model for degradation data,” Technometrics, vol. 52, no. 2, pp. 188–197, 2010. View at: Publisher Site | Google Scholar | MathSciNet
5. Z.-S. Ye, L.-P. Chen, L. C. Tang, and M. Xie, “Accelerated degradation test planning using the inverse gaussian process,” IEEE Transactions on Reliability, vol. 63, no. 3, pp. 750–763, 2014. View at: Publisher Site | Google Scholar
6. Z. S. Ye and N. Chen, “The inverse Gaussian process as a degradation model,” Technometrics, vol. 56, no. 3, pp. 302–311, 2014. View at: Publisher Site | Google Scholar
7. W. Peng, Y. F. Li, Y. J. Yang, H. Z. Huang, and M. J. Zuo, “Inverse Gaussian process models for degradation analysis: a Bayesian perspective,” Reliability Engineering and System Safety, vol. 130, no. 1, pp. 175–189, 2014. View at: Publisher Site | Google Scholar
8. C.-Y. Peng, “Inverse Gaussian processes with random effects and explanatory variables for degradation data,” Technometrics, vol. 57, no. 1, pp. 100–111, 2015. View at: Publisher Site | Google Scholar | MathSciNet
9. W. Xu and W. Wang, “RUL estimation using an adaptive inverse Gaussian model,” in Proceedings of the IEEE Conference on Prognostics and System Health Management, vol. 33, pp. 331–336, IEEE, 2013. View at: Google Scholar
10. H. Wang, G. J. Wang, and F. J. Duan, “Planning of step-stress accelerated degradation test based on the inverse gaussian process,” Reliability Engineering and System Safety, vol. 154, pp. 97–105, 2016. View at: Publisher Site | Google Scholar
11. H. Wang, K. Teng, and W. J. Xi, “Accelerated degradation modeling method based on Inverse Gaussian processes with random parameters,” Journal of Beijing University of Aeronautics and Astronautics, vol. 42, no. 9, pp. 1843–1850, 2016. View at: Google Scholar
12. Z. Y. Wu, X. P. Liu, B. Guo, D. J. Cui, and L. J. Zhang, “Multi-objective optimization design of constant stress accelerated degradation test using inverse gaussian process,” IEEE Access, vol. 7, pp. 28757–28768, 2019. View at: Publisher Site | Google Scholar
13. X. S. Si, C. H. Hu, Q. Zhang, F. H. He, and T. Zhou, “Estimating remaining useful life under uncertain degradation measurements,” Acta Electronica Sinica, vol. 43, no. 1, pp. 30–35, 2015. View at: Google Scholar
14. S. Li, Z. Chen, and E. Pan, “Step-stress accelerated degradation test plan for generalized inverse gaussian process,” Journal of Shanghai Jiaotong University, vol. 51, no. 2, pp. 186–192, 2017. View at: Google Scholar
15. J. Li, Z. Wang, X. Liu, Y. Zhang, H. Fu, and C. Liu, “A Wiener process model for accelerated degradation analysis considering measurement errors,” Microelectronics Reliability, vol. 65, pp. 8–15, 2016. View at: Publisher Site | Google Scholar
16. X.-P. Liu, L.-J. Zhang, K.-K. Shen, and Q. Gao, “Step stress accelerated degradation test modeling and remaining useful life estimation in consideration of measuring error,” Acta Armamentarii, vol. 38, no. 8, pp. 1586–1592, 2017. View at: Google Scholar
17. H. Wang, Y. Zhao, X. Ma, and H. Wang, “Optimal design of constant-stress accelerated degradation tests using the M-optimality criterion,” Reliability Engineering & System Safety, vol. 164, pp. 45–54, 2017. View at: Publisher Site | Google Scholar
18. M. Cai, D. Yang, K. Tian et al., “Step-stress accelerated testing of high-power LED lamps based on subsystem isolation method,” Microelectronics Reliability, vol. 55, no. 9-10, pp. 1784–1789, 2015. View at: Publisher Site | Google Scholar
19. L. Wang, R. Pan, X. Li, and T. Jiang, “A Bayesian reliability evaluation method with integrated accelerated degradation testing and field information,” Reliability Engineering & System Safety, vol. 112, pp. 38–47, 2013. View at: Publisher Site | Google Scholar
20. G. B. Yang, “System reliability evaluation and allocation,” in Life Cycle Reliability Engineering, pp. 94-95, Wiley, New Jersey, NJ, USA, 2007. View at: Google Scholar
21. D. Pan, Y. Wei, H. Fang, and W. Yang, “A reliability estimation approach via Wiener degradation model with measurement errors,” Applied Mathematics and Computation, vol. 320, pp. 131–141, 2018. View at: Publisher Site | Google Scholar

Copyright © 2019 Xiaoping Liu 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.