Inference from Accelerated Degradation and Failure Data Based on Gaussian Process Models

An important problem in reliability and survival analysis is that of modeling degradation together with any observed failures in a life test. Here, based on a continuous cumulative damage approach with a Gaussian process describing degradation, a general accelerated test model is presented in which failure times and degradation measures can be combined for inference about system lifetime. Some specific models when the drift of the Gaussian process depends on the acceleration variable are discussed in detail. Illustrative examples using simulated data as well as degradation data observed in carbon-film resistors are presented.     

Modelling Accelerated Degradation Data Using Wiener Diffusion With A Time Scale Transformation

Engineering degradation tests allow industry to assess the potential life span of long-life products that do not fail readily under accelerated conditions in life tests. A general statistical model is presented here for performance degradation of an item of equipment. The degradation process in the model is taken to be a Wiener diffusion process with a time scale transformation. The model incorporates Arrhenius extrapolation for high stress testing. The lifetime of an item is defined as the time until performance deteriorates to a specified failure threshold. The model can be used to predict the lifetime of an item or the extent of degradation of an item at a specified future time. Inference methods for the model parameters, based on accelerated degradation test data, are presented. The model and inference methods are illustrated with a case application involving self-regulating heating cables. The paper also discusses a number of practical issues encountered in applications. This revised version was published online in July 2006 with corrections to the Cover Date.     

Optimal design of accelerated degradation tests based on Wiener process models

Optimal accelerated degradation test (ADT) plans are developed assuming that the constant-stress loading method is employed and the degradation characteristic follows a Wiener process. Unlike the previous works on planning ADTs based on stochastic process models, this article determines the test stress levels and the proportion of test units allocated to each stress level such that the asymptotic variance of the maximum-likelihood estimator of the qth quantile of the lifetime distribution at the use condition is minimized. In addition, compromise plans are also developed for checking the validity of the relationship between the model parameters and the stress variable. Finally, using an example, sensitivity analysis procedures are presented for evaluating the robustness of optimal and compromise plans against the uncertainty in the pre-estimated parameter value, and the importance of optimally determining test stress levels and the proportion of units allocated to each stress level are illustrated.     

Inverse gaussian accelerated test models based on cumulative damage

The failure of a system under environmental stress often can be described by an accelerated test model which incorporates the environmental variable L. Here, the failure of such a system at environmental level L is modeled as the first passage of accumulated damage to a critical threshold value. Assuming a discrete additive damage model leads to a Birnbaum–Saunders-type distribution for the failure time which can be closely approximated by an inverse Gaussian-type model. However, if a continuous damage model based on a Gaussian process is assumed, a more general family of inverse Gaussian accelerated test models is obtained. Three sets of failure data are discussed to illustrate the usefulness of this general family.     

Failure Inference and Optimization for Step Stress Model Based on Bivariate Wiener Model

In this article, we consider the situation under a life test, in which the failure time of the test units are not related deterministically to an observable stochastic time varying covariate. In such a case, the joint distribution of failure time and a marker value would be useful for modeling the step stress life test. The problem of accelerating such an experiment is considered as the main aim of this article. We present a step stress accelerated model based on a bivariate Wiener process with one component as the latent (unobservable) degradation process, which determines the failure times and the other as a marker process, the degradation values of which are recorded at times of failure. Parametric inference based on the proposed model is discussed and the optimization procedure for obtaining the optimal time for changing the stress level is presented. The optimization criterion is to minimize the approximate variance of the maximum likelihood estimator of a percentile of the products’ lifetime distribution.     

A general class of cumulative damage models for materials failure

A generalized cumulative damage approach is presented which yields a large family of accelerated test inverse Gaussian-type models for strength of materials that incorporate the size effect as the acceleration variable. Previous models are generalized here in three aspects: the cumulative damage model is more general and can include damage functions other than the additive and multiplicative damages; the strength reduction function due to initial damage existing in the material is taken to be a very general function; and the initial damage process is a more general stochastic process that includes those previously assumed as special cases. The approach taken here is therefore the most general cumulative damage model obtained to date and yields a large number of potentially more useful accelerated test models for material strength. Estimation of model parameters by maximum likelihood methods is discussed, and two examples using real tensile strength data for carbon micro-composites and single carbon fibers are presented, illustrating the improvement of the new approach over previous models.     

