International Journal of Photoenergy

Volume 2015, Article ID 243648, 8 pages

http://dx.doi.org/10.1155/2015/243648

## LED Lighting System Reliability Modeling and Inference via Random Effects Gamma Process and Copula Function

^{1}Department of Industrial Engineering, Southeast University, Nanjing 211189, China^{2}Department of Mathematics, Hubei Engineering University, Xiaogan 432100, China

Received 13 December 2014; Accepted 15 January 2015

Academic Editor: Ahmad Umar

Copyright © 2015 Huibing Hao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Light emitting diode (LED) lamp has attracted increasing interest in the field of lighting systems due to its low energy and long lifetime. For different functions (i.e., illumination and color), it may have two or more performance characteristics. When the multiple performance characteristics are dependent, it creates a challenging problem to accurately analyze the system reliability. In this paper, we assume that the system has two performance characteristics, and each performance characteristic is governed by a random effects Gamma process where the random effects can capture the unit to unit differences. The dependency of performance characteristics is described by a Frank copula function. Via the copula function, the reliability assessment model is proposed. Considering the model is so complicated and analytically intractable, the Markov chain Monte Carlo (MCMC) method is used to estimate the unknown parameters. A numerical example about actual LED lamps data is given to demonstrate the usefulness and validity of the proposed model and method.

#### 1. Introduction

In recent years, significant advances have been achieved in the manufacturing of light emitting diode (LED). At present, LED manufacturers have successfully designed LED as an indicator for many devices in the lighting systems due to the benefit of its low energy and long lifetime [1]. For its longer lifetime, higher reliability, and different failure mechanisms, it is difficult to evaluate and predict the reliability of LED. In order to study this, an understanding of the possible failure mechanisms of LED lamp should be examined.

In general, the failure mechanisms of a system can be categorized into catastrophic failures and degradation failures [2]. Catastrophic failures are characterized as being both complete and sudden. A failure is considered as a catastrophic failure which means the system ceases to operate. Degradation failures are characterized as being both partial and gradual. If a system is considered to have a degradation failure, the system still functions but does not meet the performance target. It is well known that the light output form power LED is highest when new and declines gradually over time; hence, that is a degradation failure. From the result in [3], we know that the common failure modes of LED can be divided into three failure modes, including catastrophic failure, intensity degradation failure, and chromatic change degradation failure. Because the probability of catastrophic failure under standard aging test conditions is very low [4], then, for the LED, the intensity degradation and chromatic change degradation are the domination failure modes.

Degradation (e.g., wear, erosion, and fatigue) is a common phenomenon for most mechanical systems or components. Degradation can be mathematically described with a continuous process in terms of time [5]. Because the stochastic process model can provide flexibility to describe the failure generating mechanisms and the operating environment characteristics, many articles use the stochastic process approach to model the degradation path, such as Markov chain, Wiener processes, and Gamma processes [6–8].

Most of the previous research which focuses on degradation analysis only considers one performance characteristic (PC) or component failure mechanism. In practice, modern products usually have complex structure and more functions. This means that modern highly reliable products may have multiple degradation failure mechanisms. A motivating example is a lighting system consisting of many light emitting diode (LED) lamps for different purposes of lighting. LED lamp has attracted increasing interest in the field of lighting systems due to its low energy and long lifetime. For different functions (i.e., illumination and color), it may have two or more PCs. In such situations, multivariate or at least bivariate degradation model is needed to estimate the reliability of products.

There are some works dealing with estimating system reliability of bivariate or multivariate degradation data [9–14]. These works use either independence assumption of the multiple PCs, multivariate normal distribution, or modeling with covariates and consider modification to single failure classifications. However, these assumptions may not match the engineering practice. In [15, 16], Sari et al. answered the question as to how one could quantify the reliability of LED lighting system which had two or more degraded PCs, and the dependence of these PCs was described through a copula function. But in Sari’s work, they modeled the degradation data with generalized linear regression model with population average approach. Compared with stochastic process model, the regression model ignores the temporal uncertainty of a degradation process, which results in limiting its applications. Hence, based on Sari’s work, Pan et al. [17, 18] used the Wiener process and the copula function to model multiple degradation failure mechanism, and a numerical example about fatigue cracks was presented.

From the definition of Wiener process, it is known that the degradation path is not a strictly increasing function. Generally speaking, a gamma process is better suited for modeling a degradation behavior that has a strictly increasing pattern. Crowder and Lawless [19] used a gamma process to illustrate their single-inspection policy for the maintenance of automobile brake pads. van Noortwijk [20] recently surveyed the application of gamma processes in maintenance. Tsai et al. [21] discussed the misspecification effect on the prediction of product’s MTTF when the degradation model is fitted by gamma and Wiener degradation processes. Bagdonavicius and Nikulin [22] discussed a joint model including gamma degradation process and catastrophic failures. Wang et al. [23] used bivariate gamma degradation process to model the degradation data and an adaptive method of residual life was obtained.

