□ 謝雪妍

(華南理工大學(xué) 經(jīng)濟(jì)與貿(mào)易學(xué)院，廣東 廣州 510006)

1 引言

對(duì)于服裝企業(yè)來(lái)說(shuō)，需要進(jìn)行產(chǎn)品設(shè)計(jì)與研發(fā)、面料質(zhì)量檢驗(yàn)、產(chǎn)品生產(chǎn)加工、衣服款式確定、產(chǎn)品銷售等環(huán)節(jié)。由于服裝行業(yè)的特殊性，從產(chǎn)品設(shè)計(jì)研發(fā)到最后送到消費(fèi)者的手中，大約需要一年的時(shí)間，而生產(chǎn)環(huán)節(jié)前期需要投入大量資金購(gòu)置機(jī)器設(shè)備，且后期的場(chǎng)地租金和人工成本也較大，因此，企業(yè)會(huì)將生產(chǎn)環(huán)節(jié)外包給供應(yīng)商完成，而將核心競(jìng)爭(zhēng)力放在產(chǎn)品研發(fā)和銷售環(huán)節(jié)。

本文研究的D公司采取的是生產(chǎn)外包模式，把生產(chǎn)面料及成衣的任務(wù)交給供應(yīng)商完成，把自己的優(yōu)勢(shì)資源集中在設(shè)計(jì)研發(fā)和終端銷售兩大模塊。生產(chǎn)制造是供應(yīng)鏈中重要的一環(huán)，因此，根據(jù)外包供應(yīng)商的成本、質(zhì)量和交貨能力等指標(biāo)對(duì)其進(jìn)行生產(chǎn)訂單的分配任務(wù)顯得尤為重要，能夠達(dá)到降低生產(chǎn)環(huán)節(jié)成本、提高產(chǎn)品質(zhì)量和準(zhǔn)時(shí)交貨比率、提高客戶滿意度的效果。

2 相關(guān)數(shù)據(jù)分析

本文選用了D服裝公司及其生產(chǎn)外包供應(yīng)商2018年7月-12月的相關(guān)數(shù)據(jù)，主要討論泳衣、運(yùn)動(dòng)上衣和運(yùn)動(dòng)內(nèi)衣3款產(chǎn)品及4個(gè)供應(yīng)商之間的訂單分配情況。

表1 供應(yīng)商信息匯總表

表2 D服裝公司對(duì)各產(chǎn)品的生產(chǎn)需求

3 實(shí)際訂單分配

在實(shí)際訂單分配任務(wù)中，D服裝公司會(huì)根據(jù)工業(yè)采購(gòu)中心計(jì)劃小組計(jì)劃員的經(jīng)驗(yàn)進(jìn)行分配，主要將供應(yīng)商分為長(zhǎng)期合作供應(yīng)商和新供應(yīng)商兩類，與D服裝公司合作超過(guò)兩年的供應(yīng)商為長(zhǎng)期合作供應(yīng)商，不超過(guò)兩年的為新供應(yīng)商，其中供應(yīng)商4(HF)為長(zhǎng)期合作供應(yīng)商，供應(yīng)商1(DR)、供應(yīng)商2(DJ)和供應(yīng)商3(HL)為新供應(yīng)商。

長(zhǎng)期供應(yīng)商其所分配的訂單量大約是新供應(yīng)商的2倍，同一類型的供應(yīng)商平均分配訂單量。由表2得2018年7月-12月期間3類產(chǎn)品的需求量，其中品類1(運(yùn)動(dòng)內(nèi)衣)為14500、品類2(運(yùn)動(dòng)上衣)為38900，品類3(連體泳衣)為11400。表3為各供應(yīng)商實(shí)際訂單分配情況。

表3 供應(yīng)商實(shí)際訂單分配情況

根據(jù)表1的產(chǎn)品報(bào)價(jià)、固定成本、產(chǎn)品不合格率、交貨延遲率等指標(biāo)數(shù)據(jù)，結(jié)合表3的訂單分配情況，可計(jì)算出D公司實(shí)際訂單分配過(guò)程中各品類的生產(chǎn)報(bào)價(jià)總成本、固定總成本、不合格產(chǎn)品數(shù)量和延遲交貨產(chǎn)品數(shù)量。計(jì)算公式為：

總成本=生產(chǎn)報(bào)價(jià)*分配任務(wù)量+固定成本

不合格產(chǎn)品數(shù)量=產(chǎn)品不合格率*分配任務(wù)量

延遲交貨產(chǎn)品數(shù)量=交貨延遲率*分配任務(wù)量

4 模型構(gòu)建與求解

4.1 模型影響因素及目標(biāo)函數(shù)權(quán)重確定

影響D公司對(duì)生產(chǎn)供應(yīng)商訂單選擇的因素有很多，通過(guò)問(wèn)卷調(diào)查方式選取了7類因素，分別是成本、質(zhì)量、交貨能力、產(chǎn)能、供應(yīng)均衡、服務(wù)、信息化水平。通過(guò)D公司負(fù)責(zé)對(duì)接生產(chǎn)供應(yīng)商的廣州工業(yè)采購(gòu)中心的40名員工對(duì)7類因素進(jìn)行評(píng)分，結(jié)合因素能否量化和企業(yè)內(nèi)部數(shù)據(jù)獲取的難易程度，選擇模型的參數(shù)和目標(biāo)函數(shù)。

圖1 訂單分配影響因素評(píng)分

通過(guò)調(diào)查問(wèn)卷的統(tǒng)計(jì)結(jié)果選擇平均得分排名前三的成本、質(zhì)量和交貨能力三個(gè)因素作為多目標(biāo)規(guī)劃模型的目標(biāo)函數(shù)，產(chǎn)能和供應(yīng)均衡兩個(gè)因素作為模型參數(shù)在約束條件中體現(xiàn)，服務(wù)和信息化水平兩個(gè)因素因?yàn)檩^難量化、且得分排名靠后，因此不作考慮。

接下來(lái)同樣利用問(wèn)卷調(diào)查法對(duì)質(zhì)量、成本、交貨能力相對(duì)重要程度進(jìn)行評(píng)價(jià)，整理后的結(jié)果為成本相比質(zhì)量的重要程度為2.55，成本相比交貨能力的重要程度為2.4，質(zhì)量相比交貨能力的重要程度為1.75，質(zhì)量相比成本的重要程度為1/2.55=0.39，交貨能力相比成本的重要程度為1/2.4=0.42，交貨能力相比質(zhì)量的重要程度為1/1.75=0.57。整理形成以下判斷矩陣：

Matlab軟件運(yùn)行結(jié)果：w=[0.5490,0.2645,0.1866]T，CR=0.0379<0.1，因此矩陣通過(guò)一致性檢驗(yàn)?？傻贸杀尽①|(zhì)量、交貨能力的權(quán)重分別為0.5490,0.2645和0.1866，此數(shù)據(jù)可用在后續(xù)多目標(biāo)規(guī)劃模型確定各目標(biāo)的占比。

4.2 多目標(biāo)隨機(jī)規(guī)劃優(yōu)化模型建立

通過(guò)構(gòu)建隨機(jī)需求下的多目標(biāo)混合整數(shù)規(guī)劃模型，對(duì)選出的4個(gè)優(yōu)質(zhì)供應(yīng)商進(jìn)行3類產(chǎn)品的訂單任務(wù)分配，其中目標(biāo)函數(shù)為生產(chǎn)報(bào)價(jià)加上固定成本得到的總成本最小，以不合格產(chǎn)品數(shù)量為指標(biāo)使質(zhì)量最優(yōu)，以延遲交貨產(chǎn)品數(shù)量為指標(biāo)使交貨能力最強(qiáng)為三個(gè)目標(biāo)函數(shù)。

以下定義與模型相關(guān)的假設(shè)條件、符號(hào)、決策變量和相關(guān)參數(shù)。

模型假設(shè)：

①不同供應(yīng)商都可以生產(chǎn)每種產(chǎn)品；

②產(chǎn)品的需求服從正態(tài)分布；

③供應(yīng)商有供應(yīng)能力限制，即在滿交貨期條件下的最大供應(yīng)量；

④生產(chǎn)環(huán)節(jié)以成本、質(zhì)量、交貨能力三個(gè)指標(biāo)為目標(biāo)，以其他因素作為約束進(jìn)行訂單分配的研究。

符號(hào)定義：

i：供應(yīng)商序號(hào)，i=1,2，…n：

j：產(chǎn)品序號(hào)，j=1,2，…m；

n：供應(yīng)商的總數(shù)量；

m：產(chǎn)品的種類數(shù)；

決策變量：

yi：企業(yè)是否選擇供應(yīng)商i負(fù)責(zé)產(chǎn)品的生產(chǎn)，若是，yi=1；否則，yi=0

xij：供應(yīng)商i是否生產(chǎn)產(chǎn)品j，若是，xij=1；否則，xij=0

aij：供應(yīng)商i生產(chǎn)產(chǎn)品j的任務(wù)量

相關(guān)參數(shù)：

pij：供應(yīng)商i對(duì)產(chǎn)品j的單位生產(chǎn)報(bào)價(jià)；

fi：供應(yīng)商i每個(gè)生產(chǎn)周期的總固定成本；

bij：供應(yīng)商i供應(yīng)產(chǎn)品j的供應(yīng)能力，即最大供應(yīng)量；

Dj：產(chǎn)品j的階段需求量(服從正態(tài)分布N(μ，σ2)，均值為μ、方差為σ2)；

dij：供應(yīng)商i對(duì)產(chǎn)品j的交貨延遲率；

qij：從供應(yīng)商i訂購(gòu)供應(yīng)產(chǎn)品j的不合格率；

hij：供應(yīng)商i對(duì)產(chǎn)品j的缺額交貨率；

確定目標(biāo)函數(shù)：

①總成本最?。嚎偝杀?生產(chǎn)報(bào)價(jià)成本+固定成本最小化；

②產(chǎn)品質(zhì)量最優(yōu)：所有產(chǎn)品不合格數(shù)量最小化；

③交貨能力最強(qiáng)：產(chǎn)品交貨延遲數(shù)量最小化。

約束條件：

①達(dá)到一定比例的產(chǎn)品準(zhǔn)時(shí)交付(準(zhǔn)時(shí)交貨率不小于90%)；

②各供應(yīng)商對(duì)每款產(chǎn)品的足額交貨率達(dá)到一定比例(缺額交貨率不超過(guò)8%)；

hij≤0.08

③保持訂單的均衡分配(單個(gè)供應(yīng)商生產(chǎn)總量小于等于總需求量的75%)；

④供需平衡(各供應(yīng)商對(duì)于某產(chǎn)品的總生產(chǎn)供應(yīng)量≥該產(chǎn)品的需求量)；

⑤滿足供應(yīng)商產(chǎn)能負(fù)荷(最大供應(yīng)量≥生產(chǎn)供應(yīng)量)。

bij≥aij

4.3 模型的求解

本文的訂單分配為隨機(jī)需求下的多目標(biāo)混合整數(shù)規(guī)劃模型，因此，首先要利用機(jī)會(huì)約束處理隨機(jī)需求，使得模型變?yōu)橐话愕亩嗄繕?biāo)規(guī)劃模型。接著，利用最大滿意度法確定各目標(biāo)的隸屬度函數(shù)，將多目標(biāo)函數(shù)轉(zhuǎn)化為單目標(biāo)函數(shù)進(jìn)行求解。此模型目標(biāo)函數(shù)中的成本、延遲交貨數(shù)量和產(chǎn)品不合格數(shù)量的單位不相同，因此，在求解多目標(biāo)函數(shù)時(shí)會(huì)出現(xiàn)問(wèn)題。運(yùn)用最大滿意度法確定隸屬度函數(shù)能夠解決多目標(biāo)規(guī)劃模型中各目標(biāo)函數(shù)單位不一致的問(wèn)題，最終能綜合三個(gè)目標(biāo)函數(shù)進(jìn)行多目標(biāo)函數(shù)的求解。

4.3.1 利用機(jī)會(huì)約束規(guī)劃處理隨機(jī)需求

機(jī)會(huì)約束規(guī)劃是指當(dāng)約束條件中含有隨機(jī)變量時(shí)，隨機(jī)變量出現(xiàn)不滿足約束條件的可能，這時(shí)需要使約束條件成立的概率不小于某一置信水平。D公司產(chǎn)品的需求量為隨機(jī)需求，本文把它簡(jiǎn)化成服從正態(tài)分布，利用機(jī)會(huì)約束可將隨機(jī)變量等價(jià)轉(zhuǎn)換為確定性變量：

其中，當(dāng)落在置信區(qū)間的概率為90%時(shí)，對(duì)應(yīng)概率函數(shù)的反函數(shù)?-1(β)為1.28。

帶入軟件計(jì)算時(shí)取等號(hào)。

4.3.2 計(jì)算各目標(biāo)函數(shù)最大值和最小值

將表1和2數(shù)據(jù)代入Lingo軟件運(yùn)行：

maxZ1=1884682,minZ1=1851483,maxZ2=1351.310,minZ2=1080.532,maxZ3=1734,minZ3=1063.130

4.3.3 確定隸屬度函數(shù)

隸屬度函數(shù)是模糊綜合評(píng)價(jià)函數(shù)中用來(lái)綜合評(píng)價(jià)多個(gè)影響因素的方法，利用理想向量確定多目標(biāo)規(guī)劃函數(shù)中各目標(biāo)的隸屬度函數(shù)。

將各目標(biāo)函數(shù)的最大值和最小值代入上述公式可得：

4.3.4 多目標(biāo)函數(shù)轉(zhuǎn)化為單目標(biāo)函數(shù)求解

確定了各目標(biāo)的隸屬度函數(shù)后可利用最大滿意度法將多目標(biāo)函數(shù)轉(zhuǎn)化為單目標(biāo)函數(shù)。定義α為最大滿意度，則最大滿意度滿足：

α=max min{u[Z1(x)],u[Z2(x)],u[Z3(x)]}

目標(biāo)函數(shù)為maxα，約束條件為：

α≤u[Zl(x)],i=1,2,…
g(x)≤0

4.3.5 目標(biāo)函數(shù)加權(quán)處理

利用4.1得出的各目標(biāo)重要程度權(quán)重對(duì)最大滿意度進(jìn)行加權(quán)可得：

maxα=max{0.5490u[Z1(x)]+0.2645u[Z2(x)]+0.1866u[Z3(x)]}

圖2 Lingo軟件運(yùn)行結(jié)果

表4 訂單分配模型結(jié)果

4.4 模型結(jié)果分析

利用多目標(biāo)混合整數(shù)隨機(jī)規(guī)劃模型對(duì)D公司訂單分配任務(wù)進(jìn)行求解，相比于原來(lái)僅憑經(jīng)驗(yàn)分配訂單，使得訂單分配有多個(gè)方面的優(yōu)化效果。

4.4.1 優(yōu)化生產(chǎn)環(huán)節(jié)的總成本、產(chǎn)品質(zhì)量和產(chǎn)品交貨能力

表5 優(yōu)化前后的成本、質(zhì)量、交貨能力對(duì)比

其中，總成本降低了2.3%，產(chǎn)品不合格數(shù)量和交貨延遲數(shù)量分別降低了48.6%和51.4%，證明本文所運(yùn)用的供應(yīng)商選擇方法和訂單分配模型對(duì)D公司生產(chǎn)環(huán)節(jié)具有優(yōu)化效果。

