An Optimal Design of Accelerated Degradation Tests Based on Degradation Performance

Show more

1. Introduction

With the development of technology, more and more products have higher reliability and longer service life. In many cases, accelerated degradation testing is superior to accelerated life testing in assessing product life. As time goes on, accelerated degradation test has been widely used in reliability evaluation of various products.

There are many factors involved in accelerated degradation test. One of the important aspects is the optimal arrangement or design of the test. In recent years, in the research of experimental optimal design, Boulanger, Escobar [1] and Yu, Tseng [2] gave the cost function of the experiment for the general case. Park and Yum [3] used precision limits for general degradation tests to optimize the tests. Tseng and Yu [4] proposed a method to optimize the distribution of degradation test duration. Yang [5] and Wu, Chang [6] obtained the optimal test design by minimizing the asymptotic variance of quantile life estimation under the restriction of service cost for the constant stress accelerated degradation test. Assuming that the degradation rate obeys the inverse Weibull distribution, Yu and Tseng [7] designed a constant stress degradation test. Park and Yum [8] designed a step-by-step stress degradation test in which destructive degradation detection was used. Wang, Xu [9] and Ye, Chen [10] systematically introduced the degradation model of inverse Gaussian process and compared it with gamma process model. Sung and Yum [11] studied the optimal design of accelerated degradation test based on the Wiener process.

In this paper, our main purpose is to find a method based on degradation performance to plan the accelerated degradation tests. We discuss the relationship between the optimal design of a kind of accelerated degradation test and the selection of its degradation of critical value. On condition that there is a linear relationship between the parameter $\mu $ of distribution and the critical values C. A special pattern of accelerated degradation tests is presented by taking the critical values as accelerated variable. And the optimal test plan is obtained by minimizing the asymptotic variance of the MLE of the parameter of distribution.

The present paper is structured as follows. In Section 2, we will discuss the relationship between the parameters of failure distribution and the critical value level. In Section 3, we will establish the optimization design model based on degradation performance. Finally, a numerical example is presented to illustrate the procedures of the test plan in Section 4, and Section 5 concludes this paper.

2. The Critical Values

For some products with high reliability and long service life, when evaluating their reliability indexes, even if the environmental stress is increased, the failure cannot be or rarely can be observed in the normal time interval. At this time, we use degradation test. Select a performance attribute of the sample, which degrades slowly with time. When its index value exceeds (higher or lower than) a preset limit value, the sample is judged to be invalid. The preset limit is called the degradation threshold. According to reference [12], more failure data can be obtained by selecting the critical value level, and the test time can be shortened.

In this way, a key problem is how to evaluate the failure law of products with the failure data obtained at different critical value levels. This means that the relationship between the parameters of failure distribution and the critical value level should be found. According to literature [13] [14], for some products, the degradation path has a linear form:

${D}_{ik}\left(t\right)={\alpha}_{1k}+{\alpha}_{2k}t+{\alpha}_{3k}{S}_{i}+{\epsilon}_{ik}(t)$

where, ${D}_{ik}\left(t\right)$ is the degradation value of the kth sample at the ith acceleration stress to time t.

${\alpha}_{1k},{\alpha}_{2k},{\alpha}_{3k}$ are the coefficient of degradation, ${\epsilon}_{ik}\left(t\right)$ is the measurement error. Now we take the critical value of degradation as ${C}_{j}$, form ${D}_{ik}\left(t\right)={C}_{j}$, We can work out the failure time. Because of ${C}_{j}={\alpha}_{1k}+{\alpha}_{2k}t+{\alpha}_{3k}{S}_{i}+{\epsilon}_{ik}\left(t\right)$, we

obtained ${t}_{ijk}=\frac{{C}_{j}-{\alpha}_{1k}-{\alpha}_{3k}{S}_{i}}{{\alpha}_{2k}}$. This represents the failure time of the kth sample at the ith stress level at the jth critical value. Thus, in case of the numbers of invalid samples are m, the average life of the sample at the jth critical value and the ith stress level shall be:

$\begin{array}{c}{\mu}_{ij}=\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}{t}_{ijk}=\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}\frac{{C}_{j}-{\alpha}_{1k}-{\alpha}_{3k}{S}_{i}}{{\alpha}_{2k}}}}\\ =-\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}\frac{{\alpha}_{1k}}{{\alpha}_{2k}}}-\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}\frac{{\alpha}_{3k}}{{\alpha}_{1k}}{S}_{i}}+\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}\frac{1}{{\alpha}_{1k}}{C}_{j}}\end{array}$

Set $\alpha =-\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}\frac{{\alpha}_{1k}}{{\alpha}_{2k}}}$, $\beta =-\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}\frac{{\alpha}_{3k}}{{\alpha}_{1k}}}$, $\gamma =\frac{1}{m}{\displaystyle \underset{k=1}{\overset{m}{\sum}}\frac{1}{{\alpha}_{1k}}}$.

Thus, we have ${\mu}_{ij}=\alpha +\beta \cdot {S}_{i}+\gamma \cdot {C}_{j}.$

Therefore, the average life is related to stress and critical value as follows:

$\mu =\alpha +\beta \cdot S+\gamma \cdot C$ (1)

In this way, as mentioned above, more failure data can be obtained by selecting the critical value level, and the test time can be shortened, so as to accelerate the test.

3. A special Optimization Design

3.1. Test Description and Assumption

Test description

The number of test samples is n. All samples are subjected to constant stress accelerated degradation test at q different but constant stress levels. The acceleration stress level (or its deformation) is ${S}_{1},{S}_{2},\cdots ,{S}_{q}$ ( ${S}_{0}<{S}_{1}<{S}_{2}<\cdots <{S}_{q}$, ${S}_{0}$ is the normal stress level， ${S}_{q}$ is the highest stress level). The ratio of samples is ${\pi}_{i}$ to the stress level ${S}_{i}$. So the Sample number is ${n}_{i}=n\cdot {\pi}_{i}$ at the stress level ${S}_{i}$. The test adopts timing truncation. The closing time is ${\tau}_{i}$. D represents the degradation amount of sample performance characteristics. We take m critical values for the degradation. Respectively ${C}_{1},{C}_{2},\cdots ,{C}_{m}$, where ${C}_{m}$ is an extreme threshold, ${C}_{0}$ represents the critical value of degradation under normal stress.

For discussion purposes, we first standardize the stress levels as follows: Set ${X}_{i}=\frac{{S}_{i}-{S}_{q}}{{S}_{0}-{S}_{q}},i=0,1,2,\cdots ,q$, thus, we obtain $1={X}_{0}>{X}_{1}>\cdots >{X}_{q}=0$. The

advantage is that the range of stress change is limited to the interval $\left[0,1\right]$. The level of the threshold is then standardized as follows:

${Y}_{j}=\frac{{C}_{j}-{C}_{m}}{{C}_{0}-{C}_{m}},j=1,2,\cdots ,m$, then ${Y}_{0}=1$, ${Y}_{m}=0$. And all ${Y}_{j}$ are between 0 and 1.

Test assumption

1) Failure time follows lognormal distribution $LN\left(\mu ,{\sigma}^{2}\right)$ ;

2) Parameter $\mu $ is a linear function of test stress and degradation critical value (from Equation (1)), hence:

$\mu =\mu \left({X}_{i},{Y}_{j}\right)=\alpha +\beta {X}_{i}+\gamma {Y}_{j}$ (2)

3) The parameters $\sigma $ do not change with the change of test stress and degradation critical value.

3.2. Optimization Design Model

In the above case, we hope to arrange the experiment based on some standards, so that the experiment can be optimized to a certain extent. Here, the experiment is arranged by minimizing the estimated asymptotic variance $AVar\left(\stackrel{^}{\mu}\left(1,1\right)\right)$ of the parameters $\mu $ in the normal environment (the stress level is taken ${X}_{0}=1$, while the degradation critical value is taken ${Y}_{0}=1$).

Maximum likelihood estimation of model parameters:

According to the experimental description and assumption, at a sample observation value, its log likelihood function is

${L}_{i,j,k}={I}_{i,j,k}\left(-\frac{1}{2}\mathrm{ln}\left(2\pi \right)-\mathrm{ln}\sigma -\frac{1}{2}{z}_{i,j,k}^{2}\right)+\left(1-{I}_{i,j,k}\right)\mathrm{ln}\left(1-\Phi \left({\eta}_{i,j}\right)\right)$

