JIT看板生產(chǎn)系統(tǒng)WITNESS仿真建模和優(yōu)化實(shí)現(xiàn)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

20、平均等待時(shí)間。這兩項(xiàng)指標(biāo)與看板數(shù)量和維修人員數(shù)量之間具有很強(qiáng)的相關(guān)性。因此為了研究各種可變因子變化對(duì)系統(tǒng)績效的影響，需要建立一個(gè)統(tǒng)一的、優(yōu)化的模型參數(shù)環(huán)境，需要對(duì)模型的相關(guān)參數(shù)進(jìn)行優(yōu)化取值，涉及的參數(shù)有看板安全系數(shù)和維修工的數(shù)量，對(duì)其進(jìn)行優(yōu)化選擇的準(zhǔn)則是：在制品庫存盡可能低、定單平均等待時(shí)間盡可能低、維修工數(shù)量盡可能少。將參數(shù)和維修工的數(shù)量Qmen的可能取值設(shè)置如表5.8所示，從表上可以看出有13種可選取值，Qmen有6個(gè)可能取值，累積有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進(jìn)行優(yōu)化仿真運(yùn)行，將實(shí)驗(yàn)結(jié)果導(dǎo)入到Excel中，按照定單平均等待時(shí)間AWT為第一排序字段、在制品庫存WIP為第二排序字段，遞增排列。將排列結(jié)果的前9個(gè)相對(duì)最優(yōu)配置方案及實(shí)驗(yàn)結(jié)果列入表5.9。表5.9 模型仿真參數(shù)優(yōu)化結(jié)果數(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兩個(gè)數(shù)值得出的系統(tǒng)績效相同；再比較WIP數(shù)值，最低值為4794；維修工數(shù)量直到第7條配置時(shí)，才為3人，少于第一條配置。因此，選定第一條配置為系統(tǒng)的最優(yōu)配置方案。在下一章進(jìn)行的仿真實(shí)驗(yàn)研究中，如果沒有特別的說明，模型中的和Qmen將確定為這兩個(gè)數(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