RELIABILITY ANALYSIS OF A FUEL SUPPLY SYSTEM IN AN AUTOMOBILE ENGINE

Transcription

1 International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN , Vol. 9 No. III (September, 215), pp RELIABILITY ANALYSIS OF A FUEL SUPPLY SYSTEM IN AN AUTOMOBILE ENGINE R. K. AGNIHOTRI 1, CHITRANJAN SHARMA 2, PRADEEP K. JOSHI 3 & SANJAY JAIN 4 1 Professor and Head, Department of Statistics, St. John.s College, Agra-2822, India 2 Associate Professor, Department of Mathematics, Holkar Science PG College, Indore, India 3 Associate Professor, School of Computers and Elex, IPS ACADEMY, Indore, India pradeepkjoshi1@gmail.com 4 Assistant Professor, Department of Statistics, St. John.s College, Agra-2822, India Abstract The present paper deals with the analysis of a fuel supply system in an automobile engine of a four wheeler which is having both the option of fuel i.e. PETROL and CNG. Since CNG is cheaper than petrol so the priority is given to consume CNG as compared to petrol. An automatic switch is used to start petrol supply at the time of failure of CNG supply. Using regenerative point technique with Markov renewal process, the reliability characteristics which are useful to system designers are obtained. Key Words : Reliability, Redundancy, Repair Time, Transition Probability, Regenerative Points, Markov Renewal Process. 21 AMS Subject Classification : Primary 9B25, 62N5, 68M15, Secondary : 6K1, 6K2, 6K15. c http: // 125

2 126 R. K. AGNIHOTRI, C. SHARMA, P. K. JOSHI & S. JAIN 1. Introduction Many authors working in the field of system reliability have analysed various engineering systems by using different sets of assumptions. Several authors developed and analysed many hypothetical systems, but there is a need to work on the systems which are used very commonly and related to the real practical situation. [5] provided studied a system with intermittent repair and inspection under abnormal weather and also provided cost benefit analysis of the system. [6] calculated MTSF and availability analysis of a complex system composed of two sub-systems in series subject to random shocks with single repair facility. [8] showed how to optimize the parallel systems subject to two modes of failure and repair provision. [1] studied and carried out reliability analysis of a system of boiler used in Readymade garment industry. [3] studied a reliability model of 2-out-of-3 redundant system subject to degradation after repair. [2] presented reliability analysis of a system of wheels in a four wheeler. [7] introduced reliability analysis of a natural gas pressure regularity installation. [4] presented analysis of a parallel redundant system with different failure mode. Keeping the above in view, we in the present paper analysed a system of fuel supply system in an automobile engine of a four wheeler which is having both the option of fuel i.e. PETROL and CNG. Since CNG is cheaper than petrol so the priority is given to consume CNG as compared to petrol. An automatic switch is used to start petrol supply at the time of failure of CNG supply. Using regenerative point technique with Markov renewal process, the following reliability characteristics which are useful to system designers are obtained. (1) Transition and steady state transition probabilities (2) Mean Sojourn times in various states (3) Mean time to system failure (MTSF) (4) Point wise and Steady state availability of the system (5) Expected Busy period of the repairman in (, t] (6) Expected number of visits by the repairman in (, t].

3 RELIABILITY ANALYSIS OF A FUEL SUPPLY Model Description and Assumptions (1) The fuel system of an automobile engine consists of both the facilities i.e. CNG and petrol supply. Since CNG is cheaper than petrol so the priority will be given to CNG for both use and repair, to run the four wheeler. (2) Initially the four wheeler starts operation with CNG and the petrol facility is kept as standby. Both CNG and petrol has only two modes i.e. normal and failure. (3) Whenever CNG system fails petrol starts operation instantaneously with the help of an automatic switch, provided that the switch is good at the time of need, otherwise the petrol supply will wait until the switch is repaired. Probability that the automatic switch will be in good position at the time of need is fixed. (4) Since the consumption of petrol is too costly as compared to CNG, so there is a provision of rest to petrol supply after continuous working for a random amount of time. (5) A single repair facility is available in the system to repair the jobs related to CNG, petrol supply and automatic switch. (6) The failure time distributions of CNG and petrol supply are exponential with different parameters while the repair time distributions of jobs related to CNG, petrol and automatic switch are arbitrary. Also the distribution of time after which rest will be provided to petrol and its rest time are exponential with different parameter. 3. Notation and Symbols C NS : CNG is in normal supply P NS : Petrol is in normal supply P S : Petrol supply is in standby position C fr : Failed CNG supply under repair C fwr : Failed CNG supply waiting for repair P rest : Petrol is in rest position P fr : Failed petrol supply under repair P fwr : Failed petrol supply waiting for repair AS r : Automatic switch under repair

4 128 R. K. AGNIHOTRI, C. SHARMA, P. K. JOSHI & S. JAIN α : Constant failure rate of CNG supply β : Constant failure rate of petrol supply γ : Constant rate of completing rest for petrol supply δ : Constant rate of time after which rest is to be provide to petrol supply g(.), G(.) : pdf and cdf of repair time distribution of CNG supply h(.), H(.) : pdf and cdf of repair time distribution of petrol supply k(.), K(.) : pdf and cdf of repair time distribution of automatic switch p(= 1 q) : Probability that automatic switch will be in good position at the time of need. Using the above notation and symbols the following are the possible states of the fuel system of an automobile engine system: Up States : S (C NS, P S ) S 1 (C fr, P NS ) S 5 (C NS, P rest ) S 6 (C NS, P fr ). Down States : S 2 (C fwr, P S, AS r ) S 3 (C fr, P rest ) S 4 (C fr, P fwr ). The transitions between the above states are shown in Figure 1.

5 RELIABILITY ANALYSIS OF A FUEL SUPPLY Transition Probabilities Let T (= ), T 1, T 2, be the epochs at which the system enters the states S i E. Let X n denotes the state entered at epoch T n+1 i. e. just after the transition of T n. Then {T n, X n } constitutes a Markov-renewal process with state space E and Q ik (t) = P r[x n+1 = S k, T n+1 T n t X n = S i ] (1) is semi Markov-Kernal over E. The stochastic matrix of the embedded Markov chain is P = p ik = lim t Q ik = Q( ). (2) Thus P 1 = p p 2 = q p 1 = G(β + δ) p 13 = δ β+δ G(β + δ) p 14 = p (y) 16 = β β+δ [1 G(β + δ)] p 21 = p 46 = 1 p (1) 3 = p (1) 33 = p (1,4) 36 = γ [tildeg(β δ) G(γ)] γ β δ δ (γ β δ)(β + δ) [(γ β δ) γ G(β + δ) + (β + δ) G(γ)] β (γ β δ)(β + δ) [(γ β δ) γ G(β + δ) + (β + δ) G(γ)] p 35 = G(γ) p 5 = γ α+γ p 53 = α α+γ p 6 = 1 H(α) p 64 = H(α) (3-16) 5. Mean Sojourn Times The mean time taken by the system in a particular state S i before transiting to any other state is known as mean sojourn time and is defined as µ i = P [T > t]dt T is the time of stay in state S i by the system.

6 13 R. K. AGNIHOTRI, C. SHARMA, P. K. JOSHI & S. JAIN To calculate mean sojourn time µ 1 in state S i, we assume that so long as the system is in state S i, it will not transit to any other state. Therefore; µ = 1 α µ 1 = 1 β+δ [1 G(β + δ)] µ 2 = K(t)dt = m 1 (say) µ 3 = 1 γ [1 G(γ)] µ 4 = G(t)dt = m 2 (say) µ 5 = 1 α+γ µ 6 = 1 α [1 H(α)]. Contribution to Mean Sojourn Time (17-23) For the contribution to mean sojourn time in state S i E and non-regenerative state occurs, before transiting to S j E, i. e., Therefore, m ij = t q ij (t) dt = q ij (). m 1 = p t αe αt dt m 2 = q t αe αt dt m 1 = t e (β+δ)t g(t)dt m 13 = δ t e (β+δ)t G(t)dt m 14 = t βe (β+δ)t G(t)dt m 21 = t k(t)dt m 35 = t e γt g(t)dt m 46 = t g(t)dt m 5 = t γe (α+γ)t dt m 53 = t αe (α+γ)t dt m 6 = t e αt h(t)dt m 64 = t αe αt H(t)dt m (4) 16 = m (1) 3 = m (1) 33 = m (1,4) 36 = β t (1 e (β+δ)t )dg(t) β + δ γ t (e (β+δ)t e γt )dg(t) γ β δ δ γ t (e (β+δ)t e γt )G(t)dt γ β δ β t [(γ β δ) e (β+δ)t (β + δ)e γt ]dg(t). (24 39) (β + δ)(γ β δ)

