# A New Variable-Censoring Control Chart Using Lifetime Performance Index under Exponential and Weibull Distributions

Authors: Muhammad Aslam; P. Jeyadurga; S. Balamurali; Rehan Ahmad Khan Sherwani; Mohammed Albassam; Chi-Hyuck Jun
Journal: Computational Intelligence and Neuroscience
 (2021)

---

## Abstract

In reliability theory or life testing, exponential distribution and Weibull distribution are frequently considered to model the lifetime of the components or systems. In this paper, we design a control chart based on the lifetime performance index using Type II censoring for exponential and Weibull distributions. Average run length helps to measure the performance of the proposed control chart. The optimal values of the number of failure items and decision criteria used to decide whether the process is in-control or out-of-control based on the sample results are determined such that the in-control average run length is as close as to the specified average run length values. We simulate the data to illustrate the performance of the proposed control chart.


---

## Body

## 1. Introduction

In all production processes, the most important task for each manufacturer is to maintain the process stability because products with uniform quality can be achieved only when the process is stable. Process stability is defined as a state in which a process has shown a certain degree of consistency in history and is expected to continue the same state in the future. To maintain process stability, the cause of process variation should be detected and controlled because the process may be affected by either common causes or by assignable causes. Under this situation, the most important tool of statistical process control, namely, the control chart, helps the manufacturer. Specifically, a control chart is a graphical representation that shows the process variation over a period of time. Control charts can be used to achieve and maintain process stability via monitoring, analyzing, and understanding the process variations. There are two types of control charts available such as attribute control chart and variable control chart. The former is applicable where the quality characteristic of the sample is nonmeasurable, and the latter is used to monitor the measurable quality characteristic. One can find more applications of control charts in the production field (see, for example, Bersimis et al. [1]; Zhu et al. [2]; Chong et al. [3]; Jeyadurga et al. [4]; Ali et al. [5]). It is to be mentioned that statistical distributions play an important role in designing control charts. In particular, normal distribution was frequently used when designing variable control charts in earlier years. But, in recent years, numerous distributions such as exponential, Weibull, gamma, etc. are involved in designing control charts (see, for example, Aslam et al. [6, 7]; Ali and Pievatolo [8]; Ali et al. [5, 9]). In this study, we consider the designing of the control chart under exponential and Weibull distributions.The lifetime is considered as the most essential characteristic of a product among other quality characteristics. For this reason, the producers concentrate on the lifetime of the products during the manufacturing process and a life test is conducted to determine the reliability of the products. However, the testing on the entire product is not possible due to time and cost constraints, whereas the test must be conducted to inspect the reliability. Under this situation, a time-censored life test (Type I) or failure-censored life test (Type II) is performed. In Type I censoring, the test termination time is prescribed in advance and it is observed that how many sample items failed until the time is attained. But, in Type II censoring, the number of failures is specified rather than the time. That is, the life test is terminated immediately if the specified number of failed items is attained. In this work, we consider Type II censoring in which a random sample ofn items is placed in the life test and the test is terminated at the failure of the sth sample item.Process capability indices measure the ability of an in-control process to produce the products within the specification limits. One can obtain the relationship between the actual performance of the process and the specification limit by such measure. These indices are considered process improvement techniques. Similarly, Montgomery [10] introduced the lifetime performance index CL, which is utilized to evaluate the lifetime performance of a process with respect to a lower specification limit L. Some authors investigated the uniformly minimum variance unbiased estimator and hypothesis testing procedure of CL under various distributions (see, for example, Tong et al. [11]; Wu et al. [12]; Lee et al. [13]; Lee et al. [14]; Lee et al. [15]; Lee [16]; Ahmadi et al. [17], Jafari and Bafekri [18]). In general, a normal distribution is considered as suitable to model the measurable quality characteristic of the products. But, the normal distribution is not appropriate to model the random variable that represents the time-to-failure (i.e., lifetime T) of electronic components since such random variable takes only positive values. In this circumstance, exponential or Weibull distribution is preferred as a suitable lifetime model rather than the normal distribution. Therefore, in this study, we consider exponential and Weibull distributions to evaluate the lifetime performance of a process. Numerous authors designed acceptance sampling plans using the lifetime performance index under exponential, gamma, and Weibull distributions (see, for example, Aslam et al. [19]; Wu et al. [20]; Lee et al. [21]). But, the study on designing control charts for exponential and Weibull distributions based on CL under Type II censoring is not available in the literature. Therefore, we attempt to design such a control chart to monitor the lifetime performance of a process. The performance of the proposed control chart is measured by average run length (ARL) and is defined as the number of plotted points in the control chart before the out-of-control indication. In this work, the simultaneous minimization of out-of-control ARL and the number of failures to be observed under inspection are considered. The performance of the proposed control chart is explained by using simulated data.


