SUPPLY CHAIN APERIODICITY, BULLWHIP AND

SUPPLY CHAIN APERIODICITY, BULLWHIP AND STABILITY ANALYSIS WITH JURY’S INNERS STEPHEN M. DISNEY

Abstract. Jury (1974) presents a novel method for the analysis of dynamical systems based on matrices of the coefficients of the systems transfer function and its ”Inners”. Here we exploit his procedure for an analysis of a supply chain replenishment or ordering decision known as the Order-Up-To policy. We study the discrete time case and generalise the classical Order-Up-To policy by the addition of two independent proportional controllers in the policy’s feedback loops. The addition of the proportional controllers is well-known to allow the Order-Up-To policy to eliminate the bullwhip problem and we quantify this herein using Jury’s Inners approach. However, care has to be taken with the use of independent controllers as they can introduce stability problems. This is because the roots of the characteristic equation become complex, and they may even move out of the unit circle in the z-plane. We identify the conditions of stability using Jury’s Inners approach. We also investigate further the root distribution in the characteristic equation to identify the conditions under which the Order-Up-To policy is aperiodic. An aperiodic system has only a limited number of maxima and minima in its dynamic response. Thus aperiodicity is an important characteristic of a supply chain replenishment policy as it will not induce rogue seasonality.

1. Introduction ˚ Astr¨ om, Jury and Agniel (1970) present a matrix-based procedure for evaluating complex integrals of the form I 1 B(z)B(z −1 ) dz (1.1) I= 2πı A(z)A(z −1 ) z along the unit circle in a positive direction in the complex plane when A and B are polynomials with real coefficients of z. This approach was refined further by Jury (1974) and ˚ Astr¨ om (1970) and their two methodologies are surprisingly powerful and easy to use. Complex integrals of the form of 1.1 often appear in the analysis of the ”bullwhip effect” in supply chains when control theory is exploited, Dejonckheere, Disney, Lambrecht and Towill (2003). The bullwhip effect is a common term used to describe how the long-run variance of the orders amplifies as it passes through each decision making point in a supply chain, Lee, Padmanabhan and Whang, (2004, 1997a and b), Metters (1997), Disney and Towill (2003). Thus, several echelons up a supply chain, production and distribution orders bear very little resemblance Date: July 5th, 2007. Key words and phrases. Supply chains, aperiodicity, bullwhip, stability, z-transforms, orderup-to policy. 1

Disney, S.M., (2008) “Supply chain aperiodicity, bullwhip and stability analysis with Jury’s Inners”, IMA Journal of Management Mathematics, Vol. 19, No. 2, pp101-116. DOI: 10.1093/imaman/dpm033.

2

STEPHEN M. DISNEY

to the end consumers demand. The bullwhip effect leads to unnecessary costs in supply chains and thus the phenomena has received a lot of academic attention recently. This contribution is concerned with utilizing Jury’s (1974) important work in an investigation of a common supply chain replenishment rule for bullwhip effects. We will show that the evaluation of the complex integral, I, is closely related to the issue of stability, another important supply chain characteristic that has been studied recently, see Warburton, Disney, Towill and Hodgson (2003), Riddalls and Bennett (2002) and Disney and Towill (2002). Stability is a fundamental characteristic of a supply chain ordering policy. An unstable policy will react to any demand pattern with a response that oscillates with ever increasing amplitude and is thus inherently undesirable. A stable supply chain ordering policy on the other hand will react to any finite demand signal and after a finite period of time return to steady state conditions. Obviously this is a necessary property of a practical supply chain replenishment decision. We further consider the issue of aperiodicity. This is the condition that, after a finite demand, a system will return to rest without an undue number of oscillations. This is an important criterion as it will contribute to the reduction of the tendency of a supply chain to resonate and introduce rogue seasonality or a fake business cycle, Forrester (1958). As these demand oscillations are very costly (we assume supply chain oscillations will cause capacity losses and labor scheduling and training problems), an aperiodic response by a supply chain replenishment decision is highly desirable. Aperiodicity is effectively a root distribution problem that we investigate via the Hurwitz matrix. 2. Preliminaries We are following the work of Jury (1974) and refer readers to Jury (1974) for any required proof to his approach. Let us be concerned with the exploitation of his important contribution and begin with some notation. Let z be the z-transform P∞ −t operator, F (z) = f (t)z . F (z) usually manifests itself as a ratio of two t=0 B(z) polynomials in z, F (z) = A(z) , where; (2.1)

