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 (LSL) and upper specification limit (USL). The cost of producing items within the accepted range, LSL<π‘₯<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 time 𝑑=0 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 𝑑<𝑑0; 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 𝑑<𝑑0, 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:πœ‡π‘‚=π›Ύπœ‡πΌ,𝛾>0.(1) 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:𝛼=𝛼1+𝛼2=ξ€œπ‘₯=LSLπ‘₯=0π‘žξ€·π‘₯,πœ‡πΌξ€Έξ€œπ‘‘π‘₯+π‘₯=∞π‘₯=USLπ‘žξ€·π‘₯,πœ‡πΌξ€Έπ‘‘π‘₯,𝛽=𝛽1+𝛽2=ξ€œπ‘₯=LSLπ‘₯=0π‘žξ€·π‘₯,πœ‡π‘‚ξ€Έξ€œπ‘‘π‘₯+π‘₯=∞π‘₯=USLπ‘žξ€·π‘₯,πœ‡π‘‚ξ€Έπ‘‘π‘₯.(2) The expected time while the process is in-control is given by𝑇in=ξ€œπ‘‘π‘œπ‘‘=0ξ€œπ‘‘πœ=0+ξ€œπœβ‹…β„Ž(𝜏)⋅𝑓(𝑑)β‹…π‘‘πœβ‹…π‘‘π‘‘βˆžπ‘‘=π‘‘π‘œξ€œπ‘‘π‘œπœ=0πœβ‹…β„Ž(𝜏)⋅𝑓(𝑑)β‹…π‘‘πœβ‹…π‘‘π‘‘.(3) The expected time while the process is out-of-control is given by𝑇out=ξ€œπ‘‘π‘œπ‘‘=0ξ€œπ‘‘πœ=0(+ξ€œπ‘‘βˆ’πœ)β‹…β„Ž(𝜏)⋅𝑓(𝑑)β‹…π‘‘πœβ‹…π‘‘π‘‘βˆžπ‘‘=π‘‘π‘œξ€œπ‘‘π‘œπœ=0ξ€·π‘‘π‘œξ€Έβˆ’πœβ‹…β„Ž(𝜏)⋅𝑓(𝑑)β‹…π‘‘πœβ‹…π‘‘π‘‘.(4) The total cost of nonconforming items in a production cycle is given by𝐢𝑁=𝑝⋅𝑇in⋅𝑐𝐿𝛼1+π‘π‘ˆπ›Ό2ξ€»+𝑝⋅𝑇out⋅𝑐𝐿𝛽1+π‘π‘ˆπ›½2ξ€».(5) We quantify the impact of imprecision on the quality characteristic 𝑋 by using Taguchi loss function:⎧βŽͺβŽͺ⎨βŽͺβŽͺβŽ©π‘πΏ(π‘₯)=𝐿π‘₯βˆ’πœ‡πΌξ€Έ2ξ€·LSLβˆ’πœ‡πΌξ€Έ2,ifLSL≀π‘₯β‰€πœ‡πΌ,π‘π‘ˆξ€·π‘₯βˆ’πœ‡πΌξ€Έ2ξ€·USLβˆ’πœ‡πΌξ€Έ2,ifπœ‡πΌβ‰€π‘₯≀USL.(6) The total loss cost in a production cycle is given by𝐢Loss=𝑝⋅𝑇inβ‹…ξ‚Έξ€œπœ‡πΌπ‘₯=LSL𝐿(π‘₯)β‹…π‘žπ‘₯,πœ‡πΌξ€Έ+ξ€œβ‹…π‘‘π‘₯USLπ‘₯=πœ‡πΌξ€·πΏ(π‘₯)β‹…π‘žπ‘₯,πœ‡πΌξ€Έξ‚Ήβ‹…π‘‘π‘₯+𝑝⋅𝑇outβ‹…ξ‚Έξ€œπœ‡πΌπ‘₯=LSL𝐿(π‘₯)β‹…π‘žπ‘₯,πœ‡π‘‚ξ€Έ+ξ€œβ‹…π‘‘π‘₯USLπ‘₯=πœ‡πΌξ€·πΏ(π‘₯)β‹…π‘žπ‘₯,πœ‡π‘‚ξ€Έξ‚Ή.⋅𝑑π‘₯(7) 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𝐢cycle=π‘π‘œ+𝑐1β‹…ξ€œπ‘‘π‘œπ‘‘=0ξ€œβˆžπ‘™1=0𝑙1⋅𝑔1𝑙1⋅𝑓(𝑑)⋅𝑑𝑙1⋅𝑑𝑑+𝑐2β‹…ξ€œβˆžπ‘‘=π‘‘π‘œξ€œβˆžπ‘™2=0𝑙2⋅𝑔2𝑙2⋅𝑓(𝑑)⋅𝑑𝑙2+𝑐⋅𝑑𝑑𝑖𝑝(π‘βˆ’π‘‘)β‹…ξ‚Έξ€œ2π‘‘π‘‘π‘œπ‘‘=0𝑑2⋅𝑓(𝑑)⋅𝑑𝑑+𝑑2π‘œβ‹…ξ€œβˆžπ‘‘=π‘‘π‘œξ‚Ήπ‘“(𝑑)⋅𝑑𝑑+π‘π‘ ξ€œπ‘‘β‹…π‘‘π‘œπ‘‘=0ξ€œβˆžπ‘™1=(π‘βˆ’π‘‘)𝑑/𝑑𝑙1βˆ’(π‘βˆ’π‘‘)𝑑𝑑⋅𝑔1𝑙1⋅𝑓(𝑑)⋅𝑑𝑙1⋅𝑑𝑑+π‘π‘ ξ€œπ‘‘β‹…βˆžπ‘‘=π‘‘π‘œξ€œβˆžπ‘™2=(π‘βˆ’π‘‘)π‘‘π‘œ/𝑑𝑙2βˆ’(π‘βˆ’π‘‘)π‘‘π‘œπ‘‘ξ‚Όβ‹…π‘”2𝑙2⋅𝑓(𝑑)⋅𝑑𝑙2⋅𝑑𝑑+𝑝⋅𝑇in⋅𝑐𝐿𝛼1+π‘π‘ˆπ›Ό2ξ€»+𝑝⋅𝑇out⋅𝑐𝐿𝛽1+π‘π‘ˆπ›½2ξ€»+𝑝⋅𝑇inβ‹…ξ‚Έξ€œπœ‡πΌπ‘₯=LSL𝐿(π‘₯)β‹…π‘žπ‘₯,πœ‡πΌξ€Έ+ξ€œβ‹…π‘‘π‘₯USLπ‘₯=πœ‡πΌξ€·πΏ(π‘₯)β‹…π‘žπ‘₯,πœ‡πΌξ€Έξ‚Ήβ‹…π‘‘π‘₯+𝑝⋅𝑇outβ‹…ξ‚Έξ€œπœ‡πΌπ‘₯=LSL𝐿(π‘₯)β‹…π‘žπ‘₯,πœ‡π‘‚ξ€Έ+ξ€œβ‹…π‘‘π‘₯USLπ‘₯=πœ‡πΌξ€·πΏ(π‘₯)β‹…π‘žπ‘₯,πœ‡π‘‚ξ€Έξ‚Ή.⋅𝑑π‘₯(8) And the expected length of a production-inventory cycle [15] is given by𝑇cycle=ξ€œπ‘‘π‘œπ‘‘=0ξ€œπ‘™(π‘βˆ’π‘‘)𝑑/𝑑1=0𝑝𝑑𝑑⋅𝑔1𝑙1⋅𝑓(𝑑)⋅𝑑𝑙1+ξ€œβ‹…π‘‘π‘‘π‘‘π‘œπ‘‘=0ξ€œβˆžπ‘™1=(π‘βˆ’π‘‘)𝑑/𝑑𝑑+𝑙1⋅𝑔1𝑙1⋅𝑓(𝑑)⋅𝑑𝑙1+ξ€œβ‹…π‘‘π‘‘βˆžπ‘‘=π‘‘π‘œξ€œ(π‘βˆ’π‘‘)π‘‘π‘œπ‘™/𝑑2=0ξ‚΅π‘π‘‘π‘œπ‘‘ξ‚Άβ‹…π‘”2𝑙2⋅𝑓(𝑑)⋅𝑑𝑙2+ξ€œβ‹…π‘‘π‘‘βˆžπ‘‘=π‘‘π‘œξ€œβˆžπ‘™2=(π‘βˆ’π‘‘)π‘‘π‘œ/π‘‘ξ€·π‘‘π‘œ+𝑙2⋅𝑔2𝑙2⋅𝑓(𝑑).⋅𝑑𝑙2⋅𝑑𝑑.(9) From the renewal reward theorem, the average cost per unit time is given byξ€·π‘‘π΄πΆπ‘œ,πœ‡πΌξ€Έ=𝐢cycleξ€·π‘‘π‘œ,πœ‡πΌξ€Έπ‘‡cycleξ€·π‘‘π‘œξ€Έ.(10) 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; ⎧βŽͺ⎨βŽͺ⎩1β„Ž(𝜏)=𝑑,π‘‘β‰€π‘‘π‘œ,1π‘‘π‘œ,𝑑β‰₯π‘‘π‘œ.(11) In order to avoid mathematical intractability, we assume exponential time to failure with failure rateπœ†: 𝑓(𝑑)=πœ†π‘’βˆ’πœ†π‘‘.(12) 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: 𝑔1𝑙1ξ€Έ=1𝑏1,0≀𝑙1≀𝑏1,𝑔2𝑙2ξ€Έ=1𝑏2,0≀𝑙2≀𝑏2.(13) And we restrict the decision variable π‘‘π‘œ to assume values within the following interval; 0β‰€π‘‘π‘œβ‰€4.(14) We use the following values: 𝑝=130, 𝑑=100, π‘π‘œ=300, 𝑐1=1000, 𝑐2=200,𝑏1=3, 𝑏2=1, 𝑐𝑖=8, 𝑐𝑠=40, 𝑐𝐿=30, π‘π‘ˆ=20, LSL=250, USL=260, 𝜎=2, 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, 𝑐1, 𝑐2, 𝑐𝑖 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𝐿𝑁=π‘šπ‘‘π‘œ,(15) where πΏπ‘š is the expected life of the machine. The standard deviation cost in a production cycle can be expressed as𝐢𝜎=𝐾𝑒(𝜎maxβˆ’πœŽ)𝑁=πΎπ‘‘π‘œπ‘’(𝜎maxβˆ’πœŽ)πΏπ‘š.(16) Here, 𝐾 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𝐢cycle=π‘π‘œ+𝑐1β‹…ξ€œπ‘‘π‘œπ‘‘=0ξ€œβˆžπ‘™1=0𝑙1⋅𝑔1𝑙1⋅𝑓(𝑑)⋅𝑑𝑙1⋅𝑑𝑑+𝑐2β‹…ξ€œβˆžπ‘‘=π‘‘π‘œξ€œβˆžπ‘™2=0𝑙2⋅𝑔2𝑙2⋅𝑓(𝑑)⋅𝑑𝑙2+𝑐⋅𝑑𝑑𝑖𝑝(π‘βˆ’π‘‘)β‹…ξ‚Έξ€œ2π‘‘π‘‘π‘œπ‘‘=0𝑑2⋅𝑓(𝑑)⋅𝑑𝑑+𝑑2π‘œβ‹…ξ€œβˆžπ‘‘=π‘‘π‘œξ‚Ήπ‘“(𝑑)⋅𝑑𝑑+π‘π‘ ξ€œπ‘‘β‹…π‘‘π‘œπ‘‘=0ξ€œβˆžπ‘™1=(π‘βˆ’π‘‘)𝑑/𝑑𝑙1βˆ’(π‘βˆ’π‘‘)𝑑𝑑⋅𝑔1𝑙1⋅𝑓(𝑑)⋅𝑑𝑙1⋅𝑑𝑑+π‘π‘ ξ€œπ‘‘β‹…βˆžπ‘‘=π‘‘π‘œξ€œβˆžπ‘™2=(π‘βˆ’π‘‘)π‘‘π‘œ/𝑑𝑙2βˆ’(π‘βˆ’π‘‘)π‘‘π‘œπ‘‘ξ‚Όβ‹…π‘”2𝑙2⋅𝑓(𝑑)⋅𝑑𝑙2⋅𝑑𝑑+𝑝⋅𝑇in⋅𝑐𝐿𝛼1+π‘π‘ˆπ›Ό2ξ€»+𝑝⋅𝑇out⋅𝑐𝐿𝛽1+π‘π‘ˆπ›½2ξ€»+𝑝⋅𝑇inβ‹…ξ‚Έξ€œπœ‡πΌπ‘₯=LSL𝐿(π‘₯)β‹…π‘žπ‘₯,πœ‡πΌξ€Έ+ξ€œβ‹…π‘‘π‘₯USLπ‘₯=πœ‡πΌπΏξ€·(π‘₯)β‹…π‘žπ‘₯,πœ‡πΌξ€Έξ‚Ήβ‹…π‘‘π‘₯+𝑝⋅𝑇outβ‹…ξ‚Έξ€œπœ‡πΌπ‘₯=LSL𝐿(π‘₯)β‹…π‘žπ‘₯,πœ‡π‘‚ξ€Έ+ξ€œβ‹…π‘‘π‘₯USLπ‘₯=πœ‡πΌπΏξ€·(π‘₯)β‹…π‘žπ‘₯,πœ‡π‘‚ξ€Έξ‚Ή+⋅𝑑π‘₯πΎπ‘‘π‘œπ‘’(𝜎maxβˆ’πœŽ)πΏπ‘š.(17) And the cycle length is given by (9).

And again by the renewal theory, the average cost per unit time is given by𝑑ACπ‘œ,πœ‡πΌξ€Έ=𝐢,𝜎cycleξ€·π‘‘π‘œ,πœ‡πΌξ€Έ,πœŽπ‘‡cycleξ€·π‘‘π‘œξ€Έ.(18)

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𝐾=5000, πΏπ‘š=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 π‘‘π‘œ and 𝜎 to assume values within the following intervals: 0β‰€π‘‘π‘œβ‰€4, 0β‰€πœŽβ‰€πœŽmax.

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
𝑙1:Nonenegative random variable denoting corrective repair time
𝑔1(𝑙1):Corrective repair time probability density function p.d.f.
𝑏1:The upper bound on 𝑙1 when it is uniformly distributed
𝑙2:None-negative random variable denoting preventive repair time
𝑔2(𝑙2):Preventive repair time probability density function p.d.f.
𝑏2:The upper bound on 𝑙2 when it’s uniformly distributed
𝑇in:The expected time while process is in-control
𝑇out:The expected time while process is out of control
𝑑(>0):Demand rate
𝑝(>𝑑):Production rate
𝑐0(>0):Set up cost for each production run
𝑐1(>0):Corrective repair cost per unit time
𝑐2(<𝑐1):Preventive repair cost per unit time
𝑐𝑖(>0):Inventory holding cost per unit product per unit time
𝑐𝑠(>0):Shortage penalty cost per unit product
𝑐𝐿(>0):Cost incurred due to production of a nonconforming item with𝑋≀LSL
π‘π‘ˆ(>0):Cost incurred due to production of a nonconforming item with𝑋β‰₯USL
𝐢𝑁:Total cost of nonconforming items per cycle
𝐿(π‘₯):Taguchi loss function
𝐢Loss:Total loss cost per cycle
𝐢cycle:Expected total cost per production-inventory cycle
𝑇cycle:Expected length of a production-inventory cycle
AC: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
𝜎max: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
LSL:Lower specification limit on quality characteristic 𝑋
USL:Upper specification limit on quality characteristic 𝑋
𝛼1:Proportion of nonconforming while process is in-control with𝑋≀LSL
𝛼2Proportion of nonconforming while process is in-control with𝑋β‰₯USL
𝛽1:Proportion of nonconforming while process is out-of-control with𝑋≀LSL
𝛽2:Proportion of nonconforming while process is out-of-control with𝑋β‰₯USL.