# System Dynamics Model for VMI&TPL Integrated Supply Chains

Authors: Guo Li; Xiaojing Wang; Zhaohua Wang
Journal: Discrete Dynamics in Nature and Society
 (2013)

---

## Abstract

This paper establishes VMI-APIOBPCS II model by extending VMI-APIOBPCS model from serial supply chain to distribution supply chain. Then TPL is introduced to this VMI distribution supply chain, and operational framework and process of VMI&TPL integrated supply chain are analyzed deeply. On this basis VMI-APIOBPCS II model is then changed to VMI&TPL-APIOBPCS model and VMI&TPL integrated operation mode is simulated. Finally, compared with VMI-APIOBPCS model, the TPL’s important role of goods consolidation and risk sharing in VMI&TPL integrated supply chain is analyzed in detail from the aspects of bullwhip effect, inventory level, service level, and so on.


---

## Body

## 1. Introduction

Under vendor-managed inventory (VMI) operation mode, many suppliers outsource their logistics to third-party logistics (TPL) due to their poor logistics capabilities. So far, that TPL participates in VMI has been widely used in many industries. For example, Dell and Lenovo both chose Burlington Company to help them operate VMI service, and Wuhan Shenlong Automobile Company in China allows GEFCO to provide VMI service with components supply. On the one hand, this integrated operational model combining VMI with TPL ensures that the supply chain information is shared fully on one central platform. On the other hand, it can take full advantage of TPL and reduce the total operational cost of the supply chain.As to the research of VMI&TPL integrated replenishment and delivery, Çetinkaya et al. [1] take Dell as an example, which outsources VMI business to Burlington Logistics, and analyze TPL replenishment and delivery strategies. They do not only consider the optimal delivery strategies about logistics outsourcing but also find out differences of the optimal delivery strategies before and after outsourcing. Based on the above study, Çetinkaya and Lee [2] consider the time-based delivery policy and obtain the optimal delivery time structure with transportation lot constraints and capability limitations while the demand of retailers obeys the Poisson distribution. Lee et al. [3] assume that the replenishment and delivery can be started at the beginning of each period with determined demand and finite horizon and that the lead time for delivery of the replenishment is zero. They consider the problem of inventory and transportation integration, which is similar to Çetinkaya and Lee [2]. In order to achieve economies of transport scale, TPL implements goods consolidation strategy. As a result, products may be delivered to retailers in an earlier or later time, which would lead to the inventory cost or shortage cost. Their research shows that the problem is NP hard and that the demand in each period must be satisfied in just one delivery, with transportation constraints. Furthermore, they point out that if each delivery can meet demands in several consecutive cycles and the demand in the first and last periods may be met by two deliveries then the optimal replenishment and delivery policy exists. They even propose a polynomial algorithm to solve the above problem of optimal replenishment and delivery. Çetinkaya et al. [4] provide a foundation for a comparison of the impact of time-based versus quantity-based consolidation in the context of integrated inventory and transportation decisions. Numerical and analytical results verify that quantity-based consolidation is superior to the time-based version in terms of the resulting cost. Furthermore, several easily implementable TQ-based policies are proposed and their impacts on cost and service are compared to those of time-based and quantity-based versions via simulation. Mutlu et al. [5] extend the results in Çetinkaya et al. [4] and develop an analytical model for computing the expected long-run average cost of a consolidation system implementing a TQ-based policy. The presented analytical results prove that (i) the optimal TQ-based policy outperforms the optimal time-based policy and (ii) the optimal quantity-based policy is superior to the other two (i.e., optimal time-based and TQ-based) policies in terms of cost. Considering the expected maximum waiting time as a measure of timely delivery performance, however, they numerically demonstrate that the TQ-based policies improve on the quantity-based policies significantly with only a slight increase in the cost.Besides, some scholars consider replenishment strategies of TPL in VMI mode under different conditions, such as Çetinkaya et al. [6], Dejonckheere et al. [7],  Wikrom et al. [8–11], Lee [12], Hwang [13], and Howard and Marklund [14]. Çetinkaya et al. [6] consider different delivery and replenishment strategies under two kinds of transportation modes. One is self-transportation which is the same as in the above literature. The other one is outsourcing transportation. As transportation providers will take the discount policy to encourage suppliers to transport more, transportation cost may be a piecewise function in this situation. They propose two kinds of delivery and replenishment strategies based on time and quantity under different transportation modes. Dejonckheere et al. [7] investigate the utilization of a linear (Type II) or quadratic (Type III) instead of a constant (Type I) exponential smoothing forecasting mechanism in the continuous-time APIOBPCS model. Mustafa et al. [15] consider the impact of coordinated replenishment and shipment in inventory/distribution systems and analyze a system with multiple retailers and one outside supplier. They present a centralized ordering policy that orders for all retailers and some other well-known policies like (a) can-order policy, (b) echelon inventory policy, and (c) fixed-replenishment interval policy. Leopoldo et al. [16] present an alternative heuristic algorithm to solve the vendor-managed inventory system with multiproduct and multiconstraint based on EOQ with backorders considering two classical backorders costs: linear and fixed. Sadeghia et al. [17] develop a constrained multivendor multiretailer single-warehouse (MV-MR-SW) supply chain, in which both the space and the annual number of orders of the central warehouse are limited. Since the problem is formulated into an integer nonlinear programming model, the metaheuristic algorithm of particle swarm optimization (PSO) is presented to find an approximate optimum solution of the problem. In the proposed PSO algorithm, a genetic algorithm (GA) with an improved operator, namely, the boundary operator, is employed as a local searcher to turn it to a hybrid PSO. Moreover, Harigaa et al. [18] consider a supply chain composed of a single vendor and multiple retailers operating under a VMI contract that specifies limits on retailers' stock levels. They address the problem of synchronizing the vendor's cycle time with the buyers' unequal ordering cycles by developing a mixed integer nonlinear program that minimizes the joint relevant inventory costs under storage restrictions.As to logistics optimization based on system dynamics, Towill [19] establishes a new inventory-and-order-based production and control system (IOBPCS) by extending production inventory control (PIC) [20] and then optimizes the system by using the coefficient plane model. Sterman [21] constructs a general inventory management model, making use of the system dynamics, and points out that different complexity of feedback in supply chain inventory management system and the pressure of time usually lead decision makers to misunderstand the feedback information and thus make irrational decisions. John et al. [22] introduce WIP feedback control mechanism into IOBPCS model and expand the IOBPCS model into APIOBPCS. Later, Mason-Jones et al. [23] analyze the function of WIP feedback control mechanism in the models by comparing IOBPCS and APIOBPCS. Disney and Towill [24, 25] construct VMI-APIOBPCS model and analyze VMI strategy’s effects on the supply chain bullwhip effect, customer service level and inventory costs with the assumption that the enterprises face the obvious fluctuations of demand. Besides, they optimize the VMI-APIOBPCS model and obtain the optimal parameters after considering different weights of production adjustment costs, different proportions of inventory costs, and different coefficients of safety stock. Disney and Towill [24, 25] study a simple vendor-managed inventory (VMI) supply chain consisting of one production unit and one distributor. In VMI systems all supply points in the chain have access to stock positions for setting production and distribution targets. The discrete-time APIOBPCS model is used to describe the dynamics of the manufacturing unit. Pure delay is initially utilized to model the production delay. The only difference to the APIOBPCS structure presented previously is that instead of the demand signal CONS the manufacturing facility receives a “virtual” consumption signal. This is caused by adding in each time period the demand signal received by the distributor to the difference between the current time period and the previous period reorder-point. The system is checked for stability. The stability criteria that are produced are also valid for the standard APIOBPCS model, since the distributor’s policy described previously is a stable feed-forward element. One year later, Disney and Towill [26–28] analyze deeply how VMI affects the bullwhip effect in the supply chain and compare the VMI supply chain’s expected performance with that of a traditional supply chain. VMI strategy shows that it has better reactions when demand is not steady, and this kind of instability may be caused by discounts available for orders or price’s changes. Besides, the restoration of inventory level will be improved dramatically by VMI strategy. Moreover, Disney and Towill [26–28] concentrate on VMI strategy’s effects on transport operations in supply chain, especially the batch problem in transportation strategy. By using system dynamics, they establish three different kinds of models—the traditional one, the internal integrated one, and the VMI one. The simulation case shows that VMI model can reduce transportation frequency by adopting a larger batch without influencing the dynamic performance of the entire supply chain. Wikner [29] presents a methodology that introduces structure dependencies of MLMS systems in the IOBPCS production control framework. The methodology uses matrix representation to account for multiple informational channels. It is shown that for a single-level single-stage system the model is reduced to the standard IOBPCS format. The extended model has the capability to describe the dynamics of both pull-driven (base stock, kanban) and push-driven (MRP) policies.In the other field, through STELLA/iThink software platform, Chen et al. [30] construct a system dynamics model of inventory management, analyze system structures and operational mechanisms of VMI inventory management and traditional inventory management, and finally compare their operational performance. Yang and Liu [31] extend VMI-APIOBPCS model from one supplier-one retailer supply chain to one supplier-two retailers supply chain and then construct VMI-APIOBPCS II. With the integration of TPL, they establish VMI&TPL-APIOBPCS model and the simulation shows that TPL can help reduce bullwhip effect in the supply chain availably. Cho and Lazaro [32] extend PID controller for just-in-time production scheduling. Lin et al. [33] develop a fuzzy system dynamic to simulate vendor-managed inventory, automatic pipeline, and inventory-and-order-based production control system (VMI-APIOBPCS) model based on fuzzy difference equations, and these operators of difference equations adopt the weakest t-norm (TW) operators. The results of fuzzy VMI-APIOBPCS model can provide the whole extended information regarding the system behavior uncertainties for the decision makers with fuzzy interval.After that Darya and Martin [34] address the steady-state optimization of a supply chain model that belonged to the class of vendor-managed inventory, automatic pipeline, and inventory-and-order based production control systems (VMI-APIOBPCS). They optimize the supply chain with the so-called normal vector method, which has specifically been developed for the economic optimization of uncertain dynamical systems with constraints on dynamics. Kristianto et al. [35] propose an adaptive fuzzy control application to support a vendor-managed inventory (VMI). The methodology applies fuzzy control to generate an adaptive smoothing constant in the forecast method, production, and delivery plan to eliminate, for example, the rationing and gaming or the Houlihan effect and the order batching effect or the Burbidge effects and finally the bullwhip effect. In order to improve the level of integration in all aspects of supply chain reconfiguration, Kristianto et al. [36] construct an optimum supply chain network by combining optimization at the strategic and tactical level. A system dynamic based computer simulation model is used to validate the operations of the supply chain. The performance of the system is measured in terms of backorders and inventory level. The results and analysis indicate that fewer stockholding points and a shorter review period of demand can improve performance in this respect.Our work differs from these studies in important aspects. First of all, these models are mainly constructed based on two echelon supply chains, that is, Disney and Towill [24–28], Lin et al. [33], Darya and Martin [34], and Kristianto et al. [35, 36]. In our model, however, VMI&TPL integrated operational model is a relatively complex three-echelon supply chain system. The TPL role integrated with VMI should be explored. Furthermore, we also model VMI&TPL integrated supply chain from aspects of bullwhip effect, inventory levels, service level with stochastic demand, and other uncertainties, compared with VMI-APIOBPCS model. Since system dynamics suit researches of complex systems very well, VMI&TPL-APIOBPCS model is constructed based on system dynamics to simulate and analyze the performance of VMI&TPL integrated operational model.


