The bullwhip effect in capacitated supply chains ... - Semantic Scholar

25 downloads 1994 Views 821KB Size Report
Dec 20, 2011 - Game''—a simulation based teaching tool to explain the economic .... ios. For ease of presentation, we approximate different phases of.
Int. J. Production Economics 136 (2012) 318–331

Contents lists available at SciVerse ScienceDirect

Int. J. Production Economics journal homepage: www.elsevier.com/locate/ijpe

The bullwhip effect in capacitated supply chains with consideration for product life-cycle aspects Bimal Nepal a,n, Alper Murat b, Ratna Babu Chinnam b a b

Department of Engineering Technology and Industrial Distribution, Texas A&M University, College Station, TX 77843, USA Department of Industrial and Systems Engineering, Wayne State University, Detroit, MI 48202, USA

a r t i c l e i n f o

a b s t r a c t

Article history: Received 7 March 2011 Accepted 10 December 2011 Available online 20 December 2011

This paper presents an analysis of the bullwhip effect and net-stock amplification in a three-echelon supply chain considering step-changes in the production rates during a product’s life-cycle demand. The analysis is focused around highly complex and engineered products (e.g., automobiles), that have relatively long production life-cycles and require significant capital investment in manufacturing. Using a simulation approach, we analyze three stages of the product life-cycle including low volumes during product introduction, peak demand, and eventual decline toward the end of the life-cycle. Parts of the simulation model have been adopted by a major North-American automotive OEM as part of a scenario analysis tool for strategic supply network design and analysis. The simulation results show that performance of a system as a whole deteriorates when there is a step-change in the life-cycle demand. While restriction in production capacity does not significantly impact the bullwhip effect, it increases the net stock amplification significantly for the supply chain setting under consideration. Furthermore, a number of important managerial insights are presented based on sensitivity analysis of interaction effect of capacity constraints with other supply chain parameters. & 2011 Elsevier B.V. All rights reserved.

Keywords: The bullwhip effect Life-cycle demand Capacity constraints Net-stock amplification

1. Introduction The bullwhip effect is one of the most widely investigated phenomena in the modern day supply chain management research. It is the tendency to see an increase in variability in the replenishment orders with respect to true demand due to distortion in the demand information as we move upstream in the supply chain. While Lee et al. (1997) first introduced the term ‘‘bullwhip effect’’ to explain this phenomenon, it was first described by Forrester (1961) to demonstrate the demand and variance amplification in an industrial system. His idea has been studied further and illustrated through the ‘‘Beer Distribution Game’’—a simulation based teaching tool to explain the economic dynamics of stock management problem (Sterman, 1989). Lee et al. (1997) identified the following four reasons for the bullwhip problem: demand signal processing, the rationing game, order batching, and price variations. Since then, there have been a significant number of studies on this problem with respect to all the major causes of the bullwhip effect (Chen et al., 2000a; Dejonckheere et al., 2003; Disney and Towill, 2003; Moyaux et al., 2007; Boute and Lambrecht, 2007). Recently, third-party warehousing has also

n

Corresponding author. Tel.: þ1 979 845 2230; fax: þ1 979 845 4980. E-mail addresses: [email protected] (B. Nepal), [email protected] (A. Murat), [email protected] (R. Babu Chinnam). 0925-5273/$ - see front matter & 2011 Elsevier B.V. All rights reserved. doi:10.1016/j.ijpe.2011.12.018

been cited as one of the causes of the bullwhip effect (Duc et al., 2010). Classical inventory management models for multi-echelon supply chains require that the product demand process be fairly smooth in order to make it mathematically tractable (Williams, 1982). In comparison, simulation models are well suited to study complex supply chains with non-smooth demand process (e.g., lumpy demand during life-cycle) and allows for transient analysis. It is therefore widely used in the bullwhip effect analysis as well (Disney et al., 2004a; Wanphanich et al., 2010; Coppini et al., 2010). While simulation based games like ‘‘Beer Game’’ have helped researchers and practitioners understand dynamics of order and inventory fluctuations in a supply chain (SC) system, very few examples can be found that incorporate the life-cycle demand aspects into those dynamics (Disney et al., 2004b; Reddy and Rajendran, 2005). It is a well-known fact that every product has its own life-cycle demand curve—slow market growth at introduction, rapid growth during peak, and sluggish demand during saturation (maturity) or decline phase (Mahajan and Muller, 1979). Kaipa et al. (2006) discuss the nervousness of demand planning and its impact on bullwhip in an electronic SC in the face of changing demand at various life-cycle phases. Hoberg et al. (2007) also argue that the conventional approach of using a lower smoothing constant in forecasting will take a significantly long time to detect step-changes in demand. While the notion of a product life-cycle is not new and somewhat witnessed in all industries (e.g., see Kaipa et al., 2006; Berry

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

and Towill, 1982), the ‘‘signature’’ of the production life-cycle is unique to industries that produce complex engineered products (e.g., automobiles). During the first few months of introduction, the vehicle is typically produced at low volume to address any production quality and supply issues. Given the complexity of the assembly and other production facilities, it is not practical to change the production volume continuously due to the need to ‘‘balance’’ the assembly line and the supply chain. Instead, the production volume typically undergoes step changes at distinct epochs as the demand picks up for the product in the market place, e.g., in the form of additional production shifts and step changes in volume per shift (Kisiel, 2008). For example, most original equipment manufacturer (OEM) automotive products in the North-American market undergo a four to five year life-cycle (with some product ‘‘refreshing’’ every year), with demand typically waning after couple of years due to introduction of more competitive products in the market place with better functions, features, and option content. These life-cycles however tend to be different from industry to industry, with complex engineered product typically experiencing longer life-cycles due to complexity, costs, and risks associated with product development and launch. Another characteristic, at least typical of the North-American automotive OEMs, is that they predominantly operate in a ‘‘build to stock’’ production mode rather than a ‘‘build to order’’ mode. This is also attributable to complexity of the product, production facilities, challenges associated with coordination of the supply chains to support the ‘‘take rates’’ for the different options/features, and the long order-todelivery lead-times. Accordingly, the companies rely heavily on dynamic pricing strategies (in the form of dealer or customer incentives) and marketing to influence demand while adjusting supply in the long-run to match demand. Hence, this study focuses more on modeling and control of the bullwhip effect in the supply networks as a function of fluctuations in the OEM final-assembly (FA) line production volume rather than end customer ‘‘demand’’. Fig. 1 shows a typical production pattern life-cycle for an automobile. In such a case, adopting a ‘‘one size fits all’’ production and inventory management policy results in a chaotic situation especially in multi-echelon SC settings. It is necessary to investigate how inventory and order variances propagate as a product passes through different phases of the life-cycle. Capturing such transient trajectories of different SC performance measures will be very helpful especially in designing or configuring a SC network for future products, such as the automotive industry example addressed in this paper. The unique contributions of this paper are several-fold. First, we extend the work of Chen et al. (2000b) and present analytical expressions of the bullwhip effect and net-stock amplification for OEM using a different sequence of events in a replenishment period. These results confirm the bullwhip effect results developed in the literature for slightly different ordering policy. The

Production

Peak Introduction

2nd Shift Added

Ramp-up (Step-Shifts)

End-of-Life Decline (Step-Shifts) Life Cycle

Fig. 1. Typical production pattern life-cycle for an automobile with introduction, peak, and end-of-life stages.

319

analysis is then extended via simulation to a three-echelon SC consisting of OEM, Tier 1, and Tier 2 suppliers. Through simulation based models, we investigate the transient and steady-state impact of step-changes in OEM production volume on the bullwhip effect and net stock amplification. These analyses are performed for both with and without capacity constraint scenarios. For ease of presentation, we approximate different phases of the product life-cycle with three stationary phases (labeled, ‘‘introduction’’, ‘‘peak’’, and ‘‘end-of-life’’ stages) corresponding to the three stages of Fig. 1. Various sensitivity analyses are performed to explore the impact of interaction between capacity constraint and other SC parameters on the bullwhip effect and inventory variance (or net-stock amplification). Several insights of managerial importance are drawn based on the sensitivity analyses. The remainder of the paper is organized as follows. A brief review of related literature on the bullwhip effect analysis is presented in Section 2. Section 3 describes the model and the undertaken SC policies. In Section 4, we present the analytical expression for inventory variance and bullwhip effect (from the literature) for OEM. Section 5 describes the simulation model and sensitivity analysis results for a three-echelon SC with and without capacity constraints. Finally, conclusions and directions for future research are outlined in Section 6.

