CN114089699A - 一种基于细胞型膜计算模型求解混合流水车间调度的方法 - Google Patents
一种基于细胞型膜计算模型求解混合流水车间调度的方法 Download PDFInfo
- Publication number
- CN114089699A CN114089699A CN202111273590.1A CN202111273590A CN114089699A CN 114089699 A CN114089699 A CN 114089699A CN 202111273590 A CN202111273590 A CN 202111273590A CN 114089699 A CN114089699 A CN 114089699A
- Authority
- CN
- China
- Prior art keywords
- time
- membrane
- processing
- workpiece
- workpieces
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000012528 membrane Substances 0.000 title claims abstract description 114
- 238000000034 method Methods 0.000 title claims abstract description 114
- 238000004364 calculation method Methods 0.000 title claims abstract description 31
- 230000001413 cellular effect Effects 0.000 title claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 68
- 230000008569 process Effects 0.000 claims abstract description 52
- 238000004519 manufacturing process Methods 0.000 claims abstract description 37
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 31
- 238000005457 optimization Methods 0.000 claims abstract description 14
- 238000003754 machining Methods 0.000 claims description 18
- 108091026890 Coding region Proteins 0.000 claims description 13
- 238000003491 array Methods 0.000 claims description 8
- 239000013598 vector Substances 0.000 claims description 8
- 239000010410 layer Substances 0.000 claims description 4
- 239000002344 surface layer Substances 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 229910000831 Steel Inorganic materials 0.000 abstract description 14
- 239000010959 steel Substances 0.000 abstract description 14
- 238000005516 engineering process Methods 0.000 abstract 1
- 238000004891 communication Methods 0.000 description 10
- 210000004027 cell Anatomy 0.000 description 6
- 230000035772 mutation Effects 0.000 description 6
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000005315 distribution function Methods 0.000 description 4
- 230000002068 genetic effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 2
- 210000000170 cell membrane Anatomy 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000002028 premature Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000003850 cellular structure Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000001617 migratory effect Effects 0.000 description 1
- 230000003278 mimic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41865—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by job scheduling, process planning, material flow
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32252—Scheduling production, machining, job shop
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/40—Minimising material used in manufacturing processes
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于细胞型膜计算模型求解混合流水车间调度的方法，涉及钢材加工生产调度技术优化领域，解决现有传统调度算法并行能力不足、算法收敛速度慢、可移植性低、不可用于大规模计算等技术问题，包括n个工件、包含c个工序的流水线上进行加工，各工件不同工序的加工时间已知。每个工件按约定的生产工艺逻辑进行加工。对于每个工序至少存在一台具备执行该工序能力的机器，并且至少有一个工序包含多台并行机器。本发明是基于细胞型膜计算，主要的算法流程除应用选择、交叉、变异规则之外，在各层膜之间应用交流规则，大大提高算法的并行性，搜索可行方案速度更快，同时算法收敛速度提高，可以应用于较大规模的计算的优点。
Description
技术领域
本发明涉及钢材加工生产调度技术优化领域，具体涉及一种基于细胞型膜计算模型求解混合流水车间调度的方法。
背景技术
钢材是钢锭、钢坯等通过压力加工制成所需要的各种形状尺寸和性能的材料，其应用广泛、品种繁多。随着我国钢铁工业的持续发展，钢铁产品逐渐个性化、多样化。而在当前这种离散的生产方式下，产品生产工序多，每道生产工序可采用多台机器进行加工，加工过程中的信息复杂且不易控制，很容易造成经济效益低等问题。如何在生产约束条件下安排大量的钢材加工任务是企业提高生产效率，生产利润的重要手段。现目前，生产调度问题的研究较多，而针对钢材加工生产调度问题的研究相对较少。
许多钢铁行业研究机构、学者与工程师都致力于钢铁生产调度问题的优化。针对钢铁行业的混合流水车间问题，当前主要应用的算法包括动态规划、改进候鸟优化算法、遗传算法、分布式车间流水方法、改进粒子群优化算法等，这些算法存在复杂度高、控制参数过多且复杂、可移植性低、目标优化函数收敛慢等问题。对于由不同钢材生产场景的差异导致的算法参数微调代价过大，不能满足实际生产运用中的各种不确定性的因素。
现在的遗传算法是解决混合流水车间问题的一大重要方法，是一种模拟自然选择和生物进化过程的智能优化算法，算法的主要步骤是求出个体适应度，按照一定的规则进行选择、交叉、变异的过程，最终选择符合要求的的种群，即获得解决问题的较优解，但此算法目前存在着不能及时利用网络的反馈信息，获得目标方案训练需较多时间、算法并行机制能力不做，搜索速度慢、不能很好的解决大规模的计算问题。
发明内容
本发明的目的在于：为了解决上述传统调度算法并行能力不足、算法收敛速度慢、可移植性低、不可用于大规模计算等技术问题，本发明提供了一种基于细胞型膜计算模型求解混合流水车间调度的方法。
本发明采用的技术方案如下：
一种基于细胞型膜计算模型求解混合流水车间调度的方法，包括以下步骤，
步骤1：从车间生产任务表中获取待加工工件的数量n和加工工序c、以及各工件工序的具体加工时间p；
步骤2：获取车间生产线的相关信息，设置生产信息表，包括生产线适应的工序、每工序对应的并行机器数量、生产线的任务序列、已加工的统计时间；
步骤3：设置细胞型膜计算优化算法的相关参数，定义适应度函数；
步骤4：编码操作，即产生各膜的初始对象；已知目前有n个工件待加工，将工件按照1到n顺序编号，编码方式采用工件编号随机全排列，编码代表工件被处理的优先级，当有空闲机器时且满足加工条件时，优先考虑加工编码位置考前的工件；即将随机编码产生的这个序列作为一个向量，为各膜的初始化对象，可以根据模型设定的膜对象数量多次随机编码；
步骤5：解码操作，用于求解初始化编码序列、各膜应用进化规则后的编码序列、各膜应用交流规则编码序列，按照先完工先加工的规则，求出最后一个工件最后一道工序的加工时间，再求出每个膜对象的适应度，最后通过解码求出每个膜对象完成工件的最大加工时间，从而求出每层膜的适应度函数值，用于后续规则；
步骤6：各膜应用进化规则；所有膜都单独应用进化规则，此进化规则包括选择规则、交叉变异规则，达到指定的代数；
步骤7：各膜应用交流规则；所有膜都应用于进化规则后，各膜应用交流规则将膜内最好或者较好的对象送入包含该膜的外层膜，直到达到指定的交流代数；
步骤8：表层输出最好的对象，得到调度顺序，调度结束。
进一步地，所述步骤3具体包括，
步骤3.1：随机指定细胞型膜计算每层膜中包含的膜对象个数，最大进化的代数，最大交流的代数，最大以及最小的变异概率；
步骤3.2：适应度函数定义为Fit＝1/Makespan，Makespan为加工最后一个工序完工时间的最小值，即适应度函数为最大完工时间的倒数值，从定义可以看出，完工时间越少，适应度值越高，该膜对象选为下一代的概率越大。
进一步地，所述步骤5具体包括，
步骤5.1：定义存储该膜对象最大完工时间的变量，即将一个编码序列作为所有待加工工件的加工优先级别，按照先加工先完工规则加工所需最大加工时间totaltime，定义三个与加工工件数量相等的一维数组temp1，temp2，temp3，定义二维数组machine，用于存储第几道工序第几台并行机的统计时间；
数组temp1用于暂时存储当前编码序列，即存储当前的所有待加工工件的加工优先级别，数组temp2用于存储所有工件加工一道工序后的完工时间，temp3用于存储将temp2中各工件加工时间按照从小到大排序的数组；
步骤5.2：循环遍历工序，利用三个数组的依次更新，更新规则为：当要进行下一道加工工序时，temp2中存储的即相当于上一道工序所有工件的完工时间，按照先完工先加工规则，将temp2中完工时间从小到大排序，利用中间数组使得temp1存储的加工顺序进行相应的变化，将完工时间小的工件编码存在前面，同时将排序完后的temp2值赋值给temp3，从而模拟工件的依次加工过程，temp2继续用于存储下一道工序所有工件完工时间，使用二维数组machine记录执行所有工序对应的并行器的加工时间，按照机器开始时间加上加工该工件工序的时间规则累加；
步骤5.3：定义totaltime，赋值为0，用于存储按照该序列加工完成需要的最大完工时间；
步骤5.4：循环遍历所有并行机器，比较最后一道工序机器的累计时间，最大时间就是该流程的加工时间；
步骤5.5：将记录工件和机器加工时间的开始时间、结束时间、机器统计时间的三个数组进行清零，用于记录下一次的膜对象的加工时间统计；
步骤5.6：将总时间赋值给Makespan变量，利用Makespan取倒数，求得该膜对象得适应度值；存储在fit数组的相应位置。
进一步地，所述步骤5.2具体包括，
步骤5.2.1：循环遍历所有工件，比较所有机器的当前工作时间，找出最空闲的机器，便于进行新的工件生产，获得每个工件该工序的完工时间，存储在temp2中；
步骤5.2.2利用中间数组，将temp1数组的内容赋值到中间数组；
步骤5.2.3根据temp2存储的工件完工的时间从小到大进行排序，同时同步变换中间数组的序列；
步骤5.2.4利用中间数组更新temp1的值，将temp2的数据按照从小到大复制到temp3中；
步骤5.2.5继续循环，进行下一道工序。
进一步地，所述步骤5.2.1具体包括，
步骤5.2.1.1记录工序1的第一台并行机器的当前工作时间，与其他执行工序1的机器相互比较，找出时间最小的可执行工序1的并行机器，通过两层循环，也将找出可执行工序1到i的最小时间的并行机器；
步骤5.2.1.2按照顺序提取temp1中的工件号，对工件进行加工；
步骤5.2.1.3定义当前待加工工件的当前工序的开始加工的时间为可执行该工序的并行机器的当前时间值和该工件的上一道工序完工时间的最大值；
步骤5.2.1.4累计开始时间和该工件该道工序所需要花费的时间，将该值赋值给存储工件完工的数组；
步骤5.2.1.5将每个工件的完工时间赋值给数组temp2；
步骤5.2.1.6继续循环，找出用于加工该工序空闲的机器，加工下一个工件。
本发明的有益效果如下：
1.膜计算系统是从类似生物组织、器官或神经网络等细胞群的结构、功能和协作中抽象出来的优化理论计算模型，是一类具有分布式数据处理能力的动态并行计算系统。随着膜计算的逐步发展，越来越多的成熟的算法与之结合，大幅度提高原有算法的性能，应用于更多的生活场景。将其中经典的细胞型膜计算优化算法用于钢铁行业的混合流水车间作业调度问题，控制简单，节约人力物力，能够帮助车间在一定时间范围内确定工件加工顺序的较好安排，获得局部最优解。
2.本调度方法是基于细胞型膜计算，主要的算法流程除应用选择、交叉、变异规则之外，在各层膜之间应用交流规则，大大提高算法的并行性，搜索可行方案速度更快，同时算法收敛速度提高，可以应用于较大规模的计算。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，本说明书附图中的各个部件的比例关系不代表实际选材设计时的比例关系，其仅仅为结构或者位置的示意图，其中：
图1是本发明的流程示意图；
图2是本发明中CMCOA的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明，即所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
下面结合图1、图2，对本发明作详细说明。
实施例1
一种基于细胞型膜计算模型求解混合流水车间调度的方法，步骤1：从车间生产任务表中获取待加工工件的数量n和加工工序c、以及各工件工序的具体加工时间p；步骤2：获取车间生产线的相关信息，设置生产信息表，包括生产线适应的工序、每工序对应的并行机器数量、生产线的任务序列、已加工的统计时间；步骤3：设置细胞型膜计算优化算法的相关参数，定义适应度函数；步骤4：编码操作，即产生各膜的初始对象；已知目前有n个工件待加工，将工件按照1到n顺序编号，编码方式采用工件编号随机全排列，编码代表工件被处理的优先级，当有空闲机器时且满足加工条件时，优先考虑加工编码位置考前的工件；即将随机编码产生的这个序列作为一个向量，为各膜的初始化对象，可以根据模型设定的膜对象数量多次随机编码；步骤5：解码操作，用于求解初始化编码序列、各膜应用进化规则后的编码序列、各膜应用交流规则编码序列，按照先完工先加工的规则，求出最后一个工件最后一道工序的加工时间，再求出每个膜对象的适应度，最后通过解码求出每个膜对象完成工件的最大加工时间，从而求出每层膜的适应度函数值，用于后续规则；步骤6：各膜应用进化规则；所有膜都单独应用进化规则，此进化规则包括选择规则、交叉变异规则，达到指定的代数；步骤7：各膜应用交流规则；所有膜都应用于进化规则后，各膜应用交流规则将膜内最好或者较好的对象送入包含该膜的外层膜，直到达到指定的交流代数；步骤8：表层输出最好的对象，得到调度顺序，调度结束。
本发明的工作原理/工作过程为：考虑在一条流水线上进行生产。n个工件包含c个工序的流水线上进行加工，各工件不同工序的加工时间已知。每个工件按约定的生产工艺逻辑(工序顺序)进行加工。对于每个工序至少存在一台具备执行该工序能力的机器，并且至少有一个工序包含多台并行机器(能够同时执行某一相同工序任务的多台机器)。基于细胞型膜计算优化算法的车间调度方法思想包括：随机产生工件加工顺序的组合，通过该算法先产生对象(即产生工件的加工顺序序列)；根据该算法的选择规则、交叉变异规则以及交流规则共三个顺序计算过程，对工件的加工顺序以及工件工序在机器上的分配情况进行优化确定，使得最大完工时间极小化。
实施例2
进一步地，本发明模型前提需按照一定的规则给所有待加工工件产生多种编码序列，根据编码序列确定工件加工优先级别，设按照级别加工的当前工件有c道工序(工序依次从小到大按数字进行编号)，首先加工工序1，在空闲机器上按照编码中工件优先级别依次加工相应的工件；加工工序2时，按照先完工先加工的原则，在空闲机器上优先加工当前可用工件(已知该工件在上一工序完成加工)。依次类推直到最后一个工件在最后一个工序结束。以上流程将以下列递推过程和目标函数进行详细描述。
本发明的工作原理/工作过程为：符号释义如下：
递推过程如下：
c(i1,1)＝p(1,1) (式1)
c(ik,1)＝c(ik-1,1)+p(k,1) (式2)
c(i1,j)＝c(i1,j-1)+p(1,j) (式3)
c(ik,j)＝max{c(ik-1,j),c(ik,j-1)}+p(k,j) (式4)
式1表示工件1的工序1的加工时间等于该工件在可执行工序1的并行机器i上的完工时间(工件1代表是第一个被加工的工件)。
式2表示工件1的工序k完工的时间等于该工件前一工序(工序k-1)完工时间加上当前工序的加工时间；
式3表示工件j的工序1的完工时间等于同一机器上的前一已加工工件(工件j-1)的工序1的完工时间加上工件j的工序1的加工时间；
式4表示工件j在工序k的完工时间等于工件j的前一工序(工序k-1)的完工时间或同一机器上的前一已加工工件j-1的完工时间中的最大值加上工件j在工序k的加工时间。
目标函数为：Makespan＝min{c(mc,n)} (式5)
式5表示求最后一个工件最后一个工序完工时间的最小值。
1.本发明算法-细胞型膜计算优化算法(CMCOA)
细胞型膜计算优化算法是一种人工智能方法，主要构建类似生物体细胞膜的结构并抽象其内部物质交换等行为模型，最终选出膜内最优的对象。首先说明算法依赖的模型结构及规则。
(1)CMCOA的结构
CMCOA的结构是由各种膜相互嵌套组成类似于生物体中的细胞结构。此结构主要包括表层膜、中间膜和基本膜。表层膜主要与外界环境接触，同时将最优的对象输出到外部环境中；基本膜内部没有任何膜，基本膜之间是相互独立的；中间膜在结构上是位于表层膜和基本膜之间，其内部可由中间膜本身嵌套或者只有基本膜。各层膜可以包含不同的对象数目，同时不同类型膜之间可以通过交流规则交换对象，同类型膜可以通过进化规则选出该膜较优的对象。
CMCOA的结构可以用“[]”嵌套表示，每个“[]”代表一个膜，其具体表示如下：
μ＝[1[2[4]4[5]5[6]6]2[3[7]7[8]8]3]1 (式6)
式6中μ为CMCOA的结构，数字代表膜的编号。此CMCOA的结构μ可以如图2表示：
膜1为表层膜，膜2、3都是中间膜，膜4、5、6、7、8都是基本膜，表层膜1之外为环境，每个膜内部为区域。
(2)CMCOA的对象
已知目前有n个工件待加工，将工件按照1到n顺序编号，编码方式采用工件编号随机全排列，编码代表工件被处理的优先级，当有空闲机器时且满足加工条件时，优先考虑加工编码位置靠前的工件。即将随机编码产生的这个序列作为一个向量，为各膜的初始化对象，可以根据模型设定的膜对象数量多次随机编码。
编码方式注解：例如目前有7个工件，采用随机全排列的方式，可以产生的组合有7！＝20160个，若现在产生的一个序列为1465273，将该序列看作一个向量[1465273]进行后面的运算，那么当有空闲机器且满足加工条件时，优先考虑加工1号工件。根据细胞型膜计算设置的初始化膜对象个数产生相等的全排列序列。
(3)适应度函数
适应度函数为Fit函数，Fit函数表达式为：
Fit＝1/Makespan (式7)
Makespan为加工最后一个工序完工时间的最小值，即适应度函数为最大完工时间的倒数值，从定义可以看出，完工时间越少，适应度值越高，该膜对象选为下一代的概率越大。
(4)CMCOA的进化规则
CMCOA的进化规则包括选择规则、交叉变异规则。选择规则主要是从膜内对象集中选择适应度值较高的对象作为下一代，其主要目的是增强收敛性以及提高计算效率，而交叉变异规则主要在膜内产生新的对象，以提高膜内对象的多样性，避免出现过早收敛的现象。
1)选择规则
该项规则指下一代的个体(对象)按照一定机制从其父代种群(对象集)选取，在父代种群中每个个体(对象)选取的概率与该个体的适应度值有关，即若个体的适应度越大，则此个体被选为下一代种群的概率就越大，反之概率越小。选择规则能够避免遗传基因信息的丢失，并增强全局收敛性以及提高计算效率。CMCOA的选择规则具体表示如下：
FatherGeneration:
f:Prand≤Psj,
式8中okj为膜k中的对象j，f为对象okj被选入下一代的条件，prand为一个[0，1]的随机值，Psj为对象osj被选取的概率值，fitr、fiti(1≤r≤j，1≤i≤Nk)分别为对象okr、对象oki的适应度值，Nk为膜k中对象总数。
2)交叉变异规则
在CMCOA中，交叉变异规则主要是在膜内产生新的个体(对象)。交叉指两个对象对应的向量进行算术交叉，是产生新个体的主要方法，即两个个体通过交叉规则后将产生两个新个体，其主要目的是提高CMCOA中所有膜内种群的多样性。变异指单个对象对应的向量发生数值变化，其主要目的是提高局部搜索最优解的能力，同时避免CMCOA过早陷入局部最优解。由于每个对象都对应一个n(n代表工件总数)维的实数向量，同时为了方便计算，利用线性变换的方法实现CMCOA的交叉变异规则。CMCOA的交叉和变异具体表示如下：
o'ki＝αoki+(1-α)okjo'kj＝αokj+(1-α)okio'kh＝okh+β(1,1,1,1,…,…) (式9)
式9中oki、okj、okj分别为膜k中已有的对象i、对象j、对象h，o′ki、o′kj、o′kh分别为膜k中新生成的对象，fc，fm分别为膜k内个体进行交叉、变异的条件，Pci、pcj分别为膜k中对象oki、okj对象进行交叉的概率，pmh为膜k中对象okh进行变异的概率，Pcrand、pmrand分别为随机生成数并且取值范围为[0，1]，α、β都是给定常数并且取值范围也在[0，1]内。个体的交叉概率和变异概率直接影响到CMCOA的收敛速度，这些概率值可以依据种群(对象集)的多样性和进化代数进行自适应调整，这种调整策略能够避免CMCOA的过早地达到局部最优解，而且改善了CMCOA的后期收敛速度。式9中交叉概率pci、pcj以及变异概率pmh具体计算表达式如下：
式中pcg为膜k中对象oki的交叉概率pci或者对象okj的交叉概率pcj，pmh膜k内对象okh的变异概率，rc为计算膜k内所有对象的交叉概率的参数并且参数rc的取值范围在[pcmin,pcmax]，pcmax、pcmin分别为膜k内对象进行交叉的最大、最小概率并且此两个概率值为给定常数，rm为计算膜k内所有对象的变异概率的参数并且参数rm的取值范围在[pmmin,pmmax]，pmmax、pmmin分别为膜k内对象进行变异的最大、最小概率并且此两个概率值也为给定常数，fitmax为膜k内所有对象的适应度值的最大值，fitavg为膜k内所有对象的适应度值的平均值，fitg为膜k内对象oki或对象okj的适应度值，fith为膜k内对象okh的适应度值，H为膜k内种群(对象集)的信息熵并且其取值范围在[0,log2Nk]，M为在进化代数t时膜k内对象的类别数，相同的对象都归为同一类，t为膜k内当前进化的代数，nd为类别d中对象总数，Nk为膜k内所有对象的总数，Pd为类别d出现的概率并且此概率由类别d的频率值估算，fitl1≤l≤Nk为膜k内对象okl的适应度值。
(5)CMCOA的交流规则
区别于遗传思想，CMCOA除选择、交叉、变异以外，另具备交流规则，具体指将膜内较好的对象送到包含其的外层膜。在CMCOA中，所有膜先各自利用进化规则选出较好对象，然后基本膜将其膜内较好对象送入中间膜(基本膜的外层膜)，接着中间膜将其膜内较好的对象送入表层膜(中间膜的外层膜)，最终表层膜将其膜内最优的对象送到外部环境。特别地，基本膜和外层膜是模拟生物细胞总是选择能量最高的微粒进出细胞膜的过程，而中间膜的规则是根据能量分布函数原理选择较好对象进出此中间膜。交流规则具体表示如下：
式中okj为膜k中的对象j，f1为表层膜和基本膜内的对象进行交流的条件，fitj、fiti分别为对象okj、对象oki的适应度值，NK为膜k内所有对象的总数，此条件说明表层膜和基本膜只将最好的对象(对应的适应度值最大)送入包含其的外部环境或者外层膜，f2为中间膜内的对象进行交流的条件，Prand为随机生成数并且范围在[0，1]内，prj为在中间膜k中对象orj的交流概率，即对象orj被送出的概率值，此概率值参照物理学中能量分布函数进行计算，fitold为上一代膜k内所有对象中最大的适应度值，此条件表明中间膜不仅将最好的对象送入表层膜，同时将该膜内较好的对象也送入表层膜。
在CMCOA的交流规则中，表层膜、基本膜的交流规则比较简单，只需将膜内最大适应度值所对应的对象送出膜外。但是，中间膜的交流规则比较复杂，除了将此膜内最优对象(对应的适应度值最大)送出膜外，还要将较好的对象也送出膜外。这些较好对象主要由三个值来判定，这些值分别为此对象自身的适应度值、上一代该膜内最大的适应度值以及此对象被送出的概率值。在式11中，上述三个值分别对应fitj、fitold、prj，前两个适应度值可以根据式5计算获得，而prj主要参照物理学中能量分布函数进行计算。此能量分布函数具体表示如下：
Pr＝exp((Ef-E)/(kBT)) (式12)
式12中pr为电子占有能量E的概率，Ef为费米能级，KB为Boltzmann常数，T为热力学温度。此外，式12左边能量分布概率当作对象不被送出的概率，同时利用对象的适应度值替换能量值并且适应度越大，能量值就越大，利用当前进化代数替换热力学温度，从而可以得出中间膜k中对象okj的交流概率prj，具体表达式如下：
Prj＝1-exp((fitj-fitold)/g) (式13)
式13中fitj为中间膜k中对象okj的当前适应度值，fitold为上一代膜k内所有对象中最大的适应度值，g为当前的交流代数。在式11中，交流条件f2表明不仅将膜k内最好的对象送出膜外，而且若对象的适应度值大于上一代中最大适应度值，则此对象将被送出膜外，或者若随机概率值Prand小于或等于对象的交流概率，则此对象也将被送到膜外，反之此对象将被束缚在膜内。此交流规则提高了膜内对象的多样性，同时避免CMCOA过早陷入局部最优解。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于细胞型膜计算模型求解混合流水车间调度的方法，其特征在于：包括以下步骤，
步骤1：从车间生产任务表中获取待加工工件的数量n和加工工序c、以及各工件工序的具体加工时间p；
步骤2：获取车间生产线的相关信息，设置生产信息表，包括生产线适应的工序、每工序对应的并行机器数量、生产线的任务序列、已加工的统计时间；
步骤3：设置细胞型膜计算优化算法的相关参数，定义适应度函数；
步骤4：编码操作，即产生各膜的初始对象；已知目前有n个工件待加工，将工件按照1到n顺序编号，编码方式采用工件编号随机全排列，编码代表工件被处理的优先级，当有空闲机器时且满足加工条件时，优先考虑加工编码位置靠前的工件；即将随机编码产生的这个序列作为一个向量，为各膜的初始化对象，可以根据模型设定的膜对象数量多次随机编码；
步骤5：解码操作，用于求解初始化编码序列、各膜应用进化规则后的编码序列、各膜应用交流规则后的编码序列，按照先完工先加工的规则，求出最后一个工件最后一道工序的加工时间，再求出每个膜对象的适应度，最后通过解码求出每个膜对象完成工件的最大加工时间，从而求出每层膜的适应度函数值，用于后续规则；
步骤6：各膜应用进化规则；所有膜都单独应用进化规则，此进化规则包括选择规则、交叉变异规则，达到指定的代数；
步骤7：各膜应用交流规则；所有膜都应用于进化规则后，各膜应用交流规则将膜内最好或者较好的对象送入包含该膜的外层膜，直到达到指定的交流代数；
步骤8：表层输出最好的对象，得到调度顺序，调度结束。
2.根据权利要求1所述的一种基于细胞型膜计算模型求解混合流水车间调度的方法，其特征在于，所述步骤3具体包括，
步骤3.1：随机指定细胞型膜计算每层膜中包含的膜对象个数，最大进化的代数，最大交流的代数，最大以及最小的变异概率；
步骤3.2：适应度函数定义为Fit＝1/Makespan，Makespan为加工最后一个工序完工时间的最小值，即适应度函数为最大完工时间的倒数值，从定义可以看出，完工时间越少，适应度值越高，该膜对象选为下一代的概率越大。
3.根据权利要求1所述的一种基于细胞型膜计算模型求解混合流水车间调度的方法，其特征在于，所述步骤5具体包括，
步骤5.1：定义存储该膜对象最大完工时间的变量，即将一个编码序列作为所有待加工工件的加工优先级别，按照先加工先完工规则加工所需最大加工时间totaltime，定义三个与加工工件数量相等的一维数组temp1，temp2，temp3，定义二维数组machine，用于存储第几道工序第几台并行机的统计时间；
数组temp1用于暂时存储当前编码序列，即存储当前的所有待加工工件的加工优先级别，数组temp2用于存储所有工件加工一道工序后的完工时间，temp3用于存储将temp2中各工件加工时间按照从小到大排序的数组；
步骤5.2：循环遍历工序，利用三个数组的依次更新，更新规则为：当要进行下一道加工工序时，temp2中存储的即相当于上一道工序所有工件的完工时间，按照先完工先加工规则，将temp2中完工时间从小到大排序，利用中间数组使得temp1存储的加工顺序进行相应的变化，将完工时间小的工件编码存在前面，同时将排序完后的temp2值赋值给temp3，从而模拟工件的依次加工过程，temp2继续用于存储下一道工序所有工件完工时间，使用二维数组machine记录执行所有工序对应的并行器的加工时间，按照机器开始时间加上加工该工件工序的时间规则累加；
步骤5.3：定义totaltime，赋值为0，用于存储按照该序列加工完成需要的最大完工时间；
步骤5.4：循环遍历所有并行机器，比较最后一道工序机器的累计时间，最大时间就是该流程的加工时间；
步骤5.5：将记录工件和机器加工时间的开始时间、结束时间、机器统计时间的三个数组进行清零，用于记录下一次的膜对象的加工时间统计；
步骤5.6：将总时间赋值给Makespan变量，利用Makespan取倒数，求得该膜对象得适应度值；存储在fit数组的相应位置。
4.根据权利要求3所述的一种基于细胞型膜计算模型求解混合流水车间调度的方法，其特征在于，所述步骤5.2具体包括，
步骤5.2.1：循环遍历所有工件，比较所有机器的当前工作时间，找出最空闲的机器，便于进行新的工件生产，获得每个工件该工序的完工时间，存储在temp2中；
步骤5.2.2利用中间数组，将temp1数组的内容赋值到中间数组；
步骤5.2.3根据temp2存储的工件完工的时间从小到大进行排序，同时同步变换中间数组的序列；
步骤5.2.4利用中间数组更新temp1的值，将temp2的数据按照从小到大复制到temp3中；
步骤5.2.5继续循环，进行下一道工序。
5.根据权利要求4所述的一种基于细胞型膜计算模型求解混合流水车间调度的方法，其特征在于，所述步骤5.2.1具体包括，
步骤5.2.1.1记录工序1的第一台并行机器的当前工作时间，与其他执行工序1的机器相互比较，找出时间最小的可执行工序1的并行机器，通过两层循环，也将找出可执行工序1到i的最小时间的并行机器；
步骤5.2.1.2按照顺序提取temp1中的工件号，对工件进行加工；
步骤5.2.1.3定义当前待加工工件的当前工序的开始加工的时间为可执行该工序的并行机器的当前时间值和该工件的上一道工序完工时间的最大值；
步骤5.2.1.4累计开始时间和该工件该道工序所需要花费的时间，将该值赋值给存储工件完工的数组；
步骤5.2.1.5将每个工件的完工时间赋值给数组temp2；
步骤5.2.1.6继续循环，找出用于加工该工序空闲的机器，加工下一个工件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111273590.1A CN114089699B (zh) | 2021-10-29 | 2021-10-29 | 一种基于细胞型膜计算模型求解混合流水车间调度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111273590.1A CN114089699B (zh) | 2021-10-29 | 2021-10-29 | 一种基于细胞型膜计算模型求解混合流水车间调度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114089699A true CN114089699A (zh) | 2022-02-25 |
| CN114089699B CN114089699B (zh) | 2024-05-28 |
Family
ID=80298283
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111273590.1A Active CN114089699B (zh) | 2021-10-29 | 2021-10-29 | 一种基于细胞型膜计算模型求解混合流水车间调度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114089699B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114912610A (zh) * | 2022-05-26 | 2022-08-16 | 安徽理工大学 | 一种基于膜计算的单细胞膜最优解计算优化方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106408133A (zh) * | 2016-10-11 | 2017-02-15 | 山东师范大学 | 受生物膜启发dna遗传算法的铁路应急调度优化方法及系统 |
| CN111273562A (zh) * | 2020-01-15 | 2020-06-12 | 安徽理工大学 | 一种基于膜计算实现井下机器人路径优化方法 |
| NL2025034B1 (en) * | 2020-03-03 | 2021-04-20 | Univ Anhui Sci & Technology | Method for realizing underground robot path optimization based on membrane calculation |
-
2021
- 2021-10-29 CN CN202111273590.1A patent/CN114089699B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106408133A (zh) * | 2016-10-11 | 2017-02-15 | 山东师范大学 | 受生物膜启发dna遗传算法的铁路应急调度优化方法及系统 |
| CN111273562A (zh) * | 2020-01-15 | 2020-06-12 | 安徽理工大学 | 一种基于膜计算实现井下机器人路径优化方法 |
| NL2025034B1 (en) * | 2020-03-03 | 2021-04-20 | Univ Anhui Sci & Technology | Method for realizing underground robot path optimization based on membrane calculation |
Non-Patent Citations (1)
| Title |
|---|
| 刘希玉 等: "膜计算研究综述", 《山东师范大学学报（自然科学版）》, vol. 33, no. 2, 30 June 2018 (2018-06-30), pages 127 - 183 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114912610A (zh) * | 2022-05-26 | 2022-08-16 | 安徽理工大学 | 一种基于膜计算的单细胞膜最优解计算优化方法 |
| CN114912610B (zh) * | 2022-05-26 | 2024-04-09 | 安徽理工大学 | 一种基于膜计算的单细胞膜最优解计算优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114089699B (zh) | 2024-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107862411B (zh) | 一种大规模柔性作业车间调度优化方法 | |
| CN104572297B (zh) | 一种基于遗传算法的Hadoop作业调度方法 | |
| CN101901425A (zh) | 一种基于多种群协同进化的柔性作业车间调度方法 | |
| CN106855957A (zh) | 基于相似日和最小二乘支持向量机的工厂母线负荷预测 | |
| CN105740953B (zh) | 一种基于实数编码量子进化算法的不规则排样方法 | |
| CN103729694B (zh) | 基于多色集合层次结构的改进ga求解柔性车间调度的方法 | |
| CN106250583A (zh) | 基于双种群基因表达式编程的动态车间调度规则优化方法 | |
| CN106991442A (zh) | 混合蛙跳算法的自适应核k‑means方法与系统 | |
| CN106610652A (zh) | 改进编码方式的遗传算法解分布式柔性作业车间调度问题 | |
| CN114565239B (zh) | 用于产业园区的综合低碳能源调度方法及系统 | |
| CN107133088A (zh) | 一种基于粒子群算法的多核系统任务调度方法 | |
| CN114089699B (zh) | 一种基于细胞型膜计算模型求解混合流水车间调度的方法 | |
| CN112348323A (zh) | 一种多目标的能源供应与作业柔性排程方法 | |
| CN114021934A (zh) | 基于改进spea2求解车间节能调度问题的方法 | |
| CN104598657B (zh) | 一种基于memetic算法的基因模体重构方法 | |
| Zhang et al. | Furnace-grouping problem modeling and multi-objective optimization for special aluminum | |
| CN113112171B (zh) | 一种基于轮盘赌和遗传算法的批调度方法 | |
| CN112731888B (zh) | 一种批量流混合流水车间调度问题的改进候鸟优化方法 | |
| CN108805463A (zh) | 一种支持削峰型电力需求响应的生产调度方法 | |
| CN111932021A (zh) | 一种再制造系统调度方法 | |
| CN106970840A (zh) | 一种结合任务调度的软硬件划分方法 | |
| Li et al. | Flexible job shop scheduling based on genetically modified neighborhood hybrid algorithm | |
| CN106611215A (zh) | 一种新的布谷鸟搜索算法解决作业车间调度问题 | |
| CN112668864A (zh) | 一种基于狮群算法的车间生产排产方法及系统 | |
| CN115169754B (zh) | 能源调度方法、装置、电子设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |