趙　軍　劉　飛

摘要：利用系統(tǒng)動力學(xué)的方法和仿真軟件Vensim，對Y公司K產(chǎn)品的供應(yīng)鏈上的牛鞭效應(yīng)進(jìn)行了建模與仿真。針對該供應(yīng)鏈上牛鞭效應(yīng)顯著的實(shí)際情況，提出了縮短訂單延遲時間，延長庫存調(diào)節(jié)時間，采用供應(yīng)商管理庫存的綜合解決方案，仿真結(jié)果表明該解決方案能在一定程度上減輕牛鞭效應(yīng)的影響，研究結(jié)果可為Y公司的實(shí)際生產(chǎn)組織提供決策依據(jù)。

關(guān)鍵詞：系統(tǒng)動力學(xué)；供應(yīng)鏈；牛鞭效應(yīng)；仿真；供應(yīng)商管理庫存

中圖分類號：F273.7文獻(xiàn)標(biāo)識碼：A

Abstract: This paper adopts the system dynamics and its simulation software Vensim to modeling and simulating the bullwhip effect in the supply chain of the Y company's k product. Based on the reality situation of the bullwhip effect in this supply chain, it brought forward shortening delay time of order, prolonging adjust time of inventory, adopting vendor managed inventory to solve this problem, the simulation results proved that these method could alleviate the effect of the bullwhip effect in some extend, the research results could enrich the decision about the reality product organization for company Y.

Key words: system dynamics; supply chain; bullwhip effect; simulation; vendor managed inventory

牛鞭效應(yīng)指供應(yīng)鏈中需求信息從下游企業(yè)向上游企業(yè)傳遞時，信息被扭曲并逐級放大的現(xiàn)象[1]。系統(tǒng)動力學(xué)出現(xiàn)于1956年, 該方法在20世紀(jì)50年代后期主要被用于企業(yè)管理，解決諸如原材料供應(yīng)、生產(chǎn)、庫存、銷售、市場等一系列問題[1]。劉媛媛等[2]應(yīng)用系統(tǒng)動力學(xué)研究了供應(yīng)鏈節(jié)點(diǎn)變化對牛鞭效應(yīng)的影響，利用系統(tǒng)動力學(xué)專用軟件VENSIM，對“啤酒游戲”進(jìn)行仿真分析，研究當(dāng)供應(yīng)鏈上節(jié)點(diǎn)增加時“牛鞭效應(yīng)”的變化情況；李穩(wěn)安等[3]用系統(tǒng)動力學(xué)原理研究了供應(yīng)鏈中牛鞭效應(yīng)產(chǎn)生的原因及相應(yīng)的緩解對策，建立了多節(jié)點(diǎn)供應(yīng)鏈系統(tǒng)的動力學(xué)模型，給出了系統(tǒng)動力學(xué)的迭代表達(dá)式和框圖，并定性地說明了牛鞭效應(yīng)產(chǎn)生的內(nèi)在機(jī)理；王磊等[4]對研究牛鞭效應(yīng)時經(jīng)常討論到的三種模型：系統(tǒng)動力學(xué)模型、自回歸分析的AR（1）模型和Kalman濾波器（KF）模型進(jìn)行了分析和比較。楊詩華[5]在其碩士論文中基于系統(tǒng)動力學(xué)原理，應(yīng)用Powersim仿真軟件對啤酒分銷模型的三級供應(yīng)鏈中牛鞭效應(yīng)現(xiàn)象進(jìn)行了仿真設(shè)計，提出緩解牛鞭效應(yīng)的對策；苗興東[6]在其碩士論文中采用系統(tǒng)動力學(xué)研究牛鞭效應(yīng)的機(jī)理，利用系統(tǒng)動力學(xué)專用軟件VENSIM，首先研究極端情況下的牛鞭效應(yīng)（包括曲線和數(shù)據(jù)計算），然后通過調(diào)整參數(shù)來減弱牛鞭效應(yīng)，從而得出結(jié)論。本文以一家生產(chǎn)化妝品系列的外貿(mào)出口企業(yè)Y公司為背景，首先運(yùn)用系統(tǒng)動力學(xué)的方法對該公司K產(chǎn)品供應(yīng)鏈上的牛鞭效應(yīng)進(jìn)行建模，其次應(yīng)用系統(tǒng)動力學(xué)仿真軟件Vensim PLE 5.4b對該牛鞭效益進(jìn)行模擬仿真，最后針對該公司K產(chǎn)品的供應(yīng)鏈上存在的問題，提出改進(jìn)措施，以減輕牛鞭效應(yīng)的影響。

1Y公司介紹

Y公司是一家生產(chǎn)化妝品系列的外貿(mào)出口型企業(yè)，主要承接英國、法國及美國等國家的生產(chǎn)訂單，按照客戶的要求進(jìn)行生產(chǎn)。其產(chǎn)品品種多，顏色豐富，款式也非常多，甚至曾經(jīng)生產(chǎn)過的某些產(chǎn)品，其款式無一相同。Y公司K產(chǎn)品的供應(yīng)鏈由消費(fèi)者c，零售商d、批發(fā)商W和Y公司所組成。該供應(yīng)鏈上的參加人員各自扮演不同角色，且完全獨(dú)立地做出決策，他們唯一的目標(biāo)是使各自決策的利潤最大化，即各個物流環(huán)節(jié)的訂貨策略都是：Qij=SIij-Iij/Ti+Nj，其中：Qij, SIij, Iij分別為第i個物流環(huán)節(jié)第j周的訂貨量、安全庫存量和期末庫存量（在缺貨情況下取值為0）；Ti為第i個物流環(huán)節(jié)庫存的調(diào)節(jié)時間；Nj為第j周期下游物流單元的需求量。已知消費(fèi)者的需求是一個階躍函數(shù)，第一周市場需求率是4件/周，從第二周開始每周市場需求率都是8件/周。零售商的安全庫存量為12件，庫存調(diào)節(jié)時間為4周，期末庫存量初始值為12件，累計缺貨量初始值為0件，與消費(fèi)者間沒有時間延遲。批發(fā)商的安全庫存量為12件，庫存調(diào)節(jié)時間為4周，期末庫存量初始值為12件，累計缺貨量初始值為0件，與零售商間有4周的時間延遲。制造商的安全庫存量為12件，庫存調(diào)節(jié)時間為4周，期末庫存量初始值為12件，累計缺貨量初始值為0件，與批發(fā)商商間有4周的時間延遲，生產(chǎn)有2周的時間延遲。

2系統(tǒng)動力學(xué)建模與仿真

2.1K產(chǎn)品供應(yīng)鏈的因果圖

在K產(chǎn)品的供應(yīng)鏈中，有三個物流單元，包括：制造商Y公司、批發(fā)商w、零售商d，而各個物流單元之間用通訊和運(yùn)輸相聯(lián)系。根據(jù)各個環(huán)節(jié)的關(guān)系，畫出K產(chǎn)品供應(yīng)鏈的因果關(guān)系圖，如圖1所示。

2.2K產(chǎn)品供應(yīng)鏈的系統(tǒng)動力學(xué)流圖

根據(jù)圖1的因果關(guān)系圖，可進(jìn)一步得出K產(chǎn)品供應(yīng)鏈的系統(tǒng)動力學(xué)流圖，如圖2所示。這里確定六個水平變量：制造商期末庫存量，批發(fā)商期末庫存量，零售商期末庫存量，制造商累計缺貨量，批發(fā)商累計缺貨量，零售商累計缺貨量；八個速率變量：制造商的產(chǎn)出率，制造商生產(chǎn)需求率，制造商發(fā)貨率，批發(fā)商訂貨率，批發(fā)商發(fā)貨率，零售商訂貨率，市場銷售率，市場需求率。

2.3K產(chǎn)品供應(yīng)鏈的系統(tǒng)動力學(xué)模型