7 RELIABILITY ANALYSIS OF A FUEL SUPPLY Hence using (24-39) the following relations can be verified as follows. m 1 + m )2 = p = m 1 + m 13 + m (4) 16 = m 21 = + β β + δ t αe αt dt + q t αe αt dt = 1 α = µ t e (β+δ)t dg(t) + δ t αe αt dt t e (β+δ)t G(t)dt t (1 e (β+δ)t )dg(t) = n 1 t k(t)dt = µ 2 = m 1 (say). m 35 + m (1) 3 + m(1) 33 + m(1,4) + γ + γ β δ + δ γ γ β δ β (β + δ)(γ β δ) 36 = t e γt dg(t) t (e (β+δ)t e γt )dg(t) t (e (β+δ)t e γt )G(t)dt t [(γ β δ) e (β+δ)t (β + δ)e γt ]dg(t) = n 2 (say) m 46 = m 5 + m 53 = m 6 + m 64 = t g(t)dt = µ 4 = m 2 t γe (α+γ)t dt + t e αt dh(t) + t αe (α+γ)t dt = 1 α + γ = µ 5 t αe αt H(t)dt = 1 α [1 H(α)] = µ 6. (4 46) 6. Reliability and Mean Time to System Failure To obtain the reliability of the system we regard all the failed states as absorbing states. If T i be the time to system failure when it starts operation from state S i then the reliability of the system is given by R i (t) = P r[t i > t].

8 132 R. K. AGNIHOTRI, C. SHARMA, P. K. JOSHI & S. JAIN Using the probabilistic arguments in view of reliability the following recursive relations can be easily obtained. R (t) = Z (t) + q 1 (t) c R 1 (t) + q 2 (t) c R 2 (t) R 1 (t) = Z 1 (t) + q 1 (t) c R (t) + q 13 (t) c R 3 (t) R 2 (t) = Z 2 (t) + q 21 (t) c R 1 (t) R 3 (t) = Z 3 (t) + q (1) 3 (t) c R (t) + q (1) 33 (t) c R 3(t) + q 35 (t) c R 5 (t)r 5 (t) R 5 (t) = Z 5 (t) + q 5 (t) c R (t) + q 53 (t) c R 3 (t) R 6 (t) = Z 6 (t) + q 6 (t) c R (t) (47 52) Z (t) = e αt Z 1 (t) = e (β+δ)t G(t) Z 2 (t) = K(t) Z 3 (t) = γ γ β δ [γe (β+δ)t (β + δ)e γt ]G(t) (53 58) Z 5 (t) = e (α+γ)t Z 6 (t) = e αt H(t). Taking Laplace transform of the above equations (47-52) and solving them for R (s) by omitting the argument s for brevity, we get R (s) = N 1 (s)/d 1 (s) (59) N 1 (s) = (1 q (1) 33 q 35 q 53)(Z + Z 1q 1 + Z 1q 2q 21 +(Z 2q 2) + q 13(Z 3 + Z 5q 35)(q 1 + q 2q 21) (6) and D 1 (s) = (1 q (1) 33 q 35q 53)(1 q 1q 1 q 2q 21q 1) q 13(q (1) 3 + q 35q 5)(q 1 + q 2q 21). (61) Now the reliability of the system can be obtained by taking inverse Laplace transform of (59) for well known failure time distributions. Also the Mean time to system failure (MTSF) when the system starts from the state S can be obtained as E(T ) = lim s R (s) = N 1 ()/D 1 (). (62)

9 RELIABILITY ANALYSIS OF A FUEL SUPPLY Now, Z () = Z 1() = = e αt dt = 1 α = µ e (β+δ)t G(t)dt 1 β + δ [1 G(β + δ)] = µ 1 Z2() = Z3() = Z 5 (t) = Z 6 (t) = K(t)dt = m 1 γ γ β δ [γe (β+δ)t (β + δ)e γt )]G(t)dt = L e (α+γ)t dt = 1 α + γ = µ 5 (say) e αt H(t)dt = 1 α [1 H(α)] = µ 6. (63 68) Therefore N 1 () = (1 p (1) 33 p 35p 53 )(µ + µ 1 + m 1 p 2 ) + p 13 (L + µ 5 p 35 ) (69) and D 1 () = p 14 (1 p (1) 33 p 35p 53 ) + p 13 p (1,4) 36. (7) Hence, E(T ) = N 1 /D 1 (71) N 1 and D 1 are same as in (69) and (7). 7. Mean UP-Time of the System Let A i (t) be the probability that starting from state S i the system will remains up at epoch t. Using probabilistic arguments the following recursive relations can be easily

