Reliability Analysis of an Extended Shock Model and Its Optimization Application in a Production Line

Abstract

Firstly, a two-unit cold standby shock model with multiple adaptive vacations is introduced, in which the startup and replacement of repair facility are also considered. Secondly, using supplementary variable method and Laplace transform, some important reliability indices are derived, such as availability, failure frequency, mean vacation period, mean renewal cycle, mean startup period, and replacement frequency. Finally, a production line controlled by two cold-standby computers is modeled to present numerical illustration and its optimal part-time job policy at a maximum profit.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Academic Editor:Xian-Jun Long

1, School of Mathematics and Statistics, Chongqing University of Technology, Chongqing 400054, China

Received 10 February 2014; Accepted 21 March 2014; 13 May 2014

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.

1. Introduction

It is well known that the shock model is used to study the external causes which may make a system fail. For example, a computer system may fail due to the invasion of some virus or an attack from a raider. Many authors have investigated various shock models based on different assumptions. Among them, Esary et al. [1], Barlow and Proschan [2], Ross [3], and Fagiuoli and Pellerey [4] dealt with Poisson shock models. Eryilmaz [5, 6] discussed discrete-time shock model and its life behavior. Also, Gottlieb [7], Aven and Gaarder [8], Wang and Zhang [9], Lam and Zhang [10], and Tang and Lam [11] obtained the optimal replacement policies for several shock models. Recently, from the reliability viewpoint, Li and Zhao [12] studied a δ -shock model consisting of n components, and Q. T. Wu and S. M. Wu [13] analyzed a two-unit cold standby shock model with single vacation.

However, in the existing shock literature, the repair facility is assumed to be always available when a failed unit occurs, although this assumption is evidently unrealistic. In fact, in many practical situations the repair facility generally needs a startup time with random length for its preparatory work before starting repair. Furthermore, the busy repair facility is typically subject to lengthy and unpredictable breakdowns and has to be replaced (see [14-16]). On the other hand, to utilize the repairman's idle time effectively and increase profit, the system manager can assign some secondary jobs to the idle repairman. But the repairman's additional tasks will reduce system availability and sometimes yield huge economic losses. Therefore, it is very important that the system manager knows how to assign the idle repairman optimal additional jobs based on a maximum profit and high availability level. In this paper, the period when the repairman undertakes additional jobs is represented by the repairman's vacation time. A comprehensive and excellent study on the vacation models can be found in Tian and Zhang's book [17].

Based on the above facts, in this paper we present an extended shock model for two identical unit cold-standby systems. Here, cold standby means that the redundant unit cannot fail at its standby state. Our study differs from previous work [1-13] in that (i) it considers the startup and breakdown of repair facility and (ii) it introduces the multiple adaptive vacation policy (MAVP). The MAVP first proposed by Tian and Zhang [17] is more general than single vacation, multiple vacations and variant vacations, which is useful for high availability and profit optimization of the system (see Section 5); (iii) some new reliability indices are derived, such as mean renewal cycle, mean vacation period, mean startup period, mean idle period, and mean busy period; (iv) as an application of our model, a production line controlled by two identical cold-standby computers is modeled to analyze its optimal part-time job policy.

The paper is organized as follows. Section 2 presents the assumptions. In Sections 3 and 4 we obtain the solutions to state probability equations and main reliability indices by the supplementary variable method. In Section 5 the numerical illustration and optimal part-time job policy for a production line are given. Conclusions are drawn in Section 6.

2. Assumptions

The extended shock model we consider here consists of two identical cold-standby units, a repair facility and a repairman. The model assumptions are as follows.

(1) The external shocks arrive according to a Poisson process with the rate λ(>0) . The magnitude of each shock, χ , is independent with common distribution function T . Shocks only influence the operating unit. The operating unit will fail if χ outstrips a threshold [varsigma] , where [varsigma] is assumed to be nonnegative with a distribution function [varphi] .

