JIT看板生产系统WITNESS仿真建模和优化实现

1、JIT看板生产系统WITNESS真建模和优化实现I<<>>1.1生产流程描述实例系统生产制造单一类型的产品一一某部件，其生产原材料为各种类型的钢板，通过机械加工，然后组装成一个空轮（没有安装轴承座、轴承和轴）。组成某部件的部件主要为5类，分别为零件1、零件2、零件3、零件4和零件5,以数字1,2,3,4,5表示；其生产过程包括12个主要加工单元，分别为剪板切割（WS1、打坡口（WS2、卷圆（WS3、轮圈焊接（WS4、找圆（WS5、车断面（WS6、数控切割（WS7、调平（WS8、油漆（WS9、钻孔（WS10、内轮焊接（分装工作站）和轮子焊接（总装工作站）。零件1原材料为特

2、定型号的钢板，加工过程依次为两块零件1钢板通过卷板切割机切割成适合大小，通过坡口机将钢板两端结合处内外打出坡度，通过卷圆机将钢板卷成轮圈，经过轮圈焊接工段将轮圈接口处焊接起来，通过找圆机将轮圈找圆，然后进入轮子焊接工段，与内轮和零件5焊接成轮子。零件2原材料也为特定型号的钢板，加工过程为每次三块钢板通过卷板切割机切割成适合大小的板材，然后打坡口、卷圆、焊接、找圆，同零件1加工过程一样。在经过找圆工段之后，内轮圈再经过数控车床，进行端面对车，经过钻孔工段钻出工艺孔，到内轮焊接工段与零件3和零件4焊接成内轮。零件3、零件4和零件5原材料也为特定型号的钢板，首先经过数控切割机切割成型，然后调平、油漆

3、，零件4和零件5需要经过钻口工段，钻制工艺孔，然后，零件4到内轮焊接工段进行与零件2和零件3的焊接，零件5到轮子焊接工段与内轮和零件1焊接成轮子。1.2基本生产单元的分解该生产/库存系统包括四条申行线，分别为:(1) 零件1加工申行线：剪板切割、打坡口、卷圆、轮圈焊接、找圆；零件2加工申行线：剪板切割、打坡口、卷圆、轮圈焊接、找圆、车端面、钻孔；零件3加工申行线：数控切割、调平、油漆；零件4、零件5加工申行线：数控切割、调平、油漆、钻孔。该生产/库存系统包括两个并行加工模块，分别为：(1) 零件2、零件3、零件4焊接为内轮的焊接工段，即分装工作站WS11零件1、内轮、零件5焊接为轮子的轮子焊接