## 2. IOBPCS Model Family

### 2.1. Definitions of Parameters and Variables

Table1 is the definitions of parameters and variables used in this paper.Table 1
Definitions of parameters and variables.
T
i: inventory adjustment time
T
a: demand smooth time
T
p: production delay time
T
p
^: estimated value of production delay time
CON: demand rate
ACON: demand rate after forecast
ORT: production rate
CRT: production fulfillment rate
TINV: target inventory level
INV: actual inventory level
EINV: inventory deviation
WIP: work-in-process
TWIP: target work-in-process
EWIP: work-in-process deviation


### 2.2. IOBPCS Model

As shown in Figure1, inventory levels of work-in-process and finished products can be controlled by production order rate in IOBPCS system, and customers’ demand is met by inventory.Figure 1
IOBPCS model of production and inventory system.Figure2 is the causal relationship in IOBPCS model, including 4 main parts, which are demand forecast feed-forward loop, production delays, inventory feedback loop, and target inventory. Production delays refer to the time from production orders to production fulfillment rate. In series IOBPCS models, it is assumed that production process in production lead time meets a certain order, and orders keep a sequence of events. Demand forecast feed-forward loop refers to the demand forecasting mechanism which is used to predict the demand in and after production lead time. Inventory feedback loop actually is a kind of inventorydeviation adjustment mechanism. It is necessary to produce more goods to adjust inventory deviation when the actual inventory level and target inventory level differ greatly.Figure 2
IOBPCS model causal relationship.The productivity in IOBPCS model is decided by the demand forecasting mechanism and the inventory deviation adjustment mechanism, while the inventory deviation adjustment mechanism is decided by the inventory adjustment time and the target inventory level, and the demand forecasting mechanism is decided by the demand smooth time. Therefore, IOBPCS system optimization includes the definition of two basal parameters, such as demand smooth time and inventory deviation adjustment time. When designing the best production control strategy, cost from two sides should be balanced, including production cost due to production fluctuation and the inventory cost (or shortage cost) as the inventory level changes.


### 2.3. IOBPCS Expansion Model

After many scholars amending and improving the model based on IOBPCS model, now it has been turned into an IOBPCS model family consisting of five parts as Figure3 shows.Figure 3
The main parts of expanded IOBPCS model.


#### 2.3.1. Production Delay

The lead time of production delays can be regarded as production rhythm smooth time, which describes the speed of adjusting production rhythm to production order changes (ORATE). The production delay is one of the system characteristics, which cannot be controlled at random by the system’s designers, but designing different delay models will have an important effect on the entire system’s performance. Formula (1) is the dynamic behavior of the three delay models:
(1)Gp(s)=1((Tp/n)s+1)n,wheren=1, first-order delay; n=3, third-order delay; n=8, pure delay; Tp is production delay time.


#### 2.3.2.   Target Inventory Level (DINV)

In IOBPCS model, target inventory level is a fixed value or the integer multiple of the demand forecasting number (ACON) after smoothing. Target inventory level is variable in VIOBPCS. Compared with IOPBPCS, the width of ORATE is larger, but the inventory adjustment response time is shorter. The only difference between APIOBPCS model and APVIOBPCS model lies in the setting of target inventory. In APIOBPCS model, target inventory levelDINV=k*ACON, where k is a positive integer.


#### 2.3.3. Demand Forecasting Mechanism

Demand forecasting mechanism is an important part of the feed-forward loop. Demand forecasting mechanism is measured by exponential smoothing method in most literatures, because exponential smoothing method comprehensively includes all the historical information and is easy to use and to formulate a model. Exponential forecasting method (e.g., single exponential smoothing method, double exponential smoothing method, and triple exponential smoothing method) makes the steady-state error of system inventory in phase step and slope demand zero, but the steady-state error of system inventory becomes larger and larger when the demand function is a parabola.Single exponential smoothing transfer function ins region is
(2)Ga(s)=1Tas+1.Double exponential smoothing transfer function ins region is
(3)Ga(s)=2Tas+1Ta2s2+2Tas+1.Triple times exponential smoothing transfer function ins region is
(4)Ga(s)=3Ta2s2+3Tas+1Ta3s3+3Ta2s2+3Tas+1.


#### 2.3.4. Inventory Deviation Adjustment Mechanism

Inventory deviation adjustment mechanism is an inventory feedback loop which controls inventory deviation by controlling productivity. Inventory adjustment mechanism needs to consider production delay effect which means that only after a regular period of time can the controlling decision about productivity adjust the inventory level. The purpose of inventory adjustment is to reach target inventory level in a period of time (Ti). When the adjustment time is shorter, the inventory can adjust faster. As a result, shortage cost and risks become smaller. However, it requires manufacturers for a higher production capacity, because the system needs to correct the early inventory deviation in a short time by adjusting productivity, and this can lead to higher cost of production and work-in-process inventory.


#### 2.3.5. WIP Inventory Deviation Adjustment Mechanism

WIP inventory deviation is caused by actual work-in-process inventory in contradiction with target work-in-process inventory when demand changes. WIP adjustment mechanism adjusts inventory deviation by controlling productivity so that it can reach the target value in a period of time (Tw). Therefore, WIP inventory deviation adjustment rate is one of the three parts of the productivity control mechanisms. It is necessary to study and analyze the production process and get a relatively exact estimate value by making statistical analysis of production delay time before designing the productivity control mechanism. If the observed production delay time is different from actual time, inventory level in steady-state will not be in accord with target inventory level, which will cause more risk of inventory or shortage.


### 2.4. APIOBPCS Model

Figure4 is a block diagram of APIOBPCS model lying in s region. The purpose of system designing is to find a proper target inventory level and design three optimal control mechanisms (demand forecasting mechanism, inventory deviation adjustment, and WIP inventory deviation adjustment mechanism), in order to minimize system cost, including production cost and inventory cost.Figure 4
Block diagram of APIOBPCS model.In general, inventory dynamic fluctuation is measured by inventory rising time and adjustment time and overshoot, and productivity dynamic change is analyzed by frequency response method.Here is the main control mechanism of APIOBPCS model in phase-step demand.(1)
Forecasting mechanism. Formula (5) is the transfer function in s region with single exponential smoothing:
(5)Ga(s)=1Tas+1.(2)
Target inventory setting. TINV = 0.(3)
Production process. Formula (6) is the transfer function in s region with first-order delay:
(6)Gp(s)=1Tps+1.As a result, two important transfer functions about productivity and changes in inventory level can be obtained as follows:(7)F1(s)=ORT(s)CON(s)=TpTis+Tw+(Ta+Ti)Tws(1+Tas)[TiTwTps2+Tis(Tp+Tw)+Tw],F2(s)=INV(s)CON(s)=Ti(Tp^-Tp)-TiTpTws-TiTpTas-TiTwTas-TaTiTwTps2(1+Tas)[TiTwTps2+Ti(Tp+Tw)s+Tw].


## 2.1. Definitions of Parameters and Variables

Table1 is the definitions of parameters and variables used in this paper.Table 1
Definitions of parameters and variables.
T
i: inventory adjustment time
T
a: demand smooth time
T
p: production delay time
T
p
^: estimated value of production delay time
CON: demand rate
ACON: demand rate after forecast
ORT: production rate
CRT: production fulfillment rate
TINV: target inventory level
INV: actual inventory level
EINV: inventory deviation
WIP: work-in-process
TWIP: target work-in-process
EWIP: work-in-process deviation


## 2.2. IOBPCS Model

As shown in Figure1, inventory levels of work-in-process and finished products can be controlled by production order rate in IOBPCS system, and customers’ demand is met by inventory.Figure 1
IOBPCS model of production and inventory system.Figure2 is the causal relationship in IOBPCS model, including 4 main parts, which are demand forecast feed-forward loop, production delays, inventory feedback loop, and target inventory. Production delays refer to the time from production orders to production fulfillment rate. In series IOBPCS models, it is assumed that production process in production lead time meets a certain order, and orders keep a sequence of events. Demand forecast feed-forward loop refers to the demand forecasting mechanism which is used to predict the demand in and after production lead time. Inventory feedback loop actually is a kind of inventorydeviation adjustment mechanism. It is necessary to produce more goods to adjust inventory deviation when the actual inventory level and target inventory level differ greatly.Figure 2
IOBPCS model causal relationship.The productivity in IOBPCS model is decided by the demand forecasting mechanism and the inventory deviation adjustment mechanism, while the inventory deviation adjustment mechanism is decided by the inventory adjustment time and the target inventory level, and the demand forecasting mechanism is decided by the demand smooth time. Therefore, IOBPCS system optimization includes the definition of two basal parameters, such as demand smooth time and inventory deviation adjustment time. When designing the best production control strategy, cost from two sides should be balanced, including production cost due to production fluctuation and the inventory cost (or shortage cost) as the inventory level changes.


## 2.3. IOBPCS Expansion Model

After many scholars amending and improving the model based on IOBPCS model, now it has been turned into an IOBPCS model family consisting of five parts as Figure3 shows.Figure 3
The main parts of expanded IOBPCS model.


### 2.3.1. Production Delay

The lead time of production delays can be regarded as production rhythm smooth time, which describes the speed of adjusting production rhythm to production order changes (ORATE). The production delay is one of the system characteristics, which cannot be controlled at random by the system’s designers, but designing different delay models will have an important effect on the entire system’s performance. Formula (1) is the dynamic behavior of the three delay models:
(1)Gp(s)=1((Tp/n)s+1)n,wheren=1, first-order delay; n=3, third-order delay; n=8, pure delay; Tp is production delay time.


### 2.3.2.   Target Inventory Level (DINV)

In IOBPCS model, target inventory level is a fixed value or the integer multiple of the demand forecasting number (ACON) after smoothing. Target inventory level is variable in VIOBPCS. Compared with IOPBPCS, the width of ORATE is larger, but the inventory adjustment response time is shorter. The only difference between APIOBPCS model and APVIOBPCS model lies in the setting of target inventory. In APIOBPCS model, target inventory levelDINV=k*ACON, where k is a positive integer.


### 2.3.3. Demand Forecasting Mechanism

Demand forecasting mechanism is an important part of the feed-forward loop. Demand forecasting mechanism is measured by exponential smoothing method in most literatures, because exponential smoothing method comprehensively includes all the historical information and is easy to use and to formulate a model. Exponential forecasting method (e.g., single exponential smoothing method, double exponential smoothing method, and triple exponential smoothing method) makes the steady-state error of system inventory in phase step and slope demand zero, but the steady-state error of system inventory becomes larger and larger when the demand function is a parabola.Single exponential smoothing transfer function ins region is
(2)Ga(s)=1Tas+1.Double exponential smoothing transfer function ins region is
(3)Ga(s)=2Tas+1Ta2s2+2Tas+1.Triple times exponential smoothing transfer function ins region is
(4)Ga(s)=3Ta2s2+3Tas+1Ta3s3+3Ta2s2+3Tas+1.


### 2.3.4. Inventory Deviation Adjustment Mechanism

Inventory deviation adjustment mechanism is an inventory feedback loop which controls inventory deviation by controlling productivity. Inventory adjustment mechanism needs to consider production delay effect which means that only after a regular period of time can the controlling decision about productivity adjust the inventory level. The purpose of inventory adjustment is to reach target inventory level in a period of time (Ti). When the adjustment time is shorter, the inventory can adjust faster. As a result, shortage cost and risks become smaller. However, it requires manufacturers for a higher production capacity, because the system needs to correct the early inventory deviation in a short time by adjusting productivity, and this can lead to higher cost of production and work-in-process inventory.


### 2.3.5. WIP Inventory Deviation Adjustment Mechanism

WIP inventory deviation is caused by actual work-in-process inventory in contradiction with target work-in-process inventory when demand changes. WIP adjustment mechanism adjusts inventory deviation by controlling productivity so that it can reach the target value in a period of time (Tw). Therefore, WIP inventory deviation adjustment rate is one of the three parts of the productivity control mechanisms. It is necessary to study and analyze the production process and get a relatively exact estimate value by making statistical analysis of production delay time before designing the productivity control mechanism. If the observed production delay time is different from actual time, inventory level in steady-state will not be in accord with target inventory level, which will cause more risk of inventory or shortage.


## 2.3.1. Production Delay

The lead time of production delays can be regarded as production rhythm smooth time, which describes the speed of adjusting production rhythm to production order changes (ORATE). The production delay is one of the system characteristics, which cannot be controlled at random by the system’s designers, but designing different delay models will have an important effect on the entire system’s performance. Formula (1) is the dynamic behavior of the three delay models:
(1)Gp(s)=1((Tp/n)s+1)n,wheren=1, first-order delay; n=3, third-order delay; n=8, pure delay; Tp is production delay time.


## 2.3.2.   Target Inventory Level (DINV)

In IOBPCS model, target inventory level is a fixed value or the integer multiple of the demand forecasting number (ACON) after smoothing. Target inventory level is variable in VIOBPCS. Compared with IOPBPCS, the width of ORATE is larger, but the inventory adjustment response time is shorter. The only difference between APIOBPCS model and APVIOBPCS model lies in the setting of target inventory. In APIOBPCS model, target inventory levelDINV=k*ACON, where k is a positive integer.


## 2.3.3. Demand Forecasting Mechanism

Demand forecasting mechanism is an important part of the feed-forward loop. Demand forecasting mechanism is measured by exponential smoothing method in most literatures, because exponential smoothing method comprehensively includes all the historical information and is easy to use and to formulate a model. Exponential forecasting method (e.g., single exponential smoothing method, double exponential smoothing method, and triple exponential smoothing method) makes the steady-state error of system inventory in phase step and slope demand zero, but the steady-state error of system inventory becomes larger and larger when the demand function is a parabola.Single exponential smoothing transfer function ins region is
(2)Ga(s)=1Tas+1.Double exponential smoothing transfer function ins region is
(3)Ga(s)=2Tas+1Ta2s2+2Tas+1.Triple times exponential smoothing transfer function ins region is
(4)Ga(s)=3Ta2s2+3Tas+1Ta3s3+3Ta2s2+3Tas+1.


## 2.3.4. Inventory Deviation Adjustment Mechanism

Inventory deviation adjustment mechanism is an inventory feedback loop which controls inventory deviation by controlling productivity. Inventory adjustment mechanism needs to consider production delay effect which means that only after a regular period of time can the controlling decision about productivity adjust the inventory level. The purpose of inventory adjustment is to reach target inventory level in a period of time (Ti). When the adjustment time is shorter, the inventory can adjust faster. As a result, shortage cost and risks become smaller. However, it requires manufacturers for a higher production capacity, because the system needs to correct the early inventory deviation in a short time by adjusting productivity, and this can lead to higher cost of production and work-in-process inventory.


## 2.3.5. WIP Inventory Deviation Adjustment Mechanism

WIP inventory deviation is caused by actual work-in-process inventory in contradiction with target work-in-process inventory when demand changes. WIP adjustment mechanism adjusts inventory deviation by controlling productivity so that it can reach the target value in a period of time (Tw). Therefore, WIP inventory deviation adjustment rate is one of the three parts of the productivity control mechanisms. It is necessary to study and analyze the production process and get a relatively exact estimate value by making statistical analysis of production delay time before designing the productivity control mechanism. If the observed production delay time is different from actual time, inventory level in steady-state will not be in accord with target inventory level, which will cause more risk of inventory or shortage.


## 2.4. APIOBPCS Model

Figure4 is a block diagram of APIOBPCS model lying in s region. The purpose of system designing is to find a proper target inventory level and design three optimal control mechanisms (demand forecasting mechanism, inventory deviation adjustment, and WIP inventory deviation adjustment mechanism), in order to minimize system cost, including production cost and inventory cost.Figure 4
Block diagram of APIOBPCS model.In general, inventory dynamic fluctuation is measured by inventory rising time and adjustment time and overshoot, and productivity dynamic change is analyzed by frequency response method.Here is the main control mechanism of APIOBPCS model in phase-step demand.(1)
Forecasting mechanism. Formula (5) is the transfer function in s region with single exponential smoothing:
(5)Ga(s)=1Tas+1.(2)
Target inventory setting. TINV = 0.(3)
Production process. Formula (6) is the transfer function in s region with first-order delay:
(6)Gp(s)=1Tps+1.As a result, two important transfer functions about productivity and changes in inventory level can be obtained as follows:(7)F1(s)=ORT(s)CON(s)=TpTis+Tw+(Ta+Ti)Tws(1+Tas)[TiTwTps2+Tis(Tp+Tw)+Tw],F2(s)=INV(s)CON(s)=Ti(Tp^-Tp)-TiTpTws-TiTpTas-TiTwTas-TaTiTwTps2(1+Tas)[TiTwTps2+Ti(Tp+Tw)s+Tw].


## 3. VMI-APIOBPCS Model

### 3.1. Definitions of Parameters and Variables

Table2 is the definitions of parameters and variables used in this paper.Table 2
Definitions of parameters and variables.
Parameters and variables of suppliers
VCON: demand rate of suppliers
AVCON: demand factor of suppliers after smoothing
VINV: inventory level of suppliers
TINV: system target inventory level
WIP: work-in-process inventory level
TWIP: target work-in-process inventory level
EINV: inventory deviation
EWIP: deviation of work-in-process inventory
ORT: productivity of suppliers
CRT: production fulfillment rate of suppliers
T
i: adjustment time of system inventory deviation
T
a: demand rate of suppliers smooth time
T
w: work-in-process inventory deviation adjustment time
T
p: production delay time of suppliers
AEWIP: work-in-process inventory deviation regulation factor
AEINV: deviation adjustment rate of system inventory
Parameters and variables of retailers
CON: demand rate of retailers
ACON: demand rate of retailers after forecasting
RINP: inventory level of retailers (including inventory on the way)
SS: safety inventory level
G: safety inventory factor of retailers
ROP: reorder-point of distributors
SRT: delivering rate to retailers
DSS: reorder-point of retailers variation
GIT: inventory of distributors on the way
RINV: actual inventory level of retailers
L: transportation time from suppliers to retailers
ETQ: economic order quantityIn VMI operation mode, retailers share inventory information and sales information with suppliers dynamically and determine the customer service level together with suppliers. According to the fixed customer service level, suppliers choose quantity-based delivering model, which means that vehicle shipped way is chosen in order to guarantee economical efficiency of transportation when the total inventory level is lower than the reorder-point. Figure5 is VMI-APIOBPCS system dynamics model.Figure 5
VMI-APIOBPCS system dynamics model.Here are the relational formulas in VMI-APIOBPCS model. They are as follows.


#### 3.1.1. Suppliers Production or Replenishment Mechanism

(1)
Demand forecasting mechanism:(8)AVCONt=AVCONt-1+VCONt-AVCONt-11+Ta.
(2)
Target work-in-process inventory level:(9)TWIPt=AVCONt×TP.
(3)
Work-in-process inventory level:(10)WIPt=WIPt-1+ORTt-CRTt.
(4)
Finished goods inventory level:(11)VINVt=VINVt-1+CRTt-SRTt.
(5)
Target inventory level:(12)TINVt=AVCONt×Ts.
(6)
System inventory level:(13)SINVt=VINVt+RINVt.
(7)
Finished goods fulfillment rate:(14)CRTt=delay{ORTt,Tp}.
(8)
Productivity:(15)ORTt=AVCONt+AEINVt+AEWIPt.
(9)
Inventory deviation adjustment rate:(16)AEINVt=TINVt-SINVtTi.
(10)
Work-in-process inventory deviation adjustment rate:(17)AEWIPt=TWIPt-WIPtTw.


#### 3.1.2. Suppliers Delivering Mechanism

(1)
Reorder-point of retailers:(18)ROPt=ROPt-1+SSt-ROPt-11+Tq.
(2)
Inventory levels of retailers:(19)RINPt=RINVt+GITt.
(3)
Order arrival rate:(20)ARTt=delay{SRTt,L}.
(4)
Safety inventory setting of retailers:(21)SSt=CONt×G.


#### 3.1.3. Information Sharing Mechanism in VMI Supply Chain

(1)
Suppliers can get customers’ demand in time and obtain the actual total customer demand downstream through terminal customer information, including the terminal customer demand and changes of reorder-points of downstream retailers:(22)VCONt=CONt+dSSt=CONt+ROTt-ROTt-1.
(2)
Suppliers can check the inventory level of downstream retailers so that they can get the total supply chain inventory level which can optimize the total supply chain inventory decision:(23)SINVt=RINPt+VINVt.


### 3.2. VMI Distribution Supply Chain

In the actual operational process, suppliers can adopt VMI model on multiple downstream retailers. According to the aforementioned VMI serial supply chain when suppliers provide multiple retailers downstream with VMI service, they can get VMI distribution supply chain, as Figure6 shows, including one supplier and two retailers.Figure 6
Operation mode of VMI distribution supply chain.VMI-APIOBPCS II model as shown in Figure7 is constructed based on the operation mode in Figure 6.Figure 7
VMI-APIOBPCS II system dynamics model.


## 3.1. Definitions of Parameters and Variables

Table2 is the definitions of parameters and variables used in this paper.Table 2
Definitions of parameters and variables.
Parameters and variables of suppliers
VCON: demand rate of suppliers
AVCON: demand factor of suppliers after smoothing
VINV: inventory level of suppliers
TINV: system target inventory level
WIP: work-in-process inventory level
TWIP: target work-in-process inventory level
EINV: inventory deviation
EWIP: deviation of work-in-process inventory
ORT: productivity of suppliers
CRT: production fulfillment rate of suppliers
T
i: adjustment time of system inventory deviation
T
a: demand rate of suppliers smooth time
T
w: work-in-process inventory deviation adjustment time
T
p: production delay time of suppliers
AEWIP: work-in-process inventory deviation regulation factor
AEINV: deviation adjustment rate of system inventory
Parameters and variables of retailers
CON: demand rate of retailers
ACON: demand rate of retailers after forecasting
RINP: inventory level of retailers (including inventory on the way)
SS: safety inventory level
G: safety inventory factor of retailers
ROP: reorder-point of distributors
SRT: delivering rate to retailers
DSS: reorder-point of retailers variation
GIT: inventory of distributors on the way
RINV: actual inventory level of retailers
L: transportation time from suppliers to retailers
ETQ: economic order quantityIn VMI operation mode, retailers share inventory information and sales information with suppliers dynamically and determine the customer service level together with suppliers. According to the fixed customer service level, suppliers choose quantity-based delivering model, which means that vehicle shipped way is chosen in order to guarantee economical efficiency of transportation when the total inventory level is lower than the reorder-point. Figure5 is VMI-APIOBPCS system dynamics model.Figure 5
VMI-APIOBPCS system dynamics model.Here are the relational formulas in VMI-APIOBPCS model. They are as follows.


### 3.1.1. Suppliers Production or Replenishment Mechanism

(1)
Demand forecasting mechanism:(8)AVCONt=AVCONt-1+VCONt-AVCONt-11+Ta.
(2)
Target work-in-process inventory level:(9)TWIPt=AVCONt×TP.
(3)
Work-in-process inventory level:(10)WIPt=WIPt-1+ORTt-CRTt.
(4)
Finished goods inventory level:(11)VINVt=VINVt-1+CRTt-SRTt.
(5)
Target inventory level:(12)TINVt=AVCONt×Ts.
(6)
System inventory level:(13)SINVt=VINVt+RINVt.
(7)
Finished goods fulfillment rate:(14)CRTt=delay{ORTt,Tp}.
(8)
Productivity:(15)ORTt=AVCONt+AEINVt+AEWIPt.
(9)
Inventory deviation adjustment rate:(16)AEINVt=TINVt-SINVtTi.
(10)
Work-in-process inventory deviation adjustment rate:(17)AEWIPt=TWIPt-WIPtTw.


### 3.1.2. Suppliers Delivering Mechanism

(1)
Reorder-point of retailers:(18)ROPt=ROPt-1+SSt-ROPt-11+Tq.
(2)
Inventory levels of retailers:(19)RINPt=RINVt+GITt.
(3)
Order arrival rate:(20)ARTt=delay{SRTt,L}.
(4)
Safety inventory setting of retailers:(21)SSt=CONt×G.


### 3.1.3. Information Sharing Mechanism in VMI Supply Chain

(1)
Suppliers can get customers’ demand in time and obtain the actual total customer demand downstream through terminal customer information, including the terminal customer demand and changes of reorder-points of downstream retailers:(22)VCONt=CONt+dSSt=CONt+ROTt-ROTt-1.
(2)
Suppliers can check the inventory level of downstream retailers so that they can get the total supply chain inventory level which can optimize the total supply chain inventory decision:(23)SINVt=RINPt+VINVt.


## 3.1.1. Suppliers Production or Replenishment Mechanism

(1)
Demand forecasting mechanism:(8)AVCONt=AVCONt-1+VCONt-AVCONt-11+Ta.
(2)
Target work-in-process inventory level:(9)TWIPt=AVCONt×TP.
(3)
Work-in-process inventory level:(10)WIPt=WIPt-1+ORTt-CRTt.
(4)
Finished goods inventory level:(11)VINVt=VINVt-1+CRTt-SRTt.
(5)
Target inventory level:(12)TINVt=AVCONt×Ts.
(6)
System inventory level:(13)SINVt=VINVt+RINVt.
(7)
Finished goods fulfillment rate:(14)CRTt=delay{ORTt,Tp}.
(8)
Productivity:(15)ORTt=AVCONt+AEINVt+AEWIPt.
(9)
Inventory deviation adjustment rate:(16)AEINVt=TINVt-SINVtTi.
(10)
Work-in-process inventory deviation adjustment rate:(17)AEWIPt=TWIPt-WIPtTw.


## 3.1.2. Suppliers Delivering Mechanism

(1)
Reorder-point of retailers:(18)ROPt=ROPt-1+SSt-ROPt-11+Tq.
(2)
Inventory levels of retailers:(19)RINPt=RINVt+GITt.
(3)
Order arrival rate:(20)ARTt=delay{SRTt,L}.
(4)
Safety inventory setting of retailers:(21)SSt=CONt×G.


## 3.1.3. Information Sharing Mechanism in VMI Supply Chain

(1)
Suppliers can get customers’ demand in time and obtain the actual total customer demand downstream through terminal customer information, including the terminal customer demand and changes of reorder-points of downstream retailers:(22)VCONt=CONt+dSSt=CONt+ROTt-ROTt-1.
(2)
Suppliers can check the inventory level of downstream retailers so that they can get the total supply chain inventory level which can optimize the total supply chain inventory decision:(23)SINVt=RINPt+VINVt.


## 3.2. VMI Distribution Supply Chain

In the actual operational process, suppliers can adopt VMI model on multiple downstream retailers. According to the aforementioned VMI serial supply chain when suppliers provide multiple retailers downstream with VMI service, they can get VMI distribution supply chain, as Figure6 shows, including one supplier and two retailers.Figure 6
Operation mode of VMI distribution supply chain.VMI-APIOBPCS II model as shown in Figure7 is constructed based on the operation mode in Figure 6.Figure 7
VMI-APIOBPCS II system dynamics model.


## 4. VMI&TPL-APIOBPCS Model

### 4.1. VMI Operational Process Based on TPL

#### 4.1.1. VMI&TPL Integrated Operation Mode