(2) Suppose that shocks are the only cause of unit failure. The system breaks down if and only if two units fail. When the operating unit fails, the standby one begins to operate if there is one (the switch is instantaneous and perfect). The failed units are repaired in order of failures. The repairman can repair only one failed unit at a time by means of repair facility. The repaired unit is as good as new and immediately goes into standby or operating state. The repair times are independent and identically distributed (i.i.d.) with common distribution function G(t) (t...5;0) , density function g(t) (t...5;0) , hazard rate function μ(t) (=g(t)/(1-G(t))) , and finite mean E(Y) , respectively.

(3) After completing repair, the repairman will shut down repair facility and take a multiple adaptive vacation policy. Under this vacation policy, the length of each vacation is i.i.d. random variable, which follows a general distribution function V(t) (t...5;0) with density function v(t) (t...5;0) hazard rate function r(t) (=v(t)/(1-V(t))) , and mean vacation time E(V) , respectively. The maximum vacation number of the repairman, denoted by H , has an arbitrary distribution ^{{_ci }i=1∞} with probability generating function H(z)=^∑i=1∞ ..._ci^zi , |z|<1 . At each vacation completion instant, the repairman checks the system and decides the action to take according to the system state. There are three possible cases: (A) if there is any failed unit in the system, he will immediately spend a startup time Q to turn repair facility on and then start his repair until there are no failed units; (B) if there is no failed unit in the system and the total number of vacations is still less than H , he will take another vacation; (C) if there is no failed units in the system and the total number of vacations is equal to H , he will remain idle in the system until the first failed unit appears, which induces a startup time Q and subsequent repair. We assume that the startup time Q has distribution function Q(t) (t...5;0) , density function q(t)(t...5;0) , hazard rate function θ(t) (=q(t)/(1-Q(t))) , and mean startup time E(Q) , respectively.

(4) The repair facility may break down with a Poisson rate α in the process of repair. The broken facility is immediately replaced by the repairman. The replacement times are assumed to be i.i.d. random variables having a general distribution function B(t)(t...5;0) , density function b(t)(t...5;0) , hazard rate function β(t)(=b(t)/(1-B(t))) , and mean replacement time E(B) , respectively. After replacement, the repair facility continues its remaining repair. The repair time of the failed unit is cumulative.

(5) Initially, both units are new (one is operating and the other is in cold standby), and the repairman is idle. After the first busy period is completed, he begins to take a multiple adaptive vacation policy. All random variables are mutually independent.

Remark 1.

From Assumptions (1) and (5), the probability that a shock causes the operating unit to fail is given by [figure omitted; refer to PDF] where Pr(Ω) is the probability of event Ω .

3. The State Probability Equations and Solutions

We define the possible states of the system as follows: (0) osv¯ , (1) oo¯ v¯¯ , (2) o¯ o¯ v¯¯ , (3) oo¯¯¯ , (4) o¯ o¯¯¯ , (5) oo¯¯ v¯¯¯ , (6) o¯ o¯¯ v¯¯¯ , (7ij) os_vij ,j=1,2,...,i ; i=1,2,... , (8ij) oo¯_vij , j=1,2,...,i ; i=1,2,... , (9ij) o¯ o¯_vij , j=1,2,...,i ; i=1,2,... , where o , s , o¯ , o¯¯ , and o¯¯¯ represent that one unit is operating, in cold standby, waiting for repair, waiting for remaining repair (preserving the time spent in repair), and under repair, and v¯,v¯¯ , v¯¯¯ , and _vij represent that the repairman is idle, turning the repair facility on, replacing the repair facility, and taking the j th vacation under the condition that the maximum vacation number is i , j=1,2,...,i ; i=1,2,... , respectively. By definition, the system states (0) , (1) , (3) , (5) , (7ij) , and (8ij) are operable and (2) , (4) , (6) , and (9ij) are inoperable.

