Exploring the Impact of Network Structure and Demand Collaboration on the Dynamics of a Supply Chain Network Using a Robust Control Approach

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 collaborationsmight 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.


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][3][4][5], bullwhip effect [6][7][8][9], and even chaos [10][11][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 scalefree and small-world features [14].Note that the small-world feature and scaleless feature have been widely explored in complex network theory [15][16][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][19][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 2 Mathematical Problems in Engineering 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][24][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  ∞ 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.

Modeling the Supply Chain Network
Consider a two-layer supply chain network constituted with  distributors and  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: (i) 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.
(ii) 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.
(iii) 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.
(iv) Both the distributors and retailers have fixed replenishment lead times.
Denote symbols , 1 ≤  ≤ , and , 1 ≤  ≤ , as the th distributor and the th retailer, respectively.Then the supply relationship between retailers and distributors could be defined by a supply relationship matrix  = (  ) × ∈  × .
If   = 1, the distributor  and the retailer  establish the supply relationship between them.Otherwise, if   = 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  has only one element with value 1, then each retailer has only one immediate supplier.In contrast, if   = 1 (∀, ), 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   () be the demand received by retailer .Then the amount   1  2 × [  1 () −   2 ()] of the received demand of retailer  1 will be satisfied by retailer  2 if   1 () >   2 ().Note that the parameter   1  2 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  ≜ (  1  2 ) × ∈  × to exactly describe the demand diversion among retailers in the supply chain network.Assume that the demand of retailer  1 is   1 () and the demand of retailer  2 is   2 ().The essential demand of retailer  1 after demand collaboration becomes 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   (), 1 ≤  ≤ , denote the inventory position of the distributor , let   () denote the placed order by distributor  from external suppliers, and let   () denote the received order of distributor  from retailer  (1 ≤  ≤ ).Then, the resulting dynamics of the distributor  can be represented by 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   (), 1 ≤  ≤ , denote the inventory position of retailer ; then its balanced inventory equation could be obtained as 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  and the collaboration matrix .By doing this, we can incorporate the network structure and collaboration among retailers into the dynamic model.Let   = [ 1 ,  2 , . . .,   ]  , 1 ≤  ≤ , define the th column of the matrix .Then we could define the following two matrixes: where   = diag( 1 ,  2 , . . .,   ), 1 ≤  ≤ .We see that  1 ∈  × and  2 ∈  × .Based on the collaboration matrix , we could obtain the following matrix: The demand rate of retailer  is assumed to take the form where    is a known constant and   () 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   () =    ,  = 1, 2, . . ., , then the inventory positions and order quantities of retailers and distributors will return to constant or steady values after finite periods.Define    and    as the steady values of inventory position of distributor  and retailer , respectively.Let    and    denote the steady values of the amount of order of the distributor  placed from external suppliers and the amount of order of the retailer  placed from distributor , respectively.In such case, we can have lim 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: , where The input vector can be represented by Before defining the control vector, we will now introduce the operator Vec().
Definition 1. Assume that  ∈  × is a matrix with  rows and  columns and  = ( 1 ,  2 , . . .,   ), where   is the th column of .Then we can produce a single column vector using the following operator, which is given by For the supply chain network, the control vector ũ() can be represented as ũ() = [ũ  () ũ ()] ∈  + , where where Define ỹ() ∈  2++ as system output and let   denote the unit matrix with  dimensions.We then can develop the following unified state space model: where in which   and   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   = 1 and   = 0. Otherwise, the setting of   = 0 and   = 1 can be used to mitigate bullwhip effect or order fluctuations.Both the inventory and order dynamics will be considered if   = 1 and   = 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  ∞ control method.

The Robust Performance Metric and 𝐻 ∞ 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  ∞ 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: We see that the metric   is the amplification metric of the worst-case root mean square errors (RMSE) when the customer demand is unknown.Therefore, minimizing the value of   can effectively improve the dynamics and robustness of supply chain system.It is worthwhile to mention that this metric, also called  2 gain or  ∞ norm, has been widely used in the field of control engineering.Although it is not an economic metric, the metric   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   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   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   = 0 and   = 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   can take both order fluctuations and inventory fluctuations into account by setting   = 1 and   = 1.In this case, the exact formula of   for our supply chain network model becomes As emphasized, minimizing the value of the metric   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 ũ() =  x().Then the robust  ∞ control problem of our supply chain network model is to seek a feedback controller ũ() =  x() that satisfies the following optimization problem: The inventory rule to minimize the value of the robust metric   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  and a matrix  satisfies the following optimization problem: then system (11) is quadratically stable and the minimal value of the robust metric   can be approximated with the minimized  1/2 , with the state-feedback controller: ũ() =  −1 x().
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   , or the  ∞ 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.

Numerical Analysis
As introduced, the proposed robust  ∞ 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:

Mathematical Problems in Engineering
Var (r  ()) , (18) in which Var(⋅) is used to calculate the variances of time sequences.It is quite apparent to see that the metrics  1 and  2 , respectively, represent the order fluctuations at distributor layer and retailer layer.The latter two metrics  1 and  2 are used to measure the inventory oscillations for retailers and distributors.

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 Assume that the demand of each retailer follows the firstorder autoregressive process: in which  is the autoregressive parameter and   is the white noise.Setting different values for   and   , we can consider three different control ways: (i) control way 1: to reduce inventory fluctuation (  = 1,   = 0); (ii) control way 2: to reduce order fluctuation (  = 0,   = 1); (iii) control way 3: to reduce both order and inventory fluctuations (  = 1,   = 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   = 1,   = 1.As we consider a supply chain network composed of 3 retailers and 3 distributors, we could have  =  = 3.Then, the matrices  1 and  2 can be obtained as Because we do not consider demand collaboration in this section, the matrix   will be calculated as   = diag{1, 1, 1}.
From the value of  1 ,  2 , and   , we can further easily calculate the matrices ,   , , and .The inventory rule, a state-feedback controller, can be solved using Theorem 2 with LMI Toolbox in MATLAB, which is represented as 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,  1 ,  2 ,  1 , and  2 , which are summarized in Table 1.
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 tradeoff 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  ∞ control method is highly effective in improving dynamics and robustness.

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 , and then consider the change of the number of nodes.Assume that   = 1 and   = 1 to consider both inventory and order fluctuations.To concentrate our attention on network structure, we will not take demand collaboration into account.

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  can be varied to represent different structures.The supply relationship matrixes for the four types of structures are given by where  1 ,  2 ,  3 , and  4 are the supply relationship matrices for structure 1, structure 2, structure 3, and structure 4, respectively.
Using the optimization method introduced in Theorem 2, we could accordingly compute the minimal values of the robust metric   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.

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: (i) Case 1: increasing the number of distributors.(ii) 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.
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.

Impact of Demand Collaboration on Robust Performance.
As introduced, the collaboration mechanism among retailers is characterized by the collaboration matrix .The retailer  1 and retailer  2 form their collaboration relationship with the parameter   1  2 , 0 ≤   1  2 ≤ 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 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  are 1.We also assume that   = 1 and   = 1 in order to consider overall performance.In this study, we will discuss three different collaborating scenarios: (i) scenario 1:  12 =  34 =  56 = 0 (without cooperation); (ii) scenario 2:  12 =  34 =  56 = 0.25 (mild cooperation); (iii) scenario 3:  12 =  34 =  56 = 0.5 (close cooperation).
The minimal values of the robust metric   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,  1 ,  2 ,  1 , and  2 , 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.
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  ∞ control.

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  ∞ 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.

Figure 1 :
Figure 1: The structure of a general two-layer supply chain network.

Figure 2 :
Figure 2: Dynamic curves of inventory under three different control ways.
Order curves under the first control wayOrder dynamics for retailers Order dynamics for distributors Order curves under the third control way

Figure 3 :
Figure 3: Dynamic curves of order under three different control ways.

4 Figure 5 :
Figure 5: Order dynamics under different structures of supply chain network.

Figure 6 :
Figure 6: Impact of the number of nodes on robust performance.

Table 1 :
System performance under different control ways.

Table 2 :
System performance under different collaboration scenarios.