## 2. Designing of the Control Chart under Lifetime Distributions

In this section, we discuss the designing of the control charts under exponential and Weibull distributions based on the lifetime performance index.


### 2.1. Designing of the Proposed Chart under Exponential Distribution

Suppose that the lifetime of the products is defined by the random variableT and μT and σT are the mean and standard deviation of lifetimes, respectively. The index CL measures the performance of the process and is defined by (see Tong et al. [11])(1)CL=μT−LσT,where L is a lower lifetime limit. Suppose that the random variable T follows an exponential distribution with parameter λ > 0 and have the probability density function (pdf) as follows:(2)ft;λ=λe−λt,t≥0.For exponential distribution, the mean and standard deviation are same. The exponential distribution with pdf given in equation (2) has mean 1/λ (i.e., μT=1/λ) and standard deviation 1/λ (i.e., σT=1/λ). Hence, by the substitution of µT and σT in equation (1), CL can be simplified and rewritten as follows:(3)CL=1/λ−L1/λ=1/λ1/λ−L1/λ=1−λL,−∞<CL<1.The lifetime-nonconforming rate is defined by(4)p=PT≤L=1−e−λL=1−eCL−1.In practice,λ is usually unknown. In Type II censoring, s failures t1,t2,…,ts are noted from t1,t2,…,tn. Wu et al. [20] provided the estimator of λ under Type II right censoring which is given by(5)C^L=1−λ^L=1−s−1L∑i=1sti+n−sts.Suppose thatW=∑i=1sti+n−sts. For this case, s is predefined so the statistic W is sufficient for λ which leads that 2λW follows a chi-squared distribution with 2s degrees of freedom (i.e., 2λW∼χ2s2).When designing a control chart under exponential and Weibull distributions, we need two parameters such asH1 and H2 instead of control limits used in traditional control charts. That is, instead of finding control limit coefficients in the traditional designing procedure, the optimal values of H1 and H2 are found so that the in-control ARL is as close as to the specified ARL and the statistic obtained from the process is compared with these optimal values of such parameters. Then, the process is said to be in-control if the statistic lies between H1 and H2; otherwise, the process is declared as out-of-control. The construction of the proposed control chart is as follows.Step 1.
Select a random sample of sizen from the production process for the life test. Record first sn≥s failure times t1,t2,…,ts.Step 2.
CalculateC^L=1−s−1L/∑i=1sti+n−sts (see Wu et al. [20]). Declare the process is out-of-control if either C^L≥H2 or C^L≤H1. Declare the process is in-control state if H1≤C^L≤H2.
The probability that the process is to be declared as out-of-control is given by(6)Pout=PC^L≥H2+PC^L≤H1=1−PC^L≥H1+PC^L≥H2,Pout=1−P1−s−1L∑i=1sti+n−sts≥H1+P1−s−1L∑i=1sti+n−sts≥H2.
By following Wu et al. [20], Pout can be written as(7)Pout=1−Pχ2s2≥2s−11−CL1−H1+Pχ2s2≥2s−11−CL1−H2.
The probability that the process is declared to be out-of-control, when the process is atCL=CL0=1−λ0L, is given by(8)Pout0=1−Pχ2s2≥2s−11−CL01−H1|λ=λ0+Pχ2s2≥2s−11−CL01−H2|λ=λ0.
The in-control ARL of the proposed control chart whenCL=CL0 is denoted by ARL0 and given as follows:(9)ARL0=1Pout0.
Suppose that the manufacturing process has shifted fromλ0 to λ1 = kλ0, where k is a shift constant. The probability that the process is declared to be out-of-control when the process is at λ = λ1 is given by(10)Pout1=1−Pχ2s2≥2s−11−CL11−H1|λ=λ1+Pχ2s2≥2s−11−CL11−H2|λ=λ1,where CL1=1−λ1L=1−kλ0L=1−k1−CL0. (When preparing tables, specify CL0 and k. Then, CL1 can be obtained.)
The ARL for the shifted process is denoted by ARL1 and given by(11)ARL1=1Pout1.
To construct the proposed control chart for exponential and Weibull distributions based on the lifetime performance index under Type II censoring, we determine the optimal parameterss (number of failures to be observed), H1, and H2 such that the in-control ARL (ARL0) is very close to the specified ARL (say, r0). In addition, we minimize both ARL1 and s while finding the optimal parameters since the control chart is considered as effective if it detects the process shift quickly and uses the minimum number of observations to detect the shift. The following optimization problem can be utilized in determining the optimal parameters:(12)minimizeARL1andssubject toARL≥r0,n≥s,H2>H1.
The optimal parameters are determined for two specified in-control ARLs such as 300 and 370, and four different values for the lifetime performance index are considered asCL0 = 1.33, 1.50, 1.67, and 2.0 when the process is in-control. For various shift constant values, the out-of-control ARL (i.e., ARL1) is calculated using the optimal parameters and the shift constant k is taken from 0.9 to 0.1. Tables 1 and 2 provide the out-of-control ARL values along with the optimal parameters, and it should be noted that the performance of the process will decrease when there is a shift in the process. In particular, such performance decreases if there is a decrement in shift. The decreasing trend of out-of-control ARLs is observed from Tables 1 and 2 when the shift value decreases.Table 1
ARLs of the proposed control chart under exponential distributions whenr0 = 300.
ks = 5, H1 = 1.000, H2 = 1.100s = 3, H1 = 1.000, H2 = 1.102s = 3, H1 = 0.108, H2 = 1.137s = 3, H1 = 0.932, H2 = 1.204CL0 = 1.33CL0 = 1.5CL0 = 1.67CL0 = 2.01.0308.85307.58301.88307.580.9121.01139.32137.05139.320.849.3364.3563.4464.350.721.1030.4330.0730.430.69.5914.8314.6914.830.54.707.517.457.510.42.554.003.984.000.31.572.292.282.290.21.151.451.451.450.11.011.081.081.08Table 2
ARLs of the proposed control chart under exponential distribution whenr0 = 370.
ks = 5, H1 = 1.000, H2 = 1.098s = 4, H1 = 0.579, H2 = 1.127s = 5, H1 = 1.000, H2 = 1.199s = 3, H1 = 1.000, H2 = 1.199CL0 = 1.33CL0 = 1.5CL0 = 1.67CL0 = 2.01.0375.51377.22374.95376.280.9143.34154.12143.16166.390.856.9464.8656.8875.030.723.7428.3123.7234.650.610.5212.9310.5116.490.55.046.265.038.160.42.663.272.664.250.31.621.901.612.380.21.161.271.161.480.11.011.031.011.09


