A supply chain network might exhibit complex dynamics in the face of increasingly volatile and
uncertain environment. The impact of network structure and collaboration on the dynamics and robustness of supply chain network, however, remains to be explored. In this paper, a unified state space model for a two-layer supply chain network composed of multiple distributors and multiple retailers is developed. A robust control algorithm is advocated to reduce both order and demand fluctuations for unknown demand. Numerical simulations demonstrate that the robust control approach has the advantage to reduce both inventory and order fluctuations. In the simulation experiment, it is interesting to notice that complex network structure and collaborations might contribute to the reduction of inventory and order oscillations. This paper yields new insights into the overestimated bullwhip effect problem and helps us understand the complexities of supply chain networks.
1. Introduction
A supply chain network is referred to as a complex network of organizations that synchronizes a series of interrelated business processes, such as procurement, manufacturing, distribution, and retailing, to create values to final customers in the form of products and services [1]. In reality, supply chain network is facing more and more fluctuating market environment and the customer demand is usually badly uncertain due to a lot of factors, such as technology revolution, short product cycles, and promotions. These uncertainties might propagate and further lead to unexpected dynamic behaviors, for example, instability [2–5], bullwhip effect [6–9], and even chaos [10–12], which can further cause huge extra costs due to order and inventory fluctuations [13]. Therefore, it is increasingly important to understand the dynamics and improve the robustness before designing a supply chain network.
It has been found that supply chain network has scale-free and small-world features [14]. Note that the small-world feature and scaleless feature have been widely explored in complex network theory [15–17]. These features indicate that a handful of nodes might dominate the whole network and potentially enhance the capability to respond to uncertainties. In addition, determining the number of upstream suppliers and extending distribution channels to gain more market shares are also important topics in supply chain management. All these facts demonstrate that the topology of supply chain network has substantial impacts on overall performance. In fact, the supply chain network studied in this paper is also a networked control system. Similar to the literature stream of networked control system [18–20], the network structure is a significant factor that affects the control performance. In addition, the central task for the supply chain network or the networks studied in the previous literature is how to design an appropriate controller to ensure stability and satisfy specific requirements. It is worthwhile to notice that a supply chain network is usually composed of heterogeneous nodes, for example, retailers, distributors, manufacturers, and raw material suppliers. In a supply chain network, distributors usually exchange information, for example, placing and receiving orders, with nodes inside the network. Unlike distributors, retailers, however, must fulfill customer demands from external volatile environment. For this reason, the dynamics of heterogeneous nodes should be treated or described differently. Furthermore, there exist a lot of interaction activities among nodes in supply chain network, for example, cooperations and competitions. As a result, all these features, complex topologies caused by complex supply relationships, heterogeneous nodes, and interactions, make supply chain network increasingly complex, hard to describe, and challenging for exploring dynamics and robustness.
It is not surprising that the dynamics and robustness of supply chain system have received considerable attention over the past decade. Control engineering technologies are frequently used to analyze stability [2, 21, 22] and study the bullwhip effect [7, 8, 23–25], which refers to the amplification of order fluctuations as one moves up a supply chain from downstream to upstream [6, 26]. In particular, robust control methods are well used in analyzing the bullwhip effect for known inventory rules [9, 27], or designing inventory rules for unknown demand using feedback control method [28]. Some “fluid models” and simulation models are also proposed to study the dynamic properties of supply chain systems [3, 10, 29, 30].
The majority of the previous studies focus on serial supply chains. Most supply chain systems, however, have a network topology in which each entity’s decision may be influenced directly or indirectly by other nodes. More importantly, demand may arise from many customers at different locations of the supply chain network, and the interactions among nodes should be addressed. Unfortunately, few studies are undertaken to address the dynamics, robustness, and the impact of collaborations on dynamics in the context of supply chain network. For example, Helbing [31] proposes a rather general model of supply chain network and connected it to queueing theory and macroeconomics. Scholz-Reiter et al. [32] presented an approach to calculate the stability condition by mathematical system theory to guarantee stability for production networks, to identify a stability region, and to refine the region by discrete event simulations. Ouyang and Li [33] study the bullwhip effect propagation in a supply chain network with linear and time-invariant inventory management policies using a frequency domain control method. They obtained conditions to predict the existence of bullwhip effect for unknown demand. Recently, Dominguez et al. [34] study the bullwhip effect of a divergent supply chain system with simulation experiment. They begin to notice that complex structure of supply chain might yield different bullwhip effect conclusions.
From the above, the dynamics of supply chain network has attracted the attention in the very recent years. The impact of network structure and collaboration on supply chain dynamics still remain unclear. The contributions of this paper are as follows. A unified state space model is developed to model the dynamic for a supply chain network composed of heterogeneous nodes. To improve the robustness for unknown demand, a robust H∞ control approach is advocated to obtain optimal inventory rules. The time domain robust control method facilitates us to obtain the optimal replenishment policy under well-defined performance measures, in order to reduce both order and inventory fluctuations. There is a trade-off between reducing inventory variations and mitigating the bullwhip effect [23]. The effectiveness and superiority of the proposed method are numerically validated for several demand processes. Furthermore, we incorporate the network topology and demand collaboration between retailers into the dynamic model, which enables us to numerically study the impacts of network structure and demand collaboration on dynamics improvement. We show that both demand collaboration and complex supply relationship might contribute to dynamics improvement. The results obtained in this research provide new insights into supply chain management, especially for the overestimated bullwhip effect problem [35, 36], via considering complex supply relationship and demand collaboration.
The structure of this paper is arranged as follows. Section 2 describes the model and Section 3 introduces a robust metric and the optimization algorithm to calculate the minimal value of the robust metric. In Section 4, we numerically validate the effectiveness of the robust control method and also study the influences of network structure and demand collaboration on dynamics and robustness.
2. Modeling the Supply Chain Network
Consider a two-layer supply chain network constituted with m distributors and n retailers. The structure of supply chain network is shown in Figure 1. Without loss of generality, we impose the following assumptions on the supply chain network:
All the facilities, including retailers and distributors, have the same events order in each period, in order to construct a discrete-time state space model. They place replenishment orders at the beginning of each period and fulfill downstream orders or demands during the remaining time.
Each retailer might choose to place orders from many distributors and each distributor might satisfy the order from multiple retailers. All the distributors place orders from external suppliers.
The customer demands of retailers that come from external market are uncertain. All the nodes in the supply chain network follow the same way to deal with shortages; namely, unfulfilled orders or demands will be backlogged.
Both the distributors and retailers have fixed replenishment lead times.
The structure of a general two-layer supply chain network.
Denote symbols i,1≤i≤m, and j,1≤j≤n, as the ith distributor and the jth retailer, respectively. Then the supply relationship between retailers and distributors could be defined by a supply relationship matrix L=(lij)m×n∈Rm×n. If lij=1, the distributor i and the retailer j establish the supply relationship between them. Otherwise, if lij=0, there is no relationship between them. In fact, the supply relationship matrix reflects the network topology. For example, if each column of the matrix L has only one element with value 1, then each retailer has only one immediate supplier. In contrast, if lij=1(∀i,j), then each distributor supplies goods to every retailer, which represents the most complex network structure.
Suppose that the retailers in the same market can collaborate with each other, and consider the collaboration mechanism via demand diversion in [33]. Let dj(t) be the demand received by retailer j. Then the amount δj1j2×[dj1(t)-dj2(t)] of the received demand of retailer j1 will be satisfied by retailer j2 if dj1(t)>dj2(t). Note that the parameter δj1j2 can be negotiated by each pair of retailers by considering some factors such as geographic location and demand characteristics. As a result, we can define a collaboration matrix H≜(δj1j2)n×n∈Rn×n to exactly describe the demand diversion among retailers in the supply chain network. Assume that the demand of retailer j1 is dj1(t) and the demand of retailer j2 is dj2(t). The essential demand of retailer j1 after demand collaboration becomes dj1(t)+δj1j2[dj2(t)-dj1(t)]. The demand of retailer j2 is changed to dj2(t)+δj2j1[dj1(t)-dj2(t)]. Then we must have δj1j2=δj2j1 without losing any customer demand. In other words, the collaboration matrix is symmetric in order to balance the shifted demand between each pair of retailers.
As the supply chain network is composed of heterogeneous nodes, the dynamic equations for the retailers and distributors should be described separately. We will now first consider the inventory dynamics for distributors. Let wi(t), 1≤i≤m, denote the inventory position of the distributor i, let ui(t) denote the placed order by distributor i from external suppliers, and let uij(t) denote the received order of distributor i from retailer j(1≤j≤n). Then, the resulting dynamics of the distributor i can be represented by(1)wit+1=wit+uit-∑j=1nlijuijt.
All the retailers not only place orders from their upstream distributors but also fulfill external customer demand. In addition, a retailer might divert its received demand to its collaboration partner. Let rj(t), 1≤j≤n, denote the inventory position of retailer j; then its balanced inventory equation could be obtained as(2)rjt+1=rjt+∑i=1mlijuijt-djt-∑j1=1nδjj1dj1t-djt.It is worthwhile to note that as we select inventory position as the state variable, the lead time disappeared in the dynamics equation. In fact, inventory position is a comprehensive variable that incorporates on-hand inventory, backorder, and on-road inventory simultaneously.
In order to develop a unified state space model, we will now transform the supply relationship matrix L and the collaboration matrix H. By doing this, we can incorporate the network structure and collaboration among retailers into the dynamic model. Let Lj=[l1j,l2j,…,lmj]′, 1≤j≤n, define the jth column of the matrix L. Then we could define the following two matrixes:(3)M1≜diagL1′,L2′,…,Ln′,M2≜G1G2⋯Gn,where Gj=diag(l1j,l2j,…,lmj), 1≤j≤n. We see that M1∈Rn×mn and M2∈Rm×mn. Based on the collaboration matrix H, we could obtain the following matrix:(4)HL≜1-∑j≠1δ1jδ12⋯δ1nδ211-∑j≠2δ2j⋯δ2n⋯⋯⋯⋯δn1δn2⋯1-∑j≠nδnj.
The demand rate of retailer j is assumed to take the form(5)djt=dje+ξjt,where dje is a known constant and ξj(t) denotes the energy bounded disturbance with mean 0. The demand in Ouyang and Daganzo [9] and Ouyang and Li [33] also has the characteristics of the above demand pattern. If the supply chain network is stable and dj(t)=dje, j=1,2,…,n, then the inventory positions and order quantities of retailers and distributors will return to constant or steady values after finite periods. Define wie and rje as the steady values of inventory position of distributor i and retailer j, respectively. Let uie and uije denote the steady values of the amount of order of the distributor i placed from external suppliers and the amount of order of the retailer j placed from distributor i, respectively. In such case, we can have limt→∞wi(t)=wie,limt→∞rj(t)=rje,limt→∞ui(t)=uie, and limt→∞uij(t)=uije.
It is well-known that inventory rules play significant roles in determining supply chain dynamics. In fact, inventory rules are essentially feedback controllers using state information, for example, inventory position. Define the following state vector: x~(t)=w~(t)r~(t)′, where (6)w~t=w~1t,w~2t,…,w~mt=w1t-w1e,w2t-w2e,…,wmt-wme,r~t=r~1t,r~2t,…,r~nt=r1t-r1e,r2t-r2e,…,rnt-rne.The input vector can be represented by (7)d~t=d~1t,d~1t,…,d~nt′=d1t-d1e,d2t-d2e,…,dnt-dne′∈Rn.Before defining the control vector, we will now introduce the operator Vec(A).
Definition 1.
Assume that A∈Rm×n is a matrix with m rows and n columns and A=(A1,A2,…,An), where Ai is the ith column of A. Then we can produce a single column vector using the following operator, which is given by(8)VecA=A1′,A2′,…,An′′∈Rmn.
For the supply chain network, the control vector u~(t) can be represented as u~(t)=u~w(t)u~r(t)′∈Rm+mn, where (9)u~wt=u~1t,u~2t,…,u~mt=u1t-u1e,u2t-u2e,…,umt-ume,u~rt=VecMRt′,where (10)MRt=u~11tu~12t⋯u~1ntu~21tu~22t⋯u~2nt⋯⋯⋯⋯u~m1tu~m2t⋯u~mnt=u11t-u11eu12t-u12e⋯u1nt-u1neu21t-u21eu22t-u22e⋯u2nt-u2ne⋯⋯⋯⋯um1t-um1eum2t-um2e⋯umnt-umne∈Rm×n.
Define y~(t)∈R2m+n+mn as system output and let Ix denote the unit matrix with x dimensions. We then can develop the following unified state space model:(11)x~t+1=x~t+Bu~t+Bwd~t,y~t=Cx~t+Du~t,where(12)B=Im-M20M1∈Rm+n×m+mn,Bw=0-HL∈Rm+n×n,C=χiIm+n000∈R2m+n+mn×m+n,D=000χoImn+1∈R2m+n+mn×m+mn,in which χi and χo are variables with either value 1 or value 0 to define system output according to practical situations. For example, if managers only care for inventory fluctuations, they can set χi=1 and χo=0. Otherwise, the setting of χi=0 and χo=1 can be used to mitigate bullwhip effect or order fluctuations. Both the inventory and order dynamics will be considered if χi=1 and χo=1. Based on the state space model, we can design inventory rules to improve the dynamics and robustness of the supply chain network for unknown demand with robust H∞ control method.
3. The Robust Performance Metric and <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M107"><mml:mrow><mml:msub><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>∞</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> Control Method
The bullwhip effect, one kind of typical dynamic behavior, is of great significance and has been investigated extensively. The robustness of bullwhip effect, however, has been addressed recently only in a few literatures [9, 27, 28]. Meanwhile, most of the existing studies on the bullwhip effect pay major attention to order fluctuations but inventory fluctuations are neglected. To this end, this research focuses on improving the dynamics and robustness of supply chain system for unknown demand with a robust H∞ control approach with the aim to reduce both inventory and order fluctuations.
Based on the state space model (11), we will now consider the following robust metric:(13)WI≜sup∀d~t≠0∑k=0∞y~t′y~t∑k=0∞d~t′d~t1/2.We see that the metric WI is the amplification metric of the worst-case root mean square errors (RMSE) when the customer demand is unknown. Therefore, minimizing the value of WI can effectively improve the dynamics and robustness of supply chain system. It is worthwhile to mention that this metric, also called L2 gain or H∞ norm, has been widely used in the field of control engineering. Although it is not an economic metric, the metric WI is intimately related to fluctuations of order and inventory in supply chain systems. For example, Boukas [28] pointed out that reducing the value of the metric WI increases the speed of a simple production control system to respond to customer demand. Ouyang and Daganzo [9] directly used this metric to quantify the bullwhip effect for unknown demand.
From the form of formula (13), we see that reducing the value of metric WI reduces the fluctuation of system output for unknown input. Thus, we can reduce order oscillations if we select placed orders as system output by setting χi=0 and χo=1. However, in most cases, inventory fluctuations and order fluctuations should be balanced. Unlike the traditional bullwhip effect metric, which mainly focuses on dampening order fluctuations, the metric WI can take both order fluctuations and inventory fluctuations into account by setting χi=1 and χo=1. In this case, the exact formula of WI for our supply chain network model becomes(14)WI=sup∀d~t≠0∑t=1∞∑i=1mw~it2+u~it2+∑j=1nr~jt2+∑i=1mu~ijt2∑t=1∞∑j=1nd~jt21/2.
As emphasized, minimizing the value of the metric WI helps improving dynamics for unknown demand. Here we focus on designing appropriate inventory rules to achieve this goal. In this research, the inventory rule for the supply chain network can be represented in the form of state feedback controller, which is given by u~(t)=Kx~(t). Then the robust H∞ control problem of our supply chain network model is to seek a feedback controller u~(t)=Kx~(t) that satisfies the following optimization problem:(15)min∀u~t=Kx~tWI.The inventory rule to minimize the value of the robust metric WI for the state space model (11) can be obtained by Theorem 2. The proof process refers to Kim and Park [37].
Theorem 2.
Consider the discrete-time system (11). If there exists a positive-definite matrix Q and a matrix M satisfies the following optimization problem:(16)minγ1/2subject to(17)-QQ+BMBw0Q+M′B′-Q0QC′+M′D′Bw′0-γIn00CQ+DM0-I2m+n+mn<0,γ>0,then system (11) is quadratically stable and the minimal value of the robust metric WI can be approximated with the minimized γ1/2, with the state-feedback controller: u~(t)=MQ-1x~(t).
The importance of Theorem 2 lies in two aspects. Firstly, the stability of the whole supply chain network is warranted if the optimization problem in Theorem 2 is solvable. Secondly, it produces an inventory rule that minimizes the metric WI, or the H∞ norm, for unknown demand to reduce output oscillations. It also should be noted that the optimization problem can be readily solved using the standard numerical software, MATLAB, which enables us to numerically study the impacts of demand collaborations and network structures on dynamics and robustness.
4. Numerical Analysis
As introduced, the proposed robust H∞ control method is expected to reduce both order and inventory fluctuations, which will be numerically validated in the following section. Meanwhile, as the network topology and collaboration among retailers have been incorporated into the unified state space model, it enables us to study the influences of network structures and collaborations on the improvement of dynamics and robustness.
Traditionally, inventory and order fluctuations haven been measured with variance or standard deviation. For example, the bullwhip effect can be quantified with the ratio of variance between demand and order [26]. Here we consider the following four metrics in the numerical analysis:(18)Wb1=∑i=1mVaru~it,Wb2=∑i=1m∑j=1nVaru~ijt,Ws1=∑i=1mVarw~it,Ws2=∑j=1nVarr~jt,in which Var(·) is used to calculate the variances of time sequences. It is quite apparent to see that the metrics Wb1 and Wb2, respectively, represent the order fluctuations at distributor layer and retailer layer. The latter two metrics Ws1 and Ws2 are used to measure the inventory oscillations for retailers and distributors.
4.1. Effectiveness Validation for the Robust Control Method
Consider a supply chain network composed of 3 retailers and 3 distributors. Without loss of generality, we assume that there exists no collaboration among retailers and each retailer places order from all the three distributors. The supply relationship matrix and collaboration matrix can be obtained as(19)L=111111111,H=000000000.
Assume that the demand of each retailer follows the first-order autoregressive process:(20)d~jt+1=ρd~jt+1-ρεt,j∈1,2,3,in which ρ is the autoregressive parameter and εt is the white noise. Setting different values for χi and χo, we can consider three different control ways:
control way 1: to reduce inventory fluctuation (χi=1, χo=0);
control way 2: to reduce order fluctuation (χi=0, χo=1);
control way 3: to reduce both order and inventory fluctuations (χi=1, χo=1).
Using the optimization algorithm in Theorem 2, we can obtain three different inventory rules for three control ways. Take the third control way as an example. In such a parameter setting, we have χi=1,χo=1. As we consider a supply chain network composed of 3 retailers and 3 distributors, we could have n=m=3. Then, the matrices M1 and M2 can be obtained as (21)M1=111000000000111000000000111,M2=100100100010010010001001001.Because we do not consider demand collaboration in this section, the matrix HL will be calculated as HL=diag{1,1,1}. From the value of M1, M2, and HL, we can further easily calculate the matrices B, Bw, C, and D. The inventory rule, a state-feedback controller, can be solved using Theorem 2 with LMI Toolbox in MATLAB, which is represented as (22)u~t=-0.4140-0.1911-0.1911-0.2656-0.2656-0.2656-0.1911-0.4140-0.1911-0.2656-0.2656-0.2656-0.1911-0.1911-0.4140-0.2656-0.2656-0.26560.1450-0.0779-0.0779-0.3201-0.0058-0.0058-0.07790.1450-0.0779-0.3201-0.0058-0.0058-0.0779-0.07790.1450-0.3201-0.0058-0.00580.1450-0.0779-0.0779-0.0058-0.3201-0.0058-0.07790.1450-0.0779-0.0058-0.3201-0.0058-0.0779-0.07790.1450-0.0058-0.3201-0.00580.1450-0.0779-0.0779-0.0058-0.0058-0.3201-0.07790.1450-0.0779-0.0058-0.0058-0.3201-0.0779-0.07790.1450-0.0058-0.0058-0.3201x~t.The inventory and order dynamics of the supply chain network with ρ=0.6 are shown in Figures 2 and 3. With a simulation length of 15000 periods and by using a standard numerical software, we can calculate the values of the four metrics, Wb1, Wb2, Ws1, and Ws2, which are summarized in Table 1.
System performance under different control ways.
Wb1
Wb2
Ws1
Ws2
Control way 1
3.5782×103
6.2078×103
1.9027×10-27
0.7570
Control way 2
9.8163×10-5
2.4236×10-4
0.4375×103
1.1454×103
Control way 3
0.2147
0.0810
0.0141
0.7890
Dynamic curves of inventory under three different control ways.
Inventory curves under the first control way
Inventory curves under the second control way
Inventory curves under the third control way
Dynamic curves of order under three different control ways.
Order curves under the first control way
Order curves under the second control way
Order curves under the third control way
From Table 1 and Figures 2 and 3, we see that the first control way has an obvious effect on reducing inventory fluctuations. However, it inevitably results in huge order fluctuations. The second control way, in contrast, behaves well in reducing order fluctuations at the expense of high inventory fluctuations. We see that there is an obvious trade-off between improving inventory dynamics and damping order fluctuations. In fact, balancing inventory fluctuations and order fluctuation has been recognized as an important problem in the literature [23].
In practice, both inventory and order fluctuations greatly enhance difficulties for supply chain management. For example, order fluctuations bring difficulties for arranging labor force and making full use of production capacity. Meanwhile, inventory fluctuations should be rigidly controlled for distributors because high costs will incur if items cannot be sold out timely. It will be better if we can find inventory rules that reduce both inventory and order fluctuations. From the third row in Table 1, we see that the third control way performs well in reducing both inventory and order fluctuations which have been decreased. Figures 2(c) and 3(c) further confirm the advantage of the third control way with the demand process (20). In this sense, it reveals that the robust H∞ control method is highly effective in improving dynamics and robustness.
4.2. Impact of Network Structure on Dynamics and Robust Performance
The structure of supply chain network can be altered as a result of a lot of activities, such as adding new suppliers and extending distribution channels. These activities certainly impact the dynamics and robustness of supply chain network. For this purpose, the following study will be conducted with two concerns. We firstly consider the supply relationship, which is described by the supply relationship matrix L, and then consider the change of the number of nodes. Assume that χi=1 and χo=1 to consider both inventory and order fluctuations. To concentrate our attention on network structure, we will not take demand collaboration into account.
4.2.1. The Impact of Supply Relationship on System Performance
It is assumed that each retailer has the flexibility to choose its suppliers. Each retailer can place orders from multiple distributors; each retailer can also choose a single distributor as its supplier. Here we consider a supply chain network composed of 3 retailers and 6 distributors. Four types of supply are considered as shown in Figure 4. The first structure is the simplest one, in which each retailer has only one supplier. The second structure is a distribution system with a major supplier. The latter two structures represent supply chain network with more complex supply relationships. The fourth structure, in which each retailer places orders from all the distributors, is the most complex structure in this study. The supply relationship matrix L can be varied to represent different structures. The supply relationship matrixes for the four types of structures are given by (23)L1=110000001100000011,L2=100000111111000001,L3=101101110110011011,L4=111111111111111111,where L1, L2, L3, and L4 are the supply relationship matrices for structure 1, structure 2, structure 3, and structure 4, respectively.
Four types of supply chain network structures.
Using the optimization method introduced in Theorem 2, we could accordingly compute the minimal values of the robust metric WI for the supply chain network under four different structures; they are 2.2396, 2.2691, 1.9519, and 1.8361. It is obvious that the structure of supply chain network has significant influences on dynamics and robustness.
We use the demand process (20) with ρ=0.6 to generate demand sequences for retailers to observe the dynamics of order and inventory. The order dynamics of the three distributors are illustrated in Figure 5. Figure 5 shows that, compared with structure 1 and structure 2, structure 3 and structure 4 improve the dynamics for the three distributors. It means that complex supply chain network structure might help reduce order dynamics, namely, reducing the bullwhip effect. However, in the practice, we should consider other factors, such as distance between two nodes and information sharing costs. As the distance is increased, the transportation cost might be greatly increased. The second structure, which reflects a multiple-retailer system, results in high order fluctuations. In fact, too many downstream nodes can increase the risk of shortage, which has been recognized as one of the main causes of bullwhip effect [6], in the face of limited supply resources. Furthermore, since all the retailers should satisfy uncertain demands from external environment, too many retailers largely increase the uncertainties for the whole network and further lead to poor performance.
Order dynamics under different structures of supply chain network.
Order dynamics of distributors under structure 1
Order dynamics of distributors under structure 2
Order dynamics of distributors under structure 3
Order dynamics of distributors under structure 4
4.2.2. The Impact of the Number of Nodes on Robust Performance
As the supply chain network is hierarchical, we can change the number of upstream nodes or downstream nodes. Here we consider two cases:
Case 1: increasing the number of distributors.
Case 2: increasing the number of retailers.
We assume that there exist two retailers in the first case and two distributors in the second case. The first case characterizes the supply chain system with multiple suppliers as the number of distributors increases. The second case represents distribution system with multiple retailers. It is obvious that increasing the number of distributors, as shown in Figure 6(a), has a positive effect on robustness improvement. This result can be explained where increasing the number of suppliers provides new purchasing channels or opportunities and thus the possibility of shortages will be lowered. From Figure 6(b), we see that increasing the number of retailers results in poor robust performance. This result confirms the result obtained before that a distribution system usually behaves poor on robustness performance.
Impact of the number of nodes on robust performance.
Impact of distributor number on system performance
Impact of retailer number on system performance
From the above, it reveals that the structure of supply chain network is closely related to dynamics and robustness and thus seeking for optimal structures is essential for supply chain networks to achieve high economic benefits. Most importantly, our study adds new insights into the overestimated bullwhip effect problem [35, 36]. A recent study of US industry level data found that demand volatility does not increase as one moves up the supply chain [35]. The author observed that, in general, manufacturers do not have substantially greater demand volatility than retailers. However, traditional researches on the bullwhip effect believe that order fluctuations will be amplified stage by stage [6, 26]. In other words, manufacturers, probably the most upstream nodes in a supply chain network, should experience the worst order fluctuations. It seems that the findings in [35] contradict the results in traditional researches. In this research, we found that a supply chain network with complex structures behaves well in dynamics and robustness compared to simple networks. We also found that increasing the number of distributors contributes to the bullwhip mitigation. For this reason, we argue that the existing approaches that aim at quantifying the bullwhip effect neglect the network structure of supply chains. We believe that if we begin to study the bullwhip effect of supply chain network with the consideration of complex supply relationships and collaborations, not just serial supply chains, new results about the bullwhip effect might emerge.
4.3. Impact of Demand Collaboration on Robust Performance
As introduced, the collaboration mechanism among retailers is characterized by the collaboration matrix H. The retailer j1 and retailer j2 form their collaboration relationship with the parameter δj1j2, 0≤δj1j2≤1, which reflects the strength of collaboration. We will now consider a supply chain network with 3 retailers and 6 retailers. In many circumstances, one retailer tends to collaborate with a single retailer geographically close to it. To facilitate analysis, we assume that each retailer has only one collaboration partner. Then the corresponding collaboration matrix is represented as(24)H=0δ120000δ2100000000δ340000δ4300000000δ560000δ650,in which δ12, δ34, and δ56 are the collaboration parameters established by each pair of retailers. In addition, we assume that each retailer places order from all the suppliers. In other words, the values of all the elements of matrix L are 1. We also assume that χi=1 and χo=1 in order to consider overall performance. In this study, we will discuss three different collaborating scenarios:
scenario 1: δ12=δ34=δ56=0 (without cooperation);
scenario 2: δ12=δ34=δ56=0.25 (mild cooperation);
scenario 3: δ12=δ34=δ56=0.5 (close cooperation).
Assume that retailer j1 collaborates with retailer j2, ∀(j1,j2)∈{(1,2),(3,4),(5,6)}; then their demand processes are assumed to be(25)d~j1t+1=ρd~j1t+1-ρεj1t+ηj1ɛt,d~j2t+1=ρd~j2t+1-ρεj2t+ηj2ɛt,where ρ=0.2, ηj1, ηj2∈{-1,0,1}, with white noises εj1(t) and εj2(t) with mean 0 and variance 1. A variety of auto- and cross-correlations among the demands can be represented by the selection of parameters. For example, the two demands may be independent of each other (ηj1=ηj2=0), positively correlated (ηj1=ηj2=1), or negatively correlated (ηj1=-ηj2=1).
The minimal values of the robust metric WI under the above three different collaboration scenarios can be computed by Theorem 2, which are 1.8361, 1.8686, and 1.8751. It shows that the collaboration behaviors have little impacts on the robustness improvement. However, we will see that the collaboration behaviors among retailers may mitigate the bullwhip effect for some special demand processes, as demonstrated in the following simulations. Under the demand processes formulated in (25), we can obtain the values of the four metrics, Wb1, Wb2, Ws1, and Ws2, via simulation with a length of 15000 periods, which are shown in Table 2. From Table 2, we see that the demand collaboration between retailers, or demand collaboration, has obviously reduced order and inventory fluctuations. Interestingly, when the values of collaboration parameters increase, the fluctuations of inventory and order decrease, and the reduction is the most significant for the case of negative cross-correlation. In addition, the reduction of inventory and order fluctuations for retailers is more obvious than distributors. This can be explained by the fact that the collaborations are established among retailers, not distributors.
System performance under different collaboration scenarios.
Wb1
Wb2
Ws1
Ws2
Scenario 1
ηj1=ηj2=0
1.0654
0.3993
0.0432
4.1398
ηj1=ηj2=1
4.3605
1.0635
0.1767
10.5216
ηj1=-ηj2=1
1.0423
1.0057
0.0422
10.5494
Scenario 2
ηj1=ηj2=0
1.0352
0.2482
0.0436
2.5674
ηj1=ηj2=1
4.3898
0.9023
0.1849
9.0700
ηj1=-ηj2=1
1.0253
0.3944
0.0431
4.2528
Scenario 3
ηj1=ηj2=0
1.0089
0.2084
0.0445
2.0619
ηj1=ηj2=1
4.1481
0.8425
0.1827
8.4610
ηj1=-ηj2=1
1.0093
0.2055
0.0445
2.0680
In summary, collaborations among retailers might not contribute to robustness improvement; however, collaborations can reduce inventory and order fluctuations for some specific demand processes, especially for retailers. When [33] discussed the collaboration among retailers, they obtained similar results, in which they focus on order fluctuations in term of bullwhip effect. In this sense, our results confirm and also enrich the results in [33] using the robust H∞ control.
5. Conclusions
This paper has investigated the dynamics and robustness of a two-layer supply chain network constituted with multiple retailers and multiple distributors. A unified state space model is developed to characterize some basic characteristics of the supply chain network. Specifically, we describe different supply relationship with supply relationship matrix and consider one kind of demand collaboration mechanism among retailers, namely, demand diversion. As the nodes of supply chain network are heterogeneous, the dynamic equations for retailers and distributors are developed separately. Considering these characteristics enables us to investigate the influences of network structure and collaborations on dynamics and robustness.
It is well known that reducing order fluctuations usually increases inventory fluctuations [23], which tends to result in shortages and high levels of safety stock. This research, however, numerically showed that the robust H∞ control method based on the state space model has the advantage to reduce both inventory and order fluctuations for supply chain network. The simulation results show that complex supply relationships and demand collaborations among retailers contribute to the reduction of inventory and order fluctuations for unknown demand or specific demand. Collaborations among retailers might not improve the robustness of a supply chain network. However, collaborations can reduce inventory and order oscillations for some specific demand processes, especially for retailers. We argue that if we begin to study the bullwhip effect of supply chain network with the consideration of complex supply relationships and collaborations, new results about the bullwhip effect might emerge.
Finally, it is worthwhile to point out that the investigation of the dynamics and robustness of supply chain network, as we believe, is at the beginning stage. A lot of problems deserve future research. For example, if we consider more complex inventory rules, time delay problem will emerge. The relationship between system costs and dynamics is also an important topic because a supply chain network with complex supply relationships might incur high transportation costs.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (no. 71401181 and no. 71502076) and the National Natural Science Foundation for Distinguished Scholars of China (no. 71125001). This work was also supported by MOE (Ministry of Education in China) Project of Humanities and Social Sciences (no. 14YJC630136 and no. 15YJC630008).
LiG.JiP.SunL. Y.LeeW. B.Modeling and simulation of supply network evolution based on complex adaptive system and fitness landscapeDisneyS. M.TowillD. R.A discrete transfer function model to determine the dynamic stability of a vendor managed inventory supply chainDaganzoC. F.On the stability of supply chainsWeiY. C.WangH. W.QiC.On the stability and bullwhip effect of a production and inventory control systemWeiY.WangH.QiC.The impact of stock-dependent demand on supply chain dynamicsLeeH. L.PadmanabhanV.WhangS.Information distortion in a supply chain: the bullwhip effectDejonckheereJ.DisneyS. M.LambrechtM. R.TowillD. R.Measuring and avoiding the bullwhip effect: a control theoretic approachDejonckheereJ.DisneyS. M.LambrechtM. R.TowillD. R.The impact of information enrichment on the bullwhip effect in supply chains: a control engineering perspectiveOuyangY.DaganzoC.Characterization of the bullwhip effect in linear, time-invariant supply chains: some formulae and testsStermanJ. D.Modeling managerial behavior: misperceptions of feedback in a dynamic decision making experimentHwarngH. B.XieN.Understanding supply chain dynamics: a chaos perspectiveWangH. W.LiuH. X.YangJ.-B.Dynamic analysis of a two-stage supply chain—a switched system theory approachGravesS. C.A single-item inventory model for a nonstationary demand processSunH.WuJ.Scale-free characteristics of supply chain distribution networksAlbertR.BarabásiA.-L.Statistical mechanics of complex networksBarabásiA.-L.AlbertR.JeongH.Mean-field theory for scale-free random networksStrogatzS. H.Exploring complex networksSuH.JiaG.ChenM. Z. Q.Semi-global containment control of multi-agent systems with intermittent input saturationSuH. S.ChenM. Z. Q.ChenG. R.Robust semi-global coordinated tracking of linear multi-agent systems with input saturationSuH.-S.ChenM. Z. Q.Multi-agent containment control with input saturation on switching topologiesLalwaniC. S.DisneyS. M.TowillD. R.Controllable, observable and stable state space representations of a generalized order-up-to policySipahiR.DeliceI. I.Stability of inventory dynamics in supply chains with three delaysDisneyS. M.TowillD. R.Van de VeldeW.Variance amplification and the golden ratio in production and inventory controlHobergK.BradleyJ. R.ThonemannU. W.Analyzing the effect of the inventory policy on order and inventory variability with linear control theoryHobergK.ThonemannU. W.BradleyJ. R.Analyzing the effect of inventory policies on the nonstationary performance with transfer functionsChenF.DreznerZ.RyanJ. K.Simchi-LeviD.Quantifying the bullwhip effect in a simple supply chain: the impact of forecasting, lead times, and informationOuyangY.The effect of information sharing on supply chain stability and the bullwhip effectBoukasE. K.Manufacturing systems: LMI approachNagataniT.HelbingD.Stability analysis and stabilization strategies for linear supply chainsCannellaS.CianciminoE.On the bullwhip avoidance phase: supply chain collaboration and order smoothingHelbingD.Modelling supply networks and business cycles as unstable transport phenomenaScholz-ReiterB.DashkovskiyS.GörgesM.NaujokL.Stability analysis of autonomously controlled production networksOuyangY.LiX.The bullwhip effect in supply chain networksDominguezR.CannellaS.FraminanJ. M.On bullwhip-limiting strategies in divergent supply chain networksCachonG. P.RandallT.SchmidtG. M.In search of the bullwhip effectSuckyE.The bullwhip effect in supply chains—an overestimated problem?KimJ. H.ParkH. B.H∞ state feedback control for generalized continuous/discrete time-delay system