where ${z}_{i,j,k}={t}_{i,j,k}-\alpha -\beta {X}_{i}-\gamma {Y}_{j}$, ${\eta}_{i,j}={\tau}_{i}-\alpha -\beta {X}_{i}-\gamma {Y}_{j}$.

When ${t}_{i,j,k}\le {\tau}_{i}$, ${I}_{i,j,k}=1$ ; otherwise ${I}_{i,j,k}=0$. $\Phi (\cdot )$ is the distribution function of standard normal distribution. According to the description of the test, each sample has m degradation critical value to judge whether it is invalid. Therefore, $m\cdot n$ observations will be obtained for n samples. From their independence, the log likelihood function of all observations is:

$L={\displaystyle \underset{i=1}{\overset{q}{\sum}}{\displaystyle \underset{j=1}{\overset{m}{\sum}}{\displaystyle \underset{k=1}{\overset{{n}_{i}}{\sum}}{L}_{i,j,k}}}}$

The maximum likelihood estimation $\stackrel{^}{\alpha},\stackrel{^}{\beta},\stackrel{^}{\gamma},\stackrel{^}{\sigma}$ of parameters $\alpha ,\beta ,\gamma ,\sigma $ can be obtained by maximizing L. At this time, the analytical expression cannot be obtained. The numerical results can be obtained by using the nonlinear programming method in reference [15].

Taking $i=j=0$ in the Equation (2), and let ${X}_{0}=1$, ${Y}_{0}=1$, we obtain

$\mu \left(1,1\right)=\mu \left({X}_{0},{Y}_{0}\right)=\alpha +\beta {X}_{0}+\gamma {Y}_{0}=\alpha +\beta +\gamma $

Thus, the estimate of $\mu $ is

$\stackrel{^}{\mu}\left(1,1\right)=\stackrel{^}{\alpha}+\stackrel{^}{\beta}+\stackrel{^}{\gamma}$ (3)

In order to obtain the progressive variance of $\stackrel{^}{\mu}\left(1,1\right)$, we first find the Fisher information matrix of parameter maximum likelihood estimation.

According to the above form of log likelihood function ${L}_{i,j,k}$, We can get the Fisher information matrix of an observation value as follows

$\begin{array}{c}{F}_{i,j}=\left[\begin{array}{cccc}E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{{\partial}^{2}\alpha}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \alpha \partial \beta}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \alpha \partial \gamma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \alpha \partial \sigma}\right)\\ E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \alpha \partial \beta}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{{\partial}^{2}\beta}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \beta \partial \gamma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \beta \partial \sigma}\right)\\ E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \alpha \partial \gamma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \beta \partial \gamma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{{\partial}^{2}\gamma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \gamma \partial \sigma}\right)\\ E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \alpha \partial \sigma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \beta \partial \sigma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{\partial \gamma \partial \sigma}\right)& E\left(-\frac{{\partial}^{2}{L}_{i,j,k}}{{\partial}^{2}\sigma}\right)\end{array}\right]\\ =\frac{1}{{\sigma}^{2}}\left[\begin{array}{cccc}{A}_{i,j}& {X}_{i}{A}_{i,j}& {Y}_{j}{A}_{i,j}& {B}_{i,j}\\ {X}_{i}{A}_{i,j}& {X}_{i}^{2}{A}_{i,j}& {X}_{i}{Y}_{j}{A}_{i,j}& {X}_{i}{B}_{i,j}\\ {Y}_{j}{A}_{i,j}& {X}_{i}{Y}_{j}{A}_{i,j}& {Y}_{j}^{2}{A}_{i,j}& {Y}_{j}{B}_{i,j}\\ {B}_{i,j}& {X}_{i}{B}_{i,j}& {Y}_{j}{B}_{i,j}& {C}_{i,j}\end{array}\right]\end{array}$

where,

${A}_{i,j}=\Phi \left({\eta}_{i,j}\right)-\varphi \left({\eta}_{i,j}\right)\left({\eta}_{i,j}-\frac{\varphi \left({\eta}_{i,j}\right)}{1-\Phi \left({\eta}_{i,j}\right)}\right)$,

${B}_{i,j}=-\varphi \left({\eta}_{i,j}\right)\left\{1+{\eta}_{i,j}\left({\eta}_{i,j}-\frac{\varphi \left({\eta}_{i,j}\right)}{1-\Phi \left({\eta}_{i,j}\right)}\right)\right\}$

${C}_{i,j}=2\Phi \left({\eta}_{i,j}\right)-{\eta}_{i,j}\varphi \left({\eta}_{i,j}\right)\left(1+{\eta}_{i,j}^{2}-\frac{{\eta}_{i,j}\varphi \left({\eta}_{i,j}\right)}{1-\Phi \left({\eta}_{i,j}\right)}\right)$,

$\varphi (\cdot )$ is the distribution function of standard normal distribution.

Thus, the Fisher information matrix of all observations is $F={\displaystyle \underset{i=}{\overset{q}{\sum}}{\displaystyle \underset{j=1}{\overset{m}{\sum}}n{\pi}_{i}{F}_{i,j}}}$.

So, the covariance matrix of the MLE $\stackrel{^}{\alpha},\stackrel{^}{\beta},\stackrel{^}{\gamma},\stackrel{^}{\sigma}$ of the parameter $\alpha ,\beta ,\gamma ,\sigma $ can be calculated as

$\sum =\left[\begin{array}{cccc}Var\left(\stackrel{^}{\alpha}\right)& Cov\left(\stackrel{^}{\alpha},\stackrel{^}{\beta}\right)& Cov\left(\stackrel{^}{\alpha},\stackrel{^}{\gamma}\right)& Cov\left(\stackrel{^}{\alpha},\stackrel{^}{\sigma}\right)\\ Cov\left(\stackrel{^}{\alpha},\stackrel{^}{\beta}\right)& Var\left(\stackrel{^}{\beta}\right)& Cov\left(\stackrel{^}{\beta},\stackrel{^}{\gamma}\right)& Cov\left(\stackrel{^}{\beta},\stackrel{^}{\sigma}\right)\\ Cov\left(\stackrel{^}{\alpha},\stackrel{^}{\gamma}\right)& Cov\left(\stackrel{^}{\beta},\stackrel{^}{\gamma}\right)& Var\left(\stackrel{^}{\gamma}\right)& Cov\left(\stackrel{^}{\gamma},\sigma \right)\\ Cov\left(\stackrel{^}{\alpha},\stackrel{^}{\sigma}\right)& Cov\left(\stackrel{^}{\beta},\stackrel{^}{\sigma}\right)& Cov\left(\stackrel{^}{\gamma},\stackrel{^}{\sigma}\right)& Var\left(\stackrel{^}{\sigma}\right)\end{array}\right]={F}^{-1}$

According to the theory of maximum likelihood estimation and (3), we have

${h}^{\text{T}}=\left(\frac{\partial \stackrel{^}{\mu}\left(1,1\right)}{\partial \alpha},\frac{\partial \stackrel{^}{\mu}\left(1,1\right)}{\partial \beta},\frac{\partial \stackrel{^}{\mu}\left(1,1\right)}{\partial \gamma},\frac{\partial \stackrel{^}{\mu}\left(1,1\right)}{\partial \sigma}\right)=\left(1,1,1,0\right)$

Thus, we obtain

$AVar\left(\stackrel{^}{\mu}\left(1,1\right)\right)={h}^{\text{T}}\cdot \Sigma \cdot h=\left(1,1,1,0\right)\cdot \Sigma \cdot {\left(1,1,1,0\right)}^{\text{T}}$ (4)

In general, before the test, the maximum stress that the sample can withstand and the extreme critical value of degradation have been determined. Therefore, when the sample size n is determined, we can select the value of ${X}_{i}$, ${Y}_{j}$, ${\pi}_{i}$ to minimize the above-mentioned progressive variance $AVar\left(\stackrel{^}{\mu}\left(1,1\right)\right)=\left(1,1,1,0\right)\cdot \Sigma \cdot {\left(1,1,1,0\right)}^{\text{T}}$. So the optimal model can be described as follows