### 2.2. Designing of the Proposed Control Chart under Weibull Distribution

Suppose that the random variableT follows the Weibull distribution with shape parameter δ and scale parameter θ. The Weibull distribution is the generalization of the exponential distribution. The pdf of the Weibull distribution is given by(13)ft;θ,δ=δθtθδ−1exp−tθδ,t≥0,θ>0,δ>0.The mean (µT) and variance (σT) of the Weibull distribution are obtained as(14)μT=θΓ1+1δ,σT2=θ2Γ1+2δ−Γ21+1δ.The lifetime index for the Weibull distribution is given as(15)CL=μT−LσT=1AΓ1+1δ−Lθ.where −∞<CL<Γ1+1/δ/A and A=Γ1+2/δ−Γ21+1/δ.The lifetime-conforming rater is given by(16)r=PT≥L=exp−Lθδ=exp−Γ1+1δ−ACLδ=1−p,where p is the lifetime-nonconforming rate. In Type II censoring, s failures t1,t2,…,ts are noted from t1,t2,…,tn. Wu et al. [20] provided the estimator for CL when δ is known under Type II right censoring which is given by(17)C^L=1AΓ1+1δ−LΓsD1/δΓs−1/δ,where D=∑i=1sn−i+1tiδ−ti−1δ, and according to Wu et al. [20], for Type II censored data, we have 2∑i=1sZi=2θ−δD∼χ2s2. The probability that the process is declared to be out-of-control is given as in equation (6). Therefore,(18)Pout=1−P1AΓ1+1δ−LΓsD1/δΓs−1/δ≥H1+P1AΓ1+1δ−LΓsD1/δΓs−1/δ≥H2,Pout=1−Pχ2s2≥2Γ1+1/δ−ACLδΓδsΓ1+1/δ−AH1δΓδs−1/δ+Pχ2s2≥2Γ1+1/δ−ACLδΓδsΓ1+1/δ−AH2δΓδs−1/δ.The probability that the process is declared to be out-of-control when the process is atCL=CL0=1/AΓ1+1/δ−L/θ0 is given by(19)Pout0=1−Pχ2s2≥2Γ1+1/δ−ACL0δΓδsΓ1+1/δ−AH1δΓδs−1/δ+Pχ2s2≥2Γ1+1/δ−ACL0δΓδsΓ1+1/δ−AH2δΓδs−1/δ.The in-control ARL of the proposed control chart under Weibull distribution is obtained by using equation (10). Suppose that the manufacturing process has shifted from θ0 to θ1 = kθ0, where k is a shift constant. Then, the probability that the process is declared to be out-of-control when the process is at θ = θ1 is given by(20)Pout1=1−Pχ2s2≥2Γ1+1/δ−ACL1δΓδsΓ1+1/δ−AH1δΓδs−1/δ+Pχ2s2≥2Γ1+1/δ−ACL1δΓδsΓ1+1/δ−AH2δΓδs−1/δ.where CL=CL0=1/AΓ1+1/δ−L/θ0 and CL1=1/AΓ1+1/δ−L/θ1=1/AΓ1+1/δ−L/kθ0. The value of CL1 can be obtained fromCL0 as follows: CL1=1/AΓ1+1/δ1−1/k+CL0/k.One can obtain the ARL for the shifted process under Weibull distribution using equation (12). Also, the optimization problem given in equation (13) is used to determine the optimal parameters for the proposed control chart construction under Weibull distribution. The optimal parameters are determined under Weibull distribution for the same two specified in-control ARLs such as 300 and 370, and three different values CL0 = 1.33, 1.50, and 1.67 are considered as the case of exponential distribution. Here, the out-of-control ARL is calculated for shift constant k = 0.95(0.05)0.5. The optimal parameters and the corresponding out-of-control ARL values are reported in Tables 3 and 4. From these tables, it can be understood that the out-of-control ARLs decrease when the shift value decreases. In addition, it is observed that the ARL decreases if there is an increment in the required number of failure items.Table 3
ARLs of the proposed control chart under Weibull distribution with shape parameterδ = 2 when r0 = 300.
ks = 5, H1 = 0.676, H2 = 1.612s = 4, H1 = 0.927, H2 = 1.714s = 3, H1 = 1.214, H2 = 1.808CL0 = 1.33CL0 = 1.5CL0 = 1.671.00302.64303.42306.860.95219.49237.78258.950.90143.77167.78198.440.8591.24114.05146.140.8056.9976.05105.400.7535.1649.9174.720.7021.4732.2652.060.6513.0320.5635.600.607.9112.9423.890.554.848.0815.720.503.045.0510.15Table 4
ARLs of the proposed control chart under Weibull distribution with shape parameterδ = 2 when r0 = 370.
ks = 5, H1 = 0.644, H2 = 1.614s = 6, H1 = 1.088, H2 = 1.690s = 7, H1 = 1.457, H2 = 1.777CL0 = 1.33CL0 = 1.5CL0 = 1.671.00373.51370.50376.730.95271.45250.85236.850.90177.09153.11135.960.85111.7690.9176.450.8069.3653.3042.650.7542.4630.9723.730.7025.7017.9213.260.6515.4210.397.510.609.236.094.380.555.563.672.680.503.422.321.77


