特殊包裝材料供應(yīng)鏈綜合計劃

1、目錄問題描述： . .錯誤 ! 未定義書簽。特殊包裝材料制造公司供應(yīng)鏈綜合計劃建模與分析. .錯誤 ! 未定義書簽。一、模型的建立 . .錯誤 ! 未定義書簽。1公用倉庫模型 . .錯誤 ! 未定義書簽。二、 LINGO編程求解 .錯誤 ! 未定義書簽。1只租用公用倉庫 . .錯誤 ! 未定義書簽。2只租用私人倉庫 . .錯誤 ! 未定義書簽。3綜合利用公、私倉庫 . .錯誤 ! 未定義書簽。三、結(jié)果分析 . .錯誤 ! 未定義書簽。四、影響決策的各類因素 . .錯誤 ! 未定義書簽。五、需求預(yù)測中的誤差項處理 . .錯誤 ! 未定義書簽。問題描述：特殊包裝材料制造公司案例分析當(dāng)朱麗葉

2、3;威廉姆斯（ Julie Williams ）離開特殊包裝材料制造公司（ SPC）的會議室時，她想了很多東西。她的部門經(jīng)理告訴她，她的生產(chǎn)廠要負(fù)責(zé)加工過程中的原材料的管理，所有與這些原材料存儲有關(guān)的成本要記入工廠下一年度的成本。作為設(shè)施生產(chǎn)計劃經(jīng)理，她的責(zé)任是使生產(chǎn)總成本最小化這包括在加工的原材料的儲存。作為計劃制定者，朱麗葉的目標(biāo)是使工廠的生產(chǎn)效率最大化，同時不需考慮庫存的需要。新結(jié)構(gòu)要求她為在加工原材料的管理和儲存制定一個低成本的規(guī)劃。1特殊包裝材料制造公司SPC將聚苯乙烯樹脂轉(zhuǎn)化為食品工業(yè)使用的可回收的容器。聚苯乙烯是以樹脂小球這種商品的形式購入的。樹脂從大型鐵路集裝箱或公路拖車卸載下

3、來后，就裝入儲藏的密室。制作食品容器有兩個步驟。第一，樹脂被傳送到印模壓制機(jī)，用壓制機(jī)將聚苯乙烯做成一卷。塑料以兩種形式出現(xiàn)：透明的和黑色的。一卷塑料或者直接用來制作容器或者儲藏起來。第二，卷形塑料被裝入熱熔定型器中，它將薄片做成容器，然后修正容器。這兩個步驟如圖1 所示。第一步第二步印模熱熔儲存儲存乙壓制機(jī)定型樹脂烯薄片圖 1特殊包裝材料制造公司的生產(chǎn)流程在過去的 5 年，塑料包裝行業(yè)增長很快。對透明塑料做成的容器的需求來自食品店、燒烤店和餐館。對黑色塑料的需求來自于那些用它們作包裝和托盤的餐館和食品店。對透明塑料的需求在夏天達(dá)到高峰。印模壓制機(jī)的生產(chǎn)能力不足以滿足塑料薄片在旺季的需求。結(jié)果

4、，生產(chǎn)線只能被迫按照預(yù)期的需求將各種薄片儲存起來。朱麗葉和她的計劃小組討論了今后 3 年的需求預(yù)測和誤差預(yù)測問題。假定工廠現(xiàn)在要對所有產(chǎn)品的庫存和管理工作收費，因此，她必須決定用公共倉庫還是私人倉庫。如果采用私人倉庫，她還必須決定要租賃或修建的倉庫容量有多大。SPC公司目前每季度有 63 個工作日，每個工作日包括 8 小時和任何安排好的加班時間。2今后 3 年的需求預(yù)測朱麗葉的計劃小組已經(jīng)利用歷史數(shù)據(jù)，對透明和黑色塑料的季度需求進(jìn)行了預(yù)測。從 2001-2003 年的需求預(yù)測見下表。年份季度黑色塑料需求預(yù)測透明塑料需求預(yù)測200166507462457618250629388941377740