2. Review of related literature The bullwhip effect related research in supply chains has a long tradition which can be broadly divided into three streams. The first stream of research focuses on determining the impact of forecasting techniques employed by SC players on the bullwhip effect. Chen et al. (2000a) have analytically derived a quantitative measure of the bullwhip effect and demonstrated via simulation the impact of forecasting, lead-time and information. The authors used simple moving average (MA) forecasting techniques to determine the lower bounds of the bullwhip effect in simple single and multi-echelon supply chains. They later extended their work by studying impact of other forecasting techniques, in particular exponential smoothing (ES) under both auto-correlated demand and demand with linear trend (Chen et al., 2000b). Zhang (2004) compared the bullwhip effect under three different forecasting methods for a simple inventory system with a first-order autoregressive, AR(1), demand process. Zhang (2004) found the minimum mean square error (MMSE) to be the optimal forecasting method for stable AR(1) demand processes. In comparison, the MA or ES methods were more flexible and adapted better to the demand shifts over time. Disney et al. (2006) have developed exact methods to quantify the bullwhip effect for different demand processes including AR, MA, and auto-regressive moving average (ARMA). Holweg et al. (2005) discussed the challenges in implementing a Build-to-Order (BTO) production system given current SC network structures, in light of the bullwhip effect and scheduling issues especially in automotive industry. Hosada and Disney (2006) presented a control theoretic approach to quantify the bullwhip effect under the MMSE forecasting. The authors also derived an analytical expression for another important performance measure of supply chains, the ‘‘net-stock amplification’’. It is defined as the ratio of the net-stock variance over the variance of demand, which, just like the bullwhip effect, gets worse as we move up the chain. In addition to AR(1), attempts have been made to determine an upper bound for the bullwhip effect for two stage SC with first and second order autoregressive demand process (Luong, 2007; Luong and Phien, 2007). Duc et al. (2008) analyzed on a two-stage SC including one supplier and one retailer with a mixed autoregressive-moving average model, ARMA(1, 1) for

320

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

demand forecasting and base-stock inventory policy. They found that the bullwhip effect depended (although not always) on the autoregressive coefficients. For more comprehensive comparison of the bullwhip effects due to different forecasting techniques in a single stage supply chain, we refer our readers to Wang et al. (2010). The other two streams of research in the bullwhip effect analysis include an examination of impact of operations management parameters (such as ordering policy, inventory management policy, and production variation and batching) and supply chain dynamics (like information sharing) on the bullwhip effect. Dejonckheere et al. (2003) demonstrated that the bullwhip effect is imminent under an order-up-to (OUT) policy regardless of the forecasting technique. Disney and Towill (2003) have presented analytical expressions for bullwhip and net-stock (NS) amplification under the OUT policy and MMSE forecasting technique for a single echelon problem. Hoberg et al. (2007) have derived a closedform expression for a single stage supply chain to describe the transient behavior of order and inventory fluctuations when the demand process is non-stationary. A number of authors have investigated and quantified the value of information sharing to reduce the bullwhip effect in a two stage supply chain (Lee and Wang, 2000; Lee et al., 2000; Raghunathan, 2001). The prior research shows that the bullwhip effect cannot be completely eliminated even with complete information sharing. In fact, the lead-time reduction is more beneficial than the sharing of information in terms of reducing the bullwhip effect (Agrawal et al., 2009). Whereas much has been done in terms of understanding and mitigating the bullwhip effect in a two stage supply chain, prior models largely assume a fairly smooth process for OEM production volume. Hence, these models do not adequately capture the impact of step-changes in OEM production volume on the performance of the overall SC system. In this paper, an attempt has been made to model the bullwhip effect and NS amplification in a three-stage automotive SC system considering capacity constraint and step-changes in OEM production volume due to product life-cycle considerations. We also extend the existing work of Chen et al. (2000b) by developing analytical expressions for both the bullwhip effect and the NS amplification for single stage supply chain with different ordering sequence in a replenishment cycle. While the analysis is focused on automotive supply chains, the results can be applicable to other SC systems that support a fairly complex product and undergo discrete step changes in production volume over the product life-cycle phases.

3. Modeling the three-echelon SC In this paper, we consider a three-echelon SC consisting of a single OEM, a single Tier 1 (T1) supplier and a single Tier 2 (T2) supplier. The sequence of events at any echelon in any given period t is as follows:

fulfillment. We have chosen this ordering since in the automotive industry, OEMs first observe the demand (e.g., dealers place their orders with the OEM) and then the OEM makes the production and ordering decisions with their Tier-1 suppliers. There is a fixed lead-time (L) between the order placement and its receipt at each echelon, e.g., an order placed at the end of period t is received at the beginning of period t þL. We adopt a periodic review inventory policy given its popularity in the industry. Prior research also shows that it is optimal when there is no fixed ordering cost, and inventory holding cost and backlog cost are proportional to netstock or shortage quantity (Zipkin, 2000). We assume that raw material inventory level at each echelon is reviewed in every period. Hence, the lead-time L is equal to Tþ1, where T is the physical lead-time including production time plus a single review period. 3.1. Demand process We assume that the final assembly rate realized by the original equipment manufacturer (OEM), (a.k.a. assembly facility) is an AR(1) process.1 There are some recent studies that use AR(1) in modeling demand process specific to the automotive industry (Cachon et al., 2007; Cachon and Olivares, 2010). This setting is equivalent to the case of a retailer facing random customer demand and the AR(1) has been successfully employed in earlier studies (Chen et al., 2000a; Hosada and Disney, 2006; Lee et al., 2000). The demand observed by the OEM at period t (Dt) is a random variable of the form: Dt ¼ d þ rDt1 þ et

ð1Þ

where d is a non-negative constant, Dt  1 is the demand faced by the OEM at time period t  1, r is the correlation parameter with 1o r o1., and et represents the error terms that are i.i.d. with mean 0 and variance se 2 . Further, without loss of generality, we assume that each OEM assembled unit contains a single unit or sub-system supplied by the Tier 1 supplier and the bill of materials (BOM) ratio is one-to-one between the tiers, driving the demand for the suppliers. 3.2. Forecasting method We assume that the OEM and its suppliers employ an exponential smoothing forecasting technique to estimate future demand. This is the most common forecasting technique used in practice (Chen et al., 2000b). This is the most widely used forecasting technique in practice primarily due to its ease of implementation. Also, the policy is very popular in the automotive industry. Although the suppliers often have a general sense for the OEM’s planned production during the planning horizon (e.g., forecasts shared in advance through MRP and/or EDI systems), they do not know the true demand process given the randomness in actual inventory consumption at the OEM. The OEM’s exponential smoothing forecast for supply consumption ^ t Þ is given by during period t ðD

1. Replenishment: shipment from the upstream echelon is received and placed in raw material inventory. 2. Demand fulfillment: demand from downstream echelon, Dt, is observed and filled from finished product inventory (if inventory is available) or backlogged. 3. Production: production decision is made based on the finished goods inventory and depends upon observed demand, available capacity and raw material inventory on hand. 4. Ordering: overall inventory position is reviewed and a replenishment order, Ot, is placed to the upstream echelon for raw material inventory.

where 0o a1 o1 and 0o n1 o1 are OEM’s exponential smoothing constants for mean and variance, respectively. The OEM’s leadtime supply consumption beginning with period t is DLt,OEM ¼ Stt0þ¼L1t þ 1 Dt0 ,OEM , where L1 is the lead-time for OEM’s replenishment

The sequence of events above is different from those used in Chen et al. (2000b) in that ordering precedes the demand

1 The non-stationary aspects of supply consumption due to product life-cycle aspects are discussed later.

^ t1 ^ t ¼ a1 Dt1 þ ð1a1 ÞD D

ð2Þ

s^ 2t,OEM ¼ n1 ðD^ t Dt1 Þ2 þ ð1n1 Þs^ 2t1,OEM

ð3Þ

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

by the Tier 1 supplier. Accordingly, the mean and the standard deviation of OEM’s lead-time supply consumption forecast in pffiffiffiffiffi ^ t þ 1 and s ^L ^ ^ D period t are D ¼ L ¼ L s (see Boute 1 L,t,OEM 1 t þ 1,OEM t,OEM and Lambrecht, 2007). Similar expressions for the mean and the standard deviation of T1 and T2 lead-time forecasts at period t ^ L and and calculated similar to that of OEM’s are denoted as D t,j s^ L,t,j for j¼T1,T2. Note that we are differentiating between the parameters for OEM, T1 and T2 through subscripts {OEM,T1,T2} or {1,2,3}, respectively, e.g., L1 is the lead-time for OEM. 3.3. Production

PRt,OEM ¼ minfct,OEM ,Dt ,NSt,OEM g

is returned without any penalty (see Chen et al., 2000a,b; Lee et al., 2000). kj is a constant determined based on the predetermined service level at stage j. We outline the ordering policy for each player separately. 3.4.1. Ordering policy for OEM According to (7), the order placed by the OEM during the period t, net-stock and work-in-process (WIP) levels are given by Ot,OEM ¼ Dt,OEM þ ðSt,OEM St1,OEM Þ L

In this paper, we consider capacity constraint at each echelon. As presented in Lee et al. (2000), we also assume that an upstream echelon ships the required order placed by its immediate downstream customer in a given lead-time. In case the order cannot be fulfilled due to the limited capacity or insufficient on-hand inventory of raw materials, the upstream echelon will resort to alternative sources to fulfill the customer order. The alternative sources may include running overtime production or adding a shift temporarily etc. The cost of ‘‘alternative source’’ has been represented by the backlog (or expediting) cost in the model. Considering the available capacity at time t (ct,OEM), on-hand inventory of modules (i.e., OEM’s raw material), and the target supply consumption quantity (equivalent to target production), the OEM’s actual production (PR) for period t is determined as ð4Þ

