Int. J. Mathematics in Operational Research, Vol. 3, No. 4, 2011
An integrated cost model for production scheduling and perfect maintenance Laith A. Hadidi* and Umar M. Al-Turki Department of Systems Engineering, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia E-mail:
[email protected] E-mail:
[email protected] E-mail:
[email protected] *Corresponding author
M. Abdur Rahim Faculty of Business Administration, University of New Brunswick, Fredericton, NB, Canada, E3B 5A3 E-mail:
[email protected] Abstract: Production scheduling deals with scheduling production jobs on a machine (single or multiple) in order to optimise a specific objective such as total weighted completion times or total weighted tardiness. The assumption that machines are always available for processing jobs is generally used in the production scheduling literature. In reality, machines often are unavailable due to preventive maintenance activities or machine failure. Production scheduling and preventive maintenance planning are interrelated, but are most often treated separately. This interdependency seems to be overlooked in the literature. This work integrates, simultaneously, the decisions of preventive maintenance and job order sequencing for a single machine. The objective is to find the job order sequence and maintenance decisions that would minimise the expected cost. Keywords: production scheduling; maintenance; deteriorating production process; integrated models; unreliable machine. Reference to this paper should be made as follows: Hadidi, L.A., Al-Turki, U.M. and Rahim, M.A. (2011) ‘An integrated cost model for production scheduling and perfect maintenance’, Int. J. Mathematics in Operational Research, Vol. 3, No. 4, pp.395–413. Biographical notes: Laith A. Hadidi received his BSc in Mechanical Engineering and MSc in Industrial Engineering from the University of Jordan in 2002 and 2005, respectively. Currently, he is a PhD candidate in the Department of Systems Engineering at King Fahd University of Petroleum and Minerals, Dhahran, Saudi Arabia. His areas of interest include preventive maintenance planning, operations sequencing and scheduling, production systems design and simulation. Copyright © 2011 Inderscience Enterprises Ltd.
395
396
L.A. Hadidi et al. Umar M. Al-Turki is an Associate Professor in the Department of Systems Engineering at King Fahd University of Petroleum and Minerals, Dhahran, Saudi Arabia. He obtained his PhD in Decision Sciences and Engineering Systems from Rensselaer Polytechnic Institute, Troy, New York, in 1994. He obtained a MSc in Industrial Engineering and Management from Oklahoma State University in 1989. He is currently involved in doing research in the areas of scheduling and quality control and is involved in consultation in strategic planning. M. Abdur Rahim is a Professor of Quantitative Methods in the Faculty of Business Administration, University of New Brunswick, Fredericton, Canada. He received his PhD in Industrial Engineering from the University of Windsor, MSc in Systems Theory from the University of Rome, and MSc in Statistics from the University of Dhaka. His research interests are in statistical process control, inventory control, maintenance policy and total quality management. He serves on the editorial boards of several international journals. He is a Visiting Professor at King Fahd University of Petroleum and Minerals in the Department of Systems Engineering.
1 Introduction The integration between production elements has received attention in recent years. This integration is strongly justified by significant savings in operations costs; however, integrated models are not easy to solve. This difficulty is due the fact that integrated models deal with multiple conflicting objectives. In this work, production scheduling and maintenance operations, for a single machine, are integrated at the shop floor level. It is a common practice to schedule both of them independently. This independent planning is done through separate functional teams. The resulting plans of a specific function may disrupt the other function plans. For example, the maintenance function assigns scheduled shut-down intervals. These intervals will be communicated to the production unit. The suggested maintenance intervals may maximise the machine availability, but they will affect production plans. Similarly, production schedulers may have the tendency to utilise machines to their full capacity to meet demand. Under this condition, productivity may increase, but machine availability will decrease, due to having more breakdowns. This work integrates, simultaneously, the decisions of preventive maintenance and job order sequencing for a single machine. The motivation for this work is encouraged by the need of many real life applications, for example, the automotive industry (Aksoy and Ozturk, 2010). In addition, this integration is expected to provide a reduction in the total expected cost. Each job order i consists of serving Qi work pieces on a single machine. The manufacturer is expected to deliver n job orders with different processing times Pi . The machine is subjected, upon failure, to minimal repair action where each repair will cost the manufacturer cm and restores the machine to work with no improvement in its condition. The manufacturer can pick a preventive maintenance major action only before the start of serving a job order. This action will cost the manufacturer cp and will restore
An integrated cost model for production scheduling
397
the machine to an ‘as good as new’ condition reducing the chance of machine breakdown during operation. All job order work pieces are released to the shop floor at time 0, hence; the manufacturer should consider the holding cost during the scheduling horizon. The Holding Cost (HC), Minimal Repair Cost (MRC) and Preventive Maintenance Cost (PMC) compromise the expected cost. Each job order sequence will change the average expected cycle costs. The objective is to find the job order sequence and maintenance decisions that would minimise the expected cost. In the literature, some work has been done to investigate the integrated optimisation problem of preventive maintenance and production scheduling. Cassady and Kutanoglu (2005) investigated this integration, for a single machine, by conducting an extensive experimental study using small scheduling instances. The integrated objective was to minimise the total weighted expected completion times for a single machine. Similarly, Cassady and Kutanoglu (2003) investigated the benefits of integration with the objective of minimising total weighted expected tardiness for a single machine. In Yulan et al. (2008), five objectives, including minimising maintenance cost, makespan, total weighted completion time of jobs, total weighted tardiness, and maximising machine availability, are simultaneously considered to optimise the integrated problem of preventive maintenance and production scheduling introduced by Cassady and Kutanoglu (2005). In order to minimise the total tardiness, Kuo and Chang (2007) proposed a solution method to find the optimal integrated problem for a single machine under a cumulative damage process. As observed in Cassady and Kutanoglu (2003), integrated models had proved an expected reduction in cost that might reach up to a 30% savings, over treating the two problems independently. However, no analytical solutions were provided for these integrated models. Instead, a full enumeration for small size problems was suggested to solve the integrated problem (Cassady and Kutanoglu, 2003, 2005). Hence, the need for more analytical or near optimal solutions is highly appreciated. The previous models do not capture the direct costs associated with machine downtime (e.g., labour, parts, etc.) of machine failures and machine maintenance. Hence, the need for cost based models of the integrated problem is justified (Gribkovskaia et al., 2010; Khanra et al., 2010; Goswami et al., 2010; Sun-Lee and Yoon, 2010; Karamatsoukis and Kyriakidis, 2010). The use of mathematical modelling for the purpose of production scheduling or preventive maintenance planning is well established in the literature (Diaby, 2010a, 2010b). Typically, preventive maintenance planning models are stochastic models designed to either maximise equipment availability or minimise equipment maintenance costs. The focus of this work is to formulate a cost model that simultaneously considers, maintenance, minimal repair and holding costs for several production jobs with the objective to minimise the expected total costs. The maintenance is assumed to be perfect and restores the machine to an ‘as good as new’ condition. Hence, after each maintenance action a new cycle will commence without being affected by the machine condition before maintenance. The expected cycle length includes: average holding time, average maintenance time, and average minimal repair time. Similarly, total expected cycle cost includes: expected maintenance cost, expected minimal repair cost and expected holding cost. The organisation of the rest of this paper is as follows. Section 2 defines, in detail, the integrated problem and its assumptions.
398
L.A. Hadidi et al.
Sections 3 and 4 provide notations and formulation of the integrated model as a mathematical program. Section 5 presents and solves an example for the integrated problem. Section 6 derives independent maintenance model and compares it with the integrated solution found in Section 5. Finally, concluding remarks and future directions for research are shown in Section 7.
2 Problem definition Consider n job orders to be processed on a single machine that is subject to a Preventive Maintenance (PM) requirement. These job orders are available at time zero with no precedence constraints. Each job order consists of processing Qi work pieces. Moreover, the machine is available continuously along the time horizon unless a machine breakdown occurs during job processing. If so, a minimum repair is conducted that restores the machine to its condition before breakdown. Also, the job that was interrupted by the breakdown should resume the remaining portion of the job, after machine repair. Each breakdown will delay completion time of successive jobs by the needed time to repair tr (assumed to be constant). To reduce the chance for machine breakdown, a PM activity can be performed before starting any job. If so, it will restore machine condition to ‘as good as new’. This PM activity will delay successive jobs by the time of PM tp (assumed to be constant). The machine may or may not fail, causing the completion time for each job to be stochastic. The PM decisions affect the stochastic process governing machine failure. Hence, change the expected value of job completion time E(ci ). Each maintenance action costs a fixed PM cost cp . Similarly, each break down will cost a fixed minimal repair cost cm . It is assumed that all jobs are released to the ¯ i is the shop floor ready to process at the start of the schedule. Hence, given that Q average work-in-process inventory, each job order i waiting on the shop floor will ¯ i till the time when the machine completes job incur a holding cost per unit time hQ n ¯ order i, that is, E(ci ) and the total holding cost will be hE(cn ) i=1 Q i . The details of these costs will be discussed in Section 4. The problem would be to identify a set of PM decisions, as well as a set of job sequencing decisions in a way that reduces the total weighted expected completion times. Prior to the job starting, a decision has to be made whether to perform PM or not. If so, PM will take a constant time that will delay consequent jobs by such time tp . If not, the job will start at the completion time of the previous job. However, the single machine can fail during job processing. The number of failures, during job processing, is strongly affected by the machine age i.e., when the machine ages it has a higher probability to fail. Furthermore, failures are randomly distributed over machine operation. Let N (τ ) be the number of machine failures in τ time units of machine operation. Given that τPi is the time units of machine operation over Pi , then, E[N (τPi )] is the expected number of failures during machine operation for job i. For job i, the PM will reduce E[N (τPi )], but it will delay the start of job i by tp . This situation can be represented by a binary variable yi . Let 1 if PM is performed before job i yi = i = 1, 2, . . . , n. 0 Otherwise
(1)
An integrated cost model for production scheduling
399
Without loss of generality, Figure 1 represents a Gantt chart for job sequence j1 , j2 , . . . , jn−1 , jn with hatched area representing PM decisions yi tp . Figure 1 Expected cycle costs for the integrated problem (see online version for colurs)
The following assumptions are considered in the problem 1
Jobs cannot be preempted for PM
2
Jobs interrupted by failure can be resumed after repair without additional time penalty
3
Machine has increasing hazard rate
4
Upon failure, minimum repair is conducted and machine will resume with the same age
5
PM restores the machine to a ‘as good as new’ condition
6
Repair times are deterministic and known in advance
7
The number of machine breakdowns is unknown (random variable)
8
The number of breakdowns does not depend on job type
9
Machine breakdowns are independent
10
Raw material for all job orders are released at the start of the schedule.
3 Notations PM n Pj xij P[j] c[i] β η
Preventive maintenance Number of jobs to be scheduled Processing time of job j Job sequencing decision variable Processing time of the ith job in the sequence Completion time of ith job in the sequence (deterministic case) Weibull shape parameter for probability distribution of T Weibull scale parameter for probability distribution of T
400
L.A. Hadidi et al.
z(τ ) a[0] a[i] a(i−1) tr tp τ N (τ ) m(τ ) yi cp cm Qi hi
Hazard function of T Age of the machine prior to job sequencing-PM planning Age of the machine after the ith job in the sequence Age of the machine immediately prior to ith job in the sequence (after PM) Time required to repair the machine Time required to perform PM on the machine PM interval to the machine Number of machine failures in τ time units of machine operation Expected number of failures PM binary decision variable PM cost Minimal repair cost Job order size Holding cost for one working piece of job order.
4 Formulation Consider a single machine in a manufacturing system that is required to process a set of n jobs, and suppose that preempting one job for another is not permitted. The purpose of production scheduling is to choose an optimal sequence for the jobs. Let 1 if the ith job performed is job j xij = (2) 0 Otherwise P[i] processing time of ith job in the sequence then P[i] =
n
(Pj xij ) ∀ i = 1, . . . , n.
(3)
j=1
Two logical sets of constraints come out; the first set of constraints states that job i cannot seize two positions at the same time i.e., n
xij = 1 ∀ i = 1, . . . , n.
(4)
j=1
The second set of constraints states that one position cannot hold more than one job i.e., n
xij = 1 ∀ j = 1, . . . , n.
(5)
i=1
Example: Suppose a single machine has two jobs needed to be scheduled at with the following parameters: P1 = 2 and P2 = 4 then the single machine can have two possible sequences: either J1 follows J2 or J2 follows J1 .
An integrated cost model for production scheduling •
Sequence 1-2 where J1 follows J2 is expressed with following variables: x11 = 1, x12 = 0, x21 = 0, x22 = 1P[1] = 2 P[2] = 4 (see Figure 2)
•
Sequence 2-1 where J2 follows J1 is expressed with following variables: x11 = 0, x12 = 1, x21 = 1, x22 = 0 P[1] = 4 P[2] = 2 (see Figure 3).
401
Figure 2 Sequence 1-2 x11 = 1, x12 = 0, x21 = 0, x22 = 1
Figure 3 Sequence 2–1 x11 = 0, x12 = 1, x21 = 1, x22 = 0
In addition to choosing a job sequence, it is also necessary to decide whether or not to perform PM prior to each job. The integrated problem is further complicated by the fact that completion times for the jobs are stochastic, because the machine may or may not fail during each job, and PM decisions change the stochastic process governing machine failure. The completion time of a job is a random variable that depends on •
the age of the machine prior to processing the job
•
the completion time for previous jobs
•
the time to complete PM, and the PM decision
•
the job’s processing; the minimal repair time
•
the number of machine failures during the job.
PM is assumed to restore the machine an ‘as good as new’ condition i.e., machine age after PM will reduce to 0. Let a[i] be the machine age at the completion of ith job then a[i] = (1 − yi )a[i−1] + P[i] i = 1, 2, . . . , n.
402
L.A. Hadidi et al.
Suppose the machine used to process the jobs is subject to failure, and the time to failure for the machine, is governed by a Weibull probability distribution, having shaping parameter β greater than 1. When the machine fails, it is assumed that the machine is minimally repaired i.e., the machine is restored to an operating condition, but machine age is not altered. This implies that, upon machine failure, the machine operator does just enough maintenance to resume machine function. The assumption that PM restores the machine to ‘as good as new’ implies that PM is a more comprehensive action than repair that may include the replacement of many key parts in the machine. Also, the operation and maintenance of the machine (between two successive PM’s) can be modelled as a renewal process. Due to the minimal repair, the occurrence of failures during each ‘cycle’ of the renewal process can be modelled using a non homogeneous Poisson process. Given that the job i starts at age (1 − yi )a[i−1] and ends at age a[i] then a[i] z(t)dt E[N (τP[i] )] = (1−yi )a[i−1] a[i]
=
(1−yi )a[i−1]
β β−1 t dt ηβ
= m(a[i] ) − m((1 − yi )a[i−1] ) (1 − yi )a[i−1] β a[i] β − = η η
(6)
where z(t) corresponds to the hazard function to the underlying Weibull probability β τ τ distribution and m(τ ) = 0 z(t)dt = 0 ηββ tβ−1 dt = τη . Figures 4 and 5 present the effect of maintenance decision yi on the E[N (τP[i] )], where it will be a[i] reduced from (1−y z(t)dt = m(a[i] ) − m((1 − yi )a[i−1] ) (if yi = 0) down to i )a[i−1] a[i] z(t)dt = m(P[i] ) (if yi = 1). However, the PM decision, yi = 1, will delay 0 completion time of job i by tp time units. The objective is to provide jobs sequence, as well as PM schedules to minimise the expected cost E(Cost). Figure 4 Job i with PM, yi = 1 (see online version for colours)
403
An integrated cost model for production scheduling Figure 5 Job i with PM, yi = 0
The expected cycle cost includes: expected Holding Cost (HC), expected Minimal Repair Cost (M RC), and PM cost (P M C). Each job order will represent processing one batch of size Qi . Raw materials for all job orders are assumed to be released to the shop floor at time 0. Figure 6 shows the inventory level for job order i. The average inventory for a job order i is a sequence dependent. Hence Q[i] =
n
(Qj xij ) ∀ i = 1, . . . , n.
(7)
j=1
The average holding quantity is equal to the sum of areas I and II (in Figure 6) over ¯ [i] = I+II the E(C[i] ). i.e., Q E(c[i] ) ¯ [i] = Q
Q[i] E(c[i] ) −
= Q[i]
[P[i] +tr [m(a[i] )−m((1−yk )a[i−1] )]] 2
E(c[i] ) [P[i] + tr [m(a[i] ) − m((1 − yi )a[i−1] )]] 1− 2E(c[i] )
.
(8)
Figure 7 shows the expected completion times for job 1, 2, . . . , n. The expected completion times for (Job i = 1, 2, . . . , n) can be found by substituting in a recursive manner as follows
E[c[1] ] = tp y1 + P[1] + tr m(a[1] ) − m((1 − y1 )a[0] )
E[c[2] ] = E[c2 ] + tp y2 + P[2] + tr m(a[2] ) − m((1 − y2 )a[1] ) .. .
(9) E[c[n] ] = E[c[n−1] ] + tp yn + P[n] + tr m(a[n] ) − m((1 − yn )a[n−1] ) .
404
L.A. Hadidi et al.
A general formula for E[c[i] ] can be found as follows E[c[i] ] =
i
tp yk + P[k] + tr E[N (τP[k] )] i = 1, 2, . . . , n
k=1
=
i
tp yk + Pk + tr m(a[k] ) − m((1 − yk )a[k−1] ) i = 1, 2, . . . , n.
k=1
(10) Figure 6 Inventory level in the shop floor for job order i
Figure 7 Expected completion times for job order i (see online version for colours)
Hence, ¯ [i] = Q
Q[i] E(c[i] ) −
[P[i] +tr [m(a[i] )−m((1−yi )a[i−1] )]] 2
E(c[i] )
= Q[i]
1−
2
i k=1
[P[i] + tr [m(a[i] ) − m((1 − yi )a[i−1] )]] (tp yk + P[k] + tr (m(a[k] ) − m((1 − yk )a[k−1] ))
.(11)
405
An integrated cost model for production scheduling
The holding cost is usually estimated by holding cost unit h, which is expressed in monetary unit per workpiece unit per time unit. Let E(c[n] ) be the expected time to complete all n jobs. Then, the expected holding cost will be HC = h × E(c[n] ) ×
n
¯ [i] . Q
(12)
i=1
Given that cm is the expected minimal repair cost per failure. Then expected minimal repair will be MRC = cm
n
[m(a[i] ) − m((1 − yi )a[i−1] )].
(13)
i=1
Given that cp is the cost for each PM. Then expected PM cost PMC = cp
n
(14)
yi .
i=1
The expected total cost will be E(Cost) = HC + MRC + PMC n n ¯ [i] + cm = hE(c[n] ) [m(a[i] ) − m((1 − yi )a[i−1] )] Q i=1
+ cp
n
i=1
(15)
yi .
i=1
Finally, the resulting mathematical program for the integrated problem is given as follows n n n ¯ Min hE(c[n] ) [m(a[i] ) − m((1 − yi )a[i−1] )] + cp yi Q[i] + cm i=1
i=1
Subject to P[i] = Q[i] =
n j=1 n
(Pj xij ) i = 1, 2, . . . , n (Qj xij ) i = 1, . . . , n
j=1
a[0] = Initial machine age a[i] = (1 − yi )a[i−1] + P[i] i = 1, 2, . . . , n n
E(c[n] ) = tp yk + P[k] + tr m(a[k] ) − m((1 − yk )a[k−1] ) k=1
i=1
406
L.A. Hadidi et al. ¯ [i] = Q[i] Q
1−
[P[i] + tr [m(a[i] ) − m((1 − yi )a[i−1] )]]
i 2 k=1 (tp yk + P[k] + tr (m(a[k] ) − m((1 − yk )a[k−1] )) i = 1, 2, . . . , n
n j=1 n
xij = 1 i = 1, 2, . . . , n xij = 1 i = 1, 2, . . . , n
i=1
xij yi
binary i = 1, 2, . . . , n j = 1, 2, . . . , n binary i = 1, 2, . . . , n.
5 Solving the integrated problem The formulation shown in Section 4 can be solved through one of the mathematical programming languages. In the following example GAMS language was used to input the model and the BARON solver was used to reach the optimal solution. The BARON solver is a computational system designed for solving non-convex NLP optimisation problems to global optimality. When β > 1, it may be practical to perform preventive maintenance on the machine in order to reduce the increasing risk of machine failure. Table 1 considers processing 3 job orders consisting of Q1 = Q2 = Q3 = 500 work piece. Job order 1 needs 6 min for each work piece. Job order 2 needs 3 min for each work piece. Job order 3 needs 2 min for each work piece. The machine age is a[0] = 88 h. The preventive maintenance time is tp = 5 h. The machine failure rate follows a Weibull distribution with the following parameters β = 2, η = 100. Upon failure, a minimal repair is conducted with a repair time tr = 15 h. For the β t Weibull distribution m(t) = 0 1 ηββ t(β−1) dt = tη1 . h = 1.50 $/work piece/hour, cm = $1500 and cp = $500. Table 1 Example parameters
Job order
Job order size (work piece)
Processing time per work piece (minute)
Processing time per job order (hour)
500 500 500
6 3 2
50 25 16.66
1 2 3
Min hE(c[3] )
3 i=1
¯ [i] + cm Q
3 i=1
Subject to P[1] = 53x11 + 26x12 + 18.67x13
[m(a[i] ) − m((1 − yi )a[i−1] )] + cp
3 i=1
yi
An integrated cost model for production scheduling
407
P[2] = 53x21 + 26x22 + 18.67x23 P[3] = 53x31 + 26x32 + 18.67x33 Q[1] = 500x11 + 500x12 + 500x13 Q[2] = 500x21 + 500x22 + 500x23 Q[3] = 500x31 + 500x32 + 500x33 a[0] = 88 a[1] = P[1] + (1 − y1 )88
a[2] = P[2] + (1 − y2 ) P[1] + (1 − y1 )88
a[3] = P[3] + (1 − y3 ) P[2] + (1 − y2 )P[1] + (1 − y1 )88 E(c[3] ) =
3
tp yk + P[k] + tr m(a[k] ) − m((1 − yi )a[k−1] )
k=1
[P[1] + tr [m(P[1] + (1 − y1 )88) − m((1 − y1 )88)]] ¯ [1] = Q[1] 1 − Q 2(tp y1 + P[1] + tr (m(P[1] + (1 − y1 )88) − m((1 − y1 )88))) [P[2] + tr [m(a[2] ) − m((1 − y2 )a[1] )]] ¯ Q[2] = Q[2] 1 − 2 2 k=1 (tp yk + P[k] + tr (m(a[k] ) − m((1 − yk )a[k−1] )) [P + t [m(a ) − m((1 − y )a )]] r 3 [3] [3] [2] ¯ [3] = Q[3] 1 − Q 3 2 k=1 (tp yk + P[k] + tr (m(a[k] ) − m((1 − yk )a[k−1] )) x11 + x12 + x13 = 1 x21 + x22 + x23 = 1 x31 + x32 + x33 = 1 x11 + x21 + x31 = 1 x12 + x22 + x32 = 1 x13 + x23 + x33 = 1 binary variables
x11 , x12 , x13 , x21 , x22 , x23 , x31 , x32 , x33 , y1 , y2 , y3
Solution: x13 = x22 = x31 = 1 and x11 = x12 = x21 = x23 = x32 = x33 = 0 and y1 = y3 = 1 and y2 = 0 with a total expected cost that is equal to $178,030. Hence, the optimal job and PM sequences are: Job order sequence J3 − J2 − J1 and PM decisions y1 = 0, y2 = 1, y3 = 0 will give the minimal cost. This means that to get the minimum expected cost, production should start with Job 3, PM, Job 2 then Job 1. The effect of six parameters (β, tp , tr , cm , cp , h) over the expected cost is shown in the two Tables 2 and 3. A 26 = 64 factorial design was used to generate 64 trials (32 trials for β = 2 and another 32 trials for β = 3). Table 2 shows that the cost will increase with cp , cm and h increase. The need for more PM can be seen with tr = 30 h compared to tr = 15 h (trial 1 and 9 in Table 2). The machine will have more failures with the increase of β, hence, PM is more economically justified with higher machine failures. This can be shown by comparing the costs with β = 3 and β = 2 where costs are less in Table 3. The use of mathematical modelling for the purpose of production scheduling or preventive maintenance planning is well established in the literature. Typically,
408
L.A. Hadidi et al.
preventive maintenance planning models are stochastic models designed to either maximise equipment availability or minimise equipment maintenance costs. The next section will show how the integrated solution compares to the independent solution. The PM problem will be treated independently of the production sequence. The optimal PM interval is derived such that the machine availability will be maximised. The PM interval will be then superimposed on the production schedule that was found in Example 1. The resulted cost will be compared to the integrated scheduling cost.
Table 2 Parameters effect on expected cost (β = 2) Parameters Trial
β
tp
tr
cm
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
15 15 15 15 15 15 15 15 30 30 30 30 30 30 30 30 15 15 15 15 15 15 15 15 30 30 30 30 30 30 30 30
1500 1500 1500 1500 3000 3000 3000 3000 1500 1500 1500 1500 3000 3000 3000 3000 1500 1500 1500 1500 3000 3000 3000 3000 1500 1500 1500 1500 3000 3000 3000 3000
Optimal sequence cp 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000
h
PM
Production jobs
Cost
1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3
010 010 010 010 010 010 010 010 101 101 101 101 101 101 101 101 010 010 010 010 010 010 010 010 010 010 010 010 010 010 010 010
321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321
178,030 354,080 178,530 354,580 179,510 355,560 180,010 356,060 191,380 381,040 192,380 382,040 192,100 381,760 193,100 382,760 188,310 374,640 188,810 375,140 189,790 376,120 190,290 376,620 211,900 421,810 212,400 422,310 213,380 423,290 213,880 423,790
409
An integrated cost model for production scheduling Table 3 Parameters effect on expected cost (β = 3) Parameters Trial
β
tp
tr
cm
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
15 15 15 15 15 15 15 15 30 30 30 30 30 30 30 30 15 15 15 15 15 15 15 15 30 30 30 30 30 30 30 30
1500 1500 1500 1500 3000 3000 3000 3000 1500 1500 1500 1500 3000 3000 3000 3000 1500 1500 1500 1500 3000 3000 3000 3000 1500 1500 1500 1500 3000 3000 3000 3000
Optimal sequence cp 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000 500 500 1000 1000
h
PM
Production jobs
Cost
1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3 1.5 3
101 101 101 101 101 101 101 101 101 101 101 101 101 101 101 101 010 010 010 010 010 010 010 010 101 101 101 101 101 101 101 101
321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321 321
175,700 350,040 176,700 351,040 176,060 350,400 177,060 351,400 180,650 359,950 181,650 360,950 181,010 360,300 182,010 361,300 190,910 379,790 191,410 380,290 192,450 381,330 192,950 381,830 206,780 412,200 207,780 413,200 207,130 412,560 208,130 413,560
6 Solution analysis In Section 4, it was shown how maintenance decision yi can affect the E[cost]. This section will try to show and investigate the benefits of integrating the job scheduling and PM decisions. This can be done by solving independently for the optimal PM such that the availability of the machine is maximised and see how does it compare to the optimal integrated job sequence.
410
L.A. Hadidi et al.
Machine availability A(τ ) is defined as the ratio of machine uptime [0, τ ] to the total time that includes: uptime, repair time and PM time. uptime uptime + maintenance time + repair time τ = τ + tp + tr N (τ ) τ = β
. τ + tp + tr τη − 0
A(τ ) =
(16)
Hence, to find τ ∗ , A(τ ) is differentiated with respect to τ and equal it to zero. β r (β−1) τ + tp + tr τη − τ 1 + βt τ dA(τ ) ηβ = 2 dτ τ + tp + tr ( τη )β =0 β β τ τ τ + tp + t r =0 − τ + βtr η η β1 tp τ =η . tr (β − 1) 5 12 For example 1, τ ∗ = 100 15 = 57.7 h. This indicates that the machine should be maintained at age 57.7 h. Since the machine initial age is 88 that is higher than τ ∗ then PM is scheduled at time 0 then it will start processing job 3 at age 0. The machine will finish job 3 at age 16.66 and job 2 at age 41.66. Now, if the machine reaches the age of 57.7 h during job 1. In this case, there are two options: •
Schedule PM after job 1. The machine will finish job 1 at age 91.66 h (see Figure 8)
•
Schedule another PM before job 1. The machine will finish job 1 at age 53.75 h (see Figure 9)
Figure 8 1-0-0 PM, for production sequence J3 − J2 − J1
An integrated cost model for production scheduling
411
Figure 9 1-0-1 PM, for production sequence J3 − J2 − J1
For the given sequence (J3 − J2 − J1 ) there are 23 = 8 PM decisions that are shown in the following table. In Table 4, the cost for the schedule in Figure 8 is $182, 480 that is |182,480−178,030| = 2.50% higher than the optimal integrated solution found. The 178,030 cost for the schedule in Figure 9 is $181, 250 that is |181,250−178,030| = 1.8% higher 178,030 than the optimal integrated solution found. Table 4 PM effect on expected cost for job order sequence J3 − J2 − J1 , in Example 1
1 2 3 4 5 6 7 8
y1
y2
y3
Cost
0 0 0 1 1 1 1 0
0 1 1 0 0 1 1 0
1 0 1 0 1 0 1 0
181,930 178,030 181,130 182,480 181,250 187,300 190,470 202,220
Independent planning may provide optimal performance at the level of a specific function. Management usually looks at the production system as a whole and separate optimal solutions may not provide optimal solution for the whole system. Usually, there is a global optimal that includes all major functions in the production system. This global optimal can only be achieved by integrating models for all different functions. Integrated production models are expected to deal with multiple objectives with a conflicting nature. Hence, planning these elements independently will cause conflicts between functions. This disturbance can be avoided through integrated modelling.
412
L.A. Hadidi et al.
7 Conclusions In this work, production scheduling and maintenance operations, for a single machine, were integrated at the shop floor level. The objective was to find the job order sequence and maintenance decisions that would minimise the expected cost. The contribution of this work was to study the effect of production and maintenance scheduling on WIP inventory. The problem was formulated as a mixed integer program. The integrated solution provided, simultaneously, the production and maintenance schedule. Integrated modelling is expected to provide better savings over independent models. Hence, integrated modelling has recently gained momentum. However, integrated models are sophisticated and not easy to solve. Research in integrated modelling still has great potential to contribute; justified by the expected savings provided by integrated modelling. In this work, the repair and maintenance times are assumed to be constant. As an extension to this work, this assumption can be generalised where repair and preventive maintenance times are random variables following a general distribution. The model can be studied under the effect of different repair and preventive maintenance time random distributions. Another extension also can be suggested; it was assumed that the failure hazard rate is not dependent on the production jobs. However, in many cases the machine is drastically affected, by the type of job it performs. This effect can be investigated in future work.
Acknowledgements The authors would like to acknowledge King Fahd University of Petroleum and Minerals for facilitating this project. Financial assistance from the Natural Science and Engineering Research Council of Canada is greatly acknowledged. The comments and insights of the anonymous referees which improved the content of this paper are highly appreciated. The assistance of Kim Wilson for proof reading and editing the manuscript is greatly appreciated.
References Aksoy, A. and Ozturk, N. (2010) ‘Simulated annealing approach in scheduling of virtual cellular manufacturing in the automotive industry’, International Journal of Vehicle Design, Vol. 52, Nos. 1–4, pp.82–95. Cassady, C. and Kutanoglu, E. (2003) ‘Minimizing job tardiness using integrated preventive maintenance planning and production scheduling’, IIE Transactions, Vol. 35, No. 2, pp.503–513. Cassady, C. and Kutanoglu, E. (2005) ‘Integrating preventive maintenance planning and production scheduling for a single machine’, IEEE Transactions on Reliability, Vol. 54, No. 2, pp.304–309. Diaby, M. (2010a) ‘Linear programming formulation of the set partitioning problem’, International Journal of Operational Research, Vol. 8, No. 4, pp.399–427.
An integrated cost model for production scheduling
413
Diaby, M. (2010b) ‘Linear programming formulation of the vertex colouring problem’, International Journal of Mathematics in Operational Research, Vol. 2, No. 3, pp.259–289. Goswami, A., Mahata, G. and Prakash, O. (2010) ‘Optimal retailer replenishment decisions in the EPQ model for deteriorating items with two level of trade credit financing’, International Journal of Mathematics in Operational Research, Vol. 2, No. 1, pp.17–39. Gribkovskaia, I., Kovalev, S. and Werner, F. (2010) ‘Batching for work and rework processes on dedicated facilities to minimize the makespan’, Omega, Vol. 38, No. 6, pp.522–527. Karamatsoukis, C. and Kyriakidis, E. (2010) ‘Optimal maintenance of two stochastically deteriorating machines with an intermediate buffer’, European Journal of Operational Research, Vol. 207, No. 1, pp.297–308. Khanra, S., Sana, S.S. and Chaudhuri, K. (2010) ‘An EOQ model for perishable item with stock and price dependent demand rate’, International Journal of Mathematics in Operational Research, Vol. 2, No. 3, pp.320–335. Kuo, Y. and Chang, Z-A. (2007) ‘Integrated production scheduling and preventive maintenance planning for a single machine under a cumulative damage failure process’, Naval Research Logistics, Vol. 54, No. 6, pp.602–614. Sun-Lee, I. and Yoon, S. (2010) ‘Coordinated scheduling of production and delivery stages with stage-dependent inventory holding costs’, Omega, Vol. 38, No. 6, pp.509–521. Yulan, J., Zuhua, J. and Wenrui, H. (2008) ‘Multi-objective integrated optimization research on preventive maintenance planning and production scheduling for a single machine’, International Journal Advanced Manufacturing Technology, Vol. 39, Nos. 9–10, pp.954–964.