1、JIT看板生產系統(tǒng)WITNESS仿真建模和優(yōu)化實現(xiàn)1.1 生產流程描述實例系統(tǒng)生產制造單一類型的產品振動輪，其生產原材料為各種類型的鋼板，通過機械加工，然后組裝成一個空輪（沒有安裝軸承座、軸承和軸）。組成振動輪的部件主要為5類，分別為外圈、內圈、加強筋、內隔板和封口板，以數(shù)字1，2，3，4，5表示；其生產過程包括12個主要加工單元，分別為剪板切割（WS1）、打坡口（WS2）、卷圓（WS3）、輪圈焊接（WS4）、找圓（WS5）、車斷面（WS6）、數(shù)控切割（WS7）、調平（WS8）、油漆（WS9）、鉆孔（WS10）、內輪焊接（分裝工作站）和輪子焊接（總裝工作站）。其生產流程如圖5.1所示。11,5

2、產成品6需求314,53,4,53,4,51,21,21,2WS1WS21,2WS3WS4WS5WS7WS8WS9WS12WS10WS112物料流看板流WS：工作站圖5.1 實例生產/庫存系統(tǒng)生產流程圖4WS6外圈原材料為特定型號的鋼板，加工過程依次為兩塊外圈鋼板通過卷板切割機切割成適合大小，通過坡口機將鋼板兩端結合處內外打出坡度，通過卷圓機將鋼板卷成輪圈，經過輪圈焊接工段將輪圈接口處焊接起來，通過找圓機將輪圈找圓，然后進入輪子焊接工段，與內輪和封口板焊接成輪子。內圈原材料也為特定型號的鋼板，加工過程為每次三塊鋼板通過卷板切割機切割成適合大小的板材，然后打坡口、卷圓、焊接、找圓，同外圈加工過程

3、一樣。在經過找圓工段之后，內輪圈再經過數(shù)控車床，進行端面對車，經過鉆孔工段鉆出工藝孔，到內輪焊接工段與加強筋和內隔板焊接成內輪。加強筋、內隔板和封口板原材料也為特定型號的鋼板，首先經過數(shù)控切割機切割成型，然后調平、油漆，內隔板和封口板需要經過鉆口工段，鉆制工藝孔，然后，內隔板到內輪焊接工段進行與內圈和加強筋的焊接，封口板到輪子焊接工段與內輪和外圈焊接成輪子。1.2 基本生產單元的分解該生產/庫存系統(tǒng)包括四條串行線，分別為：（1）外圈加工串行線：剪板切割、打坡口、卷圓、輪圈焊接、找圓；（2）內圈加工串行線：剪板切割、打坡口、卷圓、輪圈焊接、找圓、車端面、鉆孔；（3）加強筋加工串行線：數(shù)控切割、調

4、平、油漆；（4）內隔板、封口板加工串行線：數(shù)控切割、調平、油漆、鉆孔。該生產/庫存系統(tǒng)包括兩個并行加工模塊，分別為：（1）內圈、加強筋、內隔板焊接為內輪的焊接工段，即分裝工作站WS11；（2）外圈、內輪、封口板焊接為輪子的輪子焊接工段，即總裝工作站WS12。生產流程中有資源共享和競爭的作業(yè)，如：（1）內、外圈在進行切割、打坡口、卷圓、焊接等工序時，使用的都是相同的設備；（2）加強筋、內隔板、封口板在進行切割、調平工序時，也是競爭使用相同的設備。1.3 系統(tǒng)績效指標設計系統(tǒng)績效指標考慮系統(tǒng)產出量、在制品庫存、設備稼動率和訂單滿意度四項。在進行Witness仿真程序設計的時候，將通過如下的數(shù)據(jù)模型

5、和相關的程序進行統(tǒng)計計算獲得。（1）系統(tǒng)產出量（TPUT，ThroughPUT）:統(tǒng)計在特定時間長度內，生產/庫存系統(tǒng)所產出的成品數(shù)量。該指標在Witness中可以直接在成品緩沖區(qū)的“action on input”中進行變量的累計。每當一個成品進入緩沖區(qū)，執(zhí)行式（5-1）的計算程序，即產出量增加了一個。TPUT= TPUT+1 （5-1）在仿真模型中，該統(tǒng)計任務由訂單隊列DMD.b_order的“action on input”和函數(shù)UF.tput來實現(xiàn)，具體程序參看附錄B部分。（2）在制品庫存（WIP）：用來統(tǒng)計某一時點上，生產/庫存系統(tǒng)中所擁有的原材料、半成品或成品的量，一般用價值來衡量

6、。由于生產線上具有多種物料，如外圈、內圈、加強筋、封口板、內隔板，每種類型物料所具有的價值不一樣，需要統(tǒng)計各種物料的實際價值。由于本實例并不是對實際的問題進行研究，僅僅是提供一種研究方法，所以為了便于對各種方案進行經濟性比較，根據(jù)每件產品或部件的價值，給予它們附上對應的價值權重，然后進行統(tǒng)計計算。假設加強筋最便宜，以它為1，每件物料經過一道工序加工后，每加工10分鐘，價值增加一個單位，不足10分鐘的按10分鐘計，得出實例系統(tǒng)中物料的價值權重系數(shù)如表5.1，表中下劃線前的字符為物料名稱中文首字母。表5.1 物料價值權重系數(shù)表*名稱Jqj_WS7Fkb_WS7Ngb_WS7Jqj_WS8Fkb_W

7、S8Ngb_WS8Jqj_WS9fkb_WS9ngb_WS9權重122244355名稱Wq_WS1Nq_WS1Wq_WS2Nq_WS2Wq_WS3Nq_WS3Wq_WS4Nq_WS4Wq_WS5權重303035343837444247名稱Nq_WS5Nq_WS6Nq_WS10Fkb_WS10Ngb_WS10WS11WS12權重4449528883151系統(tǒng)中的在制品庫存WIP計算公式如式（5-2）： （5-2）式中：第i種物料在各工序后的緩沖區(qū)中的數(shù)量；第i種物料在各工序機器上等待加工或正在加工的數(shù)量；第i種物料在各工序間運輸工具上的數(shù)量； 第i種物料的價值權重。在該系統(tǒng)的仿真模型中，某一仿真

8、時點上的QFi、QMi、QTi都可以通過Witness提供的函數(shù)直接得出。該指標的統(tǒng)計由目標函數(shù)UF.wip來計算，具體設計參看附錄B部分。（3）定單平均等待時間（AWT，Average Waiting Time）：用于統(tǒng)計沒能夠及時被滿足的訂單在系統(tǒng)中的平均等待時間，用來揭示隨機生產/庫存系統(tǒng)對訂單的反應靈敏度。如果某一訂單在提前期之前得到滿足，則等待時間為零；如果某一訂單在提前期之后得到滿足，則等待時間為當前時間與提前期的差值。該指標數(shù)值越低，表示顧客滿意度越高；該指標數(shù)值越高，表示顧客滿意度越低。計算公式如式（5-3）所示。WTi=Ti-Ti0TiTiTiTi （5-3）式中：n仿真期間

9、訂單總數(shù)；Ti第i批訂單實際發(fā)運時間；Ti第i批訂單預定發(fā)運時間；WTi第i批訂單的等待時間；SWT仿真期間訂單總等待時間。 在仿真模型中，這兩項統(tǒng)計任務統(tǒng)一在訂單處理單元中執(zhí)行。執(zhí)行程序體參見附錄B中DMD.demand_meet的“Actions on finish”部分。（4）設備平均稼動率（OR，Operation Rate）：用于統(tǒng)計生產系統(tǒng)中機器利用程度，仿真系統(tǒng)對機器類型的元素提供了標準狀態(tài)統(tǒng)計函數(shù)，可以統(tǒng)計機器的空閑、運行、故障和維修等狀態(tài)所占時間的百分比。在進行生產系統(tǒng)優(yōu)化過程中，需要使用的機器時間包括下列幾種： 最大操作時間：指的是設備可用的最大時間，若設備本身為廠內自購，

10、而非租借，且可完全由廠內自主使用，則最大操作時間一般為日歷時間。 負荷時間：為機器設備可稼動的時間，是最大操作時間扣除停機損失，停機損失包括計劃上的休止時間，如休假、教育訓練、保養(yǎng)等。 稼動時間：是負荷時間減去停機時間，而停機時間包括批次轉換、設備異常停止、修理、待料等時間。機器稼動率=稼動時間/負荷時間。故障間隔時間：機器每運轉多長時間就會出現(xiàn)一次故障，并需要維修。 故障維修時間：機器出現(xiàn)故障，修復所需的時間。因此，本系統(tǒng)的設備平均稼動率計算公式如式(5-4)。 (5-4)式中：m系統(tǒng)中工作單元的數(shù)量；ORj系統(tǒng)中第j個工作單元中機器設備的稼動率；在仿真模型中，該統(tǒng)計任務通過函數(shù)來實現(xiàn)，程序

11、設計參看附錄B中的UF.OR1函數(shù)體。1.4 生產線平衡在進行系統(tǒng)仿真之前，需要將每個工作單元中機器設備的設計稼動率加以平衡，這樣仿真過程不會因為設備能力不平衡，導致某些工序過快成為瓶頸工序，使前后工序的機器設備因為等待而導致利用不足。為了防止這一問題，將根據(jù)該系統(tǒng)的設計能力2500臺/年，每年250個工作日來平衡生產線。由于機器對每一部件的加工時間均值確定，只能通過設計每個加工單元中機器的數(shù)量來對生產線進行設計，使得生產線能夠很好的滿足設計的生產能力。對每個工作單元機器數(shù)量的計算過程涉及到如下的變量符號，首先對它們加以說明。TWTTotalWorkTime：每個工作單元中機器一年中的負荷時間

12、；QU QuantityUnit：每個工作單元所擁有的機器數(shù)量；TQWTotalQuantityofWork：每個工作單元所完成的產品套數(shù)；TC TimeCycle：每個單元完成一套產品的生產平均周期時間，例如：WS3設備完成的一套產品為分別加工一個外圈、一個內圈，則平均周期時間為完成一個外圈的平均周期時間（30）+完成一個內圈的平均周期時間（30），共60分鐘；TNBTotalNumofBreakdown：全年機器發(fā)生故障的次數(shù)，每月1次，均值為12次；RTRepairTime：單次機器維修所需時間，均值為半天，4小時，240分鐘；BRBreakdownRate：故障率，均值為240*12/

13、240/480=2.5%；BIBreakdownInterval：機器故障發(fā)生時間間隔；OROperating Rate：機器稼動率，假設為80%。通過該產品的BOP，得出該系統(tǒng)生產單件最終產品，在各個生產單元上需要的平均周期時間TC如表5.2所示。表5.2 生產系統(tǒng)中生產單元TC表單元WS1WS2WS3WS4WS5WS6WS10WS7WS8WS9WS11WS12時間5090601204050130838203045通過系統(tǒng)的設計生產能力和機器的生產周期，同時考慮故障比率較低，在計算機器臺數(shù)時將其忽略不記，設計計算公式如（5-5）所示。 （5-5）式（5-5）中涉及四項變量，其中由于設備的故障

14、、維修等所占時間比例忽略不記，所以機器全年負荷時間為TWT=250天*480分鐘/天=120000分鐘；假設稼動率產品套數(shù)就是生產線的年設計能力，即TQW=2500；從表5.2可以得出這12個工作單元的單套產品生產平均周期時間TCi。所以根據(jù)式（5-5）可以得出每個工作站的機器數(shù)量計算公式，如式（5-6）。 （5-6）計算式（5-6）可以得出每個加工單元所需的設備數(shù)量，如果結果是小數(shù)，取不小于結果的最小整數(shù)。通過計算得出每個單元需要配備的機器數(shù)量如表5.3所示。表5.3 加工單元設備臺數(shù)表加工單元WS1WS2WS3WS4WS5WS6WS10WS7WS8WS9WS11WS12設備數(shù)量232422

15、411112根據(jù)設備故障比率，可以計算出設備故障平均時間間隔，計算公式如式（5-7）：BI=TWT（OR-BR）/TNB=1200000.775/12=7750 （5-7）即每臺設備加工時間累計達到7750分鐘，就會發(fā)生故障，需要進行平均時間為240分鐘的維修后，才可以繼續(xù)運轉。1.5 看板數(shù)量的確定生產系統(tǒng)循環(huán)過程中看板的數(shù)量決定了生產、使用或轉讓的在制品數(shù)量，決定了系統(tǒng)運作效率的高低。建立看板控制系統(tǒng)需要確定所需要的看板卡的數(shù)量。對于本文的單看板生產系統(tǒng)需要確定生產看板的套數(shù)，以它來控制加工單元的開始生產或停止生產。得到多數(shù)研究文獻和實際應用接收的計算方式有如下兩種：文獻9的最優(yōu)看板計算公

16、式（5-8）。 （5-8）式中：k看板卡套數(shù)；D一段時期所需產品的平均數(shù)量；L補充訂貨的提前期；S安全庫存；C容器容量。文獻10的看板數(shù)量計算公式（5-9）。 （5-9）式中：R組件利用速度，如15件/分鐘；Tp單件生產時間，如：1.5分鐘/件；Td箱子的運輸和傳送時間，如：2.5分鐘；C箱子的單位容量，如：45件；X需求速率變化的百分比，如10%。使用式（5-8）來確定看板數(shù)量，由于該公式沒有應對動態(tài)需求這一項的變量，有效運用該公式的重要一步就是確定安全庫存S的數(shù)值，而一旦動態(tài)需求的隨機特征有所變化，該公式不能相應的進行反饋；式（5-9）考慮到動態(tài)需求，在公式中設定了需求速率變化X這一項，但

17、是它沒有考慮到安全庫存，計算出來的數(shù)值通常偏小，引起生產系統(tǒng)較嚴重的缺貨現(xiàn)象。本文研究生產/庫存系統(tǒng)在需求和生產服務具有隨機波動性時，隨機變量的變異系數(shù)對系統(tǒng)績效的影響，所以根據(jù)式（5-8）和式（5-9）進行修正，在式（5-9）添加安全系數(shù)項，以保持看板數(shù)量能夠對隨機變量做出反應；同時，考慮批處理作業(yè)，添加批處理數(shù)據(jù)項B，保持系統(tǒng)的有效運行。對于生產看板數(shù)量的修正公式如（5-10）式。 （5-10）式中：k生產看板數(shù)量；取大于表達式數(shù)值的最小整數(shù)；R組件利用速度；B下道工序對該部件的批處理數(shù)量；Tp單件生產時間；Td箱子的運輸和傳送時間；C箱子的單位容量；X需求速率變化的百分比；安全系數(shù)。同一

18、產品生產線的利用速率R為每天設計產出量除以每天的時間。在仿真實例中，產出量為每天10套，即R=10/480min=1/48（套/分鐘）；B和Tp各個部件各不相同，見流程數(shù)據(jù)表附錄A；由于不考慮工序間運輸時間對系統(tǒng)績效的影響，取Td為0，C為1；X就是需求隨機分布函數(shù)的變異系數(shù)CV，如果是需求確定性實驗，則X=0,如果為需求隨機性實驗，則X有需求的變異系數(shù)確定；安全系數(shù)通過仿真優(yōu)化，暫取為1.7。根據(jù)看板數(shù)量計算公式（5-9），對各個部件的各工序生產看板數(shù)量計算程序見附錄B的Initial Actions和看板計算函數(shù)UF.KB部分。1.6 仿真模型的WITNESS實現(xiàn)通過對該隨機生產/庫存系統(tǒng)

19、流程的分析和系統(tǒng)配置參數(shù)的設定，設計出系統(tǒng)的Witness仿真模型，其可視化界面如圖5.4所示。圖5.4 振動輪生產/庫存系統(tǒng)Witness仿真界面其中仿真模型中變量數(shù)據(jù)的初始化設置設計于Witness仿真模型的Initialize程序中，函數(shù)的邏輯結構、加工單元、原材料、半成品以及庫存等實體元素的輸入輸出關系等都設計在這些元素對應的action或rule編輯框內，具體程序參考附錄部分。1.7 模型參數(shù)優(yōu)化在計算機模型驗證部分，為了證明計算機模型同理論模型是相符的，只考慮了計算機模型在確定型和隨機型兩種情況下的設備稼動率和故障率，以及產品產出數(shù)量，沒有考慮系統(tǒng)的另外兩個績效指標：庫存水平和定單

20、平均等待時間。這兩項指標與看板數(shù)量和維修人員數(shù)量之間具有很強的相關性。因此為了研究各種可變因子變化對系統(tǒng)績效的影響，需要建立一個統(tǒng)一的、優(yōu)化的模型參數(shù)環(huán)境，需要對模型的相關參數(shù)進行優(yōu)化取值，涉及的參數(shù)有看板安全系數(shù)和維修工的數(shù)量，對其進行優(yōu)化選擇的準則是：在制品庫存盡可能低、定單平均等待時間盡可能低、維修工數(shù)量盡可能少。將參數(shù)和維修工的數(shù)量Qmen的可能取值設置如表5.8所示，從表上可以看出有13種可選取值，Qmen有6個可能取值，累積有78種組合，現(xiàn)在尋求仿真模型在這78種方案的哪一種方案配置下，系統(tǒng)績效最優(yōu)。表5.8 優(yōu)化參數(shù)的可能取值表變量名稱可能取值累積1.2，2.4 step=0.1

21、13Qmen1，6 step=178通過Witness仿真優(yōu)化工具包optimize進行優(yōu)化仿真運行，將實驗結果導入到Excel中，按照定單平均等待時間AWT為第一排序字段、在制品庫存WIP為第二排序字段，遞增排列。將排列結果的前9個相對最優(yōu)配置方案及實驗結果列入表5.9。表5.9 模型仿真參數(shù)優(yōu)化結果數(shù)據(jù)表WIPQmenAWTORTPUT479441.7052.0122520479441.8052.0122520503161.7052.0542520503161.8052.0542520514051.7052.0532520514051.8052.0532520523031.9052.061

22、2520523032052.0612520523032.1052.0612520從表5.9可以看出，在確定性情況下，取1.7和1.8兩個數(shù)值得出的系統(tǒng)績效相同；再比較WIP數(shù)值，最低值為4794；維修工數(shù)量直到第7條配置時，才為3人，少于第一條配置。因此，選定第一條配置為系統(tǒng)的最優(yōu)配置方案。在下一章進行的仿真實驗研究中，如果沒有特別的說明，模型中的和Qmen將確定為這兩個數(shù)值（1.7和4）。附錄B Witness 仿真程序MODEL INFORMATION:Model Name:production and inventory system simulationModel Title:opti

23、mize system parametersModel Author:jiannywangCreation Date: Thu Mar 31 11:29:17 2005Model Notes:Industrial Engineering DepartmentBusiness Administration SchoolUniversity of JiangsuInitial Actions:!pm_var.mypns = 100!pm_var.cv = 0Uvar.kbsafe = 1.7Uvar.PK_wl_qg = UF.KB (25)Uvar.PK_wl_dpk = UF.KB (50)U

24、var.PK_wl_jy = UF.KB (30)Uvar.PK_wl_hj1 = UF.KB (60)Uvar.PK_wl_zy = UF.KB (25)Uvar.PK_nl_qg = UF.KB (20)Uvar.PK_nl_dpk = UF.KB (40)Uvar.PK_nl_jy = UF.KB (25)Uvar.PK_nl_hj1 = UF.KB (50)Uvar.PK_nl_zy = UF.KB (20)Uvar.PK_nl_cdm = UF.KB (50)Uvar.PK_nl_zk = UF.KB (30)Uvar.PK_jqj_qg = UF.KB (12 * 4)Uvar.P

25、K_jqj_tp = UF.KB (14 * 4)Uvar.PK_jqj_yq = UF.KB (60 * 4)Uvar.PK_ngb_qg = UF.KB (15 * 2)Uvar.PK_ngb_tp = UF.KB (25 * 2)Uvar.PK_ngb_yq = UF.KB (60 * 2)Uvar.PK_ngb_zk = UF.KB (50 * 2)Uvar.PK_fkb_qg = UF.KB (18 * 2)Uvar.PK_fkb_tp = UF.KB (25 * 2)Uvar.PK_fkb_yq = UF.KB (60 * 2)Uvar.PK_fkb_zk = UF.KB (50

26、* 2)Uvar.PK_hj2 = UF.KB (30)Uvar.PK_hj3 = UF.KB (45)Uvar.breakinterval = 7750Uvar.repairtime = 240_ELEMENT NAME: DMD.b_orderElement Type:BufferQuantity: 1Capacity:1000Input Option:RearOutput Option:FirstSearch From:FrontActionsInput:Uvar.totaldemand = Uvar.totaldemand + DMD.batchsizeOutput:Uvar.outp

27、art = Uvar.outpart + DMD.batchsize_ELEMENT NAME: DMD.demandElement Type:PartType:Variable attributesGroup number:1Inter Arrival Time:UF.DIT () + 1First Arrival at:310.0Maximum Arrivals:UnlimitedInput / Output RulesOutput:PUSH to DMD.b_orderActionsCreate:DMD.arrivaltime = TIMEDMD.batchsize = UF.DBS (

28、) + 1DMD.leadtime = UF.DLT () + 1_ELEMENT NAME: DMD.demand_meetElement Type:MachineQuantity: 1Priority:LowestInput / Output RulesInput:SEQUENCE /Wait DMD.b_order#(1),WS12.store#(DMD.batchsize)Output:PUSH to SHIPActionsFinish:IF DMD.arrivaltime + DMD.leadtime + 20 TIME Uvar.Nunsatisfy = Uvar.Nunsatis

29、fy + 1Uvar.waittime = Uvar.waittime + TIME - DMD.arrivaltime - DMD.leadtime - 20ENDIF_ELEMENT NAME: pm_part.fengkoubanElement Type:PartType:Variable attributesGroup number:1Maximum Arrivals:0ActionsCreate:Uatt.tpnum = 4Uatt.qgnum = 8ICON = 106Uatt.skqgtime = 25Uatt.tptime = 18Uatt.zktime = 25_ELEMEN

30、T NAME: pm_part.neigebanElement Type:PartType:Variable attributesGroup number:1Maximum Arrivals:0ActionsCreate:Uatt.tpnum = 4Uatt.qgnum = 8ICON = 107Uatt.skqgtime = 25Uatt.tptime = 15Uatt.zktime = 25_ELEMENT NAME: pm_part.neilunbanElement Type:PartType:Variable attributesGroup number:1Maximum Arriva

31、ls:0ActionsCreate:Uatt.jbqgnum = 3Uatt.jbqgtime = 60Uatt.dpktime = 40Uatt.jytime = 25Uatt.hj1time = 50_ELEMENT NAME: pm_part.neiquanElement Type:PartType:Variable attributesGroup number:1Maximum Arrivals:0ActionsCreate:Uatt.zytime = 20Uatt.cdmtime = 50Uatt.zktime = 30_ELEMENT NAME: pm_part.strengthb

