Microsoft Word - 03_MTK_ZAHEDI - ASHADI - KP-- 3.docx Integrating Preventive Maintenance … (Zahedi; Ashadi Salim) 105  INTEGRATING PREVENTIVE MAINTENANCE SCHEDULING AS PROBABILITY MACHINE FAILURE AND BATCH PRODUCTION SCHEDULING Zahedi1; Ashadi Salim2 1,2Mathematics and Statistics Department, School of Computer Science, Bina Nusantara University Jln. K.H. Syahdan No. 9 Palmerah, Jakarta Barat, 11480 1zahedizahedi@binus.ac.id; 2ashadisalim@binus.ac.id ABSTRACT This paper discusses integrated model of batch production scheduling and machine maintenance scheduling. Batch production scheduling uses minimize total actual flow time criteria and machine maintenance scheduling uses the probability of machine failure based on Weibull distribution. The model assumed no non- conforming parts in a planning horizon. The model shows an increase in the number of the batch (length of production run) up to a certain limit will minimize the total actual flow time. Meanwhile, an increase in the length of production run will implicate an increase in the number of PM. An example was given to show how the model and algorithm work. Keywords: integrated model, batch scheduling, maintenance scheduling, actual flow time, machine failure probability INTRODUCTION The independence of production in scheduling and maintenance scheduling can implicate to several conditions. First, production scheduling that does not pay attention to maintenance aspect will cause the machine continues to be used, even though the machine has to be for maintenance. If maintenance is not done, it will result in machine break down, so it will interfere production scheduling. Second, the maintenance activities do not pay attention to the production. Thus, the busy machine must be stopped for maintenance. This will interfere the production schedules that have been prepared previously. Research that integrate batch scheduling on one machine with probabilistic maintenance time determined by the optimization method as Lee and Rosenblat (1987), Wang and Sheu (2001), Tseng (1996), Ben-Daya and Noman (2006), Lin and Hou (2005), Chelbi et al. (2008), Makhdoum (1996), El-Ferric and Ben-Daya (2010) as well as Zahedi (2013). Lee and Rosenblat (1987) discussed the problem of joint control of the production cycle and maintenance scheduling as inspection of parts from the machine. This study developed a relationship between maintenance and inspection showed that the optimal inspection interval is equally-spaced. Simultaneous solution of the model is a lot EMQ and schedule inspections are completed with the cost function approach. Wang and Sheu (2001) developed an integrated model to determine the lot size and the optimal inspection policy to minimize the expectations of the average cost per unit of time in the deteriorated production system. The trade-off between production cost and maintenance cost analyzed in this model to get a production lot size and optimal time inspection. Tseng (1996) proposed an optimal policy of preventive maintenance for deteriorated production system. The proposed policy 106  ComTech Vol. 7 No. 2 June 2016: 105-112 model involves the production process that the deterioration status becomes out of control, the number of items defect, and to optimize the cost of failure includes the rework cost and warranty cost. Ben-Daya and Noman (2006) developed a model that provides integrated decision optimal inventory levels, production run length and preventive maintenance schedule and warranty decisions simultaneously. The criteria are maximizing profit per unit of product, i.e. the expectation value of revenue minus the expectation cost of manufacturing cost, the cost of PM and inspection as well as warranty cost and CM cost. Lin and Hou (2005) studied the determination of the length of the production run for the model EMQ with deteriorated production process and involved the cost of CM. The objective is to find the model of optimal intervals that minimize the total cost per unit of time consisting of setup costs, inventory cost, the cost of CM and rework cost. Chelbi et al. (2008) proposed a model of determining lot size Q and time T optimal to implement preventive maintenance based on deteriorated production system life cycle. The system criteria are minimizing the average total cost per unit of time by engaging in a machine inventory during the maintenance period. Makhdoum (1996) developed models of integrated maintenance, quantity of production, inspection schedules and design quality control policy that manipulates preventive maintenance activity to be more realistic. An average x control chart is used to monitor the production process, and to earn the effect of maintenance on the quality and quantity of product. El- ferric and Ben-Daya (2010) discussed the integration of production and maintenance, where the time required for maintenance and replacement cannot be negligible. Preventive maintenance in the case is imperfect preventive maintenance, where the products can worsen after an imperfect PM. The problem is to find the length of the production run that minimizes the total actual flow time. Zahedi (2013) discussed the integration of batch scheduling and maintenance where machine time maintenance is a proportion of the length of production run with considering a common due date. This paper proposes the integration batch production scheduling and preventive maintenance (PM) determined by the probability of machine failure. The objective of the model is to minimize the total actual flow time. METHODS Supposed a set of q parts of the same item as an order will be processed on a machine with unavailability machine as the probability of machine failure. Each part requires a process to complete the operation (single stage). The decision variables are number and schedule of batches in each production run, number and schedule of PM. The objective of the model is minimization total actual flow time. The parameters of the model are the processing time per part t, setup time between two consecutive batches s, the number of parts that will be scheduled q, and the time of delivery of the all parts d (a common due date). The machine is a deteriorating machine with Weibull distribution function with probability density function as follows. exp , where , 0 and 0, (1) The cumulative distribution function for an interval (0, p') is ′ ′ 1 exp ′ . (2) Integrating Preventive Maintenance … (Zahedi; Ashadi Salim) 107  From Eq. (2) will be derived t' as ′ ln 1 ′ . (3) The preventive maintenance actions will be done if the probability had reached 0.9. The conceptual model of batch scheduling with minimizing the total actual flow time criteria (Halim, 1993) with two productions runs as Figure (1) below. Planning horizon Production run-2 Production run-1 0 B[N2max]2 B[22] B[12 APM BPM B[N1max] B[21] B[11] d : unavailability machine interval Figure 1 Gantt-chart Single Item Single Machine with Two Production Runs in Backward Approach If there is a planning horizon with g-1unavailability intervals of the PM on the machine, then the formulation of total actual flow time can be formulated as follows. 1 (4) If there are only two production runs with single unavailability interval of the PM on the machine, then the formulation of total actual flow time is 1 , (5) is the number of parts in ith batch and kth production run in backward sequencing. Xjk is binary variable will be 1 if the batch is not empty and the value 0 if the batch is empty. tPM is the length of PM time. Some constraints of the problem can be explained as follow. First, the material balance of parts in all batches scheduled the number of parts in all batches will be equal to the total number of parts that are scheduled, with perfect process without defect parts, s L[21] s L[11] s L[22] s L[12] s L[N2max]2 s L[N1max]1 108  ComTech Vol. 7 No. 2 June 2016: 105-112 (6) Second, each batch scheduled assumed arrive just in time to be processed in the floor and must be tight to the due date, or it can be written as , 1,2, … , 1,2 (7) Third, for setting the sequencing between batches will be used variable binary that has 0 or 1 value. 1if ith batch in jth production run precedes kth batch in lth production run in backward sequencing and 0 otherwise. For all batches scheduled can be formulated as , 1,2, … , 1 , 1,2 (8) Where M is a large number to guarantee sequencing, in practice it can take as M=q(s +t). Fourth, unavailability interval of machine for the PM along tPM will be in a position as: ′ ln 1 ′ (9) (10) (11) (12) ′ (13) Fifth, the maximum amount of batch is calculated by the equation 1 , (14) then Nmax= 1 (15) Integrating Preventive Maintenance … (Zahedi; Ashadi Salim) 109  Sixth, there are necessary conditions of the negativity of decision variables and binary variables as follows 1 (16) 0 (17) q (18) 1 (19) RESULTS AND DISCUSSIONS Using those constraints and objective function, the integrated batch production and preventive maintenance scheduling as the probability of machine failure to minimize total actual flow time on a deteriorating machine can be expressed as a binary-non-linear programming as follows: Model 1 (20) Constraints (21) , 1,2, … , and 1,2 (22) , 1,2, … , 1 , 1,2 (23) ′ ln 1 ′ (24) (25) (26) (27) ′ (28) (29) 110  ComTech Vol. 7 No. 2 June 2016: 105-112 1 (30) 1 (31) 0 (32) q (33) 1 (34) To solve the model, an algorithm to make the model operated is created. Algorithm Step-1. Calculate Tmin = q.t Step-2. The problem is feasible if and only if Tmin + ≤ d. Go to Step-3. If Tmin + > d then the problem is infeasible, stop. Step-3. Calculate N(max) by Equation (29). Step-4. Calculate ′by Equation (24) Step-5. Substitute values of N with N = ⎣ Nmax⎦, q, t, s, d, , and ′ into the MODEL. Step-6. Set = 1, jika ij < kl, ∀ i, j, i≠j, and Yij = 0 otherwise. Step-7. Set F11a = n2t + 1 Step-8. For i = 1, j=1,set Xij= 1, i=1,2,…, Nj(max), j = 1,2, and Xij= 0 otherwiswe. Step-9. Set i = 1, j=1, k = 1, l=2, i=1,2,…, Nj(max), j = 1,2, k=i+1, l= 1,2,…, Nj{max). Step-10. Solve Model. Step-11. Evaluate B[kl] ≥ 0, - If B[kl] ≥ 0, Write Fkla, - Evaluate Fkla< F(k-1)ja, - if Fkla< F(k-1)ja,setk = k + 1, go to step-9. - if Fkla≥ F(k-1)ja, go to step-12. -If B[kl] <0, go to step-12. Step-12. Write all values of decision variables and objective function. The following example would be given to show how the algorithm works. Example Supposed an integrated batch production scheduling and preventive maintenance scheduling as the probability of machine failure with parameters as follows. Number of parts q = 200, setup time s Integrating Preventive Maintenance … (Zahedi; Ashadi Salim) 111  = 30, processing time t = 20, preventive maintenance time tPM = 60, delivery time or due date d = 5000, a machine with Weibull distribution function with β = 1,2 and α = 0,001. Solution Step-1 to Step-3 yield Tmin= 4000, 4000+60 < 5000, and Nmax=34, then the problem is feasible for the model. Step-4 yields ′ ′ 2300. Step-5 to Step-11 yield the best solution. The principle of the steps is an increase in the number of the batch (length of production run) up to a certain limit will minimize the total actual flow time. The best solution is given in Table 1. Table 1 The Best Solution for the Example of the Problem Qij Bij APM BPM F a Q11=34,3 B11=4313,3 472.752,8 Q21=32,8 B21=3626,7 Q31=31,3 B31=2970,0 2880,0 2940,0 Q12=17,8 B12=2534,4 Q22=15.8 B22=2188,9 Q32=14,3 B32=1873,3 Q42=12,8 B42=1587,8 Q52=11,3 B52=1332,2 Q62=9,8 B62=1106,7 Q72=8,3 B72=911,1 Q82=6,8 B82=745,6 Q92=5,3 B92=610,0 The best solution in Gantt-chart is given in Figure 2 as follows. Planning horizon Production run-2 Production run-1 ... … 0 B[92] B[22] B[12] APM BPM B[31] B[21] B[11] d=5000 610,0 2188,9 2534,4 2880 2970,0 3626,7 4313,3 : Unavailability machine Figure 2 The Best Solution for the Example of the Problem s L[21] s L[11] s L[22] s L[12] s L[92] s L[31] 112  ComTech Vol. 7 No. 2 June 2016: 105-112 CONCLUSIONS The model integrates batch scheduling and maintenance scheduling to the total actual flow time. The problem in the model is divided into two, i. e., to determine batch production schedule and the second is to determine the best schedule for preventive maintenance that minimizes total actual flow time. The solution is to accommodate a trade off in the following two things. An increase in the number of the batch (length of production run) up to a certain limit will minimize the total actual flow time. Meanwhile, an increase in the length of production run will implicate an increase in the number of PM. This research can be continued by including aspects of deterioration of the machine by including the shifting of the process from in-control state to out-of-control state that follows a function of increasing failure rate Weibull, and then calculate the number of corrective maintenance to minimize the total actual flow time. REFERENCES Ben-Daya, M. & Noman, S. A. (2006). Lot sizing, preventive maintenance and warranty decisions for imperfect production systems. Journal of Quality in Maintenance Engineering, 12(1), 68-80. Chelbi, A., Rezg, N., & Radhoui, M. (2008). Simultaneous determination of production lot size and preventive maintenance schedule forunreliable production system. Journal of Quality in Maintenance Engineering, 14(2), 161-176. El-Ferric, S., & Ben-Daya, M. (2010). Integrated Production Maintenance Model under Imperfect Age-Based Maintenance Policy and Non-Negligible Maintenace Times. Asia-Pacific Journal of Operational Research, 27(4), 539-558. Halim, A. H. (1993). Batch Scheduling for Production Systems under Just in Time Environment, Doctorate Disertation. University Osaka Perfecture, Japan. Lee, H. L., & Rosenblat, M. J. (1987). Simultaneous Determination of Production Cycle and Inspection Schedules in a Production System. Management Science, 33, 1125-1136. Lin, L. C., & Hou, K. L. (2005). EMQ Model with Maintenance Actions for DeterioratingProduction System. Information and Management Sciences, 16(1), 53-65. Makhdoum, M. A. A. (1996). Integrated production, quality and maintenance models under various preventive maintenance policies. ProQuest Dissertation and Thesis. Tseng, S.T. (1996). Optimal Preventive Maintenance Policy for Deteriorating Production Systems. IIE Transactions, 28, 687-694. Wang, C. H., & Sheu, S. H. (2001). Simultaneous Determination of the Optimal Production Inventory and Product Inspection Policies for a Deteriorating Production System. Computers & Operations Research, 28, 1093-1110. Zahedi. (2013). Model Penjadwalan Batch dengan Waktu Ketaktersediaan Mesin sebagai Proporsi Production Run dengan Kriteria Minimasi Total Actual Flow Time. Jurnal MatStat, 13(2), 112-120.