首页 / 专利库 / 人工智能 / 进化算法 / 遗传算法 / 一种大规模柔性作业车间调度优化方法

一种大规模柔性作业车间调度优化方法

阅读:455发布:2020-05-11

专利汇可以提供一种大规模柔性作业车间调度优化方法专利检索,专利查询,专利分析的服务。并且大规模柔性作业车间调度方法,将大规模的生产任务重组以降低规模,然后利用自适应改进遗传 算法 求解优化。具体步骤如下:(1)首先将加工工艺类似、 工件 尺寸在同一范围内、且毛坯材质相同的工件进行聚类组批,从而减小问题求解规模;(2)设置算法初始参数,采用三层基因编码技术、OBX交叉方式和一定的变异策略,结合仿真实验选择交叉长度,并利用自适应改进 遗传算法 优化求解。该方法可以减小问题求解规模,提高求解速度;并且减小工件完工时间以及拖延期。,下面是一种大规模柔性作业车间调度优化方法专利的具体信息内容。

1.一种大规模柔性作业车间调度优化方法,其特征在于:建立柔性作业车间调度问题的数学模型,基于工件组批调度方法降解问题规模,利用自适应遗传算法优化求解,形成一套求解大规模柔性作业车间调度优化方法,步骤如下:
第一、基于工件特性相似组批
该方法首先将加工工艺类似、工件尺寸在同一范围内且毛坯材质相同的零件进行聚类组批,随机组成不同的批次,每个工件都有各自的交付期,组批之后按照本批当中最早交付期进行调度;
第二、构建大规模柔性作业车间调度问题的数学模型
大规模柔性作业车间调度问题详细描述为,在m台设备(M={Mk|M1,M2...Mm,k=1,2,...m})上加工n件工件(J={Jl|J1,J2,...Jn,l=1,2,...n}),每个工件包含N个事先确定加工顺序的工序,每个工序可以在多台设备上加工;组批之后形成r类工件R={Ri|R1,R2,...Rr,i=1,2,...r},第i类工件的第j道工序在设备k上加工时间为Tijk,Tpi表示第i类工件的第p批次最后一道工序的完工时间,Ds表示工件s的交付期,Dpi表示第i类工件的第p批次的交付期;调度目标是使最大完工时间最小及总拖延期最小;
目标函数:min f=α1f1'+α2f2'    (1)
其中:f1=min(maxTpi)p=1,2...Bi;i=1,2......r    (2)
式(2)f1表示最小最大完工时间;式(3)f2表示总拖延期最小;
约束条件:
式(5)表示第i类工件第p批次的第j道工序在设备k上加工时间等于本批次所有工件加工时间之和
Sip(j-1)k+Tip(j-1)k+Wip(j-1)k≤Sipjk'p=1,2...Bi;i=1,2...r;j=1,2...N;k,k'=1,
2...m    (6)
式(6)表示同一批次同一个工件的后一道工序的开工时间Sipjk′大于等于前道工序的完工时间
Sipjk+Tipjk+Wipjk≤Si'p'j'k p,p'=1,2...Bi;i,i'=1,2...r;j,j'=1,2...N;k=1,
2...m    (7)
式(7)表示占用同一台机器,后一个批次工件的开工时间Si′p′j′k大于等于前一个批次工件的完工时间
式(8)表示每类工件所有批次工件数量总和等于所有工件数总和
Dpi=min Ds i=1,2...r;s=1,2......Cip;p=1,2...Bi     (9)
式(9)表示第i类工件的第p批次的交付期等于本批中最早交付期,具体符号含义如下:
α1表示函数f1的权重;
α2表示函数f2的权重;
f1'表示f1归一化处理后的函数值;
f2'表示f2归一化处理后的函数值;
Sipjk表示第i类工件的第p批次第j道工序在设备k上开始加工时间;
Wipjk表示第i类工件第p批次第j道工序在设备k上的调整时间;
Xipjk表示决策变量,当i类工件的第p批次第j道工序在设备k上加工时,取值为1,否则为
0;
tsjk表示工件Js的第j道工序在设备k上加工时间;
Tipjk表示i类工件的第p批次第j道工序在设备k上加工时间;
Bi表示第i类工件组批形成的批次数;
Cip表示第i类工件第p批次的批量数;
第三、采用三层编码进行遗传算法程序中的基因编码
采用三层编码方式,第一层表示工件组批后形成的批次数量,第二层表示工件工序的加工顺序,第三层表示工件工序加工所选择的机器;
第四、遗传算法参数初始化
按照初始参数选用的原则和依据,设立种群规模为40,交叉概率0.1,变异概率0.04,最大遗传代数为300代;
第五、生成初始种群
按照编码方式,采用随机的方式产生初始种群,并且为了保持初始种群的合法性,在产生初始种群时添加约束,即在工序基因段约束工件号出现的次数应等于工件工序数;
第六、计算种群适应度值
利用公式(1)-(3)计算每个个体的目标函数值,分配适应度值,选择适应度较高进入下一代的交叉变异操作;
第七、选择
选择算子采用轮盘赌注的方法,在此方法中个体被选择的概率和其自身的适应度值相关,适应度值越大被保留下来的概率就越大;
第八、选择OBX交叉算子进行遗传算法交叉
种群经交叉洗牌后,随机选择两个个体,并截取其工序基因交叉,采用基于Order-Based Crossover(OBX)的改进交叉算子进行交叉;并通过仿真探索染色体交叉长度与算法求解精度及运算速度之间的关系,选择合适的染色体交叉长度;
第九、选择变异策略,以使变异后的基因合法化
变异分为批次变异和机器变异,每次批次变异后工序基因和机器基因都需要重新解码形成新的染色体,以保证染色体的合法性;
第十、采用自适应遗传算法进行优化求解
按照公式(10)及式(11)的正弦自适应遗传算法进行优化,以提高算法的收敛速度及收敛精度,以保证算法搜索的全局性和精确性;
交叉概率:
变异概率:
式中fmax为当前种群中最大适应度值;favg为当前种群中平均适应度值;f′为两个待交叉个体中较大的适应度值;f为待变异个体的适应度值;pc1、pm1为系数,在(0,1)区间内取值;
第十一、解码
解码为编码的逆操作,由染色体的工序基因段得出各类工件每个批次的工序加工顺序,机器基因段得出其相应的加工机器,并且根据式(4)-(9)计算出每道工序的开始及结束时间。
2.根据权利要求1所述的一种大规模柔性作业车间调度优化方法,其特征在于:所述步骤第三中的三层编码方式具体为:在有两类工件,每个工件有两道工序情况下,第一层表示工件组批后形成的批次数量,采用整数编码Chrom(1,1)表示第1类工件分为2批;第二层表示工件工序的加工顺序,采用基于工件的编码方式,其中101表示第1类工件的第一个批次,
102表示第1类工件的第二个批次,101第一次出现在工序基因内表示其第一道工序,依次类推,工序基因表示的加工顺序为101→201→102→101→201→102;第三层表示工件工序加工所选择的机器,采用整数编码,Chrom(1,9)表示工件101的第一道工序在可选机器集Mij内选择第一台机器2。
3.根据权利要求1所述的一种大规模柔性作业车间调度优化方法,其特征在于:所述步骤第八中基于OBX的改进交叉算子的交叉步骤如下:在父代P1中随机选择待交叉基因,形成交叉基因池Pos(201,101,102),由于在Chrom2的批次基因中第一类工件批次为1,因此在Chrom2的工序基因中没有102工件,故从交叉基因池中删除102,最终用于交叉的基因有Pos(201,101);在父代P2中找到基因201,由于有多道工序所以在P2中会有多个201,则随机选择一个,然后删除P2中被选中的基因,交叉基因池中其他的基因同理,最后用P1中被选中的基因顺次填补P2中空缺的基因形成子代C2;而P1中则保留备选中的基因以及不能用于交叉的基因,然后用P2中剩余的基因并去除P1中没有的个体202即第二类工件的第二批工件依次填补P1中的空缺形成子代C1。
4.根据权利要求1所述的一种大规模柔性作业车间调度优化方法,其特征在于:所述步骤第八中合适的染色体长度的选择方式具体为:选择染色体长度的10%到100%进行交叉,并各自反复进行仿真运算求其平均值,然后求加权值汇总后,发现交叉长度在占染色体长度为10%时最优;为了找出更优的目标值,又选择染色体长度的10%到20%进行交叉,并反复仿真运算求其平均值,然后求加权值,发现选择染色体长度12%进行交叉可以获得更精确的解和提高运算速度,所以选择12%的交叉长度。
5.根据权利要求1所述的一种大规模柔性作业车间调度优化方法,其特征在于:所述步骤第十中按照正弦自适应遗传算法进行优化时,在优化过程的初始阶段选取较大的交叉率和变异率以获得较快的收敛速度且保持种群的多样性,经多次迭代后,为避免最优解被破坏,此时选取较小的交叉率和变异率进行细致搜索。

