特殊包裝材料供應鏈綜合計劃

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

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

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

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

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

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

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

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

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

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

11、來制訂。你認為什么因素會影響行動和決策比如，你以為轉(zhuǎn)包商的出價會影響租用 私人空間的面積嗎朱麗葉還要考慮怎樣處理需求預測中的誤差問題，你對此有什么建議嗎特殊包裝材料制造公司供應鏈綜合計劃建模與分析一、模型的建立1 公用倉庫模型決策變量：A1t ：第 t 季度印模壓制機使用臺數(shù)t=1-12 ） （ 整數(shù)約束 ）B1tC1tD1t第 t 季度初印模壓制機雇傭人數(shù) （ 整數(shù)約束 ） 第 t 季度初印模壓制機解雇人數(shù) （整數(shù)約束 ） 第 t 季度印模壓制機員工數(shù)量 （整數(shù)約束 ）O1t：第 t 季度印模壓制機工人加班小時數(shù)F1tG1tH1tH2tH3tA2tB2tC2tD2tO2t第 t 季度塑料薄片

12、生產(chǎn)數(shù)量 （千磅）第 t 季度塑料薄片庫存數(shù)量 （千磅）第 t 季度從倉庫運往熱熔定型壓力機的塑料薄片數(shù)量（千磅）第 t 季度生產(chǎn)的并運往倉庫的塑料薄片數(shù)量 （千磅） 第 t 季度生產(chǎn)的并運往熱熔壓力機的塑料薄片數(shù)量 （千磅） 第 t 季度熱熔定型壓力機使用臺數(shù) （ 整數(shù)約束 ） 第 t 季度初熱熔定型壓力機雇傭人數(shù) （ 整數(shù)約束 ） 第 t 季度初熱熔定型壓力機解雇人數(shù) （ 整數(shù)約束 ） 第 t 季度熱熔定型壓力機員工數(shù)量 （ 整數(shù)約束 ） 第 t 季度熱熔定型壓力機工人加班小時數(shù)F2t ：第 t 季度初塑料容器生產(chǎn)數(shù)量 （ 千磅）It ：第 t 季度轉(zhuǎn)包數(shù)量 （千磅 ） Jt ：第 t 季

13、度初原材料采購量 （ 千磅） Kt：第 t 季度初季末原材料庫存 （ 千磅）Q1t：第 t 季度運往倉庫的成品Q2t：第 t 季度從熱熔定型壓力機運出的成品S2t ：第 t 季度運出倉庫的成品S1t ：第 t 季度成品的庫存約束條件：員工與設備臺數(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)度分配的約束

14、：（1）原材料庫存： 初始庫存和最終庫存為 0： K（0）=K（1）-J （1）+F1（1）=0 K（12）=0 Kt= Kt-1+Jt-F1t（2）塑料薄片庫存： 初始庫存和最終庫存為 0： G1（0）=G1（1）-H2（t）+H1（t ）-H3（t ） G1（ 12）=0 G1t=G1t-1+H2t-H1t+H3t（3）成品庫存： 初始庫存和最終庫存為 0： S1（0）=S1（1）-F2（1）+Q2（1）+S2（1）=0S1（12）=0 S1t=S1t-1+It+Q1t-S2tQ1t +Q2t= F2tH1t +H3t= F2tH2t +H3t= F1t（4）確保先使用倉庫里的庫存（先進先

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

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

17、第 t 季度熱熔定型壓力機工人加班小時數(shù)F2t ：第 t 季度初塑料容器生產(chǎn)數(shù)量 （ 千磅）It ：第 t 季度轉(zhuǎn)包數(shù)量 （千磅 ）Jt ：第 t 季度初原材料采購量 （ 千磅）Kt：第 t 季度初季末原材料庫存 （ 千磅）Q1t：第 t 季度運往倉庫的成品Q2t：第 t 季度從熱熔定型壓力機運出的成品S2t ：第 t 季度運出倉庫的成品S1t ：第 t 季度成品的庫存Wt：第 t 季度私人倉庫的租用面積（租用期限為三年）約束條件：員工與設備臺數(shù)、員工數(shù)量與解雇和雇傭的關(guān)系：6 A1t= D1tA2t= D2tD1t= D1t-1+ B1t- C1tD2t= D2t-1+ B2t- C2tA1

18、t<=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）+F1（1）=0 K（12）=0Kt= Kt-1+Jt-F1t（2）塑料薄片庫存： 初始庫存和最終庫存為 0： G1（0）=G1（1）-H2（t）+H1（t ）-H3（t ）G1（ 12）=0 G1t=G1t-1+H2t-H1t+H3t（3）成品庫存： 初始庫存和最終庫存為 0：S1（0）=S1（1）-F2（1）+Q2

19、（1）+S2（1）=0 S1（12）=0S1t=S1t-1+It+Q1t-S2tQ1t +Q2t= F2tH1t +H3t= F2tH2t +H3t= F1t（4）確保先使用倉庫里的庫存（先進先出）H1t= G1t-1（5）每季度需用的倉庫有效面積：加班約束：O1t <=60D1tO2t <=60D2t成本：（1）勞動力成本：（2）加班成本：（3）雇傭成本：（4）解雇成本：（5）原材料成本：（6）轉(zhuǎn)包成本：（7）運輸成本：（8）租金：(9) 目標函數(shù):minZC1t ：第 t 季度初印模壓制機解雇人數(shù)t=1-12 ） （ 整數(shù)約束 ）（ 整數(shù)約束 ）（ 整數(shù)約束 ）3 綜合利用公、

