CN103116324A - 基于指标预测和在线学习的微电子生产线调度方法 - Google Patents
基于指标预测和在线学习的微电子生产线调度方法 Download PDFInfo
- Publication number
- CN103116324A CN103116324A CN2012105437877A CN201210543787A CN103116324A CN 103116324 A CN103116324 A CN 103116324A CN 2012105437877 A CN2012105437877 A CN 2012105437877A CN 201210543787 A CN201210543787 A CN 201210543787A CN 103116324 A CN103116324 A CN 103116324A
- Authority
- CN
- China
- Prior art keywords
- scheduling
- fuzzy
- gene
- operations
- stage
- 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
- 238000000034 method Methods 0.000 title claims abstract description 105
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 39
- 238000004377 microelectronic Methods 0.000 title claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 30
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 9
- 230000007246 mechanism Effects 0.000 claims abstract description 4
- 108090000623 proteins and genes Proteins 0.000 claims description 61
- 238000012549 training Methods 0.000 claims description 43
- 238000012545 processing Methods 0.000 claims description 33
- 238000003754 machining Methods 0.000 claims description 14
- 238000005457 optimization Methods 0.000 claims description 9
- 230000010076 replication Effects 0.000 claims description 9
- YCBOYOYVDOUXLH-UHFFFAOYSA-N 1,2-Diethylhydrazine Chemical compound CCNNCC YCBOYOYVDOUXLH-UHFFFAOYSA-N 0.000 claims description 7
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012163 sequencing technique Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000000717 retained effect Effects 0.000 description 2
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 239000010936 titanium Substances 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
技术领域
本发明属于先进制造、自动化和信息领域,具体涉及具有可重入特征的微电子生产线的调度方法。
背景技术
针对较大规模具有可重入特征的微电子生产线调度问题,目前有如下几类工件排序方法:基于模型的方法、智能优化方法、基于启发式规则的方法。基于模型的排序方法需要对调度问题附加一些较强且不切实际的假设条件,因而该类方法主要用来对一些典型且规模较小的微电子生产线调度问题进行性能分析;由于研究的问题一般具有较大的规模,基于计算智能的方法也难以直接应用;基于启发式规则的排序方法实施简单,但该类方法考虑的因素少,对规模较大的微电子生产线调度问题一般很难获得令人满意的解;
发明内容
为解决以最小化平均流经时间和最大化瓶颈机器利用率为调度目标的较大规模微电子生产线调度问题,本发明提出一种基于指标预测和在线学习的微电子生产线调度方法(简称为DPL)。本发明首先考虑到较大规模的问题特征,所研究的调度问题被动态分解为规模相对较小的当前阶段调度子问题和用于全局指标预测的当前阶段剩余调度子问题,当前阶段剩余调度子问题本质上也是一个具有可重入特征的微电子生产线调度问题,用来获得当前阶段调度子问题较好的全局调度性能及资源冲突程度特征值;每个当前阶段的调度子问题采用提出的自适应差分进化方法(简称为SDEH)进行求解,其中每个个体的复制概率基于预测获得的资源冲突程度进行自适应调整;当前阶段剩余调度子问题通过多模糊规则(简称为MFRs)进行求解以实现全局调度性能和资源冲突程度值的预测。在每一迭代阶段,根据之前几个阶段调度子问题求解过程所获得的相关数据采用所提出的两阶段增量式在线学习方法(简称为TILM)对多模糊规则进行在线调整,以增强多模糊规则的自适应能力。将本发明应用于以最小化平均流经时间和最大化瓶颈机器利用率为调度目标的微电子生产线可产生较好的调度效果。
1.基于指标预测和在线学习的微电子生产线调度方法,其特征在于,该方法是针对一类以最小化平均流经时间和最大化瓶颈机器利用率为调度目标的较大规模微电子生产线的一种 调度方法,且所述方法是在计算机上依次按如下步骤实现的:
步骤(1):初始化,设定如下基本变量
设定问题变量:
N:微电子生产线中的工序总数
Cn:第n个工序的机器组总数,n=1,2,…,N
μnl:第n个工序第l个机器组包含的机器总数,n=1,2,…,N;l=1,2,…,Cn
L:产品类型总数
Vl:属于第l类产品的lot总数,l=1,2,…,L
Jq:第q个lot,q=1,2,…,Q
Gq:lot Jq所包含的操作总数
Uqg:操作Oqg的可加工机器组集合
Tqg:操作Oqg的可加工机器组总数,即Tqg=||Uqg||
bnc:第n个工序第c个机器组缓存
设定算法相关参数:
W:用于迭代式分解的时间窗口长度
κ:每个当前阶段调度子问题需固定加工开始时间和加工机器的操作与当前阶段调度子问题包含操作的比例
M:用于多模糊规则MFRs学习的训练数据规模
M1:每次进行多模糊规则在线学习时需增加和移除的数据规模
minFS:模糊规则的最小支撑度
步骤(2):采集包括上述工序数、品种数、各Lot的工艺路径和所包含的lot数、各工序的机器组数、各机器组中的机器个数、各操作在各机器组上的加工时间在内的调度相关信息并存储至调度数据库中,并形成待求解的微电子生产线调度问题实例。
步骤(3):微电子生产线调度问题迭代分解
在当前调度时刻,基于时间窗口长度W及各lot的工艺路径按如下方法将整个待调度问题的操作分解为已调度操作集B、当前阶段调度子问题的操作集H和当前阶段剩余调度子问题对应的操作集R。
步骤(3.1):确定已调度操作集B
所有已被确定加工机器及加工开始时间的操作构成操作集B。
步骤(3.2):确定当前阶段调度子问题的操作集H
设hq为lotJq在操作集B中的最大操作号,即
hq=max{g|Oqg∈B}q=1,2,…,Q
操作集H根据下式确定
其中,
步骤(3.3):确定用于预测的当前阶段剩余调度子问题对应的操作集R
操作集R由调度问题中不属于操作集B或H的操作组成,即
步骤(4):采用基于资源冲突程度特征值的自适应差分进化方法SDEH求解当前阶段多目标调度子问题
针对上述微电子生产线调度问题,采用如下SDEH算法求解当前阶段多目标调度子问题,形成当前阶段全局调度解;
步骤(4.1):个体编码和种群初始化
个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
采用启发式规则:最短剩余加工时间优先(Shortest Remaining Processing Time,SRPT)、最先进入队列的优先(First Come First Served,FCFS),最短加工时间优先(Shortest ProcessingTime,SPT)、最长剩余加工时间优先(Maximum Remaining Processing Time,MRPT)及上述规则的随机加权规则,求解当前阶段调度子问题,根据下述方法获得负载较大的前几个机器,并形成由Z个个体组成的初始种群P:
设Ψ为预先给定的个体基因长度的最大值,为当前迭代阶段第n个工序第c个机器组缓存bnc中等待加工操作的总数,是当前调度子问题中在第n个工序第c个机器组上加工操作的总数;根据的大小对bnc进行排序,,设排列结果为bη1,bη2,…,bηΓ,其中,Γ为缓存的总数;
然后,令
则个体由当前调度子问题中在与缓存bη1,bη2,…,bηΓ对应的机器组上加工的操作串组成;每个基因串中的基因表示对应操作在相应机器组上的加工优先级,越靠前其加工优先级越高。
步骤(4.2)初始化
设总迭代步数为Λ,令θ=1,初始复制概率为Rep。
步骤(4.3):交叉
采用顺序交叉方法对随机配对的Z/2个体对进行交叉操作,形成子代种群P1,其中,对配对的两个个体的对应基因串单独进行顺序交叉。
步骤(4.4):变异
基于变异概率λ,对需变异的父代个体基因串S,从当前种群随机选择3个不同个体,与S对应的基因串分别设为S1、S2和S3,产生随机整数r∈{1,2,...,m},m为基因串S的基因长度;然后,根据下述流程生成子代个体对应的基因串S′:
步骤(4.4.1):生成取值在区间[0,1]内的m个随机数。
步骤(4.4.2):若上述第j(j=1,2,...,m)个随机数不小于基因复制概率Repij,且j≠r,则复制基因串S中的第j个基因作为子代个体相应基因串S′的第j个基因,设被复制基因的总数为Υ,将所有Υ个被复制基因从S删除,同时删除S1、S2和S3中对应位置的基因,那么基因串S、S1、S2和S3的长度均为m-Υ;Repij为对应基因的复制概率,根据下式自适应调整:
Repij=Rep×βij
其中,βij为相应基因通过预测获得的资源冲突程度值,其定义为:设与基因xij对应的操作为Ouv,Δ为通过求解前一迭代阶段剩余调度子问题获得的调度问题解的总数,其中在δuv个解中,操作Ouv与其它操作间存在资源冲突,则
βij=δuv/Δ
步骤(4.4.3):设迭代次数i=1。
步骤(4.4.6):若i<m-Υ,则令i=i+1,执行步骤(4.4.4);否则,执行步骤(4.5)。
步骤(4.5):轮盘赌选择,其中,对当前调度子问题中的每个个体,基于预测机制,针对当前阶段剩余调度子问题,采用所提出的MFRs对当前阶段剩余调度子问题求解,获得每个个体对应的全局调度目标函数值。
步骤(4.6):迭代终止条件判别
若θ≤Λ,则执行第5.3步,否则执行步骤(4.7);
步骤(4.7):构成当前阶段的全局调度解
由第四步所得到的当前阶段调度子问题的优化解F1和当前阶段剩余调度子问题相应的优化解F2共同构成当前阶段全局调度解F=F1∪F2。
第五步:固定当前阶段调度子问题中部分操作的加工机器和加工开始时间
步骤(5.1):若当前阶段全局调度解F优于迄今已获得的全局优化解Fo,则令Fo=F。
步骤(5.2):根据迄今已获得的全局优化解Fo和预先给定的待固定操作的比例κ,基于时间轴依次固定当前阶段调度子问题中部分操作的加工机器和加工开始时间。
步骤(6):采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习
在当前阶段调度子问题求解完成后,按下述方法利用相关调度数据采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习。
步骤(6.1):生成训练数据
按下述方法,根据当前阶段调度子问题具有最好全局调度性能的解F1生成用于MFRs自适应调整的M1个数据对,若当前阶段为第一阶段,则生成M个数据对;第m个数据对具有如下形式:
步骤(6.1.2):对上述M个决策时刻两个操作的每个个属性值按如下方法进行归一化处理,获得用于MFRs自适应调整的数据对:
m=1,2,…,M;k=1,2,…,K;j=1,2
步骤(6.2):若为第一阶段,则根据生成的M个数据对生成有效的MFRs;否则根据下述提出的调整方法TILM对用于求解剩余调度子问题的MFRs进行在线调整,其中,MFRs中每条规则具有如下形式:
…,
上述模糊规则具有多个前件和一个后件,是k个条件属性,k=1,2,…,K,为的第ikp个语言变量值,Yq是属性的语言变量值总数;当p=1,对应于需要排序的两个操作中的第一个操作;当p=2,对应于第二个;是由同一属性两个操作的不同取值对应的2维模糊网格。D是结论属性,表示两个操作那个具有更高的调度优先级,从而,分类标签只有两个取值,iR=1表示第一个操作具有更高的调度优先级,iR=2表示第二个操作具有更高的调度优先级;CF(R)为规则R的确信度。
步骤(6.3):初始化
计算min FSa
min FSa=(min FS×M-M1)/M
步骤(6.4):数据移除
按下述方法更新保留的模糊网格。
步骤(6.4.1):从M个训练数据中移除最先生成的M1个训练数据。
步骤(6.4.2):对由M个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格, 计算其在剩余M-M1个训练数据下的模糊支撑度FSr。
步骤(6.4.3):若FSr≥minFSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.5):数据新增
按下述方法更新保留的模糊网格。
步骤(6.5.1):将求解当前阶段调度子问题过程中新生成的M1个训练数据增加到剩余的训练数据中。
步骤(6.5.2):对由剩余的M-M1个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,计算其在加入新生成数据后的M个训练数据下的模糊支撑度FSa。
步骤(6.5.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.6):增加新的模糊网格
根据新生成的M1个训练数据,增加新的模糊网格。
步骤(6.6.1):根据新生成的M1个训练数据,生成模糊支撑度大于或等于min FS的模糊网格。
步骤(6.6.2):对于生成的模糊网格,基于M个训练数据,按步骤(6.5.2)计算其模糊支撑度FSa。
步骤(6.6.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.7):获得有效的MFRs
将所保留的模糊网格生成对应的多模糊规则MFRs。
步骤(6.7.1):根据上述保留的模糊网格获得有效的MFRs;其中,每条模糊分类规则R的模糊支撑度FS(R)作为其确信度,即:
CF(R)=FS(R)
步骤(6.7.2):考虑到两个操作在输入顺序改变的情形下,其调度优先级应该是不变的,从而,在多模糊规则中,每条规则R均生成一条逆规则R′与其对应,逆规则具有如下形式:
…,
其中,CF(R′)为规则R′的确信度,本发明中令CF(R′)=CF(R)
步骤(7):判定算法是否结束
若待调度问题的所有操作均已调度完成,则结束算法;否则,执行步骤(3)。
根据上述基于指标预测和在线学习的微电子生产线调度方法,本发明做了大量的仿真实验,从仿真结果可看出,本发明对较大规模具有可重入特征的微电子生产线调度问题在保证较高瓶颈机器利用率下显著缩短了平均流经时间。
附图说明
图1:基于指标预测和在线学习的微电子生产线调度硬件系统结构示意图。
图2:基于指标预测和在线学习的微电子生产线调度方法流程示意图。
图3:多模糊规则在线学习方法结构示意图。
图4:传统差分进化方法(DE)和DPL算法的平均流经时间优化曲线。
具体实施方式
本发明调度方法依赖于相关数据采集系统,有调度系统客户端和调度服务器实现。在实际微电子生产线调度中应用本发明的软硬件架构示意图如图1所示,本发明的实施方式如下。步骤(1):采集上述微电子生产线工序数、品种数、各品种的工艺路径和所包含的lot数、各工序的机器组数、各机器组中的机器个数、各操作在各机器组上的加工时间在内的调度相关信息并存储至调度数据库中,并形成待求解的微电子生产线调度问题实例。
步骤(2):设定算法相关参数
W:用于迭代式分解的时间窗口长度;
κ:每个当前阶段调度子问题需固定加工开始时间和加工机器的操作与当前阶段调度子问题包含操作的比例;
M:用于多模糊规则MFRs学习的训练数据规模;
M1:每次进行多模糊规则在线学习时需增加和移除的数据规模;
min FS:模糊规则的最小支撑度。
步骤(3):微电子生产线调度问题迭代分解
在当前调度时刻,基于时间窗口长度W及各lot的工艺路径按如下方法将整个待调度问题的操作分解为已调度操作集B、当前阶段调度子问题的操作集H和当前阶段剩余调度子问题对应的操作集R,具体过程可参加“发明内容”相关部分。
步骤(4):采用基于资源冲突程度特征值的自适应差分进化方法SDEH求解当前阶段多目标调度子问题
针对上述微电子生产线调度问题,采用如下SDEH算法求解当前阶段多目标调度子问题,形成当前阶段全局调度解。
步骤(4.1):个体编码和种群初始化
个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
采用SRPT、FCFS、SPT和MRPT启发式规则及上述规则的随机加权规则,求解当前阶段调度子问题,形成由Z个个体组成的初始种群P,其中,其负载最大的机器由所提出的相 应方法确定。上述规则的含义见“发明内容”个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
步骤(4.2)初始化
设总迭代步数为Λ,令迭代次数θ=1,初始复制概率为Rep。
步骤(4.3):交叉
采用顺序交叉方法对随机配对的Z/2个体对进行交叉操作,形成子代种群P1,其中,对配对的两个个体的对应基因串单独进行顺序交叉。
步骤(4.4):变异
基于变异概率λ,对需变异的父代个体基因串S,从当前种群随机选择3个不同个体,与S对应的基因串分别设为S1、S2和S3,产生随机整数r∈{1,2,...,m},m为基因串S的基因长度;然后,根据所提出方法生成子代个体对应的基因串S′。
步骤(4.5):轮盘赌选择,其中,对当前调度子问题中的每个个体,基于预测机制,针对当前阶段剩余调度子问题,采用所提出的MFRs对其对应的全局调度目标函数值进行预测。
步骤(4.6):迭代终止条件判别
若θ≤Λ,则执行第5.3步,否则执行步骤(4.7)。
步骤(4.7):构成当前阶段的全局调度解
由第四步所得到的当前阶段调度子问题的优化解F1和当前阶段剩余调度子问题相应的优化解F2共同构成当前阶段全局调度解F=F1∪F2;
第五步:固定当前阶段调度子问题中部分操作的加工机器和加工开始时间
步骤(5.1):若当前阶段全局调度解F优于迄今已获得的全局优化解Fo,则令Fo=F。
步骤(5.2):根据迄今已获得的全局优化解Fo和预先给定的待固定操作的比例κ,基于时间轴依次固定当前阶段调度子问题中部分操作的加工机器和加工开始时间。
步骤(6):采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习
在当前阶段调度子问题求解完成后,按下述方法利用相关调度数据采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习。
步骤(6.1):生成训练数据
按“发明内容”中方法,根据当前阶段调度子问题具有最好全局调度性能的解F1生成用于MFRs自适应调整的M1个数据对,若当前阶段为第一阶段,则生成M个数据对。
步骤(6.2):若为第一阶段,则根据生成的M个数据对生成有效的MFRs;否则根据”发明 内容”提出的调整方法TILM对用于求解剩余调度子问题的MFRs进行在线调整,其中,MFRs中每条规则具有“发明内容”中给出的形式。
步骤(6.3):初始化
计算min FSa
min FSa=(minFS×M-M1)/M
步骤(6.4):数据移除
按下述方法更新保留的模糊网格。
步骤(6.4.1):从M个训练数据中移除最先生成的M1个训练数据。
步骤(6.4.2):对由M个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,
计算其在剩余M-M1个训练数据下的模糊支撑度FSr,详见“发明内容”。
步骤(6.4.3):若FSr≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.5):数据新增
按下述方法更新保留的模糊网格。
步骤(6.5.1):将求解当前阶段调度子问题过程中新生成的M1个训练数据增加到剩余的训练数据中。
步骤(6.5.2):对由剩余的M-M1个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,计算其在加入新生成数据后的M个训练数据下的模糊支撑度FSa,详见“发明内容”。
步骤(6.5.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.6):增加新的模糊网格
根据新生成的M1个训练数据,增加新的模糊网格。
步骤(6.6.1):根据新生成的M1个训练数据,生成模糊支撑度大于或等于minFS的模糊网格。
步骤(6.6.2):对于生成的模糊网格,基于M个训练数据,按步骤(6.5.2)计算其模糊支撑度FSa。
步骤(6.6.3):若FSa≥minFSa,则保留相应的模糊网格;否则,移除相应的模糊网格。
步骤(6.7):获得有效的MFRs
将所保留的模糊网格生成对应的多模糊规则MFRs,详见“发明内容”.
步骤(7):判定算法是否结束
若待调度问题的所有操作均已调度完成,则结束算法;否则,执行步骤(3);
基于指标预测和在线学习的微电子生产线调度方法的参数选择如下:
●用于迭代分解过程中的算法参数:
时间窗口长度W:设为所加工产品中所有操作平均加工时间的10倍与最大加工时间的最大值;
每个调度子问题需要固定加工开始时间和加工机器的操作的比例κ=80%。
●用于多模糊规则在线学习的算法参数:
用于多模糊规则MFRs学习的训练数据规模M=200;
每次训练数据增加和移除的规模M1=20;
模糊规则的最小支撑度min FS=0.4;
考虑到上述多模糊规则的学习过程需要较长的计算时间,本算法在每隔5个迭代阶段学习一次。
●用于求解当前阶段调度子问题的自适应差分进化算法参数:
种群规模Z=20;
总迭代步数为Λ=20;
初始复制概率Rep=0.5;
变异概率λ=0.2。
根据上述所提出的基于指标预测和在线学习的微电子生产线调度方法,本发明做了大量的仿真试验,运行的硬件环境为:P4 2.8GHz CPU,1024M RAM,操作系统为Windows、UNIX。
由于篇幅所限,仅列出部分实验结果,表1给出了3个产品在相应机器组上的加工时间,其中,∞表示相应机器组具有足够的加工能力,“-”表示对应产品不需要在相应机器组上加工;表2-4分别给出了3个产品的工艺路径:
表1工艺路径参数
表2产品1的工艺路径
表3产品2的工艺路径
表4产品3的工艺路径
实验一:传统差分进化算法(DE)与DPL算法的比较:表5给出了DE与DPL的数值比较结果。其中,P表示四个不同的问题实例集,每个集合包括包括25个问题实例,4个集合内的实例分别具有40,50,60和70个lot,但不同的实例为上述3种产品的不同组合。Ct表示相应算法的平均计算时间,Obj1和Obj2分别表示平均流经时间和瓶颈机器利用率。Ip表示DPL相对DE的最大平均改进百分比。
针对实验一给出的问题实例,图4给出了DE和DPL算法的平均流经时间优化曲线。从中可以看出,DPL优于DE。
表5DE与DPL性能比较
实验二:DPL与启发式规则的比较;针对与实验二相同的问题实例集,表6给出了DPL与各种启发式规则的性能比较,从中可看出,在平均流经时间上DPL均优于启发式规则。
表6DPL与各种启发式规则的性能比较
实验三:基因复制概率自适应调整对DPL算法的影响:表7给出了DPL与无自适应调整的DPL(简称为DPLO)的性能比较。从中可以看出基因复制概率自适应调整对DPL算法有一定的影响。
表7DPL与DPLO的性能比较
Claims (1)
1.基于指标预测和在线学习的微电子生产线调度方法,其特征在于,该方法是针对一类以最小化平均流经时间和最大化瓶颈机器利用率为调度目标的较大规模微电子生产线的一种调度方法,且所述方法是在计算机上依次按如下步骤实现的:
步骤(1):初始化,设定如下基本变量
设定问题变量:
N:微电子生产线中的工序总数
Cn:第n个工序的机器组总数,n=1,2,…,N
μnl:第n个工序第l个机器组包含的机器总数,n=1,2,…,N;l=1,2,…,Cn
L:产品类型总数
Vl:属于第l类产品的lot总数,l=1,2,…,L
Jq:第q个lot,q=1,2,…,Q
Gq:lot Jq所包含的操作总数
Oqg:lot Jq的第g个操作,g=1,2,…,Gq
Uqg:操作Oqg的可加工机器组集合
Tqg:操作Oqg的可加工机器组总数,即Tqg=||Uqg||
bnc:第n个工序第c个机器组缓存
设定算法相关参数:
W:用于迭代式分解的时间窗口长度
κ:每个当前阶段调度子问题需固定加工开始时间和加工机器的操作与当前阶段调度子问题包含操作的比例
M:用于多模糊规则MFRs学习的训练数据规模
M1:每次进行多模糊规则在线学习时需增加和移除的数据规模
min FS:模糊规则的最小支撑度
步骤(2):采集包括上述工序数、品种数、各Lot的工艺路径和所包含的lot数、各工序的机器组数、各机器组中的机器个数、各操作在各机器组上的加工时间在内的调度相关信息并存储至调度数据库中,并形成待求解的微电子生产线调度问题实例;
步骤(3):微电子生产线调度问题迭代分解
在当前调度时刻,基于时间窗口长度W及各lot的工艺路径按如下方法将整个待调度问题 的操作分解为已调度操作集B、当前阶段调度子问题的操作集H和当前阶段剩余调度子问题对应的操作集R;
步骤(3.1):确定已调度操作集B
所有已被确定加工机器及加工开始时间的操作构成操作集B;
步骤(3.2):确定当前阶段调度子问题的操作集H
设hq为lot Jq在操作集B中的最大操作号,即
hq=max{g|Oqg∈B}q=1,2,…,Q
操作集H根据下式确定
其中,
步骤(3.3):确定用于预测的当前阶段剩余调度子问题对应的操作集R
操作集R由调度问题中不属于操作集B或H的操作组成,即
步骤(4):采用基于资源冲突程度特征值的自适应差分进化方法SDEH求解当前阶段多目标调度子问题
针对上述微电子生产线调度问题,采用如下SDEH算法求解当前阶段多目标调度子问题,形成当前阶段全局调度解;
步骤(4.1):个体编码和种群初始化
个体由在负载较大的前几个机器上加工的若干个操作串组成,每个操作串对应一个基因串,即个体由若干个基因串组成。
采用启发式规则:最短剩余加工时间优先(Shortest Remaining Processing Time,SRPT)、最先进入队列的优先(First Come First Served,FCFS),最短加工时间优先(Shortest ProcessingTime,SPT)、最长剩余加工时间优先(Maximum Remaining Processing Time,MRPT)及上述规则的随机加权规则,求解当前阶段调度子问题,根据下述方法获得负载较大的前几个机器,并形成由Z个个体组成的初始种群P:
设Ψ为预先给定的个体基因长度的最大值,为当前迭代阶段第n个工序第c个机器组缓存bnc中等待加工操作的总数,是当前调度子问题中在第n个工序第c个机器组上加工操作的总数;根据的大小对bnc进行排序,,设排列结果为bη1,bη2,…,bηΓ,其中,Γ为缓存的 总数;
然后,令
则个体由当前调度子问题中在与缓存bη1,bη2,…,bηγ对应的机器组上加工的操作串组成;每个基因串中的基因表示对应操作在相应机器组上的加工优先级,越靠前其加工优先级越高;
步骤(4.2)初始化
设总迭代步数为Λ,令θ=1,初始复制概率为Rep;
步骤(4.3):交叉
采用顺序交叉方法对随机配对的Z/2个体对进行交叉操作,形成子代种群P1,其中,对配对的两个个体的对应基因串单独进行顺序交叉;
步骤(4.4):变异
基于变异概率λ,对需变异的父代个体基因串S,从当前种群随机选择3个不同个体,与S对应的基因串分别设为S1、S2和S3,产生随机整数r∈{1,2,...,m},m为基因串S的基因长度;然后,根据下述流程生成子代个体对应的基因串S′:
步骤(4.4.1):生成取值在区间[0,1]内的m个随机数;
步骤(4.4.2):若上述第j(j=1,2,...,m)个随机数不小于基因复制概率Repij,且j≠r,则复制基因串S中的第j个基因作为子代个体相应基因串S′的第j个基因,设被复制基因的总数为Υ,将所有Υ个被复制基因从S删除,同时删除S1、S2和S3中对应位置的基因,那么基因串S、S1、S2和S3的长度均为m-Υ;Repij为对应基因的复制概率,根据下式自适应调整:
Repij=Rep ×βij
其中,βij为相应基因通过预测获得的资源冲突程度值,其定义为:设与基因xij对应的操作为Ouv,Δ为通过求解前一迭代阶段剩余调度子问题获得的调度问题解的总数,其中在δuv个解中,操作Ouv与其它操作间存在资源冲突,则
βij=δuv/Δ
步骤(4.4.3):设迭代次数i=1;
步骤(4.4.5):将填入到S′最左边的空白位置;
步骤(4.4.6):若i<m-Υ,则令i=i+1,执行步骤(4.4.4);否则,执行步骤(4.5);
步骤(4.5):轮盘赌选择,其中,对当前调度子问题中的每个个体,基于预测机制,针对当前阶段剩余调度子问题,采用所提出的MFRs对当前阶段剩余调度子问题求解,获得每个个体对应的全局调度目标函数值;
步骤(4.6):迭代终止条件判别
若θ≤Λ,则执行第5.3步,否则执行步骤(4.7);
步骤(4.7):构成当前阶段的全局调度解
由第四步所得到的当前阶段调度子问题的优化解F1和当前阶段剩余调度子问题相应的优化解F2共同构成当前阶段全局调度解F=F1∪F2;
第五步:固定当前阶段调度子问题中部分操作的加工机器和加工开始时间
步骤(5.1):若当前阶段全局调度解F优于迄今已获得的全局优化解Fo,则令Fo=F;
步骤(5.2):根据迄今已获得的全局优化解Fo和预先给定的待固定操作的比例κ,基于时间轴依次固定当前阶段调度子问题中部分操作的加工机器和加工开始时间;
步骤(6):采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习
在当前阶段调度子问题求解完成后,按下述方法利用相关调度数据采用两阶段增量式学习方法TILM对多模糊规则MFRs进行在线学习;
步骤(6.1):生成训练数据
按下述方法,根据当前阶段调度子问题具有最好全局调度性能的解F1生成用于MFRs自适应调整的M1个数据对,若当前阶段为第一阶段,则生成M个数据对;第m个数据对具有如下形式:
步骤(6.1.2):对上述M个决策时刻两个操作的每个个属性值按如下方法进行归一化处理,获得用于MFRs自适应调整的数据对:
m=1,2,…,M;k=1,2,…,K;j=1,2
步骤(6.2):若为第一阶段,则根据生成的M个数据对生成有效的MFRs;否则根据下述提出的调整方法TILM对用于求解剩余调度子问题的MFRs进行在线调整,其中,MFRs中每条规则具有如下形式:
…,
上述模糊规则具有多个前件和一个后件,是k个条件属性,k=1,2,…,K,为的第ikp个语言变量值,Yq是属性的语言变量值总数;当p=1,对应于需要排序的两个操作中的第一个操作;当p=2,对应于第二个;是由同一属性两个操作的不同取值对应的2维模糊网格。D是结论属性,表示两个操作那个具有更高的调度优先级,从而,分类标签只有两个取值,iR=1表示第一个操作具有更高的调度优
先级,iR=2表示第二个操作具有更高的调度优先级;CF(R)为规则R的确信度;
步骤(6.3):初始化
计算min FSa
min FSa=(minFS×M-M1)/M
步骤(6.4):数据移除
按下述方法更新保留的模糊网格;
步骤(6.4.1):从M个训练数据中移除最先生成的M1个训练数据:
步骤(6.4.2):对由M个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格, 计算其在剩余M-M1个训练数据下的模糊支撑度FSr;
步骤(6.4.3):若FSr≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格;
步骤(6.5):数据新增
按下述方法更新保留的模糊网格;
步骤(6.5.1):将求解当前阶段调度子问题过程中新生成的M1个训练数据增加到剩余的训练数据中:
步骤(6.5.2):对由剩余的M-M1个训练数据生成的每一个模糊支撑度大于或等于min FSa的模糊网格,计算其在加入新生成数据后的M个训练数据下的模糊支撑度FSa;
步骤(6.5.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格;
步骤(6.6):增加新的模糊网格
根据新生成的M1个训练数据,增加新的模糊网格:
步骤(6.6.1):根据新生成的M1个训练数据,生成模糊支撑度大于或等于min FS的模糊网格:
步骤(6.6.2):对于生成的模糊网格,基于M个训练数据,按步骤(6.5.2)计算其模糊支撑度FSa;
步骤(6.6.3):若FSa≥min FSa,则保留相应的模糊网格;否则,移除相应的模糊网格;
步骤(6.7):获得有效的MFRs
将所保留的模糊网格生成对应的多模糊规则MFRs;
步骤(6.7.1):根据上述保留的模糊网格获得有效的MFRs;其中,每条模糊分类规则R的模糊支撑度FS(R)作为其确信度,即:
CF(R)=FS(R)
步骤(6.7.2):考虑到两个操作在输入顺序改变的情形下,其调度优先级应该是不变的,从而,在多模糊规则中,每条规则R均生成一条逆规则R′与其对应,逆规则具有如下形式:
…,
其中,CF(R′)为规则R′的确信度,本发明中令CF(R′)=CF(R)
步骤(7):判定算法是否结束
若待调度问题的所有操作均已调度完成,则结束算法;否则,执行步骤(3)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210543787.7A CN103116324B (zh) | 2012-12-17 | 2012-12-17 | 基于指标预测和在线学习的微电子生产线调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210543787.7A CN103116324B (zh) | 2012-12-17 | 2012-12-17 | 基于指标预测和在线学习的微电子生产线调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103116324A true CN103116324A (zh) | 2013-05-22 |
CN103116324B CN103116324B (zh) | 2014-06-04 |
Family
ID=48414735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210543787.7A Expired - Fee Related CN103116324B (zh) | 2012-12-17 | 2012-12-17 | 基于指标预测和在线学习的微电子生产线调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103116324B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324175A (zh) * | 2013-06-09 | 2013-09-25 | 广东豪美铝业股份有限公司 | 一种铝型材生产线的管控系统 |
CN104536412A (zh) * | 2014-12-23 | 2015-04-22 | 清华大学 | 基于指标预报和解相似度分析的光刻工序动态调度方法 |
CN104698838A (zh) * | 2014-12-23 | 2015-06-10 | 清华大学 | 基于论域动态划分和学习的模糊调度规则挖掘方法 |
CN108255141A (zh) * | 2018-01-17 | 2018-07-06 | 北京理工大学 | 一种装配调度信息生成方法及系统 |
CN110955197A (zh) * | 2018-09-27 | 2020-04-03 | 财团法人资讯工业策进会 | 产线自动配置装置及其方法 |
WO2021043007A1 (zh) * | 2019-09-06 | 2021-03-11 | 北京北方华创微电子装备有限公司 | 物料加工路径选择方法及装置 |
WO2021068274A1 (zh) * | 2019-10-10 | 2021-04-15 | 清华大学 | 基于操作完工时间快速预测的集成电路生产线调度方法 |
CN117215249A (zh) * | 2023-11-01 | 2023-12-12 | 北京航空航天大学 | 可重入制造系统的动态调控方法及相关产品 |
Citations (7)
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 |
US20030195794A1 (en) * | 2002-04-15 | 2003-10-16 | Yasuda To-Oru | Manufacture managing method |
CN101216710A (zh) * | 2007-12-28 | 2008-07-09 | 东南大学 | 一种由计算机实现的自适应选择动态生产调度控制系统 |
CN101782769A (zh) * | 2010-03-08 | 2010-07-21 | 清华大学 | 一种基于指标补偿的平均流经时间快速预测方法 |
CN101788819A (zh) * | 2010-03-08 | 2010-07-28 | 清华大学 | 大规模生产过程一种基于迭代式分解和流松驰的调度方法 |
CN101794115A (zh) * | 2010-03-08 | 2010-08-04 | 清华大学 | 一种基于规则参数全局协调优化的调度规则智能挖掘方法 |
JP2011034432A (ja) * | 2009-08-04 | 2011-02-17 | Kobe Steel Ltd | スケジューリングシステム及びスケジューリング方法、並びにスケジューリングプログラム |
-
2012
- 2012-12-17 CN CN201210543787.7A patent/CN103116324B/zh not_active Expired - Fee Related
Patent Citations (7)
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 |
US20030195794A1 (en) * | 2002-04-15 | 2003-10-16 | Yasuda To-Oru | Manufacture managing method |
CN101216710A (zh) * | 2007-12-28 | 2008-07-09 | 东南大学 | 一种由计算机实现的自适应选择动态生产调度控制系统 |
JP2011034432A (ja) * | 2009-08-04 | 2011-02-17 | Kobe Steel Ltd | スケジューリングシステム及びスケジューリング方法、並びにスケジューリングプログラム |
CN101782769A (zh) * | 2010-03-08 | 2010-07-21 | 清华大学 | 一种基于指标补偿的平均流经时间快速预测方法 |
CN101788819A (zh) * | 2010-03-08 | 2010-07-28 | 清华大学 | 大规模生产过程一种基于迭代式分解和流松驰的调度方法 |
CN101794115A (zh) * | 2010-03-08 | 2010-08-04 | 清华大学 | 一种基于规则参数全局协调优化的调度规则智能挖掘方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324175A (zh) * | 2013-06-09 | 2013-09-25 | 广东豪美铝业股份有限公司 | 一种铝型材生产线的管控系统 |
CN104536412A (zh) * | 2014-12-23 | 2015-04-22 | 清华大学 | 基于指标预报和解相似度分析的光刻工序动态调度方法 |
CN104698838A (zh) * | 2014-12-23 | 2015-06-10 | 清华大学 | 基于论域动态划分和学习的模糊调度规则挖掘方法 |
CN104698838B (zh) * | 2014-12-23 | 2017-03-29 | 清华大学 | 基于论域动态划分和学习的模糊调度规则挖掘方法 |
CN108255141A (zh) * | 2018-01-17 | 2018-07-06 | 北京理工大学 | 一种装配调度信息生成方法及系统 |
CN108255141B (zh) * | 2018-01-17 | 2019-11-26 | 北京理工大学 | 一种装配调度信息生成方法及系统 |
CN110955197A (zh) * | 2018-09-27 | 2020-04-03 | 财团法人资讯工业策进会 | 产线自动配置装置及其方法 |
WO2021043007A1 (zh) * | 2019-09-06 | 2021-03-11 | 北京北方华创微电子装备有限公司 | 物料加工路径选择方法及装置 |
US11703839B2 (en) | 2019-09-06 | 2023-07-18 | Beijing Naura Microelectronics Equipment Co., Ltd. | Material processing path selection method and device |
WO2021068274A1 (zh) * | 2019-10-10 | 2021-04-15 | 清华大学 | 基于操作完工时间快速预测的集成电路生产线调度方法 |
CN117215249A (zh) * | 2023-11-01 | 2023-12-12 | 北京航空航天大学 | 可重入制造系统的动态调控方法及相关产品 |
CN117215249B (zh) * | 2023-11-01 | 2024-02-13 | 北京航空航天大学 | 可重入制造系统的动态调控方法及相关产品 |
Also Published As
Publication number | Publication date |
---|---|
CN103116324B (zh) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103116324B (zh) | 基于指标预测和在线学习的微电子生产线调度方法 | |
CN104536412B (zh) | 基于指标预报和解相似度分析的光刻工序动态调度方法 | |
Lin et al. | An efficient deep reinforcement learning model for urban traffic control | |
Peng et al. | A multi-start variable neighbourhood descent algorithm for hybrid flowshop rescheduling | |
CN108268638A (zh) | 一种基于Spark框架的生成对抗网络分布式实现方法 | |
CN107862179A (zh) | 一种基于相似性和逻辑矩阵分解的miRNA‑疾病关联关系预测方法 | |
CN109215740A (zh) | 基于Xgboost的全基因组RNA二级结构预测方法 | |
CN112381343B (zh) | 一种基于遗传-骨干粒子群混合算法的柔性作业车间调度方法 | |
CN112987664A (zh) | 一种基于深度强化学习的流水车间调度方法 | |
US20240185025A1 (en) | Flexible Parameter Sharing for Multi-Task Learning | |
Kureichik et al. | Generation of bioinspired search procedures for optimization problems | |
Gu et al. | A discrete particle swarm optimization algorithm with adaptive inertia weight for solving multiobjective flexible job-shop scheduling problem | |
CN115640898A (zh) | 一种基于ddqn算法的大规模柔性作业车间调度方法 | |
Rad et al. | GP-RVM: Genetic programing-based symbolic regression using relevance vector machine | |
Cui et al. | Improved artificial bee colony algorithm with dynamic population composition for optimization problems | |
CN105787088A (zh) | 一种基于分段编码遗传算法的文本信息分类方法 | |
CN107590538A (zh) | 一种基于在线序列学习机的危险源识别方法 | |
CN104698838B (zh) | 基于论域动态划分和学习的模糊调度规则挖掘方法 | |
Ma et al. | A hybrid particle swarm optimization and simulated annealing algorithm for job-shop scheduling | |
Ouerfelli et al. | The genetic algorithm with two point crossover to solve the resource-constrained project scheduling problems | |
Busa-Fekete et al. | Ranking by calibrated AdaBoost | |
Ho et al. | Residual scheduling: A new reinforcement learning approach to solving job shop scheduling problem | |
Ding et al. | Multi-policy deep reinforcement learning for multi-objective multiplicity flexible job shop scheduling | |
CN109074348A (zh) | 用于对输入数据集进行迭代聚类的设备和迭代方法 | |
Li et al. | Exploiting reuse in pipeline-aware hyperparameter tuning |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140604 Termination date: 20181217 |
|
CF01 | Termination of patent right due to non-payment of annual fee |