10 134 R. K. AGNIHOTRI, C. SHARMA, P. K. JOSHI & S. JAIN obtained A (t) = M (t) + q 1 (t) c A 1 (t) + q 2 (t) c A 2 (t) A 1 (t) = M 1 (t) + q 1 (t) c A (t) + q 13 (t) c A 3 (t) + q (4) 16 (t) c A 6(t) A 2 (t) = q 21 (t) c A 1 (t) A 3 (t) = M 3 (t) + q (1) 3 (t) c A (t) + q (1) 33 (t) c A 3(t) + q 35 (t) c A 5 (t) + q (1,4) 36 (t) c A 6 (t) A 4 (t) = q 46 (t) c A 6 (t) A 5 (t) = q 5 (t) c A (t) + q 53 (t) c A 3 (t) A 6 (t) = M 6 (t) + q 6 (t) c A (t) + q 64 (t) c A 4 (t) (72 78) M (t) = e αt M 1 (t) = e (β+δ)t G(t) M 3 (t) = γ γ β δ e (β+δ+γ)t G(t) M 6 (t) = e αt H(t). (79 82) Taking laplace transform of above equations (72-78) and solving for pointwise availability A (s), by omitting the arguments s. for brevity, one gets A (s) = N 2(s) D 2 (s) (83) N 2 (s) = (1 q 46q 64)(1 q (1) 33 q 35q 53)[M + M 1 (q 1 + q 2q 21)] + q 13(q 1 + q 2q 21) and [(1 q 46q 64)q 35M 5 + q (1,4) 36 M 6 ] + M 6 q (4) 16 (q 1 + q 2q 21)(1 q (1) 33 q 35q 53) (84) D 2 (s) = (1 q 46q 64)(1 q (1) 33 q 35q 53)[1 q 1(q 1 + q 2q 21)] Thus q 13(q 1 + q 2q 21)[(1 q 46q 64)(1 q (1) 33 q 35q 53) q (1,4) 36 (1 q 6 q 46q 64)] (q 1 + q 2q 21)(1 q (1) 33 q 35q 53)q (4) 16 q 6. (85) D 2() = p 6 (1 p (1) 33 p 35p 53 )(µ +n 1 +m 1 p 2 )+[p 13 p (1,4) 36 +p (4) 16 (1 p(1) 33 p 35p 53 )(µ 6 +m 2 p 64 ). (86)

11 RELIABILITY ANALYSIS OF A FUEL SUPPLY Also M () = M1 () = M3 () = M6 () = and then e αt dt = 1/α = µ e (β+δ)t G(t)dt = µ 1 γ γ β δ (e (β+δ)t e γt ) G(t)dt = (µ 1 µ 3 )γ (γ β δ) = L 1 (say) e αt H(t)dt = µ 6 (87 9) N 2 () = (1 p (1) 33 p 35p 53 )[p 6 (µ + µ 1 ) + p (4) 16 µ 6] + p 13 [p (1,4) 36 µ 6 + p 6 (p 35 µ 5 + L 1 )] (91) Therefore A = lim t U (t) = lim s s U (s) = lim s N 2 (s)/d 2 (s) = N 2 /D 2 (92) N 2 and D 2 are same as in equations (92) and (86) respectively. 8. Mean Down-time of the System Let B i (t) be the probability that starting from regenerative down state S j the system will remains down at epoch t. Using probabilistic arguments the following recursive relations among B i (t) ss can be easily obtained B (t) = q 1 (t) c B 1 (t) + q 2 (t) c B 2 (t) B 1 (t) = q 1 (t) c B (t) + q 13 (t) c B 3 (t) + q (4) 16 (t) c B 6(t) B 2 (t) = M 2 (t) + q 21 (t) c B 1 (t) B 3 (t) = M 3 (t) + q 3 (t) c B (t) + q (1) 33 (t) c B 3(t) + q 35 (t) c B 5 (t) + q (1,4) 36 (t) c B 6 (t) B 4 (t) = q 46 (t) c B 6 (t) B 5 (t) = q 5 (t) c B (t) + q 53 (t) c B 3 (t) B 6 (t) = q 6 (t) c B (t) + q 64 (t) c B 4 (t) (93 99) M 2 (t) = K(t) M 3 (t) = e γt G(t) (1 11)

12 136 R. K. AGNIHOTRI, C. SHARMA, P. K. JOSHI & S. JAIN Taking Laplace transform of above equations (93-99) and solving them for B (s), by omitting the arguments s for brevity, one gets B (s) = N 3(s) D 3 (s) N 3 (s) = (1 q 46q 64)(1 q (1) 33 q 35q 53)M 2 q 2 and D 3 (s) is same as D 2 (s) in equation (85). Now, Then Theerefore, (12) +q 13(q 1 + q 2q 21)(1 q 46q 64)M 3 (13) M 2 () = M 3 () = K(t)dt = µ 2 = m 1 e γt G(t)dt = µ 3. (14 15) N 3 () = p 6 [(1 p (1) 33 p 33 p 35 p 53 )p 2 m 1 = p 13 µ 3 ]. (16) B = lim t A (t) = lim s s B (s) = lim s N 3 (s)/d 3 (s) = N 3 ()/D 3() = N 3 /D 3 (17) N 3 is same as in (16) and D 3 is as D 2 in (86). 9. Expected Number of Repairs of Failed Unit Let we define, V i (t) as the expected number of repairs by the repairman in (, t] given that the system initially started from regenerative state S i at t =. Then following recurrence relations among V i (t) s can be obtained as; V (t) = Q 1 (t)$v 1 (t) + Q 2 (t)$v 2 (t) V 1 (t) = Q 1 (t)$[1 + V (t)] + Q 13 (t)$v 3 (t) + Q (4) 16 (t)$[1 + V 6(t)] V 2 (t) = Q 21 (t)$v 1 (t) V 3 (t) = Q (1) 3 (t)$[1 + V (t)] + Q (1) 33 (t)$v 3(t) + Q 35 (t)$[1 + V 5 (t)] + Q (1,4) 36 (t)$[1 + V 6 (t)] V 4 (t) = Q 46 (t)$[1 + V 6 (t)] V 5 (t) = Q 5 (t)$v (t) + Q 53 (t)$v 3 (t) V 6 (t) = Q 6 (t)$[1 + V (t)] + Q 64 (t)$v 4 (t) (18 114)

13 RELIABILITY ANALYSIS OF A FUEL SUPPLY Taking laplace stieltjes transform of the above equations (18-114) and the solution of Ṽ (s) may be expressed by omitting the argument s for brevity as Ṽ (s) = N 4 (s)/d 4 (s) (115) N 4 (s) = ( )(1 ( )(1 (1) [( (1) )( 1 + (4) 16 )( ) + (1) )( ) (1,4) 36 )( ) + (4) 16 (1,4) 36 ( )] (116) and D 4 (s) = ( )(1 (1) )[1 1 ( )] 13 ( )[( )(1 (1) (1,4) 36 )( )] ( )(1 (4) ) 16 6 (117) (1) 33 Now, N 4 () = (1 p (1) 33 p 35p 53 )(p 6 (1 p 13 ) + p (4) 16 ) + p 1p 13 [p 6 (1 p (1) 33 ) + p(1,4) 36 ]. (118) Therefore, in steady state the expected number of repairs for failed unit is V = lim t [V (t)/t] = lim s Ṽ (s) = N 4 /D 4 (119) N 4 is same as (118) and D 4 is same as D 2 in (86). 1. Expected Number of Reparis of Automatic Switch Let we define, V i (t) as the expected number of repairs for automatic switch by the repairman in (, t] given that the system initially started from regenerative state S i at

14 138 R. K. AGNIHOTRI, C. SHARMA, P. K. JOSHI & S. JAIN t =. Then following recurrence relations among V i (t) s can be obtained as; V (t) = Q 1 (t)$v 1(t) + Q 2 (t)$v 2(t) V 1(t) = Q 1 (t)$[1 + V (t)] + Q 13 (t)$v 3(t) + Q (4) 16 (t)$v 6(t)] V 2(t) = Q 21 (t)$v 1(t) V 3(t) = Q (1) 3 (t)$v (t) + Q (1) 33 (t)$v 3(t) + Q 35 (t)$v 5(t)] + Q (1,4) 36 (t)$v 6(t) V 4(t) = Q 46 (t)$v 6(t) V 5(t) = Q 5 (t)$v (t) + Q 53 (t)$v 3(t) V 6(t) = Q 6 (t)$[v (t) + Q 64 (t)$v 4(t) (12 126) Taking laplace stieltjes transform of the above equations (12-126) and the solution of Ṽ (s) may be expressed as by omitting the argument s for brevity is Ṽ (s) = N 5 (s)/d 5 (s) (127) N 5 (s) = ( )(1 (1) ) 2 21 (128) and D 5 (s) is same as in D 4 (s) in (117). Now, N 5 () = p 6 p 2 (1 p (1) 33 p 35p 53 ). (129) Therefore, in steady state the expected number of repairs of automatic switch is V = lim t [V (t)/t] = lim s Ṽ (s) = N 5 /D 5 (13) N 5 is same as (129) and D 4 is same as D 2 in (86). 11. Conclusion The paper provides the reliability characteristics such as MTSF, Availability, Busy Period Analysis, etc. by using regenerative point technique with Markov renewal process, of a fuel supply system in an automobile engine of a four wheeler which is having both the option of fuel i. e. PETROL and CNG.

15 RELIABILITY ANALYSIS OF A FUEL SUPPLY References [1] Agnihotri R. K., Khare A. and Jain S., Reliability analysis of a system of boiler used in Readymade garment industry, J. Reliab. and Stat. Studies, 1 (28), [2] Agnihotri R. K., Sharma G. and Jain S., Reliability analysis of a system of wheels in a four wheeler : Int. J. Agricult. Stat. Sci., 6(1) (21), [3] Chander S. and Singh M., Reliability modeling of 2-out-of-3 redundant system subject to degradation after repair, Journal of Reliability and Statistical Studies, 2(2) (29), [4] EL-Damcese M. A. and Tamraz N. S., Analysis of a parallel redundant system with different failure mode, JRSS, 5(1) (212), [5] Goel L. R., Sharma G. C. and Gupta P., Cost benefit analysis of a system with intermittent repair and inspection under abnormal weather, Microelectron Reliab., 25 (1985), [6] Joorel J. P. S., MTSF and availability analysis of a complex system composed of two sub-systems in series subject to random shocks with single repair facility, Microelectron Reliab., 3 (199), [7] Marko Gerbec, A reliability analysis of a natural gas pressure regularity installation, Reliab. Engg. & System Safety, 95(11) (21), [8] Reddy D. R. and Rao D. V. B., Optimisation of parallel systems subject to two modes of failure and repair provision, OPSEARCH, 3 (1993), [9] Goel L. R. and Gupta P., Stochastic analysis of a two-unit parallel system with partial and catastrophic failures and preventive maintenance, Microelectron Reliability, 24(55) (1984), 881. [1] Goel L. R. and Sharma G. C. and Gupta P., Reliability analysis of a system with preventive maintenance and inspection and two types of repair, Microelectron Reliability, 26(3) (1986), [11] Goel L. R. and Shrivastava P., Comaprison of reliability characterstics of two (double component) systems with bivariate exponential lifetimes, in Journal Systems Sci, 23(1) (1992), [12] Gopalan M. N. and Nagarwalla H. E., Cost benefit analysis of a one-server two-unbit cold stand by system with repair and preventive maintenance, Microelectron. Reliab., 25(2) (1985),