整合供應(yīng)商選擇的多層級選址-庫存模型

稅文兵，何保紅，何 民

昆明理工大學(xué) 交通工程學(xué)院，昆明 650500

整合供應(yīng)商選擇的多層級選址-庫存模型

稅文兵，何保紅，何 民

昆明理工大學(xué) 交通工程學(xué)院，昆明 650500

1 引言

傳統(tǒng)的“先選址，再庫存”的順序決策方法，從總成本的角度來看，并不是最優(yōu)的[1]。近年來，將選址分配決策問題和庫存控制決策問題進(jìn)行整合，實(shí)現(xiàn)選址成本、運(yùn)輸成本和庫存成本的平衡，從而進(jìn)一步降低總成本，成為了物流系統(tǒng)優(yōu)化的發(fā)展方向之一?；谶x址分配決策和庫存控制決策的相互影響關(guān)系，從系統(tǒng)整體最優(yōu)的角度，建立同時(shí)考慮選址成本、運(yùn)輸成本和庫存成本的選址-庫存模型，對設(shè)施的位置、數(shù)量、服務(wù)范圍和庫存控制參數(shù)同時(shí)進(jìn)行優(yōu)化決策，在文獻(xiàn)中被稱為選址-庫存問題（Location-Inventory Problems，LIP）。

早在1958年，Baumol就提出應(yīng)該在選址模型的目標(biāo)函數(shù)中包含庫存成本[2]。文獻(xiàn)[1]建立了確定環(huán)境下基于經(jīng)濟(jì)訂貨批量（EOQ）的連續(xù)選址-庫存模型，并采用微分方法求解，但是該模型沒有考慮設(shè)施的建立成本。文獻(xiàn)[3]研究了需求和提前期隨機(jī)變化下的選址-庫存問題，證明了當(dāng)配送中心采用連續(xù)(S-1，S)補(bǔ)貨策略時(shí)，系統(tǒng)的總安全庫存與總配送中心數(shù)存在線性關(guān)系，并利用這種關(guān)系建立了線性的選址-庫存模型。然而，由于(S-1，S)庫存策略一般用于需求較少、提前期短的情況。對于需求量大的重要產(chǎn)品，往往采用(Q，R)庫存策略。因此，(Q，R)庫存策略下的選址-庫存問題成為了重要的研究方向。

文獻(xiàn)[4]針對由一個(gè)供應(yīng)商、多個(gè)分銷中心和多個(gè)零售商的供應(yīng)鏈，基于Axsater[5]的庫存成本計(jì)算方法，建立了第一個(gè)(Q，R)庫存策略下的選址-庫存模型。該模型較好地體現(xiàn)了選址分配決策和庫存控制決策之間的相互影響關(guān)系，從而成為了這一研究方向的基礎(chǔ)模型，得到了很多學(xué)者的進(jìn)一步研究。如文獻(xiàn)[6]建立了基于場景規(guī)劃的模型，文獻(xiàn)[7]建立了考慮設(shè)施能力限制的模型，文獻(xiàn)[8]建立了考慮工廠選址的三層級選址模型等。這些研究從不同的角度完善了Shen的模型，然而，它們都只考慮了一個(gè)庫存點(diǎn)，沒有考慮多級庫存的情況。

除(Q，R)庫存策略下的LIP模型，其他庫存策略下的模型也得到了研究。文獻(xiàn)[9]針對連續(xù)檢查(Q，R)策略下的選址-庫存模型，會(huì)使每一個(gè)配送中心的訂貨間隔期都不相同，從而導(dǎo)致供應(yīng)商在供貨時(shí)不能夠協(xié)調(diào)供貨，增加車隊(duì)的調(diào)度使用成本的問題，建立了配送中心采用定期檢查(T，S)策略下的選址-庫存模型。文獻(xiàn)[10]針對單供應(yīng)商、配送中心和顧客構(gòu)成的供應(yīng)鏈，建立了配送中心采用(s，S)策略、顧客采用(R，S，Q)策略的多級選址-庫存模型，然而，模型沒有考慮不同庫存點(diǎn)之間的相互影響關(guān)系。文獻(xiàn)[11]基于Roundy[12]的多級庫存控制方法，建立了分銷中心和顧客都有庫存的選址-庫存模型，但仍然僅限于一個(gè)供應(yīng)商的三層分銷網(wǎng)絡(luò)。文獻(xiàn)[13]基于EOQ模型假設(shè)和單點(diǎn)庫存控制策略，建立了允許供應(yīng)商可以直接向顧客供貨的多級選址-庫存模型，然而，該模型沒有明確考慮空間距離對配送中心選址的影響。

上述研究有三個(gè)特點(diǎn)：（1）從供應(yīng)鏈階段來看，大部分都只針對由工廠、配送中心和顧客構(gòu)成的三級分銷網(wǎng)絡(luò)，沒有考慮零部件供應(yīng)商的選擇；（2）從選址對象來看，都只對配送中心進(jìn)行選址分配決策，沒有對工廠進(jìn)行選址決策；（3）從所建立的模型特征來看，都是將選址分配模型和庫存控制模型融合為一個(gè)模型。前兩個(gè)特點(diǎn)決定了所建立的模型不能用于由供應(yīng)商、工廠、配送中心和顧客組成的完整供應(yīng)鏈網(wǎng)絡(luò)設(shè)計(jì)。最后一個(gè)特點(diǎn)使得所建立的模型呈現(xiàn)出高度的非線性，采用傳統(tǒng)的精確算法難以在較短時(shí)間內(nèi)求解較大規(guī)模的問題。

在一個(gè)完整的供應(yīng)鏈中，供應(yīng)商的選擇、工廠和配送中心的選址、供應(yīng)商與工廠以及工廠與配送中心之間的運(yùn)輸量以及工廠和配送中心的庫存量都存在相互影響的關(guān)系。在進(jìn)行這類供應(yīng)鏈網(wǎng)絡(luò)設(shè)計(jì)和優(yōu)化時(shí)，如果忽視了其中某個(gè)影響因素，所得到的設(shè)計(jì)方案從系統(tǒng)的角度看并不是最優(yōu)的。本文針對由供應(yīng)商、工廠、配送中心和顧客構(gòu)成的四層級供應(yīng)鏈網(wǎng)絡(luò)，構(gòu)建同時(shí)對供應(yīng)商的選擇、工廠和配送中心的選址和庫存控制參數(shù)進(jìn)行決策的LIP模型。其中，利用梯級庫存概念，處理供應(yīng)鏈中上下游庫存點(diǎn)之間的相互影響關(guān)系，建立工廠和配送中心的庫存成本計(jì)算模型。考慮到所建模型的復(fù)雜性，設(shè)計(jì)了基于部分編碼的粒子群優(yōu)化算法求解，并通過多個(gè)算例驗(yàn)證模型的有效性和算法的可行性。

2 模型建立

本文所研究的供應(yīng)鏈網(wǎng)絡(luò)由供應(yīng)商、工廠、配送中心和顧客組成，供應(yīng)商向工廠提供原材料，工廠根據(jù)配送中心的訂單進(jìn)行生產(chǎn)，產(chǎn)品生產(chǎn)出來之后通過配送中心運(yùn)往顧客。配送中心通過成品庫存滿足若干顧客的確定需求，工廠分別通過原材料庫存和成品庫存滿足生產(chǎn)和配送中心的需要。在已知顧客的需求、單位運(yùn)輸成本、候選的供應(yīng)商、工廠和配送中心位置等參數(shù)下，對供應(yīng)商的選擇、工廠和配送中心的選址、節(jié)點(diǎn)之間的運(yùn)輸量，以及設(shè)施的庫存量進(jìn)行決策。

2.1 模型參數(shù)和變量

基本參數(shù)：I為顧客的集合，某一個(gè)顧客用i表示；J為配送中心的集合，某一個(gè)配送中心用 j表示；P為工廠的集合，某一個(gè)工廠用 p表示；S為供應(yīng)商的集合，某一個(gè)配送中心用s表示。di為顧客的需求（單位/年）；DDj為配送中心的需求（單位/年）；DPp為工廠成品倉庫的需求（單位/年）；DMp為工廠原材料倉庫的需求（單位/年）；FDj為配送中心的打開成本（元/年）；FPp為工廠的打開成本（元/年）；FSs為供應(yīng)商的選擇成本（元/年）；ODj為配送中心訂貨成本（元/次）；OPp為工廠成品訂貨成本（元/次）；OMp為工廠原材料訂貨成本（元/次）；HDj為配送中心單位產(chǎn)品每年的庫存持有成本（元/單位）；HPp為工廠單位產(chǎn)品每年的庫存持有成本（元/單位）；HMp為工廠單位原材料每年的庫存持有成本（元/單位）；EDj為配送中心單位產(chǎn)品每年的梯級庫存持有成本（元/單位）；EPp為工廠單位產(chǎn)品每年的梯級庫存持有成本（元/單位）；EMp為工廠單位原材每年的料梯級庫存持有成本（元/單位）；MPp為工廠單位產(chǎn)品生產(chǎn)成本（元/單位）；MHp為工廠單位原材料處理成本（元/單位）；PHp為工廠單位成品處理成本（元/單位）；DHj為配送中心單位成品處理成本（元/單位）；PSs為供應(yīng)商的單位原材料售價(jià)（元/單位）；TDji為從配送中心到顧客的單位運(yùn)輸成本（元/單位）；TPpj為從工廠到配送中心的單位運(yùn)輸成本（元/單位）；TSsp為從供應(yīng)商到工廠的單位運(yùn)輸成本（元/單位）；BOM為成品和原材料之間的數(shù)量關(guān)系；B為大的正數(shù)。

決策變量：XDj為如果配送中心打開取1，否則取0；XPp為如果工廠打開取1，否則取0；XSs為如果供應(yīng)商被選用取1，否則取0；YDji為如果配送中心和顧客之間建立服務(wù)關(guān)系取1，否則取0；YPpj為如果工廠和配送中心之間建立服務(wù)關(guān)系取1，否則取0；YSsp為如果供應(yīng)商和工廠之間建立服務(wù)關(guān)系取1，否則取0；QPpj為工廠到配送中心的運(yùn)輸量（單位/年）；QSsp為供應(yīng)商到工廠的運(yùn)輸量（單位/年）；OQDj為配送中心訂貨批量（單位/次）；OQMp為工廠原材料訂貨批量（單位/次）；OQPp為工廠成品訂貨批量（單位/次）。

2.2 模型數(shù)學(xué)描述

整合供應(yīng)商選擇的多層級LIP模型，以設(shè)施的打開成本、原材料和成品庫存成本、原材料采購和運(yùn)輸成本、成品的生產(chǎn)和運(yùn)輸成本之和最小為目標(biāo)，同時(shí)對供應(yīng)商的選擇、工廠和配送中心的選址、配送中心和顧客之間的分配關(guān)系、層級之間的運(yùn)輸量以及工廠和配送中心的訂貨批量進(jìn)行決策，模型如下：

其中，式（1）是模型的目標(biāo)函數(shù)，前三項(xiàng)分別是供應(yīng)商的選擇成本、工廠和配送中心的打開成本；第四項(xiàng)是原材料的采購和運(yùn)輸成本，第五項(xiàng)是產(chǎn)品的生產(chǎn)和從工廠到供應(yīng)商的運(yùn)輸成本；第六項(xiàng)是從配送中心到顧客的運(yùn)輸成本；第七項(xiàng)是配送中心的庫存持有成本；第八和第九項(xiàng)是工廠的庫存持有成本。式（2）表示只有當(dāng)配送中心打開后，才能為顧客提供服務(wù)；式（3）和（4）表示只有當(dāng)供應(yīng)商、工廠和配送中心都打開后，它們之間才能建立服務(wù)關(guān)系；式（5）表示一個(gè)配送中心就能完全滿足一個(gè)顧客的需求；式（6）和（7）分別是進(jìn)出配送中心和工廠的產(chǎn)品數(shù)量相等；式（8）和（9）表示供應(yīng)商、工廠和配送中心之間運(yùn)輸量的取值限制；式（10）是配送中心的梯級庫存持有成本；式（11）和（12）分別是工廠成品和原材料的梯級庫存持有成本；式（13）和式（14）是廠和配送中心訂貨批量關(guān)系限制；式（15）～式（23）是決策變量的取值限制。

3 算法設(shè)計(jì)

所建立的四層級LIP模型屬于混合整數(shù)非線性規(guī)劃模型，同時(shí)也是典型的NP困難問題，傳統(tǒng)的精確算法難以在有限時(shí)間內(nèi)求解較大規(guī)模的問題，常用啟發(fā)式算法獲得在可接受時(shí)間花費(fèi)下的可行解[14]。啟發(fā)式算法總體上可以分為個(gè)體搜索算法和群體搜索算法。個(gè)體搜索算法從單個(gè)個(gè)體開始，逐步往最優(yōu)解靠攏，如貪心算法、松弛算法等。群體搜索算法從解的集合開始，采用并行搜索機(jī)制，接近最優(yōu)解，如遺傳算法、粒子群優(yōu)化算法（Particle Swarm Optimization，PSO）等。實(shí)踐證明，采用并行機(jī)制的群體搜索算法在效率上優(yōu)于個(gè)體搜索算法。相對于遺傳算法，粒子群優(yōu)化算法沒有交叉和變異操作，只是通過內(nèi)部速度進(jìn)行更新，因此原理更簡單、參數(shù)更少、實(shí)現(xiàn)更容易。因此，PSO被廣泛應(yīng)用于各個(gè)領(lǐng)域，是求解非線性連續(xù)優(yōu)化、混合整數(shù)非線性優(yōu)化、組合優(yōu)化問題的有力工具[15]。所以，本文采用PSO算法求解所建立的模型。同時(shí)，為了對PSO算法的求解能力進(jìn)行評價(jià)，將PSO算法的求解結(jié)果與基于分支定界算法的Lingo軟件求解結(jié)果進(jìn)行比較。

3.1 粒子群初始化

粒子群初始化包括粒子的編碼和粒子群的生成。上述LIP模型的決策變量包含0-1整數(shù)變量和連續(xù)非負(fù)變量兩大類，經(jīng)過試驗(yàn)發(fā)現(xiàn)，由于模型過于復(fù)雜，無論是全部采用二進(jìn)制編碼，還是采用混合編碼，都很難得到滿足所有約束條件的可行解。因此，本文采用部分編碼的策略，僅對其中的選址變量進(jìn)行編碼（二進(jìn)制編碼）。一旦選址變量確定后，運(yùn)輸量通過求解最小成本網(wǎng)絡(luò)流問題得到。由于運(yùn)輸量的有無反映了節(jié)點(diǎn)之間的網(wǎng)絡(luò)拓?fù)潢P(guān)系，因此不需要對其進(jìn)行專門的編碼。訂貨批量通過求解EOQ模型得到，進(jìn)而得到總的運(yùn)輸成本和總的庫存成本。

在上述編碼規(guī)則下，每一個(gè)粒子的長度為候選供應(yīng)商、工廠和配送中心數(shù)量之和，粒子群采用隨機(jī)的方法生成。

3.2 適應(yīng)度函數(shù)的確定

在選址變量確定后，可以得到總的選址成本LocaC。通過求解如下的最小成本網(wǎng)絡(luò)流問題得到總的運(yùn)輸成本TraficC。

模型中，QDj是配送中心到顧客的運(yùn)輸量，其他變量含義同前。采用最小費(fèi)用增廣路徑算法求解上述網(wǎng)絡(luò)流問題[16]，從而得到層級之間的運(yùn)輸量和TraficC。在確定選址變量和運(yùn)輸量的結(jié)果后，通過EOQ模型計(jì)算最初訂貨批量，再結(jié)合約束條件進(jìn)行調(diào)整，進(jìn)而得到滿足約束的訂貨批量和總的庫存成本InvetC。

通過上述分布求解的方式，有效地消除了所有的約束。因此，適應(yīng)度函數(shù)為：

3.3 速度和位置的更新

適應(yīng)度函數(shù)確定后，計(jì)算各個(gè)粒子的適應(yīng)度值，并對個(gè)體極值和群體極值更新。然后，對粒子速度和位置進(jìn)行更新。在離散空間下，速度和位置的更新公式如下[17]：

其中，ω稱作慣性權(quán)重，c1、c2為常數(shù)，稱作學(xué)習(xí)因子（也叫加速系數(shù)），r1、r2是兩個(gè)0到1之間的隨機(jī)數(shù)。

4 算例仿真

粒子群算法程序用Matlab 7.10中的M語言編寫，測試用的計(jì)算機(jī)處理器為Intel?CoreTMi5-2410 CPU，單核頻率為2.30 GHz，2.0 GB的內(nèi)存，操作系統(tǒng)為Windows 7家庭普通版。為了對PSO進(jìn)行評估，同時(shí)用優(yōu)化求解軟件Lingo 9.0求解每一個(gè)算例。

4.1 算例數(shù)據(jù)

根據(jù)問題規(guī)模隨機(jī)生成20個(gè)算例，每一個(gè)算例有不同的節(jié)點(diǎn)數(shù)量，總的節(jié)點(diǎn)數(shù)量從50變化到170個(gè)。例如，算例“5-10-15-20”表示該算例中候選的供應(yīng)商有5個(gè)，工廠有10個(gè)，配送中心有15個(gè)，顧客有20個(gè)。

各算例的數(shù)據(jù)產(chǎn)生規(guī)則如下：供應(yīng)商的固定選擇成本，工廠和配送中心的打開成本和服從[5 000，15 000]之間的均勻分布；供應(yīng)商到工廠的單位運(yùn)輸成本服從[5，15]之間的均勻分布，工廠到配送中心的單位運(yùn)輸成本服從[20，50]之間的均勻分布，配送中心到顧客的單位運(yùn)輸成本服從[100，200]之間的均勻分布，供應(yīng)商的原材料采購成本服從[200，300]之間的均勻分布，工廠的單位生產(chǎn)成本服從[1 000，1 500]之間的均勻分布；工廠成品一次訂貨成本服從[1 500，2 000]之間的均勻分布，工廠原材料一次訂貨成本服從[2 000，2 500]之間的均勻分布，工廠成品的單位庫存持有成本服從[150，200]之間的均勻分布，工廠單位原材料庫存持有成本服從[50，100]之間的均勻分布，配送中心一次訂貨成本服從[3 000，4 000]之間的均勻分布，配送中心的單位庫存持有成本服從[500，800]之間的均勻分布；顧客的需求服從[1 000，3 000]之間的均勻分布。對于所有的算例，BOM都取2。

另外，設(shè)定3 600 s的時(shí)間限制作為Lingo算法的終止條件，粒子群算法的終止條件為連續(xù)10次迭代，適應(yīng)度函數(shù)沒有明顯的改進(jìn)。

4.2 算法比較

上述20個(gè)算例分別用Lingo 9.0和PSO算法求解，結(jié)果如表1所示。第一列“測試算例”表明了每一個(gè)算例中候選供應(yīng)商、工廠、配送中心和顧客的數(shù)量；第二列和第三列是Lingo 9.0求解后得到的結(jié)果，分別是目標(biāo)函數(shù)最小值（“cost”單位是元）和運(yùn)算所用時(shí)間（“T”單位是s）；第四列和第五列是PSO算法求解后的結(jié)果，同樣包含目標(biāo)函數(shù)最小值和運(yùn)算所用時(shí)間。從表1可以看出：有55%的算例lingo9.0不能在3 600 s內(nèi)得到求解結(jié)果，尤其是當(dāng)總節(jié)點(diǎn)數(shù)超過110個(gè)以后，全部都不能得到最小成本，而PSO算法對于所有的算例在較短時(shí)間內(nèi)都能得到最小成本。這表明PSO的求解能力優(yōu)于優(yōu)化軟件。

表1 算例求解結(jié)果

從兩種算法都能得到最小成本的算例來看（圖1），Lingo 9.0所得到的解都要大于PSO。這說明基于分支定界算法的Lingo 9.0在有限的時(shí)間內(nèi)只能獲得局部最優(yōu)解，而且這種局部最優(yōu)解比PSO的解質(zhì)量差。

圖1 PSO和Lingo的求解結(jié)果比較

從兩種算法的求解時(shí)間來看（圖2），在節(jié)點(diǎn)數(shù)小于70個(gè)的時(shí)候，Lingo 9.0的獲得最小成本的時(shí)間要少于PSO；當(dāng)節(jié)點(diǎn)數(shù)超過70個(gè)時(shí)，情況剛好相反。而且，隨著節(jié)點(diǎn)數(shù)的增加，Lingo 9.0的求解時(shí)間增加速度非常快，超過了PSO求解時(shí)間的增加速度。

圖2 PSO和Lingo的求解時(shí)間比較

5 結(jié)論

供應(yīng)鏈網(wǎng)絡(luò)設(shè)計(jì)是供應(yīng)鏈管理的戰(zhàn)略決策問題，對供應(yīng)鏈的績效有重要的影響。合理設(shè)計(jì)供應(yīng)鏈網(wǎng)絡(luò)，要求建立選址和庫存的集成決策模型。本文針對由供應(yīng)商、工廠、配送中心和顧客構(gòu)成的供應(yīng)鏈，考慮工廠和配送中心庫存決策的相互影響，運(yùn)用梯級庫存的概念，建立了整合供應(yīng)商選擇的多層級的LIP模型。針對模型的特點(diǎn)，設(shè)計(jì)了PSO算法；通過對20個(gè)算例的求解，可以看出本文所設(shè)計(jì)的PSO算法無論是在得到的最小總成本，還是運(yùn)算速度上都明顯優(yōu)于數(shù)學(xué)規(guī)劃求解軟件Lingo，尤其是當(dāng)總節(jié)點(diǎn)數(shù)較大時(shí)。算例也表明，本文所建立的模型是有效的，設(shè)計(jì)的算法是可行的，可以用于多層級的一體化供應(yīng)鏈網(wǎng)絡(luò)設(shè)計(jì)。在未來的研究中，可以考慮需求的不確定性和更為復(fù)雜的多級庫存控制模型。

[1]Mccann P.The logistics cost location-production problem[J]. Journal of Regional Science，1993，33（4）：503-516.

[2]Baumol W J，Wolfe P.A warehouse-location problem[J].Operations Research，1958，6（2）：252-263.

[3]Nozick L K，Turnquist M A.Integrating inventory impacts into a fixed-charge model for locating distribution centers[J]. Transportation Research Part E，1998，34（3）：173-186.

[4]Shen Z J M.Efficient algorithms for various supply chain problems[D].Evanston，IL：Northwestern University，2000.

[5]Axsater S.Using the deterministic EOQ formula in stochastic inventory control[J].Management Science，1996，42（6）：830-834.

[6]Snyder L V，Daskin M S，Teo C P.The stochastic location model with risk pooling[J].European Journal of Operational Research，2007，3（179）：1221-1238.

[7]Miranda P A，Garrido R A.Incorporating inventory control decisions into a strategic distribution network design model with stochastic demand[J].Transportation Research Part E，2004，40：183-207.

[8]Park S，Lee T T，Sung C S.A three-level supply chain network design modelwith risk-pooling and lead times[J]. Transportation Research Part E，2010，46：563-581.

[9]Berman O，Krass D，Tajbakhsh M M.A coordinated locationinventory model[J].European Journal of Operational Research，2012，217：500-508.

[10]譚凌，高峻峻，王迎軍.基于庫存成本優(yōu)化的配送中心選址問題研究[J].系統(tǒng)工程學(xué)報(bào)，2004，19（l）：59-66.

[11]Romeijn H E，Shu J，Teo C P.Designing two-echelon supply networks[J].European Journal of Operational Research，2007，178：449-462.

[12]Roundy R O.98%effective integer-ratio lot-sizing for one warehouse multi-retailer systems[J].ManagementScience，1985，31：1416-1430.

[13]Tancrez J S，Lange J C，Semal P.A location-inventory model for large three-level supply chains[J].Transportation Research Part E，2012，48：485-502.

[14]邢文訓(xùn)，謝金星.現(xiàn)代優(yōu)化計(jì)算方法[M].北京：清華大學(xué)出版社，2003：13-23.

[15]楊維，李歧強(qiáng).粒子群優(yōu)化算法綜述[J].中國工程科學(xué)，2004，6（5）：87-94.

[16]Ahuja R，Magnanti T，Orlin J.Network flows：theory，algorithms，and applications[M].New York：Prentice Hall，1993.

[17]Kennedy J，Eberhart R C.Discrete binary version of the particle swarm algorithms[C]//Proceedings of the IEEE InternationalConferenceon Systems，Man，Cybernetics. Orlando：IEEE Press，1997：4104-4108.

SHUI Wenbing,HE Baohong,HE Min

Faculty of Transportation Engineering,Kunming University of Science&Technology,Kunming 650500,China

To deal with the integrated decisions of location and inventory for the supply chain network design with multi-echelons, a new location-inventory model based on echelon stock policy is developed in which the selection of suppliers is considered. The objective of the model is to determine the number and locations of suppliers,plants and Distribution Centers（DCs）,the assignment of each customer to a DC,the transportation quantity between suppliers,plants and distribution centers,and the optimal order quantity in the opening plants and DCs,which minimizes the system-wide location,production,purchase,transportation,and inventory costs.The particle swarm optimization algorithm based on partial encoding is designed to solve the proposed model.The results of 20 numerical examples of different scale show that the model is effective and can be used to design some four echelons supply chain network,and the PSO algorithm designed performs better than mathematic optimization software of Lingo 9.0 whatever the precision or the speed,especially when there are many nodes in the supply chain.

location-inventory model;selection of suppliers;particle swarm optimization algorithm;partial encoding

針對多級供應(yīng)鏈網(wǎng)絡(luò)設(shè)計(jì)中選址和庫存一體化決策問題，基于梯級庫存策略，建立了整合供應(yīng)商選擇的多層級選址-庫存模型。模型以網(wǎng)絡(luò)中供應(yīng)商的選擇成本、工廠和配送中心的打開成本、層級之間的運(yùn)輸成本、庫存成本、采購成本和生產(chǎn)成本之和最小為目標(biāo)，同時(shí)對供應(yīng)商的選擇、工廠和配送中心的選址、配送中心對顧客的分配、層級之間的運(yùn)輸量、工廠和配送中心的訂貨批量進(jìn)行決策。為了求解所建立的模型，設(shè)計(jì)了基于部分編碼的粒子群優(yōu)化算法。20個(gè)不同規(guī)模的算例測試表明：所建立的模型是有效的，能用于多層級供應(yīng)鏈網(wǎng)絡(luò)的設(shè)計(jì)；所設(shè)計(jì)的算法無論是在求解精度，還是在運(yùn)算速度上都明顯優(yōu)于數(shù)學(xué)優(yōu)化軟件Lingo 9.0，尤其是當(dāng)供應(yīng)鏈網(wǎng)絡(luò)中總節(jié)點(diǎn)數(shù)較大時(shí)。

選址-庫存模型；供應(yīng)商選擇；粒子群優(yōu)化算法；部分編碼

A

F252

10.3778/j.issn.1002-8331.1305-0245

SHUI Wenbing,HE Baohong,HE Min.Multi-echelons location-inventory model with selection of suppliers.Computer Engineering and Applications,2013,49（23）：1-5.

國家自然科學(xué)基金項(xiàng)目（No.51068013）；云南省科技廳應(yīng)用基礎(chǔ)研究計(jì)劃項(xiàng)目（No.2011FZ054）。

稅文兵（1982—），男，博士，講師，研究領(lǐng)域?yàn)槲锪飨到y(tǒng)設(shè)計(jì)與優(yōu)化；何保紅（1973—），女，博士，副教授，研究領(lǐng)域?yàn)榻煌ㄟ\(yùn)輸規(guī)劃與管理；何民（1972—），男，博士，教授，研究領(lǐng)域?yàn)榻煌ㄟ\(yùn)輸規(guī)劃與管理。E-mail:swbzzq@126.com

2013-05-20

2013-07-18

1002-8331（2013）23-0001-05

CNKI出版日期：2013-07-29 http://www.cnki.net/kcms/detail/11.2127.TP.20130729.1105.005.html