## 2.1. Designing of the Proposed Chart under Exponential Distribution

Suppose that the lifetime of the products is defined by the random variableT and μT and σT are the mean and standard deviation of lifetimes, respectively. The index CL measures the performance of the process and is defined by (see Tong et al. [11])(1)CL=μT−LσT,where L is a lower lifetime limit. Suppose that the random variable T follows an exponential distribution with parameter λ > 0 and have the probability density function (pdf) as follows:(2)ft;λ=λe−λt,t≥0.For exponential distribution, the mean and standard deviation are same. The exponential distribution with pdf given in equation (2) has mean 1/λ (i.e., μT=1/λ) and standard deviation 1/λ (i.e., σT=1/λ). Hence, by the substitution of µT and σT in equation (1), CL can be simplified and rewritten as follows:(3)CL=1/λ−L1/λ=1/λ1/λ−L1/λ=1−λL,−∞<CL<1.The lifetime-nonconforming rate is defined by(4)p=PT≤L=1−e−λL=1−eCL−1.In practice,λ is usually unknown. In Type II censoring, s failures t1,t2,…,ts are noted from t1,t2,…,tn. Wu et al. [20] provided the estimator of λ under Type II right censoring which is given by(5)C^L=1−λ^L=1−s−1L∑i=1sti+n−sts.Suppose thatW=∑i=1sti+n−sts. For this case, s is predefined so the statistic W is sufficient for λ which leads that 2λW follows a chi-squared distribution with 2s degrees of freedom (i.e., 2λW∼χ2s2).When designing a control chart under exponential and Weibull distributions, we need two parameters such asH1 and H2 instead of control limits used in traditional control charts. That is, instead of finding control limit coefficients in the traditional designing procedure, the optimal values of H1 and H2 are found so that the in-control ARL is as close as to the specified ARL and the statistic obtained from the process is compared with these optimal values of such parameters. Then, the process is said to be in-control if the statistic lies between H1 and H2; otherwise, the process is declared as out-of-control. The construction of the proposed control chart is as follows.Step 1.
Select a random sample of sizen from the production process for the life test. Record first sn≥s failure times t1,t2,…,ts.Step 2.
CalculateC^L=1−s−1L/∑i=1sti+n−sts (see Wu et al. [20]). Declare the process is out-of-control if either C^L≥H2 or C^L≤H1. Declare the process is in-control state if H1≤C^L≤H2.
The probability that the process is to be declared as out-of-control is given by(6)Pout=PC^L≥H2+PC^L≤H1=1−PC^L≥H1+PC^L≥H2,Pout=1−P1−s−1L∑i=1sti+n−sts≥H1+P1−s−1L∑i=1sti+n−sts≥H2.
By following Wu et al. [20], Pout can be written as(7)Pout=1−Pχ2s2≥2s−11−CL1−H1+Pχ2s2≥2s−11−CL1−H2.
The probability that the process is declared to be out-of-control, when the process is atCL=CL0=1−λ0L, is given by(8)Pout0=1−Pχ2s2≥2s−11−CL01−H1|λ=λ0+Pχ2s2≥2s−11−CL01−H2|λ=λ0.
The in-control ARL of the proposed control chart whenCL=CL0 is denoted by ARL0 and given as follows:(9)ARL0=1Pout0.
Suppose that the manufacturing process has shifted fromλ0 to λ1 = kλ0, where k is a shift constant. The probability that the process is declared to be out-of-control when the process is at λ = λ1 is given by(10)Pout1=1−Pχ2s2≥2s−11−CL11−H1|λ=λ1+Pχ2s2≥2s−11−CL11−H2|λ=λ1,where CL1=1−λ1L=1−kλ0L=1−k1−CL0. (When preparing tables, specify CL0 and k. Then, CL1 can be obtained.)
The ARL for the shifted process is denoted by ARL1 and given by(11)ARL1=1Pout1.
To construct the proposed control chart for exponential and Weibull distributions based on the lifetime performance index under Type II censoring, we determine the optimal parameterss (number of failures to be observed), H1, and H2 such that the in-control ARL (ARL0) is very close to the specified ARL (say, r0). In addition, we minimize both ARL1 and s while finding the optimal parameters since the control chart is considered as effective if it detects the process shift quickly and uses the minimum number of observations to detect the shift. The following optimization problem can be utilized in determining the optimal parameters:(12)minimizeARL1andssubject toARL≥r0,n≥s,H2>H1.
The optimal parameters are determined for two specified in-control ARLs such as 300 and 370, and four different values for the lifetime performance index are considered asCL0 = 1.33, 1.50, 1.67, and 2.0 when the process is in-control. For various shift constant values, the out-of-control ARL (i.e., ARL1) is calculated using the optimal parameters and the shift constant k is taken from 0.9 to 0.1. Tables 1 and 2 provide the out-of-control ARL values along with the optimal parameters, and it should be noted that the performance of the process will decrease when there is a shift in the process. In particular, such performance decreases if there is a decrement in shift. The decreasing trend of out-of-control ARLs is observed from Tables 1 and 2 when the shift value decreases.Table 1
ARLs of the proposed control chart under exponential distributions whenr0 = 300.
ks = 5, H1 = 1.000, H2 = 1.100s = 3, H1 = 1.000, H2 = 1.102s = 3, H1 = 0.108, H2 = 1.137s = 3, H1 = 0.932, H2 = 1.204CL0 = 1.33CL0 = 1.5CL0 = 1.67CL0 = 2.01.0308.85307.58301.88307.580.9121.01139.32137.05139.320.849.3364.3563.4464.350.721.1030.4330.0730.430.69.5914.8314.6914.830.54.707.517.457.510.42.554.003.984.000.31.572.292.282.290.21.151.451.451.450.11.011.081.081.08Table 2
ARLs of the proposed control chart under exponential distribution whenr0 = 370.
ks = 5, H1 = 1.000, H2 = 1.098s = 4, H1 = 0.579, H2 = 1.127s = 5, H1 = 1.000, H2 = 1.199s = 3, H1 = 1.000, H2 = 1.199CL0 = 1.33CL0 = 1.5CL0 = 1.67CL0 = 2.01.0375.51377.22374.95376.280.9143.34154.12143.16166.390.856.9464.8656.8875.030.723.7428.3123.7234.650.610.5212.9310.5116.490.55.046.265.038.160.42.663.272.664.250.31.621.901.612.380.21.161.271.161.480.11.011.031.011.09