Let S(t) be the system state at time t . For t...5;0 , we define X(t) as the elapsed startup time of repair facility at time t , Y(t) the elapsed repair time of the failed unit at time t , Z(t) the elapsed replacement time of the broken repair facility at time t , and γ(t) the elapsed vacation time of the repairman at time t . Then, {S(t),X(t),Y(t),Z(t),γ(t),t...5;0} is a vector Markov process. Define the state probabilities at time t as follows: [figure omitted; refer to PDF] where x , y , z , and w are the values taken by X(t) , Y(t) , Z(t) , and γ(t) , respectively.

In steady state, we define [figure omitted; refer to PDF]

Let ρ=pλ , since the process {S(t),X(t),Y(t),Z(t),γ(t),t...5;0} is a vector Markov process in continuous time, one can write the equations of the process in the usual way by considering the transitions occurring in t and t+Δt . For example, we have [figure omitted; refer to PDF] Letting Δt tend to zero yields [figure omitted; refer to PDF] By taking the limit t[arrow right]∞ in (^{4[variant prime]} ) , we can obtain [figure omitted; refer to PDF] In the same way, we readily get the following steady-state equations for state probabilities: [figure omitted; refer to PDF] with the boundary conditions [figure omitted; refer to PDF] and the normalization condition [figure omitted; refer to PDF]

In order to derive important reliability indices, we define the Laplace transform of a nonnegative function F(t) as ^F* (s)=^∫0∞ ...^e-st F(t)dt , and we also denote F¯(t)=1-F(t) , ^F¯* (s)=^∫0∞e-st F¯(t)dt. Define [figure omitted; refer to PDF] then we have the following.

Lemma 2.