B(z) = bn z n + bn−1 z n−1 + ... + b1 z + b0

(2.2)

A(z) = an z n + an−1 z n−1 + ... + a1 z + a0 , an > 0.

Thus F (z) is the discrete time, z-domain, transfer function of the dynamic system described by f (t), a function of the discrete (time) variable t. We will use lower case letters to denote variables in the time domain and corresponding upper case letters in the frequency domain. In practice however, it is usually not necessary to start from the time domain, as a block diagram most easily defines a system. This is easily constructed from a systems description using common control engineering techniques (for example see Nise, 1995) and allows the derivation of the systems transfer function directly in the z-domain. The block diagram of a well-known supply chain ordering or replenishment algorithm is shown below in Figure 1. We refer readers to Dejonckheere, Disney, Lambrecht and Towill (2003 and 2004) for its detailed derivation where they show that the policy is a generalisation of the OUT policy. Although, here follows a brief


SUPPLY CHAIN APERIODICITY, BULLWHIP AND STABILITY ANALYSIS

3

b introduction. First we can see the OUT policy has three feed-forward loops; D, b is a forecast of the demand process and is added directly to the DWIP and TNS. D replenishment order, O. Here the forecast is generated by exponential smoothing, 1 , where Ta is the average age of the data with a smoothing constant of α = 1+T a in the exponential smoothing forecast. Another feed-forward loop is used to set b The Target Net Stock gain (a) the Target Net Stock (TNS) position, TNS= aD. is used to achieve a given stock availability such as the fill-rate. Disney, Farasyn, Lambrecht, Towill and van de Velde (2006) investigate the link between bullwhip and the fill-rate for the ARMA(1,1) demand process. Interestingly, bullwhip reduction may not necessarily increase inventory holding costs or reduce customer service. Finally the forecast of demand is also used to generate a dynamic Work In b Progress target, DWIP by multiplying the forecast by T p , that is DWIP= T p D. T p is our estimate of the production, distribution or replenishment lead-time, Tp , which we require to determine a target for the ”orders placed but not yet received” or the WIP levels. The replenishment delay Tp , is the time between placing an order and receiving it in stock. The replenishment delay in Figure 1 is represented by the delay operator z −(Tp +1) , where the additional unit of delay is required to ensure the correct sequence of events. That is, we receive stock and satisfy demand at the beginning of the replenishment period and observe the inventory position and placed orders at the end of the period. Thus even if our physical lead-time, (Tp ) is zero, we would not observe the received stock from the suppliers or the production system until the end of the next period, hence the additional unit delay. Disney and Towill (2005) have investigated the consequences of getting the estimated delay (T p ) wrong and highlight a solution to the problem that involves calculating the WIP in a specific manner. We adopt their approach herein. Figure 1 also shows that the OUT policy has two feedback loops; the difference between the Target Net Stock (TNS) and the actual Net Stock (NS), and the difference the between Desired WIP (DWIP) and the actual WIP. Usually in inventory models these differences, or errors, are accounted for in full in the replenishment orders, thus traditionally Ti = Tw = 1. Herein we will consider the case of incorporating two proportional controllers into the feedback loops; a net stock proportional controller, T1i , and a WIP proportional controller, T1w . This generalisation of the OUT policy allows superior economic performance (Disney and Grubbstr¨om (2004) and Chen and Disney (2003)) and greater flexibility in shaping its dynamic response when balancing bullwhip, inventory and customer service targets, Disney, Farasyn, Lambrecht, Towill and van de Velde (2006). However the addition of the proportional controllers also introduces dynamic problems such as stability. Thus the benefits previously reported may be lost if care is not taken when designing replenishment algorithms. Furthermore, the generalisation is not new, it has a long history, for example see; Magee (1956), Deziel and Eilon (1967), Towill (1982), Bertrand (1986), Sterman (1989) to name a few. In this paper we are concerned with a single echelon of a supply chain. However, recently a number of important multi-echelon supply chain papers have appeared. Hoberg, Bradley and Thonemann (2007) investigate stability issues and bullwhip effects with a z-transform based approach for a multi-echelon supply chain under three operating policies; a base-stock, echelon stock and a replenishment rule with only an inventory feedback controller (rather than an inventory and WIP feedback


4

STEPHEN M. DISNEY

Figure 1. Block diagram of the generalized Order-Up-To policy controller). Gaalman and Disney (2007) have recently derived an optimal multiechelon replenishment rule for minimising the total supply inventory costs. The policy is optimal in the sense that it is the best possible linear policy and was derived using state space methods. Hosoda and Disney (2006) study a three echelon supply chain with AR(1) demand at the end customer. The evolution of the demand process as it passes through the ”order-up-to” policy is derived and conditional expectation is used to create minimal mean squared error forecasts. 3. Transfer function and Jury’s matrices of coefficients We are interesting in the dynamic (time-varying) performance of the Demand rate and the Order rates. The Demand rate is a discrete time series of the customer demand or requirements. The Order rate is a discrete time series of the replenishment orders placed on the supplier (or the production requests placed on a factory). We may manipulate the block diagram (Figure 1), to obtain following transfer function that relates the Order rate to the Demand rate. This transfer function completely describes the linear dynamic behaviour of our system. (3.1)