Bivariate Constant-Stress Accelerated Degradation Model and Inference

To assess the reliability of highly reliable products that have two or more performance characteristics (PCs) in an accurate manner, relations between the PCs should be taken duly into account. If they are not independent, it would then become important to describe the dependence of the PCs. For many products, the constant-stress degradation test cannot provide sufficient data for reliability evaluation and for this reason, accelerated degradation test is usually performed. In this article, we assume that a product has two PCs and that the PCs are governed by a Wiener process with a time scale transformation, and the relationship between the PCs is described by the Frank copula function. The copula parameter is dependent on stress and assumed to be a function of stress level that can be described by a logistic function. Based on these assumptions, a bivariate constant-stress accelerated degradation model is proposed here. The direct likelihood estimation of parameters of such a model becomes analytically intractable, and so the Bayesian Markov chain Monte Carlo (MCMC) method is developed here for this model for obtaining the maximum likelihood estimates (MLEs) efficiently. For an illustration of the proposed model and the method of inference, a simulated example is presented along with the associated computational results.     

A class of accelerated means regression models for recurrent event data

In this article, we propose a general class of accelerated means regression models for recurrent event data. The class includes the proportional means model, the accelerated failure time model and the accelerated rates model as special cases. The new model offers great flexibility in formulating the effects of covariates on the mean functions of counting processes while leaving the stochastic structure completely unspecified. For the inference on the model parameters, estimating equation approaches are developed and both large and final sample properties of the proposed estimators are established. In addition, some graphical and numerical procedures are presented for model checking. An illustration with multiple-infection data from a clinic study on chronic granulomatous disease is also provided.     

A tampered Brownian motion process model for partial step-stress accelerated life testing

In partial step-stress accelerated life testing, models extrapolating data obtained under more severe conditions to infer the lifetime distribution under normal use conditions are needed. Bhattacharyya (Invited paper for 46th session of the ISI, 1987) proposed a tampered Brownian motion process model and later derived the probability distribution from a decay process perspective without linear assumption. In this paper, the model is described and the features of the failure time distribution are discussed. The maximum likelihood estimates of the parameters in the model and their asymptotic properties are presented. An application of models for step-stress accelerated life test to fields other than engineering is described and illustrated by applying the tampered Brownian motion process model to data taken from a clinical trial.     

Mis-specification analyses of gamma and Wiener degradation processes

Degradation models are widely used these days to assess the lifetime information of highly reliable products if there exist some quality characteristics (QC) whose degradation over time can be related to the reliability of the product. In this study, motivated by a laser data, we investigate the mis-specification effect on the prediction of product's MTTF (mean-time-to-failure) when the degradation model is wrongly fitted. More specifically, we derive an expression for the asymptotic distribution of quasi-MLE (QMLE) of the product's MTTF when the true model comes from gamma degradation process, but is wrongly assumed to be Wiener degradation process. The penalty for the model mis-specification can then be addressed sequentially. The result demonstrates that the effect on the accuracy of the product's MTTF prediction strongly depends on the ratio of critical value to the scale parameter of the gamma degradation process. The effects on the precision of the product's MTTF prediction are observed to be serious when the shape and scale parameters of the gamma degradation process are large. We then carry out a simulation study to evaluate the penalty of the model mis-specification, using which we show that the simulation results are quite close to the theoretical ones even when the sample size and termination time are not large. For the reverse mis-specification problem, i.e., when the true degradation is a Wiener process, but is wrongly assumed to be a gamma degradation process, we carry out a Monte Carlo simulation study to examine the effect of the corresponding model mis-specification. The obtained results reveal that the effect of this model mis-specification is negligible.     

Reliability modeling of systems with two dependent degrading components based on gamma processes

Abstract