## 2.2. Designing of the Proposed Control Chart under Weibull Distribution

Suppose that the random variableT follows the Weibull distribution with shape parameter δ and scale parameter θ. The Weibull distribution is the generalization of the exponential distribution. The pdf of the Weibull distribution is given by(13)ft;θ,δ=δθtθδ−1exp−tθδ,t≥0,θ>0,δ>0.The mean (µT) and variance (σT) of the Weibull distribution are obtained as(14)μT=θΓ1+1δ,σT2=θ2Γ1+2δ−Γ21+1δ.The lifetime index for the Weibull distribution is given as(15)CL=μT−LσT=1AΓ1+1δ−Lθ.where −∞<CL<Γ1+1/δ/A and A=Γ1+2/δ−Γ21+1/δ.The lifetime-conforming rater is given by(16)r=PT≥L=exp−Lθδ=exp−Γ1+1δ−ACLδ=1−p,where p is the lifetime-nonconforming rate. In Type II censoring, s failures t1,t2,…,ts are noted from t1,t2,…,tn. Wu et al. [20] provided the estimator for CL when δ is known under Type II right censoring which is given by(17)C^L=1AΓ1+1δ−LΓsD1/δΓs−1/δ,where D=∑i=1sn−i+1tiδ−ti−1δ, and according to Wu et al. [20], for Type II censored data, we have 2∑i=1sZi=2θ−δD∼χ2s2. The probability that the process is declared to be out-of-control is given as in equation (6). Therefore,(18)Pout=1−P1AΓ1+1δ−LΓsD1/δΓs−1/δ≥H1+P1AΓ1+1δ−LΓsD1/δΓs−1/δ≥H2,Pout=1−Pχ2s2≥2Γ1+1/δ−ACLδΓδsΓ1+1/δ−AH1δΓδs−1/δ+Pχ2s2≥2Γ1+1/δ−ACLδΓδsΓ1+1/δ−AH2δΓδs−1/δ.The probability that the process is declared to be out-of-control when the process is atCL=CL0=1/AΓ1+1/δ−L/θ0 is given by(19)Pout0=1−Pχ2s2≥2Γ1+1/δ−ACL0δΓδsΓ1+1/δ−AH1δΓδs−1/δ+Pχ2s2≥2Γ1+1/δ−ACL0δΓδsΓ1+1/δ−AH2δΓδs−1/δ.The in-control ARL of the proposed control chart under Weibull distribution is obtained by using equation (10). Suppose that the manufacturing process has shifted from θ0 to θ1 = kθ0, where k is a shift constant. Then, the probability that the process is declared to be out-of-control when the process is at θ = θ1 is given by(20)Pout1=1−Pχ2s2≥2Γ1+1/δ−ACL1δΓδsΓ1+1/δ−AH1δΓδs−1/δ+Pχ2s2≥2Γ1+1/δ−ACL1δΓδsΓ1+1/δ−AH2δΓδs−1/δ.where CL=CL0=1/AΓ1+1/δ−L/θ0 and CL1=1/AΓ1+1/δ−L/θ1=1/AΓ1+1/δ−L/kθ0. The value of CL1 can be obtained fromCL0 as follows: CL1=1/AΓ1+1/δ1−1/k+CL0/k.One can obtain the ARL for the shifted process under Weibull distribution using equation (12). Also, the optimization problem given in equation (13) is used to determine the optimal parameters for the proposed control chart construction under Weibull distribution. The optimal parameters are determined under Weibull distribution for the same two specified in-control ARLs such as 300 and 370, and three different values CL0 = 1.33, 1.50, and 1.67 are considered as the case of exponential distribution. Here, the out-of-control ARL is calculated for shift constant k = 0.95(0.05)0.5. The optimal parameters and the corresponding out-of-control ARL values are reported in Tables 3 and 4. From these tables, it can be understood that the out-of-control ARLs decrease when the shift value decreases. In addition, it is observed that the ARL decreases if there is an increment in the required number of failure items.Table 3
ARLs of the proposed control chart under Weibull distribution with shape parameterδ = 2 when r0 = 300.
ks = 5, H1 = 0.676, H2 = 1.612s = 4, H1 = 0.927, H2 = 1.714s = 3, H1 = 1.214, H2 = 1.808CL0 = 1.33CL0 = 1.5CL0 = 1.671.00302.64303.42306.860.95219.49237.78258.950.90143.77167.78198.440.8591.24114.05146.140.8056.9976.05105.400.7535.1649.9174.720.7021.4732.2652.060.6513.0320.5635.600.607.9112.9423.890.554.848.0815.720.503.045.0510.15Table 4
ARLs of the proposed control chart under Weibull distribution with shape parameterδ = 2 when r0 = 370.
ks = 5, H1 = 0.644, H2 = 1.614s = 6, H1 = 1.088, H2 = 1.690s = 7, H1 = 1.457, H2 = 1.777CL0 = 1.33CL0 = 1.5CL0 = 1.671.00373.51370.50376.730.95271.45250.85236.850.90177.09153.11135.960.85111.7690.9176.450.8069.3653.3042.650.7542.4630.9723.730.7025.7017.9213.260.6515.4210.397.510.609.236.094.380.555.563.672.680.503.422.321.77


