Stochastic Modeling and Analysis of Multiple Nonlinear Accelerated Degradation Processes through Information Fusion

Accelerated degradation testing (ADT) is an efficient technique for evaluating the lifetime of a highly reliable product whose underlying failure process may be traced by the degradation of the product’s performance parameters with time. However, most research on ADT mainly focuses on a single performance parameter. In reality, the performance of a modern product is usually characterized by multiple parameters, and the degradation paths are usually nonlinear. To address such problems, this paper develops a new s-dependent nonlinear ADT model for products with multiple performance parameters using a general Wiener process and copulas. The general Wiener process models the nonlinear ADT data, and the dependency among different degradation measures is analyzed using the copula method. An engineering case study on a tuner’s ADT data is conducted to demonstrate the effectiveness of the proposed method. The results illustrate that the proposed method is quite effective in estimating the lifetime of a product with s-dependent performance parameters.


Introduction
Modern products are developed to have good quality and high reliability. For some safety-critical components and systems, they are even designed to last for an extremely long time to avoid the catastrophic consequences of potential failures. It is a big challenge to obtain sufficient amount of time-to-failure data by testing such products under the normal operating environments and sometimes even under harsher conditions [1]. Obviously, the traditional failure-time data analysis and testing methods are unsuitable for reliability analysis of such highly reliable products.
In engineering applications, many failure mechanisms can be traced to underlying degradation processes (e.g., cumulative wear, crack growth, corrosion, fatigue, material aging, etc.) [2]. As a result, the degradation analysis method has been introduced to handle reliability modeling problems based on the products' degradation information obtained from historical data or degradation tests [3]. Especially, accelerated degradation testing (ADT) has been used to collect the degradation data by exposing the test specimens to severer-than-normal conditions. In the last two decades, ADT has been intensively studied as an effective tool for reliability verification and lifetime evaluation of modern products. Successful applications of ADT include the reliability analyses of batteries [4], super luminescent diode (SLD) [5], smart electricity meter [6], metal oxide semiconductor field effect transistors (MOSFETs) [7], etc. Due to these advantages, the copulas method has attracted much attention in reliability engineering. Sari et al. [26] introduced a copula function to describe the correlation between two performance parameters, and combined it with a generalized linear model for bivariate constant-stress degradation data. Pan et al. [27,28] discussed the bivariate degradation modeling approaches based on Wiener processes and Copulas, respectively, under the normal operating stress and constant-stress accelerating conditions as well. Similarly, Peng et al. [18], Liu et al. [29] and Hao et al. [30,31] also adopted copulas to characterize the dependency between two performance parameters, and developed their own bivariate degradation models. Wang and Pham [19] applied time-varying copulas for handling the s-dependent relationship among competing degradation processes, and provided a numerical example with two degradation processes. Hong et al. [32] used copulas to investigate the influence of degradation dependency of multiple components on the optimal maintenance decisions. Xi et al. [33] developed a copula-based sampling method for data-driven prognostics. However, most of the previous work focuses on degradation processes under normal operating conditions, and scarce studies have been conducted on the dependency of multiple degradation processes under accelerated conditions. This paper is aimed at making an early attempt to model s-dependent multivariate ADT data using general Wiener process and copulas. First, the general Wiener process is used to model nonlinear univariate accelerated degradation processes. Then, the copula method is adopted to find the joint probability distribution of multiple degradation processes. Next, parameter estimation of the proposed model is performed using the Inference Functions for Margins (IFM) method, and the Akaike Information Criterion (AIC) is employed to compare the goodness-of-fit of candidate models. Finally, a real-world case study on the tuner's CSADT data is utilized to demonstrate the usefulness of the proposed method.
The remainder of this paper is organized as follows: Section 2 presents the univariate ADT model based on the general Wiener process and its parameter estimation method. Section 3 elaborates on the copula-based multivariate ADT modeling method, including introduction to copulas, multivariate dependent accelerated degradation model, and statistical inference. Section 4 provides the case study, and Section 5 concludes the paper.

