Abstract
The paper considers a production facility that might deteriorate suddenly at some point during the production run time; after deterioration, nonconforming items are produced in a greater rate compared to the rate before deterioration. Moreover, the production facility may ultimately break down; consequently, the production lot is aborted before completion. If breakdown happens, corrective action is started immediately; otherwise, the production lot is completed and preventive repair is implemented at the end of the production cycle to enhance system reliability. The mathematical model is formulated under general distributions of failure, corrective, and repair times, while the numerical examples are solved under exponential failure and uniform repair times. The formulated model successfully determines the optimal lot size in addition to the optimal process parameters (mean and standard deviation) simultaneously.
1. Introduction
In real-life situations, most production systems are unreliable to a significant degree, and process deterioration occurs at some point in the production run time. A process shift from in-control state to out-of-control state will cause the production of more nonconforming items. This process may go worse and ultimately machine breakdown happens causing production interruption, and consequently plans for meeting demand are severely affected. Building on the previous argument, the need for more realistic modeling of EMQ (Economic Manufacturing Quantity) is rising in the manufacturing field; such modeling should take into consideration system variables which affect performance in a random manner; some of those variables include the following.(1)Time to shift from in-control state to out-of control state.(2)Time to breakdown.(3)Corrective maintenance time.(4)Preventive maintenance time.(5)Rate of producing nonconforming items before and after process shift; this can be handled in various ways.(a)Assuming nonconforming items are produced only after deterioration and they are produced at constant rate.(b)Assuming nonconforming items are produced before and after a sudden deterioration; accordingly, the process mean shifts to a new value and stays unchanged. The rate of producing nonconforming items is obviously greater after deterioration, and in both cases, the rate depends on the value of the process mean.(c)Or deterioration might be considered to happen gradually following a specific pattern (e.g., linear, quadratic, or exponential).(6)Usual costs related to maintenance, production, and inventory such as corrective and preventive repair costs, inventory holding cost, shortage penalty cost, cost of producing nonconforming items, and set-up cost.
Obviously, the two factors that have a great influence on EMQ decisions are process deterioration and machine breakdown. During the last few years, researchers from all over the world have done a tremendous amount of work investigating those two issues and their effect on EMQ.
Dagpunar [1] examined the lot-sizing problem with machine time to failure following a Weibull distribution; the machine is minimally repaired until the interrupted lot is completed; at the end of the production cycle, the machine is restored to as-good-as-new condition and new cycle is started. Kuhn [2] suggested a stochastic dynamic programming model to determine the optimal lot sizing decision where the equipment is subject to stochastic breakdowns. The analysis considered two cases; first it is assumed that, after the machine breakdown, the setup is totally lost and new setup cost is incurred; the second case considers the situation in which the cost of resuming the production run after a failure might be substantially lower than the production set-up cost. Ben-Daya [3] dealt with an integrated model for the joint determination of economic production quantity and preventive maintenance (PM) level for an imperfect process having a general deterioration distribution with increasing hazard rate; the effect of PM activities on the deterioration pattern of the process is modeled using the imperfect maintenance concept. Chung [4] showed that the long-run average cost function per unit of time for the case of exponential breakdowns is unimodal but neither convex nor concave, and he obtained an approximation for lower and upper bounds on lot sizing under this condition.
Giri and Dohi [5] considered the net present value (NPV) approach to determine the economic manufacturing quantities for an unreliable production system over an infinite planning horizon. The NPV of the expected total cost was obtained under general breakdown time and general repair time. The criteria for the existence and uniqueness of the optimal production time (lot size) were derived under exponential breakdown and constant/zero repair time.
Giri and Yun [6] considered an economic manufacturing quantity problem for an unreliable manufacturing system where the machine is subject to random breakdown and at most two failures can occur in a production cycle. Upon first failure, the repair action is started immediately and the demand is met first from the on-hand inventory. The shortages, if they occur due to longer repair time, are backlogged partially by resuming the production run after machine repair. While backlogging, if failure occurs again then the accumulated shortages until completion of the second repair are assumed to be lost. The model was formulated under general breakdown and general repair time distributions.
Hsieh and Lee [7] considered two economic manufacturing quantity models with unrepairable and repairable standby key modules. They determined the economic production run length and the economic number of standbys in a deteriorating production process. Jaber [8] investigated the lot-sizing problem for reduction in setups, with reworks, and interruptions to restore process quality after process deterioration takes place. He assumed the rate of producing defectives to benefit from any changes for eliminating the defects and thus reduces with each quality restoration action.
Chiu et al. [9] studied the optimal lot-sizing decision for a production system with rework, random scrap rate and service level constraint. El-Ferik [10] studied the joint determination of both economic production quantity and preventive maintenance (PM) schedules under the realistic assumption that the production facility is subject to random breakdown and the maintenance is imperfect. The manufacturing system was assumed to deteriorate while in operation, with an increasing failure rate. The system undergoes PM either upon failure or after having reached a predetermined age, whichever of them occurs first. Chelbi et al. [11] proposed an integrated production-maintenance model for unreliable production system producing conforming and nonconforming items. The model is designed to simultaneously determine the optimal values of the lot size and the age at which preventive maintenance must be performed in order to minimize the total cost.
Chakraborty et al. [12] developed integrated production, inventory, and maintenance models for a deteriorating production system in which the production facility may shift from an βin-controlβ state to an βout-of-controlβ state and break down at any random point in time during a production run. In case of machine breakdown, production of the interrupted lot is aborted and a new production lot is started when the on-hand inventory is depleted and after corrective repair is completed. The process is inspected during each production run to examine the state of the production process.
Pentico et al. [13] investigated EOQ with partial and full backordering production systems. Hu and Zong [14] proposed an extended product inspection policy for a deteriorating production system. Based on their model, they showed that there exists a production run time and a corresponding unique inspection policy such that the expected total cost per item per cycle is minimized.
Up to our knowledge, no one has considered process targeting and production lot sizing simultaneously under process deterioration and machine breakdown conditions. In this paper we are going to study the joint effect of process deterioration, machine breakdown, and random repair times (corrective and preventive) on the optimal lot-sizing decisions in addition to process parameters (mean and standard deviation). We will consider process deterioration to happen suddenly, meaning that the process mean immediately shifts from original value to a new value .
While determining the EMQ (or production run time) is important for production planning, one can question the economic benefits that can be drawn from a model that aims to find the optimal values of the process parameters: mean and standard deviation. For a quality characteristic, the performance of the process is quantified using the specification limits, lower specification limit ( and upper specification limit . The cost of producing items within the accepted range, , is evaluated by the Taguchi loss function which is dependent on the process mean. Also the cost associated with producing items out of the specification limits is entirely affected by the process mean. Process targeting problems are found in practical situations such as filling problems where one has to control the overfilling and underfilling costs by finding the optimal process mean which lies within the specification limits. On the other hand, the process standard deviation reflects the machineβs accuracy and precision. A lower standard deviation means more expensive machines with higher accuracy levels. Hence, finding the optimal value for the standard deviation which can achieve the lowest cost is essential in minimizing the overall cost of the process.
2. Model Formulation
The model under consideration is describing a production facility which may shift from in-control-state to out-of-control state at any random time during a production run; in both states, nonconforming items are produced in different rates depending on the process mean value. The production facility may also break down at any time during a production cycle. Once a shift to the out-of-control state has occurred, it is assumed that the production process will stay in that state until the whole lot has been produced or machine breakdown occurs. If the machine breakdown happens during the production run, then the interrupted lot is aborted and a new lot is started after corrective maintenance when all available inventory is depleted (no resumptionβNRβpolicy).
Assuming that the production process starts in an βin-controlβ state at time to produce an amount , if the machine breakdown does not occur until the time , then preventive maintenance (PM) is carried out to bring the machine to the as-good-as-new condition before the start of the next production run. However, if the machine breaks down before producing the required quantity , then the corrective repair action is started immediately. During the machine repair (corrective and preventive), the demand is met first from the accumulated on-hand inventory. Shortages may occur due to longer corrective/preventive repair times. If shortages occur, they are not delivered after the machine repair.
Figures 1β4 show all possible scenarios during a production cycle. In Figure 1, machine failure happens at time ; consequently, the production lot is aborted and corrective repair is started immediately and completed before inventory is totally depleted causing no shortages. Figure 2 shows the case where machine failure happens at time , and shortages occur due to long corrective repair time. In Figures 3 and 4, no machine failure happens, accordingly preventive repair started immediately at the end of the production cycle, and again shortages may occur due to longer repair times.




2.1. Lot Sizing and Targeting Process Mean
When deterioration happens after time , process mean immediately shifts from to. The in-control and out-of-control process means are assumed to be related as in the following equation: The probability density functions for the quality characteristic before and after process shift are and , respectively, and they are assumed to be normal with known and constant standard deviation . The fraction nonconforming before and after the process shift is given by: The expected time while the process is in-control is given by The expected time while the process is out-of-control is given by The total cost of nonconforming items in a production cycle is given by We quantify the impact of imprecision on the quality characteristic by using Taguchi loss function: The total loss cost in a production cycle is given by The total expected cost per production-inventory cycle is composed of setup cost, repair costs (corrective and preventive), inventory holding cost, shortage cost [15], cost of nonconforming items, and the loss cost; it is given by And the expected length of a production-inventory cycle [15] is given by From the renewal reward theorem, the average cost per unit time is given by Our objective is to determine the optimal production run time and the optimal process mean simultaneously.
Example 1. Now we solve an example by assuming the following.
Time to shift from in-control to out-of-control state is following a uniform distribution;
In order to avoid mathematical intractability, we assume exponential time to failure with failure rate:
Assuming exponential time to failure can be justified by noticing that in the proposed model we implement either corrective or preventive maintenance after each production run, and we assume that machine is restored to as-good-as-new condition; accordingly aging effect is minimized and process can be considered as βmemory lessβ with time independent hazard rate .
Also we assume the corrective and preventive repair times to follow uniform distribution functions:
And we restrict the decision variable to assume values within the following interval;
We use the following values: , , , , , , , , , , , , , with different combinations of failure rate and deterioration factor .
From the results shown in Table 1, we notice that, under a fixed value of failure rate, the optimal production run time decreases when deterioration factor increases in order to reduce the out-of-control interval and hence reduce the number of nonconforming products, while the optimal process mean decreases and then starts to increase after a specific value of which reflects the conflict between the different terms in the cost function.
On the other hand, under a fixed value of deterioration factor , the optimal production run time increases when failure rate increases which is quite strange but can be explained by noticing that not all terms in the cost function have the same preference regarding the length of the production run time under higher failure rates. For example, from a corrective repair cost perspective, it is better to have shorter run times, while from a preventive repair cost perspective, it is better to have longer run times; and hence it is a matter of conflict resolution between the cost function terms, and accordingly the optimal run time is highly sensitive to model parameters, for example, , , and .
And we notice reasonably that the average cost per unit time increases when any of the two parameters (and) is increased.
2.2. Lot Sizing, Targeting Process Mean and Standard Deviation
In what follows we consider a model to determine the optimal production run time , the optimal process mean , and the optimal standard deviation simultaneously. Lowering the standard deviation incurs additional cost because it requires more accurate machines which are expensive.
Here we are assuming that process mean is not machine-related characteristic while the standard deviation is machine related. Moreover we assume that we have a wide range of machines available in the market varying in their accuracy (standard deviation); we need to determine the one with the optimal standard deviation which will minimize the total cost. Also we assume that all available machines are sharing the same breakdown distribution since they are all new and provided by qualified suppliers and they only differ in their accuracy. For instance, two identical CNC machines can differ in their accuracy due to different toolkits used.
The number of production cycles the new machine can serve is given by where is the expected life of the machine. The standard deviation cost in a production cycle can be expressed as Here, is the cost incurred to reduce by one unit, and is the maximum allowed value for .
The total cost in this case is given by And the cycle length is given by (9).
And again by the renewal theory, the average cost per unit time is given by
Example 2. Now we solve an example to find the optimal production cycle length, the optimal process mean, and the optimal process standard deviation simultaneously. We set, , with different combinations of failure rate and deterioration factor , and the other parameters are same as in Example 1. And we restrict the decision variables and to assume values within the following intervals: , .
From the results shown in Table 2, we notice that, under a fixed value of,(a)decreases when is increased in order to reduce nonconformities by reducing the out of control period;(b) decreases then starts to increase after a specific value of due to the conflict between the terms in the cost function;(c) increases then starts to decrease after a specific value of .
On the other hand, under a fixed value of we notice the following:(a) increases when is increased;(b) almost does not change when is increased;(c) increases when is increased.
And finally we notice that the average cost per unit time increases when any of the two parameters (and) is increased.
3. Conclusion
In this paper, we have investigated the joint effect of process deterioration and machine breakdown on production lot-sizing and process-targeting decisions. The model is formulated under general failure and corrective repair times, while numerical examples are solved under exponential failure time and uniform repair (preventive and corrective) time distributions. We have considered two models: the first one for determining the optimal production cycle length and the optimal process mean simultaneously, while in the second model we introduce a new cost resulting from lowering the process standard deviation ; accordingly we determine the optimal standard deviation in addition to and . In future research, deterioration can be considered to happen gradually following some pattern (e.g., linear or exponential) rather than being sudden.
Notations
: | Nonenegative random variable denoting time to machine breakdown |
: | Time to breakdown probability density function pdf |
: | Failure rate (parameter for) when itβs exponentially distributed |
: | Random variable denoting the time taken by the machine to shift from βin-controlβ state to βout-of-controlβ state |
: | p.d.f. of the shift time (from in-control to out-of-control state) |
: | Production cycle length |
: | Nonenegative random variable denoting corrective repair time |
: | Corrective repair time probability density function p.d.f. |
: | The upper bound on when it is uniformly distributed |
: | None-negative random variable denoting preventive repair time |
: | Preventive repair time probability density function p.d.f. |
: | The upper bound on when itβs uniformly distributed |
: | The expected time while process is in-control |
: | The expected time while process is out of control |
: | Demand rate |
: | Production rate |
: | Set up cost for each production run |
: | Corrective repair cost per unit time |
: | Preventive repair cost per unit time |
: | Inventory holding cost per unit product per unit time |
: | Shortage penalty cost per unit product |
: | Cost incurred due to production of a nonconforming item with |
: | Cost incurred due to production of a nonconforming item with |
: | Total cost of nonconforming items per cycle |
: | Taguchi loss function |
: | Total loss cost per cycle |
: | Expected total cost per production-inventory cycle |
: | Expected length of a production-inventory cycle |
: | The average cost per unit time |
: | Cost of process standard deviation |
: | Expected number of production cycles a new machine can serve |
: | Expected life of the machine |
: | The cost incurred to reduce by one unit |
: | The maximum allowed value for |
: | Random variable denoting the quality characteristic of process under consideration |
: | p.d.f. of the quality characteristic |
: | In-control process mean |
: | Out-of-control process mean |
: | Process deterioration factor |
: | Process standard deviation |
: | Lower specification limit on quality characteristic |
: | Upper specification limit on quality characteristic |
: | Proportion of nonconforming while process is in-control with |
Proportion of nonconforming while process is in-control with | |
: | Proportion of nonconforming while process is out-of-control with |
: | Proportion of nonconforming while process is out-of-control with. |