RELIABILITY AND AVAILABILITY ANALYSIS OF AN 8-STEP AUTO UNIT MANUFACTURING PROCESS HAVING A FAULT IN MAINTENANCE
Surabhi Sengar
G. B. Pant University of Agriculture and Technology, Pantnagar (India)
sursengar@gmail.com
Abstract
Human plays a vital role in the manufacturing process of a product during the planning, design, assembly, production, and maintenance phase. This paper presents a systematic method of determining the reliability and availability of an 8-step auto unit manufacturing system taking into consideration the case of occurring human error during maintenance. The whole process involves eight major units as Combustion engine, power train unit, fuel feed unit, fuel injection unit, drain or exhaust unit, engine cooling unit, brake unit, and body frame. With the failure of any of these units, the whole process can fail. Also, a constant failure rate and a general repair time are taken into the consideration for each operative unit that makes up the process. The integral differential equations were generated based on Markov modeling of the process and solution derived by considering repair time distribution using the Laplace transform. Calculations of different reliability aspects like process availability, mean time to failure steady-state nature, and profit analysis are done using supplementary variable technique and copula methodology.
Keywords: Reliability, Availability, Failure probability, Supplementary variable technique, Markov model, Human error.
1. Introduction
For a very long period, a good amount of work has been done to find the reliability of various manufacturing and industrial systems. Many researchers applied different methods like the Hidden Markov models, state space method, K map approach, probabilistic rational model method, etc. to analyze the reliability and steady-state nature of some realistic manufacturing systems like Aircraft and sea-going ships, bufferless production system, Aircraft commutators, automobile assembly process, fiber plant unit, etc. [1, 4, 5, 12, 13, and 14]. Here the discussion is focused on the reliability and availability analysis of an 8-step auto unit manufacturing process.
The automobile sector has played a significant role in the economic growth and development of any country. Because of the steady growth in the GDP of a country due to the contribution of the auto industry, it is also called an "industry of industries". Although the industry started in Germany and France now it has become one of the global all over the world.
The automobile manufacturing process is a multifaceted process involving various steps starting from designing, building, quality checks, and shipping. The first phase of the process is to design, develop, and product analysis. Once this phase is completed the next step is fulfill the requirements of tools and apparatus for the bulk production of the vehicle. This phase includes everything right from engine assembly, welding and painting to stamping machines. Following this
Surabhi Sengar RT&A, No 2 (73) RELIABILITY AND AVAILABILITY ANALYSIS OF A 8-STEP AUTO UNIT_Volume 18, June 2023
phase, is production planning which decides the quantity needs to manufacture. Before launching
the product, it has to undergo different quality and safety checks. The purpose of these checks is to
ensure that the product is designed by considering all safety standards and will perform failure-free
operation.
In this competitive scenario, to sustain itself in the market and to fulfill customer expectation it is required to increase the product performance, topographies, and quality, as a result, the complexity and mechanization of the product are increased and it has resulted in several issues related to maintenance and repair [8, 9]. When we talk about manufacturing lines, a repairman plays a vital role in this field during the design, installation, production, and maintenance. The probability of occurrence of human error is more in this maintenance phase only. Human error in repair is a concern, that has not gotten the proper attention of the researchers. It can be defined as failure to perform a particular operation that could result in a delay in the production or could damage a machine or equipment. These errors may occur due to various reasons like insufficient facilities and services, less technical skills, design errors or improper planning, etc. These errors may be classified into different categories like Design, operational, installation, quality check, maintenance, etc. Maintenance error is defined as an unintentional failure occurred during the maintenance of the product because of improper repair or proactive measures. The probability of occurrence of maintenance error increases as the product gets older [6, 10].
By considering all the above facts the present study proposes a methodology to assess the reliability and performance of a manufacturing process having a fault in maintenance due to human error. The assumed process involves eight major units as Combustion engine, power train unit, fuel feed unit, fuel injection unit, drain or exhaust unit, engine cooling unit, brake unit, and body frame. The main focus of the study is on failure during maintenance due to human error. The following assumptions have been made for the whole process:
• The combustion engines can fail due to operational errors.
• The power train unit contains: a clutch, transmission gear, differential and final drive
> It is assumed that the process can fail due to inappropriate fitting of the power train and transmission gear.
> The process can be in a degraded state because of inappropriate working of differential.
> Also, the process can completely fail due to final drive failure.
• It is considered that the process can fail due mistake in fuel injection pump timings during an inspection in the fuel feed unit.
• Furthermore, it is also assumed that the process can fail due to maintenance errors and mistakes in quality checks in the fuel injection unit because of misfiring and disturbed pressure levels.
• Improper assembly of the exhaust units can cause process failure and overheating problems in the engine cooling units can take the process into a degraded state.
• The process can fail due to design error in break unit.
• Improper installation of axle and chassis can lead to process failure in the build unit.
The whole process can fail due to the failure of any of these units. A joint repair policy is applied to repair the system in the power train when the failure occurs due to failure of the differential and final drive.
• Fuel injection unit when the failure occurs due to misfiring and disturbed pressure level. Here, the joint probability distribution is applied with the help of the Gumbel-Hougaard
Copula methodology [7, 11]. Also, failures follow exponential time distribution while general time distribution is applied for repairs. To help the production industry some parametric investigations for process reliability, availability, mean time to failure, and profit analysis has been made [2, 3]. The state transition figure is shown in figure-1.
2. Notations
Different notations used in the model are given in table 1 and the state specification of the model is described in table 2.
Table 1: Notations
Indicates the probability that initially the whole system is fully operational, denoted by state S0
The probability that the process is in the failed state due to failure of the ith unit at any time, here, i=1,2,4,5,6,7,9,10,11 and Elapsed repair time j= x, y, v, w, m, r, n, v, k
Combustion engine failure rate due to operational error
Failure rates of transmission gear and fuel injection pump in fuel supply Differential unit failure rate and final drive failure rate
The failure rate of misfiring and disturbed pressure level to be low pressure of fuel injection unit Exhaust unit failure rate
Failure rate due to overheating
The cooling unit failure rate Break unit failure rate Build unit failure rate
Shows the repair rate of kth failure in the time interval (j, j+A), where, k = EO, TI, D, DS, FSS, M, LP, EA, OVH, CS, BR, BU, and j = x, y, v, w, m, r, n, v, k Profit cost and service cost per unit of time respectively
u = em, U2 = +M(m) andX1 = **, X2 = </>D+DS(w), then joint probability is given by the expression,
$lp +m (m) = exp[ m° + (log^P +m (m)f )l'd], fo+ds (w) = exp№ + (logfe+ds (w))^)1^] using Gumbel- Hougaard copula methodology.
3. State specification
The state specification of the process is given by the table 2.
Table 2: Process state specification
States Description System
State
50 State in which the process is fully operational G
51 State in which the process gets failed due to operational error in combustion FR engine
52 State in which the process is failed due to inappropriate fitting of the power FR train and transmission gear
Po(t) :
P(j,t) :
aeo :
aTI /aFss :
ad / ads : aM ! a LP :
aea :
aovh :
acs :
abr :
aBU :
Pk(j) : K1, K2 :
Consider that,
53 State in which the process is in a degraded state because of inappropriate D working of differential
54 State in which the process is completely failed due to final drive failure FR
55 State in which the process fails due mistakes in fuel injection pump timings FR during an inspection in the fuel feed unit
56 State in which the process fails due to maintenance error and mistake in FR quality check in fuel injection unit because of misfiring and disturbed pressure level
57 State in which process gets fail due to improper assembly of the exhaust unit FR
58 State in which the process gets degraded because of overheating problem in D the engine cooling unit
59 State in which process is failed due to failure of the cooling unit FR
510 State in which process gets failed due to design error in break unit FR
511 State in which the process gets failed due to improper installation of axle and FR chassis in the build unit
Note: G= Good state; FR= Failed state under repair; D = Degraded state
K, (k)
Si
P11 (k, t)
p.. w
r
P10 (v, t)
aEO
- N 'eo aEA
S7 1 r pEA
P7 (r, t)
pFSS(m)
a
BU
S
10
am
S
1
Figure 1: State transition diagram for 8-step auto manufacturing process
4. Mathematical formulation of the model
The differential transition state probabilities for different states of the model are given by the following equations
d + a +a +a +a + (a +a \+a +a +a +a P(t) =
dt EO TI D FSS M LP) EA OVH BR BU 0
u u u
J(eo (x)P (x, t)dx J ( (y)P (y, t)dy J$D+DS (w)P (w, t)dw
+
0
+
+
J (v)P (v, t)dv + J (m)P (m, t)dm + J ^ (r)P (r, t)dr +
u u
J(s (n)P (n, t)dn + J(r (l)P0 (l, t)dl + J$BU (k)Pi (k, t)dk
0
8 8 , , N
--1---+ (b ( x)
dt dx EO
0
8 8 , , x
— + — + ( ( y)
8t 8y TI
d
--+ a
dt DS
8 8
--1---+ ( (w)
8t 8w D+DS
8 8 ■ — + — +( (v)
8t 8v FSS
8 8 , . . --1---+ (b (m)
8t 8m LP+M
" 8 8 " --1---+( (r )
8t 8r EA
d
~z + acs dt
8 8 , . . --1---+( (n)
8t 8n CS
8 8 , " — + — + ( (v)
8t 8v BR
8 + ^ + ( (k)
8t 8k BU
P1( x, t ) = 0 P2( y, t) = 0
P3(t) = aD P0(t)
P4(w, t ) = 0 P5(v, t) = 0
P6(m, t ) = 0 P7(r, t) = 0 P (t ) = aovhP0(t) P9(n, t ) = 0
Pl0(v, t ) = 0 P\\(k, t ) = 0
Boundary Conditions:
P (0, t ) = akPQ (t ) P (0,t) = a P3(t)
4 V ' DS 3 y y
P (0, t) = a P (t)
C \ / Q V y
9 v y CS 8
Where, ¡=1,2,5,6,7,10,11 andk=EO, TI, D, FSS, M, LP, EA, OVH, BR, BU
(1) (2)
(3)
(4)
(5)
(6)
(7)
(8) (9)
(10) (11) (12)
(13)
(14)
0
0
0
0
Initial Condition:
Po(0) 1 , otherwise zero. (16)
For finding the solution of the mathematical model, Taking Laplace transforms of equation (1) to (12), subject to the initial condition (16), and then solving them, we get the following transition state probabilities of the model:
Po(s) = 1
B(s)
— ak x D (s)
P ( s ) = -i-
) B ( s )
where i=1,2, 5, 6, 7,10,11 and k= EO, TI, FSS, LP+M, EA, BR, BU
P3(s) = x- 1
[s + ads ] B(ss)
= ■ adads „ dpd+ds (s)
P8(s) =
[s + ads ] b(s) aovh ^ 1
[s + aCs ] B(s)
aCs xaOVH x Dpcs (s)
P9(s) = - PCS
(17) (1 8)
(19)
(20) (21) (22)
(s + acs) x B(s)
where,
b(s) = s + aeo + ati +ad + afss + (am + alp ) + aea + aovh + abr + abu
- -aTI x Spn (s) - adads x SpD+os (s) - -
aeo x SPEO (s) [s + ads ] afss x SPfss (s)
_ _ acsaovh x S (s) _
-alp+m x SPlr+M (s) -aea x SPea (s) (s + acs ) CS -abr x SPBR (s)
abu x sPbu (s) plp +m (m) = exP[m9 + (1ogplp + m (m))9)l'9]
pd+ds o) = exp[w11 + (\ogplp + m (ww )l'1]
Also, up-state and do wn-state probabilities of the system are given by: Pup(s) = Po (s) + P3 (s) + P8 (s)
(23)
(24)
(25)
1
-x
B(s)
2 ! aD } aOVH
[s + aDs ] [s + aCs ]
(26)
Pdown(s) = P1 (s) + P2 (s) + P4 (s) + P5 (s) + P6 (s) + Pi (s) + P9 (s) + P10 (s) + P11 (s) 1 x [
B(s)
aeo xdPeo (s) + ati x DPT1 (s) + + a ] dPd+dS (s) + afss xdPfss (s)
[s + ads ]
+ aLP +m x Dplp+M (s) + aEA x DpM (s) + a<CS XaomDpCRSG +
(s + aCS )
abr x d0br (s) + abu x d'bu (s)
4.1. Steady-state behavior of the system
(27)
By using Abel's lemma,
lim{sF (s)} = lim F (t) We get following up and down time independent operational probabilities:
Pup =
1
B(0)
1
I \ aD ! aOVH aDS aCS
ad ad x ads
Pdown = -— x [ aeo x mpeo + atix mpn (s) +- —;-r"
B(0) ads + [ads ]
x M„
+ afss x m'pss + alp + m x mp„+u + aea x m0a +
(28)
aovh
a
where,
Mi = lim • s^0
Sp(s)=
+ aC0HMpCS +abr xMpBR +abu xMpBU ] acs
lim
B (0) = s^0B(s)
[1 - Si (s) 1
I s I i = peo ,pti,pd + ds ,pfss,plp +m ,pea ,pcs ,pbr ,pbu
P,
7 S + p , i = peo ,pti ,pd+ds ,pfss ,plp +m ,pea ,pcs ,pbr ,pbu
4.2. Special case
when repair rates follow exponential time distribution then,
S Pi (s) =
P.
Let,
s + P
, Vi
SPp+m (s) = ■
exp
m6 +
s + exp
[logpLP + M (m)f\
/6
m6 +
[logpLP + M (m)f\
/6
,(s) =
exp
w6 +
[logpD + DS (w)]l]L
1
Pup(s) =
s + exp 1
w6 +
[logpD + DS Mpf^
in equations (17) to (22) we get,
B1(s)
^ ! aD j _ aOVH
[s + aDs ] [s + aCs ]
- . . 1 r aEO aTI aD aDaDS 1 Pdown(s) = ——x [-—+-1J— +-D— + - D DS
B1(s) [s + pEO] [s + pTI] [s + pDS] [s + pDS][s + pD + DS] afss alp + m aea aovh acsaovh
[s + pfss] [s + plp + m] [s + pea] [s + acs] [s + pcs][s + acs]
| abr { abu
[s + pbr] [s + pbu]
cs
(29)
(30)
(31)
(32)
and
x
x
]
Surabhi Sengar RT&A, No 2 (73) RELIABILITY AND AVAILABILITY ANALYSIS OF A 8-STEP AUTO UNIT_Volume 18, June 2023
where,
aEO x PeO
Bl \s) = s + aEO + aTI +aD + aFSS + (aM + aLP ) + aEA + aOVH + aBR + aBU - "T °
[s + pEO ]
ati x &ti ad x ads pd + ds afss x pfss alp +m xplp + m aea x pea [s + pi ] [S + ads ] + pd+ds ] + pfss ] + pp +m ] + pea ]
acs x aovh x pcs abr x pbr abu xpbu
[s + pcs][s + acs] [s + pbr] [s + pbu] (33)
5. Numerical Computation
Considering different values of failure rates like aEO = 0.006, a = 0.009, ad = 0.008, aDS = 0.01,
aFSS = 0.07, aM = 0.008, aLP =0.006, aEA =0.002, aOVH =0.005, a<CS =0.003, abr =0.009, aBU =0.004, ®i = 1 for i= EO, TI, D, DS, FSS, M, LP, EA, OVH, CS, BR, BU, 9 = ^ = 1, and x= y= w= v= m= r= n= v=k = 1. Consider that all repair follows exponential time distribution then by substituting these values in equation (26) and taking inverse Laplace transformation, we have
Pup (t) = -0.002992831292 e^003000000000 ^ +0.09089745085 e(-1 10098171 -0.1045972416 1q(-7) e(-9969905879t) -0.04625367733e(-0 2043987365 t) -0.009328505569e(-000471367337 t) +0.9676775738 (34)
By, putting t=0,1, 2...,10 in equation (34) we obtain the variation of availability for time shown in figure 2. Similarly, by considering different numerical values for failure rates we get the graphs of reliability and mean time to failure given in figures 3, 4 and 5.
5.1. Cost Analysis
If it is considered that the service facility is always available, then the expected profit function in the interval (0, t] is given by
t
Ep (t) = Ki \ Pup (t)dt - K21 0 F
where, K1 and K2 are the revenue and service cost per unit of time respectively, then E P (t) = K1 [0.2602109240 e(-0.5000000000 t) + 0.003420894509 e(-0.9100000000 t) +0.3707166241 e(-1.577110005 t) +0.2282704474 e(-0.8601482323 t) cos(0.007394121216 t) +0.2901464146 e(-0.8601482323 t) sin(0.007394121216 t) -0.05867072106 e(-0.7862665523 t) -0.2283382234 e(-0.7285677930 t)-0.003350319082 e(-0.6042513821 t)-9.234941008 e(-0.1268411365 t) +8.996326343] -K2t (35)
Keeping K1 = 1 and varying K2 at 0.1, 0.2, 0.3, 0.4, 0.5 in equation (35), one can obtain Figure 6.
6. Results & Conclusion
For the more concrete behavior analysis of the process, a numerical calculation of availability, reliability, mean time to failure for various failure rates and cost function have been made. The following conclusions may be drawn based on the study conducted in the present paper.
Figure 2 shows a rapid decrease in the availability of the system for time initially because of an error in maintenance but later on, it becomes stable. Although, the reliability of the process shown in figure 3 has a constant decline and stabilizes at 0.4. Figures 4 and 5 show the variations in the meantime to failure (MTTF) for different failure rates. As the failure rates increase, the MTTF decreases and also it is observed that it is highest for fuel supply unit failure shown in figure 4 and 5. Now, cost analysis reveals that an increase in service cost results in decreased profit as shown in figure 6.
Figure 2: Time vs. Availability
time
Figure 3: Time vs. Reliability
33
25 H-1-1-1-1-1-1-1-1-1-1
0.001 0.002 0.003 0.004 0.005 0.006 0.007 O.OOB 0.009 0.01
Figure 4: MTTF vs. aEO, alp, am , abr
Figure 5: MTTF vs.
ads aFSS
Figure 6: Time vs. cost
References
[1] Bao, Y., Mays, L. W. (1990). Model for water distribution system reliability. Journal of Hydraulic Engineering, 116, 1119-1137.
[2] Gnedenko B. and Igor A.U. Probabilistic Reliability Engineering,1995.
[3] Gnedenko B, Igor V.P., Igor A. U. and Chakravarty S. Statistical Reliability Engineering, first edition,1999.
[4] Firas S. and Walid A. K. (2018). Reliability Assessment of Buffer less Production System: A Boolean Algebra and K-Map Approach, Proceedings of the International Conference on Industrial Engineering and Operations Management Washington DC, USA.
[5] Henryk T., Mariusz Z., Mariusz W. (2016). A method for reliability assessment of structural components of aircraft and sea going ships with taking into account given failure generation model, POLISH MARITIME RESEARCH, 2 (90), 23: 83-90.
[6] Ioannis S. T. and Chalikias M. (2021). On Reliability Structures with Two Common Failure Criteria Under Age-Based Maintenance Policy. Reliability: theory & applications, 2(64): 27-34.
[7] Ram, M. and Singh, S. B. (2010). Availability, MTTF, and cost analysis of the complex system under pre-emptive-repeat repair discipline using Gumbel-Hougaard family copula. International Journal of Quality & Reliability Management, 27(5): 576-595.
[8] Sengar, S. and Singh, S.B. (2012). Operational behavior and reliability measures of a viscose staple fiber plant including deliberate failures. International journal of reliability and applications,13 (1): 1-17.
[9] Sengar, S. (2020). Sensitivity analysis of a complex manufacturing system with the application of wait-in-line theory. Reliability: theory & applications, 4(59): 27-34.
[10] Sengar, S., Mangay, R., & Yigit, K. (2021). Stochastic analysis of complex redundant system having problem of waiting line in repair using copula methodology. Reliability: theory & applications, 16:383-391.
[11] Sengar S. Reliability analysis and cost optimization in manufacturing systems, System reliability Engineering, Walter de-Gruyter publication, 2021.
[12] Sengar S. and Mangey R. (2022). Reliability and Performance Analysis of a Complex Manufacturing System with Inspection facility using Copula Methodology. Reliability: Theory & Applications. 4(71): 494-508.
[13] Wawrzynski,W., Zieja, M., Tomaszewski, J., Michalski, M. (2021). Reliability Assessment of Aircraft Commutators. Energies, 14, 7404.
[14] Wasiu A. A. and Olalekan O. (2014). Reliability Analysis of a 3-Machine Power Station Using State Space Approach. International Journal of Engineering Research and Applications, 4(7): 208217.