O(z) =

z 1+Tp (T p Ti (z − 1) − (a + Ta + Ti )Tw + (1 + a + Ta + Ti )Tw z) (Ta (z − 1) + z)(Tw + Ti ((1 + Tw (z − 1))z Tp − 1))

Jury’s Inners approach exploits matrices that are a function of the order (n) of the transfer function. Thus we have to specify a value of Tp here. We choose Tp = 2 as the case of Tp = 1 can be readily solved using time domain techniques such as those exploited by Disney and Towill (2003) and Disney, Towill and van de Velde (2004). Herein we also assume the demand pattern is a stationary independently and identically distributed (i.i.d.) stochastic process. Although this is a rather restrictive assumption we do know this type of demand process presents at least one product in Procter and Gamble’s product range, (Disney, Farasyn, Lambrecht, Towill and Van de Velde (2006). The best possible forecast of an i.i.d. process (that is, the one with the minimum mean squared error) is the unconditional mean of the demand process. If demand is not i.i.d. then other types of forecasting mechanisms may be more desirable. For example for some integrated demand processes then exponential smoothing (with −0.5 < Ta < ∞) will be optimal, Chatfield et al, (2001). If demand is of



5

an ARIMA type, the Box and Jenkins (1970) approach will be helpful. In any case, the general goal is to create two forecasts; one is a forecast of the demand over the replenishment lead-time, the other is a forecast in the period after the replenishment lead-time. The variance of the inventory (net stock) levels will equal the sum of the variances of the forecast errors of demand over the lead-time and review period. This was proved, for a single echelon of a supply chain by Vassian (1955) and extended to the multi-echelon case by Hosoda and Disney (2006). Thus, it is always important to create the best possible forecasts from an inventory cost perspective. We may also set the mean demand to zero without loss of generality as we are considering a linear system. The linear system assumption also implies that all unmet demand is fully backlogged (or available from another source). To achieve the constant forecast of the unconditional mean we set Ta = ∞. The gain in the inventory feed-forward loop, a, and the Work-In-Progress (WIP) feedforward gain, T p , will simplify out of the system’s transfer function when Ta = ∞. This is the optimal exponential forecasting constant for the demand process we have assumed. The transfer function then becomes B(z) Tw z 3 = . A(z) Tw + Ti (z − 1)(1 + z + Tw z 2 ) This transfer function can be expressed in the form of 2.1 and 2.2 as coefficients of z with

(3.2)

O(z) =

(3.3)

b0 = b1 = b2 = 0, b3 = Tw

and (3.4)

a0 = Tw − Ti , a1 = 0, a2 = Ti (1 − Tw ), a3 = Ti Tw .

An important sub-set of 3.2 occurs when Ti = Tw . In this case the transfer function simplifies down to z (3.5) O(z) = . 1 + Ti (z − 1) Jury’s Inners method makes use of the determinants of matrices whose elements consist of the coefficients of the transfer function. We require four matrices here. They are;   an an−1 an−2 . . . a0  0 an an−1 . . . a1     0 0 an . . . a2  (3.6) Xn+1 =    .. .. .. ..  ..  . . .  . . 0 0 0 . . . an

(3.7)

Yn+1



0

0

...

0

   =   

0 .. .

0

...

a0

a0 .. .

0 a0 a1

a0 a1 a2

... ... ...

an−2 an−1 an

0 a0

       


6

STEPHEN M. DISNEY

 (3.8)

Xn−1

(3.9)

   =  

Yn−1



an 0 0 .. .

an−1 an 0 .. .

an−2 an−1 an .. .

... ... ... .. .

a2 a3 a4 .. .

0

0

0

0

an



0

0

...

0

   =   

0 .. .

0

...

a0

a0 .. .

0 a0 a1

a0 a1 a2

... ... ...

an−4 an−3 an−2

0 a0

     

       

Equations 3.6 and 3.7 are required for the bullwhip analyses, while 3.8 and 3.9 are required for the stability investigation. Here we can see the concept of the ”Inners”. Xn−1 is the first Inner of Xn+1 . 4. Stability via Jury’s Inners The question of stability is a fundamental aspect of dynamic systems. A stable system will react to a finite input and return to steady state conditions in a finite time. An unstable system will either diverge exponentially to positive or negative infinity or oscillate with ever increasing amplitude. A critically stable system will fall into a limit cycle of constant amplitude to any finite input. Oscillations in the order rates in supply chain are costly. Thus as a first step to dynamically designing a supply chain replenishment rule, we must ensure the replenishment rule is stable. Jury shows that the necessary and sufficient conditions for stability of a system is given by; A(1) > 0, (−1)n A(−1) > 0 and the matrices ∆± n−1 = Xn−1 ± Yn−1 are positive innerwise. Lets take each criteria in turn. A(1) must be greater than zero. A(1) = A(z)|z→1 , that is A(1) is given by 2.2 with the z’s replaced with 1. When the coefficients are substituted into A(1), A(1) = Tw . Thus this condition is satisfied if Tw > 0. (−1)n A(−1) must be greater than zero. In alike manner as above, A(−1) = A(x)|z→−1 , that is, Equation 2.2 with the z’s replaced with -1. (−1)n = −1 as the order, n, is 3 in our example. We substitute in the coefficients yielding (−1)n A(−1) = 2Ti Tw − Tw . Thus for stability we require 2Ti Tw − Tw > 0 ⇒ 2Ti Tw > Tw ⇒ 2Ti > 1 ⇒ Ti > 21 . ∆± n−1 = Xn−1 ±Yn−1 matrices are positive innerwise. A matrix is positive innerwise if its determinant is positive and all the determinants of its Inners are also positive. The ∆± n−1 = Xn−1 ± Yn−1 matrices here are; (4.1)

∆+ n−1

=

a3 a0

a2 + a0 a1 + a3

=

Ti Tw Tw − Ti

Tw (1 − Ti ) Ti Tw



(4.2)

∆− n−1 =

a3 −a0

a2 − a0 a3 − a1

=

Ti Tw Ti − Tw

2Ti − Tw (1 + Ti ) Ti Tw

7

We can see that both ∆± n−1 are 2x2, thus for our stability analysis here we only need to test whether the determinants of ∆± n−1 are positive as there are no Inners. The determinants are; (4.3)

|∆+ n−1 | = Tw (Ti (1 + Ti (Tw − 1) + Tw ) − Tw )

(4.4)

2 2 |∆− n−1 | = Ti Tw (3 + Ti ) − 2Ti + (Ti (Ti − 1) − 1)Tw .

Solving 4.3 for its roots yields a non-critical stability boundary of

(4.5)

Ti |nc

p −1 − Tw ± 1 − 2Tw + 5Tw2 = 2(Tw − 1),

which is completely dominated by the roots of 4.4,

(4.6)

p Tw2 − 3Tw ± Tw 1 − 2Tw + 5Tw2 Ti |c = 2(Tw − 2 + Tw2 ).

which describes the critical stability boundary. Investigation reveals that the value of Ti required for stability is given by (4.7)

(4.8)

(4.9)

Ti |(−2 0 that was given in the first step of Jury’s approach. We may also observe that when Ti → ∞ or Tw → ∞ then the numerator of 5.1 become infinite, giving us the last boundary of Figure 3, (Ti < ∞). 6. Aperiodicity via Jury’s Inners An aperiodic system has a time domain response with only a finite number of maxima or minima in its time domain response (less than n, the order of the system), Jury (1974) and Meerov and Jury (1998). Thus an aperiodic supply chain ordering policy will avoid costly long-term oscillations in the ordering patterns generated. This will obviously contribute to a reduction of any system induced cycles such as the well-known business cycle. An aperiodic condition occurs when all the roots of the characteristic equation are distinct and lie on the real axis in the interval [0, 1) in the z-plane, Jury (1974). Lets first consider the case of Equation 3.5, the generalised OUT policy when Ta = ∞ and Tw = Ti . The characteristic equation is 1 + Ti z − Ti that clearly has a single, simple root at z = TiT−1 . Thus the generalised OUT policy is aperiodic i when Tw = Ti and Ti = [1, ∞). The case of Equation 3.2, the generalised OUT policy with independent controllers, is more difficult to study as there are three roots to the characteristic



11

equation, of which two may be complex. In order to determine whether the roots are distinct, real and in the interval of [0, 1) in the z-plane we may apply the linear fractional transformation, 6.1, to the characteristic equation to map the segment [0, 1) in the z-plane onto the negative real axis in the s-plane. (6.1)

z=

s s−1

After this transformation is applied to Equation 3.2 we may use a continuous time stability test to determine the aperiodicity condition has now reduced to the requirement that the all the roots of A(s) are distinct and lie on the left-hand side of s-plane. (6.2)

A(s) = (s − 1)3 Tw + Ti (1 + s(s(2 + Tw ) − 3))

Collecting together terms surrenders the following coefficients of s; (6.3)

a0 = Ti − Tw , a1 = 3(Tw − Ti ), a2 = Ti (2 + Tw ) − 3Tw , a3 = Tw .

Jury (1974) and Meerov and Jury (1998) show us that the aperiodicity constraints are that ai > 0, ∀i = 0, .., n and that the Hurwitz matrix of A(s) must be positive innerwise. Lets first consider the signs of the coefficients ai , i = 0..n; (6.4) a0 > 0 ⇒ Ti > Tw , a1 > 0 ⇒ Tw > Ti , a2 > 0 ⇒ Ti >

3Tw , a3 > 0 ⇒ T w > 0 2 + Tw

We can show that the third order representation of the matched controller case, that is when Ti = Tw , contains two roots at the origin. From the requirement that the coefficients, ai , must be positive, we require that Ti > Tw and Tw > Ti simultaneously. This does not surprise us, as the double root is not distinct, contravening Meerov and Jury’s conditions for aperiodicity and re-enforcing the need to simplify first. The conditions for a positive a2 and a3 ensure that when Ti = [1, ∞) the distinct root of A(z) lies on the real axis of the z-plane between [0, 1). When Tw 6= Ti the double root at the origin will become either a pair of roots on the real line; one positive, one negative, or they will become complex. We may determine when the roots of A(z) become complex by considering when the Hurwitz matrix of A(s) as given by 6.5 is positive innerwise, (in fact Meerov and Jury (1998) assert that, in their experience, only the determinant of ∆5 is critical and needs to be considered and not it’s inners; we concur herein).  (6.5)

  ∆5 =   

a3 0 0 0 3a3

a2 a3 0 3a3 2a2

a1 a2 3a3 2a2 a1

a0 a1 2a2 a1 0

0 a0 a1 0 0

   .  

Substituting in the values of a0 ..an (the ai for A(s), not (A(z)) yields


12

STEPHEN M. DISNEY

Figure 5. Aperiodicity of the generalised OUT policy

(6.6)    ∆5 =   

Tw 0 0 0 3Tw

Ti (2 + Tw ) − 3Tw 3(Tw − Ti ) Ti − T w Tw Ti (2 + Tw ) − 3Tw 3(Tw − Ti ) 0 3Tw 2(Ti (2 + Tw ) − 3Tw ) 3Tw 2(Ti (2 + Tw ) − 3Tw ) 3(Tw − Ti ) 2(Ti (2 + Tw ) − 3Tw ) 3(Tw − Ti ) 0

0 Ti − Tw 3(Tw − Ti ) 0 0

and its determinant is, (6.7)

|∆5 | = 8Ti2 Tw − 6Ti Tw2 + 8Ti2 Tw2 − 12Ti Tw3 + 2Ti Tw3 .

Solving |∆5 | = 0 for Ti yields the following roots; (6.8)

Ti = 0, Ti = 0, Ti = Tw , Ti =

27Tw3 . (2 + Tw )2 (4Tw − 1)

Together with the positive innerwise requirement, 6.8, shows that the congruent 3 27Tw roots of A(z) are real if Tw > Ti > (2+Tw )2 (4T ; complex otherwise. So although, w −1) the OUT policy is not aperiodic when Ti 6= Tw as there is always one root with a neg3 27Tw ative real part, the system is ”weakly aperiodic” when Tw > Ti > (2+Tw )2 (4T . w −1) We use the term ”weakly aperiodic” to denote the case when the congruent roots of A(z) lie between (-1,1) on the real line and the simple root lies between [0,1) on the real line. We have collected these aperiodicity results graphically in Figure 5. We have verified our findings on the distribution of real roots of A(z) by plotting them in the complex z-plane for a selection of settings of parameters Ti and Tw in Figure 6. Here we can see that plots; E,G and H are aperiodic, B,C,D,I and


     


13

Figure 6. Verification of the poles of A(z)

Figure 7. Simulation of impulse responses J are weakly aperiodic and that plots A,F,K,L,M,N and O are non-aperiodic as predicted. For insights into the implications for the time domain solution of aperiodic systems we have further plotted the impulse response for the 15 sets of parameters of Figure 6 in Figure 7. From Figure 7 we can see that weakly aperiodic systems do not oscillate around zero. Strictly aperiodic systems approach zero exponentially directly from the demand impulse. Weakly periodic systems may initially diverge but soon return to exponential decay. Non-aperiodic systems possess a greater number of longer lasting oscillations and presumably these oscillations are costly in real supply chains.


14

STEPHEN M. DISNEY

Figure 8. Simulation to i.i.d. random demand

7. Simulation to a normal, independently and identically distributed demand pattern In order to gain an understanding of what our stability, bullwhip and aperiodic results might mean to a supply chain, Figure 8 shows six responses of the generalized OUT policy to an i.i.d. random demand. Plot A shows the demand pattern and the response produced by the classical OUT policy . They are the same because when Tw = Ti = 1 the OUT policy simply ”passes on orders”. This is the response that would have been obtained if a ”standard, text book” like approach was used to design a supply chain replenishment rule. Plot B shows a typical aperiodic response in the generalised OUT policy. Compared to the nearby non-aperiodic response (Plot C) we can see that there is a reduction in the fluctuations later in time. Clearly the amplitudes of the cylical demands in periods 25-30 has been attenuated. In Plot D we have set Tw = 99999 to effectively remove the WIP feedback loop. The plot shows a critically stable system. Here the order rate is stuck in an orbit around the mean demand with a large amplitude. Obviously this is a response to avoid in a real supply chain. Plot E highlights a weakly aperiodic response in the newly discovered stable region, that is when Tw < −2. Clearly this is a very damped response. The final Plot (F) demonstrates that the generalised OUT policy is also capable of level scheduling. In this policy Ti = Tw = 99999. This has the effect of removing both the inventory and the WIP feedback loops. Care must be taken with this level scheduling policy as there is no control of the inventory levels. In fact the inventory levels are simply an accumulation of the random demands and is thus a random walk. It is also possible to minimise the inventory and order variances in a trade-off as was shown in Disney, Towill and van de Velde (2004) for the case of Tw = Ti . Indeed, for the case where costs are equally related to inventory and order variance, the √ optimal Tw = Ti = 1+2 5 , the golden ratio for all lead-times. This was illustrated in plot B.



15

8. Conclusions Using Jury’s Inners approach we have studied the dynamics of a generalized OUT policy, a replenishment rule often used in supply chain to control material flow. To use the Inner’s approach we needed to specify a lead-time, we chose leadtime Tp = 2. Having to specify the lead-time appears to be the only weakness of Jury’s Inners approach. There appears to be no limit to the order of the transfer function when using Jury’s Inners. In other (time-based) methods, there is often a limit to the order of the transfer function that can be considered as inversion of a transfer function often requires us to find its roots; not an easy task and impossible in radicals for orders greater than 5. With the Inners approach we may avoid this issue altogether as we only require the determinants of matrices of the coefficients of the transfer function. This can always be done analytically, but the expressions may become large for very high order transfer functions. These high order transfer functions are generated when long lead-times or complex, multi-echelon supply chain models are present. Neither-the-less, computers are very good at dealing with matrix operations numerically if large, awkward determinants exist. A new stability test has been developed that is more direct than the traditional matrix based procedures as it simply involves solving 2an A(1)(−1)n A(−1)|∆− n−1 | = 0. This was achieved by integrating the bullwhip and stability analysis of Jury (1974). It appears not to generate redundant stability criteria. The stability analysis together with appropriate scaling of the parameter plane (by using 1/Ti and 1/Tw ), identified new areas of stable responses that were previously overlooked. We have also investigated the roots distribution of the transfer function in order to identify when aperiodic responses are generated by the OUT policy. The classical OUT policy (with Tw = Ti = 1) was found to be aperiodic. This means in a practical sense that the system does not oscillate unduly and has a damped response. An aperiodic response is also generated when Tw = Ti in the generalized OUT policy, thus confirming the importance of the so called Deziel and Eilon (1967) settings. We also found ”weakly aperiodic” responses in the parameter plane. Weakly aperiodic systems have only real roots, but one or more of them may be negative and they may have small short-term oscillations in the time domain, although they quickly dissipate and the long-term behaviour is damped. Aperiodicity is important in supply chains, as it will contribute to a reduction in rogue seasonality. The Inners approach is a very powerful tool for the supply chain analyst. Here we have used it to generate closed form bullwhip expressions for a system that is too complex to be solved using other techniques. For example, the ”squared impulse response” method of Disney and Towill (2003) is unable to solve the system considered herein. The Inners approach is also very efficient to use for both stability and bullwhip questions as many of the steps are similar in each analysis and we may gain some economies of scale. It is hoped this analytical approach will provide solutions to many more supply chain problems. In summary, the unique contributions of this paper are; the derivation of the bullwhip expression for the Tp = 2 case of the OUT policy with independent feedback controllers, the identification of stable responses with negative values of the WIP feedback controller (Tw ), a more direct stability test based on Jury’s approach, and the identification of two new supply chain performance characteristics - aperiodic responses and weakly aperiodic responses.


16

STEPHEN M. DISNEY

9. Biblography ˚ Astr¨ om, K.J. 1970. Introduction to stochastic control theory. Richard Bellman, University of Southern California. ˚ Astr¨ om, K.J., Jury, E.I., Agnail, R.G. 1970. A numerical method for the evaluation of complex integrals. IEEE Trans on Automatic Control AC-15 4 468-471 Aug 1970. Bertrand, J.W.M. 1986. Balancing production level variations and inventory variations in complex production systems. International Journal of Production Research. 24 5 1059-1074. Box, G.E. and Jenkins, G.M. 1970. Time series analysis forecasting and control. Holden-Day, San Francisco. Chatfield, C., Koehler, A.B., Ord, J.K. and Snyder, R.D. 2001. A new look at models for exponential smoothing. The Statistician 50 2 147-159. Chen, Y.F., Disney, S.M. 2003. The order-up-to policy sweet spot - Using proportional controllers to eliminate the bullwhip problem. EUROMA POMS Conference, Como Lake, Italy, 16-18 June 2 551-560. Dejonckheere, J., Disney, S.M., Lambrecht, M.R., Towill, D.R. 2003. Measuring the Bullwhip Effect: A control theoretic approach to analyse forecasting induced Bullwhip in order-up-to policies. European Journal of Operations Research 147 3 567-590. Dejonckheere, J., Disney, S.M., Lambrecht, M.R., Towill, D.R. 2004. The impact of information enrichment on the Bullwhip effect in supply chains: A control engineering perspective. European Journal of Operational Research, 153 727-750. Deziel, D.P., Eilon, S. 1967. A linear production - inventory control rule. Production Engineer 43 93-104. 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 echelon of a supply chain. European Journal of Operational Research 173 151-172. Disney, S.M., Grubbstr¨ om, R.W. 2004. The economic consequences of a production and inventory control policy. International Journal of Production Research 42 17 3419 - 3431. Disney, S.M., Towill, D.R. 2002. A discrete linear control theory model to determine the dynamic stability of Vendor Managed Inventory supply chains. International Journal of Production Research 40 1 179-204. 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 3 157-167. Disney, S.M., Towill, D.R. 2005. Eliminating inventory drift in supply chains. International Journal of Production Economics 93-94 331-344. Disney, S.M., Towill, D.R., van de Velde, W. 2004. Variance amplification and the golden ratio in production and inventory control. International Journal of Production Economics 90 3 295-309. Forrester, J. 1958. Industrial Dynamics, A major breakthrough for decision makers. Harvard Business Review, July-August 67-96. Gaalman, G. and Disney, S. M. (2007) On the echelon-order-up-to policy: A stochastic optimal control approach, Manufacturing and Service Operations Management Conference, Beijing, June 18-19.



17

Hoberg, K., Bradley, J.R. and Thonemann, U.W. 2007. Analyzing the effect of the inventory policy on order and inventory variability with linear control theory, European Journal of Operational Research 176 (3) 1620-1642. Hosoda, T. and Disney, S.M., (2006), ”On variance amplification in a threeechelon supply chain with minimum mean squared error forecasting”, OMEGA: The International Journal of Management Science 4 344-358 Jury, E.I., 1974, Inners and the stability of dynamic systems, John Wiley and Sons, New York. Lee, H.L., Padmanabhan, V., Whang, S. 1997a. The bullwhip effect in supply chains. Sloan Management Review Spring 93-102. Lee, H.L., Padmanabhan, V., Whang, S. 1997b. Information distortion in a supply chain: the Bullwhip Effect. Management Science 43 543-558. Lee, H.L., Padmanabhan, V., Whang, S. 2004. Comments on Information Distortion in a Supply Chain: The Bullwhip Effect. The Bullwhip Effects: A reflection, Management Science 50 12 1887-1893. Magee, J.F. 1956. Guides to inventory policy-II. Problems of uncertainty. Harvard Business Review, March-April, 103-116. Meerov, M.V., Jury, E.I. 1998. On periodicity robustness. International Journal of Control. 70 2 193-201. Metters, R. 1997. Quantifying the bullwhip effect in supply chains. Journal of Operations Management 15 89-100. Nise, N.S. 1995. Control systems engineering. The Benjamin Cummings Publishing Company, Inc., California. Riddalls, C.E., Bennett, S. 2002. The stability of supply chains. International Journal of Production Research 40 459-475. Sterman, J. 1989. Modelling managerial behaviour: Misperceptions of feedback in a dynamic decision making experiment. Management Science 35 3 321-339. Towill, D.R. 1982. Dynamic analysis of an inventory and order based production control system. International Journal of Production Research 20 6 671-687. Warburton, R.D.H., Disney, S.M., Towill, D.R., Hodgson, J.P.E. 2004. Further insights into The stability of supply chains. International Journal of Production Research 42 3 639-648. Vassian H.J. 1955. Application of discrete variable servo theory to inventory control. Journal of the Operations Research Society of America 3 272-282. Cardiff Business School, Cardiff University, Aberconway Building, Colum Drive, Cardiff, CF10 3EU, UK. E-mail address: [email protected]. Tel: +44(0)2920 876310. Fax: +44(0)2920 874301


SUPPLY CHAIN APERIODICITY, BULLWHIP AND

SUPPLY CHAIN APERIODICITY, BULLWHIP AND

Suggest Documents

Removing bullwhip from the Tesco supply chain

bullwhip effect and supply chain modelling and analysis ... - CiteSeerX

Supply Chain Management Practices and Scope for Bullwhip ... - ebsco

Analysis of bullwhip effect on supply chain with Q ...

The impact of the supply chain structure on bullwhip ...

Philips Electronics Synchronizes Its Supply Chain to End the Bullwhip ...

THE BULLWHIP EFFECT IN SUPPLY CHAIN ... - ERTEK PROJECTS

Information distortion in a supply chain: The bullwhip effect

SUPPLY AND PRODUCTION NETWORKS: FROM THE BULLWHIP ...

Supply Chain Vs. Supply Chain

The bullwhip-effect in the electricity supply

Supply Chain Strategy - UPS Supply Chain Solutions

Introducing Supply Chain Management - Logistics & Supply Chain ...

Supply Chain and Finance

Supply Chain

Supply-Chain

supply chain evolution and green supply chain perspective

Supply Chain Management Practices and Supply Chain Performance ...

Supply chain design and design for supply chain

The bullwhip effect in a close loop supply chain - System Dynamics ...

Study of a supply chain system with the presence of bullwhip effect

Controller Design and Reduction of Bullwhip for a General Supply ...

Purchasing and Supply Chain Management

strategy and supply chain management