## 3. Industrial Application

The contribution of control charts in producing high-reliability products through monitoring and detecting the process variation is an important one. Most of the industrialists have considered that this control chart is one of the useful tools to make the business as profitable. Suppose that the manufacturer of the electronic component wants that the lifetime of producing products should be very high and he decides to use the control chart for process monitoring. The proposed control chart under exponential distribution will be appropriate for monitoring the lifetime since it is considered that exponential distribution is the suitable model to represent the lifetime of electronic components. The performance of the proposed control chart in detecting the process shift is investigated by using simulated data, and data are given in Table5. Such data have been simulated for exponential distribution, for the following specifications: in-control ARL is 300 (i.e., r0 = 300) and sample size n = 30, and we assume that CL=CL0 = 1.67 (i.e., the process is in-control) where the lower limit is L = 0.387. The optimal values obtained for these specifications are as follows: s = 3, H1 = 0.108, and H2 = 1.137. To compute the performance index, the first three failure times of sample items in each of the 10 subgroups is considered, and such index is calculated and plotted in Figure 1. It can be observed from the figure that the performance index of each subgroup lies between decision criteria. Therefore, the process is declared to be in-control. Similarly, we can show the lifetime performance of the process by using simulated data where the lifetime follows Weibull distribution.Table 5
Simulated data for exponential distribution.
123456789100.0560.0200.0160.0100.1220.0210.0200.0760.0150.0100.1050.0370.0300.0500.1370.0500.0200.0810.0440.0510.1430.0420.0720.0560.1650.0820.1250.1020.0550.0600.2290.0540.2570.0870.1740.0980.1420.1050.0640.1040.3490.0780.2600.1470.2850.1940.1560.1120.0740.1080.3870.0810.2760.1720.2860.2060.2210.1480.1390.1110.4000.1350.2780.2570.2890.2540.2980.1810.1600.1320.4080.1660.2980.3130.2940.2800.3140.2810.1610.1660.4280.2370.3490.3390.4060.3660.3300.3290.2370.2060.4450.2750.4670.3720.4380.3970.3720.4700.2550.2150.4730.2860.5730.5330.4630.4550.4270.5090.2900.2210.5310.3160.6110.5800.5120.4740.5100.5130.3130.2860.5810.3190.6250.6170.5590.5610.6360.6180.3380.2970.6840.3270.6720.6220.7830.6330.6640.6220.4380.3120.7240.3740.6790.8180.8430.6510.6660.8080.5190.3380.7770.4240.8750.8420.8430.6520.7270.9160.5450.3540.7900.5270.8830.9880.9330.7030.8081.0360.5940.4660.7970.6811.0261.0210.9600.7150.9061.0710.6730.4980.9880.7881.0561.1571.0340.7580.9831.0770.7510.6841.0170.8161.0911.5361.0440.9791.0211.0890.9040.7151.0611.2251.1091.5671.0751.0151.0441.3940.9090.9371.1151.3681.1481.5781.3141.2131.2061.4760.9931.3061.1861.4341.1491.6341.3311.5341.2291.4811.0171.3821.3861.8231.2212.7571.3511.6351.2901.5571.1121.4881.4371.9251.2483.3431.4571.6401.3161.6281.2441.8781.6791.9801.9053.4171.7602.1481.4471.6901.6902.3531.8522.1721.9483.5411.8772.2611.7751.9042.3282.4132.3603.2113.0545.0712.2992.5712.0012.1752.3453.2352.8813.2333.4095.1342.3804.0462.0652.2792.7453.9372.9485.9874.7905.3174.5354.1492.3792.2825.0694.114C^L0.8140.3720.6250.5250.8410.6730.7810.7430.5160.555Figure 1
Performance of the proposed control chart.


## 4. Conclusion

A number of variable-control chart designing methods are available for monitoring the mean, standard deviation, etc. Similarly, some studies are available on designing the control charts for monitoring the performance of the manufacturing process using the capability index. We have also investigated the lifetime performance of the process using one of the capability indices, namely, lifetime performance index with the help of control chart. The designing of the proposed control chart for both the exponential and Weibull distributions has been discussed, and to reduce the time consumption, Type II censoring scheme is used. We have minimized the out-of-control average run length and the number of failure items which are some of the important parameters helpsto decide whether the process is in-control or out-of-control. The performance of the proposed control chart has been explained by using simulated data. It is concluded that the proposed control chart will be very helpful to the manufacturer for monitoring the performance of the process. The proposed chart has the limitation that it can be applied when the lifetime follows the Weibull distribution and cannot be applied for normally distributed data. In this study, designing of the control chart is given for monitoring the lifetime performance of the process where the lifetime follows exponential and Weibull distributions. As a result, there can be a future study on designing control charts for other lifetime distributions and also, an economic designing of the proposed control chart will also be considered.



---

*Source: 1350169-2021-12-27.xml*