4、工段，即总装工作站WS12生产流程中有资源共享和竞争的作业，如：(1)内、零件1在进行切割、打坡口、卷圆、焊接等工序时，使用的都是相同的设备；(2)零件3、零件4、零件5在进行切割、调平工序时，也是竞争使用相同的设备。1.3系统绩效指标设计系统绩效指标考虑系统产出量、在制品库存、设备稼动率和订单满意度四项。在进行Witness仿真程序设计的时候，将通过如下的数据模型和相关的程序进行统计计算获得。(1)系统产出量(TPUTThroughPU:统计在特定时间长度内，生产/库存系统所产出的成品数量。该指标在Witness中可以直接在成品缓冲区的“actiononinput”中进行变量的累计。每当一个

5、成品进入缓冲区，执行式(5-1)的计算程序，即产出量增加了一个。TPUT=TPUT+1(5-1)在仿真模型中，该统计任务由订单队歹0DMD.b_order的“actiononinput”和函数UF.tput来实现，具体程序参看附录B部分。(2)在制品库存(WR：用来统计某一时点上，生产/库存系统中所拥有的原材料、半成品或成品的量，一般用价值来衡量。由丁生产线上具有多种物料，如零件1、零件2、零件3、零件5、零件4,每种类型物料所具有的价值不一样，需要统计各种物料的实际价值。由丁本实例并不是对实际的问题进行研究，仅仅是提供一种研究方法，所以为了便丁对各种方案进行经济性比较，根据每件产品或部件的价

6、值，给予它们附上对应的价值权重，然后进行统计计算。假设零件3最便宜，以它为1,每件物料经过一道工序加工后，每加工10分钟，价值增加一个单位，不足10分钟的按10分钟计，得出实例系统中物料的价值权重系数如表5.1,表中下划线前的字符为物料名称中文首字母。表5.1物料价值权重系数表*名称Jqj_WS7Fkb_WSNgb_WS7qj_WS8Fkb_WS；Ngb_WSJqj_WS9fkb_WS9ngb_WS122244355，称Wq_WS；Nq_WS1Wq_WS2Nq_WS：2Wq_WS3Nq_WS3Wq_WSiNq_WS4Wq_WS仪虫303035343837444247名称Nq_WS5Jq_WSN

7、q_WS1Fkb_WS0Ngb_WS，0WS11WS12仪虫4449528883151系统中的在制品库存WIP计算公式如式(5-2):(5-2)式中：第i种物料在各工序后的缓冲区中的数量；第i种物料在各工序机器上等待加工或正在加工的数量；第i种物料在各工序间运输工具上的数量;第i种物料的价值权重。在该系统的仿真模型中，某一仿真时点上的QF、QM、QT都可以通过Witness提供的函数直接得出。该指标的统计由目标函数UF.wip来计算，具体设计参看附录B部分。(2) 定单平均等待时间(AWTAverageWaitingTime):用于统计没能够及时被满足的订单在系统中的平均等待时间，用来揭示随机

8、生产/库存系统对订单的反应灵敏度。如果某一订单在提前期之前得到满足，则等待时间为零；如果某一订单在提前期之后得到满足，则等待时间为当前时间与提前期的差值。该指标数值越低，表示顾客满意度越高；该指标数值越高，表示顾客满意度越低。计算公式如式(5-3)所示。(5-3)式中：n仿真期间订单总数；Ti，第i批订单实际发运时间；Ti第i批订单预定发运时间；WT第i批订单的等待时间；SWT一仿真期间订单总等待时间。在仿真模型中，这两项统计任务统一在订单处理单元中执行。执行程序体参见附录B中DMD.demand_mee®“Actionsonfinish”部分。设备平均稼动率(OROperation

9、Rate):用于统计生产系统中机器利用程度，仿真系统对机器类型的元素提供了标准状态统计函数，可以统计机器的空闲、运行、故障和维修等状态所占时间的白分比。在进行生产系统优化过程中，需要使用的机器时间包括下列几种：(i) 最大操作时间：指的是设备可用的最大时间，若设备本身为厂内自购，而非租借，且可完全由厂内自主使用，则最大操作时间一般为日历时间。(ii) 负荷时间：为机器设备可稼动的时间，是最大操作时间扣除停机损失，停机损失包括计划上的休止时间，如休假、教育训练、保养等。(iii) 稼动时间：是负荷时间减去停机时间，而停机时间包括批次转换、设备异常停止、修理、待料等时间。机器稼动率=稼动时间/负荷

10、时间。(iv) 故障间隔时间：机器每运转多长时间就会出现一次故障，并需要维修。(v) 故障维修时间：机器出现故障，修复所需的时间。因此，本系统的设备平均稼动率计算公式如式(5-4)。(5-4)式中：m系统中工作单元的数量；OR系统中第j个工作单元中机器设备的稼动率；在仿真模型中，该统计任务通过函数来实现，程序设计参看附录B中的UF.OR1函数体。1.4生产线平衡在进行系统仿真之前，需要将每个工作单元中机器设备的设计稼动率加以平衡，这样仿真过程不会因为设备能力不平衡，导致某些工序过快成为瓶颈工序，使前后工序的机器设备因为等待而导致利用不足。为了防止这一问题，将根据该系统的设计能力2500台/年，

11、每年250个工作日来平衡生产线。由于机器对每一部件的加工时间均值确定，只能通过设计每个加工单元中机器的数量来对生产线进行设计，使得生产线能够很好的满足设计的生产能力。对每个工作单元机器数量的计算过程涉及到如下的变量符号,首先对它们加以说明。TWTTotalWorkTime:每个工作单元中机器一年中的负荷时间；QUQuantityUnit:每个工作单元所拥有的机器数量；TQ府一TotalQuantityofWork:每个工作单元所完成的产品套数；TCTimeCycle:每个单元完成一套产品的生产平均周期时间，例如：WS3S备完成的一套产品为分别加工一个零件1、一个零件2,则平均周期时间为完成一个

12、零件1的平均周期时间（30'）+完成一个零件2的平均周期时间（30'）,共60分钟；TNTotalNumofBreakdown：全年机器发生故障的次数，每月1次，均值为12次；RRepairTime:单次机器维修所需时间，均值为半天，4小时，240分钟；BBreakdownRate：故障率，均值为240*12/240/480=2.5%;BIBreakdownInterval:机器故障发生时间间隔；OOperatingRate:机器稼动率，假设为80%通过该产品的BOP得出该系统生产单件最终产品，在各个生产单元上需要的平均周期时间TC如表5.2所示。表5.2生产系统中生产单元TC

13、表单元WS1WS2WS3WS4WS5WS6WS10WS7WS8WS9WS11WS12时间5090601204050130838203045通过系统的设计生产能力和机器的生产周期，同时考虑故障比率较低，在计算机器台数时将其忽略不记，设计计算公式如(5-5)所示。(5-5)式(5-5)中涉及四项变量，其中由丁设备的故障、维修等所占时间比例忽略不记，所以机器全年负荷时间为TWT=25K*480分钟/天=120000分钟；假设稼动率产品套数就是生产线的年设计能力，即TQW=2500从表5.2可以得出这12个工作单元的单套产品生产平均周期时间TC。所以根据式(5-5)可以得出每个工作站的机器数量计算公式

14、,如式(5-6)。(5-6)计算式(5-6)可以得出每个加工单元所需的设备数量，如果结果是小数，取不小丁结果的最小整数。通过计算得出每个单元需要配备的机器数量如表5.3所示。表5.3加工单元设备台数表加工单元WS1WS2WS3WS4WS5WS6WS1WS7WS8WS9WS11WS12设备数量232422411112根据设备故障比率，可以计算出设备故障平均时间间隔，计算公式如式(5-7):BI=T"(OR-BR/TNB=12000(0<0.775/12=7750(5-7)即每台设备加工时间累计达到7750分钟，就会发生故障，需要进行平均时间为240分钟的维修后，才可以继续运转。1

15、.5看板数量的确定生产系统循环过程中看板的数量决定了生产、使用或转让的在制品数量，决定了系统运作效率的高低。建立看板控制系统需要确定所需要的看板卡的数量。对丁本文的单看板生产系统需要确定生产看板的套数，以它来控制加工单元的开始生产或停止生产。得到多数研究文献和实际应用接收的计算方式有如下两种：文献9的最优看板计算公式（5-8）。（5-8）式中：k看板卡套数；D一一段时期所需产品的平均数量；L补充订货的提前期；S安全库存；C-一容器容量。文献10的看板数量计算公式（5-9）。（5-9）式中：R组件利用速度，如15件/分钟；Tp单件生产时间，如：1.5分钟/件；Td箱子的运输和传送时间，如：2.5

16、分钟；C箱子的单位容量，如：45件；X需求速率变化的白分比，如10%使用式（5-8）来确定看板数量，由于该公式没有应对动态需求这一项的变量，有效运用该公式的重要一步就是确定安全库存S的数值，而一旦动态需求的随机特征有所变化，该公式不能相应的进行反馈；式（5-9）考虑到动态需求，在公式中设定了需求速率变化X这一项，但是它没有考虑到安全库存，计算出来的数值通常偏小，引起生产系统较严重的缺货现象。本文研究生产/库存系统在需求和生产服务具有随机波动性时，随机变量的变异系数对系统绩效的影响，所以根据式（5-8）和式（5-9）进行修正，在式（5-9）添加安全系数a项，以保持看板数量能够对随机变量做出反应；

17、同时，考虑批处理作业，添加批处理数据项B,保持系统的有效运行。对于生产看板数量的修正公式如（5-10）式。（5-10）式中：k生产看板数量；取大于表达式数值的最小整数；R组件利用速度；B下道工序对该部件的批处理数量；Tp单件生产时间；Td箱子的运输和传送时间；C箱子的单位容量；X需求速率变化的白分比；a安全系数。同一产品生产线的利用速率R为每天设计产出量除以每天的时间。在仿真实例中，产出量为每天10套，即R=10/480min=1/48（套/分钟）；B和Tp各个部件各不相同，见流程数据表附录A;由于不考虑工序问运输时间对系统绩效的影响，取Td为0,C为1;X就是需求随机分布函数的变异系数CV如

18、果是需求确定性实验，贝UX=0,如果为需求随机性实验，则X有需求的变异系数确定；安全系数a通过仿真优化，暂取为1.7。根据看板数量计算公式（5-9）,对各个部件的各工序生产看板数量计算程序见附录B的InitialActions和看板计算函数UF.KB部分。1.6仿真模型的WITNESS现通过对该随机生产/库存系统流程的分析和系统配置参数的设定，设计出系统的Witness仿真模型，其可视化界面如图5.4所示。图5.4某部件生产/库存系统Witness仿真界面其中仿真模型中变量数据的初始化设置设计丁Witness仿真模型的Initialize程序中,函数的逻辑结构、加工单元、原材料、半成品以及库存

19、等实体元素的输入输出关系等都设计在这些元素对应的action或rule编辑框内，具体程序参考附录部分。1.7模型参数优化在计算机模型验证部分，为了证明计算机模型同理论模型是相符的，只考虑了计算机模型在确定型和随机型两种情况下的设备稼动率和故障率，以及产品产出数量，没有考虑系统的另外两个绩效指标：库存水平和定单平均等待时间。这两项指标与看板数量和维修人员数量之间具有很强的相关性。因此为了研究各种可变因子变化对系统绩效的影响，需要建立一个统一的、优化的模型参数环境，需要对模型的相关参数进行优化取值，涉及的参数有看板安全系数a和维修工的数量，对其进行优化选择的准则是：在制品库存尽可能低、定单平均等待时间尽可能低、维修工数量尽可能少。将参