Many engineering systems have multiple components with more than one degradation measure which is dependent on each other due to their complex failure mechanisms, which results in some insurmountable difficulties for reliability work in engineering. To overcome these difficulties, the system reliability prediction approaches based on performance degradation theory develop rapidly in recent years, and show their superiority over the traditional approaches in many applications. This paper proposes reliability models of systems with two dependent degrading components. It is assumed that the degradation paths of the components are governed by gamma processes. For a parallel system, its failure probability function can be approximated by the bivariate Birnbaum–Saunders distribution. According to the relationship of parallel and series systems, it is easy to find that the failure probability function of a series system can be expressed by the bivariate Birnbaum–Saunders distribution and its marginal distributions. The model in such a situation is very complicated and analytically intractable, and becomes cumbersome from a computational viewpoint. For this reason, the Bayesian Markov chain Monte Carlo method is developed for this problem that allows the maximum likelihood estimates of the parameters to be determined in an efficient manner. After that, the confidence intervals of the failure probability of systems are given. For an illustration of the proposed model, a numerical example about railway track is presented.     

A Bivariate Gamma Model for a Latent Degradation Process

A model is presented in this article based on a bivariate gamma process in which, the first component is latent and determines the failure time and the second represents the marker. This process is a more realistic model for a degradation process. After introducing the model, we obtain failure and survival probability distributions and discuss parametric and predictive inferences based on the Maximum Likelihood method and in a Bayesian setup.     

Optimal Design for Step-Stress Accelerated Degradation Test with Multiple Performance Characteristics Based on Gamma Processes

Step-stress accelerated degradation test (SSADT) plays an important role in assessing the lifetime distribution of highly reliable products under normal operating conditions when there are not enough test units available for testing purposes. Recently, the optimal SSADT plans are presented based on an underlying assumption that there is only one performance characteristic. However, many highly reliable products usually have complex structure, with their reliability being evaluated by two or more performance characteristics. At the same time, the degradation of these performance characteristics would be always positive and strictly increasing. In such a case, the gamma process is usually considered as a degradation process due to its independent and nonnegative increments properties. Therefore, it is of great interest to design an efficient SSADT plan for the products with multiple performance characteristics based on gamma processes. In this work, we first introduce reliability model of the degradation products with two performance characteristics based on gamma processes, and then present the corresponding SSADT model. Next, under the constraint of total experimental cost, the optimal settings such as sample size, measurement times, and measurement frequency are obtained by minimizing the asymptotic variance of the estimated 100 qth percentile of the product’s lifetime distribution. Finally, a numerical example is given to illustrate the proposed procedure.     

Modelling accelerated life test data by using a Bayesian approach

Summary. Because of the high reliability of many modern products, accelerated life tests are becoming widely used to obtain timely information about their time-to-failure distributions. We propose a general class of accelerated life testing models which are motivated by the actual failure process of units from a limited failure population with a positive probability of not failing during the technological lifetime. We demonstrate a Bayesian approach to this problem, using a new class of models with non-monotone hazard rates, the hazard model with potential scope for use far beyond accelerated life testing. Our methods are illustrated with the modelling and analysis of a data set on lifetimes of printed circuit boards under humidity accelerated life testing.     

Multiple-Steps Step-Stress Accelerated Degradation Modeling Based on Wiener and Gamma Processes

During the step-stress accelerated degradation test (SSADT) experiment, the operator usually elevates the stress level at a predetermined time-point for all tested products that had not failed. This time-point is determined by the experience of the operator and the test is carried on until the performance degradation value of the product crosses the threshold value. In fact, this mode only works when a lot of products have been used in the experiment. But in the SSADT experiment, the number of products is relatively few, and so the test control to the products should be done more carefully. Considering the differences in products, we think the time-point of elevating stress level varies randomly from product-to-product. We consider a situation in which when the degradation value crosses a pre-specified value, the stress level is elevated. Under the circumstances, the time at which the degradation path crosses the pre-specified value is uncertain, and so we may regard it as a random variable. We discuss multiple-steps step-stress accelerated degradation models based on Wiener and gamma processes, respectively, and we apply the Bayesian Markov chain Monte Carlo (MCMC) method for such analytically intractable models to obtain the maximum likelihood estimates (MLEs) efficiently and present some computational results obtained from our implementation.     

Optimal step-stress accelerated degradation test plans for inverse Gaussian process based on proportional degradation rate model