5、6420027509834951492035570569891153994507200383679235572122461781910889170214950MAD=608MAD=786表 5-1透明和黑色塑料容器的需求預(yù)測3印模壓制機(jī)印模壓制流程是資本密集型的，因為它需要設(shè)備的投資。目前生產(chǎn)線上有14 臺機(jī)器，每臺機(jī)器的生產(chǎn)能力是每小時加工3000 磅原料。機(jī)器在生產(chǎn)透明塑料和黑色塑料之間要有一個轉(zhuǎn)化的過程。由于這種轉(zhuǎn)化，公司估計要有5%的生產(chǎn)能力會損失掉。這樣，機(jī)器的實際生產(chǎn)能力就變?yōu)槊啃r 2850 磅。每臺機(jī)器需要 6 名工人，公司給每個工人的薪酬是每小時 15 美元，包括各種福利。加

6、班費是正常工資的 150%。每季度工人加班最多 60 小時。印模壓制機(jī)是很昂貴的，同時購買新的機(jī)器需要再額外雇用6 名工人。每臺新機(jī)器會帶來 80000 美元的固定成本。雇用的新工人還要被培訓(xùn)，培訓(xùn)費是每人 3000 美元。結(jié)果，公司決定短期內(nèi)不購買任何新的機(jī)器。每一個季度，機(jī)器都會有閑置。但是如果解雇工人的話，每解雇1 個工人，要花費2500 美元。如果利用閑置機(jī)器，公司在每個工人身上要花費3000 美元的培訓(xùn)費。4熱熔定型壓縮機(jī)目前生產(chǎn)線上有25 臺熱熔定型壓縮機(jī)。每臺機(jī)器需要一個操作者，機(jī)器以每小時 2000 磅的速度生產(chǎn)容器。公司給工人每小時 15 美元工資（包括各種福利）。加班費是正

7、常工資的 150%。工人每季度最多加班 60 小時。每季度機(jī)器都會閑置，解雇一個工人要花費 2500 美元，而培訓(xùn)一個新工人的成本是 3000 美元。5轉(zhuǎn)包生產(chǎn)公司可以將塑料的生產(chǎn)轉(zhuǎn)包出去，因為在開放市場上生產(chǎn)能力總是充足的。如果轉(zhuǎn)包給其他制造商，公司每1000 磅要支付 100 美元。6原材料管理樹脂每 1000 磅 10 美元。市場上樹脂不短缺，可以很容易的以上述價格購進(jìn)。公司每季度購買一次樹脂以滿足生產(chǎn)計劃的需求。這樣就必須將樹脂儲存起來。印模定型壓制機(jī)生產(chǎn)出成卷的塑料，熱熔定型壓制機(jī)接下來加工，生于的塑料片通過傳送帶運往兩個公共倉庫。當(dāng)熱熔定型壓制機(jī)需要原材料的時候，傳送帶又將塑料片運

8、送回來。公司總的運輸成本是每1000 磅 2 美元。7公共倉庫公共倉庫的收費包括原材料的管理費用和保存費。公司和地方倉庫簽訂合同以每 1000 磅為單位儲存原材料，原材料管理費用每 1000 磅從 4 美元 -6 美元不等，儲存費用是每季度末每 1000 磅付 10 美元 -12 美元不等。 SPC公司每年就各項費用和地方倉庫進(jìn)行商談。8私人倉庫私人倉庫的運營需要投資或者購買設(shè)備或者租用現(xiàn)有的設(shè)備。任何地點的租金取決于該地點的建設(shè)成本以及相對于長期租約的相對成本。租約一般是 3 年期的，但公司也可以利用自己談判的優(yōu)勢將期限縮短。特殊包裝材料公司有幾個租賃方案，所以這些都優(yōu)于建設(shè)一個新的倉庫。在

9、每個地方租金都是每平方米 4 美元。一般說來，每平方米可儲存1000 磅。私人倉庫同樣有固定成本和可變成本。這個私人倉庫是從一個“第三方物流提供者”那里租賃來的，租金是每季度每1000 磅塑料 3 美元。為了保證這個租金水平，特殊包裝材料公司必須將租約簽滿 3 年。結(jié)果，公司每季度都要付租金，盡管有時它并不用來儲存貨物。因此，公司在作決定的時候必須要考慮這個因素。公司在決定使用多大面積的倉庫時必須考慮幾個因素。可使用倉儲空間是指倉庫中可以用來儲存貨物的那部分空間。我們要考慮過道、運輸和卸貨的空間、行政辦公的空間和天花板的高度。儲存密度是另一個要考慮的因素。另外，公司還要考慮原材料運輸?shù)乃俣群痛?/p>

10、數(shù)，這是因為工人的技術(shù)水平和倉庫的外觀都取決于上述因素。比如，如果原材料很容易運回的話，那么倉庫在布局時就必須要加入更多的過道和臺階。9實際行動和決策朱麗葉和她的小組必須采取兩項措施。第一，根據(jù)表5-1的 3 年預(yù)測，制定總的生產(chǎn)計劃。第二，從下列3 種選擇中選出一個。繼續(xù)在公共倉庫儲存原材料。租賃一個私人倉庫來儲促上述原材料。將上述兩種方法結(jié)合起來。如果租用私人倉庫，朱麗葉必須考慮租用的平方米數(shù)。這個決策會在2001-2003 年實施。很明顯，這個決策必須和3 年期的總體計劃聯(lián)系起來考慮。理想條件下，由于上述兩個決策是相互聯(lián)系的，所以它們應(yīng)該結(jié)合起來制訂。你認(rèn)為什么因素會影響行動和決策比如，

11、你以為轉(zhuǎn)包商的出價會影響租用私人空間的面積嗎朱麗葉還要考慮怎樣處理需求預(yù)測中的誤差問題，你對此有什么建議嗎特殊包裝材料制造公司供應(yīng)鏈綜合計劃建模與分析一、模型的建立1 公用倉庫模型決策變量：A1t ：第 t季度印模壓制機(jī)使用臺數(shù)（ t=1-12）(整數(shù)約束 )B1t ：第 t季度初印模壓制機(jī)雇傭人數(shù)(整數(shù)約束 )C1t ：第 t季度初印模壓制機(jī)解雇人數(shù)( 整數(shù)約束 )D1t ：第 t季度印模壓制機(jī)員工數(shù)量( 整數(shù)約束 )O1t：第 t季度印模壓制機(jī)工人加班小時數(shù)F1t ：第 t季度塑料薄片生產(chǎn)數(shù)量(千磅)G1t：第 t季度塑料薄片庫存數(shù)量(千磅)H1t ：第 t季度從倉庫運往熱熔定型壓力機(jī)的塑

12、料薄片數(shù)量(千磅 )H2t ：第 t季度生產(chǎn)的并運往倉庫的塑料薄片數(shù)量(千磅)H3t ：第 t季度生產(chǎn)的并運往熱熔壓力機(jī)的塑料薄片數(shù)量(千磅)A2t ：第 t季度熱熔定型壓力機(jī)使用臺數(shù)( 整數(shù)約束 )B2t ：第 t季度初熱熔定型壓力機(jī)雇傭人數(shù)( 整數(shù)約束 )C2t ：第 t季度初熱熔定型壓力機(jī)解雇人數(shù)( 整數(shù)約束 )D2t ：第 t季度熱熔定型壓力機(jī)員工數(shù)量( 整數(shù)約束 )O2t：第 t季度熱熔定型壓力機(jī)工人加班小時數(shù)F2t ：第 t季度初塑料容器生產(chǎn)數(shù)量( 千磅)It ：第 t季度轉(zhuǎn)包數(shù)量 ( 千磅 )Jt ：第 t季度初原材料采購量 ( 千磅 )Kt ：第 t季度初季末原材料庫存(千磅)

13、Q1t：第 t季度運往倉庫的成品Q2t：第 t季度從熱熔定型壓力機(jī)運出的成品S2t ：第 t季度運出倉庫的成品S1t ：第 t季度成品的庫存約束條件：員工與設(shè)備臺數(shù)、員工數(shù)量與解雇和雇傭的關(guān)系：6 A1t= D1tA2t= D2tD1t= D1t-1+ B1t- C1tD2t= D2t-1+ B2t- C2tA1t<=14A2t<=25產(chǎn)能約束：3*8*63 D1t+3O1t = F1t2*8*63 D2t+2O2t = F2tQ2t+S2t=Zt (Zt為每季度的總需求)庫存與調(diào)度分配的約束：（ 1）原材料庫存：初始庫存和最終庫存為 0：K（ 0） =K（1） -J （ 1） +

14、F1（1） =0K（ 12） =0Kt= Kt-1+Jt-F1t（ 2）塑料薄片庫存：初始庫存和最終庫存為 0：G1（ 0） =G1（ 1） -H2（ t ） +H1（ t ） -H3（ t ）G1（ 12） =0G1t=G1t-1+H2t-H1t+H3t（ 3）成品庫存：初始庫存和最終庫存為 0：S1(0)=S1(1)-F2(1)+Q2(1)+S2(1)=0S1(12)=0S1t=S1t-1+It+Q1t-S2tQ1t +Q2t= F2tH1t +H3t= F2tH2t +H3t= F1t（4）確保先使用倉庫里的庫存（先進(jìn)先出）：H1t= G1t-1加班時間約束：O1t <=60D1t

15、O2t <=60D2t成本：（ 1）勞動力成本：（ 2）加班成本：（ 3）雇傭成本：（4）解雇成本：（ 5）原材料成本：（ 6）轉(zhuǎn)包成本：（ 7）運輸成本：（ 8）倉庫儲存費用（以季末庫存量作為平均庫存量）:（ 9）倉庫管理費用（以季末庫存量作為平均庫存量）：（ 10）目標(biāo)函數(shù)：2 私人倉庫模型決策變量：A1t ：第 t季度印模壓制機(jī)使用臺數(shù)（ t=1-12 ） ( 整數(shù)約束 )B1t ：第 t季度初印模壓制機(jī)雇傭人數(shù)(整數(shù)約束 )C1t ：第 t季度初印模壓制機(jī)解雇人數(shù)( 整數(shù)約束 )D1t ：第 t季度印模壓制機(jī)員工數(shù)量( 整數(shù)約束 )O1t：第 t季度印模壓制機(jī)工人加班小時數(shù)F1t

16、 ：第 t季度塑料薄片生產(chǎn)數(shù)量(千磅)G1t：第 t季度塑料薄片庫存數(shù)量(千磅)H1t ：第 t季度從倉庫運往熱熔定型壓力機(jī)的塑料薄片數(shù)量(千磅 )H2t ：第 t季度生產(chǎn)的并運往倉庫的塑料薄片數(shù)量(千磅)H3t ：第 t季度生產(chǎn)的并運往熱熔壓力機(jī)的塑料薄片數(shù)量(千磅)A2t ：第 t季度熱熔定型壓力機(jī)使用臺數(shù)( 整數(shù)約束 )B2t ：第 t季度初熱熔定型壓力機(jī)雇傭人數(shù)( 整數(shù)約束 )C2t ：第 t季度初熱熔定型壓力機(jī)解雇人數(shù)( 整數(shù)約束 )D2t ：第 t季度熱熔定型壓力機(jī)員工數(shù)量( 整數(shù)約束 )O2t：第 t季度熱熔定型壓力機(jī)工人加班小時數(shù)F2t ：第 t季度初塑料容器生產(chǎn)數(shù)量( 千磅)

17、It ：第 t季度轉(zhuǎn)包數(shù)量 ( 千磅 )Jt ：第 t季度初原材料采購量(千磅)Kt ：第 t季度初季末原材料庫存(千磅)Q1t：第 t季度運往倉庫的成品Q2t：第 t季度從熱熔定型壓力機(jī)運出的成品S2t ：第 t季度運出倉庫的成品S1t ：第 t季度成品的庫存Wt：第 t季度私人倉庫的租用面積（租用期限為三年）約束條件：員工與設(shè)備臺數(shù)、員工數(shù)量與解雇和雇傭的關(guān)系：6 A1t= D1tA2t= D2tD1t= D1t-1+ B1t- C1tD2t= D2t-1+ B2t- C2tA1t<=14A2t<=25產(chǎn)能約束：3*8*63 D1t+3O1t = F1t2*8*63 D2t+2

18、O2t = F2tQ2t+S2t=Zt (Zt為每季度的總需求)庫存與調(diào)度分配的約束：（ 1）原材料庫存：初始庫存和最終庫存為 0：K（ 0） =K（1） -J （ 1） +F1（1） =0K（ 12） =0Kt= Kt-1+Jt-F1t（ 2）塑料薄片庫存：初始庫存和最終庫存為 0：G1（ 0） =G1（ 1） -H2（ t ） +H1（ t ） -H3（ t ）G1（ 12） =0G1t=G1t-1+H2t-H1t+H3t（ 3）成品庫存：初始庫存和最終庫存為 0：S1(0)=S1(1)-F2(1)+Q2(1)+S2(1)=0S1(12)=0S1t=S1t-1+It+Q1t-S2tQ1t

19、+Q2t= F2tH1t +H3t= F2tH2t +H3t= F1t（4）確保先使用倉庫里的庫存（先進(jìn)先出）：H1t= G1t-1（5）每季度需用的倉庫有效面積：加班約束：O1t <=60D1tO2t <=60D2t成本：（ 1）勞動力成本：（ 2）加班成本：（ 3）雇傭成本：（ 4）解雇成本：（ 5）原材料成本：（ 6）轉(zhuǎn)包成本：（ 7）運輸成本：（ 8）租金：令（9）目標(biāo)函數(shù)：3 綜合利用公、私倉庫模型決策變量：A1t ：第 t季度印模壓制機(jī)使用臺數(shù)（ t=1-12 ） ( 整數(shù)約束 )B1t ：第 t季度初印模壓制機(jī)雇傭人數(shù)(整數(shù)約束 )C1t ：第 t季度初印模壓制機(jī)解雇

20、人數(shù)( 整數(shù)約束 )D1t ：第 t季度印模壓制機(jī)員工數(shù)量( 整數(shù)約束 )O1t：第 t季度印模壓制機(jī)工人加班小時數(shù)F1t ：第 t季度塑料薄片生產(chǎn)數(shù)量 ( 千磅 )G1t：第 t季度塑料薄片公用倉庫庫存數(shù)量(千磅 )G2t：第 t季度塑料薄片私人倉庫庫存數(shù)量(千磅 )H1t ：第 t季度從公用倉庫運往熱熔定型壓力機(jī)的塑料薄片數(shù)量(千磅)H4t ：第 t季度從私人倉庫運往熱熔定型壓力機(jī)的塑料薄片數(shù)量(千磅)H2t ：第 t季度生產(chǎn)的并運往公用倉庫的塑料薄片數(shù)量(千磅)H5t ：第 t季度生產(chǎn)的并運往私人倉庫的塑料薄片數(shù)量(千磅)H3t ：第 t季度生產(chǎn)的并運往熱熔壓力機(jī)的塑料薄片數(shù)量(千磅)A

21、2t ：第 t季度熱熔定型壓力機(jī)使用臺數(shù)( 整數(shù)約束 )B2t ：第 t季度初熱熔定型壓力機(jī)雇傭人數(shù)( 整數(shù)約束 )C2t ：第 t季度初熱熔定型壓力機(jī)解雇人數(shù)( 整數(shù)約束 )D2t ：第 t季度熱熔定型壓力機(jī)員工數(shù)量( 整數(shù)約束 )O2t：第 t季度熱熔定型壓力機(jī)工人加班小時數(shù)F2t ：第 t季度初塑料容器生產(chǎn)數(shù)量( 千磅)It ：第 t季度轉(zhuǎn)包數(shù)量 ( 千磅 )Jt ：第 t季度初原材料采購量 ( 千磅 )K1t ：第 t季度初季末原材料公用倉庫庫存(千磅 )K2t ：第 t季度初季末原材料私人倉庫庫存(千磅 )Q1t：第 t季度運往公用倉庫的成品Q3t：第 t季度運往私人倉庫的成品Q2t

22、：第 t季度從熱熔定型壓力機(jī)運出的成品S2t ：第 t季度運出公用倉庫的成品S4t ：第 t季度運出私人倉庫的成品S1t ：第 t季度成品的公用倉庫庫存S3t ：第 t季度成品的私人倉庫庫存約束條件：員工與設(shè)備臺數(shù)、員工數(shù)量與解雇和雇傭的關(guān)系：6 A1t= D1tA2t= D2tD1t= D1t-1+ B1t- C1tD2t= D2t-1+ B2t- C2tA1t<=14A2t<=25產(chǎn)能約束：3*8*63 D1t+3O1t = F1t2*8*63 D2t+2O2t = F2tQ2t+S2t+S4t=Zt (Zt為每季度的總需求)庫存與調(diào)度分配的約束：（ 1）原材料庫存：初始庫存和

23、最終庫存為 0：K1（ 0） +K2（ 0） =K1（ 1） +K2（ 1） -J （ 1） +F1（ 1） =0K1（ 12） +K2（ 12） =0K1t+K2t= K1 （ t-1 ）+ K2 （ t-1 ） +Jt-F1t（ 2）塑料薄片庫存：初始庫存和最終庫存為 0：G1（ 0） +G2（ 0） =G1（ 1） +G2（ 1） -H2（ t ） -H5（ t ） -H1（ t ） -H4（ t ） -H3（ t ）G1（ 12） =0G1t+G2t=G1（ t-1 ） +G2（ t-1 ） +H2t+H5t-H1t-H4t+H3t（ 3）成品庫存：初始庫存和最終庫存為 0：S1(0)

24、+S3(0)=S1(1)+S3(1)-F2(1)+Q2(1)+S2(1)+S4(1)=0S1(12)=0S1t+S3t=S1 （ t-1 ） +S3（ t-1 ） +It+Q1t+Q3t-S2t-S4tQ1t+Q3t +Q2t= F2tH1t+H4t +H3t= F2tH2t+H5t +H3t= F1t（4）確保先使用倉庫里的庫存（先進(jìn)先出）：H1t+H4t= G1 （ t-1 ） +G2（ t-1 ）（5）每季度需用的倉庫有效面積（假定轉(zhuǎn)包的產(chǎn)品只存在公用倉庫）：加班時間約束：O1t <=60D1tO2t <=60D2t成本：（ 1）勞動力成本：（ 2）加班成本：（3）雇傭成本：

25、（ 4）解雇成本：（ 5）原材料成本：（ 6）轉(zhuǎn)包成本：（ 7）運輸成本：（ 8）倉庫儲存費用（以季末庫存量作為平均庫存量）:（ 9）倉庫管理費用（以季末庫存量作為平均庫存量）：（ 10）私人倉庫租賃費用：令（ 11）目標(biāo)函數(shù)：二、 LINGO編程求解1 只租用公用倉庫（ 1） LINGO代碼model:sets:!12 個季度的需求量;season/1.12/:Z,A1,B1,C1,D1,O1,F1,G1,H1,H2,H3,A2,B2,C2,D2,O2,F2,I,J,K,Q1,Q2,S1, S2;endsetsdata:Z=14112 22826 15187 17841 15858 2550

