Mathematical modeling of a multi-product EMQ model with an enhanced end items issuing policy and failures in rework.

This study uses mathematical modeling to examine a multi-product economic manufacturing quantity (EMQ) model with an enhanced end items issuing policy and rework failures. We assume that a multi-product EMQ model randomly generates nonconforming items. All of the defective are reworked, but a certain portion fails and becomes scraps. When rework process ends and the entire lot of each product is quality assured, a cost reduction n + 1 end items issuing policy is used to transport finished items of each product. As a result, a closed-form optimal production cycle time is obtained. A numerical example demonstrates the practical usage of our result and confirms a significant savings in stock holding and overall production costs as compared to that of a prior work (Chiu et al. in J Sci Ind Res India, 72:435-440 2013) in the literature.

Background

Mathematical modeling is used in this study to examine a multi-product EMQ model with rework failures and an enhanced cost reduction end items issuing policy. The EMQ model made use of a mathematical technique to balance the setup and holding costs incurred in a production cycle, and derive most economic manufacturing quantity that minimizes the long run average system costs per unit time (Taft 1918). The assumption of traditional EMQ model includes a perfect manufacturing process for a single product, and a continuous finished product distribution policy. Although the assumptions of the EMQ model are simple and somehow unrealistic, its concept along with solution procedure has since been extensively applied to the fields of inventory control and production management (Hadley and Whitin 1963; Silver et al. 1998; Nahmias 2009; Battini et al. 2010a; Andriolo et al. 2014; Azzi et al. 2014; Glock et al. 2014). In order to increase machine utilization, the vendors in manufacturing sector often fabricate multiple products in sequence on a single machine. Rosenblatt and Finger (1983) considered a single machine multi-item production problem, whereas the machine was an electrochemical machining system, and its outputs are impact sockets of different sizes for power wrenches. They used a grouping procedure for various different products along with a modified version of an existing algorithm to confirm that the cycle times are the multiples of the shortest cycle time. Federgruen and Katalan (1998) examined stochastic economic batch scheduling problems with periodic base-stock policies, where all products are fabricated according to a given periodic item-order. They proposed some effective heuristics to minimize system-wide costs for such a periodic item-sequence production. Muramatsu et al. (2013) studied a multi-item multi-process dynamic lot size scheduling problem with setup time and general product structure. Various heterogeneous decision features such as lot sizing, lot sequencing, dispatching are considered. A near-optimal solution procedure was proposed to determine the decision features involved in this problem simultaneously. Caggiano et al. (2009) proposed a method for computing the channel fill rates in a multi-product, multi-echelon service parts distribution system. A simulation approach was employed to study multi-item three-echelon production- distribution system. They showed that the estimation errors are insignificant over a wide range of base stock level vectors and they also presented an enhanced approximation method to the problem. Jodlbauer and Reitner (2012) investigated a stochastic make-to-order multi-product manufacturing system under a common cycle policy. Effects of the safety stock, demand, cycle time, and setup time on the service level and on the total system costs were investigated. Papers that related to various aspects of planning and optimization issues on multi-item production can also be referred to Lin et al. (2014); Wu et al. (2014).

In real manufacturing environments, due to unpredictable factors production of defective products is inevitable. Sometimes, certain portion of nonconforming products can be reworked and repaired with additional repair expense. Agnihothri and Kenett (1995) studied the effects of defects on various system performance measures for a manufacturing process with rework. They offered management guidelines for allocating additional budget in process improvement for continuously increasing yield rate, identifying potential bottlenecks under increasing workloads, and providing extra resources to release bottlenecks. Extra studies that addressed various aspects of specific features including imperfect quality production and the rework processes can also be found elsewhere (Grosfeld-Nir and Gerchak 2002; Biswas and Sarker 2008; Battini et al. 2010b; Glock 2011a, 2012a; Chiu and Chang 2014; Wee et al. 2014; Murugan and Selladurai 2014; Hishamuddin et al. 2014; Battini et al. 2014; Chiu et al. 2014). The conventional EMQ model assumes a simplified ‘continuous product issuing policy’. However, in real supply chains environments multiple or periodic product delivery policy is commonly adopted. Banerjee (1986) examined a joint economic lot-size model for vendor and buyer, with the focus on minimizing the total joint relevant cost. He concluded that a jointly optimal ordering policy along with an appropriate price adjustment, could be economically beneficial for both buyer and vendor. Swenseth and Godfrey (2002) revealed that the freight rate functions can be combined into inventory replenishment decisions without lowering the accuracy of decisions, nor will these functions increase the complexity of the decision making process. Abdul-Jalbar et al. (2008) studied a one-vendor multi-buyer multi- echelon finite production rate system. The objectives of their study are to determine the optimal production and shipment schedule, and the most economic order size for buyers so that average total cost per unit time can be minimized. In a recent study, Chiu et al. (2013) determined optimal common production cycle time for a multi-item production system with discontinuous n fixed quantity multiple delivery policy and rework failures. Mathematical modeling and optimization techniques are employed in their study to solve the problem. As a result, a closed-form optimal common cycle time that minimizes the expected system costs is obtained. Effect of rework failure on the optimal cycle time was investigated through a numerical example. Additional studies that addressed various aspects of periodic or multi-delivery issues of vendor–buyer integrated systems can be referred to (Glock 2011b, 2012b; Katsaliaki et al. 2014; Safaei 2014; Rodger 2014; Sana et al. 2014).