4.4.2 考慮供應(yīng)商的綜合能力

由表4可得，在產(chǎn)品類別3(連體泳衣)中，供應(yīng)商DJ比供應(yīng)商DR的成本要高，但因?yàn)楣?yīng)商DJ的質(zhì)量和交貨能力方面比供應(yīng)商DR更有優(yōu)勢(shì)，因此分配給DJ的量更多。這說(shuō)明了在多目標(biāo)規(guī)劃模型所選定的是成本、質(zhì)量和交貨能力三個(gè)目標(biāo)而非單一目標(biāo)，并按各權(quán)重得到最終結(jié)果，這樣能夠使D公司在多個(gè)方面對(duì)供應(yīng)商的綜合能力做出評(píng)價(jià)，從而更有依據(jù)地完成訂單分配的任務(wù)。而在D公司實(shí)際訂單分配流程中，沒(méi)有以目標(biāo)為導(dǎo)向，簡(jiǎn)單地進(jìn)行平均分配，沒(méi)有依據(jù)性。

4.4.3 考慮各目標(biāo)的重要程度

分析表1和表4數(shù)據(jù)可發(fā)現(xiàn)，在產(chǎn)品類別1(運(yùn)動(dòng)內(nèi)衣)中，供應(yīng)商2(DJ)的交貨延遲率為0.8%，相比另外三個(gè)供應(yīng)商來(lái)說(shuō)交貨能力最有優(yōu)勢(shì)，但因?yàn)?.1小節(jié)中所算出交貨能力的權(quán)重較小，供應(yīng)商2(DJ)在產(chǎn)品1中的綜合能力排名仍然靠后，因此所分配到的訂單量較少。因此，不同目標(biāo)的權(quán)重對(duì)模型結(jié)果十分重要，在確定權(quán)重時(shí)要全方位考慮，如本文中的權(quán)重采用D公司工業(yè)采購(gòu)中心的40名內(nèi)部員工的評(píng)分?jǐn)?shù)據(jù)，他們與生產(chǎn)供應(yīng)商的接觸最密切，評(píng)價(jià)出來(lái)的分?jǐn)?shù)具有說(shuō)服力。

在實(shí)際訂單分配流程中，D公司僅憑計(jì)劃員和供應(yīng)鏈經(jīng)理等人的經(jīng)驗(yàn)和主觀意見進(jìn)行訂單分配任務(wù)。而此優(yōu)化模型綜合考慮了影響訂單分配的多方面因素，且模型數(shù)據(jù)容易獲取，能夠?qū)公司的訂單分配任務(wù)進(jìn)行改進(jìn)和優(yōu)化，且在日后能定期統(tǒng)計(jì)更新供應(yīng)商的數(shù)據(jù)，使模型結(jié)果更準(zhǔn)確。

5 總結(jié)

在訂單分配優(yōu)化模型中，舍棄了D服裝公司參考以往訂單分配方案和依賴計(jì)劃小組計(jì)劃員進(jìn)行訂單分配的主觀方法，改為采用隨機(jī)需求下的多目標(biāo)規(guī)劃模型，設(shè)置了成本最小、質(zhì)量最優(yōu)和交貨能力最強(qiáng)三個(gè)目標(biāo)，供應(yīng)商最大供應(yīng)量、D公司產(chǎn)品需求量、訂單均衡等因素為約束條件，根據(jù)目標(biāo)函數(shù)和約束條件分配生產(chǎn)任務(wù)，接著運(yùn)用機(jī)會(huì)約束和最大滿意度法求解多目標(biāo)隨機(jī)規(guī)劃模型，最后得到訂單分配結(jié)果。

通過(guò)實(shí)證分析結(jié)果表明，本文訂單分配優(yōu)化策略有效地降低成本、提高質(zhì)量和準(zhǔn)時(shí)交貨率。