where, ct,OEM is the OEM’s production capacity at time t and NSt,OEM is the net-stock of modules available with OEM at the beginning of period t and is calculated as NSt,OEM ¼ |t,OEM þ NSt1,OEM . Here, |t,OEM . the quantity of the modules received from T1 supplier at the beginning of period t, which equals the order placed by the OEM at the end of period t  L1, i.e., OtL1 ,OEM . The inventory system at the suppliers is similar to that of the OEM with backorders (Lee et al., 2000). Similarly, the production rates in period t for T1 and T2 suppliers are given by PRt,T1 ¼ minfct,T1 ,Ot,OEM ,NSt,T1 g

ð5Þ

PRt,T2 ¼ minfct,T2 ,Ot,T1 ,NSt,T2 g

ð6Þ

where NSt,T1 ¼ ft,T1 þNSt1,T1 is on-hand inventory of components (T1’s raw material) with T1 at the beginning of period t, and Ot,OEM is the order placed by OEM to T1 during period t. Likewise, NSt,T2 ¼ ft,T2 þ NSt1,T2 is on-hand inventory of raw material with T2 at the beginning of period t, and Ot,T1 is the order placed by T1 to T2 during period t. The ct,T1 and ct,T2 are the production capacities of T1 and T2 at time t, respectively. Note that PR is a non-negative quantity. In other words, a negative production quantity in a period, t, refers to zero production in that period. 3.4. Ordering policy We assume that all three echelons adopt an adaptive OUT policy for their inventory management system. This policy minimizes the total holding and shortage cost over the long run (Kahn, 1987; Chen et al., 2000a). In the OUT policy, the general expression for order quantity in period t at stage j (j¼OEM,T1,T2) is expressed as Ot,j ¼ Dt,j þ ðSt,j St1,j Þ

ð7Þ

^ L þkj s ^ L,j,t St,j ¼ D t,j

ð8Þ

where, St,j is the base stock or OUT level at time period t and stage j, ^ L,j,t is an estimate of the standard deviation of demand over and s the lead-time at stage j in time t. Further, it may be noted that this policy allows Ot,j to be negative, which means the excess inventory

321

ð9Þ

^ ^ St,OEM ¼ D t,OEM þ kOEM sL,t,OEM

ð10Þ

NSt,OEM ¼ NSt1,OEM þ OtL OEM PRt,OEM

ð11Þ

WIPt,OEM ¼ WIPt1,OEM þOt1,OEM ft,OEM

ð12Þ

where, NSt,OEM and WIPt,OEM are the net-stock and WIP inventory of modules with the OEM at the beginning of period t. Using similar logic, we can write the ordering policy equations for T1 and T2 as follows. 3.4.2. Ordering policy for Tier 1 supplier Ot,T1 ¼ Ot,OEM þðSt,T1 St1,T1 Þ L

ð13Þ

^ ^ St,T1 ¼ D t,T1 þ kT1 sL,t,T1

ð14Þ

NSt,T1 ¼ NSt1,T1 þ OtL ,T1 PRt,T1

ð15Þ

WIPt,T1 ¼ WIPt1,T1 þ Ot1,T1 ft,T1

ð16Þ

3.4.3. Ordering policy for Tier 2 supplier Ot,T2 ¼ Ot,T1 þðSt,T2 St1,T2 Þ

ð17Þ

^ L þ kT2 s ^ L,t,T2 St,T2 ¼ D t,T2

ð18Þ

NSt,T2 ¼ NSt1,T2 þ OtL ,T2 PRt,T2

ð19Þ

WIPt,T2 ¼ WIPt1,T2 þ Ot1,T2 ft,T2

ð20Þ

4. The bullwhip and NS amplification effects under exponential smoothing forecasting In this section we consider only the OEM echelon of the SC and derive analytical expressions for the bullwhip effect and inventory amplification in the absence of capacity constraints. Since the closed-form analytical expressions for the capacitated case is not available, we analyze the impact of capacity on the bullwhip effect and net-stock amplification through empirical simulation study in the next section. 4.1. The bullwhip effect Chen et al. (2000b) consider an OUT policy in which, at every period, the retailer first places an order based on the inventory level and then observes and fulfills the demand. In contrast, we herein consider that first the demand is observed and fulfilled, and then the order is placed. Since the ordering policy is different in this paper, we re-visit the analysis of the bullwhip effect by accounting for this difference. Let B1(a1, L1) be the bullwhip effect for the OEM with exponential smoothing constant a1 and lead-time L1. We first derive the expression for the variance of the lead-time supply

322

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331 2

^ L,t,OEM . For this we note consumption forecast error for the OEM s ^ 2 L,t,OEM depends on the continuthat, as in Chen et al. (2000b), s ously updated estimate for the variance of the forecast error using (3). Alwan et al. (2003) note that these estimations based on ‘‘large samples’’ lead to asymptotically unbiased and consistent estimates and demonstrated through Monte Carlo simulations that the results using estimated from (3) are in very close approximation with the asymptotic case. Hence, we consider ^ 2 L,t,OEM as follows: the asymptotic case and express the s L

s^ 2 L,t,OEM ¼ VðDLt,OEM D^ t,OEM Þ L

L

L ^ ^ ¼ VðDLt,OEM Þ þ VðD t,OEM Þ2 CovðDt,OEM , Dt,OEM Þ:

ð21Þ

Using DLt,OEM ¼ Stt0þ¼L1t þ 1 Dt0 ,OEM , we can show the equivalence of the first term as follows:   VðDt,OEM Þ 2rð1rL1 Þ L1 ð1þ rÞ ð22Þ VðDLt,OEM Þ ¼ ð1rÞ ð1rÞ where VðDt,OEM Þ ¼ s2e =ð1-rÞ. The second term can be expressed ^L ^ using D t,OEM ¼ L1 Dt þ 1 and the asymptotic behavior of the forecasting system from Chen et al. (2000b, Eq. (4)), as follows: 0 1 1 X ^L ^ t,OEM Þ ¼ ðL1 Þ2 V @a1 VðD Þ ¼ VðL1 D ð1a1 Þj1 Dtj þ 1,OEM A t,OEM

0 ¼ ðL1 a1 Þ2 @

j¼1 1 X

ð1a1 Þ2ðj1Þ VðDtj þ 1,OEM Þ

j¼1

þ2

1 X 1 X

1 j1

ð1a1 Þ

i1

ð1a1 Þ

CovðDtj þ 1,OEM ,Dti þ 1,OEM ÞA

j ¼ 1 i4j

0 ¼ ðL1 a1 Þ VðDt,OEM Þ@ 2

1 X

ð1a1 Þ2ðj1Þ

j¼1

þ2

1 X 1 X

j þ i2

ð1a1 Þ

1 i1 A

r

Note that, by multiplying each side of Eq. (2) and using ^L ^ D t,OEM ¼ L1 Dt þ 1,OEM , we obtain the following: ^L ^ D t,OEM ¼ L1 a1 Dt,OEM þ L1 ð1a1 ÞðDt,OEM Þ ^ t,OEM þ L1 D ^ t,OEM ¼ L1 a1 Dt,OEM L1 a1 D ^L ^ t,OEM þ D ¼ L1 a1 Dt,OEM L1 a1 D t1,OEM and ^L ^ ^L D t,OEM Dt1,OEM ¼ L1 a1 ðDt,OEM Dt,OEM Þ:

ð25Þ

When we substitute (25) into (24) and express the variance ^L ^L ^ Ot,OEM ¼ Dt,OEM þ ðD t,OEM Dt1,OEM Þ ¼ Dt,OEM þ L1 a1 ðDt,OEM Dt,OEM Þ ^ t,OEM Þ VðOt,OEM Þ ¼ VðDt,OEM Þ þðL1 a1 Þ2 VðDt,OEM D ^ t,OEM Þ: þ 2ðL1 a1 ÞCovðDt,OEM ,Dt,OEM D Following the steps of the third result in Zhang (2004) for deriving expressions in VðOt,OEM Þ, we obtain the bullwhip effect expression as ( ! ) 2a2 L2 VðOt,OEM Þ 1r ¼ 1 þ 2a1 L1 þ 1 1 ð26Þ VðDt,OEM Þ 2a1 1b1 r where b1 ¼1 a1. We note that while the bullwhip effect expression (26) is identical to that derived in Alwan et al. (2003) and Zhang (2004) and in agreement with Chen et al. (2000b), it is based on the different sequence of events during a replenishment period. Note ^ 2 L,t,OEM is not constant (e.g., due to finite that in case where error s sample), then the right-hand-side of (26) will be a lower bound on the bullwhip effect. The above formula pertains to a situation without capacity constraint, as well as smooth supply consumption process, that is, without considering lifecycle demand. 4.2. Net-stock amplification effect

:

j ¼ 1 i4j