$\text{Min}:AVar\left(\stackrel{^}{\mu}\left(1,1\right)\right)=\left(1,1,1,0\right)\cdot \Sigma \cdot {\left(1,1,1,0\right)}^{\text{T}}$ (5)

where, the decision variables are: ${X}_{i}$, ${Y}_{j}$, ${\pi}_{i}$, $i=1,2,\cdots ,q-1$ ; $j=1,2,\cdots ,m-1$, and $0\le {X}_{i},{Y}_{j},{\pi}_{i}\le 1$.

3.3. A Special Case of Optimization Model

As a special case of optimization model, we set $q=2,m=2$ to calculate $AVar\left(\stackrel{^}{\mu}\left(1,1\right)\right)$.

In this case, the Fisher information matrix is $F={\displaystyle \underset{i=}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}n{\pi}_{i}{F}_{i,j}}}=\frac{n}{{\sigma}^{2}}{F}_{0}$, ${F}_{0}={\left({f}_{i,j}\right)}_{4\times 4}$ is a fourth-order matrix, where:

${f}_{1,1}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{\pi}_{i}{A}_{i,j}}}$, ${f}_{2,2}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{X}_{i}^{2}{\pi}_{i}{A}_{i,j}}}$, ${f}_{3,3}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{Y}_{j}^{2}}}{\pi}_{i}{A}_{i,j}$,

${f}_{4,4}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{\pi}_{i}{C}_{i,j}}}$, ${f}_{1,2}={f}_{2,1}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{X}_{i}{\pi}_{i}{A}_{i,j}}}$, ${f}_{1,3}={f}_{3,1}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{Y}_{j}{\pi}_{i}{A}_{i,j}}}$,

${f}_{1,4}={f}_{4,1}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{\pi}_{i}{B}_{i,j}}}$, ${f}_{2,3}={f}_{3,2}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{X}_{i}{Y}_{j}{\pi}_{i}{A}_{i,j}}}$,

${f}_{2,4}={f}_{4,2}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{X}_{i}{\pi}_{i}{B}_{i,j}}}$, ${f}_{3,4}={f}_{4,3}={\displaystyle \underset{i=1}{\overset{2}{\sum}}{\displaystyle \underset{j=1}{\overset{2}{\sum}}{Y}_{j}{\pi}_{i}{B}_{i,j}}}$,

Substituting the above results into (5), we obtain

$AVar\left(\stackrel{^}{\mu}\left(1,1\right)\right)=\left(1,1,1,0\right)\cdot \Sigma \cdot {\left(1,1,1,0\right)}^{\text{T}}=\left(1,1,1,0\right)\cdot {F}_{0}^{-1}\cdot {\left(1,1,1,0\right)}^{\text{T}}=\frac{n}{{\sigma}^{2}}V$

where, $V=\frac{{E}_{1}+{E}_{2}+{E}_{3}+2\left({E}_{4}+{E}_{5}+{E}_{6}\right)}{\left|{F}_{0}\right|}$, $\left|{F}_{0}\right|$ is the determinant of matrix ${F}_{0}$.

${E}_{1}={f}_{2,2}{f}_{3,3}{f}_{4,4}+2{f}_{2,3}{f}_{3,4}{f}_{4,2}-{f}_{2,4}^{2}{f}_{3,3}-{f}_{3,4}^{2}{f}_{2,2}-{f}_{2,3}^{2}{f}_{4,4}$

${E}_{2}={f}_{1,4}{f}_{2,4}{f}_{3,3}+{f}_{1,2}{f}_{3,4}^{2}+{f}_{1,3}{f}_{2,3}^{2}{f}_{4,4}-{f}_{1,2}{f}_{3,3}{f}_{4,4}-{f}_{1,4}{f}_{2,3}{f}_{3,4}-{f}_{1,3}{f}_{2,4}{f}_{3,4}$

${E}_{3}={f}_{1,2}{f}_{2,3}{f}_{4,4}+{f}_{1,4}{f}_{2,2}{f}_{3,4}+{f}_{1,3}{f}_{2,4}^{2}-{f}_{1,4}{f}_{2,3}{f}_{2,4}-{f}_{1,2}{f}_{2,4}{f}_{3,4}-{f}_{1,3}{f}_{2,2}{f}_{4,4}$