In addition, similar to Sari’s work, Pan et al. [17, 18] also only use the population average approach and do not consider the random effect in their model. In fact, considering that each product possibly experiences different sources of variations during its operation, for a degradation model to be realistic, it is more appropriate to incorporate unit to unit variability in the degradation process. The degradations of such products can be described by random effects models, where the random effects can capture the unit to unit differences. Particularly for the actual LED lighting system, there are few researches using the copula method with the random effects Gamma process to assess the products reliability.

From the above analysis, it is clear that the reliability assessment of the system with two PCs has not been studied thoroughly. In this paper, we assume that the LED lighting system has two PCs (i.e., illumination and color), and each PC is governed by a random effects Gamma process. Moreover, we assume that the two PCs are dependent and their dependency can be characterized by a copula function. Since the likelihood function in such a situation is quite complicated, the Markov chain Monte Carlo (MCMC) algorithm is used to obtain the unknown parameters.

The rest of the paper is organized as follows. In Section 2, some assumptions and copula basic are described. Then, the bivariate degradation model based on random effects Gamma process and the copula function is introduced in Section 3. In Section 4, the estimation of unknown parameters based on the MCMC algorithm is obtained. A numerical example about actual LED lamps data is given in Section 5. Finally, some conclusions are made in Section 6.

#### 2. Some Basic Assumptions and Copula Brief Introduction

##### 2.1. Basic Assumptions

To analyze the problem statistically, some assumptions are used for the reliability modeling in this paper. The details of each assumption are explained in the corresponding sections.(1)The samples are independent, and no catastrophic failures occur during the degradation process.(2)The marginal degradation processes can be modeled as Gamma process with random effect.(3)For a specific product, the degradation measurements on the two PCs are observed at the same time (balanced data).(4)A product is supposed to have failed if one of the PCs exceeds the corresponding failure threshold for the first time. Here, the failure threshold vector is denoted by .(5)The two PCs are dependent on each other, and the dependency can be characterized by a Frank copula function.

##### 2.2. Copulas and Their Properties

###### 2.2.1. Definition and Base Properties

Copula function is a powerful tool to model the dependence structure among multiple PCs since they (see Nelsen [24]) give the dependence structure which relates the known marginal distributions of each PC to their multivariate joint distribution.

A two-dimensional copula is conventionally defined as a bivariate cumulative distribution function with uniform margins. A probabilistic way to define the copula is provided by the theorem of Sklar (see Nelsen [24]).

Theorem 1. *Let and be random variables with continuous distributions and , respectively, and let be the two-dimensional cumulative distribution function. Then, there exists a two-dimensional copula such that, for all ,
**If and are continuous, then is unique; conversely, if is a copula and and are distribution functions, then the function defined by (1) is a joint distribution function with margins and .*

*2.2.2. Measures of Association*

*We will consider here the standard dependence measures, Spearman’s and Kendall’s . These measures are related to the copula since the latter is an expression of the stochastic relationship between and within the entire range of values the variables can take. It is not difficult to show that
where the range of and can be shown to be []. For further properties of and , see Nelsen [24].*

*2.2.3. Archimedean Copula*

*Archimedean copulas have a wide range of applications, because(a)they can be constructed easily;(b)they have many nice properties;(c)a lot of copulas belong to this class;(d)they can be easily extended from 2-dimension to -dimension when satisfying some conditions.*

*In this paper, the Frank copula which belongs to Archimedean copula family is used to depict the dependence among multiple PCs as follows:
where is the Frank copula parameter and . And the relationship between Kendall’s τ and the Frank copula parameter is given by
where is a Debye function of the first kind (see Nelsen [24]).*

*3. Reliability Model Based on Gamma Process with Random Effects*

*3.1. Marginal Reliability Model Based on Random Effects Gamma Process*

*As mentioned earlier, a gamma process is more suitable for describing a monotone increasing degradation path. A well-adopted form for the regular Gamma process can be expressed as
where is a Gamma distribution with shape parameter and scale parameter , and the corresponding probability density function (PDF) is
where is Gamma function, and the Gamma process has the following properties:(1) with probability one;(2) has independent increments; that is, and are independent, ;(3), .*

*Assume that the degradation path of a product is governed by (5). Given the threshold value , the product’s lifetime is defined as
*

*Since the Gamma process has a monotone path, then the cumulative distribution function (CDF) of the lifetime can be expressed as
*

*From the properties of the Gamma process, we know that the average performance degradation characteristic in (5) is a linear function about . When it is not linear, Whitmore and Schenkelberg [25] used a transformed time scale to make it linear. Generally, the transformation is denoted by
*

*Then, a nonstationary Gamma process can be written as
*

*Considering that each item possibly experiences different sources of variations during its operation, for a degradation model to be realistic, it is more appropriate to incorporate item to item variability in the degradation process. In this paper, for simplicity, we assume that and are the fixed parameters which are common to all products; is a random parameter representing the heterogeneity among different products. Then, the random effects Gamma process can be written as
where is a Gamma distribution with shape parameter and scale parameter , and the corresponding probability density function (PDF) of is
*