^L Note that the above equation VðD t,OEM Þ does not depend on t, L which was also the case for VðDt,OEM Þ. Similar expression can be ^L derived for CovðDLt,OEM , D t,OEM Þ as follows: ^L CovðDLt,OEM , D t,OEM Þ 0 1 tX þ L1 1 X j1 Dt0 ,OEM ,L1 a1 ð1a1 Þ Dtj þ 1,OEM A ¼ Cov@ j¼1 t’ ¼ t þ 1 0 1 tX þ L1 1 X 2 j1 ¼ ðL1 a1 Þ ð1a1 Þ Cov@ Dt0 ,OEM ,Dtj þ 1,OEM A j¼1 t’ ¼ t þ 1 1 X  ¼ ðL1 a1 Þ2 ð1a1 Þj1 CovðDt þ 1,OEM ,Dtj þ 1,OEM Þ

Lemma 1. The inventory variance amplification ratio experienced by the OEM under OUT policy, AR(1) demand (in other words, supply consumption) process and exponential smoothing forecasting technique is given by the expression ð1rL1 Þð2a1 L1 þ a1 L1 rÞ VarðNSt,OEM Þ L1 ð1 þ rÞ ¼ r VarðDt,OEM Þ ð1rÞ ð1rÞ2 ð1b1 rÞ

j¼1

 þ    þCovðDt þ L1 ,OEM ,Dtj þ 1,OEM Þ 1 X ¼ ðL1 a1 Þ2 VðDt,OEM Þ ð1a1 Þj1 ½rj þ rj þ 1 þ    þ rj þ L1 1 :

We derive the expressions for net-stock amplification without capacity constraint under exponential smoothing forecasting. Hosada and Disney (2006) present a control theoretic approach to derive the expression for net-stock amplification for a three echelon SC with MMSE forecasting. While there are a number of prior researches that have developed the expressions for the bullwhip effect considering exponential smoothing constant, we are not aware of any prior work on the expressions for variance amplification under exponential smoothing forecasting.

þ

ð23Þ

j¼1

a1 L1 2 ð1þ b1 rÞ ð1 þ b1 Þð1b1 rÞ

ð27Þ

L

^ Note that CovðDLt,OEM , D t,OEM Þ is also independent of time as the ^ 2 L,t,OEM . Hence, the standard deviation of leadfirst two terms in s ^ 2 L,t,OEM remains constant under the asymptotic time forecast error s ^ 2 L,t,OEM ¼ s ^ 2 L,t1,OEM . Using this observation and Eqs. analysis, e.g., s (9) and (10), one can obtain Ot,OEM ¼ Dt,OEM þðSt,OEM St1,OEM Þ L

L

^ ^ ^ ^ ¼D t,OEM þ kOEM sL,t,OEM ðDt1,OEM þkOEM sL,t1,OEM Þ þDt,OEM ^L ^L ¼ Dt,OEM þ ðD t,OEM Dt1,OEM Þ þ kOEM ^L ^L ¼ Dt,OEM þ ðD t,OEM Dt1,OEM Þ:

s^ L,t,OEM kOEM s^ L,t1,OEM ð24Þ

Proof. Let NSt(L1,a1) be the net-stock of module inventories with the OEM in time period t with exponential smoothing constant a1 and lead-time L1. It has been shown that, for OUT policy, the variance of net inventory levels and the variance of forecast errors over lead-time are equal (see Hosada and Disney, 2006). Therefore ^ L1 ^ 2OEM,L1,t ¼ VarðDL1 VarðNSt,OEM Þ ¼ s t,OEM Dt,OEM Þ 1 ^ L1 Þ2CovðDL1 , D ^ L1 Þ Þ þ VarðD ¼ VarðDLt,OEM t,OEM t,OEM t,OEM

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

where, DL1 t,OEM is lead-time demand for the OEM at time period t. Furthermore, from Zhang (2004), we obtain h i ð1 þ rÞ 2rð1rL1 Þ 1 VarðDLt,OEM Þ ¼ L1 1 r  ð1rÞ2 VarðDt,OEM Þ ^ L1 Þ ¼ VarðD t,OEM

  ða1 L1 Þ2

1 a1 ð2a1 Þ



 1þ ð1a1 Þr VarðDt,OEM Þ 1ð1a1 Þr

and,    rrL1 þ 1 1 ^ L1 Þ ¼ ða1 L1 Þ VarðDt,OEM Þ: CovðDLt,OEM ,D t,OEM ð1rÞð1b1 rÞ

" VarðNSt,OEM Þ L1 ð1 þ rÞ 2rð1rL1 Þ ¼  3 VarðDt,OEM Þ 1r ð1rÞ2    1 1 þ b1 rÞ þða1 L1 Þ2 a1 ð2a1 Þ 1b1 rÞ   rrL1 þ 1 ða1 L1 Þ ð1rÞð1b1 rÞ ¼

ð1rL1 Þð2a1 L1 þ a1 L1 rÞ L1 ð1þ rÞ r ð1rÞ ð1rÞ2 ð1b1 rÞ þ

a1 L1 2 ð1 þ b1 rÞ ð1þ b1 Þð1b1 rÞ

Scenario 4: step-change in supply consumption rate due lifecycle demand and with capacity constraints. Lastly, Eqs. (29) and (30) depict the formulae for the bullwhip effect and net-stock amplifications at a any node. Bullwhip effect ¼

Order variance at an echelon : Variance of demand at OEM

NS amplification ¼

Putting them all together, we get the net-stock amplification ratio as " L1 ð1þ rÞ 2rð1rL1 Þ VarðNSt,OEM Þ ¼  1r ð1rÞ2    1 1 þ b1 r Þ þða1 L1 Þ2 a1 ð2a1 Þ 1b1 rÞ   rrL1 þ 1 ða1 L1 Þ VarðDt,OEM Þ ð28Þ ð1rÞð1b1 rÞ

&

323

Net stock variance at an echelon : Variance of demand at OEM

ð29Þ

ð30Þ

5.1. Supply chain settings and baseline parameters for the simulation model The initial values for SC parameters have been drawn from the industry interviews. However, the actual company data are disguised to protect the confidentiality of participants. The bill of materials (BOM) ratio is one-on-one between the OEM and the two tiers. The production rate or supply consumption process at the OEM is assumed to be first-order autoregressive with a mean of 5400 units per period and auto-correlation coefficient of 0.2. It is assumed that the suppliers are located close to the OEM’s assembly facility. The simulation is run for 500 periods and all the results are reported based on 30 replication runs. Table 1 presents the other SC settings and parameters for simulation models for all four scenarios. For Scenarios 3 and 4, the 500 period simulation run consists of an introduction phase (periods 1–90), a peak supply consumption phase (periods 91–300), and a decline phase (periods 301–500). We assume that this phase distribution is a representative of lifecycle demand of a typical vehicle model. 5.2. Scenario 1: three echelon uncapacitated SC network with smooth supply consumption (base case)

5. Simulation experiments and results Simulation models extend the analytical formulas presented in Section 4 for the bullwhip effect and NS amplification for OEM to a three echelon SC (OEM, T1, and T2). The three-echelon automotive SC includes the following four scenarios: Scenario 1: without capacity constraints and smooth consumption process. Scenario 2: with capacity constraints and smooth consumption process. Scenario 3: step-change in supply consumption rate due to lifecycle demand and without capacity constraints; and

In this scenario, we studied a base case with three-stage SC with the OUT inventory policy under exponential smoothing forecasting. We did not consider any restriction on capacity. There were no step-changes in the supply consumption rate. The objective of running the base case analysis was to compare the results obtained using the analytical formulas with those from simulation. We have found that the bullwhip effect for OEM using (26) and that of the simulation model to be 2.26 and 2.25, respectively. Similarly, the values of net-stock amplification ratio calculated using our analytical formula (27) and using simulation

Table 1 SC settings and parameters used in simulation. SC parameters

Smooth supply consumption Without capacity constraints

Step-changes in supply consumption With capacity constraints

With capacity constraints

Peak demand: 5400 units per period Introduction: 60% of peak demand Decline: 80% of peak demand

Mean supply consumption rate

5400 units per period

Standard deviation of error term

OEM: 15% of mean supply consumption rate; T1: 20% of mean; T2: 25% of mean

Exponential smoothing forecasting parameters

Without capacity constraints

a ¼ 0.25; b ¼ 0.2 for all three echelons

Total (production þtransportation) lead-time

OEM ¼1 day; T1¼ 2 day; and T2 ¼2 day

Production capacity

Not applicable

OEM: 15% more than the expected demand per period T1: 10% more than OEM T2: 10% more than T1

Peak demand: a ¼ 0.3; b ¼0.25 for all 3 echelons Introduction: a ¼ 0.25; b ¼0.2 for all 3 echelons Decline: a ¼ 0.25; b ¼0.2 for all 3 echelons

Not applicable

OEM: 15% more than the expected demand per period T1: 10% more than OEM T2: 10% more than T1

324

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

models were found to be 3.02 and 2.99, respectively. As expected, both the bullwhip effect and NS amplification increase exponentially as we move up the chain (see Table 2).

5.3. Scenario 2: three echelon capacitated SC network with smooth supply consumption In the second scenario, we studied the same three stage SC but with a restriction on capacity at each echelon. As in Scenario 1, every order was fulfilled and shipped out to customer without any supply shortage while some orders were backlogged because of the limited capacity. As shown in Table 2, capacity constraint causes a drastic increase in net-stock amplification in raw materials at each echelon compared to that in the first scenario with infinite production capacity (e.g., for OEM, see 385.41 vs. 2.95). The significant rise in inventory variance in the upstream players can be attributed to: (1) in the finite production case, the production is limited by the actual order while raw materials are ordered based on the demand forecast; (2) the lead-times and service levels are higher in the given SC setting as we move up the chain. Furthermore, it has also been found that the finished goods inventory is higher for all three players in the infinite capacity

Table 2 Results for scenario analysis 1 and 2. Capacitated and uncapacitated SC with smooth supply consumption

Mean

Bullwhip effect OEM Capacitated Uncapacitated Tier 1 Capacitated Uncapacitated Tier 2 Capacitated Uncapacitated NS amplification OEM Capacitated Uncapacitated Tier 1 Capacitated Uncapacitated Tier 2 Capacitated Uncapacitated

Standard deviation

2.25 2.24

0.04 0.05

7.62 7.58

0.27 0.32

29.49 29.34

1.54 1.83

395.41 2.95

99.66 0.15

394.51 12.30

102.08 0.95

2248.01 50.99

458.28 15.88

case whereas there is almost no finished goods inventory at the suppliers’ nodes in the finite production capacity case (Table 3). A close look at uncapacitated and capacitated cases in Table 2 reveals that there is no evidence of any significant impact of the capacity constraint on the bullwhip effect in the given SC setting. In order to further explore the relationship between capacity restriction and the bullwhip effect, we formulated and tested several hypotheses as presented in Table 4. It should be noted here that the demand-capacity ratios differ in the different hypotheses tests. For example, the hypothesis tests H1 and H5 represent the most common setting where capacity at OEM is assumed to be 15% higher than the mean demand (e.g., to buffer against the demand variability) and the suppliers will have another 10% additional or ‘‘default’’ capacity as compared to their immediate customer (this is very typical in the automotive industry and often required by contracts). In comparison, hypothesis tests H2 and H6 (H4 and H8) represent a uniform 30% decrease (increase) in capacity across the board. Lastly, H3 and H7 represent the cases of capacity matching the demand. These percentages were chosen arbitrarily to demonstrate a general framework for testing the impact of capacity on the bullwhip effect and the net-stock amplification. In all cases, we compare the capacitated performance measures with their uncapacitated counterparts. The p-values for the hypotheses tests in Table 4 show that the NS amplifications are higher and statistically significant for all three echelons when capacity is limited. However, there was no sufficient evidence (see p-values in Table 4) to establish the impact of capacity constraint on the bullwhip effect at all three tiers of the simulated SC setting, especially when there was no short shipment of the orders. The above mentioned hypothesis tests H1–H8 represent the range of most practical cases. However, for the sake of completeness, we also ran two additional scenarios representing the extreme settings of capacity vs. demand. Those scenarios included: (1) when capacity was five times lower than the mean demand faced by OEM (represented by H9 and H10), and (2) when capacity was five times higher than the mean demand faced by OEM (represented by H11 and H12). The hypothesis tests on these extreme scenarios are presented in the last two rows of Table 4. The test results have shown that when the capacity was extremely low compared to the mean demand, it led to a significant increase in inventory variance amplification. However, the capacity had little or no impact on the inventory variance when the capacity was extremely (e.g., five times) higher than the mean demand except for OEM. It is worth noting here that this situation corresponds to the infinite capacity scenario. Yet again,

Table 3 Change in finished goods and raw materials inventory cost due to capacity constraint. Infinite capacity

Finished goods inventory Mean FGI cost (T1)

sT1 FG$ Mean FGI cost (T2)

sT2 FG$ Raw materials inventory Mean RMI cost (OEM)

sOEM RM$ Mean RMI cost (T1)

sT1 RM$ Mean RMI cost (T2)

sT2 RM$

Finite capacity

Introduction

Ramp up

Saturation

Introduction

Peak

Decline

$683 $51

$1568 $242

$1572 $232

$3 $0

$0 $0

$0 $0

$649 $61

$1667 $557

$1490 $272

$2 $0

$1 $1

$0 $0

$4830 $3855

$52,853 $12,411

$10,638 $3921

$487 $117

$2052 $293

$4269 $536

$5593 $6642

$77,663 $24,374

$17,870 $7966

$699 $80

$2096 $220

$3727 $435

$5006 $5639

$53,150 $19,714

$6267 $3661

$641 $80

$2704 $224

$4823 $413

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

325

Table 4 Hypothesis testing for impact of capacity constraints on bullwhip and NS amplification effects. Demand vs. capacity

Bullwhip effect

NS amplification

Hypothesis

Capacity ¼ 1.15  demand (baseline) Capacity o 0.7  demand (capacity o demand) Capacity ¼ demand Capacity 41.3  demand (capacity 4demand) Capacity ¼ 0.2  demand Capacity ¼ 5  demand

60.00

p-Values

H1 H2 H3 H4 H9 H10

Hypothesis

OEM

T1

T2

0.337 0.396 0.278 1.000 0.36 1.00

0.120 0.459 0.368 0.562 0.08 0.55

0.120 0.353 0.276 0.411 0.17 0.59

H5 H6 H7 H8 H11 H12

p-Values OEM

T1

T2

0.004 0.000 0.000 0.000 0.00 0.00

0.000 0.000 0.000 0.000 0.00 0.51

0.000 0.000 0.000 0.000 0.00 0.70

Bullwhip Effect

NS Amplification

Bullwhip Analysis

Inventory Variation Analysis

Bullwhip effect

Uncapacitated

Tier-1 Tier-2

40.00 30.00 20.00 10.00

Net Stock Amplification

OEM

50.00

0.00 Introduction

Peak

90.00 80.00 70.00 60.00 50.00 40.00 30.00 20.00 10.00 0.00

Decline

Bullwhip Analysis

Net Stock Amplification

Bullwhip effect

Peak

Decline

Inventory Variation Analysis (Note: Log Scale for NS Amplification)

50.00

Capacitated

Tier-2

Product Lifecycle Phase

60.00

40.00 30.00 20.00 10.00 0.00 Peak

Tier-1

Introduction

Product Lifecycle Phase

Introduction

OEM

Decline

Product Lifecycle Phase

1000.00

100.00

10.00

1.00 Introduction

Peak

Decline

Product Lifecycle Phase

Fig. 2. Bullwhip and NS amplification effects when considering life-cycle demand aspects.

regardless of capacity/demand ratio, for given SC settings, there was not enough evidence to make any strong conclusion about the relationship between the capacity restriction and the bullwhip effect. 5.4. Scenarios 3 and 4: three echelon capacitated and uncapacitated SC network with step changes in supply consumption In these scenarios, we studied a three-echelon SC with the OUT inventory policy under exponential smoothing forecasting technique and step-changes in OEM production volume (supply consumption rate) due to life-cycle demand. The impact of stepchanges in demand during life cycle phases for uncapacitated and capacitated scenarios on the bullwhip effect and NS amplification are shown in Fig. 2. For both the scenarios, that is, with and without capacity constraints, we observed that the performance of the SC deteriorated

in all measures during the peak demand. The bullwhip effect and net-stock amplification were significantly higher during the peak demand periods than those during slow demand or decline phase. Interestingly, the magnitudes of order variances and inventory variances were directly proportional to the size of the mean production rate for the uncapacitated scenario. As expected, the SC performance worsened on all accounts as we moved up the supply chain. Tier 2 supplier was the most vulnerable of all, in terms of both the bullwhip effect and the inventory variance, under step-changes in demand. As in the case of smooth supply consumption, Fig. 2 indicated that the capacity constraint did not have a notable effect on the bullwhip effect. We also noted that the inventory variance in the uncapacitated scenario was the highest during the ramp-up phase. However, the effect of capacity constraints was higher on net-stock amplification in the saturation phase than the introduction phase for all SC players (OEM, T1 and T2). This observation lead to the conclusion that OEM and suppliers, while adjusting their capacity

326

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

to level the production, should also beware of the impact of step changes on the inventory variance. In order to further explore the influence of step-changes in OEM production volumes on the performance of the SC, several additional hypotheses were formulated and tested using additional simulation results from an uncapacitated SC network scenario (see Table 5). The objective of these tests was to determine whether the mean bullwhip effects at different phases of the life-cycle demand are significantly different (statistically). The test results showed that for all echelons (OEM, T1, and T2), the mean bullwhip effects during peak phase were significantly higher than those during the introduction and decline phases (Table 5). The same trend was also seen with respect to the net-stock or inventory variance amplification.

Table 5 Hypothesis testing for the impact of life-cycle demand on bullwhip and NS amplification effects. Impact of life- Bullwhip effect cycle stages Hypothesis p-Values OEM Peak vs. introduction Peak vs. decline Decline vs. introduction

T1

NS amplification Hypothesis p-Values T2

OEM

T1

T2

H13

0.000 0.00 0.00 H16

0.000 0.00 0.00

H14

0.000 0.00 0.00 H17

0.006 0.00 0.00

H15

0.997 1.00 1.00 H18

0.004 0.17 0.98

10000 demand orders

Bullwhip Effect

8000 6000

f

4000 2000

Intro

0 1

26

51

Peak 76

101

126

151

176

201

Decline 226

251

276

301

326

351

