# Production Lot Sizing and Process Targeting under Process Deterioration and Machine Breakdown Conditions

Authors: Muhammad Al-Salamah; Ali Abudari
Journal: Modelling and Simulation in Engineering
 (2012)

---

## 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.


---

## Body

## 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μI to a new value μO.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 (LSL) and upper specification limit (USL). The cost of producing items within the accepted range, LSL<x<USL, 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 timet=0 to produce an amount Q, if the machine breakdown does not occur until the time to=Q/p, 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 Q, 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.Figures1–4 show all possible scenarios during a production cycle. In Figure 1, machine failure happens at time t<t0; 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 t<t0, 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.Figure 1
Breakdown/no shortage.Figure 2
Breakdown/shortage.Figure 3
No Breakdown/no shortage.Figure 4
No Breakdown/shortage.


### 2.1. Lot Sizing and Targeting Process Mean

When deterioration happens after timeτ, process mean immediately shifts from μI toμO. The in-control and out-of-control process means are assumed to be related as in the following equation:(1)μO=γμI,γ>0.
The probability density functions for the quality characteristic X before and after process shift are q(x,μI) and q(x,μO), 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:(2)α=α1+α2=∫x=0x=LSLq(x,μI)dx+∫x=USLx=∞q(x,μI)dx,β=β1+β2=∫x=0x=LSLq(x,μO)dx+∫x=USLx=∞q(x,μO)dx.
The expected time while the process is in-control is given by(3)Tin=∫t=0to∫τ=0tτ⋅h(τ)⋅f(t)⋅dτ⋅dt+∫t=to∞∫τ=0toτ⋅h(τ)⋅f(t)⋅dτ⋅dt.
The expected time while the process is out-of-control is given by(4)Tout=∫t=0to∫τ=0t(t-τ)⋅h(τ)⋅f(t)⋅dτ⋅dt+∫t=to∞∫τ=0to(to-τ)⋅h(τ)⋅f(t)⋅dτ⋅dt.
The total cost of nonconforming items in a production cycle is given by(5)CN=p⋅Tin⋅[cLα1+cUα2]+p⋅Tout⋅[cLβ1+cUβ2].
We quantify the impact of imprecision on the quality characteristic X by using Taguchi loss function:(6)L(x)={cL(x-μI)2(LSL-μI)2,ifLSL≤x≤μI,cU(x-μI)2(USL-μI)2,ifμI≤x≤USL.
The total loss cost in a production cycle is given by(7)CLoss=p⋅Tin⋅[∫x=LSLμIL(x)⋅q(x,μI)⋅dx+∫x=μIUSLL(x)⋅q(x,μI)⋅dx]+p⋅Tout⋅[∫x=LSLμIL(x)⋅q(x,μO)⋅dx+∫x=μIUSLL(x)⋅q(x,μO)⋅dx].
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(8)Ccycle⁡=co+c1⋅∫t=0to∫l1=0∞l1⋅g1(l1)⋅f(t)⋅dl1⋅dt+c2⋅∫t=to∞∫l2=0∞l2⋅g2(l2)⋅f(t)⋅dl2⋅dt+cip(p-d)2d⋅[∫t=0tot2⋅f(t)⋅dt+to2⋅∫t=to∞f(t)⋅dt]+csd⋅∫t=0to∫l1=(p-d)t/d∞{l1-(p-d)td}⋅g1(l1)⋅f(t)⋅dl1⋅dt+csd⋅∫t=to∞∫l2=(p-d)to/d∞{l2-(p-d)tod}⋅g2(l2)⋅f(t)⋅dl2⋅dt+p⋅Tin⋅[cLα1+cUα2]+p⋅Tout⋅[cLβ1+cUβ2]+p⋅Tin⋅[∫x=LSLμIL(x)⋅q(x,μI)⋅dx+∫x=μIUSLL(x)⋅q(x,μI)⋅dx]+p⋅Tout⋅[∫x=LSLμIL(x)⋅q(x,μO)⋅dx+∫x=μIUSLL(x)⋅q(x,μO)⋅dx].
And the expected length of a production-inventory cycle [15] is given by(9)Tcycle⁡=∫t=0to∫l1=0(p-d)t/d(ptd)⋅g1(l1)⋅f(t)⋅dl1⋅dt+∫t=0to∫l1=(p-d)t/d∞(t+l1)⋅g1(l1)⋅f(t)⋅dl1⋅dt+∫t=to∞∫l2=0(p-d)to/d(ptod)⋅g2(l2)⋅f(t)⋅dl2⋅dt+∫t=to∞∫l2=(p-d)to/d∞(to+l2)⋅g2(l2)⋅f(t).⋅dl2⋅dt.
From the renewal reward theorem, the average cost per unit time is given by(10)AC(to,μI)=Ccycle⁡(to,μI)Tcycle⁡(to).
Our objective is to determine the optimal production run time to* and the optimal process meanμI* 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;(11)h(τ)={1t,t≤to,1to,t≥to.
In order to avoid mathematical intractability, we assume exponential time to failure with failure rateλ:
(12)f(t)=λe-λt.
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:(13)g1(l1)=1b1,0≤l1≤b1,g2(l2)=1b2,0≤l2≤b2.
And we restrict the decision variable to to assume values within the following interval;
(14)0≤to≤4.
We use the following values: p=130, d=100, co=300, c1=1000, c2=200,b1=3, b2=1, ci=8, cs=40, cL=30, cU=20, LSL=250, USL=260, σ=2, with different combinations of failure rate λ and deterioration factor γ.
From the results shown in Table1, we notice that, under a fixed value of failure rateλ, the optimal production run time to* 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 μI* 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 to* 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, c1, c2, ci and cs.
And we notice reasonably that the average cost per unit time increases when any of the two parameters (λandγ) is increased.Table 1
Dependent of the optimal policy (to* & μI*) on parameters λ and γ.
γto*μI*AC(to*,μI*)λ=0.041.0052.626254.4961013.531.012.612253.71511091.022.558252.721428.281.032.501254.1411730.861.042.495255.3991759.87λ=0.061.0052.654254.4961083.911.012.638253.7151178.081.022.582252.721492.011.032.522254.1411789.711.042.515255.3991818.25λ=0.081.0052.682254.4961152.191.012.665253.7151244.871.022.606252.721553.821.032.542254.1411846.791.042.536255.3991874.88


### 2.2. Lot Sizing, Targeting Process Mean and Standard Deviation

In what follows we consider a model to determine the optimal production run timeto*, the optimal process mean μI*, 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(15)N=Lmto,
where Lm is the expected life of the machine. The standard deviation cost in a production cycle can be expressed as(16)Cσ=Ke(σmax⁡-σ)N=Ktoe(σmax⁡-σ)Lm.
Here, K is the cost incurred to reduce σ by one unit, andσmax is the maximum allowed value for σ.The total cost in this case is given by(17)Ccycle⁡=co+c1⋅∫t=0to∫l1=0∞l1⋅g1(l1)⋅f(t)⋅dl1⋅dt+c2⋅∫t=to∞∫l2=0∞l2⋅g2(l2)⋅f(t)⋅dl2⋅dt+cip(p-d)2d⋅[∫t=0tot2⋅f(t)⋅dt+to2⋅∫t=to∞f(t)⋅dt]+csd⋅∫t=0to∫l1=(p-d)t/d∞{l1-(p-d)td}⋅g1(l1)⋅f(t)⋅dl1⋅dt+csd⋅∫t=to∞∫l2=(p-d)to/d∞{l2-(p-d)tod}⋅g2(l2)⋅f(t)⋅dl2⋅dt+p⋅Tin⋅[cLα1+cUα2]+p⋅Tout⋅[cLβ1+cUβ2]+p⋅Tin⋅[∫x=LSLμIL(x)⋅q(x,μI)⋅dx+∫x=μIUSLL(x)⋅q(x,μI)⋅dx]+p⋅Tout⋅[∫x=LSLμIL(x)⋅q(x,μO)⋅dx+∫x=μIUSLL(x)⋅q(x,μO)⋅dx]+Ktoe(σmax⁡-σ)Lm.
And the cycle length is given by (9).And again by the renewal theory, the average cost per unit time is given by(18)AC⁡(to,μI,σ)=Ccycle⁡(to,μI,σ)Tcycle⁡(to).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 setK=5000, Lm=500, σmax=3 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 to and σ to assume values within the following intervals: 0≤to≤4, 0≤σ≤σmax.From the results shown in Table2, we notice that, under a fixed value ofλ,(a)
to*decreases when γ is increased in order to reduce nonconformities by reducing the out of control period;(b)
μI* 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 γ.Table 2
Dependent of the optimal policy (to*,μI*,σ*) on parameters λ and γ.
γto*μI*σ*AC(to*,μI*,σ*)λ=0.041.0052.575254.2541.4641119.391.012.558253.4761.4831226.581.022.499252.5371.5811574.041.032.454254.0312.3191880.181.042.440255.3921.9411914.69λ=0.061.0052.581254.2631.4801194.871.012.563253.4861.5021300.151.022.503252.5581.6251641.171.032.464254.022.3471939.331.042.444255.3951.9611974.51λ=0.081.0052.584254.2721.4961268.161.012.566253.4961.5211371.571.022.507252.5781.6681706.181.032.473254.0082.3751996.691.042.448255.3971.9822032.53On the other hand, under a fixed value ofγ we notice the following:(a)
to* increases when λ is increased;(b)
μI* 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.


## 2.1. Lot Sizing and Targeting Process Mean

When deterioration happens after timeτ, process mean immediately shifts from μI toμO. The in-control and out-of-control process means are assumed to be related as in the following equation:(1)μO=γμI,γ>0.
The probability density functions for the quality characteristic X before and after process shift are q(x,μI) and q(x,μO), 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:(2)α=α1+α2=∫x=0x=LSLq(x,μI)dx+∫x=USLx=∞q(x,μI)dx,β=β1+β2=∫x=0x=LSLq(x,μO)dx+∫x=USLx=∞q(x,μO)dx.
The expected time while the process is in-control is given by(3)Tin=∫t=0to∫τ=0tτ⋅h(τ)⋅f(t)⋅dτ⋅dt+∫t=to∞∫τ=0toτ⋅h(τ)⋅f(t)⋅dτ⋅dt.
The expected time while the process is out-of-control is given by(4)Tout=∫t=0to∫τ=0t(t-τ)⋅h(τ)⋅f(t)⋅dτ⋅dt+∫t=to∞∫τ=0to(to-τ)⋅h(τ)⋅f(t)⋅dτ⋅dt.
The total cost of nonconforming items in a production cycle is given by(5)CN=p⋅Tin⋅[cLα1+cUα2]+p⋅Tout⋅[cLβ1+cUβ2].
We quantify the impact of imprecision on the quality characteristic X by using Taguchi loss function:(6)L(x)={cL(x-μI)2(LSL-μI)2,ifLSL≤x≤μI,cU(x-μI)2(USL-μI)2,ifμI≤x≤USL.
The total loss cost in a production cycle is given by(7)CLoss=p⋅Tin⋅[∫x=LSLμIL(x)⋅q(x,μI)⋅dx+∫x=μIUSLL(x)⋅q(x,μI)⋅dx]+p⋅Tout⋅[∫x=LSLμIL(x)⋅q(x,μO)⋅dx+∫x=μIUSLL(x)⋅q(x,μO)⋅dx].
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(8)Ccycle⁡=co+c1⋅∫t=0to∫l1=0∞l1⋅g1(l1)⋅f(t)⋅dl1⋅dt+c2⋅∫t=to∞∫l2=0∞l2⋅g2(l2)⋅f(t)⋅dl2⋅dt+cip(p-d)2d⋅[∫t=0tot2⋅f(t)⋅dt+to2⋅∫t=to∞f(t)⋅dt]+csd⋅∫t=0to∫l1=(p-d)t/d∞{l1-(p-d)td}⋅g1(l1)⋅f(t)⋅dl1⋅dt+csd⋅∫t=to∞∫l2=(p-d)to/d∞{l2-(p-d)tod}⋅g2(l2)⋅f(t)⋅dl2⋅dt+p⋅Tin⋅[cLα1+cUα2]+p⋅Tout⋅[cLβ1+cUβ2]+p⋅Tin⋅[∫x=LSLμIL(x)⋅q(x,μI)⋅dx+∫x=μIUSLL(x)⋅q(x,μI)⋅dx]+p⋅Tout⋅[∫x=LSLμIL(x)⋅q(x,μO)⋅dx+∫x=μIUSLL(x)⋅q(x,μO)⋅dx].
And the expected length of a production-inventory cycle [15] is given by(9)Tcycle⁡=∫t=0to∫l1=0(p-d)t/d(ptd)⋅g1(l1)⋅f(t)⋅dl1⋅dt+∫t=0to∫l1=(p-d)t/d∞(t+l1)⋅g1(l1)⋅f(t)⋅dl1⋅dt+∫t=to∞∫l2=0(p-d)to/d(ptod)⋅g2(l2)⋅f(t)⋅dl2⋅dt+∫t=to∞∫l2=(p-d)to/d∞(to+l2)⋅g2(l2)⋅f(t).⋅dl2⋅dt.
From the renewal reward theorem, the average cost per unit time is given by(10)AC(to,μI)=Ccycle⁡(to,μI)Tcycle⁡(to).
Our objective is to determine the optimal production run time to* and the optimal process meanμI* 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;(11)h(τ)={1t,t≤to,1to,t≥to.
In order to avoid mathematical intractability, we assume exponential time to failure with failure rateλ:
(12)f(t)=λe-λt.
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:(13)g1(l1)=1b1,0≤l1≤b1,g2(l2)=1b2,0≤l2≤b2.
And we restrict the decision variable to to assume values within the following interval;
(14)0≤to≤4.
We use the following values: p=130, d=100, co=300, c1=1000, c2=200,b1=3, b2=1, ci=8, cs=40, cL=30, cU=20, LSL=250, USL=260, σ=2, with different combinations of failure rate λ and deterioration factor γ.
From the results shown in Table1, we notice that, under a fixed value of failure rateλ, the optimal production run time to* 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 μI* 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 to* 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, c1, c2, ci and cs.
And we notice reasonably that the average cost per unit time increases when any of the two parameters (λandγ) is increased.Table 1
Dependent of the optimal policy (to* & μI*) on parameters λ and γ.
γto*μI*AC(to*,μI*)λ=0.041.0052.626254.4961013.531.012.612253.71511091.022.558252.721428.281.032.501254.1411730.861.042.495255.3991759.87λ=0.061.0052.654254.4961083.911.012.638253.7151178.081.022.582252.721492.011.032.522254.1411789.711.042.515255.3991818.25λ=0.081.0052.682254.4961152.191.012.665253.7151244.871.022.606252.721553.821.032.542254.1411846.791.042.536255.3991874.88


## 2.2. Lot Sizing, Targeting Process Mean and Standard Deviation

In what follows we consider a model to determine the optimal production run timeto*, the optimal process mean μI*, 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(15)N=Lmto,
where Lm is the expected life of the machine. The standard deviation cost in a production cycle can be expressed as(16)Cσ=Ke(σmax⁡-σ)N=Ktoe(σmax⁡-σ)Lm.
Here, K is the cost incurred to reduce σ by one unit, andσmax is the maximum allowed value for σ.The total cost in this case is given by(17)Ccycle⁡=co+c1⋅∫t=0to∫l1=0∞l1⋅g1(l1)⋅f(t)⋅dl1⋅dt+c2⋅∫t=to∞∫l2=0∞l2⋅g2(l2)⋅f(t)⋅dl2⋅dt+cip(p-d)2d⋅[∫t=0tot2⋅f(t)⋅dt+to2⋅∫t=to∞f(t)⋅dt]+csd⋅∫t=0to∫l1=(p-d)t/d∞{l1-(p-d)td}⋅g1(l1)⋅f(t)⋅dl1⋅dt+csd⋅∫t=to∞∫l2=(p-d)to/d∞{l2-(p-d)tod}⋅g2(l2)⋅f(t)⋅dl2⋅dt+p⋅Tin⋅[cLα1+cUα2]+p⋅Tout⋅[cLβ1+cUβ2]+p⋅Tin⋅[∫x=LSLμIL(x)⋅q(x,μI)⋅dx+∫x=μIUSLL(x)⋅q(x,μI)⋅dx]+p⋅Tout⋅[∫x=LSLμIL(x)⋅q(x,μO)⋅dx+∫x=μIUSLL(x)⋅q(x,μO)⋅dx]+Ktoe(σmax⁡-σ)Lm.
And the cycle length is given by (9).And again by the renewal theory, the average cost per unit time is given by(18)AC⁡(to,μI,σ)=Ccycle⁡(to,μI,σ)Tcycle⁡(to).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 setK=5000, Lm=500, σmax=3 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 to and σ to assume values within the following intervals: 0≤to≤4, 0≤σ≤σmax.From the results shown in Table2, we notice that, under a fixed value ofλ,(a)
to*decreases when γ is increased in order to reduce nonconformities by reducing the out of control period;(b)
μI* 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 γ.Table 2
Dependent of the optimal policy (to*,μI*,σ*) on parameters λ and γ.
γto*μI*σ*AC(to*,μI*,σ*)λ=0.041.0052.575254.2541.4641119.391.012.558253.4761.4831226.581.022.499252.5371.5811574.041.032.454254.0312.3191880.181.042.440255.3921.9411914.69λ=0.061.0052.581254.2631.4801194.871.012.563253.4861.5021300.151.022.503252.5581.6251641.171.032.464254.022.3471939.331.042.444255.3951.9611974.51λ=0.081.0052.584254.2721.4961268.161.012.566253.4961.5211371.571.022.507252.5781.6681706.181.032.473254.0082.3751996.691.042.448255.3971.9822032.53On the other hand, under a fixed value ofγ we notice the following:(a)
to* increases when λ is increased;(b)
μI* 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 lengthto* and the optimal process mean μI* simultaneously, while in the second model we introduce a new cost Cσ resulting from lowering the process standard deviation σ; accordingly we determine the optimal standard deviation σ* in addition to to* and μI*. In future research, deterioration can be considered to happen gradually following some pattern (e.g., linear or exponential) rather than being sudden.



---

*Source: 393495-2012-02-02.xml*