32、anElement Type:PartType:Variable attributesGroup number:1Maximum Arrivals:0ActionsCreate:Uatt.tpnum = 8Uatt.qgnum = 16ICON = 105Uatt.skqgtime = 12Uatt.tptime = 14_ELEMENT NAME: pm_part.wailunElement Type:PartType:Variable attributesGroup number:1Maximum Arrivals:0ActionsCreate:Uatt.zytime = 25_ELEME

33、NT NAME: pm_part.wailunbanElement Type:PartType:Variable attributesGroup number:1Maximum Arrivals:0ActionsCreate:Uatt.jbqgnum = 2Uatt.jbqgtime = 50Uatt.dpktime = 50Uatt.jytime = 30Uatt.hj1time = 60_ELEMENT NAME: repairmenElement Type:LaborQuantity always available:4_ELEMENT NAME: UF.AWTNotes: THE AV

34、ERAGE WAITING TIME OF ORDER Element Type:FunctionType:RealParameters:0ActionsExecution:RETURN Uvar.waittime / UF.numD ()_ELEMENT NAME: UF.btNotes: break time distribution funtion Element Type:FunctionType:RealParameters:2mean,Realppns,IntegerActionsExecution:Uvar.funpt = GAMMA (1 / Uvar.cv * 2,mean

35、* Uvar.cv * 2,Uvar.mypns + ppns)RETURN Uvar.funpt !mean_ELEMENT NAME: UF.DBSNotes: DEMAND BATCH SIZE Element Type:FunctionType:IntegerParameters:0ActionsExecution:RETURN GAMMA (1 / Uvar.cv * 2,29 * Uvar.cv * 2,Uvar.mypns + 14)_ELEMENT NAME: UF.DITNotes: DEMAND INTER ARRIVAL TIME Element Type:Functio

36、nType:RealParameters:0ActionsExecution:RETURN GAMMA (1 / Uvar.cv * 2,1439 * Uvar.cv * 2,Uvar.mypns + 13) !1439_ELEMENT NAME: UF.DLTNotes: DEMAND LEAD TIME Element Type:FunctionType:IntegerParameters:0ActionsExecution:RETURN GAMMA (1 / Uvar.cv * 2,1439 * Uvar.cv * 2,Uvar.mypns + 15) !1439_ELEMENT NAM

37、E: UF.KBElement Type:FunctionType:IntegerParameters:1meantime,RealActionsExecution:Uvar.fun_kbsum = 1 / 48 * meantime * (1 + Uvar.cv + Uvar.kbsafe) + 1RETURN Uvar.fun_kbsum!x=the coeffection variation!a=the safe index of kanban!p=the process time multi batchsize_ELEMENT NAME: UF.numDElement Type:Fun

38、ctionType:IntegerParameters:0ActionsExecution:RETURN TOTALIN (DMD.b_order)_ELEMENT NAME: UF.OR1Element Type:FunctionType:RealParameters:0ActionsExecution:RETURN (SUTIL (WS9.mach,2) + SUTIL (WS1.mach,2) + SUTIL (WS8.mach,2) + SUTIL (WS5.mach,2) + SUTIL (WS6.mach,2) + SUTIL (WS10.mach,2) + SUTIL (WS7.

39、mach,2) + SUTIL (WS12.mach,2) + SUTIL (WS4.mach,2) + SUTIL (WS11.mach,2) + SUTIL (WS3.mach,2) + SUTIL (WS2.mach,2) / 12_ELEMENT NAME: UF.ptNotes: process time distribution funtion Element Type:FunctionType:RealParameters:2mean,Realppns,IntegerActionsExecution:Uvar.funpt = GAMMA (1 / Uvar.cv * 2,mean

40、 * Uvar.cv * 2,Uvar.mypns + ppns)RETURN Uvar.funpt !mean_ELEMENT NAME: UF.rtNotes: repair time distribution funtion Element Type:FunctionType:RealParameters:2mean,Realppns,IntegerActionsExecution:Uvar.funpt = GAMMA (1 / Uvar.cv * 2,mean * Uvar.cv * 2,Uvar.mypns + ppns)RETURN Uvar.funpt !mean_ELEMENT

41、 NAME: UF.TPUTElement Type:FunctionType:IntegerParameters:0ActionsExecution:RETURN Uvar.outpart_ELEMENT NAME: UF.WIPElement Type:FunctionType:RealParameters:0ActionsExecution:!add all finish product WIP valueUvar.valwip (1) = NPARTS (WS1.store) * 30 + NPARTS (WS2.store(1) * 35 + NPARTS (WS2.store(2)

42、 * 34 + NPARTS (WS3.store(1) * 38 + NPARTS (WS3.store(2) * 37 + NPARTS (WS4.store(1) * 44 + NPARTS (WS4.store(2) * 42 + NPARTS (WS5.store(1) * 47 + NPARTS (WS5.store(2) * 44 + NPARTS (WS6.store) * 49 + NPARTS (WS11.store) * 83 + NPARTS (WS12.store) * 151Uvar.valwip (1) = Uvar.valwip (1) + NPARTS (WS

43、7.store(1) * 1 + NPARTS (WS7.store(2) * 2 + NPARTS (WS7.store(3) * 2 + NPARTS (WS8.store(1) * 2 + NPARTS (WS8.store(2) * 4 + NPARTS (WS8.store(3) * 4 + NPARTS (WS9.store(1) * 3 + NPARTS (WS9.store(2) * 5 + NPARTS (WS9.store(3) * 5 + NPARTS (WS10.store(1) * 52 + NPARTS (WS10.store(2) * 8 + NPARTS (WS

44、10.store(3) * 8 + NPARTS (DMD.demand_meet) * 151!add all material in machine valueUvar.valwip (3) = NPARTS2 (WS1.mach,pm_part.wailunban,0) * 27 + NPARTS2 (WS1.mach,pm_part.neilunban,0) * 28 + NPARTS2 (WS2.mach,pm_part.wailunban,0) * 30 + NPARTS2 (WS1.mach,pm_part.neilunban,0) * 30 + NPARTS2 (WS3.mac

45、h,pm_part.wailunban,0) * 35 + NPARTS2 (WS3.mach,pm_part.neilunban,0) * 34 + NPARTS2 (WS4.mach,pm_part.wailunban,0) * 38 + NPARTS2 (WS4.mach,pm_part.neilunban,0) * 37 + NPARTS2 (WS5.mach,pm_part.wailun,0) * 44 + NPARTS2 (WS5.mach,pm_part.neiquan,0) * 42 + NPARTS2 (WS6.mach,pm_part.neiquan,0) * 44 + N

46、PARTS2 (WS12.mach,pm_part.fengkouban,0) * 8 + NPARTS2 (WS12.mach,pm_part.wailun,0) * 47 + NPARTS2 (WS12.mach,pm_part.neilun,0) * 83Uvar.valwip (3) = Uvar.valwip (3) + NPARTS (WS7.mach) * 6 + NPARTS2 (WS8.mach,pm_part.strengthban,0) + NPARTS2 (WS8.mach,pm_part.fengkouban,0) * 2 + NPARTS2 (WS8.mach,pm

47、_part.neigeban,0) * 2 + NPARTS2 (WS9.mach,pm_part.strengthban,0) * 2 + NPARTS2 (WS9.mach,pm_part.fengkouban,0) * 4 + NPARTS2 (WS9.mach,pm_part.neigeban,0) * 4 + NPARTS2 (WS10.mach,pm_part.fengkouban,0) * 5 + NPARTS2 (WS10.mach,pm_part.neigeban,0) * 5 + NPARTS2 (WS10.mach,pm_part.neiquan,0) * 49 + NP

48、ARTS2 (WS11.mach,pm_part.strengthban,0) * 3 + NPARTS2 (WS11.mach,pm_part.neiquan,0) * 52 + NPARTS2 (WS11.mach,pm_part.neigeban,0) * 8RETURN Uvar.valwip (1) + Uvar.valwip (2) + Uvar.valwip (3)_ELEMENT NAME: WS1.machNotes: IF NPARTS (b_finish(1) + NPARTS2 (mach,wailunban,0) = 2 PULL from b_raw(1) ELSE

49、IF NPARTS (b_finish(1) + NPARTS2 (mach,wailunban,0) NPARTS (b_finish(2) + NPARTS2 (mach,neilunban,0) AND NPARTS (b_raw(2) = 3 PULL from b_raw(2) ELSE Wait ENDIF Element Type:MachineQuantity: 2Priority:LowestType:BatchBatch Min:Uatt.jbqgnumBatch Max:Uatt.jbqgnumCycle Time:UF.pt (Uatt.jbqgtime,1)Input

50、 / Output RulesInput:IF N = 1 AND NPARTS (WS1.store(1) Uvar.PK_wl_qgPULL from pm_part.wailunban out of WORLDELSEIF N = 2 AND NPARTS (WS1.store(2) Uvar.PK_wl_qgPULL from pm_part.neilunban out of WORLDELSE WaitENDIFOutput:PUSH pm_part.wailunban to WS1.store(1),pm_part.neilunban to WS1.store(2)Labor Re

51、quirementsRepair #1:repairmen#1BreakdownDescription:Breakdown Number 1Breakdown type:Busy TimeAt Start of Cycle:YesDown Interval:UF.bt (Uvar.breakinterval,21)Repair Time:UF.rt (Uvar.repairtime,41)Scrap Part:NoSet-up on repair:No_ELEMENT NAME: WS10.machElement Type:MachineQuantity: 4Priority:LowestTy

52、pe:SingleCycle Time:UF.pt (Uatt.zktime,10)Input / Output RulesInput:IF NPARTS2 (WS10.mach,pm_part.neiquan,0) + NPARTS (WS10.store(1) Uvar.PK_nl_zk OR NPARTS2 (WS10.mach,pm_part.neigeban,0) + NPARTS (WS10.store(2) Uvar.PK_ngb_zk OR NPARTS2 (WS10.mach,pm_part.fengkouban,0) + NPARTS (WS10.store(3) 0 AN

53、D 2 * (NPARTS2 (WS10.mach,pm_part.neiquan,0) + NPARTS (WS10.store(1) = NPARTS2 (WS10.mach,pm_part.neigeban,0) + NPARTS (WS10.store(2) AND 2 * (NPARTS2 (WS10.mach,pm_part.neiquan,0) + NPARTS (WS10.store(1) 0 AND NPARTS (WS10.store(2) = NPARTS (WS10.store(3)PULL from WS9.store(2)ELSEIF NPARTS2 (WS10.m

54、ach,pm_part.fengkouban,0) + NPARTS (WS10.store(3) Uvar.PK_fkb_qgPULL from WS9.store(3)ENDIFELSE WaitENDIFOutput:PUSH pm_part.neiquan to WS10.store(1),pm_part.neigeban to WS10.store(2),pm_part.fengkouban to WS10.store(3)Labor RequirementsRepair #1:repairmen#1BreakdownDescription:Breakdown Number 1Bre

55、akdown type:Busy TimeAt Start of Cycle:YesDown Interval:UF.bt (Uvar.breakinterval,30)Repair Time:UF.rt (Uvar.repairtime,50)Scrap Part:NoSet-up on repair:No_ELEMENT NAME: WS11.machElement Type:MachineQuantity: 1Priority:LowestType:AssemblyAssembly Qty:7Cycle Time:UF.pt (30,11)Input / Output RulesInpu

56、t:IF NPARTS (WS11.store) Uvar.PK_hj2SEQUENCE /Wait WS10.store(1)#(1),WS10.store(2)#(2),WS9.store(1)#(4)ELSE WaitENDIFOutput:PUSH to WS11.storeLabor RequirementsRepair #1:repairmen#1BreakdownDescription:Breakdown Number 1Breakdown type:Busy TimeAt Start of Cycle:YesDown Interval:UF.bt (Uvar.breakinte

57、rval,31)Repair Time:UF.rt (Uvar.repairtime,51)Scrap Part:NoSet-up on repair:NoActionsFinish:CHANGE ALL to pm_part.neilun_ELEMENT NAME: WS12.machElement Type:MachineQuantity: 2Priority:LowestType:AssemblyAssembly Qty:4Cycle Time:UF.pt (45,12)Input / Output RulesInput:IF NPARTS (WS12.store) Uvar.PK_hj

58、3SEQUENCE /Wait WS5.store(1)#(1),WS11.store#(1),WS10.store(3)#(2)ELSE WaitENDIFOutput:PUSH to WS12.storeLabor RequirementsRepair #1:repairmen#1BreakdownDescription:Breakdown Number 1Breakdown type:Busy TimeAt Start of Cycle:YesDown Interval:UF.bt (Uvar.breakinterval,32)Repair Time:UF.rt (Uvar.repair

59、time,52)Scrap Part:NoSet-up on repair:NoActionsFinish:CHANGE ALL to pm_part.lunzi_ELEMENT NAME: WS2.machElement Type:MachineQuantity: 3Priority:LowestType:SingleCycle Time:UF.pt (Uatt.dpktime,2)Input / Output RulesInput:IF NPARTS2 (WS2.mach,pm_part.wailunban,0) + NPARTS (WS2.store(1) Uvar.PK_wl_dpk

60、OR NPARTS2 (WS2.mach,pm_part.neilunban,0) + NPARTS (WS2.store(2) 0 AND NPARTS2 (WS2.mach,pm_part.wailunban,0) + NPARTS (WS2.store(1) NPARTS2 (WS2.mach,pm_part.neilunban,0) + NPARTS (WS2.store(2)PULL from WS1.store(1)ELSEIF NPARTS2 (WS2.mach,pm_part.neilunban,0) + NPARTS (WS2.store(2) Uvar.PK_nl_dpkP