朝陽重型機(jī)器有限公司

朝陽重型機(jī)器有限公司

——配送中心選址目錄重心法選址模型假設(shè)條件模型結(jié)構(gòu)求解思路實證分析企業(yè)的選取與數(shù)據(jù)的調(diào)查重心法的實證模型及其數(shù)據(jù)處理假設(shè)條件重心法的應(yīng)用對象是OD(Origin-Destination)流量的交通網(wǎng)絡(luò)問題,即起點到終點的運輸流量構(gòu)成的物流網(wǎng)絡(luò)規(guī)劃問題。重心法進(jìn)行決策的依據(jù)是產(chǎn)品運輸成本的最小化,這樣就涉及到如下幾個假設(shè)前提條件:（1）運輸費用只與配送中心和配送點的直線距離有關(guān)，不考慮城市交通狀況；（2）選擇配送中心時，不考慮配送中心所處地理位置的地產(chǎn)價格；（3）運輸費率與運輸距離和運輸量呈線性關(guān)系；（4）決策各點的需求量不是地理位置上所實際發(fā)生的需求量,而是一個匯總量,這個量聚集了分散在一定區(qū)域內(nèi)眾多的需求量；（5）各配送點的需求量已知；（6）可以估計各個備選配送中心的固定費用（包括基本建設(shè)費和固定經(jīng)營費）；（7）可以估計經(jīng)營管理產(chǎn)生的可變費用，并在總費用中加以考慮。模型結(jié)構(gòu)設(shè)有n個配送點，他們各自的坐標(biāo)是（xi，yi）（i=1,2,3,…,n）配送中心的坐標(biāo)是（x0，y0）。運輸費用為E；總費用為C則有：E=aiwidi(2.1)minC(x)=1EIi+2VIi+3CIi(2.2)式中:ai表示從配送中心到配送點i每單位運量、單位運距的運輸費用；wi表示配送中心到配送點i的運輸量，也表示第i個配送點的需求量；di表示從配送中心到配送點i的直線距離；Ii表示由重心法得到的各個備選地址；Wi表示各個配送點的需求量之和；EIi表示備選地址Ii總的運輸費用；VIi表示各備選地址Ii總的可變費用；CIi表示各備選地址Ii的固定費用；表示權(quán)系數(shù)(可以根據(jù)決策者的需求來定）且，其中∈（0，1）求解思路采用迭代法計算出12個重心點和重心點的運輸成本上式中：di=[(x0-xi)2+(y0-yi)2]1/2(2.3)采用微分法，將式（2.3）代入（2.1）中，為了求出使E最小的x0,y0值，對得到的公式求偏導(dǎo)，令=aiwi(x0-xi)/di=0,=aiwi(y0-yi)/di=0(2.4)由式（2.4）可以分別求得最為合適的x0和y0，即X0=,y0=(2.5)方程式(2.5)的右邊還含有未知數(shù)(x0,y0)，如果從兩個方程式的右邊完全消除x0和y0,計算將變得很復(fù)雜，計算量也很大.因此，可以采用迭代的方法進(jìn)行計算，通過迭代，得到各個備選的配送中心Ii。用迭代方法計算的方法如下：（1）以所有需求點的重心坐標(biāo)作為配送中心的初始位置坐標(biāo)(,)；（2）利用方程式（2.1）和（2.3）計算與(,)相應(yīng)的總的運輸費用E0；(3）把(,)分別代入方程式（2.3）和（2.5）中，計算配送中心的改善地點(,)；這樣反復(fù)計算下去，直到計算出12個重心點。（4）利用方程式（2.1）和（2.3）計算各個地點相對應(yīng)的總的運輸費用E；采用多元線性回歸對總成本目標(biāo)函數(shù)的系數(shù)進(jìn)行求解設(shè)y為因變量，,為自變量，并且y=C(x),=EIi,=VIi,=CIi,則多元線性回歸模型為：=(2.6)

設(shè)分別作為參數(shù)的估計量，得樣本回歸方程為：=（i=1,2…,n）(2.7)用Excel輔助計算可得到3個待估參數(shù)的估計值。采用迭代法對優(yōu)化好的模型進(jìn)行求解用迭代方法計算的方法如下：（1）以所有需求點的重心坐標(biāo)作為配送中心的初始位置坐標(biāo)(,)；（2）利用方程式（2.1）和（2.3）計算與(,)相應(yīng)的總的運輸費用E0；（3）把(,)分別代入方程式（2.3）和（2.5）中，計算配送中心的改善地點(,)；（4）利用方程式（2.1）和（2.3）計算相對應(yīng)的總的運輸費用E1；（5）把E1和E0進(jìn)行比較，如果E1＜E0則返回（2.3）的計算，再把(,)代入方程式（2.3）和（2.5）中，計算配送中心的再改善地點。如果則說明是最優(yōu)解。這樣反復(fù)計算下去，直至求出最優(yōu)解為止。

根據(jù)上面解的情況，把求出的最優(yōu)解之前的次優(yōu)解、、以及最優(yōu)解所對應(yīng)的位置作為配送中心的備選地址，記為Ii（i=0，1，…，K）。且EIi=E=aiwidi；的值為上一節(jié)所求的值。然后，將所需要的數(shù)值代入（2.2）式直接計算即可，最小的C(x)所對應(yīng)的Ii即為最優(yōu)解。企業(yè)的選取與數(shù)據(jù)的調(diào)查朝陽重型機(jī)器有限公司是在原朝重（集團(tuán)）有限責(zé)任公司、朝陽重型機(jī)器有限責(zé)任公司、朝陽重型機(jī)器廠等三家企業(yè)改制后組成的一個全新的公司。是中國建材機(jī)械行業(yè)大型骨干企業(yè)。裝備實力、產(chǎn)品銷售、創(chuàng)新能力居中國建材機(jī)械行業(yè)領(lǐng)先地位。朝重有進(jìn)出口自營權(quán)。是ISO9001質(zhì)量體系認(rèn)證合格單位。多年來，朝重先后榮獲“國家質(zhì)量一級合格單位”、“國家質(zhì)量管理獎”、“國家節(jié)能銀牌獎”、“中國環(huán)保產(chǎn)業(yè)百強企業(yè)第一名”、“中國企業(yè)最佳信譽和中國企業(yè)最佳形象AAA級單位”等榮譽稱號。朝陽重型機(jī)器有限公司主要以生產(chǎn)、研制、開發(fā)“朝重牌”建材機(jī)械產(chǎn)品為主，年生產(chǎn)能力3萬余噸。朝陽重型機(jī)器有限公司具備提供300T∕D——4000T∕D大中型水泥廠成套裝備的設(shè)計開發(fā)、生產(chǎn)制造、質(zhì)量檢驗、吊裝運輸、安裝調(diào)試的能力。同時，還提供環(huán)保設(shè)備，墻體材料成套設(shè)備，礦山、冶金、化工、壓力容器、煤炭、糧食行業(yè)的通用、專用設(shè)備以及公路碎石生產(chǎn)線主機(jī)設(shè)備等。朝陽重型機(jī)器有限公司的供應(yīng)商遍布全國各地，其供貨時間和數(shù)量相對比較隨機(jī)，即朝陽重型機(jī)器有限公司發(fā)出訂貨通知就供貨，這樣會使得朝陽重型機(jī)器有限公司方面因需要接受各地的零件而不得不建造較大的儲存空間，而接受到的零件并不會一次馬上消耗掉，因而會造成因儲存而形成的浪費。并因為各地供貨都是小批量的，因而無法形成規(guī)模效應(yīng)，這就使得朝陽重型機(jī)器有限公司在運輸方面也需要大量的投資。在這種情況下，選擇一個配送中心作為自己供貨的暫存區(qū)就顯得尤為重要。由于朝陽重型機(jī)器有限公司供應(yīng)商以長三角地區(qū)的居多，所以配送中心的選擇以長三角地區(qū)為主。一般情況下，配送中心擔(dān)任原料的收集和成品的銷售兩個任務(wù)，但在這次選址中，單考慮原料收集任務(wù)。重心法的數(shù)據(jù)處理實證模型所需數(shù)據(jù)

本課題的數(shù)據(jù)主要是通過朝陽重型機(jī)器有限公司的內(nèi)部調(diào)查取得企業(yè)內(nèi)部生產(chǎn)數(shù)據(jù)，再對數(shù)據(jù)進(jìn)行篩選加工。主要選取該公司長三角地區(qū)的供應(yīng)商的運輸重量和單位運費，備選配送中心的固定費用和總的可變費用等數(shù)據(jù)來進(jìn)行分析。供應(yīng)商坐標(biāo)整理

根據(jù)朝陽重型機(jī)器有限公司提供的2009年的數(shù)據(jù)和在中國地圖上建立直角坐標(biāo)系，統(tǒng)計出各個供應(yīng)商的坐標(biāo)，得出表4－1。表4.1續(xù)表備選配送中心的固定費用和總的可變費用根據(jù)朝陽重型機(jī)器有限公司提供的數(shù)據(jù)計算出各個備選配送中心的固定費用和總的可變費用，得出表4－2。

表4.2單位：萬元續(xù)表重心法實證模型的求解過程采用迭代法計算出12個重心點和重心點的運輸成本采用迭代法計算出12個重心點和重心點的運輸成本，計算結(jié)果如表4－3。表4－3采用多元線性回歸對總成本目標(biāo)函數(shù)的系數(shù)進(jìn)行求解用Excel輔助計算結(jié)果如下：圖4－1應(yīng)用excel“數(shù)據(jù)分析”功能求多元線性回歸的回歸系數(shù)由圖4－1的輸出結(jié)果，可以得到本例中的回歸系數(shù)為=0.3,=0.4,=0.3。故所求回歸方程為=采用迭代法對優(yōu)化好的模型進(jìn)行求解用迭代方法計算的結(jié)果如表4.4.續(xù)表根據(jù)上面解的情況，把求出的最優(yōu)解(162.75,56.03)之前的次優(yōu)解(162.76,56.03)、(162.76,56.02)、(162.77,56.02)以及最優(yōu)解(162.75,56.03)所對應(yīng)的位置作為配送中心的備選地址，記為Ii（i=0，1，2…，K）。且EIi=E=aidiwi