${E}_{4}={f}_{1,1}{f}_{3,3}{f}_{4,4}+2{f}_{1,3}{f}_{1,4}{f}_{3,4}-{f}_{1,4}^{2}{f}_{3,3}-{f}_{1,1}{f}_{3,4}^{2}-{f}_{1,3}^{2}{f}_{4,4}$

${E}_{5}={f}_{1,4}^{2}{f}_{2,3}+{f}_{1,1}{f}_{2,4}{f}_{3,4}+{f}_{1,2}{f}_{1,3}{f}_{4,4}-{f}_{1,1}{f}_{2,3}{f}_{4,4}-{f}_{1,2}{f}_{1,4}{f}_{3,4}-{f}_{1,4}{f}_{1,3}{f}_{2,4}$

${E}_{6}={f}_{1,1}{f}_{2,2}{f}_{4,4}+2{f}_{1,2}{f}_{1,4}{f}_{2,4}-{f}_{1,4}^{2}{f}_{2,2}-{f}_{1,1}{f}_{2,4}^{2}-{f}_{1,2}^{2}{f}_{4,4}$

Because the total number n of samples and ${\sigma}^{2}$ are constant, the optimization model is:

$\text{Min}:V=\frac{{E}_{1}+{E}_{2}+{E}_{3}+2\left({E}_{4}+{E}_{5}+{E}_{6}\right)}{\left|{F}_{0}\right|}$ (6)

where, the decision variables are: ${X}_{1}$, ${Y}_{1}$, ${\pi}_{1}$, and $0\le {X}_{1},{Y}_{1},{\pi}_{1}\le 1$.

This nonlinear programming problem can be solved by the method of reference [15].

4. Numerical Example

In order to illustrate the above optimization process, we take 50 LEDs for accelerated degradation test. The degraded performance attribute is taken as its luminous intensity. The acceleration variable of the test is taken as voltage. Under normal working conditions, the working current of the sample is 40 mA. According to previous experience, the maximum allowable working current is 330 mA (without changing its failure mechanism). Here we consider accelerated degradation tests at two stress levels. The voltage level is ${S}_{1},{S}_{2}$. $40<{S}_{1}<330$, ${S}_{2}=330$. When the sample is working under normal conditions, if its luminous intensity is lower than 50% of its initial intensity, it will be judged as invalid. If the initial luminous intensity is set to 1, the critical value under normal conditions is ${C}_{0}=0.5$. Now, we accelerate the degradation threshold at the same time. Set the extreme critical value as ${C}_{2}=0.1$. We optimize the experimental arrangement according to the above model.

The estimation of parameters in the model can be based on previous experience or similar data. Here, we set its estimate as $\stackrel{^}{\alpha}=5.5$, $\stackrel{^}{\beta}=6$, $\stackrel{^}{\gamma}=3$, $\stackrel{^}{\sigma}=1$. We take the cut-off time of the test at two stress levels as ${\tau}_{1}={\tau}_{2}=8$. Put the above data into the optimization model (6), using the nonlinear programming method, we obtain ${X}_{1}=0.412$, ${Y}_{1}=0.646$, ${\pi}_{1}=0.557$

Thus, we have ${n}_{1}=n\cdot {\pi}_{1}=50\times 0.557\approx 28$, ${n}_{2}=n-{n}_{1}=22$.

And ${X}_{1}=\frac{{S}_{1}-{S}_{2}}{{S}_{0}-{S}_{2}}$, ${S}_{1}=210.5$ ; and ${Y}_{1}=\frac{{C}_{1}-{C}_{2}}{{C}_{0}-{C}_{2}}$, ${C}_{1}=0.24$.

Therefore, the optimal design of the test is: Take 28 samples and put them into 210.5 mA for accelerated degradation test. When the luminous intensity of the test sample—led is lower than 24% of the initial value, we determine that the sample is invalid once. When the luminous intensity of the test sample is lower than 50% of the initial value, we decide that it will fail again. The closing time of the test is $t=8$. In addition, 22 samples are placed at 330 mA for accelerated degradation test. When the luminous intensity of the test sample is lower than 24% of the initial value, we determine that the sample is invalid once. When the luminous intensity of the test sample is lower than 50% of the initial value, we decide that it will fail again. The closing time of the test is $t=8$.