The expressions of _πm , m=0,1,...,9 are given as follows: [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

Proof.

Solving (6) and (10), respectively, and using (19), we get [figure omitted; refer to PDF] Applying (28), (7) becomes [figure omitted; refer to PDF] which leads to [figure omitted; refer to PDF] Similarly, by the theory of first-order, linear, ordinary differential equation, we get [figure omitted; refer to PDF]

Now, we need to determine the values of _P0 , _Pm (0) , m=1,2,3,4 , and _P7ij (0) , j=1,2,...,i ; i=1,2,... .

Substituting (32) and (35) into (21) and (22), respectively, leads to [figure omitted; refer to PDF] which gives rise to [figure omitted; refer to PDF] From (5), (35), and (40), we obtain [figure omitted; refer to PDF] Substituting (42), (36), and (41) into (15) gives [figure omitted; refer to PDF] Combining (16), (37), and (41) leads to [figure omitted; refer to PDF] By (17) and (18), we get [figure omitted; refer to PDF] It follows from (45) that [figure omitted; refer to PDF] Finally, according to (47) and (43), we have [figure omitted; refer to PDF] Thus, _πm , m=0,1,...,9 are obtained by (28)-(37) and the above results. Again by the normalization condition (24), we get the value of _P3 (0) .

4. Reliability Indices

Based on the results obtained in Lemma 2, some important reliability indices are easily derived as follows.

Theorem 3.

(1) The steady-state availability of the system, that is, in steady state, the probability that the system is operating, denoted by A , is [figure omitted; refer to PDF]

(2) The steady-state failure frequency of the system, that is, in steady state, the rate of occurrence of failures of the system, is given by [figure omitted; refer to PDF]

(3) In steady state, let _PI and _PV denote the repairman's idle and vacation probabilities, and _PS and _PB are the startup and busy (repair and replacement) probabilities of repair facility, respectively; then [figure omitted; refer to PDF]

(4) In steady state, let ξ denote the unavailability of the repair facility, that is, the probability that the repair facility is under replacement, and η is the replacement frequency of the repair facility, that is, the rate of occurrence of replacements; then [figure omitted; refer to PDF] where _P3 (0) is given by Lemma 2.

Proof.

According to the frequency formula in [18], W and η are easily obtained. The rest of probabilities are derived by their definitions, Lemma 2, and direct calculations.

Noting that the time points that the repairman completes repair and begins to take vacation are regenerative ones, we define the following.

(i) System renewal cycle denoted by _LR : this is the length of time from the beginning of the last vacation period to the beginning of the next vacation period.

(ii) Vacation period denoted by _LV : this is the length of total vacation time of repairman per renewal cycle.

(iii): Idle period denoted by _LI : this is the length of total idle time of repairman per renewal cycle.

(iv) Startup period denoted by _LS : this is the length of total startup time of repair facility per renewal cycle.

(v) Busy period denoted by _LB : this is the sum of repair time and replacement time of repair facility per renewal cycle.

Theorem 4.

Denote by E(_LV ) and E(J) the mean vacation period and mean actual vacation number of repairman, respectively; then [figure omitted; refer to PDF]

Proof.

Let _L¯V (t)=Pr(_LV ...5;t) denote the complementary distribution of vacation period _LV ; then E(_LV )=^L¯V* (0) . Now, we will construct a vector Markov process with finite absorbing states to get ^L¯V* (s) .

Let all system states defined in Section 3 except (kij) , k=7,8,9 ; j=1,2,...,i ; and i=1,2,... be absorbing states; then we obtain a new vector Markov process. For this new process, using similar state probability notations to Section 3, we get the following set of equations: [figure omitted; refer to PDF] with the boundary conditions [figure omitted; refer to PDF] where δ(t)=1 if t=0 , δ(t)=0 if t...0;0 , and the initial conditions: _P7i1 (0,w)=_ci δ(w) , i=1,2,... , the rest are zero.

These equations can be solved in a similar manner to that in Section 3. The Laplace transforms of solutions are given by [figure omitted; refer to PDF] Therefore, we have [figure omitted; refer to PDF] Substituting the values of ^Pkij* (s,w) , k=7,8,9 ; j=1,2,...,i ; and i=1,2,... into the above equation and direct calculations complete the proof of Theorem 4.

Corollary 5.

Denote by E(_LR ) , E(_LI ) , E(_LS ) , and E(_LB ) the mean renewal cycle, mean idle period, mean startup period, and mean busy period, respectively; then [figure omitted; refer to PDF] where _P3 (0) is given by Lemma 2.

Proof.

The results are easily derived by noting that [figure omitted; refer to PDF]

Remark 6.

Since the time points in which the repairman completes repair and begins to take vacation are regenerative ones, by the frequency formula in [18] and _P3 (y) in (32), the steady state renewal frequency of system, that is, the rate of occurrence of renewals of system, is given by κ=^∫0∞ ...μ(y)_P3 (y)dy=_P3 (0)^g* (Λ) . Again by Proposition 3.3 in [19], we get E(_LR )=1/κ=1/(_P3 (0)^g* (Λ)) . The result derived by this method agrees with that in Corollary 5, which confirms that the deduction processes in Theorem 4 and Corollary 5 are right. More importantly, Theorem 4 presents an effective analysis technique to derive the mean vacation period and mean actual vacation number.

Remark 7 (special cases).

(1) If p=0 , then each shock does not harm the operating unit; If p=1 , then each shock will cause the operating unit to fail.

(2) If α[arrow right]0 , then our model becomes a two-unit cold standby shock system with repair facility startup and multiple adaptive vacations.

(3) If we let Q(t)=1-^e-θt (θ>0, t...5;0) , then when θ[arrow right]∞ , α[arrow right]0 , our model becomes a two-unit cold standby shock system with multiple adaptive vacations.

(4) If we let H=k and H(z)=^zk , |z|<1 , k=0,1,2,... , then our model can describe a two-unit cold standby shock system with repair facility startup and replacement under a variant vacation policy. In this special case, H=0,1,k and ∞ represent no vacation, single vacation, k vacations, and multiple vacations, respectively.

5. Numerical Illustration and Optimal Part-Time Job Policy for a Production Line

As an application example of our model and its results, we consider a production line controlled by a huge operating computer, which may fail due to the invasion of some virus. The virus arrives according to a Poisson process with rate λ (>0) . If the magnitude of a shock χ outstrips a threshold [varsigma] the operating computer will fail, where χ and [varsigma] are assumed to have distribution functions T(x) and [varphi](x) , x>0 , respectively. When the operating computer fails, an identical standby computer begins to operate. If two computers fail the production will be stopped. The repairer repaired the failed computer(s) in order of failure at a repair station. Before repair the repair station needs a random startup time Q . After startup, the repairer starts his repair until no failed computer appears. The repair time of each computer G is assumed to be a random variable. The repaired computer goes into operating or standby state. Whenever no failed computer occurs, the repairer shuts down the repair station and performs H part-time jobs (vacations), where H is a discrete random variable with generating function H(z), |z|<1 . The part-time jobs can make profits for the production line. Upon completion of each part-time job time V , the repairer checks the computer failure and decides whether or not to resume startup and subsequent repair. If at this moment there is no failed computer, a decision may be made for taking another part-time job to be performed. If the failed computer occurs, he will restart the startup of repair station before repair. If H part-time jobs are completed, the repairer will be idle and wait for the computer failure. Further, the repair of failed computer may be interrupted due to some unpredictable events, which occur according to a Poisson process with rate α (>0) . The interrupted repair is immediately recovered with a random time B . The repair will restart when the interruption is recovered.

Firstly, we numerically analyze the effects of part-time job, virus arrival rate, startup, and repair recovery time on main reliability indices of production line. These indices include the stop production frequency W , mean renewal cycle E(_LR ) , mean startup period E(_LV ) , mean part-time job period E(_LV ) , the availability A , busy probability _PB , unavailability ξ , and repair recovery frequency η .

For convenience, we assume that (1) the time of each part-time job V is exponentially distributed with parameter r (>0) , (2) the repair time G is deterministic with mean E(G)=0.25 , (3) the startup time Q obeys 2-stage Erlang distribution with parameter θ (>0) , and (4) the repair recovery time B follows hyperexponential distribution with density function b(t)=0.4β^e-βt +0.6(2β)^e-2βt , β,t>0 . Also, suppose that T(x)=1-^e-σx and [varphi](x)=1-^e-(1/2)σx , σ,x>0 ; then we get p=^∫0∞ ...[varphi](x)dT(x)=1/3 . The numerical results are reported in Tables 1-5 by varying values of H , r , λ , θ and β , where in Tables 2-5 it is assumed that H is geometrically distributed with mean 5 .

Table 1: Main reliability indices of production line for different maximum part-time job number H with (λ,r,θ,α,β)=(3,15,6,0.35,14) .

H	ξ	η	A	W	P B	E ( L R )	E ( L S )	E ( L V )
0	0.0057	0.0802	0.9171	0.2614	0.2350	1.5251	0.1667	0.0000
1	0.0057	0.0802	0.9163	0.2628	0.2348	1.5301	0.1667	0.0667
3	0.0057	0.0800	0.9148	0.2651	0.2344	1.5392	0.1667	0.1878
5	0.0057	0.0799	0.9135	0.2672	0.2341	1.5473	0.1667	0.2942
∞	0.0057	0.0791	0.9046	0.2817	0.2318	1.6055	0.1667	1.0667

Table 2: Main reliability indices of production line for different part-time job rate r with (λ,θ,α,β)=(3,6,0.35,14) .

r	ξ	η	A	W	P B	E ( L R )	E ( L S )	E ( L V )
5	0.0056	0.0783	0.8946	0.2861	0.2292	1.6434	0.1667	0.6000
8	0.0057	0.0794	0.9070	0.2751	0.2324	1.5825	0.1667	0.4327
11	0.0057	0.0797	0.9114	0.2700	0.2335	1.5592	0.1667	0.3409
14	0.0057	0.0799	0.9134	0.2673	0.2341	1.5477	0.1667	0.2820
17	0.0057	0.0800	0.9145	0.2657	0.2344	1.5412	0.1667	0.2406

Table 3: Main reliability indices of production line for different virus arrival rate λ with (r,θ,α,β)=(15,6,0.35,14) .

λ	ξ	η	A	W	P B	E ( L R )	E ( L S )	E ( L V )
1	0.0021	0.0288	0.9869	0.0404	0.0843	3.4655	0.1667	0.3067
2	0.0040	0.0577	0.9553	0.1384	0.1632	2.0062	0.1667	0.2848
3	0.0057	0.0800	0.9139	0.2667	0.2342	1.5452	0.1667	0.2667
4	0.0072	0.1013	0.8686	0.4077	0.2968	1.3327	0.1667	0.2513
5	0.0086	0.1200	0.8228	0.5508	0.3514	1.2188	0.1667	0.2381

Table 4: Main reliability indices of production line for different startup time parameter θ with (λ,r,α,β)=(3,15,0.35,14) .

θ	ξ	η	A	W	P B	E ( L R )	E ( L S )	E ( L V )
4	0.0056	0.0782	0.8937	0.2851	0.2290	1.6429	0.2500	0.2667
7	0.0057	0.0804	0.9193	0.2601	0.2356	1.5169	0.1429	0.2667
10	0.0058	0.0813	0.9287	0.2463	0.2380	1.4654	0.1000	0.2667
13	0.0058	0.0817	0.9335	0.2378	0.2392	1.4374	0.0769	0.2667
16	0.0059	0.0819	0.9363	0.2320	0.2399	1.4198	0.0625	0.2667

Table 5: Main reliability indices of production line for different repair recovery time parameter β with (λ,r,θ,α)=(3,15,6,0.35) .

β	ξ	η	A	W	P B	E ( L R )	E ( L S )	E ( L V )
5	0.0159	0.0795	0.9088	0.2692	0.2431	1.5634	0.1667	0.2667
9	0.0089	0.0798	0.9124	0.2676	0.2370	1.5508	0.1667	0.2667
13	0.0061	0.0799	0.9137	0.2668	0.2346	1.5460	0.1667	0.2667
17	0.0047	0.0800	0.9143	0.2664	0.2333	1.5434	0.1667	0.2667
21	0.0038	0.0800	0.9147	0.2661	0.2325	1.5418	0.1667	0.2667

By means of analysis results derived in Section 4, Table 1 shows that the effects of maximum part-time job number H on main reliability indices of production line for the set of parameters (λ,r,θ,α,β)=(3,15,6,0.35,14) . Here, we assume that H=0,1,3,5,∞ ; that is, suppose that the repairman takes no part-time job, single part-time job, 3 part-time jobs, 5 part-time jobs, and multiple part-time jobs, respectively. It is observed that the stop production frequency W , mean renewal cycle E(_LR ) , and mean part-time job period E(_LV ) all increase monotonously, and the availability A , busy probability _PB , and repair recovery frequency η all decrease monotonously as the value H increases. The same conclusion holds for decreasing value of part-time job rate r in Table 2.

The effects of virus arrival rate λ on main reliability indices are reported in Table 3, where we set (r,θ,α,β)=(15,6,0.35,14) . As is expected, A , E(_LR ) , and E(_LV ) decrease, and W , _PB , ξ , and η increase as λ increases. But the mean startup period E(_LS ) does not change. The effect of varying startup parameter θ is shown in Table 4 for the set of parameters (λ,r,α,β)=(3,15,0.35,14) . We see that the mean startup period increases monotonously as the value θ decreases.

Table 5 reports the effects of repair recovery time parameter β on main reliability indices for the set of parameters (λ,r,θ,α)=(3,15,6,0.35) . It should be noted that the repair recovery time does not affect the mean startup period and mean part-time job period. The trends shown by Tables 1-5 are as expected.

Next, we numerically analyze the optimal part-time job policy of production line while maintaining a maximum profit. Let us define

: _xA ...1; income of the production line per unit time when one computer is operating,

: _xV ...1; income per unit time when the repairer is performing part-time jobs,

: _xW ...1; loss of the production line when the production stops for each time,

: _xI ...1; idle cost of the repairer per unit time,

: _xS ...1; startup cost of the repair station per unit time,

: _xB ...1; repair and repair recovery costs per unit time,

: _xη ...1; loss of the production line when the repair is recovered for each time;

then the average profit of production line per unit time is given by [figure omitted; refer to PDF] where A , W , E(_LR ) , E(_LI ) , E(_LS ) , E(_LB ) , E(_LV ) , and η are given by Theorems 3-4 and Corollary 5, respectively. For optimization analysis, we let

(1) _xA =1000 , _xV =2000 , _xW =400 , _xI =16 , _xS =25 , _xB =60 , and _xη =160 ;

(2) λ=0.5 , T(x)=1-^e-σx , [varphi](x)=1-^e-(2/3)σx , and σ,x>0 ;

(3) the startup time Q obeys 4-stage Erlang distribution with mean E(Q)=1 ;

(4) the repair time G follows 3-stage Erlang distribution with mean E(G)=5/6 ;

(5) the repair is interrupted according to a Poisson process with rate α=0.2 ;

(6) the repair recovery time B is deterministic with mean E(B)=1/3 .

The following two optimal part-time job policies will be determined so as to maximize the average profit Γ and maintain an availability no less than 0.85: (a) the maximum number of part-time jobs K when the time of each part-time job V follows an exponential distribution with mean E(V)=2 ; (b) the maximum value of u when the maximum number of part-time jobs H obeys a negative binomial distribution Nb(3,4/5) and the time of each part-time job V is a positive integer u .

For case (a), we set V(t)=1-^e-0.5t ,t...5;0 and H(z)=^zK ,|z|<1 ; the optimization problem is given by [figure omitted; refer to PDF] For case (b), we set [figure omitted; refer to PDF] and H(z)=^(4z/(5-z))3 , |z|<1 ; the optimization problem is [figure omitted; refer to PDF] With Matlab 7.0 and a heuristic method, the numerical results for two optimization problems are reported in Tables 6 and 7, respectively. Table 6 shows that for case (a), a maximum average profit value per unit time of 2141.3 is achieved at ^K* =6 . It is seen from Table 7 that for case (b), the maximum average profit per unit time is 1946.6 as the maximum positive integer value of u equals 2.

Table 6: Availability A and average profit per unit time Γ(K) of production line for different maximum number K of part-time jobs.

K	1	2	3	4	5	6	7	8	9
A	0.9019	0.8821	0.8692	0.8606	0.8547	0.8506	0.8478	0.8458	0.8443
Γ ( K )	1367.1	1666.2	1860.6	1990.8	2079.7	2141.3	2184.3	2214.5	2235.9

Table 7: Availability A and average profit per unit time Γ(u) of production line for different time u of each part-time job.

u	1	2	3	4	5	6	7	8	9
A	0.9115	0.8635	0.8103	0.7593	0.7127	0.6707	0.6330	0.5991	0.5684
Γ ( u )	1623.5	1946.6	2088.8	2153.7	2183.3	2195.7	2199.3	2198.3	2194.7

6. Conclusions

In this paper, with the help of supplementary variable method and Laplace transform, we successfully obtain main steady-state reliability indices of an extended shock model, such as availability, failure frequency, mean renewal cycle, mean vacation period, and mean startup period. Some special cases are given. As an application, a production line is modeled to study numerical illustration and its optimal part-time job policy. For future research, one could consider some discrete time shock models with a repairman vacation and their optimization applications.

Acknowledgments

The author would like to thank the referees and editor for their valuable comments and suggestions. This work is supported by the Basic and Frontier Research Foundation of Chongqing of China (cstc2013jcyjA00008) and the Scientific Research Starting Foundation for Doctors of Chongqing University of Technology (2012ZD48).

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

[1] J. D. Esary, A. W. Marshall, F. Proschan, "Shock models and wear processes," The Annals of Probability , vol. 1, no. 4, pp. 627-649, 1973.

[2] R. E. Barlow, F. Proschan Mathematical Theory of Reliability , John Wiley & Sons, New York, NY, USA, 1965.

[3] S. M. Ross, "Generalized Poisson shock models," The Annals of Probability , vol. 9, no. 5, pp. 896-898, 1981.

[4] E. Fagiuoli, F. Pellerey, "Preservation of certain classes of life distributions under Poisson shock models," Journal of Applied Probability , vol. 31, no. 2, pp. 458-465, 1994.

[5] S. Eryilmaz, "Life behavior of a system under discrete shock model," Discrete Dynamics in Nature and Society , vol. 2012, 2012.

[6] S. Eryilmaz, "On the lifetime behavior of a discrete time shock model," Journal of Computational and Applied Mathematics , vol. 237, no. 1, pp. 384-388, 2013.

[7] G. Gottlieb, "Optimal replacement for shock models with general failure rate," Operations Research , vol. 30, no. 1, pp. 82-92, 1982.

[8] T. Aven, S. Gaarder, "Optimal replacement in a shock model: discrete time," Journal of Applied Probability , vol. 24, no. 1, pp. 281-287, 1987.

[9] G. J. Wang, Y. L. Zhang, " δ -shock model and its optimal replacement policy," Journal of Southeast University , vol. 31, no. 5, pp. 121-124, 2001.

[10] Y. Lam, Y. L. Zhang, "A shock model for the maintenance problem of a repairable system," Computers & Operations Research , vol. 31, no. 11, pp. 1807-1820, 2004.

[11] Y.-Y. Tang, Y. Lam, "A δ -shock maintenance model for a deteriorating system," European Journal of Operational Research , vol. 168, no. 2, pp. 541-556, 2006.

[12] Z. H. Li, P. Zhao, "Reliability analysis on the δ -shock model of complex system," IEEE Transactions on Reliability , vol. 56, no. 2, pp. 340-348, 2007.

[13] Q. T. Wu, S. M. Wu, "Reliability analysis of two-unit cold standby repairable systems under Poisson shocks," Applied Mathematics and Computation , vol. 218, no. 1, pp. 171-182, 2011.

[14] J. C. Ke, "Bi-level control for batch arrival queues with an early startup and un-reliable server," Applied Mathematical Modelling , vol. 28, no. 5, pp. 469-485, 2004.

[15] J.-C. Ke, "Modified T vacation policy for an M/G/1 queueing system with an unreliable server and startup," Mathematical and Computer Modelling , vol. 41, no. 11-12, pp. 1267-1277, 2005.

[16] G. Choudhury, K. Deka, "An M/G/1 retrial queueing system with two phases of service subject to the server breakdown and repair," Performance Evaluation , vol. 65, no. 10, pp. 714-724, 2008.

[17] N. S. Tian, Z. G. Zhang Vacation Queueing Models-Theory and Applications , Springer, New York, NY, USA, 2006.

[18] D. H. Shi, "A new method for calculating the mean failure numbers of a repairable system during (0,t]," Acta Mathematicae Applicatae Sinica , vol. 8, no. 1, pp. 101-110, 1985.

[19] S. M. Ross Applied Probability Models with Optimization Applications , Dover, New York, NY, USA, 1992.