CN101788819B - 大规模生产过程一种基于迭代式分解和流松弛的调度方法 - Google Patents
大规模生产过程一种基于迭代式分解和流松弛的调度方法 Download PDFInfo
- Publication number
- CN101788819B CN101788819B CN2010101194078A CN201010119407A CN101788819B CN 101788819 B CN101788819 B CN 101788819B CN 2010101194078 A CN2010101194078 A CN 2010101194078A CN 201010119407 A CN201010119407 A CN 201010119407A CN 101788819 B CN101788819 B CN 101788819B
- Authority
- CN
- China
- Prior art keywords
- scheduling
- workpiece
- workpieces
- subproblem
- time
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 118
- 230000008569 process Effects 0.000 title claims abstract description 49
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 16
- 238000011031 large-scale manufacturing process Methods 0.000 title claims abstract description 13
- 238000004519 manufacturing process Methods 0.000 claims abstract description 45
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 19
- 230000007246 mechanism Effects 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims description 54
- 238000004088 simulation Methods 0.000 claims description 22
- 238000003754 machining Methods 0.000 claims description 16
- 238000011156 evaluation Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 7
- 230000002068 genetic effect Effects 0.000 claims description 6
- 230000035772 mutation Effects 0.000 claims description 5
- 239000012634 fragment Substances 0.000 claims description 4
- 238000007514 turning Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 238000012854 evaluation process Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000010187 selection method Methods 0.000 claims description 2
- 238000004377 microelectronic Methods 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 238000004904 shortening Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
生产过程优化调度对制造企业缩短制造周期、提高机器利用率、降低生产成本等具有重要作用。本发明针对在微电子、机械等离散行业广泛存在的一类以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程,公开了一种基于迭代式分解和流松弛的调度方法。该方法采用基于预测机制的迭代式分解算法结构,将原调度问题迭代地分解为多个阶段进行求解,在各求解阶段,首先基于工件聚类和流松弛手段建立全局调度指标预测模型,然后,在基于上述预测模型获得的全局调度指标预测值指导下,进行调度子问题的形成和优化求解。将本发明应用于上述以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程后,可有效缩短制造周期,提高生产效率。
Description
技术领域
本发明属于自动控制、信息技术和先进制造领域,针对以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程,公开了一种基于迭代式分解和流松弛的调度方法。
背景技术
生产过程优化调度是提高制造企业生产管理和控制水平的重要手段,其目的是在满足各类资源约束和工艺约束的条件下,通过合理安排各加工任务(工件)在各机器前的加工顺序,使某项或多项生产指标达到最优。生产过程调度水平的提高对企业缩短制造周期、提高机器利用率、降低生产成本,进而提高企业生产效率、经济效益和市场竞争力均具有重要作用。目前,常见的生产过程调度方法包括启发式方法、运筹学方法、软计算方法和人工智能方法等。但由于实际生产过程调度涉及的工件和操作数量多(工件达数百至上千个,操作达数千至上万个),且生产约束较为复杂,已有方法在实际大规模生产过程调度中的应用效果不够理想。
在实际微电子、机械等行业大型制造企业中,其生产任务(工件)常具有可分类特性,即在一次调度过程中,所有工件可根据工艺路径及相应操作的加工时间的不同划分为若干类,同类工件具有较大的相似性,且同类工件数量较多。同时,最小化制造周期是上述企业常见的调度目标。本发明即是针对该类以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程,提供一种基于迭代式分解和流松弛的调度方法。在上述行业大中型企业以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程中应用本调度方法,可有效缩短制造周期、提高生产效率。
另一方面,目前上述行业大多企业均已实施了MES、ERP等系统,已具备采集订单、工艺、进度、设备等信息的基础条件,这为实施本发明提供了可能。
发明内容
本发明针对以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程,公开了一种基于迭代式分解和流松弛的调度方法。该方法通过迭代式分解机制把大规模调度问题动态分解为多个阶段进行求解(每个阶段对应一个规模较小的调度子问题),并在求解每个阶段对应的调度子问题前,首先采用一种基于MD(Macro Distance)距离的工件聚类算法对所有待调度工件进行聚类,在此基础上,基于流松弛手段构建当前阶段调度子问题的全局调度指标预测模型,其可实现对当前调度子问题解全局性能的快速评价,进而,将上述预测模型应用于当前调度子问题的求解过程中,以提高求解效果。
该调度方法的基本流程如图1所示,具体实现步骤说明如下:
步骤(1):初始化,设定如下基本变量:
工件Ji需经过ni个机器组的加工,其工艺路径Ri记为:
步骤(2):采集包括上述工件个数、机器组个数、各机器组中的机器个数、各工件的工艺路径、各操作的加工时间在内的调度相关信息并存储至调度数据库中;
步骤(3):从上述调度数据库中读取包括工件个数、机器组个数、各机器组中的机器个数、各工件的工艺路径、各操作的加工时间在内的调度相关信息,建立生产过程调度模型(,该问题简称为“原调度问题”,所建立的调度模型可表示为如下形式:
其中,Ci为工件Ji的预计完工时间;sti,k为操作Oi,k的预计开始加工时间;Ai为工件Ji工艺路径上所有相邻操作对的集合;Ek,l为在机器gk,l上待加工的所有操作对的集合;
步骤(4):基于时间分解机制将原调度问题迭代地分解为多个阶段进行求解,并在每个求解阶段,首先基于工件聚类和流松弛手段建立全局调度指标预测模型,该预测模型涉及工件加工进度特征指标预测和全局调度性能指标预测,然后,将该预测模型用于形成和求解当前阶段调度子问题过程中;该过程按如下步骤进行:
步骤(4.1):采用基于Macro Distance距离的K-均值聚类算法,将当前阶段所有待调度工件划分为c个聚类,其中,Macro Distance距离简写为“MD距离”,c为聚类数;该聚类过程采用如下流程进行:
步骤(4.1.1):计算工件距离矩阵(md(Ji,Jj))n×n,其中,md(Ji,Jj)表示工件Ji和Jj之间的距离;md(Ji,Jj)的计算方式如下:
其中,d(Ri,Rj,k)=|posk,i-posk,j|表示机器组Gk在工件Ji和Jj的工艺路径Ri和Rj中的相对位置之差,posk,i表示机器组Gk在工件Ji的工艺路径Ri中所处的相对位置,其计算方法为:
其中,Ind(Ri,Gk)为机器组Gk在工艺路径Ri中所处的绝对位置序号,即若该机器组位于Ri中的第一位,则Ind(Ri,Gk)=1,若第二位,则Ind(Ri,Gk)=2,以此类推;Len(Ri)为工件ji的总操作数;
步骤(4.1.2):设定初始聚类中心点
步骤(4.1.3):将各工件划分至各聚类
步骤(4.1.4):重新计算各聚类的聚类中心点
在将各工件分别划分至各聚类中后,重新计算各聚类中心点;聚类中心点的计算方法为:在该聚类中,若某工件与其它工件的平均MD距离最短,则该工件为该聚类的中心点;
步骤(4.1.5):若所有聚类的聚类中心点未改变,则转步骤(4.1.7);
步骤(4.1.6):k=k+1,转步骤(4.1.3);
步骤(4.1.7):若聚类后某些类中的工件数量小于给定最小值NCmin,则将该类并入最近的其它聚类中,并重新计算新聚类的中心点,其中,以聚类中心点之间的MD距离的大小来衡量各聚类之间的距离的远近;
步骤(4.2):在步骤(4.1)所得到的聚类结果基础上,将各聚类内的所有待调度工件均用其聚类中心点对应的工件取代,然后,建立当前阶段待调度问题的流松弛模型并求其最优解;该步骤的具体流程如下:
步骤(4.2.1):将各聚类内的所有工件均用其所属聚类的中心点对应的工件取代,此时,原调度问题转化为一个具有多类工件、每类工件内部各工件完全相同的调度问题;
步骤(4.2.2):对步骤(4.2.1)所形成的调度问题,基于流松弛假设,建立相应的流松弛模型;在流松弛模型中,各类工件以“流”方式存在,即工件可分解为无限小的“碎片”分别进行加工;在上述假设下,加工过程中由于工件的各“碎片”可能同时处于不同的机器组,因此,由ci个聚类工件组成的i类工件的加工进度状况可用非负实数组描述;其中,xi,k(t)为t时刻位于机器组Gk缓冲区所有前续操作已完成,当前操作尚未开始的工件中i类工件的待加工量,该加工量以相应操作的加工时间之和表示,该值为一非负实数;基于上述假设,所建立的流松弛模型具有如下形式:
其中,为机器组Gk可加工的所有操作的集合;为i类工件在机器组Gk上的平均加工时间;为i类工件在机器组Gk中任一台机器上的加工速度;Ti,k(t)为在时间区间[0,t]内i类工件在机器组Gk中已完成的加工量,用已完成操作的加工时间之和表示;1(·)为指示函数,满足:
步骤(4.2.3):上述流松弛模型所对应的一种最优解为如下形式:
xi,k(t)=0 i=1,2,…,c k=2,3,…,ni
其中,xi=xi,0为第i类工件在零时刻的总负载量,Ωk为机器组Gk中每台机器的平均负载,即:
步骤(4.3):基于上述流松弛模型的最优解,确定预测工件加工进度特征指标所需的相关参数;当前阶段对应的工件加工进度特征指标预测值f(Θ)为如下形式:
其中,Θ表示当前调度子问题的求解策略,Ti,k(tend)表示基于当前调度子问题解对原调度问题进行仿真,仿真起始时间为上一阶段调度子问题对应时间窗口的结束时刻,仿真结束时间为当前阶段调度子问题对应时间窗口的结束时刻tend,i类工件于当前阶段调度子问题对应时间窗口结束时刻tend在机器组Gk上的实际加工完成量,该完成量采用已调度工件加工时间之和表示;工件加工进度特征指标预测值f(Θ)反映了在调度子问题对应的时间窗口长度T内,基于调度子问题解获得的各类工件的加工完成量和流松弛模型最优解对应各类工件理论完成量的偏差;
步骤(4.4):选择全局调度性能指标预测所需的若干调度规则;
步骤(4.5):确定当前调度子问题对应时间窗口的长度,在此基础上采用基于问题特征的遗传算法形成并求解该调度子问题,并在调度子问题中采用基于全局调度指标预测模型得到的工件加工进度特征指标预测值和全局调度性能指标预测值对当前调度子问题解性能进行全局评价;当前阶段调度子问题对应的数学模型可表示如下:
其中,tq为第q个求解阶段开始的时刻;Jq为在时刻tq所有未完成的工件集合;Oq为在第q个求解阶段所有需要给出调度策略的操作集合;为第q个求解阶段属于工件Ji的所有相邻操作对集合;为在机器gk,l上待加工的所有相邻操作对集合;为该调度子问题的优化目标函数,为最小化预测得到的最大完工时间。
求解当前调度子问题的基于问题特征的遗传算法主要环节描述如下:
a)编码
采用基于规则的编码方法,编码形式为:
Ch={r1,r2,…,rk,…,rK}
其中,rk表示机器组Gk所采用的调度规则;
b)译码
按时间推进当前阶段调度子问题的仿真过程,并在每个仿真时刻,处理相应的调度事件;工件在各机器组上进行加工时,加工顺序由该机器组对应的调度规则确定;当存在多个可加工机器空闲时,随机选择一台机器上机;同时,遵循活动调度策略,即在存在等待加工的操作时,不允许相应的可加工机器出现空闲;在本调度方法中,调度子问题对应的时间窗口长度由下式确定:
T=max(σk)/Q,
其中,Q为时间窗口长度系数,其大小反映了该调度子问题对应的时间窗口的长度;
c)初始种群产生
对各机器组,从适用于优化制造周期目标的较好调度规则,包括:SPT、WINQ、NINQ、LRPT、FIFO、MOPNR 6种调度规则中随机选出一条规则作为该机器组对应编码位相应的调度规则;
各调度规则的含义介绍如下:
●SPT:操作加工时间最小者优先;
●WINQ:下道工序缓冲区内等待操作加工时间总和最少的工件优先;
●NINQ:下道工序缓冲区内等待操作个数最少的工件优先;
●LRPT:剩余加工时间最大者优先;
●FIFO:先到的工件优先;
●MOPNR:操作数最多的工件优先;
d)交叉和变异
按概率Pc对随机选定的两个个体采用两点交叉方法进行交叉;在对个体进行变异时,采用单点变异方法按概率Pm进行变异;
e)评价和选择
对调度子问题解,采用如下两项指标进行性能评价:
●f(Θ):工件加工进度特征指标预测值;
1)首先基于当前调度子问题解,对原调度问题进行仿真,仿真起始时间为上一阶段调度子问题对应时间窗口的结束时刻,仿真结束时间为当前阶段调度子问题对应时间窗口的结束时刻tend);
2)在上述仿真过程完成后,分别采用SPT、LRPT、WINQ规则继续对原调度问题进行仿真,仿真起始时间为当前调度子问题对应时间窗口的结束时刻,仿真结束时间为原调度问题所有操作完成加工;
3)将基于上述3条调度规则进行仿真得到的全局调度性能指标中的最小值作为全局调度性能指标预测值;
采用随机联赛选择方法对种群进行选择,联赛规模为2;
则选择值较低的解作为较好解;否则,选择工件加工进度特征指标预测值f(Θ)较低的解作为较好解;其中,φ为目标分级系数,取值在区间[0,1]中,其大小反映了各聚类内工件加工进度特征指标预测值f(Θ)在对调度子问题解的评价过程中所起的重要程度;该值越大,说明f(Θ)在其评价过程中的重要性就越大;为此,取φ=αe-βN(s),其中,N(s)表示当前所形成的调度子问题序号;
f)停止条件
算法迭代次数达到最大迭代次数;
步骤(4.6):在当前调度子问题求解完成后,固定当前调度子问题对应操作的调度策略,然后转步骤(4.1),继续求解下一阶段的调度子问题,直至原调度问题求解完成;
步骤(5):将最终所得到的原调度问题的调度方案下发执行。
在微电子、机械等行业大中型制造企业以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程中应用本调度方法,可有效缩短制造周期、提高生产效率。
附图说明
图1:基于迭代式分解和流松弛的调度方法流程示意图。
图2:全局调度指标预测模型结构示意图。
图3:本调度方法在实际制造企业实施所需的软硬件架构示意图。
具体实施方式
本发明公开的基于迭代式分解和流松弛的调度方法依赖于相关数据采集系统,由调度系统客户端和调度服务器实现。在实际制造企业大规模生产过程调度中应用本发明的软硬件架构示意图如图3所示,本发明的实施方式如下。
步骤(1):采集上述工件个数、机器组个数、各机器组中的机器个数、各工件的工艺路径、各操作的加工时间等调度相关信息并存储至调度数据库中;
步骤(2):从上述数据库中读取工件个数、机器组个数、各机器组中的机器个数、各工件的工艺路径、各操作的加工时间等调度相关信息,建立生产过程调度模型;
步骤(3):基于时间分解机制将上述调度模型迭代地分解为多个阶段进行求解,并在每个阶段,首先基于工件聚类和流松弛手段建立全局调度指标预测模型,该预测模型涉及工件加工进度特征指标预测和全局调度性能指标预测,然后,基于该预测模型,形成并求解调度子问题。具体过程可参见“发明内容”相关部分。
步骤(4):将所得到的调度方案下发执行。
为验证上述基于迭代式分解和流松弛的调度方法在以制造周期为调度目标、工件具有可分类特性的大规模生产过程调度问题上的应用效果,本发明做了大量的仿真试验。在仿真试验中,调度问题数据按如下方式产生:
1.产生15个机器组,各机器组中机器数量为[5,10]间均匀分布的随机整数;
2.产生一个工件,其操作数为[5,15]间的随机整数,各操作的加工时间为[10,50]间均匀分布的随机整数,各操作所在的加工机器组随机指定(但同一工件中不同操作的加工机器组之间不可重复);
3.对步骤2重复c次,产生c个工件;
4.从上述c个工件中随机挑选一个工件,对其工艺路径和各操作的加工时间进行随机扰动(对工艺路径采用交换相邻机器组的方法进行扰动,交换次数占总操作数的比例为η=0.2;对加工时间采用加上一个随机数的方法,扰动后的加工时间pnew=pold×(1+r),r∈[-γ,γ],扰动比例为γ=0.1),形成新的工件;
5.对步骤4重复n-c次,产生n-c个工件。
按上述方法可产生出不同规模的、工件具有可分类特性的调度问题,分别记为P3-1~P3-15。其中,问题规模采用工件数n表示。
基于迭代式分解和流松弛的调度方法所涉及的参数选择如下:
●工件聚类过程参数:
聚类个数:取c=n/50
●调度子问题形成与求解参数:
时间窗口长度系数:Q=20
迭代代数:20
种群规模:10
目标分级系数中:α=0.3,对P3-1~P3-5:β=0.25
对P3-6~P3-10:β=0.2
对P3-11~P3-15:β=0.15
交叉变异概率:pc=0.2,pm=0.05
本发明设计了用于与本方法(简称RCP-ID)进行对比的调度方法,包括:
■SPT/LRPT/WINQ:适于优化制造周期目标的启发式规则
■RBGA:基于规则组合的遗传算法
■RCP-ID(Cmax):仅采用全局调度性能指标预测值,不采用工件加工进度特征指标f(Θ)作为调度子问题优化指标的RCP-ID方法,即在计算调度子问题解的适应度值时,目标分级系数φ=0
本调度方法运行的硬件环境为:P42.8GHz CPU,512M RAM,操作系统为Windows XP,编程语言为Visual C++6.0。
表1 本调度方法与其它方法的求解性能比较
问题P3-1~P3-5:工件数n为200
问题P3-6~P3-10:工件数n为800
问题P3-11~P3-15:工件数n为2000
①:该算法中,所采用的调度规则包括SPT、WINQ、NINQ、LRPT、FIFO、MOPNR共六种调度规则,采用基于规则的编码,各机器组对应一个编码位,其种群规模、初始种群产生方法、交叉变异方法与说明书中调度子问题求解算法相同,对调度目标值指数定标后轮盘赌方式进行选择,但迭代次数通过限制算法运行时间来控制(对P3-1至P3-5算法迭代1小时,P3-6至P3-10算法迭代1.5小时,P3-11至P3-15算法迭代2小时)
②:为本方法与RCP-ID(Cmax)方法相比得到的改进率
③:为本方法与SPT/LRPT/WINQ/RBGA中最好的一组结果相比得到的改进率
表2 本方法与其它方法的运行时间比较 单位:秒
①:由于SPT/LRPT/WINQ三种启发式规则的运行时间很接近,因此不再一一列举。表中SPT/LRPT/WINQ一列为采用三种启发式规则求解原调度问题所需的最长运行时间
从表1和表2可看出,在所列的三种规模的调度问题上,本调度方法的优化效果均优于较好的启发式调度规则及基于规则组合的遗传算法。同时,在大规模调度问题(P3-11~P3-15)上的调度效果更明显优于RBGA及启发式规则(平均改进率为14.2%)。在运行时间上,本方法也在可以接受的范围内,且其性能明显优于迭代2小时的RBGA算法,这表明本方法在具有工件可分类特征的大规模生产过程调度问题上具有很好的优化效果。另外,与RCP-ID(Cmax)方法相比,本调度方法也显示出了较大的优越性。
Claims (1)
1.一种基于迭代式分解和流松弛的调度方法,其特征在于,该方法是针对以最小化制造周期为调度目标、工件具有可分类特征的大规模生产过程的一种调度方法,且所述方法在计算机上是按如下步骤实现的:
步骤(1):初始化,设定如下基本变量:
工件Ji需经过ni个机器组的加工,其工艺路径Ri记为:
所有操作的集合记为O;
步骤(2):采集包括上述工件个数、机器组个数、各机器组中的机器个数、各工件的工艺路径、各操作的加工时间在内的调度相关信息并存储至调度数据库中;
步骤(3):从上述调度数据库中读取包括工件个数、机器组个数、各机器组中的机器个数、各工件的工艺路径、各操作的加工时间在内的调度相关信息,建立生产过程调度模型,该问题简称为“原调度问题”,所建立的调度模型可表示为如下形式:
其中,Ci为工件Ji的预计完工时间;sti,k为操作Oi,k的预计开始加工时间;Ai为工件Ji工艺路径上所有相邻操作对的集合;Ek,l为在机器gk,l上待加工的所有操作对的集合;
步骤(4):基于时间分解机制将原调度问题迭代地分解为多个阶段进行求解,并在每个求解阶段,首先基于工件聚类和流松弛手段建立全局调度指标预测模型,该预测模型涉及工件加工进度特征指标预测和全局调度性能指标预测,然后,将该预测模型用于形成和求解当前阶段调度子问题过程中;该过程按如下步骤进行:
步骤(4.1):采用基于Macro Distance距离的K-均值聚类算法,将当前阶段所有待调度工件划分为c个聚类,其中,Macro Distance距离简写为“MD距离”,c为聚类数;该聚类过程采用如下流程进行:
步骤(4.1.1):计算工件距离矩阵(md(Ji,Jj))n×n,其中,md(Ji,Jj)表示工件Ji和Jj之间的距离;md(Ji,Jj)的计算方式如下:
其中,d(Ri,Rj,k)=|posk,i-posk,j|表示机器组Gk在工件Ji和Jj的工艺路径Ri和Rj中的相对位置之差,posk,i表示机器组Gk在工件Ji的工艺路径Ri中所处的相对位置,其计算方法为:
其中,Ind(Ri,Gk)为机器组Gk在工艺路径Ri中所处的绝对位置序号,即若该机器组位于Ri中的第一位,则Ind(Ri,Gk)=1,若第二位,则Ind(Ri,Gk)=2,以此类推;Len(Ri)为工件Ji的总操作数;
步骤(4.1.2):设定初始聚类中心点
步骤(4.1.3):将各工件划分至各聚类
步骤(4.1.4):重新计算各聚类的聚类中心点
在将各工件分别划分至各聚类中后,重新计算各聚类中心点;聚类中心点的计算方法为:在该聚类中,若某工件与其它工件的平均MD距离最短,则该工件为该聚类的中心点;
步骤(4.1.5):若所有聚类的聚类中心点未改变,则转步骤(4.1.7);
步骤(4.1.6):k=k+1,转步骤(4.1.3);
步骤(4.1.7):若聚类后某些类中的工件数量小于给定最小值NCmin,则将该类并入最近的其它聚类中,并重新计算新聚类的中心点,其中,以聚类中心点之间的MD距离的大小来衡量各聚类之间的距离的远近;
步骤(4.2):在步骤(4.1)所得到的聚类结果基础上,将各聚类内的所有待调度工件均用其聚类中心点对应的工件取代,然后,建立当前阶段待调度问题的流松弛模型并求其最优解;该步骤的具体流程如下:
步骤(4.2.1):将各聚类内的所有工件均用其所属聚类的中心点对应的工件取代,此时,原调度问题转化为一个具有多类工件、每类工件内部各工件完全相同的调度问题;
步骤(4.2.2):对步骤(4.2.1)所形成的调度问题,基于流松弛假设,建立相应的流松弛模型;在流松弛模型中,各类工件以“流”方式存在,即工件可分解为无限小的“碎片”分别进行加工;在上述假设下,加工过程中由于工件的各“碎片”可能同时处于不同的机器组,因此,由ci个聚类工件组成的i类工件的加工进度状况可用非负实数组 描述;其中,xi,k(t)为t时刻位于机器组Gk缓冲区的所有前续操作已完成,当前操作尚未开始的工件中i类工件的待加工量,该加工量以相应操作的加工时间之和表示,该值为一非负实数;基于上述假设,所建立的流松弛模型具有如下形式:
其中, 为机器组Gk可加工的所有操作的集合; 为i类工件在机器组Gk上的平均加工时间; 为i类工件在机器组Gk中任一台机器上的加工速度;Ti,k(t)为在时间区间[0,t]内i类工件在机器组Gk中已完成的加工量,该加工量用已完成操作的加工时间之和表示;1(·)为指示函数,满足:
步骤(4.2.3):上述流松弛模型所对应的一种最优解为如下形式:
xi,k(t)=0 i=1,2,…,c k=2,3,…,ni
其中,Ωk为机器组Gk中每台机器的平均负载,即:
步骤(4.3):基于上述流松弛模型的最优解,确定预测工件加工进度特征指标所需的相关参数;当前阶段对应的工件加工进度特征指标预测值f(Θ)为如下形式:
其中,Θ表示当前调度子问题的求解策略,Ti,k(tend)表示基于当前调度子问题解对原调度问题进行仿真,仿真起始时间为上一阶段调度子问题对应时间窗口的结束时刻,仿真结束时间为当前阶段调度子问题对应时间窗口的结束时刻tend,i类工件于当前阶段调度子问题对应时间窗口结束时刻tend在机器组Gk上的实际加工完成量,该完成量采用已调度工件加工时间之和表示;工件加工进度特征指标预测值f(Θ)反映了在调度子问题对应的时间窗口长度T内,基于调度子问题解获得的各类工件的加工完成量和流松弛模型最优解对应各类工件理论完成量的偏差;
步骤(4.4):选择全局调度性能指标预测所需的若干调度规则;
步骤(4.5):确定当前调度子问题对应时间窗口的长度,在此基础上采用基于问题特征的遗传算法形成并求解该调度子问题,并在调度子问题中采用基于全局调度指标预测模型得到的工件加工进度特征指标预测值和全局调度性能指标预测值对当前调度子问题解性能进行全局评价;当前阶段调度子问题对应的数学模型可表示如下:
其中,tq为第q个求解阶段开始的时刻;Jq为在时刻tq所有未完成的工件集合;Oq为在第q个求解阶段所有需要给出调度策略的操作集合; 为第q个求解阶段属于工件ji的所有相邻操作对集合; 为在机器Gk,l上待加工的所有相邻操作对集合; 为该调度子问题的 优化目标函数,为最小化预测得到的最大完工时间;
求解当前调度子问题的基于问题特征的遗传算法主要环节描述如下:
a)编码
采用基于规则的编码方法,编码形式为:
Ch={r1,r2,…,rk,…,rK}
其中,rk表示机器组Gk所采用的调度规则;
b)译码
按时间推进当前阶段调度子问题的仿真过程,并在每个仿真时刻,处理相应的调度事件;工件在各机器组上进行加工时,加工顺序由该机器组对应的调度规则确定;当存在多个可加工机器空闲时,随机选择一台机器上机;同时,遵循活动调度策略,即在存在等待加工的操作时,不允许相应的可加工机器出现空闲;在本调度方法中,调度子问题对应的时间窗口长度由下式确定:
其中,Q为时间窗口长度系数,其大小反映了该调度子问题对应的时间窗口的长度;
c)初始种群产生
对各机器组,从适用于优化制造周期目标的较好调度规则,包括:SPT、WINQ、NINQ、LRPT、FIFO、MOPNR共6种调度规则中随机选出一条规则作为该机器组对应编码位相应的调度规则;
各调度规则的含义介绍如下:
●SPT:操作加工时间最小者优先;
●WINQ:下道工序缓冲区内等待操作加工时间总和最少的工件优先;
●NINQ:下道工序缓冲区内等待操作个数最少的工件优先;
●LRPT:剩余加工时间最大者优先;
●FIFO:先到的工件优先;
●MOPNR:操作数最多的工件优先;
d)交叉和变异
按概率Pc对随机选定的两个个体采用两点交叉方法进行交叉;在对个体进行变异时,采用单点变异方法按概率Pm进行变异;
e)评价和选择
对调度子问题解,采用如下两项指标进行性能评价:
●f(Θ):工件加工进度特征指标预测值;
1)首先基于当前调度子问题解,对原调度问题进行仿真,仿真起始时间为上一阶段调度子问题对应时间窗口的结束时刻,仿真结束时间为当前阶段调度子问题对应时间窗口的结束时刻tend;
2)在上述仿真过程完成后,分别采用SPT、LRPT、WINQ规则继续对原调度问题进行仿真,仿真起始时间为当前调度子问题对应时间窗口的结束时刻,仿真结束时间为原调度问题所有操作完成加工;
3)将基于上述3条调度规则进行仿真得到的全局调度性能指标中的最小值作为全局调度性能指标预测值;
采用随机联赛选择方法对种群进行选择,联赛规模为2;
则选择 值较低的解作为较好解;否则,选择工件加工进度特征指标预测值f(Θ)较低的解作为较好解;其中,φ为目标分级系数,取值在区间[0,1]中,其大小反映了各聚类内工件加工进度特征指标预测值f(Θ)在对调度子问题解的评价过程中所起的重要程度;该值越大,说明f(Θ)在其评价过程中的重要性就越大;为此,取φ=αe-βN(s),其中,N(s)表示当前所形成的调度子问题序号,α=0.3,β=0.2;
f)停止条件
算法迭代次数达到最大迭代次数;
步骤(4.6):在当前调度子问题求解完成后,固定当前调度子问题对应操作的调度策略,然后转步骤(4.1),继续求解下一阶段的调度子问题,直至原调度问题求解完成;
步骤(5):将最终所得到的原调度问题的调度方案下发执行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101194078A CN101788819B (zh) | 2010-03-08 | 2010-03-08 | 大规模生产过程一种基于迭代式分解和流松弛的调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101194078A CN101788819B (zh) | 2010-03-08 | 2010-03-08 | 大规模生产过程一种基于迭代式分解和流松弛的调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101788819A CN101788819A (zh) | 2010-07-28 |
CN101788819B true CN101788819B (zh) | 2011-09-14 |
Family
ID=42532059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101194078A Expired - Fee Related CN101788819B (zh) | 2010-03-08 | 2010-03-08 | 大规模生产过程一种基于迭代式分解和流松弛的调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101788819B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963810A (zh) * | 2010-10-22 | 2011-02-02 | 北京理工大学 | 一种分批流水生产压件数量分析方法 |
CN102521442B (zh) * | 2011-12-06 | 2013-07-24 | 南京航空航天大学 | 基于特征样本的飞机结构件神经网络加工时间预测方法 |
CN102945517B (zh) * | 2012-10-19 | 2016-02-10 | 利诚服装集团股份有限公司 | 一种基于聚类分析的服装标准工时的数据挖掘方法 |
CN103116324B (zh) * | 2012-12-17 | 2014-06-04 | 清华大学 | 基于指标预测和在线学习的微电子生产线调度方法 |
CN103105837B (zh) * | 2012-12-21 | 2015-04-01 | 上海交通大学 | 基于可变时间窗实施两级混合优化批处理调度的方法 |
CN104635772B (zh) * | 2014-12-08 | 2017-02-08 | 南京信息工程大学 | 一种制造系统自适应动态调度方法 |
CN104536412B (zh) * | 2014-12-23 | 2017-02-22 | 清华大学 | 基于指标预报和解相似度分析的光刻工序动态调度方法 |
CN107247686B (zh) * | 2017-05-22 | 2020-11-17 | 电子科技大学 | 一种基于并行算法的fetd仿真模拟方法 |
CN109085803B (zh) * | 2018-07-06 | 2021-07-16 | 昆明理工大学 | 一种dram内存封装过程的优化调度方法 |
TW202013107A (zh) * | 2018-09-27 | 2020-04-01 | 財團法人資訊工業策進會 | 產線自動配置裝置及其方法 |
CN109501110A (zh) * | 2018-12-14 | 2019-03-22 | 四川长虹电器股份有限公司 | 注塑产品自动排产方法 |
CN109858515A (zh) * | 2018-12-24 | 2019-06-07 | 合肥工业大学智能制造技术研究院 | 用于对智能制造的供应链进行订单分批配置的方法及系统 |
CN110673560B (zh) * | 2019-10-10 | 2021-02-05 | 清华大学 | 基于操作完工时间快速预测的集成电路生产线调度方法 |
CN110968057B (zh) * | 2019-11-05 | 2022-01-28 | 天津大学 | 一种解决智能制造车间多目标优化的规划方法 |
CN113570118B (zh) * | 2021-07-06 | 2024-06-18 | 浙江工业大学 | 一种基于调度规则的车间排产和分析方法 |
CN117113109B (zh) * | 2023-07-26 | 2024-02-27 | 常州润来科技有限公司 | 基于机器学习的无缝内螺纹铜管制造工艺确定方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6456996B1 (en) * | 1998-06-05 | 2002-09-24 | I2 Technologies Us, Inc. | Computer implemented scheduling system and process using abstract local search technique |
CN101216710A (zh) * | 2007-12-28 | 2008-07-09 | 东南大学 | 一种由计算机实现的自适应选择动态生产调度控制系统 |
-
2010
- 2010-03-08 CN CN2010101194078A patent/CN101788819B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6456996B1 (en) * | 1998-06-05 | 2002-09-24 | I2 Technologies Us, Inc. | Computer implemented scheduling system and process using abstract local search technique |
CN101216710A (zh) * | 2007-12-28 | 2008-07-09 | 东南大学 | 一种由计算机实现的自适应选择动态生产调度控制系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101788819A (zh) | 2010-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101788819B (zh) | 大规模生产过程一种基于迭代式分解和流松弛的调度方法 | |
Yuan et al. | Research on intelligent workshop resource scheduling method based on improved NSGA-II algorithm | |
Lova et al. | An efficient hybrid genetic algorithm for scheduling projects with resource constraints and multiple execution modes | |
ElMaraghy et al. | Scheduling of manufacturing systems under dual-resource constraints using genetic algorithms | |
Yuan et al. | A co-evolutionary genetic algorithm for the two-machine flow shop group scheduling problem with job-related blocking and transportation times | |
CN103745273B (zh) | 一种半导体制造过程的多性能预测方法 | |
Wen et al. | Dynamic scheduling method for integrated process planning and scheduling problem with machine fault | |
US8983881B2 (en) | Benchmarking progressive systems for solving combinatorial problems | |
CN108694502A (zh) | 一种基于XGBoost算法的机器人制造单元自适应调度方法 | |
Zhang et al. | An integer-coded differential evolution algorithm for simple assembly line balancing problem of type 2 | |
CN113406939A (zh) | 一种基于深度q网络的无关并行机动态混合流水车间调度方法 | |
CN114580678A (zh) | 一种产品维修资源调度方法和系统 | |
CN107831740A (zh) | 一种应用于笔记本零件的分布式生产制造过程中的优化调度方法 | |
Jiang et al. | Energy-conscious flexible job shop scheduling problem considering transportation time and deterioration effect simultaneously | |
Shen et al. | An uncertain programming model for single machine scheduling problem with batch delivery | |
Wang et al. | Order acceptance and scheduling problem with outsourcing in seru production system considering lot-spitting | |
Chaudhry | Job shop scheduling problem with alternative machines using genetic algorithms | |
Gong et al. | Energy-efficient flexible job shop scheduling problem considering discrete operation sequence flexibility | |
Halim et al. | Single-machine integrated production preventive maintenance scheduling: A simheuristic approach | |
Aliahmadi et al. | Flexible flow shop scheduling with forward and reverse flow under uncertainty using the red deer algorithm | |
Lehuang et al. | A multi-objective variable neighborhood search algorithm for precast production scheduling | |
Qin et al. | Linear Disassembly Line Balancing Problem with Tool Deterioration and Solution by Discrete Migratory Bird Optimizer | |
Leißau et al. | Backward-Oriented Decision and Planning Approaches in Production Scenarios: A Systematic Literature Review and Potential Solution Approach | |
Wang et al. | Multi-order scheduling optimisation considering product operation and worker allocation in divisional seru | |
An et al. | A self-adaptive co-evolutionary algorithm for multi-objective flexible job-shop rescheduling problem with multi-phase processing speed selection, condition-based preventive maintenance and dynamic repairman assignment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110914 Termination date: 20120308 |