20、私倉庫模型 決策變量：A1t ：第 t 季度印模壓制機使用臺數(shù)B1t ：第 t 季度初印模壓制機雇傭人數(shù)D1t ：第 t 季度印模壓制機員工數(shù)量（ 整數(shù)約束 ）O1t：第 t 季度印模壓制機工人加班小時數(shù)F1tG1tG2tH1t第 t 季度塑料薄片生產(chǎn)數(shù)量 （千磅）第 t 季度塑料薄片公用倉庫庫存數(shù)量 （千磅）第 t 季度塑料薄片私人倉庫庫存數(shù)量 （千磅）第 t 季度從公用倉庫運往熱熔定型壓力機的塑料薄片數(shù)量H4t ：第 t 季度從私人倉庫運往熱熔定型壓力機的塑料薄片數(shù)量（ 千磅 ） （ 千磅 ）H2tH5tH3tA2tB2tC2tD2tO2t第 t 季度生產(chǎn)的并運往公用倉庫的塑料薄片數(shù)量（千

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

22、材料私人倉庫庫存 （千磅 ）Q1t：第 t 季度運往公用倉庫的成品Q3t：第 t 季度運往私人倉庫的成品Q2t：第 t 季度從熱熔定型壓力機運出的成品S2t ：第 t 季度運出公用倉庫的成品S4t ：第 t 季度運出私人倉庫的成品S1t ：第 t 季度成品的公用倉庫庫存S3t ：第 t 季度成品的私人倉庫庫存約束條件：員工與設備臺數(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+2O

23、2t = F2tQ2t+S2t+S4t=Zt （Zt 為每季度的總需求 ）庫存與調(diào)度分配的約束：（1）原材料庫存： 初始庫存和最終庫存為 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

24、）成品庫存： 初始庫存和最終庫存為 0：S1（0）+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）確保先使用倉庫里的庫存（先進先出） ： H1t+H4t= G1 （ t-1 ） +G2（ t-1 ）（5）每季度需用的倉庫有效面積（假定轉(zhuǎn)包的產(chǎn)品只存在公用倉庫）：加班時間約束：O1t <=60D1tO2t <=60D2t成本：（1）勞

25、動力成本：（2）加班成本：（3）雇傭成本：4）解雇成本：（5）原材料成本：（6）轉(zhuǎn)包成本：（7）運輸成本：（8）倉庫儲存費用（以季末庫存量作為平均庫存量）（9）倉庫管理費用（以季末庫存量作為平均庫存量）：（10）私人倉庫租賃費用：令11）目標函數(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

26、15858 25504 16947 19906 17602 2818218708 21971 ;enddata! 約束條件 ;!1. 員工與設備臺數(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);f

27、or(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(p)=0);for(season(t):for(season(p)|p#eq#t-

28、1:K(t)-J(t)+F1(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(

29、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) ;! 確保先使用倉庫里的庫存 ;!4. 加班約束 ;for(season(t):O1(t) <=60*D1(t);for(season(t):O2(t)

30、<=60*D2(t);! 成本 ;!1. 勞動力成本 ; M1=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(seaso

31、n(t):2*(H1(t)+H2(t)+I(t);!8. 倉庫儲存費用 ; M8=sum(season(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);! 目標函數(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(B

32、2);for(season:gin(C2);for(season:gin(D2);End(2)運行結(jié)果Global optimal solution found.Objective value: +08Objective bound: +08Infeasibilities:Reduced CostExtended solver steps: 441 Total solver iterations:12413Variable ValueM1 5329800.M2M3M4M5 1175687.M6M7M8 1394828.M9 1585898.Z( 1) Z( 2) Z( 3) Z( 4) Z(

33、5) Z( 6) Z( 7) Z( 8) Z( 9) Z( 10) Z( 11) Z( 12) A1( 1) A1( 2) A1( 3) A1( 4) 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)

34、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) 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(

35、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)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

36、)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) 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( 1

37、0) 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) 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

38、) 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) Q2( 11) Q2( 12) S1( 1) S1(

39、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, S2,W;ends

40、etsdata:Z=14112 22826 15187 17841 15858 25504 16947 19906 17602 2818218708 21971 ;enddata! 約束條件 ;!1. 員工與設備臺數(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)= D

41、2(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(p)= 0);

42、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:S1(p

43、)=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

44、):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(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(se

45、ason(t):60*I(t);!7. 運輸成本 ; M7=sum(season(t):2*(H1(t)+H2(t)+I(t);!8. 租金 ;Y=max(season(t):W);M8=12*3*Y;! 目標函數(shù) ;min=M1+M2+M3+M4+M5+M6+M7+M8;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é)

46、果Local optimal solution found.Objective value:Objective bound:Infeasibilities:Extended solver steps:Total solver iterations:VariableM1M2M3M4M5M6M7YM8Z( 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)7380320.7380320.544103334Value5382720.1173220.A1( 4) A1( 5) A1( 6) A1(

47、 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

48、)O1( 1)O1( 2)O1( 3)O1( 4)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

49、)H2( 1)H2( 2)H2( 3)H2( 4)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

50、) B2( 12) C2( 1) C2( 2) C2( 3) C2( 4) 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

51、( 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) 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)Q

52、1( 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)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)W( 1)W( 2)W( 3)W( 4)W( 5)

53、W( 6)W( 7)W( 8)W( 9)W( 10)W( 11)W( 12)3 綜合利用公、私倉庫( 1) lingo 代碼! 公共 + 私人的綜合模型 ;model:sets:!12 個季度的需求量 ;season/1.12/:Z,A1,B1,C1,D1,O1,F1,G1,H1,H2,H3,H4,H5,A2,B2,C2,D2,O2,F2,G2,I,J,K1 ,K2,Q1,Q2,Q3,S1,S2,S3,S4,W;endsetsdata:Z=14112 22826 15187 17841 15858 25504 16947 19906 17602 28182 18708 21971 ;enddata ! 約束條件 ;!1. 員工與設備臺數(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(