Optimization of the Accelerated Degradation Test Plan for Electrical Connector Contact Pairs Based on a Nonlinear Wiener Process

Accelerated degradation test is an effective method to evaluate the reliability of products with long life and high reliability. (e performance of most products fluctuates randomly in the degradation process, so it is suitable to use Wiener process. At present, the diffusion coefficient is regarded as constant in Wiener process, while the drift coefficient is related to stress. However, in practice, the amplitude of product performance fluctuation increases with the increase of stress level, which is not constant. (erefore, for the nonlinear Wiener case where both the drift coefficient and the diffusion coefficient are stress dependent, this paper studies the constant-stress accelerated degradation test theories and methods. Taking the contact pairs of electrical connectors as the research object, the minimum variance of reliable life estimate under normal stress is taken as the target. After determining the censored time at each stress level, the test stress level, the sample distribution ratio at each stress level, and the test interval at the one-third power scale of time are taken as design variables. (e test plan under 3, 4, and 5 stress levels is optimized and compared with the general test plan. (e influence of the difference between high and low stress levels on the evaluation accuracy is analyzed. Finally, the sensitivity analysis of parameters shows that the optimization plan has good robustness, and the change of stress quantity has little influence on the robustness of the plan.


Introduction
With the rapid development of modern science and technology, many products with high reliability and long life have emerged in the fields of machinery, electronics, aerospace, etc., which makes it more difficult for engineers to obtain enough failure data in a reasonable and effective time.
erefore, an accelerated degradation test (ADT) method is developed, which studies the product performance degradation process and extrapolates the reliability measure of products by statistical analysis of the degradation data and thus solves the reliability evaluation problem under lack of failure data or no failure data. At present, the research on accelerated degradation test focuses on performance degradation models, statistical analysis methods for degradation data, and optimization design for accelerated degradation test. e most widely used performance degradation models are Wiener processes and Gamma processes. Wiener process is suitable for nonmonotonic degradation of product performance and has been widely studied and applied in the analysis of degraded data [1]. Whitmore [2] established a degradation model that takes into account the randomness of performance degradation and measurement error. e randomness of performance degradation is described by Wiener process, and the randomness of measurement error is described by normal distribution. Tseng et al. [3,4] established an acceleration model of LEDs under step stress and multiple stresses based on Wiener process, respectively. Ye [5] considered the feature of high reliability of products with high degradation rate and high volatility and proposed a Wiener degradation model with linear relationship between diffusion coefficient and drift coefficient. On the basis of YE, Feng [6] proposed a degradation model, which is more consistent with product degradation data by considering the influence of measurement error. e Gamma process is used to describe the monotonic degradation of performance. Zhao [7] proposed a degradation modeling method based on Gamma process in the situation where the product degradation is caused by the continuously cumulative damages. Wang [8] studied parameter estimation method of nonhomogeneous Gamma process with random effects and proposed a semiparametric pseudolikelihood estimation method. e statistical analysis of degradation data usually includes the parameter estimation and reliability index estimation. ere are many estimation methods for unknown parameters of models, such as least squares (LSE), Bayesian (BAYES) estimation, maximum likelihood estimation (MLE), and expectation maximum (EM) algorithms. Liu [9] and Wang [10] solved unknown parameters of degradation path model with random effect by combining LSE and MLE. About the related works using maximum likelihood estimation methods, the reader can refer to [11][12][13]. Si [11], Ye [12], and Wang [13] estimated the parameters of degradation model by using MLE. Zhang [14] trained hyperparameters in the inverse problem models by using MLE. Robinson [15] and Gebraeel [16] estimated mixed effect model parameter by using the Bayesian estimation method. Li [17] proposed a grid-based EM method, which can learn online and classify signal patterns; the computational efficiency had greatly improved compared with traditional EM method. Chen [18] proposed a method that combines EM algorithm and unscented Kalman filter for estimating the accurate white noise process. e objective of the optimal design of accelerated degradation test plan is to obtain as much degradation information as possible under the shortest test time and the smallest test cost to improve the reliability evaluation accuracy. e optimal design of the accelerated degradation test plan mainly focuses on the situation that the product performance degradation model obeys the Wiener process and the Gamma process. When the performance degradation model obeys the linear Wiener process, the related work can be referred to [19,20] With the minimum cost as the optimization goal, Tang [19] proposed a step accelerated degradation test plan considering the evaluation accuracy. With the minimum asymptotic variance of the maximum likelihood estimation of the qth quantile of the lifetime distribution as goal, Heonsang et al. [20] optimized the design of constant accelerated degradation test. When the product performance degradation process is as nonlinear Wiener process, Liao and Tseng [4] provided the optimal step-stress accelerated degradation test plans by minimizing the asymptotic variance of the estimated 100pth percentile of the product's lifetime distribution under constraints on the total cost. When the performance degradation model obeys the Gamma process, some reference materials can be found in [21,22]. Tseng et al. [21] provided the optimal step-stress accelerated degradation test plans by minimizing the approximate variance of the estimated MTTF of the lifetime distribution of the product under constraints on the total cost. Tsai et al. [22] provided the optimal design plans for degradation tests based on a Gamma degradation process with random effects by minimizing the asymptotic variance of the estimate of the 100pth percentile of the lifetime distribution of the product under the constraint that the total experimental cost does not exceed a prespecified budget.
As a basic electromechanical component for transmitting signals and electric energy on the system, the various electrical connector is widely used in the system; meanwhile, its reliability directly determines the reliability of the system equipment; therefore, it is important to correctly evaluate the electrical connector reliability. With the progress of electrical connector process technology, its reliability is higher and has longer lifetime. As present, the accelerated degradation test method usually is used to evaluate electrical connector reliability index. In the test, the entire electrical connector is used as a test sample. e reliability index of the electrical connector is evaluated based on the performance degradation data obtained. If electrical connector model specification is changed, we need to carry out an accelerated degradation test of this model specification. is will consume a lot of manpower, material resources, and financial resources; therefore, the accelerated degradation test of electrical connector component is developed. No type of electrical connector will be omitted by the combination of component type spectrum. e literature [23] indicates that the lifetime of an electrical connector is primarily determined by the contact pair performance. erefore, if electrical connector contact pairs are considered as the accelerated test object, it not only solves reliability assessment problem of various types of electrical connector but also saves costs and improves efficiency. e contact resistance growth process of contact pairs of electrical connectors obeys the nonlinear Wiener process and has certain fluctuation. In addition, in most of the above accelerated test optimization design studies, only the drift coefficient is related to the stress in the nonlinear Wiener process model, and the diffusion coefficient is a constant.
e literature [23] shows that, for most electronic products, the fluctuation range of the performance parameters increases with the increase of the stress level; that is, the diffusion coefficient is not constant. erefore, the assumption of a constant diffusion coefficient will affect the evaluation accuracy of the test. As such, with the contact pair of a certain type of aerospace electrical connector as test object and the minimum variance of product reliable life maximum likelihood estimation value at normal stresses as goal, this paper studies the constant-stress accelerated degradation test theories and methods; meanwhile, this paper also optimizes the constant-stress accelerated degradation test plan of contact pairs of electrical connectors. e rest of this article is as follows. Section 2 presents the statistical model of reliability and the design criterion of accelerated degradation test plan. Section 3 describes the optimization theory of accelerated degradation test plan design. Section 4 establishes the mathematical model to optimize the accelerated degradation test plan design.
rough a case study of the contact pairs of certain type electrical connector, Section 5 illustrates the proposed optimal design method and compared the optimal design plan with traditional test plan. Section 6 analyzes the sensitivity of the parameters. Section 7 concludes this article.

Contact Performance Degradation Model for Contact Pairs of Electrical
Connectors. An aerospace electrical connector transmits electrical signals via contact spots on the contact pair surface. e contact performance gradually degrades due to temperature stress during long-term storage. Due to the fact that the growth of oxidation corrosive substances on the surface of contact spots obeys a nonlinear Wiener process, macroscopically, the degradation process of contact resistance obeys the nonlinear Wiener process [23]; namely, where r(t) represents the contact resistance at time t; r 0 represents the initial contact resistance, u is the drift coefficient, σ is the diffusion coefficient, W(τ(t)) denotes Brownian motion with mean 0 and variance τ(t), and, according to the law of oxidation corrosive substances growth obeying the cubic law and both diffusion term and drift term are related to the law of oxidation corrosive substances growth; therefore literature [24] indicates that τ(t) � t 1/3 (refer to Appendix B for the derivation process).
If the failure mechanism of the electrical connector remains unchanged at different temperatures stress, its specific degradation model can be expressed as follows: (1) e product degradation process is independent and obeys the nonlinear Wiener process (2) e relationship between the drift and diffusion coefficients in the Wiener process model and the temperature stress satisfies the following equations: where a, b, c 0 , and c 1 are unknown parameters, x � 1000/(273.15 + T) denotes the temperature stress level after conversion, φ(x) � (x − x 0 )/(x M − x 0 ), x M , x 0 , respectively, represent the highest and normal temperature stresses, and T is the ambient temperature (°C)

Lifetime Distribution Model for Contact Pairs of Electrical Connectors.
For degradation failure products, when the performance of a contact pair degrades with time and reaches the failure threshold D, the product fails. us, the time when the performance parameter is first taken to reach failure threshold D is defined as the lifetime of the product. When the degradation process follows the Wiener process, the product lifetime follows the inverse Gaussian distribution; that is, the lifetime distribution of the contact pairs of electrical connectors is

Guidelines for Test Design.
In order to facilitate the optimal design of constant-stress accelerated degradation test plan for contact pairs of electrical connectors, the following assumptions are given: (1) During the test, all samples should be tested simultaneously. e censored time for each test stress level is τ i (i � 1, · · ·, M, where M is the number of test stress levels) unless all samples fail. (2) In the degradation test, the highest stress x M and the normal stress x 0 are given in advance. e highest stress level is determined based on the criterion that the product failure mechanism does not change [25]. (3) To ensure that sufficient degradation data are obtained for each test stress level, the sample size for each stress level n i cannot be less than a certain minimum value n min . According to the literature [26], in order to ensure statistical accuracy, the sample size at each stress level should not be less than 5, so the minimum value of the sample size is taken as n min � 5. (4) In order to ensure sufficient degradation of contact performance in a certain period of time, the acceleration coefficient at the lowest test stress level cannot be smaller than a certain value.

Guidelines for Optimization of Degradation Test Plans.
In order to ensure the accuracy of the product lifetime estimation under certain test time and test cost, with the minimum variance of product pth quantile lifetime estimation value at normal stress x 0 as goal, the constant-stress accelerated degradation test plan is optimized.

Method of Lifetime Estimation.
e maximum likelihood estimation method (MLE) is a statistical method widely used to process accelerated degradation test data. Estimators obtained by the MLE have many excellent properties such as being asymptotically unbiased, asymptotically normal, and asymptotically efficient, and they are easy to be computed. For large sample sizes, the standard deviation of the MLE is the smallest compared to those of other estimation methods. is approach is also preferred to other estimation methods when the sample size is small. erefore, in the process of optimization of degradation test plan, MLE is used to calculate the pth quantile of the product lifetime distribution and its variance. In addition, prior to the degradation test plan optimization, the estimated values of the degradation model parameters a, b, c 0 , and c 1 must be known. ese can only be estimated from experience or preliminary test data. us, the test plan determined in this manner may not be optimal, but, undoubtedly, it is always better to determine the test plan based on these roughly estimated model parameters than blindly.

Standardization of Test Stress.
To facilitate the statistical analyses of the data, simplify the model, and make the conclusions more general, the test stress level x i is standardized as follows: us,

Sample Likelihood Function Based on the Nonlinear
Wiener Process. It is assumed that, at the ith test stress level ζ i (i � 1, 2, · · · , M), n i contact pairs are subjected to the accelerated degradation test. For the jth (j � 1, 2, · · · , n i ) contact pair, a contact resistance increment Δr ijk � r ijk − r ij(k−1) . k(k � 2, · · · , l) denotes the number of measurements.
As per the reliability statistical model, the contact resistance increment follows a normal distribution; that is, where us, the corresponding log-likelihood function can be obtained.
For a total of M stress levels, the log-likelihood function becomes

Information and Covariance
Matrix. According to the MLE theory, the variance and covariance matrix of the degradation model parameters a, b, c 0 , and c 1 are as follows: According to the MLE theory, the variance and covariance matrix Σ are the inverse matrix of the Fisher information matrix F; that is, Σ � F − 1 . e information matrix is obtained by taking the mathematical expectation of negative second partial derivative of the log-likelihood function L with respect to the degradation model parameters a, b, c 0 , and c 1 and is given by (refer to Appendix A for the derivation process), where

Asymptotic Variance of Function Estimation.
From equation (4), the pth quantile t p (ζ i ) of the lifetime distribution under stress ζ i is a function of parameters r 0 , D, a, b, c 0 , and c 1 . e literature [27] indicates that the pth quantile t p (ζ i ) under stress ζ i can be approximated as 4 Mathematical Problems in Engineering where τ −1 is inverse function of τ and z p is the pth quantile of standard normal distribution.
Since the maximum likelihood estimation has asymptotic normality, its asymptotic variance under a normal stress ζ 0 is where Let and then ,

Objective Function.
For an accelerated degradation test with the stress level number of M(0 < ζ 1 < ζ 2 < · · · < ζ M � 1) and a total of n contact pairs, nπ i samples are allocated to each stress level and π i is the ratio of samples input for each stress level. us, the Fisher information matrix F can be expressed as From equations (12) and (13), the asymptotic variance of the pth quantile of the MLE lifetime estimation for an electrical connector contact pairs at normal stress ζ 0 is where V M is the variance factor when the number of stress levels is M.
As per the provisions of the test method, before the plan optimization, the highest temperature stress level T M , normal temperature stress level T 0 , test censored time τ i at each stress level, and failure threshold D are given. Moreover, if the sample size n is also given, the model parameters r 0 , a, b, c 0 , and c 1 can also be obtained through preliminary tests. Consequently, under the principle that the variance of the pth quantile of the lifetime distribution at use stress level is minimum, when the number of stress levels is M, equation (17) shows that the variance factor V M is the objective function of the optimal design.

Selection of Design Variables.
Under the conditions that the model parameter estimation values r 0 , a, b, c 0 , and c 1 are known, if the test censored time τ i at each stress level, the highest stress ζ M , normal stress level ξ 0 , and the objective function are completely determined by the remaining M − 1 stress levels ξ 1 , ξ 2 , · · · , ξ M−1 and the allocation ratio of samples π 1 , π 2 , · · · , π M and test interval in the one-third power scale of time f 1 , f 2 , · · · , f M at each stress level, therefore the design variable for the optimal design is ζ 1 , ζ 2 , · · · , ζ M−1 , π 1 , π 2 , · · · , π M , where the test censored time τ i refers to the time when the test stops at each stress level, while the test interval f i (i � 1, · · · , M) refers to the interval of periodic inspection.

Determination of Constraints
(1) According to the fourth criterion of the test plan design, the standard test stress level should satisfy the following: (2) To ensure the acceleration of test and obtain more information about product degradation in the effective time, and according to the fourth criterion of the test plan design and the preliminary results, the acceleration factor τ m/0 should satisfy the following: where the acceleration factor τ m/0 is defined as the ratio of the pth quantile of the lifetime under the lowest test stress level ζ m relative to its counterpart at use stress levelζ 0 . (3) e allocation ratio of samples π 1 , π 2 , · · · , π M at each stress level must satisfy the following condition: M i�1 π i � 1 and 0 < π i < 1, (i � 1, 2, · · · , M − 1). (22) (4) According to the third criterion of the test plan design, the total sample size input for the test is n, and the sample size at each stress level should not be smaller than the minimum value n min . erefore, the sample size under each stress is the following: (5) On the one hand, in order to ensure that the sample has a certain degradation trend, the test time in the test chamber should not be less than 30 minutes; that is to say, the test interval of each stress on the scale of one-third of the time should not be less than 0.79 h; on the other hand, in order to obtain more reliability information of the product, the test time shall be no less than 6 times as much as possible.
According to the previous test results [28], the electrical connector did not fail in 1209 hours at 105°C. In order to obtain enough sample performance data before the test censored time, the censored time under each stress level of this test is set as 5000 hours. erefore, it is concluded that the test interval of each stress on the scale of one-third of time should not be higher than 2.85 h. In conclusion, the range of test interval f i at one-third time scale under each stress level is where f i is the test interval under each stress level and M is the number of stress levels. (6) In order to ensure that the contact pair of electrical connector has enough degradation in each test, the following conditions should be met in each test under each stress:

Determination of Sample
Size. Based on the MLE theory, under normal stress level, the estimation value t p (ζ 0 ) of the pth quantile at a certain time asymptotically follows a normal distribution with mean t p (ζ 0 ) and variance Var[t p (ζ 0 )]; that is, If t p (ζ 0 ) is replaced by t p (ζ 0 ), the upper and lower confidence limits, respectively, for t p (ζ 0 ) at the confidence level ϕ are where K ϕ is the (1 + ϕ)/2-th quantile of the standard normal distribution at the confidence level ϕ. Under the condition that the probability is ϕ, when t p (ζ 0 ) is required to be within the range of t p (ζ 0 ) ± W (where t p (ζ 0 ) denotes the true value), the sample size is the following:

Failure reshold and Test Stress Level.
In this paper, the contact pair of a certain type of aerospace electrical connector are considered as the test object. Test stresses must satisfy the principle that the failure mechanism at test stress levels is consistent with that at normal stress level. As per the fact that the failure of aerospace electrical connectors is contact failure (contact failure threshold, D � 5 mΩ) at a test temperature stress level of 158°C [28], we thereby conclude that the highest temperature stress level is T M � 158°C, while the corresponding stress after conversion is x M � 2.3194, and the standard stress is ζ M � 1. e results of the investigation on the actual storage profile of electrical connectors indicate that the normal temperature is T 0 � 30°C, while the converted stress is x 0 � 3.2987 and the standard stress is ξ 0 � 0.

Initial Estimation of Degradation Model Parameters.
Based on the results of the statistical analysis of the preliminary test data [24], the rough estimated values of the parameters of the reliability statistical model for the contact pairs of certain type electrical connector are obtained: Let p � 0.05; that is, reliability R � 0.95. e variance factor V M of the product lifetime maximum likelihood estimation value under a normal stress level is taken as the objective function. e lowest stress level ξ 1 and the allocation ratio of samples π 1 , π 2 , · · · , π M and test interval at the one-third power scale of time f 1 , f 2 , · · · , f M are used as the design variables; meanwhile, equal reciprocal interval of temperature method is used to calculate stress level. erefore, the intermediate stress level can be calculated according to formula (31). e accelerated degradation test plan for the contact pair of certain type of electrical connector is optimized.
e optimized test plans for stress levels of M � 3, 4, 5 are shown in Tables 1-3: If the confidence level and confidence width are defined as ϕ � 80% and 2W � 0.4, respectively, then K ϕ � 1.28155. erefore, the total numbers of samples for optimum test plans with stress levels of M � 3, 4, 5 are as follows: n M�3 � 37, n M�4 � 42, and n M�5 � 48.
From Tables 1-3, it can be found that, in the optimum test plan of stress level M � 3, 4, 5, the estimation accuracy decreases gradually with the increase of the number of stress levels, and, in each optimum test plan, the allocation ratio of sample distribution under the minimum stress level is larger than that under other stress levels. It shows that the minimum stress has greater influence on the estimation accuracy of the objective function than other stress levels.
Although the estimation accuracy of the three stress levels is the best, from the perspective of the robustness and effectiveness of the test, the engineering design can choose an optimal design with a stress level of M � 4.

Comparison of Optimal Test Plan and Traditional Test
Plan.
e optimal test plan for stress level M � 3, 4, 5 is shown in Tables 1-3

Analysis of the Influence of the Lowest Stress Level and the Difference of High and Low Stress Levels on Variance Factor
e allocation ratio of sample for the lowest stress level is higher than others in Tables 1-3; this result indicates that influence of the lowest stress level on variance factor V M is higher than other stress levels. erefore, this section analyzes how the lowest stress level T 1 influences variance factor V M . Taking the traditional test plan of four stress levels as example, this paper chooses the lowest stress level T 1 as 100°C, 105°C, 110°C, and 115°C; meanwhile, due to the preliminary test result, the highest stress level is 158°C in the above. In order to explore how the variance factor changes under different values of the highest stress level, the highest stress level T 2 is 125°C, 135°C, 145°C as a comparison; according to the equal interval formula, the middle stress level  Table 7. e results of Table 7 are as follows: Figure 1 shows that the variance factor V M increases with the increase of the lowest stress level T 1 when T 2 is constant and T 1 is increased from 100°C to 115°C (iii) Variance factor V M decreases with the increase of the highest stress level T 2 when T 1 is constant (iv) Figures 2 and 3 indicate that the variance factor V M increases with the decrease of the value of T 2 − T 1 ; therefore, the engineering test should choose the large difference value of T 2 and T 1 in reasonable range

Analysis of Parameter Sensitivity
Sensitivity is expressed as a dimensionless index; it reflects the extent to which the output value of the model changes with small changes in parameters. In view of several key parameters (a, b, c 0 , and c 1 ) of the reliability statistical model of electrical connector contact pairs, this paper explores how the change of parameters influences the object function value. e values of the other three parameters are fixed; one of the parameter values needs to be changed. e change range is taken as ±5%, ±10%, and ±15%.

Analysis of Parameter Sensitivity with Four Stress Levels.
When the number of stress levels is 4, the fluctuation results of object function value change with parameters are listed in Table 8 and Figures 4-6. τ 1 , τ 2 , τ 3 , and τ 4 represent the errors of parameters a, b, c 0 , and c 1 , respectively. e results are as follows.
According to Table 8 and Figures 4-6, when the error range of the parameter is ±5%, ±10%, and ±15%, the change  Mathematical Problems in Engineering amount of variance factor is very small, so it shows that the robustness of this optimization plan is good.

Analysis of Parameter Sensitivity with ree Stress Levels.
In order to explore whether the optimization plan still has a good robustness when the number of stresses changes, the optimization plan of 3 stress levels is explored when the error is ±10%. e results of Table 9 and Figure 7 show that the optimization plan of 3 stress levels is the same as the optimization plan of 4 stress levels, which has a good robustness. is shows that the robustness of the plan will not change greatly with the change of stress quantity.           electrical connectors subjected to nonlinear wiener process as an example, the optimization result of test plan with stress levels of 3, 4, and 5 shows that, at the same confidence level, the estimation accuracy of optimal test plan is improved compared to the traditional test plan; meanwhile, the total    number of samples is also reduced. e estimation accuracy with stress level of 3 is improved by a factor of 1.16 and the total sample size is reduced by a factor of 1.14. e estimation accuracy with stress level of 4 is improved by a factor of 1.21 and the total sample size is reduced by a factor of 1.2. e estimation accuracy with stress level of 5 is improved by a factor of 1.16 and the total sample size is reduced by a factor of 1.17. At the same estimation accuracy, the optimal design test plan time is shortened. e optimal design test plan time with stress level of 3 is reduced by a factor of 1.26. e optimal design test plan time with stress level of 4 is reduced by a factor of 1.31. e optimal design test plan time with stress level of 5 is reduced by a factor of 1.29. e effect of the lowest stress level on the variance factor V M is analyzed. e results show that the variance factor V M increases with the increase of the lowest stress level; meanwhile, the variance factor V M increases with the decrease of the difference between the highest stress and the lowest stress. erefore, the engineering test should choose the large difference value of T 2 and T 1 in reasonable range. When the error range of parameter is 5%−15%, the change of the estimation accuracy of the optimization plan is very small, the plan has good robustness, and the change of the stress quantity has little effect on the plan's robustness. e second partial derivatives of L i with respect to each parameter are given by