26、4 16947 19906 17602 28182 18708 21971 ;enddata! 約束條件 ;!1.員工與設(shè)備臺數(shù)、員工數(shù)量與解雇和雇傭的關(guān)系;for(season(t):6 *A1(t)= D1(t);for(season(t):A2(t)= D2(t);for(season(t):for(season(p)|p#eq#t-1:D1(t)+C1(t)-B1(t)= D1(p);for(season(t):for(season(p)|p#eq#t-1:D2(t)+C2(t)-B2(t)= D2(p);for(season(t):A1(t)<=14);for(season(t

27、):A2(t)<=25);!2.產(chǎn)能約束 ;for(season(t):*8*63 *D1(t)+*O1(t) = F1(t);for(season(t):2*8*63 *D2(t)+2*O2(t) = F2(t);for(season(t):Q2(t)+S2(t)=Z(t);! Zt 為每季度的總需求 ;!3.庫存與調(diào)度分配的約束;! 原材料庫存 ;for(season(p)|p#eq#1:K(p)-J(p)+F1(p)= 0);for(season(p)|p#eq#12:K(p)=0);for(season(t):for(season(p)|p#eq#t-1:K(t)-J(t)+F1

28、(t)= K(p);! 塑料薄片庫存 ;! 初始塑料薄片庫存為 0; for(season(p)|p#eq#1:G1(p)-H2(p)+H1(p)-H3(p)=0); for(season(p)|p#eq#12:G1(p)=0); for(season(t):for(season(p)|p#eq#t-1:G1(t)-H2(t)+H1(t)-H3(t)=G1(p);! 成品庫存 ;! 初始成品庫存為 0;for(season(p)|p#eq#1:S1(p)-F2(p)+Q2(p)+S2(p)=0);for(season(p)|p#eq#12:S1(p)=0);for(season(t):for(

29、season(p)|p#eq#t-1:S1(t)-I(t)-Q1(t)+S2(t)=S1(p);for(season(t):Q1(t) +Q2(t)= F2(t);for(season(t):H1(t) +H3(t)= F2(t);for(season(t):H2(t) +H3(t)= F1(t);for(season(t):for(season(p)|p#eq#t-1:! 確保先使用倉庫里的庫存 ; !4. 加班約束 ; for(season(t):O1(t) <=60*D1(t); for(season(t):O2(t) <=60*D2(t);! 成本 ;!1.勞動力成本 ;M

30、1=sum(season(t):15*8*63*(D1(t)+D2(t);!2.加班成本 ;M2=sum(season(t):*(O1(t)+O2(t);!3.雇傭成本 ;M3=sum(season(t):3000*(B1(t)+B2(t);!4.解雇成本 ;M4=sum(season(t):2500*(C1(t)+C2(t);!5.原材料成本 ;M5=sum(season(t):10*J(t);!6.轉(zhuǎn)包成本 ;M6=sum(season(t):60*I(t);!7.運輸成本 ;M7=sum(season(t):2*(H1(t)+H2(t)+I(t);!8.倉庫儲存費用;M8=sum(sea

31、son(t):10*(K(t)+G1(t)+S1(t);!9. 倉庫管理費用;M9=sum(season(t):4*(J(t)+F1(t)+H1(t)+H2(t)+I(t)+S2(t)+Q1(t);! 目標(biāo)函數(shù) ;min=M1+M2+M3+M4+M5+M6+M7+M8+M9;for(season:gin(A1);for(season:gin(B1);for(season:gin(C1);for(season:gin(D1);for(season:gin(A2);for(season:gin(B2);for(season:gin(C2);for(season:gin(D2);End（2）運行結(jié)果

32、Global optimal solution found.Objective value:+08Objective bound:+08Infeasibilities:Extended solver steps:441Total solver iterations:12413VariableValueReduced CostM15329800.M2M3M4M51175687.M6M7M81394828.M91585898.Z( 1)Z( 2)Z( 3)Z( 4)Z( 5)Z( 6)Z( 7)Z( 8)Z( 9)Z( 10)Z( 11)Z( 12)A1( 1)A1( 2)A1( 3)A1( 4)

33、A1( 5)A1( 6)A1( 7)A1( 8)A1( 9)A1( 10)A1( 11)A1( 12)B1( 1)B1( 2)B1( 3)B1( 4)B1( 5)B1( 6)B1( 7)B1( 8)B1( 9)B1( 10)B1( 11)B1( 12)C1( 1)C1( 2)C1( 3)C1( 4)C1( 5)C1( 6)C1( 7)C1( 8)C1( 9)C1( 10)C1( 11)C1( 12)D1( 1)D1( 2)D1( 3)D1( 4)D1( 5)D1( 6)D1( 7)D1( 8)D1( 9)D1( 10)D1( 11)D1( 12)O1( 1)O1( 2)O1( 3)O1( 4)

34、O1( 5)O1( 6)O1( 7)O1( 8)O1( 9)O1( 10)O1( 11)O1( 12)F1( 1)F1( 2)F1( 3)F1( 4)F1( 5)F1( 6)F1( 7)F1( 8)F1( 9)F1( 10)F1( 11)F1( 12)G1( 1)G1( 2)G1( 3)G1( 4)G1( 5)G1( 6)G1( 7)G1( 8)G1( 9)G1( 10)G1( 11)G1( 12)H1( 1)H1( 2)H1( 3)H1( 4)H1( 5)H1( 6)H1( 7)H1( 8)H1( 9)H1( 10)H1( 11)H1( 12)H2( 1)H2( 2)H2( 3)H2( 4)

35、H2( 5)H2( 6)H2( 7)H2( 8)H2( 9)H2( 10)H2( 11)H2( 12)H3( 1)H3( 2)H3( 3)H3( 4)H3( 5)H3( 6)H3( 7)H3( 8)H3( 9)H3( 10)H3( 11)H3( 12)A2( 1)A2( 2)A2( 3)A2( 4)A2( 5)A2( 6)A2( 7)A2( 8)A2( 9)A2( 10)A2( 11)A2( 12)B2( 1)B2( 2)B2( 3)B2( 4)B2( 5)B2( 6)B2( 7)B2( 8)B2( 9)B2( 10)B2( 11)B2( 12)C2( 1)C2( 2)C2( 3)C2( 4)

36、C2( 5)C2( 6)C2( 7)C2( 8)C2( 9)C2( 10)C2( 11)C2( 12)D2( 1)D2( 2)D2( 3)D2( 4)D2( 5)D2( 6)D2( 7)D2( 8)D2( 9)D2( 10)D2( 11)D2( 12)O2( 1)O2( 2)O2( 3)O2( 4)O2( 5)O2( 6)O2( 7)O2( 8)O2( 9)O2( 10)O2( 11)O2( 12)F2( 1)F2( 2)F2( 3)F2( 4)F2( 5)F2( 6)F2( 7)F2( 8)F2( 9)F2( 10)F2( 11)F2( 12)I( 1)I( 2)I( 3)I( 4)I( 5

37、)I( 6)I( 7)I( 8)I( 9)I( 10)I( 11)I( 12)J( 1)J( 2)J( 3)J( 4)J( 5)J( 6)J( 7)J( 8)J( 9)J( 10)J( 11)J( 12)K( 1)K( 2)K( 3)K( 4)K( 5)K( 6)K( 7)K( 8)K( 9)K( 10)K( 11)K( 12)Q1( 1)Q1( 2)Q1( 3)Q1( 4)Q1( 5)Q1( 6)Q1( 7)Q1( 8)Q1( 9)Q1( 10)Q1( 11)Q1( 12)Q2( 1)Q2( 2)Q2( 3)Q2( 4)Q2( 5)Q2( 6)Q2( 7)Q2( 8)Q2( 9)Q2( 10

38、)Q2( 11)Q2( 12)S1( 1)S1( 2)S1( 3)S1( 4)S1( 5)S1( 6)S1( 7)S1( 8)S1( 9)S1( 10)S1( 11)S1( 12)S2( 1)S2( 2)S2( 3)S2( 4)S2( 5)S2( 6)S2( 7)S2( 8)S2( 9)S2( 10)S2( 11)S2( 12)2 只租用私人倉庫（ 1） lingo 代碼model:! 私人倉庫模型 ; sets:!12 個季度的需求量 ;season/1.12/:Z,A1,B1,C1,D1,O1,F1,G1,H1,H2,H3,A2,B2,C2,D2,O2,F2,I,J,K,Q1,Q2,S1,

39、 S2,W;endsets data:Z=14112 22826 15187 17841 15858 25504 16947 19906 17602 28182 18708 21971 ;enddata! 約束條件 ;!1.員工與設(shè)備臺數(shù)、員工數(shù)量與解雇和雇傭的關(guān)系;for(season(t):6 *A1(t)= D1(t);for(season(t):A2(t)= D2(t);for(season(t):for(season(p)|p#eq#t-1:D1(t)+C1(t)-B1(t)= D1(p);for(season(t):for(season(p)|p#eq#t-1:D2(t)+C2(t

40、)-B2(t)= D2(p);for(season(t):A1(t)<=14);for(season(t):A2(t)<=25);!2.產(chǎn)能約束 ;for(season(t):*8*63 *D1(t)+*O1(t) = F1(t);for(season(t):2*8*63 *D2(t)+2*O2(t) = F2(t);for(season(t):Q2(t)+S2(t)=Z(t);! Zt 為每季度的總需求 ;!3.庫存與調(diào)度分配的約束;! 原材料庫存 ;for(season(p)|p#eq#1:K(p)-J(p)+F1(p)= 0);for(season(p)|p#eq#12:K(

41、p)= 0);for(season(t):for(season(p)|p#eq#t-1:K(t)-J(t)+F1(t)= K(p);! 塑料薄片庫存 ;for(season(p)|p#eq#1:G1(p)-H2(p)+H1(p)-H3(p)=0);for(season(p)|p#eq#12:G1(p)=0);for(season(t):for(season(p)|p#eq#t-1:G1(t)-H2(t)+H1(t)-H3(t)=G1(p);! 成品庫存 ;for(season(p)|p#eq#1:S1(p)-F2(p)+Q2(p)+S2(p)=0);for(season(p)|p#eq#12:

42、S1(p)=0);for(season(t):for(season(p)|p#eq#t-1:S1(t)-I(t)-Q1(t)+S2(t)=S1(p);for(season(t):Q1(t) +Q2(t)= F2(t);for(season(t):H1(t) +H3(t)= F2(t);for(season(t):H2(t) +H3(t)= F1(t);! 確保先使用倉庫里的庫存 ; for(season(t): for(season(p)|p#eq#t-1: H1(t)= G1(p) ; for(season(t):W(t) = H2(t)+I(t)+Q1(t);!4.加班約束 ;for(season(t):O1(t) <=60*D1(t);for(season(t):O2(t) <=60*D2(t);! 成本 ;!1.勞動力成本 ;M1=sum(season(t):15*8*63*(D1(t)+D2(t);!2.加班成本 ;M2=sum(seas