*Then, the PDF of in (11) is given by
*

*From (13), for any fixed t, the random variable has an distribution with and degrees of freedom, and we denote as the CDF of .*

*Then, the CDF of in (11) is given by
*

*Then, the CDF of the lifetime can be expressed as
*

*Suppose that a product has two PCs and each PC is governed by a random effects Gamma process. The th PC is defined by (11) as
where and .*

*Let be the threshold value of the th PC; the lifetime of the th PC is defined as
*

*Then, the failure time distribution of the th PC can be obtained as
and the reliability function is
*

*3.2. Degradation Model Based on Bivariate Degradation Data*

*Suppose that a product has two PCs and the measurements of degradation process are random variables at observation time point . From the basic assumptions, the product is considered to be failed if any one PC reaches its corresponding failure threshold, which is known as . Therefore, the product still works when each PC keeps below its failure thresholds. Given the failure time of the th PC, suppose that the lifetime of the system is , and . Then, the product reliability can be written as follows:
*

*If the two degradation failure mechanisms are assumed to be independent, the product reliability in (20) can be rewritten as
*

*From (19) and (21), the product reliability with two PCs can be easily obtained. However, if the degradation failure mechanisms are not independent of each other, we cannot provide accurate product reliability estimation. Therefore, the copula method is utilized to establish the dependent structure among various degradation measurements.*

*Copulas provide a very convenient way to model and measure the dependence among multiple degradation failure mechanisms. The advantage of the copula function is that the joint reliability function can be modeled directly through the univariate marginal reliability functions of the individual failure processes, and the copula has no constraints on the univariate marginal distribution.*

*Suppose that is the CDF of lifetime for each PC, and is the joint CDF of and . According to Sklar’s theorem, there exists a unique copula such that
where is the parameter vector of the copula function.*

*Then, similarly to [16, 18], the system reliability in (20) can be obtained as
*

*If the product has two PCs linked by bivariate Frank copula given in (3), then we can obtain the system reliability function as
*

*4. Statistical Inferential Methods for Unknown Parameters*

*Here, we discuss the estimation of parameters required to implement the reliability function in (24). The unknown parameters are
*

*From (18), (19), and (24), we know that the model not only has nine parameters, but also is very complicated from a computational viewpoint. For this reason, the MCMC with the Gibbs sampling techniques is employed in this study to estimate model parameters.*

*The MCMC method is a simulation technique in which the analytical posterior distribution is difficult to be computed. A Markov chain is generated by sampling the current point based on the previous one. The MCMC method works successfully in Bayesian computing. By using the MCMC method, it is possible to generate samples from the posterior distribution and to use these samples to estimate the desired features of the posterior distribution. In addition, the simulation algorithm can be easily extensible to models with a large number of parameters or high complexity. The MCMC techniques, including the Metropolis-Hastings (M-H) algorithm and the Gibbs sampler, have become very popular in recent years as methods for generating a sample from a complicated model (see the details in [26, 27]).*

*The Gibbs sampler is a special case of an MCMC algorithm. It generates a sequence of samples from the full conditional probability distributions of two or more random variables. Gibbs sampling requires decomposing the joint posterior distributions into full conditional distributions for each parameter in the model and then sampling from them. We can use the Gibbs sampler to generate a sample, and then the unknown parameters are estimated. Let denote the full conditional posterior distribution of , where and is the observed data.*

*Then, the MCMC algorithm with the Gibbs sampling can be summarized as follows. Step 1.Initialize .Step 2.Set .Step 3.Generate from conditional distribution .Step 4.Generate from conditional distribution .Step 5.Generate from conditional distribution .Step 6.Generate from conditional distribution .Step 7.Set , and repeat steps 3–6, .Step 8.Estimate desired features based on the simulate sample , where denotes the burn-in number. For example, we can use the mean of to be the estimator of .*

*Using the Bayesian software package WinBUGS (see [28]) and carrying out the Gibbs sampling, the estimator of the model parameters can be obtained.*

*5. Numerical Example: Application to a LED Lighting System*

*In recent years, LED has attracted increasing interest in the field of lighting systems owing to its high efficiency, low energy, and long lifetime. Due to its different functions (i.e., illumination and color), there may have been multiple failure mechanisms about the LED lighting system. The actual LED lamps data are taken from Chaluvadi [29], where the LED lamps are put into test under an accelerated current of 40 mA. In the original data, 12 samples are tested for LED data and the measurements are taken at the same measurement times. The measured frequency of its light intensity is 50 hours. For demonstrating the bivariate degradation model, similarly to Sari et al. [16] and Pan et al. [18], we choose 12 samples and the data will be treated as if half of it is the first PC and the other half represents the second PC. That is to say, in the LED lighting system, the performance demand of the first group LED lamps is for lighting, and the performance demand of the second group LED lamps is for chromatic change. The used data are the data measured only until 250 hours. The LED lighting system is considered to have failed if one of the two LED datasets is more than 50. Table 1 lists the LED data.*