376

401

426

251

276

301

326

351

376

401

426

451

476

501

14000 demand orders

12000

Bullwhip Effect

10000 8000 f

6000 4000 2000 0 1

26

51

76

101

126

151

176

Bullwhip Effect

Intro 26000 24000 22000 20000 18000 16000 14000 12000 10000 8000 6000 4000 2000 0

201

226

Peak

451

476

501

Decline

demand orders

f

1

26

51

Intro

76

101

126

151

176

201

226

251

276

Peak

301

326

351

376

401

Decline Lifecycle Demand

Fig. 3. Supply consumption and order fluctuations during different phases of the production life-cycle.

426

451

476

501

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

most unstable situation with the highest variability and the largest order quantity as compared to the demand faced by the OEM. We have also experimented with different levels of step change for ‘‘Introduction-to-Peak’’ and ‘‘Peak-to-Decline’’ transitions and observed that the bullwhip effect is more significant in the former transition than the latter transition Similar transient behavior can also be seen in the net-stock plots for the uncapacitated SC (Fig. 5). When there is sudden rise in demand (transition from the intro to the peak phase), the entire system experiences stock-outs as the NS levels become negative. As a response to this, all three players increase the order quantities, which eventually bring the NS into positive territory. Nevertheless, while the OEM is still facing stock-outs, the lower tiers are already building up the inventories. In the next transition (from the peak to the decline), when the demand starts declining, there is sudden drop in the net-stock as companies place smaller orders. However, there is still large accumulation of inventories available at T2 as compared to that at the OEM and T1. In summary, step-changes in the life-cycle demand process increases the nervousness of the entire SC with worst impact at the upstream echelons. As the system stabilizes over time, there are lower values of the bullwhip effect at the decline phase as a result.

Interestingly, the bullwhip effect during the decline phase was lower than that during the introduction phase for all three echelons. While the demand observed was higher during the decline phase than that during the introduction phase, this trend in the bullwhip effect might be a little surprising to SC managers. We argue (without proof) that this is attributable to the following: (1) the introduction phase is shorter (first 90 periods) than the decline phase (last 190 periods in a 500 period production simulation run) meaning that the system may be still dealing with ‘‘burn-in’’ effect or has not reached the steady state yet; (2) by the time we reach period 300 (the cut-off point for the decline phase), the system would have accumulated a large inventory, thereby, reducing the order size and variability in the order quantities. These phenomena can be observed in Fig. 3. Furthermore, it shows that order fluctuation (or the bullwhip effect) is significantly higher during peak demand (periods 90– 310) for all three echelons. 5.5. Analysis of transient behaviors of SC system during stepchanges in OEM production volume Figs. 4 and 5 provide a closer look on the propagation of the bullwhip effect and the inventory variance across the SC, especially the transient behavior of the system right after a step change in supply consumption rate at the OEM. As seen in Fig. 4, there is a sudden jump in the OEM orders (around period 90) when the lifecycle is transitioning from the introduction phase to the peak phase. It causes a huge forecast error at the upstream nodes which in turn results into monumental variations in the order quantities at both T1 and T2. On the other hand, when the life-cycle is transitioning from the peak phase to the decline phase, there is sudden reduction in the order quantity and system gets to a steady state with reduced variability (in line with demand pattern) after certain time but there is still a bullwhip effect at all levels.. As expected, the Tier 2 faces the

5.6. Sensitivity analysis In this section we present sensitivity analysis of interaction effects of capacity constraint with other supply chain and product life cycle parameters on the bullwhip effect and NS amplifications. The main goal of this exercise is to test the robustness of the results obtained through earlier scenario analyses. In particular, we are interested in the following three types of experiments: (i) Impact on the bullwhip effect and the NS amplification due to interaction between capacity constraints and other supply

50000

9000

Demand

Demand

40000

8000 6000 4000

Demand, Order

Order

Demand, Order

Demand, Order

10000

6000

3000

2000

Order

0

80

85

90

95

100

105

110

80

85

30000 20000

0

90

95

Period

100 105 110

80

Period

85

90

95

105

100

Period

80000

9000

8000

Order

10000

Demand

0

6000 4000 2000

Demand

295

300

6000

3000 Demand Order

Order

0 290

Demand, Order

Demand

Demand, Order

Demand, Order

327

305 Period

310

315

0 290 295 300 305 310 315 Period

Order

60000 40000 20000

0 290 295 300 305 310 315 Period

Fig. 4. Transient behavior of order quantity with respect to demand faced by each echelon during the transition (without capacity constraints).

110

328

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

7000 6000 5000 4000 3000 2000 1000 0 -1000 80 -2000 -3000

15000

Demand

10000

NS

Demand

50000

Demand

NS

40000

NS

30000

5000

20000

0 -5000 85

90

95

100

10000 80

85

90

95

100

105

-10000

5000

6000

4000

4000

3000

2000

2000 Demand

1000

NS

0 310

0 -2000

95

100

105

110

Period

Demand NS

8000

305

90

-30000

10000

6000

300

85

Period

7000

295

0 -10000 80 -20000

-15000 Period

290

110

290

295

300

305

310

315

18000 16000 14000 12000 10000 8000 6000 4000 2000 0 -2000 290

Demand NS

295

300

305

310

315

315

Period

Period

Period

Fig. 5. Trajectories of variations in net-stocks with respect to demand faced by each echelon during the transition (without capacity constraints).

chain parameters like forecasting constants, autocorrelation coefficient of supply consumption process, and duration of lifecycle phases. (ii) Impact of change in duration of lifecycle phases on the bullwhip effect and NS amplification. (iii) Impact of lifecycle phases with smooth supply consumption process on the bullwhip effect and NS amplification. In the first sensitivity analysis, we ran three different sets of experiments that were: (a) by varying the exponential smoothing constants for the mean and the variance of forecasted demand; (b) by changing the autocorrelation coefficient for demand (i.e., supply consumption process in this case); and (c) by changing the duration of the peak phase. In order to compare the results, these experiments were repeated for the uncapacitated case as well. The results of the sensitivity analyses along with the hypothesis tests are presented in Table 6a and b. It has been found that there was no statistically significant difference in the magnitude of the bullwhip effects between the uncapacitated and capacitated supply chains for all 3 nodes regardless of changes in the above-mentioned supply chain parameters (see Table 6a). Therefore, based on these simulation results, we can generalize our earlier conclusion in that the bullwhip effect does not depend on the capacity constraints for the given supply chain setting. On the other hand, the mean NS amplification values were significantly higher at all nodes in the capacitated SC than those in the uncapacitated SC (Table 6b). In the second sensitivity analysis, we ran two experiments one each for a capacitated three-stage supply chain with unequal and equal duration of lifecycle phases. Both the cases consist of stepchange in supply consumption process. The results of hypothesis tests comparing the two experiments are shown in Table 7. As anticipated, there is significant difference in the magnitudes of the bullwhip effect between the introduction phase and the peak phase, and between the peak and the decline phase regardless of duration of lifecycle phases. Similar results were observed for the NS amplification. Table 7 also depicts the results of the third sensitivity analysis in which we ran two experiments by varying the nature of the demand process, that is, lifecycle phases with

and without step-changes in supply consumption process. In this case, it has been found that for smooth supply consumption process the bullwhip effect during the introduction phase was significantly lower than that during the peak. However, there was no statistically significantly difference in the bullwhip effect and net stock amplification between the peak and the decline phases. We found these results to be very intuitive because of the very nature of the whole supply consumption process as the system has yet to be at steady state. On the other hand, both the bullwhip effect and the net stock amplifications were significantly different between different lifecycle phases for step-change supply consumption process. In other words, based on the results of second and third sensitivity analyses, it is safe to say that both the bullwhip effect and net stock amplifications are significantly increased due to step-change in the demand during different life cycle phases of a product.

6. Conclusions and further work The bullwhip effect has generated tremendous interest in the SC research community. Numerous simulation-based and analytical models have been developed in the literature to quantify as well as to reduce the bullwhip effect. Apart from the information sharing issue, the prior works mostly dealt with traditional problem settings such as the correlated demand process, the common forecasting techniques, and an OUT inventory policy. Most of the prior models have been developed around a single stage supply chain. Furthermore, these models have rarely considered some of the real-world issues such as capacity restriction and life-cycle demand. This paper has presented an analysis of a three echelon SC considering capacity restriction and stepchanges in supply consumption rate due to life-cycle demand phases. An analytical expression has been presented to quantify the net-stock amplifications for a single stage SC using exponential forecasting. The analysis was then extended to a three stage SC network via simulation and empirical testing. The threeechelon SC setting consisted of an OEM and Tier 1, Tier 2 suppliers.

Table 6 Sensitivity analysis results for interaction of capacity constraint with other supply chain characteristics and parameters. (a) Bullwhip effect

Bull Whip Effect

Capacitated OEM Mean Std. deviation Tier 1 Mean Std. deviation Tier 2 Mean Std. deviation

Lifecycle phases duration sensitivity

Forecasting parameters sensitivity

Auto-correlation ¼ 0.4

Auto-correlation ¼ 0.0

Shortened peak

Equal values in all phases

Higher alpha and beta at Intro