To decrease the logistics cost and avoid the delivery risks in VMI system, in practice an upstream enterprise normally prefers to outsource its purchasing business to the third-party logistics (TPL) and requires his supplier to keep the inventory in the warehouse operated by TPL. For example, BAX Global is responsible for Apple, Dell, IBM, and other IT companies with their supplies in Southeast Asia, and United Parcel Service manages goods and materials procurement for Fender overseas and achieves its integration of process in distribution. Besides, Shanghai Volkswagen and Wuhan Shenlong Automobile adopt VMI&TPL integrated operation mode to effectively support the mixed flow job shop manufacturing with JIT delivering components to the work station directly.After TPL is introduced into VMI, we consider the supply chain including one supplier (S), one TPL, and two retailers (R1 and R2). Suppliers give the rights of inventory operation and decision to TPL through a contract. TPL is responsible for replenishment and delivery in the total supply chain which means that TPL stores finished products in the warehouse near suppliers and, meanwhile, builds a district distribution center in order to meet retailers’ requirements in time. Besides, considering scale effect of transportation, TPL takes a certain delivery strategy in the district distribution center.Figure8 mainly describes VMI&TPL integrated operation mode. Similar to VMI operation mode, the information in the supply chain is shared fully. Retailers share real-time sales information with suppliers and TPL. Meanwhile, TPL shares inventory information with suppliers and suppliers provide TPL with production information. It is different from VMI operation mode that TPL is in charge of inventory management in the total supply chain, oversees the inventory level in the whole supply chain, and sends requests for replenishment to suppliers in time so that suppliers can organize capacities to produce according to orders. Besides, TPL organizes transit power and sends products to the district distribution center close to suppliers and to retailers and distributes products to retailers according to sales information and contracts in VMI&TPL integrated operation mode.Figure 8
VMI&TPL-APIOBPCS integrated operating model.


#### 4.1.2. VMI&TPL Operational Process

After carrying out TPL and VMI integrated operation, suppliers are not responsible for concrete logistics activities but give rights of inventory operation and decision to TPL through a contract. Therefore, TPL does not only undertake physical distribution business but is also responsible for orders generated in integrated logistics operation. The operational process of the whole system is illustrated in Figure9.(1)
TPL updates retailers’ inventory information everyday according to inventory information provided by retailers.(2)
TPL makes recommended orders according to retailers’ inventory level and service level and replenishment point confirmed in advance.(3)
TPL sends orders to retailers and chooses the proper distribution route according to self-inventory level and retailers’ demand after retailers confirm their orders. It is necessary to send requests for replenishment to suppliers and ask them to replenish inventory in time if self-inventory reaches replenishment level.(4)
Suppliers know well about the logistics operational situation by information sharing, then replenish inventory according to TPL’s demand, and settle accounts in time according to the orders confirmed by retailers.Figure 9
VMI&TPL integrated operation process.


### 4.2. VMI&TPL-APIOBPCS Model

The definitions of parameters and variables in VMI&TPL-APIOBPCS model are given in Table3.Table 3
Definitions of parameters and variables.
ORT: productivity
CRT: production fulfillment rate
RPT: replenishment rate
SRT: delivery rate
W-ROP: TPL replenishment point
D-ROP: TPL redelivery point
GIT: TPL transportation inventory
RINV: retailers’ temporary inventory levelAccording to the operational structure in Figure8 and operational process in Figure 9, three subsystems of VMI&TPL-APIOBPCS models, including suppliers’ production subsystem, TPL replenishment and delivery subsystem, and retailers’ sales subsystem, are analyzed as follows.


#### 4.2.1. Suppliers’ Production Subsystem

As shown in Figure10, suppliers’ production decisions are influenced by three aspects, which are demand information (terminal customer requirements), system inventory level, and work-in-process inventory level.Figure 10
Causality diagram of suppliers’ production subsystem.Difference equations of suppliers’ production operational process can be obtained according to the causality in Figure10 as shown in the following formulas.Work-in-process:(24)WIPt=WIPt-1+ORTt-CRTt.Work-in-process deviation:(25)EWIPt=DWIPt-WIPt.Productivity:(26)ORTt=AVCONt-1+EINVt-1Ti+EWIPt-1Tw.Productivity fulfillment rate:(27)CRTt=ORTt-Tp.System inventory deviation:(28)EINVt=TINV-SINVt.Target system inventory:(29)TINV=AVCON*(Tp+Tq).


#### 4.2.2. TPL Replenishment and Delivery Subsystem

