SciELO - Scientific Electronic Library Online

 
vol.18 issue4Constructive heuristic for the vertex bisection problemReal time localization solution for land vehicle application using low-cost integrated sensors with GPS author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Journal of applied research and technology

On-line version ISSN 2448-6736Print version ISSN 1665-6423

J. appl. res. technol vol.18 n.4 Ciudad de México Aug. 2020  Epub July 31, 2021

https://doi.org/10.22201/icat.24486736e.2020.18.4.1195 

Artículos

Optimal hybrid inventory replenishment runtime for a vendor- buyer coordinated system with breakdowns and rework/disposal of nonconforming items

Singa Wang Chiua 

Hua-Yao Wub 

Peng-Cheng Sungc  * 

aDept. of Business Administration, Chaoyang University of Technology, Taichung 413, Taiwan

bPhysics Dept., College of Liberal Arts and Sciences, State University of New Yorkat Oswego, NY, USA

cDept. of Industrial Engineering & Management, Chaoyang University of Technology, Taichung,


Abstract:

Most multinational/transnational manufacturers that own internal supply chains and operate in turbulent world markets apply strategies related to product quality, low-cost, and timely delivery. This study aims to assist such firms with making accurate decisions to enable their competitive strategies and cope with the realities of limited capacity and unreliable equipment. We examine a vendor-buyer coordinated system featuring batch fabrication, outsourcing, quality reassurance, discontinuous deliveries, and an unreliable machine. The system outsources a portion of a lot to reduce manufacturing uptime, and the in-house fabrication system experiences undesirable defective items and Poisson distributed breakdowns. In each cycle, corrective action and rework/disposal of defective stocks are undertaken as these incidents occur, and upon receipt of outsourced products and when the entire batch is quality ensured, it makes multiple deliveries of the end products. Using modeling, formulation, derivation, and an optimization methodology, we obtain the problem’s cost function and justify its convexity. We then apply differential calculus and propose a recursive algorithm to derive the problem’s optimal replenishment runtime. A numerical illustration is offered to show the applicability of the result that reveals various important system characteristics/ capabilities, such as the distinct and combined influences of breakdowns, outsourcing, rework, scrap, and delivery-frequency factors on various system parameters, performance, and optimal runtime. The methods proposed here can facilitate managerial operations planning and strategic decision making in an intra-supply chain setting in practice.

Keywords: Industrial engineering; Replenishment runtime; Breakdown; Outsourcing; Quality reassurance; Vendor-buyer coordinated system; Multiple deliveries

1. Introduction

This study examines a vendor-buyer coordinated system featuring batch fabrication, outsourcing, quality reassurance, discontinuous deliveries, and an unreliable machine. Most real-life manufacturing systems experience unanticipated nonconforming items and breakdowns. As such instances occur, corrective action and rework/disposal of nonconforming items must be undertaken to avoid delay in the fabrication schedule and attain the desired product quality. Vinod and Solberg (1984) examined the single- and multi-stage unreliable fabrication systems using queueing models. The authors derived the exact solution for the queueing model with a single-stage and presented two approximations for the closed network queueing model with multiple stages. Their approximation results/performances were validated/compared against the exact solution in the literature. Groenevelt, Pintelon, and Seidmann (1992) studied an unreliable fabrication facility with safety stocks and batch production, wherein a constant failure rate and random failure-repair time of the facility were assumed. Diverse bounds of service-level were examined to find their impacts on different system variables. A production control discipline was proposed to investigate the relationship between the safety stocks and the renewal process of a particular type of single server queue. The authors also showed how their approaches could be applied to broader decision makings in resource allocation fields. Dohi, Okamura, and Osaki (2001) considered an economic manufacturing quantity (EMQ) model with preventive maintenance, stochastic facility breakdowns, (PM), and safety stocks. Their purpose was to jointly decide the optimal control policies of the PM schedule and the quantities of safety stocks that keep the total cost at a minimum. Besides, the authors found out that both the safety stocks and total cost rise as breakdown rate increases. Chakraborty, Giri, and Chaudhuri (2009) studied the production lot-size problem considering breakdowns and different inspection schedules for a deteriorating process. Corrective action of breakdown situation and preventive maintenance are undertaken accordingly. The authors proposed models based on general shift and various distributions of machine failure and repair times to derive a suboptimal lot-size policy. Numerical examples were offered to show the results’ applicability and sensitivity analyses on system performances with/without inspection policy. Goerler and Voß (2016) used a mixed-integer programming approach to explore the capacitated batch-size problem with defective products and rework processes. Various numerical experiments were performed to explore the influences of changes in defective instances on the required computer times for obtaining the optimal batch-size solutions. Additional works (Al-Bahkali & Abbas, 2018; Arun, Lincon, & Prabhakaran, 2019; Ghalme, Mankar, & Bhalerao, 2017; Richter, 1996; Saari & Odelius, 2018; Sarker, Jamal, & Mondal, 2008; Shakoor, Abu Jadayil, Jaber, & Jaber, 2017; Souha, Soufien, & Mtibaa, 2018; Vujosevic, Makajic-Nikolic, & Pavlovic, 2017; Zahraee, Rohani, & Wong, 2018) studied the impact of various characteristics of unreliable facility and rework/disposal of defective products on fabrication systems and operations management.

Production managers apply an outsourcing strategy to effectively reduce fabrication uptime or release in-house facility’s workloads. Vining and Globerman (1999) presented a conceptual structure for comprehending the correct and less risky outsourcing decision. Specifically, through identifying the pre- and post-outsourcing risks and implementing certain suggested strategies to avoid or lessen those potential risks in advance (or pre-outsourcing stage). The authors referred to transaction costs in the literature to support their conceptual framework. De Fontenay and Gans (2008) considered a bargaining perception on strategic subcontracting and supply competition, wherein the subcontracting decision of a downstream company and its upstream fabrication resources is involved. The authors portrayed a downstream company has a choice to either subcontract to a reputable upstream company or a new and independent firm. Hence, it faces a trade-off between the higher resource value linked to those who could consolidate upstream capabilities and the lower input costs afforded by an independent competition. The result of their study indicates that outsourcing to an established firm is more beneficial. Rosar (2017) explored the connection between strategic subcontracting and optimal purchase policy. First, a subcontracting choice that relies on a non-cost-savings mechanism was presented and analyzed. Then, the author extended it to a cost-savings relating rationale, with the discussion of the incentives of sellers who employ in nested subcontracting policies. Additional works (Chiu, Liu, & Hwang, 2017; Chiu, Chiu, Lin, & Chang, 2019a; Mohammadi, 2017; Skowronski & Benton, 2018) investigated the impact of distinct outsourcing characteristics on the manufacturing systems and enterprise management.