Lower alpha and beta at Intro

Intro

Intro

Intro

Intro

Intro

Intro

Peak

Decline

Peak

Decline

Peak

Extended peak Decline

Intro

Peak

Decline

Peak

Decline

Peak

Decline

2.43 0.10

2.13 0.08

2.37 0.09

2.72 0.11

2.34 0.06

2.27 0.10

2.62 0.13

2.25 0.07

2.33 0.09

2.58 0.09

2.26 0.08

2.29 0.10

2.30 0.08

2.25 0.07

2.60 0.11

2.28 0.05

2.28 0.07

2.32 0.09

2.60 0.11

2.58 0.10

7.49 0.66

9.65 0.79

6.88 0.49

8.69 0.63

11.85 0.92

8.23 0.42

8.03 0.65

11.19 1.08

7.61 0.42

8.34 0.61

10.76 0.70

7.66 0.52

8.12 0.64

8.09 0.62

7.57 0.49

10.97 0.97

7.94 0.40

7.81 0.49

8.28 0.63

10.88 0.92

10.64 0.86

29.53 3.64

45.15 5.46

25.58 2.60

36.02 3.92

59.00 6.76

32.72 2.48

32.83 3.95

55.07 7.68

29.34 2.32

33.96 3.61

51.78 5.18

29.35 2.93

32.84 3.80

32.21 4.17

28.98 2.83

53.67 7.50

31.41 2.68

30.48 2.64

33.79 3.91

52.71 6.66

50.77 6.23

2.20 0.12

2.42 0.09

2.16 0.07

2.39 0.07

2.75 0.10

2.34 0.06

2.31 0.08

2.59 0.10

2.24 0.04

2.31 0.09

2.59 0.07

2.27 0.08

2.30 0.10

2.29 0.07

2.25 0.07

2.62 0.17

2.28 0.07

2.27 0.06

2.34 0.09

2.62 0.12

2.56 0.12

7.48 0.80

9.55 0.71

7.00 0.42

8.79 0.54

12.13 0.84

8.21 0.43

8.21 0.51

10.94 0.77

7.53 0.28

8.25 0.64

10.84 0.55

7.75 0.55

8.17 0.66

8.00 0.47

7.65 0.42

11.15 1.26

7.96 0.50

7.74 0.47

8.46 0.68

11.07 0.99

10.49 0.95

29.46 4.37

44.11 4.98

26.08 2.22

36.75 3.71

61.31 6.37

32.60 2.69

33.27 3.08

53.13 5.60

28.85 1.79

33.84 3.83

52.49 3.81

30.04 3.32

33.43 3.87

31.75 2.74

29.52 2.43

54.65 7.88

31.69 2.96

29.86 3.01

34.85 4.35

53.98 6.87

50.00 6.59

0.13 0.31 0.43

0.40 0.51 0.46

0.24 0.22 0.18

0.95 0.89 0.85

0.09 0.24 0.63

0.32 0.31 0.27

0.50 0.39 0.36

0.39 0.58 0.90

0.63 0.96 0.55

0.63 0.52 0.40

0.70 0.77 0.55

0.61 0.53 0.62

1.00 0.50 0.43

0.59 0.54 0.62

1.00 0.87 0.70

0.55 0.57 0.40

0.39 0.29 0.33

0.50 0.44 0.47

0.49 0.52 0.64

Net Stock Amplification Auto-regressive demand process sensitivity

Lifecycle phases duration sensitivity

Forecasting parameters sensitivity

Auto-correlation ¼0.4 Auto-correlation ¼ 0.0 Shortened peak Extended peak Intro Peak Decline Intro Peak Decline Intro Peak Decline Intro Peak

Equal values in all phases Higher alpha and beta at Intro Lower alpha and beta at Intro Decline Intro Peak Decline Intro Peak Decline Intro Peak Decline

3.18 0.30

3.42 0.27

3.31 0.33

2.66 0.26

3.21 0.21

2.85 0.20

2.89 0.30

3.59 0.29

3.13 0.23

3.09 0.29

3.36 0.22

3.32 0.37

3.03 0.27

3.34 0.23

3.18 0.24

3.12 0.30

3.37 0.23

3.17 0.28

3.04 0.37

3.47 0.24

3.33 0.28

14.10 1.99

17.55 1.67

14.32 1.95

11.19 1.34

15.90 1.66

11.29 1.01

12.12 1.74

18.23 2.15

12.85 1.34

13.30 2.29

16.69 1.38

14.13 1.93

13.28 1.68

14.52 1.60

13.23 1.26

14.99 1.84

14.51 1.30

13.05 1.41

13.29 2.24

17.25 1.71

16.17 2.06

53.44 8.38

74.79 10.03

49.94 8.30

48.39 7.31

77.77 13.16

43.36 5.74

49.53 7.79

84.27 14.96

46.54 6.14

53.21 9.92

75.35 8.89

51.50 8.43

54.06 8.86

55.08 8.30

47.13 5.49

69.60 10.11

54.25 6.82

47.49 5.65

54.46 9.39

77.72 10.77

71.96 11.80

18.10 102.31 7.79 39.73

83.95 26.26

18.05 105.84 7.85 40.08

84.57 25.53

19.58 12.54

38.12 120.45 13.13 39.77

20.64 12.24

128.24 45.85

35.12 15.50

18.89 9.87

83.23 29.40

67.28 21.75

19.90 9.03

81.28 22.53

70.68 26.74

18.77 7.58

72.95 25.58

65.27 26.58 329

Capacitated OEM Mean Std. deviation Tier 1

Decline

2.20 0.10

Hypothesis test (p-value) Uncapacitated aCapacitated OEM 1.00 0.69 Tier-1 0.96 0.61 Tier-2 0.95 0.44 (b) NS amplification

Uncapacitated OEM Mean Std. deviation Tier 1 Mean Std. deviation Tier 2 Mean Std. deviation

Peak

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

Uncapacitated OEM Mean Std. deviation Tier 1 Mean Std. deviation Tier 2 Mean Std. deviation

Auto-regressive demand process sensitivity

0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 Hypothesis test (p-value) Uncapacitated ocapacitated OEM 0.00 Tier-1 0.00 Tier-2 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

0.00 0.00 0.00

788.31 239.08 158.52 55.74 468.74 118.17 490.25 117.99 118.89 835.10 412.68 31.61 147.71 88.03

121.02 854.34 407.14 31.51 162.87 87.04

144.17 477.80 630.60 58.60 129.15 179.84

137.32 1387.75 228.20 45.57 263.51 86.51

138.98 43.87

399.23 111.53

257.04 99.55

436.21 157.81

867.02 224.05

100.60 36.96 116.57 33.47 31.72 11.99 79.85 31.69 91.49 23.64 39.78 15.69 73.34 23.59 96.16 29.79 28.90 9.49 44.68 18.34 187.62 55.11 30.26 10.87 66.78 118.32 20.37 40.45 30.37 14.02 87.49 22.87 27.77 140.93 7.26 38.61 87.89 22.52 27.60 135.57 7.37 36.92

Mean Std. deviation Tier 2 Mean Std. deviation

Auto-correlation ¼0.4 Auto-correlation ¼ 0.0 Shortened peak Extended peak Intro Peak Decline Intro Peak Decline Intro Peak Decline Intro Peak

Lifecycle phases duration sensitivity Net Stock Amplification Auto-regressive demand process sensitivity

(b) NS amplification

Table 6 (continued )

Equal values in all phases Higher alpha and beta at Intro Lower alpha and beta at Intro Decline Intro Peak Decline Intro Peak Decline Intro Peak Decline

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

Forecasting parameters sensitivity

330

Table 7 Sensitivity analysis results for impact of change in duration of lifecycle phases on the bullwhip effect and NS amplification. p-Values for hypothesis tests (2sample T-test)

Bullwhip effect OEM Intro vs. peak Peak vs. decline Tier 1 Intro vs. peak Peak vs. decline Tier 2 Intro vs. peak Peak vs. decline

Nature of demand change

Duration of lifecycle phases

Smooth changes with lifecycle phases

Stepchange with lifecycle phases

Equal duration with stepchange in demand

Unequal duration with step change in demand

0.02

0.00

0.00

0.00

0.81

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.78

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.65

0.00

0.00

0.00

0.00

0.04

0.00

0.01

0.01

0.01

0.00

0.06

0.00

0.00

0.01

0.00

0.00

0.00

0.00

0.00

0.00

0.00

NS amplification OEM Intro vs. 0.00 peak Peak vs. 0.54 decline Tier 1 Intro vs. 0.00 peak Peak vs. 0.52 decline Tier 2 Intro vs. 0.00 peak Peak vs. 0.33 decline