TPL is responsible for replenishment and delivery decisions between suppliers and retailers. On the one hand, TPL sends requests for replenishment to suppliers in order to ensure proper inventory of the TPL warehouse (TPL-W). On the other hand, TPL needs to deliver products to the TPL distribution center (TPL-DC) near retailers in order to meet retailers’ demand. Figure11 is TPL replenishment and delivery subsystem causality diagram.Figure 11
TPL replenishment and delivery subsystem causality diagram.Formulas (30)–(32) are difference equations of TPL replenishment strategy and formulas (33)–(37) are difference equations of TPL delivery strategy. TPL total inventory level:(30)TPL-INVt=DC-INVt+TPLWt.TPL-W replenishment point:(31)W-ROPt=AVCONt×Gw.TPL replenishment capacity:(32)ROTt={W-ROPt-TPL-INVtifTPL-INVt<W-ROPt0ifTPL-INVt≥W-ROPt.TPL-DC redelivery point:(33)D-ROPt=AVCONt×Gd.TPL-DC inventory level:(34)DC-INVt=GITt+TPL-DCt.TPL transportation inventory level:(35)GITt=GITt-1+SPTt-SATt.(T, S) delivery strategy:
(36)SOTt={D-ROPt-DC-INVtifDC-INVt<D-ROPt0ifDC-INVt≥D-ROPt.(R, Q) delivery strategy:
(37)SOTt={n*ETQifDC-INVt<D-ROPt0ifDC-INVt≥D-ROPt.


#### 4.2.3. Retailers’ Sales Subsystem

Retailers’ sales subsystem is relatively simple. Though TPL can adopt JIT delivery to retailers, as service level should be improved, retailers still need to keep a small quantity of safety inventory and this inventory level is related to the demand fluctuation level. Figure12 is retailers’ sales subsystem causality diagram.Figure 12
Retailers’ sales subsystem causality diagram.According to the retailer’s sales system causal diagram, the difference equations for retailers can be obtained as follows:(38)R1=AVCON1t×G1,R2=AVCON2t×G2,JIT1t=R1t-RINV1t,JIT2t=R2t-RINV2t,CONSt=CON1t+CON2t,JITt=JIT1t+JIT2t.


#### 4.2.4. VMI&TPL-APIOBPCS System Dynamics Model

According to the aforementioned analysis, VMI&TPL-APIOBPCS system dynamics model is constructed as shown in Figure13.Figure 13
VMI&TPL-APIOBPCS model causality diagram.


## 4.1. VMI Operational Process Based on TPL

### 4.1.1. VMI&TPL Integrated Operation Mode

To decrease the logistics cost and avoid the delivery risks in VMI system, in practice an upstream enterprise normally prefers to outsource its purchasing business to the third-party logistics (TPL) and requires his supplier to keep the inventory in the warehouse operated by TPL. For example, BAX Global is responsible for Apple, Dell, IBM, and other IT companies with their supplies in Southeast Asia, and United Parcel Service manages goods and materials procurement for Fender overseas and achieves its integration of process in distribution. Besides, Shanghai Volkswagen and Wuhan Shenlong Automobile adopt VMI&TPL integrated operation mode to effectively support the mixed flow job shop manufacturing with JIT delivering components to the work station directly.After TPL is introduced into VMI, we consider the supply chain including one supplier (S), one TPL, and two retailers (R1 and R2). Suppliers give the rights of inventory operation and decision to TPL through a contract. TPL is responsible for replenishment and delivery in the total supply chain which means that TPL stores finished products in the warehouse near suppliers and, meanwhile, builds a district distribution center in order to meet retailers’ requirements in time. Besides, considering scale effect of transportation, TPL takes a certain delivery strategy in the district distribution center.Figure8 mainly describes VMI&TPL integrated operation mode. Similar to VMI operation mode, the information in the supply chain is shared fully. Retailers share real-time sales information with suppliers and TPL. Meanwhile, TPL shares inventory information with suppliers and suppliers provide TPL with production information. It is different from VMI operation mode that TPL is in charge of inventory management in the total supply chain, oversees the inventory level in the whole supply chain, and sends requests for replenishment to suppliers in time so that suppliers can organize capacities to produce according to orders. Besides, TPL organizes transit power and sends products to the district distribution center close to suppliers and to retailers and distributes products to retailers according to sales information and contracts in VMI&TPL integrated operation mode.Figure 8
VMI&TPL-APIOBPCS integrated operating model.


### 4.1.2. VMI&TPL Operational Process

After carrying out TPL and VMI integrated operation, suppliers are not responsible for concrete logistics activities but give rights of inventory operation and decision to TPL through a contract. Therefore, TPL does not only undertake physical distribution business but is also responsible for orders generated in integrated logistics operation. The operational process of the whole system is illustrated in Figure9.(1)
TPL updates retailers’ inventory information everyday according to inventory information provided by retailers.(2)
TPL makes recommended orders according to retailers’ inventory level and service level and replenishment point confirmed in advance.(3)
TPL sends orders to retailers and chooses the proper distribution route according to self-inventory level and retailers’ demand after retailers confirm their orders. It is necessary to send requests for replenishment to suppliers and ask them to replenish inventory in time if self-inventory reaches replenishment level.(4)
Suppliers know well about the logistics operational situation by information sharing, then replenish inventory according to TPL’s demand, and settle accounts in time according to the orders confirmed by retailers.Figure 9
VMI&TPL integrated operation process.


## 4.1.1. VMI&TPL Integrated Operation Mode

To decrease the logistics cost and avoid the delivery risks in VMI system, in practice an upstream enterprise normally prefers to outsource its purchasing business to the third-party logistics (TPL) and requires his supplier to keep the inventory in the warehouse operated by TPL. For example, BAX Global is responsible for Apple, Dell, IBM, and other IT companies with their supplies in Southeast Asia, and United Parcel Service manages goods and materials procurement for Fender overseas and achieves its integration of process in distribution. Besides, Shanghai Volkswagen and Wuhan Shenlong Automobile adopt VMI&TPL integrated operation mode to effectively support the mixed flow job shop manufacturing with JIT delivering components to the work station directly.After TPL is introduced into VMI, we consider the supply chain including one supplier (S), one TPL, and two retailers (R1 and R2). Suppliers give the rights of inventory operation and decision to TPL through a contract. TPL is responsible for replenishment and delivery in the total supply chain which means that TPL stores finished products in the warehouse near suppliers and, meanwhile, builds a district distribution center in order to meet retailers’ requirements in time. Besides, considering scale effect of transportation, TPL takes a certain delivery strategy in the district distribution center.Figure8 mainly describes VMI&TPL integrated operation mode. Similar to VMI operation mode, the information in the supply chain is shared fully. Retailers share real-time sales information with suppliers and TPL. Meanwhile, TPL shares inventory information with suppliers and suppliers provide TPL with production information. It is different from VMI operation mode that TPL is in charge of inventory management in the total supply chain, oversees the inventory level in the whole supply chain, and sends requests for replenishment to suppliers in time so that suppliers can organize capacities to produce according to orders. Besides, TPL organizes transit power and sends products to the district distribution center close to suppliers and to retailers and distributes products to retailers according to sales information and contracts in VMI&TPL integrated operation mode.Figure 8
VMI&TPL-APIOBPCS integrated operating model.


## 4.1.2. VMI&TPL Operational Process

After carrying out TPL and VMI integrated operation, suppliers are not responsible for concrete logistics activities but give rights of inventory operation and decision to TPL through a contract. Therefore, TPL does not only undertake physical distribution business but is also responsible for orders generated in integrated logistics operation. The operational process of the whole system is illustrated in Figure9.(1)
TPL updates retailers’ inventory information everyday according to inventory information provided by retailers.(2)
TPL makes recommended orders according to retailers’ inventory level and service level and replenishment point confirmed in advance.(3)
TPL sends orders to retailers and chooses the proper distribution route according to self-inventory level and retailers’ demand after retailers confirm their orders. It is necessary to send requests for replenishment to suppliers and ask them to replenish inventory in time if self-inventory reaches replenishment level.(4)
Suppliers know well about the logistics operational situation by information sharing, then replenish inventory according to TPL’s demand, and settle accounts in time according to the orders confirmed by retailers.Figure 9
VMI&TPL integrated operation process.


## 4.2. VMI&TPL-APIOBPCS Model

The definitions of parameters and variables in VMI&TPL-APIOBPCS model are given in Table3.Table 3
Definitions of parameters and variables.
ORT: productivity
CRT: production fulfillment rate
RPT: replenishment rate
SRT: delivery rate
W-ROP: TPL replenishment point
D-ROP: TPL redelivery point
GIT: TPL transportation inventory
RINV: retailers’ temporary inventory levelAccording to the operational structure in Figure8 and operational process in Figure 9, three subsystems of VMI&TPL-APIOBPCS models, including suppliers’ production subsystem, TPL replenishment and delivery subsystem, and retailers’ sales subsystem, are analyzed as follows.


### 4.2.1. Suppliers’ Production Subsystem

As shown in Figure10, suppliers’ production decisions are influenced by three aspects, which are demand information (terminal customer requirements), system inventory level, and work-in-process inventory level.Figure 10
Causality diagram of suppliers’ production subsystem.Difference equations of suppliers’ production operational process can be obtained according to the causality in Figure10 as shown in the following formulas.Work-in-process:(24)WIPt=WIPt-1+ORTt-CRTt.Work-in-process deviation:(25)EWIPt=DWIPt-WIPt.Productivity:(26)ORTt=AVCONt-1+EINVt-1Ti+EWIPt-1Tw.Productivity fulfillment rate:(27)CRTt=ORTt-Tp.System inventory deviation:(28)EINVt=TINV-SINVt.Target system inventory:(29)TINV=AVCON*(Tp+Tq).


### 4.2.2. TPL Replenishment and Delivery Subsystem

TPL is responsible for replenishment and delivery decisions between suppliers and retailers. On the one hand, TPL sends requests for replenishment to suppliers in order to ensure proper inventory of the TPL warehouse (TPL-W). On the other hand, TPL needs to deliver products to the TPL distribution center (TPL-DC) near retailers in order to meet retailers’ demand. Figure11 is TPL replenishment and delivery subsystem causality diagram.Figure 11
TPL replenishment and delivery subsystem causality diagram.Formulas (30)–(32) are difference equations of TPL replenishment strategy and formulas (33)–(37) are difference equations of TPL delivery strategy. TPL total inventory level:(30)TPL-INVt=DC-INVt+TPLWt.TPL-W replenishment point:(31)W-ROPt=AVCONt×Gw.TPL replenishment capacity:(32)ROTt={W-ROPt-TPL-INVtifTPL-INVt<W-ROPt0ifTPL-INVt≥W-ROPt.TPL-DC redelivery point:(33)D-ROPt=AVCONt×Gd.TPL-DC inventory level:(34)DC-INVt=GITt+TPL-DCt.TPL transportation inventory level:(35)GITt=GITt-1+SPTt-SATt.(T, S) delivery strategy:
(36)SOTt={D-ROPt-DC-INVtifDC-INVt<D-ROPt0ifDC-INVt≥D-ROPt.(R, Q) delivery strategy:
(37)SOTt={n*ETQifDC-INVt<D-ROPt0ifDC-INVt≥D-ROPt.


### 4.2.3. Retailers’ Sales Subsystem

Retailers’ sales subsystem is relatively simple. Though TPL can adopt JIT delivery to retailers, as service level should be improved, retailers still need to keep a small quantity of safety inventory and this inventory level is related to the demand fluctuation level. Figure12 is retailers’ sales subsystem causality diagram.Figure 12
Retailers’ sales subsystem causality diagram.According to the retailer’s sales system causal diagram, the difference equations for retailers can be obtained as follows:(38)R1=AVCON1t×G1,R2=AVCON2t×G2,JIT1t=R1t-RINV1t,JIT2t=R2t-RINV2t,CONSt=CON1t+CON2t,JITt=JIT1t+JIT2t.


### 4.2.4. VMI&TPL-APIOBPCS System Dynamics Model

According to the aforementioned analysis, VMI&TPL-APIOBPCS system dynamics model is constructed as shown in Figure13.Figure 13
VMI&TPL-APIOBPCS model causality diagram.


## 4.2.1. Suppliers’ Production Subsystem

As shown in Figure10, suppliers’ production decisions are influenced by three aspects, which are demand information (terminal customer requirements), system inventory level, and work-in-process inventory level.Figure 10
Causality diagram of suppliers’ production subsystem.Difference equations of suppliers’ production operational process can be obtained according to the causality in Figure10 as shown in the following formulas.Work-in-process:(24)WIPt=WIPt-1+ORTt-CRTt.Work-in-process deviation:(25)EWIPt=DWIPt-WIPt.Productivity:(26)ORTt=AVCONt-1+EINVt-1Ti+EWIPt-1Tw.Productivity fulfillment rate:(27)CRTt=ORTt-Tp.System inventory deviation:(28)EINVt=TINV-SINVt.Target system inventory:(29)TINV=AVCON*(Tp+Tq).


## 4.2.2. TPL Replenishment and Delivery Subsystem

TPL is responsible for replenishment and delivery decisions between suppliers and retailers. On the one hand, TPL sends requests for replenishment to suppliers in order to ensure proper inventory of the TPL warehouse (TPL-W). On the other hand, TPL needs to deliver products to the TPL distribution center (TPL-DC) near retailers in order to meet retailers’ demand. Figure11 is TPL replenishment and delivery subsystem causality diagram.Figure 11
TPL replenishment and delivery subsystem causality diagram.Formulas (30)–(32) are difference equations of TPL replenishment strategy and formulas (33)–(37) are difference equations of TPL delivery strategy. TPL total inventory level:(30)TPL-INVt=DC-INVt+TPLWt.TPL-W replenishment point:(31)W-ROPt=AVCONt×Gw.TPL replenishment capacity:(32)ROTt={W-ROPt-TPL-INVtifTPL-INVt<W-ROPt0ifTPL-INVt≥W-ROPt.TPL-DC redelivery point:(33)D-ROPt=AVCONt×Gd.TPL-DC inventory level:(34)DC-INVt=GITt+TPL-DCt.TPL transportation inventory level:(35)GITt=GITt-1+SPTt-SATt.(T, S) delivery strategy:
(36)SOTt={D-ROPt-DC-INVtifDC-INVt<D-ROPt0ifDC-INVt≥D-ROPt.(R, Q) delivery strategy:
(37)SOTt={n*ETQifDC-INVt<D-ROPt0ifDC-INVt≥D-ROPt.


## 4.2.3. Retailers’ Sales Subsystem

Retailers’ sales subsystem is relatively simple. Though TPL can adopt JIT delivery to retailers, as service level should be improved, retailers still need to keep a small quantity of safety inventory and this inventory level is related to the demand fluctuation level. Figure12 is retailers’ sales subsystem causality diagram.Figure 12
Retailers’ sales subsystem causality diagram.According to the retailer’s sales system causal diagram, the difference equations for retailers can be obtained as follows:(38)R1=AVCON1t×G1,R2=AVCON2t×G2,JIT1t=R1t-RINV1t,JIT2t=R2t-RINV2t,CONSt=CON1t+CON2t,JITt=JIT1t+JIT2t.


## 4.2.4. VMI&TPL-APIOBPCS System Dynamics Model

According to the aforementioned analysis, VMI&TPL-APIOBPCS system dynamics model is constructed as shown in Figure13.Figure 13
VMI&TPL-APIOBPCS model causality diagram.


## 5. Simulation Analysis

According to the aforementioned VMI&TPL-APIOBPC system dynamics model, the two different conditions with phase-step and random demands are investigated, respectively, and parameter settings are as follows.(1)
Production subsystem parameters settings. Referring to that of Disney and Towill [24, 25], Ta=4, Ti=14, Tw=8, and Tp=5.(2)
Parameters settings of replenishment and delivery subsystem,Gw=8, Gd=5, ETQ=10, and Tq=5.(3)
Parameters settings of sales subsystem,G1=3, G2=3.Then simulate the two models, VMI-APIOBPCS II and VMI&TPL-APIOBPCS, using Vensim, and run the test for 100 units of time (month).


### 5.1. Phase-Step Demand Test

The demand test functions CONS1 and CONS2 are both phase-step functions,CONS1=STEP  (5,0)  (STEP({height},{stime})), CONS2=STEP  (10,10).


#### 5.1.1. Fluctuations of Production Order

As shown in Figure14, when TPL is introduced into VMI, suppliers’ productivity becomes smoother, which is mostly due to the smoother system inventory levels. In the meanwhile, the response time of productivity reaching the steady-state in the VMI&TPL-APIOBPCS model is relatively longer.Figure 14
Comparison of production fluctuations in suppliers’ production orders.


#### 5.1.2. Fluctuations of Inventory Level

(1)
System inventory levels. As shown in Figure15, in VMI&TPL-APIOBPCS model, the system inventory level significantly decreases and becomes smoother. Under phase-step demand, the two model system inventories in models gradually return to a steady-state value (TINV = 150) and are similar to the suppliers’ productivity. The response time of system inventory level reaching the steady-state in VMI&TPL-APIOBPCS model is relatively longer.
(2)
Suppliers’ inventory level. Figure16 is the suppliers’ inventory level. After introducing TPL, suppliers’ inventory level is effectively smoothed and reduced, which is because that the fluctuation of downstream replenishment batch is more consistent compared to that in VMI model; see Figure 17.
(3)
Comparison of downstream overall inventory level. In VMI-APIOBPCS model, the downstream inventory mainly includes retailers’ inventories. However, in the VMI&TPL-APIOBPCS model, the downstream inventory includes not only the retailers’ safety stock but also the inventory in TPL warehouse and distribution center. Figure18 is the downstream inventory level. After introducing TPL, the downstream inventory level (including the TPL) increases slightly, largely because the downstream structure is added with a new subject (TPL).Figure 15
The comparison of system inventory level fluctuation.Figure 16
The comparison of supplier inventory level fluctuation.Figure 17
The comparison of replenishment batch fluctuation.Figure 18
The comparison of downstream inventory level fluctuation.From the aforementioned inventory analysis, we can know that after introducing TPL into VMI supply chain, although supply chains add an echelon, the downstream inventory level increases slightly, but the replenishment batch becomes smoother after TPL’s participation, thus effectively smoothing suppliers’ productivity and reducing supplier inventory levels.


#### 5.1.3. Comparison of Service Level

What defines service level of retailers is the ratio between the retailers’ inventory level and customers’ demand. Figures19 and 20 are the service level of the two operation modes, respectively. Since the initial state of the system is zero, the initial service level of the retailers is zero, and the service level gradually improves after a period of time. In VMI operation mode, the average service levels of two retailers are 61.5% and 46.2%, respectively. In VMI&TPL integrated operation mode, the average service level of two retailers is 58.5%.Figure 19
Service levels in the VMI operation mode.Figure 20
Service levels in the VMI&TPL integrated operation mode.With the deterministic demand, service level of the first retailer falls slightly, but the second retailer's service level increases substantially, and the overall average service level of retailers improves after introducing TPL. However, the response time of the system increases, because the system service level will not be enhanced till a long period of shortage.


### 5.2. Random Demand Test

The demand test functions, such as CONS1 and CONS2, are both phase-step functions.CONS1=RANDOMNORMAL (0,10,5,20,1) (RANDOM NORMAL ({min},{max},{mean},{stdev},{seed})), CONS2=RANDOMNORMAL (0,5,3,10,1).


#### 5.2.1. Suppliers’ Productivity

Similar to the phase-step demand situation, supplier’s productivity is steadier apparently in VMI&TPL-APIOBPCS mode.


#### 5.2.2. Comparison of Inventory Levels

(1)
System inventory level. As Figure22 shows, compared with the VMI model, the system inventory level in VMI&TPL-APIOBPCS model is lower and largely benefits from the dramatical decrease of suppliers’ inventory level, as shown in Figure 23.
(2)
Suppliers’ inventory levels. As Figure23 shows, in VMI&TPL-APIOBPCS model, suppliers’ inventory levels are significantly reduced as suppliers’ productivity becomes smoother, and TPL replenishment batch is steadier compared to VMI model.
(3)
Downstream inventory level: As Figure24 shows, in VMI&TPL-APIOBPCS model, the downstream inventory level increases slightly due to an additional echelon of TPL with concentrated restocking and delivery in supply chain.


#### 5.2.3. Comparison of Service Level

Under the random demand and in the operation mode of VMI the first retailer’s average service level is 78.6% and the second retailer’s average service level is 40.4%.In the VMI&TPL integrated operation mode, the overall service levels of the two retailers are 94.8%. It is clear that after introducing TPL the two service levels of retailers are improved under random demand; this benefits from the risk-sharing effect in centralized inventory after the introduction of TPL.


### 5.3. Simulation Discussion

The complex simulations under phase-step demand and random demand show the following results as this study summarized from the follow-up interviews in practice.First, VMI&TPL integrated operation mode can smoothen suppliers’ productivity under both phase-step demand and random demand as TPL introduced into VMI makes the whole system inventory levels smoother (see Figure21). However, the response time of reaching the steady-state in the VMI&TPL-APIOBPCS model is relatively longer, which may be caused by TPL centralized operation. Particularly, the fluctuation of suppliers’ productivity under random demand is more stable than under phase-step demand, as in practice TPL places orders from suppliers periodically, and its response to fluctuation is slow.Figure 21
Fluctuation diagram of suppliers’ productivity.Figure 22
Fluctuation diagram of system inventory level.Figure 23
Fluctuation diagram of suppliers’ inventory level.Figure 24
Fluctuation diagram of lower inventory level.Second, VMI&TPL integrated operation mode can reduce the system inventory level significantly. Similarly, the response time of system inventory level reaching the steady-state in VMI&TPL-APIOBPCS model is relatively longer. Compared with under phase-step demand, system inventory level can be lower under random demand. This illustrates that the TPL centralized replenishment has a scale of economics and much risk pooling effects which can decrease the whole system inventory level [37]. Besides, suppliers’ inventory level is effectively smoothed and reduced since downstream replenishment batch is more consistent and scale of economics compared to those in VMI model. However, in the VMI&TPL-APIOBPCS model, the downstream inventory includes not only the retailers’ safety stock but also the inventory in TPL warehouse and distribution center. As a result, the downstream inventory level rises slightly after introducing TPL.Third, VMI&TPL integrated operation mode can improve average service level. Under phase-step demand, service level of the first retailer falls slightly, but the second retailer's service level increases substantially. As a contrast, service level of two retailers increases under random demand. On the whole, the service level under random demand is improved significantly than under phase-step demand (see Figures25 and 26). These may be caused by risk pooling effect, especially under random demand.Figure 25
The service level of VMI operational model under random demand.Figure 26
The service level of VMI&TPL integrated operational model under random demand.


## 5.1. Phase-Step Demand Test

The demand test functions CONS1 and CONS2 are both phase-step functions,CONS1=STEP  (5,0)  (STEP({height},{stime})), CONS2=STEP  (10,10).


### 5.1.1. Fluctuations of Production Order

As shown in Figure14, when TPL is introduced into VMI, suppliers’ productivity becomes smoother, which is mostly due to the smoother system inventory levels. In the meanwhile, the response time of productivity reaching the steady-state in the VMI&TPL-APIOBPCS model is relatively longer.Figure 14
Comparison of production fluctuations in suppliers’ production orders.


### 5.1.2. Fluctuations of Inventory Level

(1)
System inventory levels. As shown in Figure15, in VMI&TPL-APIOBPCS model, the system inventory level significantly decreases and becomes smoother. Under phase-step demand, the two model system inventories in models gradually return to a steady-state value (TINV = 150) and are similar to the suppliers’ productivity. The response time of system inventory level reaching the steady-state in VMI&TPL-APIOBPCS model is relatively longer.
(2)
Suppliers’ inventory level. Figure16 is the suppliers’ inventory level. After introducing TPL, suppliers’ inventory level is effectively smoothed and reduced, which is because that the fluctuation of downstream replenishment batch is more consistent compared to that in VMI model; see Figure 17.
(3)
Comparison of downstream overall inventory level. In VMI-APIOBPCS model, the downstream inventory mainly includes retailers’ inventories. However, in the VMI&TPL-APIOBPCS model, the downstream inventory includes not only the retailers’ safety stock but also the inventory in TPL warehouse and distribution center. Figure18 is the downstream inventory level. After introducing TPL, the downstream inventory level (including the TPL) increases slightly, largely because the downstream structure is added with a new subject (TPL).Figure 15
The comparison of system inventory level fluctuation.Figure 16
The comparison of supplier inventory level fluctuation.Figure 17
The comparison of replenishment batch fluctuation.Figure 18
The comparison of downstream inventory level fluctuation.From the aforementioned inventory analysis, we can know that after introducing TPL into VMI supply chain, although supply chains add an echelon, the downstream inventory level increases slightly, but the replenishment batch becomes smoother after TPL’s participation, thus effectively smoothing suppliers’ productivity and reducing supplier inventory levels.


### 5.1.3. Comparison of Service Level

What defines service level of retailers is the ratio between the retailers’ inventory level and customers’ demand. Figures19 and 20 are the service level of the two operation modes, respectively. Since the initial state of the system is zero, the initial service level of the retailers is zero, and the service level gradually improves after a period of time. In VMI operation mode, the average service levels of two retailers are 61.5% and 46.2%, respectively. In VMI&TPL integrated operation mode, the average service level of two retailers is 58.5%.Figure 19
Service levels in the VMI operation mode.Figure 20
Service levels in the VMI&TPL integrated operation mode.With the deterministic demand, service level of the first retailer falls slightly, but the second retailer's service level increases substantially, and the overall average service level of retailers improves after introducing TPL. However, the response time of the system increases, because the system service level will not be enhanced till a long period of shortage.


## 5.1.1. Fluctuations of Production Order

As shown in Figure14, when TPL is introduced into VMI, suppliers’ productivity becomes smoother, which is mostly due to the smoother system inventory levels. In the meanwhile, the response time of productivity reaching the steady-state in the VMI&TPL-APIOBPCS model is relatively longer.Figure 14
Comparison of production fluctuations in suppliers’ production orders.


## 5.1.2. Fluctuations of Inventory Level

(1)
System inventory levels. As shown in Figure15, in VMI&TPL-APIOBPCS model, the system inventory level significantly decreases and becomes smoother. Under phase-step demand, the two model system inventories in models gradually return to a steady-state value (TINV = 150) and are similar to the suppliers’ productivity. The response time of system inventory level reaching the steady-state in VMI&TPL-APIOBPCS model is relatively longer.
(2)
Suppliers’ inventory level. Figure16 is the suppliers’ inventory level. After introducing TPL, suppliers’ inventory level is effectively smoothed and reduced, which is because that the fluctuation of downstream replenishment batch is more consistent compared to that in VMI model; see Figure 17.
(3)
Comparison of downstream overall inventory level. In VMI-APIOBPCS model, the downstream inventory mainly includes retailers’ inventories. However, in the VMI&TPL-APIOBPCS model, the downstream inventory includes not only the retailers’ safety stock but also the inventory in TPL warehouse and distribution center. Figure18 is the downstream inventory level. After introducing TPL, the downstream inventory level (including the TPL) increases slightly, largely because the downstream structure is added with a new subject (TPL).Figure 15
The comparison of system inventory level fluctuation.Figure 16
The comparison of supplier inventory level fluctuation.Figure 17
The comparison of replenishment batch fluctuation.Figure 18
The comparison of downstream inventory level fluctuation.From the aforementioned inventory analysis, we can know that after introducing TPL into VMI supply chain, although supply chains add an echelon, the downstream inventory level increases slightly, but the replenishment batch becomes smoother after TPL’s participation, thus effectively smoothing suppliers’ productivity and reducing supplier inventory levels.


## 5.1.3. Comparison of Service Level

What defines service level of retailers is the ratio between the retailers’ inventory level and customers’ demand. Figures19 and 20 are the service level of the two operation modes, respectively. Since the initial state of the system is zero, the initial service level of the retailers is zero, and the service level gradually improves after a period of time. In VMI operation mode, the average service levels of two retailers are 61.5% and 46.2%, respectively. In VMI&TPL integrated operation mode, the average service level of two retailers is 58.5%.Figure 19
Service levels in the VMI operation mode.Figure 20
Service levels in the VMI&TPL integrated operation mode.With the deterministic demand, service level of the first retailer falls slightly, but the second retailer's service level increases substantially, and the overall average service level of retailers improves after introducing TPL. However, the response time of the system increases, because the system service level will not be enhanced till a long period of shortage.


## 5.2. Random Demand Test

The demand test functions, such as CONS1 and CONS2, are both phase-step functions.CONS1=RANDOMNORMAL (0,10,5,20,1) (RANDOM NORMAL ({min},{max},{mean},{stdev},{seed})), CONS2=RANDOMNORMAL (0,5,3,10,1).


### 5.2.1. Suppliers’ Productivity

Similar to the phase-step demand situation, supplier’s productivity is steadier apparently in VMI&TPL-APIOBPCS mode.


### 5.2.2. Comparison of Inventory Levels

(1)
System inventory level. As Figure22 shows, compared with the VMI model, the system inventory level in VMI&TPL-APIOBPCS model is lower and largely benefits from the dramatical decrease of suppliers’ inventory level, as shown in Figure 23.
(2)
Suppliers’ inventory levels. As Figure23 shows, in VMI&TPL-APIOBPCS model, suppliers’ inventory levels are significantly reduced as suppliers’ productivity becomes smoother, and TPL replenishment batch is steadier compared to VMI model.
(3)
Downstream inventory level: As Figure24 shows, in VMI&TPL-APIOBPCS model, the downstream inventory level increases slightly due to an additional echelon of TPL with concentrated restocking and delivery in supply chain.


### 5.2.3. Comparison of Service Level

Under the random demand and in the operation mode of VMI the first retailer’s average service level is 78.6% and the second retailer’s average service level is 40.4%.In the VMI&TPL integrated operation mode, the overall service levels of the two retailers are 94.8%. It is clear that after introducing TPL the two service levels of retailers are improved under random demand; this benefits from the risk-sharing effect in centralized inventory after the introduction of TPL.


## 5.2.1. Suppliers’ Productivity

Similar to the phase-step demand situation, supplier’s productivity is steadier apparently in VMI&TPL-APIOBPCS mode.


## 5.2.2. Comparison of Inventory Levels

(1)
System inventory level. As Figure22 shows, compared with the VMI model, the system inventory level in VMI&TPL-APIOBPCS model is lower and largely benefits from the dramatical decrease of suppliers’ inventory level, as shown in Figure 23.
(2)
Suppliers’ inventory levels. As Figure23 shows, in VMI&TPL-APIOBPCS model, suppliers’ inventory levels are significantly reduced as suppliers’ productivity becomes smoother, and TPL replenishment batch is steadier compared to VMI model.
(3)
Downstream inventory level: As Figure24 shows, in VMI&TPL-APIOBPCS model, the downstream inventory level increases slightly due to an additional echelon of TPL with concentrated restocking and delivery in supply chain.


## 5.2.3. Comparison of Service Level

Under the random demand and in the operation mode of VMI the first retailer’s average service level is 78.6% and the second retailer’s average service level is 40.4%.In the VMI&TPL integrated operation mode, the overall service levels of the two retailers are 94.8%. It is clear that after introducing TPL the two service levels of retailers are improved under random demand; this benefits from the risk-sharing effect in centralized inventory after the introduction of TPL.


## 5.3. Simulation Discussion

The complex simulations under phase-step demand and random demand show the following results as this study summarized from the follow-up interviews in practice.First, VMI&TPL integrated operation mode can smoothen suppliers’ productivity under both phase-step demand and random demand as TPL introduced into VMI makes the whole system inventory levels smoother (see Figure21). However, the response time of reaching the steady-state in the VMI&TPL-APIOBPCS model is relatively longer, which may be caused by TPL centralized operation. Particularly, the fluctuation of suppliers’ productivity under random demand is more stable than under phase-step demand, as in practice TPL places orders from suppliers periodically, and its response to fluctuation is slow.Figure 21
Fluctuation diagram of suppliers’ productivity.Figure 22
Fluctuation diagram of system inventory level.Figure 23
Fluctuation diagram of suppliers’ inventory level.Figure 24
Fluctuation diagram of lower inventory level.Second, VMI&TPL integrated operation mode can reduce the system inventory level significantly. Similarly, the response time of system inventory level reaching the steady-state in VMI&TPL-APIOBPCS model is relatively longer. Compared with under phase-step demand, system inventory level can be lower under random demand. This illustrates that the TPL centralized replenishment has a scale of economics and much risk pooling effects which can decrease the whole system inventory level [37]. Besides, suppliers’ inventory level is effectively smoothed and reduced since downstream replenishment batch is more consistent and scale of economics compared to those in VMI model. However, in the VMI&TPL-APIOBPCS model, the downstream inventory includes not only the retailers’ safety stock but also the inventory in TPL warehouse and distribution center. As a result, the downstream inventory level rises slightly after introducing TPL.Third, VMI&TPL integrated operation mode can improve average service level. Under phase-step demand, service level of the first retailer falls slightly, but the second retailer's service level increases substantially. As a contrast, service level of two retailers increases under random demand. On the whole, the service level under random demand is improved significantly than under phase-step demand (see Figures25 and 26). These may be caused by risk pooling effect, especially under random demand.Figure 25
The service level of VMI operational model under random demand.Figure 26
The service level of VMI&TPL integrated operational model under random demand.


## 6. Conclusion

This paper constructs the VMI&TPL-APIOBPCS model after introducing TPL into VMI distribution based on VMI-APIOBPCS system dynamics model. The system performance of VMI&TPL integrated supply chain under phase-step and random demand is considered. The simulation analysis shows that though the supply chain is turned into a three-echelon structure from a two-echelon one TPL can effectively smooth the replenishment and delivery quantity between suppliers and retailers by goods collection, thus dramatically reducing the inventory level of the suppliers and the whole system, effectively smoothing the production rhythm of suppliers and improving the service level of customers.Although system dynamics method can describe and simulate VMI&TPL integrated operational model, it lacks the optimization of TPL replenishment and delivery policy in this operation mode. As a result, it is necessary to optimize the TPL replenishment and delivery policy under various demands using mathematical programming and optimization theories, so as to enhance the study of VMI&TPL integrated operational model further.



---

*Source: 178713-2013-11-25.xml*