For the purpose of lowering producer’s inventory holding cost as well as minimizing the expected overall system cost, we extend Chiu et al.’s (2013) work by replacing their n fixed quantity end items delivery policy with an enhanced n + 1 issuing policy. Under the new policy, one extra delivery of end items is made in producer’s production uptime to satisfy customers’ demands during production uptime and rework time. Then, upon completion of the rework, additional n installments (fixed quantity) of end items are shipped at a fixed time interval. The objectives of this study are to determine optimal common cycle time that minimizes the long-run average system cost per unit time, and investigate effects of random defective rate, rework failures, and the enhanced end items issuing policy on the optimal operating cycle time as well as on the expected system costs per unit time.

Problem description and mathematical modeling

We use mathematical modeling to examine a multi-product EMQ system with a cost reduction multi-distribution policy and rework failures. Consider there are L products to be fabricated in sequence on a single machine and there is a xi portion of nonconforming items being randomly produced at a rate d1i during the production of product i (where ). All items are screened and cost of quality inspection is included in unit production cost C. Under the normal operation (which shortage are not permitted), the constant production rate P1i for product i must satisfies , where λi denotes demand rate for product i per year, and d1i = xiP1i. All nonconforming items are reworked at a rate of P2i each cycle immediately following the end of regular production process. The additional rework cost is CRi per item. A failure-in-rework rate φi exists during the reworking process, the production rate of scrap items during rework d2i = φiP2i, and those that fail in repair are discarded at a unit disposal cost CSi. A specific n + 1 multi-shipment policy is proposed in an attempt to reduce the vendor’s stock holding cost as compared to n multi-delivery policy used in Chiu et al. (2013). Under the proposed n + 1 delivery policy, the purpose of the first shipment of finished goods is to meet buyers’ product demands during vendor’s uptime and reworking time. Then upon completion of rework process, n fixed quantity installments of the end items are distributed to buyers at a fixed time interval tn (see Fig. 1).

Fig. 1

On-hand inventory of perfect quality items in the proposed multi-item production system with rework failures and a cost reduction distribution policy

Figure 2 depicts producer’s on-hand inventory level of perfect quality items of product i in the proposed n + 1 delivery model (in blue lines) and the expected reduction in vendor’s stock holding costs (yellow shaded area) in comparison with that of Chiu et al. (2013). Cost related variables used in our analysis comprise the setup cost K per cycle, unit holding cost h, unit holding cost h1 during rework, the fixed shipping cost K1 for product i per delivery, and unit shipping cost CT for each product i. Additional variables also include the following: T = common production cycle length, our decision variable, Q = batch size per cycle for product i, n = number of installments (fixed quantity) of the finished lot to be shipped to buyers per cycle, H1 = on-hand inventory in units of product i for meeting buyer’s demand during uptime t1 and reworking time t2, H2 = maximum level of on-hand inventory of product i when the regular production ends, H = maximum level of on-hand inventory in units of product i when the rework process ends, t = time required for producing enough items to meet demand of product i during vendor’s uptime t1 and reworking time t2, t1 = production uptime for product i, t2 = the reworking time for product i, t3 = the delivery time for product i, tn = fixed interval of time between each installment of finished product i being delivered during t3i, I(t) = level of on-hand perfect quality items at time t, IS(t) = level of on-hand scrap items of product i at time t, TC(Q) = total production-inventory-delivery cost per cycle for product i, E[TCU(T)] = the expected system costs per unit time for L products in the proposed system.

Fig. 2

Expected reduction in producer’s inventory holding costs (yellow shaded area) for each product i in the proposed model in comparison with that in Chiu et al. (2013)

From Fig. 1, the following equations (for ) can be obtained directly:The level of on-hand inventory of scrap items in the proposed model is depicted in Fig. 3 and the following equations (for ) can be obtained:

Fig. 3

On-hand inventory of scrap items in the proposed multi-item production system with rework failures and a cost reduction distribution policy

The variable holding cost for finished items of product i during delivery time t3 areThe fixed and variable transportation costs for product i per cycle areThe production-inventory-delivery cost per cycle for L products comprises the setup cost, the variable manufacturing, reworking, and disposal costs (Fig. 3), the fixed and variable shipping cost, the holding cost during t1i, t2i, and t3. Therefore, TC(Qi) becomes

Because the defective rate x is assumed to be a random variable with a known probability density function, in order to take the randomness of x into account, the expected value of x is used in this study. By substituting all parameters from Eqs. (1) to (13) in Eq. (14), and with further derivations the expected E[TCU(T)] can be obtained as follows:Let Then Eq. (15) becomes

Derivation of the optimal cycle time

Before derivation of the optimal common production cycle time T*, one should prove that the expected cost function E[TCU(T)] is convex. By differentiating E[TCU(T)] with respect to T gives the following first and second derivatives:It can be seen that Eq. (19) is positive, since Ki, n, K1i, and T are all positive. Because the second derivative of E[TCU(T)] > 0, one confirms that E[TCU(T)] is convex for all T different from zero. It follows that by letting the first derivative of E[TCU(T)] = 0, one can derive the optimal common production cycle time T*. LetorTherefore, one has T* as follows:

Capacity and setup time effects on the optimal cycle time

Generally speaking, setup time is relatively short in comparison with uptime. But, if the setup time becomes a factor, one has to ensure that the cycle length is long enough to account for the setup, production, and reworking times of L products (Nahmias 2009). Let S denote the production setup time for product i, the Eq. (23) must hold.Substituting Eq. (1) in Eq. (23) one hasTherefore, when setup time becomes a significant factor in production planning, one must choose the optimal common cycle length from the maximum of [T*, Tmin].

Numerical example

With the purpose of easing the comparison efforts for readers, this section uses the same example as in Chiu et al. (2013) to demonstrate the proposed research result. Reconsider a production plan of producing five different items on a single machine in sequence under a common cycle time policy. Annual demand rates λ for these five items are 3000, 3200, 3400, 3600, and 3800, respectively. The production rates P1i are 58,000; 59,000; 60,000; 61,000 and 62,000, respectively. During individual production process, there are random nonconforming rates x for each item and they follow Uniform distribution over the intervals of [0, 0.05], [0, 010], [0, 0.15], [0, 020], and [0, 0.25], respectively. All nonconforming items produced go through a rework process at the rates P2i of 1800, 2000, 2200, 2400, and 2600 items per year, respectively; with additional unit rework costs CR of $50, $55, $60, $65, and $70, respectively. During the reworking, there are failure-in-rework rates φi of 0.05, 0.10, 0.15, 0.20, and 0.25, respectively. Additional values of system parameters are given as follows: K = the setup costs are $3800, $3900, $4000, $4100, and $4200, respectively, C = production costs per item are $80, $90, $100, $110, and $120, respectively, CS = disposal costs per item are $20, $25, $30, $35, and $40, respectively, K1 = fixed costs per delivery are $1800, $1900, $2000, $2100, and $2200, respectively, CT = unit transportation costs are $0.1, $0.2, $0.3, $0.4, and $0.5, respectively, n = number of shipments per cycle, it is assumed to be a constant 3 (i.e., n + 1 = 4), h = unit holding costs are $10, $15, $20, $25, and $30, respectively, h1 = holding costs per reworked item are $30, $35, $40, $45, and $50, respectively.

The optimal common cycle time T* = 0.7279 (years) can be obtained by applying Eq. (22). Total expected system costs E[TCU(T*)] = $2,013,956 can also be obtained from computation of Eq. (15). Variation of mean defective rate and mean failure-in-rework rate effects on the expected system cost E[TCU(T)] is illustrated in Fig. 4. It is noted that as mean defective rate increases the E[TCU(T)] increases significantly, and as mean failure-in-rework rate increases the system cost E[TCU(T)] increases slightly.

Fig. 4

Variation of mean defective rate and mean failure-in-rework rate effects on the expected system cost E[TCU(T)]

As stated earlier, the proposed model aims at reducing vendor’s inventory holding cost for each product i during the production cycle. As a result from this numerical example, the percentage of overall holding cost reduction is 24.7 % (i.e., from $109,476 (Chiu et al. 2013) down to $82,431). Figure 5 demonstrates the percentage of holding cost reduction for five different products, respectively as compared to that of Chiu et al.’s work (where n-delivery policy is adopted).

Fig. 5

Percentage of reduction in producer’s inventory holding costs for five different products in example, respectively

In summary, the proposed study realizes a significant system cost savings of $56,358 (i.e., $2,070,314 − $2,013,956) or 16.09 % of other system interrelated costs (i.e., E[TCU(T)] − λC, which is the expected system cost excludes variable manufacturing cost).

Conclusions

With an aim at reducing producer’s inventory cost as well as minimizing the expected overall system costs per unit time, this study incorporating an enhanced n + 1 product issuing policy into Chiu et al.’s model (2013), and with the help of mathematical modeling and optimization method a closed-form optimal common production cycle time for the proposed multi-product EMQ model with rework failures was derived. A numerical example is given to demonstrate the applicability of our research result, reveal joint effects of random defective rate and failure-in-rework rate on the optimal policy (refer to Fig. 4), and confirm significant savings in producer’s inventory holding cost (Fig. 5). For future study, an interesting topic will be to include the machine breakdown factor into such a multi-product EMQ model.