说明书全文

一种大规模柔性作业车间调度优化方法

技术领域

[0001] 该发明涉及离散型组合问题的智能优化算法技术领域。特别是大规模柔性作业车间调度优化方法。

背景技术

[0002] 不超过20×50(机床×工件)的车间调度问题为中小规模调度问题,大规模车间调度问题则有以下几种情况【LIANG Xu,WANG Jia,H UANG Ming.New coding method for massive production scheduling problem[J].Computer Integrated Manufacturing System s,2008,10(14):1974-1982】:①工件数J>50,机器数M>20时;②当J≤50,M>20,J×M>1000时;③当J>50,M≤20,J×M>1000时。近60多年来,作业车间的生产排程与调度问题属于NP难题,不少学者对求解此类提出很多求解效果较好的优化方法,但是求解规模通常都是中小规模。大多数优化算法因其求解效率和质量问题,不能直接用于求解大规模调度问题。目前解决大规模调度的方法有:智能优化算法、拉格朗日松弛分解法、基于问题分解法。各种分解法和智能优化算法首先被应用于规模不断增大的车间调度问题。随着工件数、机器数的增加,JSP问题的求解空间变得复杂,传统的求解算法在解的质量和求解时间上不能满足要求。

发明内容

[0003] 本发明的目的是针对现有技术存在的问题而提供一种求解速度快、减小工件完工时间及降低工件拖延期的基于工件组批的大规模柔性作业车间调度方法。
[0004] 本发明的目的是这样实现的:本发明将柔性作业车间生产任务采用聚类的工件组批方法,将大规模的生产任务重组以降低规模,然后利用自适应改进遗传算法求解优化的一种排程调度方法。该方法的具体步骤如下,第一,首先将加工工艺类似、工件尺寸在同一范围内、且毛坯材质相同的工件进行聚类组批,以减小问题求解规模,随机组成不同的批次,每个工件都有各自的交付期,组批之后按照本批当中最早交付期进行调度。
[0005] 第二、构建大规模柔性作业车间调度问题的数学模型
[0006] 所讨论的大规模柔性作业车间调度问题详细描述为,在m台设备(M={Mk|M1,M2...Mm,k=1,2,...m})上加工n件工件(J={Jl|J1,J2,...Jn,l=1,2,...n}),每个工件包含N个事先确定加工顺序的工序,每个工序可以在多台设备上加工。组批之后形成r类工件R={Ri|R1,R2,...Rr,i=1,2,...r},第i类工件的第j道工序在设备k上加工时间为Tijk,Tpi表示第i类工件的第p批次最后一道工序的完工时间,Ds表示工件s的交付期,Dpi表示第i类工件的第p批次的交付期。调度目标是使最大完工时间最小及总拖延期最小。
[0007] 目标函数:minf=α1f'1+α2f'2        (1)
[0008] 其中:f1=min(maxTpi)p=1,2...Bi;i=1,2......r     (2)
[0009]
[0010] 式(2)f1表示最小最大完工时间;式(3)f2表示总拖延期最小。
[0011] 约束条件:
[0012]
[0013]
[0014] 式(5)表示第i类工件第p批次的第j道工序加工时间等于本批次所有工件加工时间之和
[0015] Sip(j-1)k+Tip(j-1)k+Wip(j-1)k≤Sipjk'p=1,2...Bi;i=1,2...r;j=1,2...N;k,k'=1,2...m      (6)
[0016] 式(6)表示同一批次同一个工件的后一道工序的开工时间大于等于前道工序的完工时间
[0017] Sipjk+Tipjk+Wipjk≤Si'p'j'k p,p'=1,2...Bi;i,i'=1,2...r;j,j'=1,2...N;k=1,2...m     (7)
[0018] 式(7)表示占用同一台机器,后一个批次工件的开工时间大于等于前一个批次工件的完工时间
[0019]
[0020] 式(8)表示每类工件所有批次工件数量总和等于所有工件数总和
[0021] Dpi=minDs i=1,2...r;s=1,2......Cip,p=1,2...Bi      (9)
[0022] 式(9)表示第i类工件的第p批次的交付期等于本批中最早交付期
[0023] 具体符号含义见表2。
[0024] 表1符号含义
[0025]
[0026] 第三、采用三层编码进行遗传算法程序中的基因编码
[0027] 本专利为适应大规模柔性作业调度问题而采用三层编码方式。假设有两类工件,每个工件有两道工序,如图2所示,第一层表示工件组批后形成的批次数量,采用整数编码表示第1类工件分为2批。第二层表示工件工序的加工顺序,采用基于工件的编码方式,其中101表示第1类工件的第一个批次,102表示第1类工件的第二个批次,101第一次出现在工序基因内表示其第一道工序,依次类推,工序基因表示的加工顺序为101→201→102→101→
201→102。第三层表示工件工序加工所选择的机器,采用整数编码,表示工件101的第一道工序在可选机器集内选择第一台机器2。
[0028] 第四、遗传算法参数初始化
[0029] 按照初始参数选用的原则和依据,设立种群规模为40,交叉概率0.1,变异概率0.04,最大遗传代数为300代。
[0030] 第五、生成初始种群
[0031] 按照编码方式,采用随机的方式产生初始种群,并且为了保持初始种群的合法性,在产生初始种群时添加约束,即在工序基因段约束工件号出现的次数应等于工件工序数。
[0032] 第六、计算种群适应度值
[0033] 利用公式(1)-(3)计算每个个体的目标函数值,分配适应度值,选择适应度较高进入下一代的交叉变异操作。
[0034] 第七、选择
[0035] 选择算子采用轮盘赌注的方法,因该方法便于操作,故被广泛地运用于选择操作。在此方法中个体被选择的概率和其自身的适应度值相关,适应度值越大被保留下来的概率就越大。
[0036] 第八、选择OBX交叉算子进行遗传算法交叉
[0037] 种群经交叉洗牌后,随机选择两个个体,如图4所示,并截取其工序基因交叉,本专利是采用基于Order-Based Crossover(OBX)的改进交叉算子,其交叉步骤见图4,如图4a所示,在父代P1中随机选择选择待交叉基因,形成交叉基因池Pos(201,101,102),但是由于在Chrom2的批次基因中第一类工件批次为1,因此在Chrom2的工序基因中没有102工件,故从交叉基因池中删除102,最终用于交叉的基因有Pos(201,101)。在父代P2中找到基因201,由于有多道工序所以在P2中会有多个201,则随机选择一个如图4b,然后删除P2中被选中的基因,交叉基因池中其他的基因同理,最后用P1中被选中的基因顺次填补P2中空缺的基因形成子代C2如图4c。而P1中则保留备选中的基因以及不能用于交叉的基因,然后用P2中剩余的基因并去除P1中没有的个体202(第二类工件的第二批工件)依次填补P1中的空缺形成子代C1如图4c。
[0038] 染色体交叉长度会影响算法性能,本专利研究染色体交叉长度与算法求解精度及运算速度之间的关系。选择染色体长度的10%到100%进行交叉,并各自反复进行仿真运算求其平均值,然后求加权值汇总后,发现交叉长度在占染色体长度为10%时最优。为了找出更优的目标值,本专利又选择染色体长度的10%到20%进行交叉,并反复仿真运算求其平均值,然后求加权值,发现选择染色体长度12%进行交叉可以获得更精确的解和提高运算速度,所以选择12%的交叉长度。
[0039] 第九、选择变异策略,以使变异后的基因合法化
[0040] 变异分为批次变异和机器变异,其中批次变异如图5,每次批次变异后工序基因和机器基因都需要重新解码形成新的染色体,以保证染色体的合法性,如第一类工件批次减少1则工序基因相应的减少一个工件102,而第二类工件的批次增加1则工序基因则增加工件202,两者的机器基因均根据工序基因选择相应的机器。机器变异如图6,如Parent(1,9)变异为2,表示选择M11(2,3)中的3号机器。
[0041] 第十、采用自适应遗传算法进行优化求解
[0042] 正弦自适应是交叉率及比变异率根据个体平均适应度值以及最大适应度值按照线性函数曲线变化的一种自适应方法。在优化过程的初始阶段选取较大的交叉率和变异率可以获得较快的收敛速度且保持种群的多样性,经多次迭代后,为避免最优解被破坏,此时选取较小的交叉率和变异率进行细致搜索,该方法所用公式详见式(10)及式(11)【杨建斌,孙树栋,刚刚,王萌.自适应遗传算法求解模糊作业车间调度问题[J].机械科学与技术,2013,(01):16-21.】。经实验验证,正弦自适应遗传算法可以提高算法的收敛速度及收敛精度,并且利用此方法可以保证算法搜索的全局性和精确性。
[0043] 交叉概率:
[0044]
[0045] 变异概率:
[0046]
[0047] 式中fmax为当前种群中最大适应度值;favg为当前种群中平均适应度值;f′为两个待交叉个体中较大的适应度值;f为待变异个体的适应度值;pc1、pm1为系数,在(0,1)区间内取值。
[0048] 第十一、解码
[0049] 解码为编码的逆操作,由染色体的工序基因段得出各类工件每个批次的工序加工顺序,机器基因段得出其相应的加工机器,并且根据式(4)-(9)计算出每道工序的开始及结束时间。
[0050] 本专利提出了一种基于工件组批的大规模作业车间调度方法,该方法首先将加工工艺类似、管径尺寸在同一范围内,且毛坯材质相同的零件进行聚类组批,随机组成不同的批次,每个工件都有各自的交付期,组批之后按照本批当中最早交付期进行调度。然后利用基于OBX交叉自适应遗传算法优化其组批的批次数及批次加工顺序。该方法可以减小问题求解规模,提高求解速度;并且减小工件完工时间以及拖延期。
[0051] 本发明的有益效果是:
[0052] (1)减小问题求解规模,提高求解速度
[0053] 根据组批原则将特性相似的工件进行组批加工,然后利用基于OBX交叉方式的自适应遗传算法进行优化其组批批次数,从而降低问题规模,并大幅度提高求解速度。
[0054] (2)减小工件完工时间及拖延期
[0055] 将同类型的工件进行组批之后可以减少工件的加工准备时间,优化工件组批批次数可以降低工件的拖延期。
[0056] 本发明设置算法初始参数,采用三层基因编码技术、OBX交叉方式和一定的变异策略,结合仿真实验选择交叉长度,并利用自适应改进遗传算法优化求解。本发明提出了一种根据工件的生产特性相似性进行组批并使用改进遗传算法求解大规模柔性作业车间的调度方法,该方法可以减小问题求解规模,提高求解速度;并且减小工件完工时间以及拖延期。这对于求解大规模车间作业调度具有重要的意义和显著的工程实际应用价值。附图说明
[0057] 图1本发明流程图
[0058] 图2染色体编码图。
[0059] 图3父代染色体序列。
[0060] 图4基于OBX的交叉流程图。
[0061] 图5批次变异流程图。
[0062] 图6机器变异流程图。

具体实施方式

[0063] 下面结合附图和具体实施方式对本发明作进一步的描述。
[0064] 实施实例1
[0065] 实例1采用230×18规模的JSP问题,即230个工件在18台机器上进行加工。每一个工件包括四个工序,加工的工序顺序为Qi1、Qi2、Qi3、Qi4。其中Qi1、Qi2和Qi3为加工操作,只能选择在加工的机器上操作;Qi4为检测操作,只能选择在检测的机器上操作。而机器的分布为机器1-15为加工机器,机器16-18为检测机器。
[0066] 实施实例2
[0067] 实例2采用460×18规模的JSP问题,即460个工件在18为加工机器,加工机器和检测机器分布同实例1。
[0068] 具体操作
[0069] 第一、基于工件特性相似组批
[0070] 该方法首先将加工工艺类似、工件尺寸在同一范围内且毛坯材质相同的零件进行聚类组批,随机组成不同的批次,每个工件都有各自的交付期,组批之后按照本批当中最早交付期进行调度。例如有10件导管工件,其主要特性见表1。管径有两个尺寸范围6mm-12mm和14mm-38mm,工件的主要工艺采用A、B、C、D、E、F、G代替,其中E和F可以使用同种机器加工,在工序顺序相同的条件下,工序相同或者不同工序可以采用同类型机器的情况都可以认为是工艺类似,如B→C→F→G和B→C→E→G则可认为两工艺流程相似。按照以上组批方法可得组批结果为:第一组(1、5、9),其交付期为2017/07/1;第二组(2、10),其交付期为2017/07/1;第三组(3、4、6、7、8),其交付期为2017/07/2。
[0071] 表1工件特性表
[0072]
[0073] 第二、构建数学模型
[0074] 按照公式(1)-(9)建立实例1的优化数学模型:
[0075] 目标函数:min f=α1f1'+α2f2'(1)
[0076] 其中:f1=min(maxTpi)i=1,2,.....48;p=1,2......Bi       (2)
[0077]
[0078] 约束方程:
[0079]
[0080]
[0081] Sip(j-1)k+Tip(j-1)k+Wip(j-1)k≤Sipjk'
[0082] (6)
[0083] i=1,2,.....48;p=1,2......Bi;j=1,2,3,4;k,k'=1,2......18
[0084] Sipjk+Tipjk+Wipjk≤Si'p'j'k
[0085] (7)
[0086] i,i'=1,2,.....r;p=1,2......Bi;j,j'=1,2,3,4;k=1,2......18
[0087]
[0088] Dpi=minDs i=1,2,.....r;s=1,2......Cip      (9)
[0089] 第三、采用OBX交叉算子的改进遗传算法进行优化
[0090] 其余实例操作方式一致,获得数学模型后,按照初始参数选用的原则和依据,设立种群规模为40,交叉概率0.1,变异概率0.04,最大遗传代数为300代。
[0091] 采用MATLAB编程,按照上述规则将工件聚类组批,随机产生初始种群,经解码之后,按照一定概率选择出优秀的个体进入下一代交叉、变异更新种群,循环迭代直到达到最大迭代代数后,输出最优完工时间及对应的调度解。
[0092] 两种不同规模的实例分别运用非组批自适应遗传算法和本专利的组批自适应遗传算法各运算十次,以进行相关数据比较,说明该发明中优化方法的有效性。计算后各规模的JSP问题的最优解及相关数据见表4。根据不同的实验数据结果,对性能进行比较,从以下数据可见最小完成时间/拖延期和运算效率等性能指标都得到了较大的提高。
[0093] 最小完成时间与搜索性能
[0094] 表4实验结果数据表
[0095]
[0096] 由表4可知,两种不同大规模的JSP问题经过组批调度之后拖延期均为0,而非组批调度方法均有不同程度的拖延;最大完工时间至少缩减21%,最多缩减25%;运算效率至少提高61%,最大提高66%。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