制造商產(chǎn)出率是其生產(chǎn)需求率的一階延遲，制造商的發(fā)貨率是批發(fā)商訂貨率的三階延遲，批發(fā)商的發(fā)貨率是零售商訂貨率的三階延遲。K產(chǎn)品供應(yīng)鏈的系統(tǒng)動力學(xué)模型具體如下：（1）零售商安全庫存量=12件；（2）零售商庫存調(diào)節(jié)時間=4周；（3）零售商庫存調(diào)節(jié)率=（零售商安全庫存量-零售商期末庫存量）/零售商庫存調(diào)節(jié)時間，件/周；（4）零售商期末庫存量=INTEG（DELAY3（批發(fā)商發(fā)貨率，4）-市場銷售率，12），件；（5）零售商期初庫存量=DELAY3（零售商期末庫存量, 1），件；（6）零售商累計缺貨量=INTEG（市場需求率-市場銷售率，0），件；（7）零售商訂貨率=零售商庫存調(diào)節(jié)率+市場需求率，件/周；（8）市場需求率=4+STEP（4 , 2），件/周；（9）市場銷售率=MIN（零售商期初庫存量+DELAY3（批發(fā)商發(fā)貨率，4），市場需求率），件/周；（10）批發(fā)商安全庫存量=12件；（11）批發(fā)商庫存調(diào)節(jié)時間=4周；（12）批發(fā)商庫存調(diào)節(jié)率=（批發(fā)商安全庫存量-批發(fā)商期末庫存量）/批發(fā)商庫存調(diào)節(jié)時間，件/周；（13）批發(fā)商期末庫存量=INTEG（DELAY3（制造商發(fā)貨率, 4）-批發(fā)商發(fā)貨率，12），件；（14）批發(fā)商期初庫存量=DELAY3（批發(fā)商期末庫存量, 1），件；（15）批發(fā)商累計缺貨量=INTEG （DELAY3（零售商訂貨率, 4）-批發(fā)商發(fā)貨率，0），件；（16）批發(fā)商訂貨率=DELAY3（零售商訂貨率, 4）+批發(fā)商庫存調(diào)節(jié)率，件/周；（17）批發(fā)商發(fā)貨率=MIN（批發(fā)商期初庫存量+DELAY3（制造商發(fā)貨率, 4）, DELAY3（零售商訂貨率, 4）+DELAY3（批發(fā)商累計缺貨量, 1）），件/周；（18）制造商安全庫存量=12件；（19）制造商庫存調(diào)節(jié)時間=4周；（20）制造商庫存調(diào)節(jié)率=（制造商安全庫存量-制造商期末庫存量）/制造商庫存調(diào)節(jié)時間，件/周；（21）制造商期末庫存量=INTEG（產(chǎn)出率-制造商發(fā)貨率，12），件；（22）制造商期初庫存量=DELAY3（制造商期末庫存量, 1），件；（23）制造商累計缺貨量=INTEG（DELAY3（批發(fā)商訂貨率, 4）-制造商發(fā)貨率，0），件；（24）生產(chǎn)需求率=DELAY3（批發(fā)商訂貨率, 4）+制造商庫存調(diào)節(jié)率，件/周；（25）制造商發(fā)貨率=MIN（制造商期初庫存量+產(chǎn)出率, DELAY3（批發(fā)商訂貨率, 4）+DELAY3（制造商累計缺貨量, 1）），件/周；（26）產(chǎn)出率=DELAY3（生產(chǎn)需求率, 2），件/周。

2.4Vensim系統(tǒng)動力學(xué)仿真

本仿真的時間范圍是0～300周, 步長為0.125，在Vensim PLE 5.4b版本上仿真運(yùn)行該模型，各個物流環(huán)節(jié)期末庫存量變化趨勢如圖3所示。從圖3中我們可以看出，由于K產(chǎn)品供應(yīng)鏈上各個環(huán)節(jié)各自為政，庫存調(diào)節(jié)時間較短，訂單延遲時間較長，導(dǎo)致了該條供應(yīng)鏈上牛鞭效應(yīng)顯著，這與Y公司K產(chǎn)品現(xiàn)階段供應(yīng)鏈實(shí)際情況相符合。

2.5解決方案