In real supply-chain environments, the transportation of goods is commonly planned using multi-shipment at specific time intervals. Thomas and Griffin (1996) examined the conventional business processes in the stages of procurement, fabrication, and distribution, and indicated the need for coordinating these stages as a supply chain. The authors suggested taking advantage of recent progress in advance communication technology to place specific emphasis on the effective management of the coordinated supply-chain model to reduce overall operating costs. Swenseth and Godfrey (2002) examined the stock refilling decisions incorporating certain transportation cost functions from the literature and showed that no unnecessary complexity was added to the decision process, nor loss in accuracy of the decision. Farsijani, Nikabadi, and Ayough (2012) employed the simulated annealing methodology to explore a multiproduct economic production quantity (EPQ) model with discrete shipping orders and space constraints. Their batch- production model also considered realistic factors such as the imperfect manufacturing process, rework of defective stocks, and allowable shortages. The LINGO package helped solve the linear examples, and the simulated annealing methodology assisted in resolving the non-linear combinatorial optimization examples. Montarelo, Glardon, and Zufferey (2017) used the Tabu searching metaheuristic to investigate a four-echelon stock management decision in a decentralized supply chain setting. The authors proposed a global simulation methodology and set different service levels to deal with the market’s random demands, to explore/optimize the four-echelon linear/nonlinear supply chains. Their result showed that there are substantial differences among echelons in crucial stock and cost parameters. The authors claimed their approach could be generalized for boarder applications. Additional works (Arabi, Dehshiri, & Shokrgozar, 2018; Bolaños, Escobar, & Echeverri, 2018; Chiu, Wu, & Tseng, 2019b; Morales, Franco, & Mendez-Giraldo, 2018; Nielsen & Saha, 2018; Paz, Granada-Echeverri, & Escobar, 2018; Puška, Kozarević, Stević, & Stovrag, 2018; Stažnik, Babić, & Bajor, 2017; Zhao, Qian, Nakamura, & Nakagawa, 2018) studied the influence of distinct features of multiple deliveries on various types of manufacturing-transportation and supply-chain systems. Few prior works have investigated the joint influence of breakdowns, outsourcing, multiple deliveries, and rework/disposal of defective stocks on the optimal batch-fabrication runtime decision, this work aims to fill the gap.

2. The proposed model

This study determines the optimal hybrid inventory replenishment runtime for a vendor-buyer coordinated system with the breakdown, outsourcing, multiple deliveries, and rework/disposal of defective items. Suppose the annual demand rate λ of a manufactured product is supplied by a vendor at a fabrication rate of P1 units per year in a vendor-buyer coordinated system. To shorten fabrication uptime of the batch production plan, the vendor decides to outsource a π portion of the batch size Q (where 0 < π < 1). Thus, Kπ and Cπ denote the fixed and unit costs relating to the vendor’s outsourcing policy. The relationship between outsourcing relevant parameters and their corresponding in-house variables is shown as follows:

Cπ=1+β2C (1)

Kπ=1+β1K (2)

where C and K represent the in-house fabrication unit and setup cost, respectively; and β2 and β1 denote the relating ratios between these variables. It is noted that when π = 1, our model turns into a “buy” rather than “make” model; in contrast, when π = 0, the proposed model becomes a purely in-house production model.

The in-house fabrication process may produce an x portion of nonconforming items randomly, at a rate d 1 (where d 1 = xP 1 and P 1 stands for the in-house manufacturing rate). To prohibit the stock-out situation, we assume that (P 1 - d 1 - λ) > 0. Careful inspection of the nonconforming items separates the rework-able from the scrap (where the scrap ratio θ1 among the nonconforming is assumed). In each batch fabrication cycle, a rework process immediately follows the regular manufacturing process, at a reworking rate of P 2 and an extra cost C R is associated with each reworked item. Also, we assume an imperfect rework process, a scrap ratio θ2 among the reworked items exists. Hence, the overall scrap rate is φ (which sums up to (θ1 + (1 - θ1) θ2) in each cycle, and all scraps are disposed with unit disposal cost C S. As to the outsourced products, we assume that their quality is guaranteed by the outside provider, and they are scheduled to be received at the end of the in-house rework process, before the beginning of the delivery time of finished goods.

Moreover, the in-house production machine is not reliable, it is subject to random failure (which follows the Poisson distribution, with β as mean per year). When a failure occurs (as shown in subsection 2.1), a specific abort/resume stock controlling policy is used. Its guideline is to instantly repair the failure and promptly resume fabrication of the interrupted/unfinished lot when the machine is restored. A constant failure repair time t r is assumed; in case that actual repair time is greater than t r, a piece of rental equipment will be put in use to avoid further delay in production. Upon completion of the fabrication and rework processes, and receipt of outsourced items, n equal-size installments of the lot are shipped to the buyer at fixed time interval t' during distribution time t' . The additional notation used in this study is listed below.

t-

time before a random failure occurs (in years),

Q-

batch size,

M-

machine repair cost,

t-

uptime in the proposed hybrid replenishment vendor-buyer coordinated system with random failure and quality assurance - the decision variable,

t'-

rework time in the failure occurrence case,

T'π-

cycle length in the failure occurrence case,

d2-

production rate of scrap items during t' ,

h-

perfect item’s unit holding cost,

h1-

reworked item’s unit holding cost,

h2-

buyer stock’s unit holding cost,

h3-

safety stock’s unit holding cost,

C1-

safety stock’s unit cost,

CT-

unit transportation cost,

K1-

fixed transportation cost,

g- tr,

fixed machine repair time,

D-

quantity per delivery,

I -

the leftover stocks in each delivery time interval,

H0-

level of perfect stocks when a failure occurs,

H1-

level of perfect stocks when the fabrication process ends,

H2-

level of perfect stocks when the rework process ends,

H-

level of perfect stocks after receipt of outsourced items,

I(t)-

level of perfect stocks at time t,

IF(t)-

level of safety stocks at time t,

Id(t)-

level of nonconforming stocks at time t,

Is(t)-

level of scrap at time t,

Ic(t)-

level of buyer’s stocks at time t,

TC(t)1 =

total system cost per cycle in the failure occurrence case,

E[TC(t)1] =

the expected total system cost per cycle in the failure occurrence case,

E[T'π] =

the expected cycle length in the failure occurrence case,

t-

rework time in the case of no failure occurrence,

t-

stock delivery time in the case of no failure occurrence,

t

time interval between any two deliveries in the case of no failure occurrence,

Tπ

cycle length in the case of no failure occurrence,

TC(t)2 =

total system cost per cycle in the case of no failure occurrence,

E[TC(t)2] =

the expected total system cost per cycle in the case of no failure occurrence,

E[TCU(t)] =

the expected system cost per unit time for the proposed system with or without failure occurrence,

E[T’π] =

the expected cycle length in the case of no failure occurrence,

t1-

uptime for the proposed system without breakdown, nor outsourcing,

t2-

rework time for the proposed system without breakdown, nor outsourcing,

t3-

delivery time for the proposed system without breakdown, nor outsourcing,

T-

cycle length for the proposed system without breakdown, nor outsourcing,

Tπ-

replenishment cycle length for the proposed system with or without failure occurrence.

The following subsections examine two distinct cases due to the random failure in the proposed model:

2.1. Case 1: A random failure occurs during fabrication uptime

2.1.1. During the fabrication process of Case 1

Figure 1 shows the level of perfect inventories in this case (i.e., t < t), wherein at the time when a failure happens, the level of inventory reaches H0 and once the failure is repaired, it continues to pile up to H1 at the end of uptime and reaches H2 at the end of rework process. Then, the outsourced items are received and the level of the perfect stock reaches to H, before the beginning of product distribution time t'.

Figure 1 Level of perfect inventories in the proposed hybrid inventory replenishment vendor- buyer coordinated system with random breakdown and quality reassurance (in brown) as compared to the proposed system without breakdown, nor outsourcing (in black). 

Figure 2 displays the on-hand level of safety stock in the proposed system. It indicates that in the failure occurrence case, the safety stock will be added to the finished batch and delivered in t' for meeting extra buyer’s demand during tr.

Figure 2 Level of safety stock in the proposed system with breakdown occurrence. 

Figures 3 and 4 illustrate the levels of nonconforming and scrap items in the proposed system with failure occurrence, respectively.

Figure 3 Level of nonconforming items in the proposed system with breakdown occurrence. 

Figure 4 Level of scrap items in the proposed system with breakdown occurrence. 

Based on the aforementioned description of the in-house fabrication process, one can observe the following straightforward equations (please refer to Figures 1 to 5):

Figure 5 Level of perfect inventories in the proposed hybrid inventory replenishment vendor- buyer coordinated system with quality reassurance, but no machine breakdown (in brown) as compared to the same system without outsourcing option (in black). 

T'π=t1π+tr+t'2π+t'3π (3)

t1π=Q1-πP1=H1P1-d1 (4)

t'2π=1-πQx1-θ1P2 (5)

t'3π=T'π-t1π+tr+t'2π (6)

H0=P1-d1t (7)

H1=P1-d1t1π (8)

H2=H1+P2-d2t'2π (9)

d1t1π=xP1t1π=x1-πQ (10)

φx1-πQ=θ1+1-θ1θ2x1-πQ. (11)

2.1.2. During the delivery time of Case 1

Total delivery quantity H at the beginning of product distribution time t' must include λtr as shown in Eq. (12). Total inventories during the product distribution time t' (Chiu et al., 2019b) is exhibited in Eq. (13).

H=H2+πQ+λtr (12)

1n2i=1n-1iHt'3π =1n2n(n-1)2Ht'3π =n-12nHt'3π (13)

2.1.3. The status of buyer’s stocks in Case 1

Total inventories at the buyer side during the cycle length T' π can be calculated (Chiu et al., 2019b) as shown in Eq. (14).

nt'nπD-λt'nπ2+nn-12It'nπ+nI2t1π+t'2π=12Ht'3πn+H-λt'3πT'π (14)

2.1.4. Total cost per cycle for Case 1

Total cost per cycle in failure occurrence case, TC(t )1 comprises both the variable and fixed outsourcing and in-house fabrication costs, machine repaired cost, safety stock relevant costs (see Fig. 2), both fixed and variable shipping costs, rework and disposal costs, and total holding costs (including perfect items, nonconforming and reworked items, and buyer’s stocks) during the entire cycle, as shown in Eq. (15).

TCt1π1=CππQ+Kπ+C1-πQ+K+M+C1λtr+h3λtrt1π+tr+t'2π+nK1+CTQ1-φx1-π+λtr+CRx1-πQ1-θ1+CSφx1-πQ+hH1+d1t1π2t1π+H0tr+d1ttr+H1+H22t'2π+n-12nHt'3π+h1P2t'2π2t'2π+h 22Ht'3πn+H-λt'3πT'π (15)

Substitute Equations (1) to (14) in Eq. (15), and use the expected value to cope with the randomness of x, the expected total system cost per cycle in the failure occurrence case E[TC(t )1] can be derived as follows:

ETCt1π1=Cππt1πP11-π+Kπ+Ct1πP1+K+M+nK1+CTt1πP11-πy0+λg+C1λg+h3λgt1π+λg2+λgExt1πP11-θ1P2+CRExt1πP11-θ1+CSExφt1πP1+hP1tg +Ex2t1π2P121-θ12P2h11-θ1-h+t1π2P122nλ1-πh2-hy0y1-y2+hgt1πP12y1-y2+ht1π2P122λ1-πy021-π+λExφ1-π-πP1+λEx1-θ11-2πP2+h2λg22+h2gt1πP12y1+y2+h2t1π2P1221-πy0y2λ+h2-hgt1πP12ny1-y2 (16)

where

y0=1-Exφ1-π, y1=11-π-Exφ, y2=λP1+λEx1-θ1P2

2.2. Case 2: No machine failure occurrence during fabrication uptime

Figure 5 displays the level of perfect inventories in this case (i.e., tt ). Since no failure occurs, the inventory level goes up to H 1 when uptime ends and it reaches H 2 when rework is completed. Upon receipt of the outsourced items, the level of perfect inventories reaches H, before the beginning of product distribution time t .

Since no machine failure occurs, the safety stock remains unused throughout the cycle length T π. The following straightforward equations for this no failure occurrence case can be directly observed:

Tπ=t1π+t2π+t3π (17)

t1π=Q1-πP1=H1P1-d1 (18)

t2π=1-πQx1-θ1P2 (19)

t3π=Tπ-t1π+t2π (20)

H1=P1-d1t1π (21)

H2=H1+P2-d2t2π (22)

H=H2+πQ (23)

Equations (10) and (11) remain valid in this case, and the inventories in distribution time t and at the buyer side in T π can be computed using formulas shown in Eqs. (24) and (25) (Chiu et al., 2019b).

1n2i=1n-1iHt3π =1n2n(n-1)2Ht3π =n-12nHt3π (24)

12Ht3πn+H-λt3πTπ (25)

2.2.1. Total cost per cycle for Case 2

Total cost per cycle in no failure occurrence case, TC(t )2 comprises both the variable and fixed outsourcing and in-house fabrication costs, safety stock holding cost (see Fig. 8), rework and disposal costs, transportation costs (both variable and fixed costs), and total holding costs (including reworked, perfect items, and nonconforming items, and buyer’s stocks) during the entire cycle, as shown in Eq. (26).

TCt1π2=CππQ+Kπ+C1-πQ+K+h3λtrTπ+CTQ1-φx1-π+nK1+CRx1-πQ1-θ1+CSφx1-πQ+h1P2t2π2t2π+hH1+d1t1π2t1π+H1+H22t2π+n-12nHt3π+h 22Ht3πn+H-λt3πTπ (26)

Substitute equations (17) to (25) and (10) to (11) in Eq. (26), and use the expected value to cope with the randomness of x, the expected total system cost per cycle for case 2, E[TC(t )2] can be derived as follows:

ETCt1π2=Cππt1πP11-π+Kπ+Ct1πP1+K+CRExt1πP11-θ1+CSExφt1πP1+nK1+CTt1πP1y01-π+Ex2t1π2P121-θ12P2h11-θ1-h+h2t1π2P1221-πy0y2λ+t1π2P12y0h2-h2nλ1-πy1-y2+ht1π2P122λ1-πy021-π+λExφ1-π-πP1+λEx1-θ11-2πP2+h3gt1πP11-πy0 (27)

3. Solution processes to the problem

Because of the assumption of Poisson distributed failure rate β per year, the time to failure obeys an Exponential distribution with f(t) = βe-βt (i.e., the density function) and F(t) = (1 - e-βt) (i.e., the cumulative density function). Also, since the scrap rate φ is random, hence, the cycle length is not constant.The renewal reward theorem is employed to cope with the variable cycle length. Therefore, E[TCU(t)] can be computed as follows:

ETCUt1π=0t1πETCt1π1ftdt+ t1πETCt1π2ftdtE[Tπ] (28)

where E[ T π], E[T' π], and E[T π] represent the following:

ETπ= 0 t1π ET'π ftdt+ t1π  ETπ ftdt (29)

E[T'π]=Q1-φEx1-π+λtrλ=t1πP111-π-φEx+λtrλ (30)

E[Tπ]=Q1-φEx1-πλ=t1πP111-π-φExλ (31)

Substitute formulas (16), (27), and (29) in formula (28), along with extra efforts in derivations, one can obtain E[TCU(t )] as follows (for details please refer to Appendix A):

ETCUt1π=λy1+λg1-e-βt1πt1πP1W0t1π+W1t1π+W2+t1πW5-hge-βt1π+W3e-βt1πt1π+W4-W4e-βt1π (32)

The first and second derivatives of E[TCU(t )] are shown in equations (B-1) and (B-2) in Appendix B. Since the first term on the right-hand side (RHS) of Eq. (B-2) is positive, it follows that the E[TCU(t )] is convex if the second term on the RHS of Eq. (B-2) is also positive. That means if γ (t ) > t > 0 holds (see Eq. (B-3) for details).

Once Eq. (B-3) is verified to be true, we can solve the optimal t * by setting the first derivative of E[TCU(t )] = 0 (refer to Eq. (B-1)). Since the first term on the RHS of Eq. (B-1) is positive, we obtain the following:

hg+W4P1y1P1βe-βt1π+W5P1y1P1-λgβe-βt1πt1π2+W3P1-y1P1βe-βt1π+W5P12λg-2λge-βt1π+hg-W2P1λgβe-βt1πt1π-W0+W1P1y1P1+λgβe-βt1π+W3P1-λgβe-βt1π-y1P1e-βt1π-hg+W4P1λg-e-2βt1π+e-βt1π-W2+W4P1λge-βt1π-1=0 (33)

Let z 0, z 1, and z 2 represent the following:

z0=hg+W4P1y1P1βe-βt1π+W5P1y1P1-λgβe-βt1πz1=W3P1-y1P1βe-βt1π+W5P12λg-2λge-βt1π+hg-W2P1λgβe-βt1πz2=-W0+W1P1y1P1+λgβe-βt1π+W3P1-λgβe-βt1π-y1P1e-βt1π-hg+W4P1λg-e-2βt1π+e-βt1π-W2+W4P1λge-βt1π-1

Then, we can rearrange Eq. (33) as follows:

z0t1π2+z1t1π+z2=0 (34)

Apply the square roots solution, t π * can be found as follows:

t1π*=-z1±z12-4z0z22z0 (35)

As the cumulative density function of Exponential distribution F(t ) = (1 - e -βt1π ) is throughout for [0.1], so does its complement e -βt1π . Moreover, Eq. (33) can be rearranged as follows:

e-βt1π=-W5t1πP1y1P1-2λg+W0+W1P12y1-W2+W4P1λghg+W4P12y1βt1π2+-W3P12y1β-W5t1πP1λgβ+hg-W2P1λgβt1π-W2+W4P1λg-hg+W4P1λg1-e-βt1π-2W5t1πP1λg-W0+W1P1λgβ+W3P1-λgβ-y1P1 (36)

To solve the optimal t *, we start with letting e-βt = 0 and e-βt = 1, then compute Eq. (35) to find the bounds for t (i.e., t 1πU and t 1πL). Next step use present t 1πU and t 1πL to compute and obtain update values of e-βt and e-βt. Re-compute Eq. (35) using the current e-βt and e-βt to obtain the update bounds t 1πU and t 1πL. If (t 1πU = t 1πL) holds, then, t * is derived (i.e., t * = t 1πU = t 1πL); otherwise, repeat the above-mentioned steps, until it holds.

4. Numerical example

The following numerical example demonstrates the applicability of our obtained result. The assumed values of system variables in this example are shown in Table 1.

Table 1 Assumed of values of system variables 

β K1 Cπ λ C C1 β2 P1 Kπ CR K CS CT h2
1 90 2.8 4000 2.0 2.0 0.4 10000 60 1.0 200 0.3 0.01 1.6
π n θ1 M θ2 h3 β1 P2 x φ g h h1

First, we verify if E[TCU(t )] is convex (i.e., whether Eq. (B-3) holds). Since e -βt1π falls within the interval of [0, 1], let e -βt1π = 0 and e -βt1π = 1, and apply Eq. (35) to gain t 1πU = 0.2875 and t 1πL = 0.0909 initially. Then, use t 1πU and t 1πL to calculate e -βt1πU and e -βt1πL . Finally, apply Eq. (B-3) with the present values of e -βt1πL , e -βt1πU , t 1πL, and t 1πU to confirm that γ(t1πL) = 0.3103 > t1πL = 0.0909 > 0 and γ(t1πU) = 0.5320 > t1πU = 0.2875 > 0, respectively. Therefore, the convexity of E[TCU(t)] is assured for β = 1.0, and optimal t* exists. Additionally, a wider range of β values have been used to test for convexity of E[TCU(t)] to demonstrate the boarder applicability of the obtained result from this study (see Table 2)

Table 2 Verification of convexity of E[TCU(t)] against different βs

β γ(t1πU) t1πU γ(t1πL) t1πL
10 0.7927 0.2844 0.0467 0.0216
8 0.6141 0.2845 0.0573 0.0263
6 0.4998 0.2847 0.0744 0.0336
5 0.4621 0.2848 0.0874 0.0389
4 0.4370 0.2850 0.1060 0.0461
3 0.4268 0.2853 0.1346 0.0561
2 0.4415 0.2858 0.1851 0.0703
1 0.5320 0.2875 0.3103 0.0909
0.5 0.7277 0.2909 0.5215 0.1044
0.01 6.0228 0.5277 5.6043 0.1200

To solve the optimal t *, we start with letting e -βt1π = 0 and e -βt1π = 1 and apply Eq. (35) to gain the bounds for t (i.e., t 1πU = 0.2875 and t 1πL = 0.0909). Next, we repeatedly use resent t 1πU and t 1πL to compute and update values of e -βt1πU and e -βt1πL , and re-compute Eq. (35) using current e -βt1πU and e -βt1πL until t 1πU = t 1πL = t *. Table 3 exhibits the step-by-step results for searching t *. Therefore, the optimal uptime for this example t * = 0.1224 and E[TCU(t *)] = $12,542.25.

Table 3 Step-by-step results for searching t*. 

Step # t1πU e-βt1πU t1πL e-βt1πL t1πU - t1πL E[TCU(t1πU)] E[TCU(t1πL)]
- - 0 - 1 - - -
1 0.2875 0.7501 0.0909 0.9131 0.1966 $13,371.17 $12,637.28
2 0.1539 0.8573 0.1151 0.8913 0.0388 $12,598.72 $12,546.23
3 0.1292 0.8788 0.1207 0.8863 0.0085 $12,545.38 $12,542.44
4 0.1239 0.8835 0.1220 0.8851 0.0019 $12,542.41 $12,542.26
5 0.1227 0.8845 0.1223 0.8849 0.0004 $12,542.26 $12,542.25
6 0.1224 0.8848 0.1224 0.8848 0.0000 $12,542.25 $12,542.25

4.1. Impact of core system feature on the problem

The convexity of E[TCU(t )] and the initial bounds for t is exhibited in Figure 6.

Figure 6 The convexity of E[TCU(t)] and the initial bounds for t

Figure 7 illustrates the impact of variations in φ along with various x values on E[TCU(t*)]. It indicates that as both φ and x rise, E[TCU(t*)] increases noticeably.

Figure 7 The impact of variations in φ along with various x values on E[TCU(t*)]. 

The influence of changes in mean-time-to-breakdown 1/β on E[TCU(t*)] is displayed in Figure 8. It shows our optimal solution E[TCU(t*)] = $12,542 (for x = 0.2 and n = 3), it also specifies that as 1/β increases to over 0.17, E[TCU(t*)] begins to decline significantly; and as 1/β rises to extremely large (e.g., 1/β ≥ 100), E[TCU(t*)] = $11,962 (i.e., the same result as what is obtained from a problem without breakdown occurrence).

Figure 8 The influence of changes in 1/β on E[TCU(t*)] 

The impact of differences in the number of deliveries n (per cycle) on E[TCU(t *)] is depicted in Figure 9. It shows our optimal solution given n = 3, it also indicates that when n = 2 we have the minimal E[TCU(t *)], and as n increases, E[TCU(t *)] goes up significantly.

Figure 9 The impact of differences in the number of deliveries n on E[TCU(t*)]. 

The effect of variations in the outsourcing portion π on utilization is demonstrated in Figure 10. It shows that utilization noticeably decreases as π increases; and for π = 0.4 (as we assumed in our example), utilization declines from 47.72% to 28.11%.

Figure 10 The effect of variations in outsourcing portion π on utilization. 

The breakup of E[TCU(t *)] of our example is exhibited in Figure 11. It reveals the sum of outsourcing relevant setup and variable costs is 37.7%; total in-house relevant costs are 51.1% (including quality and breakdown related expenses), and supply chain relevant cost (including delivery and buyer’s holding costs) adds up to 11.2%.

Figure 11 The breakup of E[TCU(t*)]. 

The influence of changes in the number of deliveries n (per cycle) on the delivery and stock holding costs is illustrated in Figure 12. It reveals that as n increases, the fixed product distribution cost goes up significantly and in-house holding cost rises accordingly (the latter is simply due to a slow stock movement from the vendor to the buyer when n increases); on the contrary, the buyer holding cost drops accordingly.

Figure 12 The influence of changes in n on the delivery and stockholding costs. 

4.2. The joint impact of the core system features on the problem

The joint impact of differences in uniformly distributed nonconforming rate x and total scrap rate φ on the optimal decision variable t* is explored and illustrated in Fig. 13. It shows that t* increases significantly as both x and φ rise.

Fig. 13 The joint impact of differences in x and φ on the optimal t*. 

The combined influence of variations in the outsourcing portion of a batch π and mean time to breakdown 1/β on the optimal decision variable t* is studied and depicted in Fig. 14. It reveals that t* decreases enormously as π increases, especially when 1/β value is less than 0.17; and as 1/β increases to over 0.17 and π < 0.45, t* declines noticeably.

Figure 14 The combined influence of variations in π and 1/β on the optimal t1π*

The joint effect of changes in the meantime to breakdown 1/β and total scrap rate φ on E[TCU(t*)] is demonstrated in Fig. 15. It shows E[TCU(t*)] decreases considerably as 1/β increases and as φ goes up, E[TCU(t*)] increases slightly. The effect of 1/β on E[TCU(t*)] is more significant than that from φ.

Figure 15 The joint effect of changes in 1/β and φ on E[TCU(t*)]. 

4.3. Discussion and limitation

This study develops the inventory replenishment model based on a case where only one or no breakdowns occur during a production cycle. Table C-1 (see Appendix C) presents the Poisson probabilities results for a machine with different mean breakdown rates per year. It specifies that for a machine in good condition, or with an average of less than one breakdown occurrence per year, this study is appropriate, as there is over 99.31% chance that only one or no breakdowns will occur (refer to Table C-1).

Besides, for a machine in fair condition, or with an average of less than or equal to two breakdown occurrences per year, our model indicates that there is over 97.27% chance of one or no breakdowns occurring (see Table C-1). However, to explore the fabrication planning for a piece of equipment having a mean breakdown rate greater than five per year, the suitability of our model will fall below 80%, therefore, a different model must be developed for this specific condition.

5. Conclusions

This study aims to assist multinational/translational manufacturing firms with making the accurate decisions in their intra-supply chain environments, enable competitive strategies (including quality, low-cost, and timely delivery), and cope with the realities of limited capacity and unreliable equipment. Therefore, we examine a vendor-buyer coordinated system featuring batch fabrication, outsourcing, multiple deliveries, rework/disposal of defective items, and Poisson distributed breakdowns. Using the modeling, formulation, derivation, and optimization procedure along with a recursive algorithm, we in sequence obtain the problem’s cost function, justify its convexity, and find the problem’s optimal replenishment runtime. A numerical illustration is offered to show the applicability of the result that reveals various key system characteristics/capabilities, such as the distinct and joint influences of breakdowns (see Figs. 8, 11, 14, and 15), outsourcing (Figs. 10, 11, and 14), rework/scrap (Figs. 7, 11, 13, and 15), and delivery-frequency factors (Figs. 9, 11, and 12) on various system parameters, performance, and optimal runtime (Fig. 6). The methods proposed here and their results can facilitate managerial operations planning and strategic decision-making in an intra-supply chain setting in practice.

Acknowledgments

This project is sponsored by the Ministry of Science and Technology of Taiwan (Funding number: MOST 107-2221-E-324-015.

References

Arabi, M., Dehshiri, M.A., & Shokrgozar, M. (2018). Modeling transportation supply and demand forecasting using artificial intelligence parameters (Bayesian model). Journal of Applied Engineering Science, 16(1), 43-49. https://doi.org/10.5937/jaes16-12829 [ Links ]

Al-Bahkali, E.A., & Abbas, A.T. (2018). Failure analysis of vise jaw holders for hacksaw machine. Journal of King Saud University Engineering Sciences, 30(1), 68-77. https://doi.org/10.1016/j.jksues.2015.12.007 [ Links ]

Arun, P., Lincon, S.A., & Prabhakaran, N. (2019). An automated method for the analysis of bearing vibration based on spectrogram pattern matching. Journal of Applied Research and Technology, 17(2), 126-136. https://doi.org/10.22201/icat.16656423.2019.17.2.805 [ Links ]

Bolaños, R.I., Escobar, J.W., & Echeverri, M.G. (2018). A metaheuristic algorithm for the multi-depot vehicle routing problem with heterogeneous fleet. International Journal of Industrial Engineering Computations, 9(4), 461-478. https://doi.org/10.5267/j.ijiec.2017.11.005 [ Links ]

Chakraborty, T., Giri, B. C., & Chaudhuri, K. S. (2008). Production lot sizing with process deterioration and machine breakdown. European Journal of Operational Research, 185(2), 606-618. https://doi.org/10.1016/j.ejor.2007.01.011 [ Links ]

Chiu, Y-S.P., Liu, C-J., & Hwang, M-H. (2017). Optimal batch size considering partial outsourcing plan and rework. Jordan Journal of Mechanical and Industrial Engineering, 11(3), 195-200. [ Links ]

Chiu, Y-S.P., Chiu, V., Lin, H-D., & Chang, H-H. (2019a). Meeting multiproduct demand with a hybrid inventory replenishment system featuring quality reassurance. Operations Research Perspectives, 6. https://doi.org/10.1016/j.orp.2019.100112 [ Links ]

Chiu, S.W., Wu, C-S., & Tseng, C-T. (2019b) Incorporating an expedited rate, rework, and a multi-shipment policy into a multi-item stock refilling system. Operations Research Perspectives, 6. https://doi.org/10.1016/j.orp.2019.100115 [ Links ]

De Fontenay, C.C., & Gans, J.S. (2008). A bargaining perspective on strategic outsourcing and supply competition. Strategic Management Journal, 29(8), 819- 839. https://doi.org/10.1002/smj.697 [ Links ]

Dohi, T., Okamura, H., & Osaki, S. (2001). Optimal control of preventive maintenance schedule and safety stocks in an unreliable manufacturing environment. International Journal of Production Economics, 74(1-3), 147-155. https://doi.org/10.1016/S0925-5273(01)00121-9 [ Links ]

Farsijani, H., Nikabadi, M.S., & Ayough, A. (2012). A simulated annealing approach to optimize multi-products EPQ model with discrete delivery orders, imperfect production processes and service level constraint. World Applied Sciences Journal, 16(8), 1142-1157. [ Links ]

Ghalme, S., Mankar, A., & Bhalerao, Y. (2017). Integrated Taguchi-simulated annealing (SA) approach for analyzing wear behaviour of silicon nitride. Journal of Applied Research and Technology , 15(6), 624-632. https://doi.org/10.1016/j.jart.2017.08.003 [ Links ]

Goerler, A., & Voß, S. (2016). Dynamic lot-sizing with rework of defective items and minimum lot-size constraints. International Journal of Production Research, 54(8), 2284-2297. https://doi.org/10.1080/00207543.2015.1070970 [ Links ]

Groenevelt, H., Pintelon, L., & Seidmann, A. (1992). Production batching with machine breakdowns and safety stocks. Operations Research, 40(5), 959-971. https://doi.org/10.1287/opre.40.5.959 [ Links ]

Mohammadi, M. (2017). The tradeoff between outsourcing and using more factories in a distributed flow shop system. Economic Computation and Economic Cybernetics Studies and Research, 51(4), 279-295. [ Links ]

Montarelo, L.A., Glardon, R., & Zufferey, N. (2017). A global simulation-optimisation approach for inventory management in a decentralised supply chain. Supply Chain Forum: An International Journal, 18(2), 112-119. https://doi.org/10.1080/16258312.2017.1305255 [ Links ]

Morales, F., Franco, C., & Mendez-Giraldo, G. (2018). Dynamic inventory routing problem: Policies considering network disruptions. International Journal of Industrial Engineering Computations, 9(4), 523-534. https://doi.org/10.5267/j.ijiec.2017.11.001 [ Links ]

Nielsen, I.E., & Saha, S. (2018). Procurement planning in a multi-period supply chain: An epiphany. Operations Research Perspectives , 5, 383-398. https://doi.org/10.1016/j.orp.2018.11.003 [ Links ]

Paz, J.C., Granada-Echeverri, M., & Escobar, J.W. (2018). The multi-depot electric vehicle location routing problem with time windows. International Journal of Industrial Engineering Computations, 9(1), 123-136. [ Links ]

Puška, A., Kozarević, S., Stević, Ž., & Stovrag, J. (2018). A new way of applying interval fuzzy logic in group decision making for supplier selection. Economic Computation and Economic Cybernetics Studies and Research, 52(2), 217-234. https:/doi.org/10.24818/18423264/52.2.18.13 [ Links ]

Richter, K. (1996). The EOQ repair and waste disposal model with variable setup numbers. European Journal of Operational Research, 95(2), 313-324. https://doi.org/10.1016/0377-2217(95)00276-6 [ Links ]

Rosar, F. (2017). Strategic outsourcing and optimal procurement. International Journal of Industrial Organization, 50, 91-130. https://doi.org/10.1016/j.ijindorg.2016.11.001 [ Links ]

Saari, J., & Odelius, J. (2018). Detecting operation regimes using unsupervised clustering with infected group labelling to improve machine diagnostics and prognostics. Operations Research Perspectives, 5, 232-244. https://doi.org/10.1016/j.orp.2018.08.002 [ Links ]

Sarker, B.R., Jamal, A.M.M., & Mondal, S. (2008). Optimal batch sizing in a multi-stage production system with rework consideration. European Journal of Operational Research, 184(3), 915-929. https://doi.org/10.1016/j.ejor.2006.12.005 [ Links ]

Shakoor, M., Abu Jadayil, W., Jaber, N., & Jaber, S. (2017). Efficiency assessment in emergency department using lean thinking approach. Jordan Journal of Mechanical and Industrial Engineering, 11(2), 97-103. [ Links ]

Skowronski, K., & Benton, W.C., Jr. (2018). The Influence of Intellectual Property Rights on Poaching in Manufacturing Outsourcing. Production and Operations Management, 27(3), 531-552. https://doi.org/10.1111/poms.12813 [ Links ]

Souha, B., Soufien, G., & Mtibaa, A. (2018). Using system generator to design a hardware implementation of a fault-tolerant control of induction motor for electrical vehicle. Journal of Engineering Research, 6(2), 138-154. [ Links ]

Stažnik, A., Babić, D., & Bajor, I. (2017). Identification and analysis of risks in transport chains. Journal of Applied Engineering Science, 15(1), 61-70. https://doi.org/10.5937/jaes15-12179 [ Links ]

Swenseth, R.S., & Godfrey, R.M. (2002). Incorporating transportation costs into inventory replenishment decisions. International Journal of Production Economics, 77(2), 113-130. https://doi.org/10.1016/S0925-5273(01)00230-4 [ Links ]

Thomas, D.J., & Griffin, P.M. (1996). Coordinated supply chain management. European Journal of Operational Research, 94, 1-15. https://doi.org/10.1016/0377-2217(96)00098-7 [ Links ]

Vining, A., & Globerman, S. (1999). A conceptual framework for understanding the outsourcing decision. European Management Journal, 17(6), 645-654. https://doi.org/10.1016/S0263-2373(99)00055-9 [ Links ]

Vinod, B., & Solberg, J.J. (1984). Performance models for unreliable flexible manufacturing systems. Omega, 12(3), 299-308. https://doi.org/10.1016/0305-0483(84)90025-2 [ Links ]

Vujosevic, M., Makajic-Nikolic, D., & Pavlovic, P. (2017). A new approach to determination of the most critical multi-state components in multi-state systems. Journal of Applied Engineering Science , 15(4), 401-405. https://doi.org/10.5937/jaes15-15936 [ Links ]

Zahraee, S.M., Rohani, J.M., & Wong, K.Y. (2018). Application of computer simulation experiment and response surface methodology for productivity improvement in a continuous production line: Case study. Journal of King Saud University - Engineering Sciences, 30(3), 207-217. https://doi.org/10.1016/j.jksues.2018.04.003 [ Links ]

Zhao, X., Qian, C., Nakamura, S., & Nakagawa, T. (2018). A summary of replacement policies with number of failures. International Journal of Mathematical, Engineering and Management Sciences, 3(2), 136-150. https://doi.org/10.33889/ijmems.2018.3.2-011 [ Links ]

Appendix A

Detailed derivations for Eq. (32) are given below.

First, the integration results for the numerator and the denominator of Eq. (28) are shown in Eqs. (A-1) and (A-2), respectively.

0t1πETC(t1π)1ftdt+t1πETC(t1π)2ftdt=Kπ+K+nK1+t1πδ1+t1π2δ2+M1-e-βt1π+CTλg1-e-βt1π-C1λg1-e-βt1π+hP1g-t1πe-βt1π-1βe-βt1π+1β+h3λg21-e-βt1π+12h2λg21-e-βt1π+g2nh2-ht1πP1y1-y21-e-βt1π+hg2t1π-P1y1-y21-e-βt1π+g2h2+2h3t1πP1y1-y21-e-βt1π (A-1)

where

δ1=CππP11-π+CP1+CTy1P1+CRExP11-θ1+CSφExP1
δ2=Ex2P121-θ12P2h11-θ1-h+P12h2-hy02nλ1-πy1-y2+h22P121-πy0y2λ+h2λP12(1-π)y02(1-π)+λP1Exφ1-π-π+λE[x](1-θ1)(1-2π)P2ETπ=t1πP111-π-Exφλ+g1-e-βt1π (A-2)

With further derivation, one obtains E[TCU(t )] as follows:

ETCUt1π=λy1+λg1-e-βt1πt1πP1W0t1π+W1t1π+W2+t1πW5-hge-βt1π+W3e-βt1πt1π+W4-W4e-βt1π (32)

where

W0=KπP1+KP1+nK1P1

W1=MP1+CTλgP1+C1λgP1+h3λg2P1+12h2λg2P1+hgβ

W2=Cππ11-π+C+CTy1+CREx1-θ1+CSφEx

W3=-MP1-CTλgP1-C1λgP1-h3λg2P1-12h2λg2P1-hgβ

W4=hg2y01-π-λP1-λEx1-θ1P2+g2nh2-hy1-y2+g2h2+2h3y1+y2

W5=Ex2P11-θ12P2h11-θ1-h+P1y12nλh2-hy1-y2+h2P1y0y22λ1-π+h2λP1(1-π)y02(1-π)+λP1Exφ1-π-π+Exλ1-θ1P2[1-2π]

Appendix B

The first and second derivatives of E[TCU(t )] are shown in equations (B-1) and (B-2) below:

dETCUt1πdt1π=λy1t1πP1+λg1-e-βt1π2-W0+W1P1y1P1+λgβe-βt1π+W3P1-y1t1πP1βe-βt1π-λgβe-βt1π-y1P1e-βt1π+W5t1πP1y1t1πP1+2λg-2λge-βt1π-t1πλgβe-βt1π-hg+W4P1-y1t1π2P1βe-βt1π-λge-2βt1π-t1πλgβe-βt1π+λge-βt1π-W2+W4P1λgt1πβe-βt1π+e-βt1π-1 (B-1)

and

d2ETCUt1πdt1π2=λy1t1πP1+λg1-e-βt1π3W0+W1P12y12P12+y1t1πP1λgβ2e-βt1π+4y1t1πP1λgβ2e-βt1π+λ2g2β2e-2βt1π+λ2g2β2e-βt1π+W3P1e-βt1πy12t1πP12β2+2y12P12+2y12t1πP12β+2y1t1πP1λgβ2+2y1t1πP1λgβ+y1t1πP1λgβ2e-βt1π+2y1t1πP1λgβe-βt1π+λ2g2β2e-2βt1π+λ2g2β2+W5P1λgy1t1π3P1λgβ2e-2βt1π+2λge-2βt1π+t1π2λgβ2e-2βt1π-4λge-βt1π+4t1π3λgβ2e-2βt1π+t1π2λgβ2e-βt1π-4t1πλgβe-2βt1π+2λg-hg+W4P1e-2βt1πy12t1π3P12β2eβt1π+2y1t1π2P1λgβ2eβt1π-2y1t1πP1λgβeβt1π+y1t1π2P1λgβ2+2λ2g2β+4y1t1πP1λgβ+t1π2λ2g2β2eβt1π+2y1P1λg-2y1P1λgeβt1π-2λ2g2βeβt1π+t1πλ2g2β2+(W2+W4)P1λgt1πλgβ2e-2βt1π+t1πλgβ2e-βt1π+2λgβ2e-βt1π+y1t1π2P1β2e-βt1π+2y1P1β2e-βt1π+2y1t1πP1βe-βt1π-2y1P1 (B-2)

Since the first term on the right-hand side (RHS) of Eq. (B-2) is positive, it follows that the E[TCU(t )] is convex if the second term on the RHS of Eq. (B-2) is also positive. That means if the following γ(t ) > t > 0 holds.

γt1π=W0+W12y12P12+4y1P1λgβe-βt1π+λ2g2β2e-2βt1π+λ2g2β2e-βt1π+W3e-βt1π2y12P12+2y1P1λgβ+λ2g2β2+2y1P1λgβe-βt1π+λ2g2β2e-βt1π+W5λg2λge-2βt1π-4λge-βt1π+2λg-hg+W4e-2βt1π2y1P1λg-2y1P1λgeβt1π-2λ2g2βeβt1π+2λ2g2β+W2+W4λg2λgβe-2βt1π-2λgβe-βt1π+2y1P1e-βt1π-2y1P1-W0+W1y1P1λgβ2e-βt1π+W3e-βt1πy1y1t1πP12β2+2y1P12β+2P1λgβ2+P1λgβ2e-βt1π+W5λgy1t1π2P1β2e-βt1π+t1πλgβ2e-2βt1π+4λgβe-2βt1π+t1πλgβ2e-βt1π-4λgβe-βt1π-hg+W4e-2βt1πy12t1π2P12β2eβt1π+4y1P1λgβ+y1t1πP1λgβ2+λ2g2β2+2y1t1πP1λgβ2eβt1π-2y1P1λgβeβt1π+λ2g2β2eβt1π+W2+W4λgλgβ2e-2βt1π+λgβ2e-βt1π+y1t1πP1β2e-βt1π+2y1P1βe-βt1π>t1π>0 (B-3)

Appendix C

Table C-1: The probabilities of various Poisson 

β t* P(x=0) P(x=1) P(x≤1) P(x<1)
5.0 0.1644 43.95% 36.13% 80.09% 19.91%
4.0 0.1480 55.32% 32.75% 88.07% 11.93%
3.0 0.1356 66.59% 27.08% 93.67% 6.33%
2.0 0.1271 77.55% 19.71% 97.27% 2.73%
1.5 0.1243 82.99% 15.47% 98.46% 1.54%
1.0 0.1224 88.48% 10.83% 99.31% 0.69%
0.5 0.1214 94.11% 5.71% 99.82% 0.18%
0.01 0.1213 99.88% 0.12% 100.00% 0.00%

Escriba aquí la ecuación.

e-βt1x*(βt1x*)xx! (C-1)

Peer Review under the responsibility of Universidad Nacional Autónoma de México.

Received: July 12, 2019; Accepted: May 27, 2020; Published: August 31, 2020

Corresponding author. E-mail address: sungpc@cyut.edu.tw (Peng-Cheng Sung).

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License