基于数学规划模型的工件加工问题的求解

1、基于工件加工问题的求解摘要对于一个加工企业而言，如何在最短时间内完成加工任务，是一个企业提高竞争力和利润的关键。本文就是一篇关于工件加工的排序优化问题，在给定的数据和符合实际生产的条件下，合理的安排工件的加工顺序，使总加工时间达到最少。对于工件加工次序模型的求解，我们可以运用许多方法来进行求解，但是考虑到3台机床加工10个零件的给定一加工顺序，所有零件通过机床的顺序是一致的; 每个零件在各机床的加工时间已知，且每台机床在同一时间只能加工一个零件。M2及M3工序上会出现等待。如果采用不同序加工,那么在M1上已加工好的零件,在M2上加工的时间会落到在M1上比其后加工的零件的后面,则其在M2上等待的

2、时间更长,同样在M2与M3工序上也是这样，要求加工时间最短的加工顺序，就必须尽量减少零件在M2及M3工序上的等待时间，由于零件必须在它们要求的时间内完工，即某零件在任务开始起到该零件加工完毕之间所用的总时间应少于该零件的规定完工时间。所以要使各个零件在车间待的总时间最短，其加工零件顺序固然只有一种。那么就要合理选择加工零件的种类及其加工的次序。本题根据已知数据，结合问题中的具体要求，我们引入0/1变量建立零件排序的数学规划模型。通过lingo得出其中的最优排序方案。使得完成这批工件加工任务所需要的总时间最省。然后我们对各个排序后的零件完成特定工序所需花费时间进行求和得到整个加工程序所需总时间。

3、总时间包括了各个零件在机床的加工时间以及加工其它零件的等待时间。最后，根据建立的模型求出某车间加工十个零件所需最短的时间为413分钟，总加工时间最短的加工顺序为D-H-G-I-J-E-A-F-C-B,具体结果如表1-1，1-2。若件加工还要满足下面条件，零件D必须在零件E之前加工；零件H与零件J的加工必须相连；机床M3加工每个零件等待时间不能超过5分钟，总等待时间不能超过30分钟。那么继续利用lingo软件求解可以得出在此条件下最优的顺序为GI-DHJEAFCB，所需最短的时间为425分钟，具体结果如表3-1，3-2。关键词： 线性规划 0/1变量 数学规划模型 lingo软件 一、 问题重述

4、车间上午8:00开始加工十个零件,这些零件必须依次通过机床M1,M2,M3，其加工时间如下表(单位：分钟)。M1M2M3A131520B102018C201615D8106E91413F192014G111612H16918I15127J13791. 试建立模型求出使总加工时间最短的加工顺序。2. 写出各零件加工起止时间表，求出各机床的等待时间。3. 若零件加工还要满足下面条件，零件D必须在零件E之前加工；零件H与零件J的加工必须相连；机床M3加工每个零件等待时间不能超过5分钟，总等待时间不能超过30分钟。试建立模型，重新回答前面两个问题。二、 问题分析零件在M1工序上的总加工时间是固定的，无

5、关乎零件加工顺序。问题关键在于零件在M2以及M3工序上回出现等待。若采取不同加工顺序，那么在M1上已加工好的零件会在M2工序上出现等待。同样在M2上加工好的零件在M3工序上也会出现这种现象。要求加工时间最短的加工顺序，就必须尽量减少零件在M2及M3工序上的等待时间，由于零件必须在它们要求的时间内完工，即某零件在任务开始起到该零件加工完毕之间所用的总时间应少于该零件的规定完工时间。若要使整个加工任务的零件总价值最大，应选择最优加工零件的种类及其加工的次序。若零件D必须在零件E之前加工且零件H与零件J的加工必须相连；机床M3加工每个零件等待时间不能超过5分钟，总等待时间不能超过30分钟。那么就需在

6、第一个建立好的模型上附上条件，得出该条件下的最优次序。三、 模型假设假设一：在后面的模型中，我们都假定了忽略零件在转换工序时的运输时间。即将整个零件加工过程简化为一个连续的过程，只考虑机床在加工零件时其他零件的等待时间。假设二：零件之间是相互独立的,从生产的角度看,先加工一个零件并不响对后面零件的加工。不象有些流水线生产那样,存在固定的加工顺序。假设三：工人都是熟练工人,零件在工序上的加工时间是固定不变的,与工人的操作水平无关。假设四：零件在三个工序上采有同顺序加工,即在工序M1上的加工顺序与在M2及M3工序上的加工顺序相同。在工序M1上的加工是连续不断的。四、 符号说明Xi (1)：i零件在

7、车床M1加工所需时间Xi (2) ：i零件在车床M2加工所需时间Xi (3) ：i零件在车床M3加工所需时间Ti (1):i零件完成在M1加工的总时间Ti (2):i零件完成在M2加工的总时间Ti (3):i零件完成在M3加工的总时间Ti-1(2):(i-1)零件完成在M2加工的总时间(i1)Ti-1(3)(i-1)零件完成在M3加工的总时间(i1)T:加工十个零件模型的总时间五、 模型的建立由问题分析可知工件i在M1工序完成的时间： Ti1=Ti-11+Ti1 (1)对于工件i在M1工序完成的时间Ti1与(i-1)工件完成在M2加工的总时间Ti-12要分两种情况分析:(1)当Ti1Ti-12

8、时，即i工件完成M1工序的总时间大于或等于(i-1)工件完成M2工序的总时间，此时i工件不需要等待(i-1)工件而立即就进入下一工序，因此i工件完成 M2工序的总时间表达式为: Ti2=Ti1+Xi2 (2)当Ti1Ti-12时，即i工件完成M1工序的总时间小于或等于(i-1)工件完成M2工序的总时间，此时i工件需要等待(i-1)工件完成M2工序才能进入M2加工。因此i工件 完成M2工序的总时间表达式为: Ti2=Ti-12+Xi2 综合以上两种情况，得到i工件完成M2工序的总时间计算公式为:-=+Ti2=maxTi1，Ti-12+Xi2 (i1) (2)同理：对于工件i在M2工序完成的时间T