5. Results and Conclusion

In this work we established an optimization design model based on degradation performance for accelerated degradation test. The model shows that the test can be accelerated by selecting the critical value level. For ease of application, we simplified the model, and presented a numerical example to illustrate the procedures of the test plan.

Acknowledgements

This research was partially supported by the PhD research startup foundation of Guizhou Normal University (Grant No. GZNUD[2017]27) & Science and Technology Foundation of Guizhou Province (LKS[2012]11), China. & Teaching Project of Guizhou Normal University in 2016: Contract No. [2016] XJ No. 09.

References

[1] Boulanger, M. and Escobar, L.A. (1994) Experimental Design for a Class of Accelerated Degradation Tests . Technometrics, 36, 260-272.

https://doi.org/10.1080/00401706.1994.10485803

[2] Yu, H.F. and Tseng, S.T. (1999) Designing a Degradation Experiment. Naval Research Logistics, 46, 689-706.

https://doi.org/10.1002/(SICI)1520-6750(199909)46:6<689::AID-NAV6>3.0.CO;2-N

[3] Park, J.I. and Yum, B.J. (1997) Optimal Design of Accelerated Degradation Tests for Estimating Mean Lifetime at the Use Condition. Engineering Optimization, 28, 199-230.

https://doi.org/10.1080/03052159708941132

[4] Tseng, S.T. and Yu, H.F. (1997) A Termination Rule for Degradation Experiments. IEEE Transactions on Reliability, 46, 130-133.

https://doi.org/10.1109/24.589938

[5] Yang, G. and Yang, K. (2002) Accelerated Degradation-Tests with Tightened Critical Values. IEEE Transactions on Reliability, 51, 463-468.

https://doi.org/10.1109/TR.2002.804490

[6] Wu, S.J. and Chang, C.T. (2002) Optimal Design of Degradation Tests in Presence of Cost Constraint. Reliability Engineering & System Safety, 76, 109-115.

https://doi.org/10.1016/S0951-8320(01)00123-5

[7] Yu, H.F. and Tseng, S.T. (2004) Designing a Degradation Experiment with a Reciprocal Weibull Degradation Rate. Quality Technology & Quantitative Management, 1, 47-63.

https://doi.org/10.1080/16843703.2004.11673064

[8] Park, S.J., Yum, B. J. and Balamurali, S. (2004) Optimal Design of Step-Stress Degradation Tests in the Case of Destructive Measurement. Quality Technology & Quantitative Management, 1, 105-124.

https://doi.org/10.1080/16843703.2004.11673067

[9] Wang, X. and Xu, D. (2010) An Inverse Gaussian Process Model for Degradation Data. Technometrics, 52, 188-197.

https://doi.org/10.1198/TECH.2009.08197

[10] Ye, Z.S. and Chen, N. (2014) The Inverse Gaussian Process as a Degradation Model. Technometrics, 56, 302-311.

https://doi.org/10.1080/00401706.2013.830074

[11] Sung, S.I. and Yum, B.J. (2016) Optimal Design of Step-Stress Accelerated Degradation Tests Based on the Wiener Degradation Process. Quality Technology & Quantitative Management, 13, 367-393.

https://doi.org/10.1080/16843703.2016.1189179

[12] Yang, K. and Yang, G. (1998) Degradation Reliability Assessment Using Severe Critical Values. International Journal of Reliability, Quality and Safety Engineering, 5, 85-95.

https://doi.org/10.1142/S0218539398000091

[13] Lu, J.C. and Yang, P.Q. (1997) Statistical Inference of a Time-to-Failure Distribution Derived from Linear Degradation Data. Technometrics, 39, 391-400.

https://doi.org/10.1080/00401706.1997.10485158

[14] Lall, P., Pecht, M. and Hakim, E. (1997) Influence of Temperature on Microelectronics and System Reliability. CRC Press, London.

[15] Bazaraa, M. and Shetty, C. (1979) Nonlinear Programming: Theory and Algorithms. John Wiley & Sons, New York.