這里我們綜合采用縮短訂單延遲時間，延長庫存調(diào)節(jié)時間，采用供應(yīng)商管理庫存的解決辦法?，F(xiàn)階段K產(chǎn)品供應(yīng)鏈上，制造商產(chǎn)出率是其生產(chǎn)需求率的一階延遲，制造商的發(fā)貨率是批發(fā)商訂貨率的三階延遲，批發(fā)商的發(fā)貨率是零售商訂貨率的三階延遲，全都改進(jìn)為一階延遲；各物流環(huán)節(jié)庫存調(diào)整時間為4天，全都改進(jìn)為延長8天；采用供應(yīng)商管理庫存的解決方案。各個物流環(huán)節(jié)的訂貨策略改變?yōu)椋篞ij=SEij-Iij/Ti+Nj，SEij為期望庫存量，對下一環(huán)節(jié)需求率做一次指數(shù)平滑預(yù)測。變化的系統(tǒng)動力學(xué)模型為：（1）零售商期望庫存量=SMOOTHI（市場需求率, 1, 零售商平滑系數(shù)），件；（2）零售商庫存調(diào)節(jié)率=（零售商期望庫存量-零售商期末庫存量）/零售商庫存調(diào)節(jié)時間，件/周；（3）批發(fā)商期望庫存量=SMOOTHI（零售商訂貨率, 1, 批發(fā)商平滑系數(shù)），單位：件；（4）批發(fā)商庫存調(diào)節(jié)率=（批發(fā)商期望庫存量-批發(fā)商期末庫存量）/批發(fā)商庫存調(diào)節(jié)時間，件/周；（5）制造商期望庫存量=SMOOTHI（批發(fā)商訂貨率, 1, 制造商平滑系數(shù)），件。平滑系數(shù)取為0.1；（6）制造商庫存調(diào)節(jié)率=（制造商期望庫存量-制造商期末庫存量）/制造商庫存調(diào)節(jié)時間，件/周。原模型中去掉零售商、批發(fā)商、制造商安全庫存量的方程，其他不變。改進(jìn)后效果如圖4所示。

3結(jié)論與評價

本文利用系統(tǒng)動力學(xué)建模，在系統(tǒng)動力學(xué)仿真軟件Vensim PLE 5.4b版本上進(jìn)行仿真，主要得到以下結(jié)論：（1）圖3表明：由于訂單延遲時間較長，庫存調(diào)節(jié)時間較短，采用傳統(tǒng)庫存管理辦法，Y公司K產(chǎn)品供應(yīng)鏈上每個物流環(huán)節(jié)各自為證，牛鞭效應(yīng)顯著；（2）從圖4中我們可以看出，綜合利用三種解決方案，可以得到顯著的效果，系統(tǒng)在60周左右處于完全平衡狀態(tài)，與單靠縮短延遲時間相比，提前了30周的時間。本文針對Y公司K產(chǎn)品的供應(yīng)鏈上牛鞭效應(yīng)顯著的情況，進(jìn)行了定性分析，并給出可行的綜合解決方案，然而延遲時間縮短為多少，庫存調(diào)節(jié)時間延長為多少，供應(yīng)商管理庫存有許多方法，對下游物流環(huán)節(jié)需求的預(yù)測，對庫存的調(diào)節(jié)率都有相當(dāng)多的理論與方法，怎樣保證改進(jìn)方案既經(jīng)濟(jì)又可行，還需具體問題具體分析，這些都是本文值得進(jìn)一步研究的地方。

參考文獻(xiàn)：

[1]彭揚(yáng), 伍蓓. 物流系統(tǒng)優(yōu)化與仿真[M]. 北京: 中國物質(zhì)出版社, 2007.

[2]劉媛媛, 王海燕. 供應(yīng)鏈節(jié)點(diǎn)變化對牛鞭效應(yīng)影響的系統(tǒng)動力學(xué)仿真研究[J]. 物流技術(shù), 2006,25(6):50-53.

[3]李穩(wěn)安, 趙林度. 牛鞭效應(yīng)的系統(tǒng)動力學(xué)分析[J]. 東南大學(xué)學(xué)報: 哲學(xué)社會科學(xué)版, 2002,33(10):96-98.

[4]王磊, 陳競先, 唐志. 供應(yīng)鏈中牛鞭效應(yīng)的模型與分析[J]. 物流技術(shù), 2004,10(1):42-46.

[5]楊詩華. 供應(yīng)鏈管理中牛鞭效應(yīng)機(jī)理分析與仿真[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2006.

[6]苗興東. 系統(tǒng)動力學(xué)在供應(yīng)鏈建模中的應(yīng)用[D]. 成都: 西南交通大學(xué), 2006.