While each node had capacity restriction, there were no short shipments meaning that back orders were allowed. We studied the impact of step-changes in the life-cycle demand (which caused sudden change in the supply consumption rates at each echelon) on the bullwhip effect and NS amplification, at three phases (introduction, peak, and decline or end-of-life) of the product life-cycle. Four scenarios were analyzed through simulation. The first scenario involved a simple three-echelon SC network with unlimited capacity, AR(1) demand process, and exponential smoothing forecasting. As anticipated, it showed that key performance metrics such as the bullwhip effect and NS amplification increased exponentially as we moved up the chain. The second scenario considered a capacitated SC network with all the others parameters remaining unchanged from the previous scenario. It was discovered that the NS amplification increased at much higher rate in the capacitated network than that in the uncapacitated network at all echelons. Obviously, the rate of the NS increase was the worst for T2, followed by T1 and OEM, respectively. Interestingly enough, although the bullwhip effect was present at all echelons, a direct correlation between the bullwhip effect and capacity level could not be established. In the third and fourth scenarios, we studied relationship between the stepchanges in production rate due to life-cycle demand and the change in the performance level of the SC system with and without capacity constraints, respectively. The simulation results

B. Nepal et al. / Int. J. Production Economics 136 (2012) 318–331

showed that performance of a system as a whole deteriorated when there was a step-change in the life-cycle demand. Particularly, key performance metrics such as the bullwhip effect and the NS amplification rose sharply during the peak phase. Lastly, the NS amplification was less adversely affected by the limited capacity when there were step changes in the supply consumption rate. In addition, a number of sensitivity analyses were performed to examine the impact on the supply chain performance due to interaction of capacity constraint with the other SC characteristics and operating parameters like forecasting, autocorrelation coefficient, and duration of life cycle phases. The sensitivity analysis results also showed that there was no evidence to establish the correlation between the capacity constraint and the bullwhip effect regardless of its interaction with the other parameters. However, there was significant difference in the magnitudes of the NS amplification between the capacitated and the uncapacitated supply chains in all the cases. Sensitivity analysis of duration of life cycle phases on the supply chain performance was also performed to ascertain the impact of a step-change on the bullwhip effect and the NS amplification. However, the results have confirmed our original hypothesis that a step-change in demand during different phases of life cycle would directly contribute to the bullwhip effect and the NS amplification regardless of duration of lifecycle phases. While this study has revealed many important insights with respect to capacity constraint and step-changes in demand, it has made some critical assumptions. For example, the paper assumed complete shipment of orders (although with penalty costs as presented in Lee et al., 2000) and common SC policies throughout the life-cycle regardless of the demand situation. Future work will involve the analysis with both incomplete shipments and target finished goods inventory. Furthermore, we used AR(1) to model demand process. However, it would be interesting to analyze if the similar conclusions on the bullwhip effect hold true in a situation where unexpected events such as product recalls can affect the demand.

Acknowledgments We would like to thank the anonymous reviewers for constructive comments and recommendations, which helped to improve the paper. References Agrawal, S., Sengupta, R.N., Shanker, K., 2009. Impact of information sharing and lead-time on bullwhip effect and on-hand inventory. European Journal of Operational Research 192 (2), 576–593. Alwan, L., Liu, J., Yao, D., 2003. Stochastic characterization of upstream demand process in a supply chain. IIE Transactions 35 (3), 207–221. Berry, D., Towill, D.R., 1982. Material flow in electronic product based supply chains. International Journal of Logistics Management 3 (2), 77–94. Boute, R.N., Lambrecht, M.R., 2007. Exploring the Bullwhip Effect by Means of Spreadsheet Simulation. Department of Decision Sciences and Information Management (KBI), Katholieke Universiteit Leuven. KBI 0706. URL: /http:// papers.ssrn.com/sol3/papers.cfm?abstract_id=1094649S. Cachon, G.P., Randall, T., Schmidt, G.M., 2007. In search of the bullwhip effect. Manufacturing and Service Operations Management 9 (4), 457–479. Cachon, G.P., Olivares, M., 2010. Drivers of finished-goods inventory in the U.S. automobile industry. Management Science 56 (1), 202–216. Chen, Y.F., Drezner, Z., Ryan, J.K., Simchi-Levi, D., 2000a. Quantifying the bullwhip effect in a simple supply chain: the impact of forecasting, lead-times, and information. Management Science 46 (3), 436–443. Chen, Y.F., Ryan, J.K., Simchi-Levi, D., 2000b. The impact of exponential smoothing forecasts on the bullwhip effect. Naval Research Logistics 47, 269–286.

331

Coppini, M., Rossignoli, C., Rossi, T., Strozzi, F., 2010. Bullwhip effect and inventory oscillations analysis using the beer game model. International Journal of Production Research 48, 3943–3956. Dejonckheere, J., Disney, S.M., Lambrecht, M.R., Towill, D.R., 2003. Measuring and avoiding the bullwhip effect: a control theoretic approach. European Journal of Operational Research 147 (3), 567–590. Disney, S.M., Towill, D.R., 2003. On the bullwhip and inventory variance produced by an ordering policy. OMEGA, The International Journal of Management Science 31, 157–167. Disney, S.M., Farasyn, I., Lambrecht, M., Towill, D., Van De Velde, W., 2004a. Dampening Variability by Using Smoothing Replenishment Rules. Available at SSRN: /http://ssrn.com/abstract=873588S. Disney, S.M., Towill, D.R., Van De Velde, W., 2004b. Variance amplification and the golden ratio in production and inventory control. International Journal of Production Economics 90, 295–309. Disney, S.M., Farasyn, I., Lambrecht, M., Towill, D.R., Van De Velde, W., 2006. Taming the bullwhip effect whilst watching customer service in a single supply chain echelon. European Journal of Operational Research 173 (1), 151–172. Duc, T.T.H., Luong, H.T., Kim, Y.-D., 2008. A measure of bullwhip effect in supply chains with a mixed autoregressive-moving average demand process. European Journal of Operational Research 187 (1), 243–256. Duc, T.T.H., Luong, H.T., Kim, Y.-D., 2010. Effect of third party warehouse on bullwhip effect and inventory cost in supply chains. International Journal of Production Economics 124 (2), 395–407. Forrester, J., 1961. Industrial Dynamics. MIT Press, Cambridge, MA. Hoberg, K., Thoneman, U., Bradley, J.R., 2007. Analyzing the effect of inventory policies on the nonstationary performance with transfer functions. IIE Transactions 39, 911–924. Hosada, T., Disney, S.M., 2006. On variance amplification in a three-echelon supply chain with minimum mean square error forecasting. OMEGA, The International Journal of Management Science 34, 344–358. Holweg, M., Disney, S.M., Hines, P., Naim, M., 2005. Towards responsive vehicle supply: a simulation-based investigation into automotive scheduling systems. Journal of Operations Management 23, 507–530. Kahn, J.A., 1987. Inventories and the volatility of production. American Economic Review 77 (4), 667–679. Kaipa, R., Korhonen, H., Hartiala, H., 2006. Planning nervousness in a demand supply network: an empirical study. The International Journal of Logistics Management 17 (1), 95–113. Kisiel, R., 2008. Honda Plant Adds Parts Center. Automotive News (November: 24). Lee, H.L., Padmanabhan, P., Whang, S., 1997. Information distortion in a supply chain: the bullwhip effect. Management Science 43 (4), 543–558. Lee, H.L., So, K.C., Tang, C.S., 2000. The Value of Information Sharing in Two-Level Supply Chain Management Science 46 (5), 626–643. Lee, H.L., Wang, S., 2000. Information sharing in supply chain. International Journal of Technology Management 20 (3), 373–387. Luong, H.T., 2007. Measure of bullwhip effect in supply chains with autoregressive demand process. European Journal of Operational Research 180 (3), 1086–1097. Luong, H.T., Phien, N.H., 2007. Measure of bullwhip effect in supply chains: the case of high order autoregressive demand process. European Journal of Operational Research 183, 197–209. Mahajan, V., Muller, E., 1979. Innovation diffusion and new product growth models in marketing. Journal of Marketing 43 (4), 55–68. Moyaux, T., Chaib-draa, B., D’Amours, S., 2007. Information sharing as a coordination mechanism for reducing the bullwhip effect in a supply chain. IEEE Transactions on Systems, Man, and Cybernetics—Part C: Applications and Reviews 37 (3), 396–408. Raghunathan, S., 2001. Information sharing in a supply chain: a note on its value when demand is nonstationary. Management Science 47 (4), 605–610. Reddy, A.M., Rajendran, C., 2005. A simulation study of dynamic order-up-to policies in a supply chain with non-stationary customer demand and information sharing. International Journal of Advanced Manufacturing Technology 25, 1029–1045. Sterman, J., 1989. Modeling managerial behavior: misperceptions of feedback in a dynamic decision making experiment. Management Science 35 (3), 321–339. Wang, J.-L., Kuo, J.-H., Chou, S.-Y., Wang, S.-Z., 2010. A comparison of bullwhip effect in a single stage supply chain for autocorrelated demands when using correct, MA, and EWMA methods. Expert Systems with Applications 37, 4726–4736. Wanphanich, P., Kara, S., Kayis, B., 2010. Analysis of the bullwhip effect in multiproduct, multi-stage supply chain systems—a simulation approach. International Journal of Production Research 48, 4501–4517. Williams, T.M., 1982. Reorder levels for lumpy demand. Journal of Operational Research Society 33, 185–189. Zhang, X., 2004. The impact of forecasting methods on the bullwhip effect. International Journal of Production Economics 88, 15–27. Zipkin, P.H., 2000. Foundations of Inventory Management. McGraw-Hill, New York.