General Univariate Accelerated Degradation Model
As mentioned earlier, ADT attempts to obtain a product's degradation information more quickly by exposing the test units to harsher-than-normal conditions. In many applications, two ADT methods have been widely used: constant-stress ADT (CSADT) and step-stress ADT (SSADT). For illustration purposes, this paper will focus on CSADT.
ADT modeling is concerned with accelerated degradation models and data analyses [37]. Generally, accelerated degradation model is composed of both stress-acceleration model and degradation model, where the stress-acceleration model describes the relationship between the levels of accelerating stress and the degradation rate, and the degradation model depicts the evolution of degradation process over time. In this section, a nonlinear univariate CSADT model based on the general Wiener process is presented.

Stress-Acceleration Model
To estimate the lifetime or degradation rate of a product at the use conditions, a stress-acceleration model should be obtained first. Usually, stress-acceleration model can be obtained from either some physical knowledge about the product or empirical observations. Some widely used models include Arrhenius model, Eyring model, and inverse power model [38]. After a proper transformation of stress S i , a stress-acceleration model for degradation rate dr(S i ) of the product under stress level S i can be expressed as [11]: where a and b are unknown constants, and ϕ(S i ) is a function of S i depending on the stress type, e.g., ϕ(S i ) = 1/S i for absolute temperature and ϕ(S i ) = lnS i for voltage.

An ADT Model Based on the General Wiener Process
The general Wiener process used for degradation modeling can be expressed as [16]: where X(t) is the degradation value at time t; X(0) is the initial degradation value (for simplicity, we assume X(0) = 0); B(.) is the standard Brownian motion; µ is the drift coefficient reflecting the rate of degradation; σ is the diffusion parameter describing the unit-to-unit variability and variability due to operational and environmental conditions; λ(t; θ) and τ(t; γ) are the non-decreasing time-scale transformation which can describe the nonlinear property of the degradation process, and θ and γ are the corresponding parameters. By further assuming that Λ pt; θq " t 0 λ pt; θq dt, we have: Clearly, the basic Wiener process model [1], as shown in Equation (4), is a special case of Equation (3) when Λ(t; θ) = τ(t; γ) = Λ(t): For clarity, Equation (3) is called model M 0 , while Equation (4) is called model M 1 in this paper.
In model M 0 , since E[X(t)] = µΛ(t; θ) and Var[X(t)] = σ 2 τ(t; γ), one can see that θ reflects the behavior of the expectation of degradation pattern while γ reflects the behavior of the variance of degradation. In model M 1 , the two parameters are assumed to be the same.
When using model M 0 for ADT analysis, the drift coefficient µ is modeled as a function of stresses, i.e., as a stress-acceleration model [1,5,10]. According to Equation (1), we have: The corresponding lifetime distribution can be obtained based on the time when the degradation path X(t) exceeds the failure threshold ω for the first time, i.e., First Passage Time (FPT). According to the property of linear Wiener process, the probability density function (PDF) of FPT follows an inverse Gaussian distribution [39]: Hence, for model M 1 considering the time-scale transformation, the PDF of FPT and the reliability function of the product can be easily obtained as: respectively. However, the analytical formula of PDF of FPT for model M 0 , where Λ(t; θ)‰τ(t; γ), cannot be directly obtained from the above equations. Let the time-scale transformation be s = τ(t; γ), so that t = τ´1(s; γ). We also define Λ(t; θ) =Λ (τ´1(s; γ), θ) = ρ(s; θ). As a result, Equation (3) becomes [16]: Xpsq " µρps; θq`σBpsq (9) We define: Under some mild assumptions, the PDF of stochastic process X(s) can be approximated by (according to Theorem 2 in Si et al. [15] p. 56): ρ ps; θq σs`κ ps; θq σ˙¨e xp˜´p ω´µρ ps; θqq 2 2σ 2 s¸( 11) Hence, the PDF of failure time when X(t) exceeds ω can be obtained by substituting s with τ(t; γ) as: Equation (12) can be further approximated as: and the corresponding Cumulative Distribution Function (CDF) of FPT is: Suppose that Λ(t; θ) = t θ and τ(t; γ) = t γ , then Equation (12) becomes: When θ = γ = 1, Equation (13) is an inverse Gaussian distribution as in Equation (6) since (13) is an inverse Gaussian distribution with time-scale transformation as in Equation (7). Thus, Equation (13) is the generalized PDF of FPT which can cover existing linear and time-scale transformation Wiener process model as its limiting cases.
Then, according to Equation (14), the reliability function that relies on S 0 can be expressed as:

Parameter Estimation of General ADT Model
In this section, we briefly provide a two-stage Maximum Likelihood Estimation (MLE) for unknown parameters in CSADT with a single performance parameter. The unknown parameters are in vector Θ = {θ, γ, σ, a, b}.
Suppose that a univariate CSADT with K stress levels has been conducted, S i , i = 1, . . . , K is indicated the ith accelerated stress level, and S 0 is normal operating condition. There are n i specimens under stress level S i . Then, X ijk is the kth degradation value of unit j at the stress level i and t ijk is the corresponding measurement time, i = 1, 2, . . . , K; j = 1, 2, . . . , n i ; k = 1, 2, . . . , m ij . Let X ij " pX ij1 , X ij2 , . . . , X ijm ij q 1 and t ij " pΛpt ij1 ; θq, Λpt ij2 ; θq, . . . , Λpt ijm ij ; θqq 1 . According to the properties of Wiener process, X ij follows a multivariate normal distribution: ( 17) where: Let µ " pµ 11 , . . . , µ 1n 1 , . . . , µ K1 , . . . , µ Kn K q. The likelihood function of CSADT data can be easily obtained, and the corresponding log-likelihood function is: Taking the first partial derivatives of Equation (19) with respect to µ ij and σ 2 and setting the resulting equations to zero, the estimatesμ ij andσ 2 relying on θ and γ can be obtained as: Substituting Equation (20) into Equation (19), the log-likelihood function is only a function of θ and γ: The estimatesθ andγ can be obtained using the two-dimensional search [40] that maximizes Equation (21). Then,μ ij andσ 2 can be easily computed from Equation (20).
Finally, given theμ ij s and the corresponding stress S i , the values ofâ andb can be easily computed from Equation (5) by the regression analysis.

Copulas
Identifying a multivariate probability distribution is difficult due to the complexity of marginal distributions and the curse of dimensionality. Copulas simplify this process by separating the learning of the marginal distributions from the learning of the multivariate dependence structure that links the marginal distributions to form a joint distribution [35], as shown in Figure 1.
The following theorem provides the necessary and sufficient condition for the existence of Copula, which elucidates the role that copulas play in describing the relationship between a multivariate distribution and the associated univariate marginal distributions. [41]): Let X = (X 1 , . . . , X N ) be a random vector with marginal distributions F 1 (x 1 ), . . . , F N (x N ), and let F be their joint distribution function. Then, there exists a copula function C such that:

Theorem 1. (Sklar's theorem
. . , N, are continuous, the copula C is unique. Conversely, if F 1 (x 1 ), . . . , F N (x N ) are univariate distributions, the function F defined by Equation (22) is the joint distribution function of margins F 1 (x 1 ), . . . , F N (x N ). According to the theorem, a multivariate Copula function can be defined as follows: [36]: An N-dimensional copula is a function C from I N = [0,1] N to I with the following properties: (1) Grounded: for every u = (u1, …, uN) in I N , C(u) = 0 if at least one coordinate of u is 0; (2) Uniform margins: if all coordinates of u are 1 except uk, then: where N(z) = #{k: zk = uk}.
The density of a copula function C is denoted by c, which may be achieved by taking partial derivative as: Furthermore, the joint density function corresponding to distribution function F(x1, …, xN) can be calculated by: where fn(·) , n = 1, …, N, is the PDF of marginal distribution Fn(·).
Copulas have several attractive properties [42]. First, copulas provide a convenient way to model the marginal distributions and the joint dependence structure separately. Second, they are invariant under increasing and continuous transformations. Third, they have no constraints on the univariate marginal distribution. Because of these advantages, copulas are chosen to emphasize and measure the dependence among multiple degradation processes in this paper. Table 1 lists several commonly used multivariate copula functions, where u1, …, uN denote univariate marginal distributions. The Akaike Information Criterion (AIC) will be employed to According to the theorem, a multivariate Copula function can be defined as follows: [36]: An N-dimensional copula is a function C from I N = [0,1] N to I with the following properties: Cpuq " Cp1, . . . , 1, u k , 1, . . . , 1q " u k (3) N-increasing: for each hyperrectangle B " where N(z) = #{k: z k = u k }.
The density of a copula function C is denoted by c, which may be achieved by taking partial derivative as: Furthermore, the joint density function corresponding to distribution function F(x 1 , . . . , x N ) can be calculated by: where f n (¨), n = 1, . . . , N, is the PDF of marginal distribution F n (¨).
Copulas have several attractive properties [42]. First, copulas provide a convenient way to model the marginal distributions and the joint dependence structure separately. Second, they are invariant under increasing and continuous transformations. Third, they have no constraints on the univariate marginal distribution. Because of these advantages, copulas are chosen to emphasize and measure the dependence among multiple degradation processes in this paper. Table 1 lists several commonly used multivariate copula functions, where u 1 , . . . , u N denote univariate marginal distributions. The Akaike Information Criterion (AIC) will be employed to implement goodness-of-fit test and quantitatively select the most suitable copula type from several candidate copulas. In particular, the AIC is defined as: where k is the number of estimated parameters in the model, and L is the maximum likelihood value for the model. Given a set of candidate copulas for the multivariate accelerated degradation processes, the preferred copula function is the one with the minimum AIC value. Table 1. Summary of some multivariate copulas.
1 Φ is the distribution function of a standard normal random variable, Φ ρ is the N-variate standard normal distribution with mean vector 0 and covariance matrix ρ; 2 t υ is a univariate t distribution with υ degrees of freedom, t ρ,υ is the multivariate Student's t distribution with a correlation matrix ρ with υ degrees of freedom.

Multivariate Dependent Accelerated Degradation Model
In this section, a multivariate accelerated degradation model is presented to estimate the system reliability under the normal operating conditions. First, the CSADT data of each performance characteristic is modeled by a general Wiener process (see Section 2.2), and general accelerated degradation models for multiple performance parameters are established respectively. Then, copulas are utilized to describe the dependence among multiple performance parameters and conduct multi-parameters degradation information fusion.
To analyze the multivariate CSADT data, the following assumptions are made:

A1.
The degradation measurements of all specimens are taken at the same time.

A2.
A specimen is considered to be failed if one of the features reaches its corresponding failure threshold for the first time.

A3.
The degradation processes of all performance parameters in CSADT can be depicted by the general Wiener process.

A4.
For a dependent system with multiple performance parameters, the dependency among the multiple performance parameters can be characterized by copulas.
Suppose that a system has N degradation measures. Let X p (t), p = 1, 2, . . . , N denote the degradation process of the pth performance parameter at normal operating condition S 0 . According to above assumptions, the system is failed if any one performance parameter reaches its corresponding failure threshold represented by D p . Let the failure time of the pth characteristic index be T p , then the system lifetime is T = min(T 1 , . . . , T N ). So, the system reliability under S 0 can be expressed as: Rptq " PpT ą tq " PpT 1 ą t, . . . , T N ą tq " PpX 1 ptq ă D 1 , . . . , X N ptq ă D N q Apparently, if the multiple degradation measures are independent, Equation (28) becomes: where R p (t), p = 1, 2, . . . , N is the reliability function associated with the pth performance parameter under S 0 as defined in Equation (16). However, the multiple performance parameters are usually dependent. In order to address the dependency among multiple degradation processes, a copula function can be applied, and the best model can be selected based on AIC. Essentially, the system reliability at time t under S 0 can be performed as [19,26,27,30]: Rptq " CpR 1 ptq, R 2 ptq, . . . , R N ptq; δq where δ is the parameter set of a family of copulas.

Statistical Inference
Consider a copula-based multivariate distribution for random vector X with PDF: ( 31) where F p (; Θ p ) and f p (; Θ p ), p = 1, . . . ,N are CDF and PDF of the pth marginal distribution with parameter set Θ p = {θ p , γ p , σ p , a p , b p }, c(;δ) is the density of copula function C(;δ) with parameter set δ.
Clearly, the parameter set Ω = {Θ 1 , . . . , Θ N , δ} can be estimated by maximizing Equation (32). However, as the dimension of variables increases, it is difficult to obtain the optimal solution.
Joe [35] proposed a computationally attractive alternative to MLE, called Inference Functions for Margins (IFM) method, to estimate the parameters in multivariate copula models. The idea is to decompose Equation (32) into two parts: the contribution, denoted by L C , from the dependence structure in data represented by copula C, and the contributions, denoted by L p , p = 1, . . . , N, from each marginal distribution. Then, the approach estimates the parameters of marginal distributions and the copula in two stages. In the first stage, the parameters of each marginal distribution are estimated from the corresponding L p . In the second stage, the copula's parameters are obtained by maximizing L C using the marginal distribution parameters estimated in the first stage.
In this paper, the IFM method is used to estimate the parameters Ω = {Θ 1 , . . . , Θ N , δ} of the copula-based multivariate CSADT model in the following two stages. For the CSADT data of each degradation measure, the parameter estimation method given in Section 2.3 can be used to estimate the parameter set Θ p = {θ p , γ p , σ p , a p , b p } of each marginal distribution by maximizing the log-likelihood function. Then, assuming that the failure threshold of the pth performance characteristic is D p , p = 1, 2, . . . , N, we can respectively get the reliability functions of all performance parameters under normal operating condition S 0 using Equation (16).

Case Study
In this section, a case study is used to demonstrate the validity and estimation performance of the proposed copula-based multivariate CSADT model and computational method.

Problem Description
The tuner is a kind of microwave electronical assembly, which is required to have a long lifetime and high reliability. Previous failure analysis shows that temperature is the major factor that causes tuner output degradation. To quickly evaluate the lifetime and reliability of the tuner, a temperature CSADT with four stress levels was conducted. Six performance parameters of each tuner were measured by a computerized measuring system, including the power gain GA, GB, and GC, the noise figure NA, NB and NC, where G denotes power gain, N means noise figure, A, B and C indicate different directions. Figure 2 shows a sample from one tuner.

Problem Description
The tuner is a kind of microwave electronical assembly, which is required to have a long lifetime and high reliability. Previous failure analysis shows that temperature is the major factor that causes tuner output degradation. To quickly evaluate the lifetime and reliability of the tuner, a temperature CSADT with four stress levels was conducted. Six performance parameters of each tuner were measured by a computerized measuring system, including the power gain GA, GB, and GC, the noise figure NA, NB and NC, where G denotes power gain, N means noise figure, A, B and C indicate different directions. Figure 2 shows a sample from one tuner.
For all performance parameters, the sampling interval was Δt = 5 h. The other settings of the test are shown in Table 2, where S1, S2, S3, and S4 are the accelerated stress levels, and the normal operating condition is S0 = 25 °C. The accelerated degradation data of all performance parameters obtained by CSADT are shown in Figures 3-6.   For all performance parameters, the sampling interval was ∆t = 5 h. The other settings of the test are shown in Table 2, where S 1 , S 2 , S 3 , and S 4 are the accelerated stress levels, and the normal operating condition is S 0 = 25˝C. The accelerated degradation data of all performance parameters obtained by CSADT are shown in Figures 3-6.        Figure 7, one can see that all pairs of degradation measures are significantly correlated. Therefore, in order to achieve more accurate and practical evaluation results of the tuner's CSADT, it is worth considering the dependence among the multiple degradation processes.

Univariate ADT Models with General Wiener Process
First of all, the original degradation data of GA, GB, and GC are transformed by: where x0 and xi are the initial and ith degradation measurements.

Univariate ADT Models with General Wiener Process
First of all, the original degradation data of GA, GB, and GC are transformed by:

Univariate ADT Models with General Wiener Process
First of all, the original degradation data of GA, GB, and GC are transformed by: where x 0 and x i are the initial and ith degradation measurements. Similarly, the original degradation data of NA, NB, and NC are also transformed by: Then, the proposed general Wiener process model, i.e., Equation (3), is used to fit the CSADT data of each performance parameter separately. For convenience and without loss of generality, the formulas Λ(t; θ) = t θ and τ(t; γ) = t γ are used in this paper according to the behavior of the degradation over time. It is important to mention that although we only focus on a combination of Λ(t; θ) = t θ and τ(t; γ) = t γ , other forms of Λ(t; θ) and τ(t; γ) can be adopted without extra theoretical difficulties, if applied [16]. The estimates of univariate ADT model parameters for each performance characteristic are obtained by the method provided in Section 2.3 and are listed in Table 3. From Table 3, it can be seen that the estimate of γ for each performance parameter is close to 1, which indicates that the proposed model could be simplified. Simplify the univariate ADT model with the assumption of γ = 1, i.e., τ(t; γ) = t. The estimates of the simplified model (SM) parameters for each degradation measure are listed in Table 4. To justify the use of the full model (FM) or the simplified model (assume γ = 1) for the univariate ADT data of each performance parameter, the likelihood ratio (LR) test [43] is implemented. The test statistic distribution has 1 degree of freedom, and the significance level is 5%. The result of the likelihood ratio test is shown in Table 5. Since the critical value is χ 2 1,0.05 " 3.84, we accept the assumption of γ = 1 except for GA. While, the p-value of GA is close to 0, which also indicates that there is strong evidence suggesting that the full model fits the GA data better than the simplified model. Hence, the proposed general Wiener process model with τ(t; γ) = t γ is used to fit the CSADT data of GA, and the simplified model (assume γ = 1) is used to fit the other performance parameters of tuner.
For the failure thresholds of degradation measures, it is assumed that the tuner fails when one of the power gains, i.e., GA, GB, and GC, drops to 60% of its initial value. Meanwhile, the tuner would also fail if one of the noise figures, i.e., NA, NB, and NC, increases to 100% of its initial value. According to Equation (16), the marginal reliability functions of six performance parameters at the normal stress can be calculated, as shown in Figure 8. From Equation (29), we also estimate the system reliability of the tuner at the normal operating condition without considering the dependence among multiple performance parameters, which is displayed in Figure 8. From Figure 8, one can see that the tuner's reliability is far less than the marginal reliabilities of individual degradation processes.

Multivariate Accelerated Degradation Model with Copulas
Based on the marginal reliability functions, the copula method is applied to fit the dependent multivariate ADT data. To select the best copula, the AIC is used when comparing candidate copulas. Five widely used multivariate copulas are selected to fit the joint distribution of multiple degradation processes. The results of parameter estimation and AIC values are summarized in Table 6, and the joint reliability curves of tuner with different copulas are shown in Figure 9.
Reliability Reliability Figure 8. Marginal reliability function for each performance parameter.

Multivariate Accelerated Degradation Model with Copulas
Based on the marginal reliability functions, the copula method is applied to fit the dependent multivariate ADT data. To select the best copula, the AIC is used when comparing candidate copulas. Five widely used multivariate copulas are selected to fit the joint distribution of multiple degradation processes. The results of parameter estimation and AIC values are summarized in Table 6, and the joint reliability curves of tuner with different copulas are shown in Figure 9.

Multivariate Accelerated Degradation Model with Copulas
Based on the marginal reliability functions, the copula method is applied to fit the dependent multivariate ADT data. To select the best copula, the AIC is used when comparing candidate copulas. Five widely used multivariate copulas are selected to fit the joint distribution of multiple degradation processes. The results of parameter estimation and AIC values are summarized in Table 6, and the joint reliability curves of tuner with different copulas are shown in Figure 9.  Figure 9. System reliability functions obtained using different multivariate copulas. Reliability Reliability Figure 9. System reliability functions obtained using different multivariate copulas. As can be seen in Table 6, the multivariate ADT model with Frank copula has the smallest AIC. In other words, Frank copula may best describe the dependence among the multiple accelerated degradation processes. On the other hand, the ADT model with the Gaussian copula has the largest AIC, and thus is not adequate to describe the dependency between the two performance parameters by linear correlation [28].
Finally, the s-dependent system reliability of tuner is compared with that obtained under an s-independent assumption. From Figure 10, one can see that the s-dependent reliability curve is higher than s-independent reliability curve. In other words, making the s-independent assumption may be under-estimate the system reliability. This is not uncommon in many engineering applications. Therefore, it is more reasonable to use the proposed method in product reliability estimation when the strong dependency among multiple performance parameters exhibits. As can be seen in Table 6, the multivariate ADT model with Frank copula has the smallest AIC. In other words, Frank copula may best describe the dependence among the multiple accelerated degradation processes. On the other hand, the ADT model with the Gaussian copula has the largest AIC, and thus is not adequate to describe the dependency between the two performance parameters by linear correlation [28].
Finally, the s-dependent system reliability of tuner is compared with that obtained under an s-independent assumption. From Figure 10, one can see that the s-dependent reliability curve is higher than s-independent reliability curve. In other words, making the s-independent assumption may be under-estimate the system reliability. This is not uncommon in many engineering applications. Therefore, it is more reasonable to use the proposed method in product reliability estimation when the strong dependency among multiple performance parameters exhibits.

Conclusions
In this paper, we have proposed a multivariate ADT model based on general Wiener process and copulas for a system with s-dependent nonlinear degradation processes. We first introduced a general Wiener process to describe univariate accelerated degradation process and provided the corresponding parameter estimation method. This model can cover the commonly used linear and time-scale transformed Wiener processes as its limiting cases. Then, the copulas were employed to develop a more flexible s-dependent multivariate ADT model, where the dependency among different degradation processes is linked by the copulas. The IFM method is applied to estimate the unknown parameters of copula and marginal distributions, and the AIC is used as the criterion to choose the best copula from several candidates. Finally, the case study validates the effectiveness and superiority of the general multivariate ADT model on solving real-world reliability engineering problems. It provides a useful tool for the reliability evaluation and the development of maintenance strategies for products with multiple performance parameters. This is also a new approach for information fusion in engineering practice. Our future research will be focused on the development of a model for analyzing multiple dependent nonlinear step-stress ADT data. Figure 10. Comparison of reliability functions of tuner and its performance parameters.

Conclusions
In this paper, we have proposed a multivariate ADT model based on general Wiener process and copulas for a system with s-dependent nonlinear degradation processes. We first introduced a general Wiener process to describe univariate accelerated degradation process and provided the corresponding parameter estimation method. This model can cover the commonly used linear and time-scale transformed Wiener processes as its limiting cases. Then, the copulas were employed to develop a more flexible s-dependent multivariate ADT model, where the dependency among different degradation processes is linked by the copulas. The IFM method is applied to estimate the unknown parameters of copula and marginal distributions, and the AIC is used as the criterion to choose the best copula from several candidates. Finally, the case study validates the effectiveness and superiority of the general multivariate ADT model on solving real-world reliability engineering problems. It provides a useful tool for the reliability evaluation and the development of maintenance strategies for products with multiple performance parameters. This is also a new approach for information fusion in engineering practice. Our future research will be focused on the development of a model for analyzing multiple dependent nonlinear step-stress ADT data.