The purpose of this paper is to address the optimal design of the step-stress accelerated degradation test (SSADT) issue when the degradation process of a product follows the inverse Gaussian (IG) process. For this design problem, an important task is to construct a link model to connect the degradation magnitudes at different stress levels. In this paper, a proportional degradation rate model is proposed to link the degradation paths of the SSADT with stress levels, in which the average degradation rate is proportional to an exponential function of the stress level. Two optimization problems about the asymptotic variances of the lifetime characteristics' estimators are investigated. The optimal settings including sample size, measurement frequency and the number of measurements for each stress level are determined by minimizing the two objective functions within a given budget constraint. As an example, the sliding metal wear data are used to illustrate the proposed model.     

Reliability statistical analysis about products of Birnbaum–Saunders fatigue life distribution for life test and progressive stress accelerated life test

Birnbaum–Saunders fatigue life distribution is an important failure model in the probability physical methods. It is more suitable for describing the life rules of fatigue failure products than common life distributions such as Weibull distribution and lognormal distribution. Besides, it is mainly applied to analytical research about fatigue failure and degradation failure of electronic product performance. The characteristic properties such as numerical characteristics and image features of density function and failure rate function are studied for generalized BS fatigue life distribution GBS(α, β, m) in this paper. Then the point estimates and approximate interval estimates of parameters are proposed for generalized BS fatigue life distribution GBS(α, β, m), and the precision of estimates are investigated by Monte Carlo simulations. Finally, when the scale parameter satisfies inverse power law model, the failure distribution model is given for the products of two-parameter BS fatigue life distribution BS(α, β) under progressive stress accelerated life test according to the time conversion idea of famous Nelson assumption, and then the points estimates of parameters are given.     

Joint modeling of degradation and failure time data

This paper surveys some approaches to model the relationship between failure time data and covariate data like internal degradation and external environmental processes. These models which reflect the dependency between system state and system reliability include threshold models and hazard-based models. In particular, we consider the class of degradation–threshold–shock models (DTS models) in which failure is due to the competing causes of degradation and trauma. For this class of reliability models we express the failure time in terms of degradation and covariates. We compute the survival function of the resulting failure time and derive the likelihood function for the joint observation of failure times and degradation data at discrete times. We consider a special class of DTS models where degradation is modeled by a process with stationary independent increments and related to external covariates through a random time scale and extend this model class to repairable items by a marked point process approach. The proposed model class provides a rich conceptual framework for the study of degradation–failure issues.     

A Bayesian semiparametric accelerated failure time model for arbitrarily censored data with covariates subject to measurement error

A flexible Bayesian semiparametric accelerated failure time (AFT) model is proposed for analyzing arbitrarily censored survival data with covariates subject to measurement error. Specifically, the baseline error distribution in the AFT model is nonparametrically modeled as a Dirichlet process mixture of normals. Classical measurement error models are imposed for covariates subject to measurement error. An efficient and easy-to-implement Gibbs sampler, based on the stick-breaking formulation of the Dirichlet process combined with the techniques of retrospective and slice sampling, is developed for the posterior calculation. An extensive simulation study is conducted to illustrate the advantages of our approach.     

Residual life estimation based on bivariate non-stationary gamma degradation process

Due to the growing importance in maintenance scheduling, the issue of residual life (RL) estimation for some high reliable products based on degradation data has been studied quite extensively. However, most of the existing work only deals with one-dimensional degradation data, which may not be realistic in some cases. Here, an adaptive method of RL estimation is developed based on two-dimensional degradation data. It is assumed that a product has two performance characteristics (PCs) and that the degradation of each PC over time is governed by a non-stationary gamma degradation process. From a practical consideration, it is further assumed that these two PCs are dependent and that their dependency can be characterized by a copula function. As the likelihood function in such a situation is complicated and computationally quite intensive, a two-stage method is used to estimate the unknown parameters of the model. Once new degradation information of the product being monitored becomes available, random effects are first updated by using the Bayesian method. Following that, the RL at current time is estimated accordingly. As the degradation data information accumulates, the RL can be re-estimated in an adaptive manner. Finally, a numerical example about fatigue cracks is presented in order to illustrate the proposed model and the developed inferential method.