9、i2与(i-1)工件完成在M3加工的总时Ti-13,有i工件完成M3工序的总时间计算公式为: Ti3=maxTi2,Ti-13+Xi3 (i1) (3)综合（1）（2）（3）可得加工十件零件需要的总时间为: T=Ti3+i=210maxTi2,Ti3+Xi3 约束条件：s.t. =1,(i=1,2,3,10) =1,(j=1,2,3,10)六、 模型求解算法流程图通过lingo求解，我们得出以下结论：顺序号零件号M1加工时间XI(1)（分钟）M2加工时间分XI(2)钟）M3加工时间Xi(3 （分钟）1D81062H169183G1116124I151275J13796E914137A13152

10、08F1920149C20161510B102018 顺序号零件号完成M1工序总时间MI(1)（分钟）完成M2工序总时间MI(2（分钟）完成M3工序总时间MI(3)（分钟）1D818242H2433513G3551634I5063705J6370796E7286997A851011218F1041241389C12414015510B134160178 表1-1序顺间时止起号序工M1M2M3D8:00-8:088:08-8:188：18-8:24H8：08-8:248:24-8:338:33-8:51G8:24-8:358:35-8:518:51-9:03I8:35-8:508:51-9:03

11、9:03-9;10J8:50-9:039:03-9:109:10-9:19E9:03-9:129:12-9:269:26-9:39A9:12-9:259:26-9:419:41-10:01F9:25-9:449:44-10:0410:01-10:18C9:44-10:0410:04-10:2410:24-10:35B10:04-10:2610:26-10:4010:40-10:58 表1-2工序号等待时间顺序号M1M2M3D000H067G020I000J000E000A000F040C005B006 顺序号零件号M1加工时间XI(1)（分钟）M2加工时间分XI(2)钟）M3加工时间Xi(3

12、（分钟）1G1116122I151273D81064H16918 5J13796E914137A1315208F1920149C20161510B102018顺序号零件号完成M1工序总时间MI(1)（分钟）完成M2工序总时间MI(2（分钟）完成M3工序总时间MI(3)（分钟）1G1127392I2639463D3249554H4858765J6168856E7084987A83991198F1021221369C12213815310B132158176表3-1序顺间时止起号序工M1M2M3D8:00-8:118:11-8:278：27-8:39H8：11-8:268:27-8:398:39-

13、8:46G8:26-8:328:39-8:498:49-8:55I8:32-8:488:49-8:588:58-9:16J8:48-9:019:01-9:099:16-9:25E9:01-9:109:10-9:249:25-9:38A9:10-9:239:24-9:399:39-9:59F9:23-9:429:42-10:0210:02-10:26C9:42-10:0210:02-10:1810:26-10:33B10:02-10:3210:32-10:3810:38-10:56表3-2工序号等待时间顺序号M1M2M3G000I000D003H003J030E010A001F033C000B0

14、145七、 模型评价 在本题中，对于第一题中构建了数学规划模型，将整批工件的加工任务拆分为在最优的排序下每个工件的实际加工情况来分析，根据各工件在加工过程中加工时间和总时间之间的联系，寻求各工件加工总时间的具体算法。再利用Lingo软件进行求解模型，得出工件的最优排序。其中逻辑严谨，论证充分，算法简洁准确。有效地提高了软件求解效率。考虑到它其实是一个有限源的“单队多服务台”模型由于运用了lingo软件进行求解，使模型本身就具有很大的可移植性，方便了其到m件零件n台机床的延伸与应用。但是我们的模型还是存在一些不足之处。模型对于无限源的情况并不适用，与现实生活并不相符合，只是用于理论的研究。但是总

15、的来说，该模型对于日常生活中的工作安排问题还是具有很大的实际运用价值，值得推广。八、 参考文献1清华大学运筹学教材编写组。运筹学（第三版）M清华大学出版社，20052洪文，吴本忠，Lingo4.0 for windows 优化软件及应用，北京大学出版社，2001.3颜文勇，数学建模，高等教育出版社，北京，2011.64姜启源 谢金星 叶俊等，数学模型（第四版），高等教育出版社，2010。5卢开澄.单目标、多目标与整数规划M.北京:清华大学出版社,1999.6张建中.线性规划,M.北京:科学出版社,19997 /view/63bfcfed9ad51f39

16、1.html九、附录model: !工件先M1后M2再M3的的排序问题;sets: gongjian/g1.g10/:M1_shijian,M2_shijian,M3_shijian; shunxu/s1.s10/:M1_time,M2_time,M3_time,M1_fintime,M2_fintime,M3_fintime; links(shunxu,gongjian): note;endsets !目标函数; min=sum(shunxu(I):M3_fintime(I);!重新排序后各工件的M1床加工时间; for(shunxu(I): che_time(I)=sum(gongjian

17、(J):M1_shijian(J)*note(I,J); );! 重新排序后各工件的M2床加工时间;for(shunxu(I): zuan_time(I)=sum(gongjian(J):M2_shijian(J)*note(I,J); );! 重新排序后各工件的M3床加工时间;for(shunxu(I): xi_time(I)=sum(gongjian(J):M3_shijian(J)*note(I,J); ); !每个顺序位只能有一个工件; for(shunxu(I): sum(gongjian(J): note(I,J)=1; ); !每个工件只能排在一个顺序位上; for(gongji

18、an(J): sum(shunxu(I): note(I,J)=1; );!从新排序后各工件的M1完工时间;for(shunxu(I): che_fintime(I)=sum(shunxu(J)|J#le#I:M1_time(J); );!从新排序后各工件在M2工序的完工时间;for(shunxu(I)| I#gt#1: M2_fintime(I)=if(M1_fintime(I)#ge#M2_fintime(I-1),M1_fintime(I),M2_fintime(I-1)+M1_time(I); ); !顺序1中的M1床工序的完工时间; M2_fintime(1)=M1_fintime(1)+M2_time(1);! 重新排序后各工件在完成M3工序的时间;for(shunxu(I)| I#gt#1: M3_fintime(I)=if(M2_fintime(I)#ge#M3_fintime(I-1),M2_fintime(I),M3_fintime(I-1)+M3_time(I); );!顺序1中的M3工序的完工时间;M3_fintime(1)=M2_fintime(1)+M3_time(1); !定义0/1变量; for(links: