CN111026051B - 一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法 - Google Patents
一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法 Download PDFInfo
- Publication number
- CN111026051B CN111026051B CN201911200006.2A CN201911200006A CN111026051B CN 111026051 B CN111026051 B CN 111026051B CN 201911200006 A CN201911200006 A CN 201911200006A CN 111026051 B CN111026051 B CN 111026051B
- Authority
- CN
- China
- Prior art keywords
- casting
- equipment
- population
- time
- heat treatment
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 227
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 55
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 39
- 238000005266 casting Methods 0.000 claims abstract description 230
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 42
- 238000005265 energy consumption Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims description 133
- 238000012545 processing Methods 0.000 claims description 129
- 238000010438 heat treatment Methods 0.000 claims description 58
- 239000000428 dust Substances 0.000 claims description 57
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 30
- 238000003754 machining Methods 0.000 claims description 26
- 239000003345 natural gas Substances 0.000 claims description 15
- 238000007528 sand casting Methods 0.000 claims description 14
- 230000003111 delayed effect Effects 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 8
- 238000003723 Smelting Methods 0.000 claims description 7
- 238000009434 installation Methods 0.000 claims description 6
- 238000000465 moulding Methods 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 229910052757 nitrogen Inorganic materials 0.000 claims description 4
- 238000004321 preservation Methods 0.000 claims description 4
- 239000003245 coal Substances 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 229910052698 phosphorus Inorganic materials 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000002922 simulated annealing Methods 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 6
- 239000000956 alloy Substances 0.000 description 3
- 229910045601 alloy Inorganic materials 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000010924 continuous production Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000269350 Anura Species 0.000 description 1
- 238000010923 batch production Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000003110 molding sand Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000009827 uniform distribution 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
- 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)
- General Factory Administration (AREA)
Abstract
本发明涉及一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法。通过对典型铸造企业的设备能耗状况和铸造工艺流程进行***分析,找出其存在的主要问题。构建碳排放约束下的柔性铸造流水车间调度模型,设计改进蛙跳算法,并通过此算法求解低碳调度模型。本发明自上而下对算法的个体编解码、种群初始化、子种群个体搜索策略、子种群个体更新方式、种群最优个体更新方式等进行全面改进以实现对考虑设备碳排放约束下的铸造车间调度问题的求解,同时引入模拟退火局部搜索机制以提高改进算法的性能。
Description
技术领域
本发明涉及柔性铸造流水车间的低碳调度方法,具体涉及基于改进蛙跳算法的低碳调度方法。
背景技术
为推动促进我国铸造行业朝着绿色化、数字化、智能化工厂生产模式转型升级,实现铸件生产中铸造设备能耗状况量化分析和铸件单件的可调度性就显得十分重要。通过对典型铸造企业的设备能耗状况和铸造工艺流程进行***分析,发现其存在的主要问题如下:
1)铸造设备能耗量化分析困难
铸造设备不同于精加工机床,铸造生产中所需加工设备、辅助设备和运输设备种类复杂,铸造设备处于待机状况、负载状态及空载状态下的碳排放量化方法均不相同。绝大部分工序具有特定加工设备和辅助设备,其碳排放量化过程较为困难。
2)铸造工艺流程复杂且繁琐
绝大多数铸件生产属于单件小批量生产模式且生产流程较长。铸造工艺存在多重约束。真实的铸造车间并非零等待作业车间,同一铸件的造型工序与制芯工序不存在前后约束限制条件且可独立并行加工。铸件相邻工序之间存在着固定的检测时间间隔和运输时间。
铸造车间排产调度方式与传统精加工车间排产调度方式差别较大且目前铸造行业普遍缺乏成熟的生产调度方案。目前绝大多数与铸造相关的调度问题主要集中在某几道工序批调度研究上,而通过优化调度方案进而降低铸造设备碳排放的研究还处于一片空白。
本所发明涉及铸造车间调度问题针对砂型铸造生产的离散工艺工部,原因如下:
1)某批次所需生产的铸件种类和数量确定时,砂型铸造生产的连续工艺工部即混砂工部所需配制的型砂和芯砂的吨量确定,熔炼工部所需熔炼的合金种类和吨量确定。当连续工艺工部的设备利用率和生产技术水平保持稳定的前提下,混砂工部和熔炼工部的加工时间、设备负荷和设备碳排放等指标确定且无法通过调度方案选择进行优化。
2)在实际砂型铸造生产中,铸件的每道离散工序都对应于不少于一台的可选加工(除尘)设备集,加工(除尘)设备集中不同设备的待机功率、空载功率和单位负载功率均不完全相同。当铸件任意离散工序选择不同铸造加工(除尘)设备时,铸件的安装时间、负载时间和卸载时间均不完全相同,铸件在任意相邻离散工序之间的运输时间由所选铸造加工(除尘)设备之间的距离决定。因此通过合理为铸件离散工序合理安排加工顺序和选择加工(除尘)设备不仅可以达到缩短最大完工时间、设备最大负荷和设备总负荷,还可以达到降低加工(除尘)设备和运输设备的碳排放量的目的。
混合蛙跳算法作为一种基于群体智能的全局寻优算法,具有概念简单、参数少、计算速度快和全局搜索寻优能力强的特征,但,SFL算法所求解的目标函数值容易陷入局部最优而不是全局最优解,且算法中青蛙的子族群周期性进行聚集并进行食物源信息交换和共享并不全面。
发明内容
针对跳蛙算法自身存在易陷入局部最优及信息共享不足的缺陷,本发明提出了一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,实现对考虑设备碳排放约束下的铸造车间调度问题的求解,同时引入模拟退火局部搜索机制以提高改进算法的性能。
本发明所设计的基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于,包括以下步骤:
1、设备碳排放约束下的柔性铸造流水车间调度模型构建
1.1对研究的柔性铸造流水车间调度模型进行必要的假设;
1.2设置参数;
1.3设置约束条件;
1.4构建模型的优化目标函数;
2、改进蛙跳算法设计与实现
2.1选择合适的方法对个体进行编码和解码;
2.2种群初始化;
2.3计算Pareto非支配排序与个体拥挤距离;
2.4针对提出的模型设计了更优引导式和突变跳跃式两种局部搜索策略;在每次全局迭代进化过程中对每个子种群进行Pareto排序,随机选择某个rank值为1的个体作为子种群最优解Xb,然后鉴于SA良好的局部搜索能力,基于SA对各子种群的Xb进行局部搜索进而提高算法的寻优能力避免算法陷入局部次优解。
2.5采用锦标赛法构建子种群;
2.6改进子种群最劣解Xw的更新方式;
2.7基于最小算术位置差的全局最优解Xg更新方式;
2.8判断种群的迭代进化次数是否达到设定次数,若达到则输出种群中所有Pareto非支配等级为1的个体并结束算法;否则循环执行步骤2.4至2.7。
本发明所构建的考虑碳排放约束的柔性铸造流水车间调度模型需要作出以下假设:
(1)在零时刻,任意设备均可被使用,任意铸件均可被加工;
(2)任何铸件的任意工序在任意时刻只可选择一台设备进行加工且每台设备在任意时刻只能加工一个铸件的一道工序,任意设备前可容纳无限量的待加工铸件;
(3)不同铸件间没有加工顺序约束,同一铸件的造型工序与制芯工序间没有加工顺序约束,其它工序之间存在加工顺序约束;
(4)每个铸件的工序加工顺序不可改动,且除最后一道工序外的其它每道工序被加工完毕后必须经过一定时间的休整方可运输至下一所选加工设备进行加工;
(5)同一铸件任意工序的加工过程均不可间断且同一铸件在相邻工序间的运输过程不可中断;
(6)铸件的任意一道工序一旦完成加工经过必要休整即被运输到下一相邻工序所选加工设备进行加工,不存在铸件成批运输的情况;
(7)对于铸件存在除尘工艺的工序的任一可选加工设备对应唯一可选除尘设备,规定加工设备与其对应除尘设备在同一时刻启动和关闭且任一加工设备与其对应除尘设备待机时间、空载时间和负载时间均相同;
(8)规定加工设备、除尘设备和运输设备的启动时间和启动状态下的能耗和碳排放忽略不计;
(9)规定任意加工(除尘)设备在首个铸件到达的安装时刻启动,在最后一个铸件完工的卸载时刻关闭。铸件在加工(除尘)设备上安装、卸载过程中,加工(除尘)设备处于空载;铸件在加工(除尘)设备上加工(除尘)和运输过程中,加工(除尘)设备处于负载,其他时间加工(除尘)设备处于待机;
(10)运输设备即行车类型相同且数量充足,行车在安装铸件时启动,卸载铸件时停止,行车运输过程保持恒速,忽略行车待机、空载状态的能耗且铸件在相邻工序间运输时不存在运输等待时间;
(11)规定单批次生产的铸件为同种材质且所需熔炼合金总吨量不超过熔炼设备单批次可熔炼合金的最大容纳吨量;
具体参数设置如下:
C:铸件集合;
E:铸造加工设备集合(含热处理设备);
H:铸造热处理设备集合;
D:铸造除尘设备集合;
TS:运输设备集合;
n:待生产铸件总量;
m:铸造加工设备总量(含热处理设备);
r:铸造热处理设备总量;
m’:铸造除尘设备总量;
y:运输设备总数;
i:铸件序号,i=1,2,3…n;
Ci:铸件集合中序号为i的铸件;
Mi:铸件Ci的理论重量;
ni:铸件Ci的工序总数;
h:铸造加工设备的序号,h=1,2,3…m(包含热处理设备);
Eh:铸造加工设备集合中序号为h的设备(包含热处理设备);
h’:铸造热处理设备序号,h’=1,2,3…r;
Hh’:铸造热处理设备集合中序号为h’的设备;
c:铸造除尘设备序号c=1,2,3…m’;
Dc:铸造除尘设备集合中序号为c的设备;
s:铸造运输设备序号s=1,2,3…y;
TSs:铸造运输设备中序号为s的设备;
Oij:铸件Ci的第j道工序;
Oij’:铸件Ci的第j道工序,该工序需要进行除尘操作;
Oi1:铸件Ci的并行工序中的造型工序;
Oi2:铸件Ci的并行工序中的制芯工序;
Oi3:铸件Ci的下芯合箱工序;
Oip:铸件Ci的浇注工序;
Oij”:铸件Ci的热处理工序;
Eij:工序Oij的可选加工设备集合;
Oijh:工序Oij在Eij中所选的加工设备h上加工;
Oij’hc:工序Oij’在Eij中所选的加工设备h所对应的除尘设备c上进行除尘操作;
Oij”h’:工序Oij”在热处理设备集合H中选用设备h’进行加工;
pijh:工序Oij在其所选加工设备h上的安装时间;
tijh:工序Oij在其所选加工设备h上的负载时间;
uijh:工序Oij在其所选加工设备h上的卸载时间;
tij’c:工序Oij’在其所选除尘设备c上的负载时间;
tij”h’:工序Oij”在其所选热处理设备h’上的负载时间;
Jij:工序Oij与其下一道相邻工序之间的固定休整时间;
tijlh:工序Oij所选加工设备l到其下道相邻工序的所选加工设备h的运输时间;
Delijh:工序Oij在其所选用加工设备h上由于加工设备资源约束而延迟的时间;
STijh:工序Oij在其所选加工设备h上的最早可开工时间;
Cijh:工序Oij在其所选加工设备h上最早完工时间;
Ci’j’h:加工设备h上加工的工序Oij的上一道相邻工序Oi’j’在加工设备h上的最早完工时间;
Psl:浇注工序Oij’的最早可加工限制时刻;
Delip:浇注工序Oip延迟加工等待时间;
Deliph:浇注工序Oip由于加工设备资源约束而延迟的时间;
Delipx:浇注工序Oip由于最早加工时刻Psl限制而延迟的时间;
Mt:砂型铸造熔炼工序的特定加工时间区间;
Iwt:砂型铸造保温电炉的理论负载时间;
Ikt:砂型铸造保温电炉的理论空载时间;
Oi’j’:加工设备h上所加工的工序Oij的上一道相邻工序;
CTi:铸件Ci的完工时间;
Psh:加工设备Eh的待机功率;
Pkh:加工设备Eh的空载功率;
Plh:加工设备Eh的负载功率;
Psc:除尘设备Dc的待机功率;
Pkc:除尘设备Dc的空载功率;
Plc:除尘设备Dc的负载功率;
Psh’:热处理设备Hh’的待机功率;
Pkh’:热处理设备Hh’的空载功率;
Plh’:热处理设备Hh’的负载功率;
Th'l:热处理设备Hh’单位负载重量与单位时间下的天然气排放量(单位:m3/(kg·h))
Pki:铸造保温电炉的空载功率;
Pli:铸造保温电炉的负载功率;
Pks:运输设备行车的空载功率;
Pls:运输设备行车的负载功率;
Ee:加工设备(包含热处理设备)总电能消耗量,单位kw.h;
Ed:除尘设备电能总消耗量,单位kw.h;
Ets:运输设备电能总消耗量,单位kw.h;
Ei:铸造保温电炉电能总消耗量,单位kw.h;
Eht:热处理设备天然气总消耗量,单位m3;
SCe:电力标准折煤系数,0.1229kgce/(kW·h);
SCt:天然气标准折煤系数,1.2143kgce/m3;
EFe:电能的碳排放系数,4.035kgCO2e/kgce;
EFt:天然气的碳排放系数,1.744kgCO2e/kgce;
CBz:铸造设备总碳排放量,单位kgCO2e;
Tt:铸造设备总负荷,单位h;
决策变量:
Xijh:0-1变量,工序Oij在加工设备h上加工时则Xijh=1,否则为0;
Xij”h’:0-1变量,工序Oij”在热处理设备h’上加工则Xij”h’=1,否则为0;
Xijd:0-1变量,工序Oij存在除尘操作则Xijd=1,否则为0;
Xijc:0-1变量,工序Oij在除尘设备c上除尘则Xijc=1,否则为0;
Xijs:0-1变量,工序Oij与其下一道相邻工序之间在运输设备c上运输则Xijs=1,否则为0;
约束条件:
1)铸件的任意工序Oij一旦开始加工就不可中断其加工过程,工序Oij的完工时间Cijh大于等于该工序在其所选加工设备h的最早可加工时间STijh与安装时间pijh,加工时间tijh和卸载时间uijh之和,具体如下式:
Cijh≥STijh+pijh+tijh+uijh
2)考虑到铸件的造型工序Oi1和制芯工序Oi2并行加工的特性,铸件的下芯合箱工序Oi3的STi3h存在两种取值结果,如公式所示:
3)同一铸件任意相邻工序之间有着严格的加工顺序约束,下式表示同一铸件的某道工序Oi(j+1)必须在前一道工序Oij结束后经过一定的固定休整时间Jij和运输时间tijhk后才可进行加工:
STi(j+1)k≥Cijh+Jij+tijhk
4)工序Oi’j’表示工序Oij之前在加工设备h上所加工的上一道相邻工序,下式表明各铸造加工设备在同一时刻只能加工一道工序。
STijh≥Ci'j'h
5)铸件的任意一道工序Oij存在多个可选加工设备,但其一次只能被可选加工设备集合Eij中的某个加工设备h单独加工。
6)当工序Oij不属于浇注工序Oip,此时Delijh具体表示安排在设备h上的工序Oij因设备资源约束而延迟的时间。如图5所示,Delij理论上存在两种不同取值结果:当铸件Ci被运输到加工设备Eh的时刻,加工设备Eh处于空闲状态,工序Oij直接开始加工。当铸件Ci被运输到加工设备Eh的时刻,工序Oi’j’正在加工设备Eh上加工,工序Oij需等待工序Oi’j’加工完成才可加工。Delijh的两种取值结果如下式所示:
7)对于铸造生产的浇注工序Oip,Delip具体表示工序Oip因设备资源约束而延迟的时间Deliph和最早可加工限制时刻Psl限制而延迟的时间Delipx之和,如公式所示:
Delip=Delipx+Deliph
如图6所示,Delip理论上存在五种不同取值结果:当STiph≤Psl时,则Delij'x=Psl-STiph,此时Deliph存在两种不同的取值结果:
当STiph>Psl且Ci'j'h≤Psl时,则此时Delijh=0且Delipx=0;当STiph>Psl且Ci'j'h>Psl时,则Delipx=0,此时Delijh存在两种不同的取值结果:
8)砂型铸造保温电炉的理论负载时间Iwt和理论空载时间Ikt受到批次铸件熔炼工序的结束时刻Psl,最早可开始浇注的时刻minSTiph,最晚浇注完工时刻maxCTiph的限制约束。对于Iwt则有:
对于Ikt则有:
目标函数的构建:
1)最小化铸件最大完工时间
min f1=max(Ci)
其中
2)最小化铸造设备碳排放量
该指标是指参与该批次铸件生产的普通铸造工艺加工设备、除尘设备、热处理设备、铸造保温电炉及运输设备的碳排放量总和,本文对铸造设备分别处于待机、空载和负载状态下的能源消耗状况进行分析:
min f2=CBz
其中
CBz=(Ee+Ed+Ets+Ei)·SCe·EFe+Eht·SCt·EFt
Ee为加工设备电力消耗量:
Ed为除尘设备电力消耗量:
Ets为运输设备电力消耗量:
Ei为铸造保温电炉电力消耗量:
Ei=Pki·Ikt+Pwi·Iwt
Eht为热处理设备负载状态下的天然气消耗量:
3)最小化铸造设备总负荷
该指标是指所有参与铸造调度的加工设备、除尘设备、运输设备、保温电炉负荷的总和:
min f3=Tt
其中
2、算法设计
1)个体的编码与解码
ISFLA算法的种群中每只青蛙个体均代表铸造车间调度问题解空间的一个解。编码是对青蛙个体进行进化操作,生成解和评估解的重要载体。柔性铸造流水车间调度问题的编码主要包括两个部分,即加工(除尘)设备选择和离散工序排序。针对这一特点,本发明采用MSOS整数编码方式。该编码方式由两段组成:加工设备选择段(Xm)和工序排序段(X0)。两段长度均为D,它们一起构成了种群青蛙个体的位置信息。
2)种群初始化
本发明结合所研究的铸造车间调度问题考虑了铸件运输时间的特征,采用基于运输时间的最短工作时间法(Short Working Machine,SWM)策略和随机选择初始化(RadomSelection,RS)加工设备分配部分,工序排序部分由随机方法生成。通过实验测试将SWM和RS两种初始化策略的比例进行设定为1:1,使其在保证初始解均匀分布和避免算法早熟的前提下进而提高了算法搜索效率。
3)Pareto非支配排序与个体拥挤距离计算
改进的Pareto排序方法由以下两部分组成,第一部分为所有解设置了两个变量np和Sp,其中P=1,2,…N,N为种群大小,np用来记录所有解中支配解P的个数,Sp用来记录所有解中被P支配的解集。第二部分根据第一部分排序后的结果进行分层,并初始化分层序号为1,首先将种群中所有np为0的个体移除种群,将当前非支配排序分层的序号赋予这些个体;然后将这些个体Sp中所对应的个体np减1,并将分层序号递增加1,如此反复循环执行直到青蛙种群中所有个体的分层序号被赋予。
拥挤距离是对同一排序分层内个体在每一个目标上和它相邻两个体距离之差进行绝对值求和。个体i在第k个目标fk上的拥挤距离为k=1,2,…m,m为目标的个数 是个体i在第k个目标上相邻两个体i-1和i+1的目标值。个体i的拥挤距离di如下式所示:
4)基于模拟退火(SA)的子种群最优解Xb局部搜索策略
本发明所采用的ISFLA算法在每次全局迭代进化过程中对每个子种群进行Pareto排序,随机选择某个rank值为1的个体作为子种群最优解Xb,然后鉴于SA良好的局部搜索能力,基于SA对各子种群的Xb进行局部搜索进而提高算法的寻优能力避免算法陷入局部次优解。
5)基于锦标赛法的子种群构建方式
ISFLA算法采用锦标赛法构建子种群,即随机选择初始种群中的两个个体X1和X2,如果X1支配X2则让X1进入第一个子种群,同时将X2放回初始种群中;若两个个体彼此非受支配关系则随机选择一个个体进入第一个种群并将未选中个体放回初始种群;随机选择两个个体重复以上过程依次构建所以子种群直到所有子种群均被构建完成。
6)子种群最劣解Xw更新方式
ISFLA算法在子种群每次迭代进化过程中进行Pareto排序,随机选择某个rank值最大的个体作为子种群最劣解Xw。然后基于子种群最优解Xb对子种群最劣解Xw的工序排序段和加工设备排序段分别进行POX交叉和RPX交叉,若交叉后产生的新解Xnew支配Xw,则将Xw用Xnew替换,否则基于全局最优解Xg对子种群最劣解Xw进行POX交叉和RPX交叉,若交叉产生的新解Xnew不被Xw支配,则将Xw用Xnew替换,否则对Xw执行翻转变异,将Xw用变异产生的新解Xnew替代。
7)基于最小算术位置差的全局最优解Xg更新方式
本发明基于青蛙个体每次跳跃后的位置相对于全局最优解Xg的位置的差值提出位置差的概念,具体实现方法是首先对各子种群青蛙个体混合并进行Pareto非支配排序,然后计算所有非支配等级为1的青蛙个体Xrank=1的位置的差,由于所研究铸造调度模型的三个目标函数都是最小化问题,因此选择一个最小算术位置差的个体作为新的全局最优解Xg+1。
本发明的优点在于:
本发明构建了低碳铸造车间生产调度问题模型,并针对该模型设计改进求解算法,为铸造企业提高生产效益和节能减排问题提供一套可选择的解决方案。
跳蛙算法用于求解设备碳排放约束下的柔性铸造流水车间调度问题的有效性和优越性。
1)本发明根据铸件实际生产过程中的工艺约束特征与铸造设备能耗特征,提出一种考虑设备碳排放约束的柔性铸造流水车间调度策略,进而实现了绿色铸造调度理论在实际铸件生产中的应用,弥补了低碳铸造车间调度相关学术研究领域上的空白。
2)本发明构建了低碳铸造车间生产调度问题模型,并针对该模型设计改进求解算法,以实际铸造企业生产数据作为实例,开发一套铸造设备碳排放优化调度管理***,以实现铸件生产过程和设备能耗精细化管控为目标,为铸造企业提高生产效益和节能减排问题提供一套可选择的解决方案。
3)改进的蛙跳算法求解设备碳排放约束下的柔性铸造流水车间调度模型时有效降低了计算的复杂度。而且本算法所设计的更优引导式和突变跳跃式两种局部搜索策略可以使算法有效的跳出局部最优解。
附图说明
图1为本发明中的改进跳蛙算法具体的实施流程图。
图2为4种算法Pareto第一前沿3维分布图。
图3为4种算法Pareto第一前沿2维分布图。
图4为设备碳排量目标值最小的解方案对应的甘特图。
图5为当工序Oij不属于浇注工序Oip,Delij理论上存在两种不同取值结果。
图6为对于铸造生产的浇注工序Oip,Delip理论上存在五种不同取值结果。
具体实施方式
本发明所设计的基于改进蛙跳(ISFLA)算法的柔性铸造流水车间低碳调度方法,包括以下步骤:
1、设备碳排放约束下的柔性铸造流水车间调度模型构建
1.1对研究的柔性铸造流水车间调度模型进行必要的假设;
1.2设置参数;
1.3设置约束条件;
1.4构建模型的优化目标函数;
2、改进蛙跳算法(ISFLA)设计与实现
2.1选择合适的方法对个体进行编码和解码;
2.2种群初始化;
2.3计算Pareto非支配排序与个体拥挤距离;
2.4针对提出的模型设计了更优引导式和突变跳跃式两种局部搜索策略;
2.5采用锦标赛法构建子种群;
2.6改进子种群最劣解Xw的更新方式;
2.7基于最小算术位置差的全局最优解Xg更新方式;
2.8判断种群的迭代进化次数是否达到设定次数,若达到则输出种群中所有Pareto非支配等级为1的个体并结束算法;否则循环执行步骤2.4至2.7。
其中,设置参数具体包括:
C:铸件集合;
E:铸造加工设备集合(含热处理设备);
H:铸造热处理设备集合;
D:铸造除尘设备集合;
TS:运输设备集合;
n:待生产铸件总量;
m:铸造加工设备总量(含热处理设备);
r:铸造热处理设备总量;
m’:铸造除尘设备总量;
y:运输设备总数;
i:铸件序号,i=1,2,3…n;
Ci:铸件集合中序号为i的铸件;
Mi:铸件Ci的理论重量;
ni:铸件Ci的工序总数;
h:铸造加工设备的序号,h=1,2,3…m(包含热处理设备);
Eh:铸造加工设备集合中序号为h的设备(包含热处理设备);
h’:铸造热处理设备序号,h’=1,2,3…r;
Hh’:铸造热处理设备集合中序号为h’的设备;
c:铸造除尘设备序号c=1,2,3…m’;
Dc:铸造除尘设备集合中序号为c的设备;
s:铸造运输设备序号s=1,2,3…y;
TSs:铸造运输设备中序号为s的设备;
Oij:铸件Ci的第j道工序;
Oij’:铸件Ci的第j道工序,该工序需要进行除尘操作;
Oi1:铸件Ci的并行工序中的造型工序;
Oi2:铸件Ci的并行工序中的制芯工序;
Oi3:铸件Ci的下芯合箱工序;
Oip:铸件Ci的浇注工序;
Oij”:铸件Ci的热处理工序;
Eij:工序Oij的可选加工设备集合;
Oijh:工序Oij在Eij中所选的加工设备h上加工;
Oij’hc:工序Oij’在Eij中所选的加工设备h所对应的除尘设备c上进行除尘操作;
Oij”h’:工序Oij”在热处理设备集合H中选用设备h’进行加工;
pijh:工序Oij在其所选加工设备h上的安装时间;
tijh:工序Oij在其所选加工设备h上的负载时间;
uijh:工序Oij在其所选加工设备h上的卸载时间;
tij’c:工序Oij’在其所选除尘设备c上的负载时间;
tij”h’:工序Oij”在其所选热处理设备h’上的负载时间;
Jij:工序Oij与其下一道相邻工序之间的固定休整时间;
tijlh:工序Oij所选加工设备l到其下道相邻工序的所选加工设备h的运输时间;
Delijh:工序Oij在其所选用加工设备h上由于加工设备资源约束而延迟的时间;
STijh:工序Oij在其所选加工设备h上的最早可开工时间;
Cijh:工序Oij在其所选加工设备h上最早完工时间;
Ci’j’h:加工设备h上加工的工序Oij的上一道相邻工序Oi’j’在加工设备h上的最早完工时间;
Psl:浇注工序Oij’的最早可加工限制时刻;
Delip:浇注工序Oip延迟加工等待时间;
Deliph:浇注工序Oip由于加工设备资源约束而延迟的时间;
Delipx:浇注工序Oip由于最早加工时刻Psl限制而延迟的时间;
Mt:砂型铸造熔炼工序的特定加工时间区间;
Iwt:砂型铸造保温电炉的理论负载时间;
Ikt:砂型铸造保温电炉的理论空载时间;
Oi’j’:加工设备h上所加工的工序Oij的上一道相邻工序;
CTi:铸件Ci的完工时间;
Psh:加工设备Eh的待机功率;
Pkh:加工设备Eh的空载功率;
Plh:加工设备Eh的负载功率;
Psc:除尘设备Dc的待机功率;
Pkc:除尘设备Dc的空载功率;
Plc:除尘设备Dc的负载功率;
Psh’:热处理设备Hh’的待机功率;
Pkh’:热处理设备Hh’的空载功率;
Plh’:热处理设备Hh’的负载功率;
Th'l:热处理设备Hh’单位负载重量与单位时间下的天然气排放量(单位:m3/(kg·h))
Pki:铸造保温电炉的空载功率;
Pli:铸造保温电炉的负载功率;
Pks:运输设备行车的空载功率;
Pls:运输设备行车的负载功率;
Ee:加工设备(包含热处理设备)总电能消耗量,单位kw.h;
Ed:除尘设备电能总消耗量,单位kw.h;
Ets:运输设备电能总消耗量,单位kw.h;
Ei:铸造保温电炉电能总消耗量,单位kw.h;
Eht:热处理设备天然气总消耗量,单位m3;
SCe:电力标准折煤系数,0.1229kgce/(kW·h);
SCt:天然气标准折煤系数,1.2143kgce/m3;
EFe:电能的碳排放系数,4.035kgCO2e/kgce;
EFt:天然气的碳排放系数,1.744kgCO2e/kgce;
CBz:铸造设备总碳排放量,单位kgCO2e;
Tt:铸造设备总负荷,单位h;
决策变量:
Xijh:0-1变量,工序Oij在加工设备h上加工时则Xijh=1,否则为0;
Xij”h’:0-1变量,工序Oij”在热处理设备h’上加工则Xij”h’=1,否则为0;
Xijd:0-1变量,工序Oij存在除尘操作则Xijd=1,否则为0;
Xijc:0-1变量,工序Oij在除尘设备c上除尘则Xijc=1,否则为0;
Xijs:0-1变量,工序Oij与其下一道相邻工序之间在运输设备c上运输则Xijs=1,否则为0。
设置的约束条件具体如下:
Cijh≥STijh+pijh+tijh+uijh;
STi(j+1)k≥Cijh+Jij+tijhk;
STijh≥Ci'j'h;
Delip=Delipx+Deliph;
所构建的目标函数包括:
1.41最小化铸件最大完工时间
min f1=max(Ci)
其中
1.42最小化铸造设备碳排放量
min f2=CBz
其中
CBz=(Ee+Ed+Ets+Ei)·SCe·EFe+Eht·SCt·EFt
Ee为加工设备电力消耗量:
Ed为除尘设备电力消耗量:
Ets为运输设备电力消耗量:
Ei为铸造保温电炉电力消耗量:
Ei=Pki·Ikt+Pwi·Iwt
Eht为热处理设备负载状态下的天然气消耗量:
1.43最小化铸造设备总负荷
minf3=Tt
其中
改进蛙跳算法的具体实现基本步骤如下:
以杭州某铸造企业A典型生产周期内的铸件生产作为实例,然后通过对铸件和相关设备实例数据进行实验测试进而设定ISFLA算法与其他几种常见对比算法的基本参数,最后通过调度实例结果分析验证了所采用的ISFLA算法在求解设备碳排放约束下的铸造调度问题方面的有效性和优越性。
企业A以生产砂型铸件为主,其订单主要具有小批量和单材质的特点,其某一典型生产周期需要生产的铸件基本参数、铸造设备相关参数、铸件加工时间、铸件工序固定休整时间和铸件运输时间等信息要素分别见表1、2、3、4、5。
表1铸件基本参数信息
表2铸造设备相关参数信息
表3铸件加工时间参数信息
表4铸件工序固定休整时间(h)
表5铸件运输时间信息要素(min)
EQX表示设备EX作为运输起点加工设备,EZX表示设备EX作为运输终点加工设备。
为评估本发明所采用ISFLA算法的性能,将其与原始的SFLA算法、多目标粒子群算法(HDPSO)以及原始的NSGA-II算法进行对比。四种算法的具体参数设置如表6所示。
表6 ISFLA、SFLA、NSGA-II和HDPSO的基本参数设置
四种算法均在Windows 10***、CPU主频四核2.3GHz、运行内存8GB的计算机中运行,实验实施条件为MATLAB R2016a on an Intel Core i7 2.30GHz PC with 8.00GB ofmemory.将每种算法每次运行后得到的Pareto第一前沿解记录下来,四种算法的Pareto第一前沿解的三维分布如图2所示,其第一前沿的二维分布如图3所示。其中,(a)4种算法的f1、f2分布图,(b)4种算法的f1、f3分布图(c)4种算法的f2、f3分布图。由图可知所提算法获得的前沿排列在其他3种算法第一前沿的前列,说明所采用的ISFLA算法较其它三种算法寻优能力更强,性能更好。
为降低铸造设备能耗、减少铸造设备碳排放和响应国家提出的绿色制造理念,选择本实验中ISFLA算法10次运行所得Pareto第一前沿解集中铸造设备碳排量目标值最小的解(39.63333,7811.227,80.38333)作为决策人员所选择最优解,其调度甘特图如图4所示。图中每一个矩形代表铸件的一道工序,矩形的长度表示该道工序的加工所需时间矩形中符号表示加工工序信息,例如符号为106的矩形表示编号为1的铸件的第6道加工工序在13号铸造加工设备上加工。
Claims (9)
1.一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于,包括以下步骤:
1、设备碳排放约束下的柔性铸造流水车间调度模型构建
1.1对研究的柔性铸造流水车间调度模型进行必要的假设;
1.2设置参数;
1.3设置约束条件;
1.4构建模型的优化目标函数;
2、改进蛙跳算法设计与实现
2.1选择合适的方法对个体进行编码和解码;
2.2种群初始化;
2.3计算Pareto非支配排序与个体拥挤距离;
2.4针对提出的模型设计了更优引导式和突变跳跃式两种局部搜索策略;在每次全局迭代进化过程中对每个子种群进行Pareto排序,随机选择某个rank值为1的个体作为子种群最优解Xb,然后鉴于SA良好的局部搜索能力,基于SA对各子种群的Xb进行局部搜索进而提高算法的寻优能力避免算法陷入局部次优解;
2.5采用锦标赛法构建子种群;
2.6改进子种群最劣解Xw的更新方式,
2.7基于最小算术位置差的全局最优解Xg更新方式;
2.8判断种群的迭代进化次数是否达到设定次数,若达到则输出种群中所有Pareto非支配等级为1的个体并结束算法;否则循环执行步骤2.4至2.7。
2.根据权利要求1所述的基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于:所述设置参数具体包括:
C:铸件集合;
E:铸造加工设备集合(含热处理设备);
H:铸造热处理设备集合;
D:铸造除尘设备集合;
TS:运输设备集合;
n:待生产铸件总量;
m:铸造加工设备总量(含热处理设备);
r:铸造热处理设备总量;
m’:铸造除尘设备总量;
y:运输设备总数;
i:铸件序号,i=1,2,3…n;
Ci:铸件集合中序号为i的铸件;
Mi:铸件Ci的理论重量;
ni:铸件Ci的工序总数;
h:铸造加工设备的序号,h=1,2,3…m(包含热处理设备);
Eh:铸造加工设备集合中序号为h的设备(包含热处理设备);
h’:铸造热处理设备序号,h’=1,2,3…r;
Hh’:铸造热处理设备集合中序号为h’的设备;
c:铸造除尘设备序号c=1,2,3…m’;
Dc:铸造除尘设备集合中序号为c的设备;
s:铸造运输设备序号s=1,2,3…y;
TSs:铸造运输设备中序号为s的设备;
Oij:铸件Ci的第j道工序;
Oij’:铸件Ci的第j道工序,该工序需要进行除尘操作;
Oi1:铸件Ci的并行工序中的造型工序;
Oi2:铸件Ci的并行工序中的制芯工序;
Oi3:铸件Ci的下芯合箱工序;
Oip:铸件Ci的浇注工序;
Oij”:铸件Ci的热处理工序;
Eij:工序Oij的可选加工设备集合;
Oijh:工序Oij在Eij中所选的加工设备h上加工;
Oij’hc:工序Oij’在Eij中所选的加工设备h所对应的除尘设备c上进行除尘操作;
Oij”h’:工序Oij”在热处理设备集合H中选用设备h’进行加工;
pijh:工序Oij在其所选加工设备h上的安装时间;
tijh:工序Oij在其所选加工设备h上的负载时间;
uijh:工序Oij在其所选加工设备h上的卸载时间;
tij’c:工序Oij’在其所选除尘设备c上的负载时间;
tij”h’:工序Oij”在其所选热处理设备h’上的负载时间;
Jij:工序Oij与其下一道相邻工序之间的固定休整时间;
tijlh:工序Oij所选加工设备l到其下道相邻工序的所选加工设备h的运输时间;
Delijh:工序Oij在其所选用加工设备h上由于加工设备资源约束而延迟的时间;
STijh:工序Oij在其所选加工设备h上的最早可开工时间;
Cijh:工序Oij在其所选加工设备h上最早完工时间;
Ci’j’h:加工设备h上加工的工序Oij的上一道相邻工序Oi’j’在加工设备h上的最早完工时间;
Psl:浇注工序Oij’的最早可加工限制时刻;
Delip:浇注工序Oip延迟加工等待时间;
Deliph:浇注工序Oip由于加工设备资源约束而延迟的时间;
Delipx:浇注工序Oip由于最早加工时刻Psl限制而延迟的时间;
Mt:砂型铸造熔炼工序的特定加工时间区间;
Iwt:砂型铸造保温电炉的理论负载时间;
Ikt:砂型铸造保温电炉的理论空载时间;
Oi’j’:加工设备h上所加工的工序Oij的上一道相邻工序;
CTi:铸件Ci的完工时间;
Psh:加工设备Eh的待机功率;
Pkh:加工设备Eh的空载功率;
Plh:加工设备Eh的负载功率;
Psc:除尘设备Dc的待机功率;
Pkc:除尘设备Dc的空载功率;
Plc:除尘设备Dc的负载功率;
Psh’:热处理设备Hh’的待机功率;
Pkh’:热处理设备Hh’的空载功率;
Plh’:热处理设备Hh’的负载功率;
Th'l:热处理设备Hh’单位负载重量与单位时间下的天然气排放量(单位:m3/(kg·h))
Pki:铸造保温电炉的空载功率;
Pli:铸造保温电炉的负载功率;
Pks:运输设备行车的空载功率;
Pls:运输设备行车的负载功率;
Ee:加工设备(包含热处理设备)总电能消耗量,单位kw.h;
Ed:除尘设备电能总消耗量,单位kw.h;
Ets:运输设备电能总消耗量,单位kw.h;
Ei:铸造保温电炉电能总消耗量,单位kw.h;
Eht:热处理设备天然气总消耗量,单位m3;
SCe:电力标准折煤系数,0.1229kgce/(kW·h);
SCt:天然气标准折煤系数,1.2143kgce/m3;
EFe:电能的碳排放系数,4.035kgCO2e/kgce;
EFt:天然气的碳排放系数,1.744kgCO2e/kgce;
CBz:铸造设备总碳排放量,单位kgCO2e;
Tt:铸造设备总负荷,单位h;
决策变量:
Xijh:0-1变量,工序Oij在加工设备h上加工时则Xijh=1,否则为0;
Xij”h’:0-1变量,工序Oij”在热处理设备h’上加工则Xij”h’=1,否则为0;
Xijd:0-1变量,工序Oij存在除尘操作则Xijd=1,否则为0;
Xijc:0-1变量,工序Oij在除尘设备c上除尘则Xijc=1,否则为0;
Xijs:0-1变量,工序Oij与其下一道相邻工序之间在运输设备c上运输则Xijs=1,否则为0。
5.根据权利要求1所述的基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于:
本发明采用MSOS整数编码方式,该编码方式由两段组成:加工设备选择段(Xm)和工序排序段(X0),两段长度均为D,它们一起构成了种群青蛙个体的位置信息。
6.根据权利要求5所述的基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于:
所述种群初始化中采用基于运输时间的最短工作时间法策略和随机选择初始化加工设备分配部分,工序排序部分由随机方法生成。
7.根据权利要求6所述的基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于:
所述计算Pareto非支配排序与个体拥挤距离具体过程如下:
Pareto排序方法由以下两部分组成,第一部分为所有解设置了两个变量np和Sp,其中P=1,2,…N,N为种群大小,np用来记录所有解中支配解P的个数,Sp用来记录所有解中被P支配的解集,第二部分根据第一部分排序后的结果进行分层,并初始化分层序号为1,首先将种群中所有np为0的个体移除种群,将当前非支配排序分层的序号赋予这些个体;然后将这些个体Sp中所对应的个体np减1,并将分层序号递增加1,如此反复循环执行直到青蛙种群中所有个体的分层序号被赋予;
拥挤距离是对同一排序分层内个体在每一个目标上和它相邻两个体距离之差进行绝对值求和,个体i在第k个目标fk上的拥挤距离为m为目标的个数是个体i在第k个目标上相邻两个体i-1和i+1的目标值,个体i的拥挤距离di如下式所示:
8.根据权利要求7所述的基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于:
采用锦标赛法构建子种群具体过程如下:
随机选择初始种群中的两个个体X1和X2,如果X1支配X2则让X1进入第一个子种群,同时将X2放回初始种群中;若两个个体彼此非受支配关系则随机选择一个个体进入第一个种群并将未选中个体放回初始种群;随机选择两个个体重复以上过程依次构建所以子种群直到所有子种群均被构建完成。
9.根据权利要求8所述的基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法,其特征在于:
在子种群每次迭代进化过程中进行Pareto排序,随机选择某个rank值最大的个体作为子种群最劣解Xw,然后基于子种群最优解Xb对子种群最劣解Xw的工序排序段和加工设备排序段分别进行POX交叉和RPX交叉若交叉后产生的新解Xnew支配Xw,则将Xw用Xnew替换,否则基于全局最优解Xg对子种群最劣解Xw进行POX交叉和RPX交叉,若交叉产生的新解Xnew不被Xw支配,则将Xw用Xnew替换,否则对Xw执行翻转变异,将Xw用变异产生的新解Xnew替代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911200006.2A CN111026051B (zh) | 2019-11-29 | 2019-11-29 | 一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911200006.2A CN111026051B (zh) | 2019-11-29 | 2019-11-29 | 一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111026051A CN111026051A (zh) | 2020-04-17 |
CN111026051B true CN111026051B (zh) | 2021-01-01 |
Family
ID=70207094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911200006.2A Active CN111026051B (zh) | 2019-11-29 | 2019-11-29 | 一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111026051B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112819342B (zh) * | 2021-02-04 | 2021-09-07 | 中国水利水电科学研究院 | 一种碳排放量核算方法、装置、设备及存储介质 |
CN113139720B (zh) * | 2021-04-19 | 2023-10-13 | 武汉理工大学 | 一种具有学习效应的铸造车间生产调度优化方法 |
CN117252401B (zh) * | 2023-11-17 | 2024-02-02 | 北京太极信息***技术有限公司 | 一种资源调度方法、装置、电子设备及存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6288727B1 (en) * | 1997-03-31 | 2001-09-11 | Peter Akemann | Method and apparatus for providing realistic movement of human characters on a video display |
CN103425840A (zh) * | 2013-08-14 | 2013-12-04 | 西北工业大学 | 一种基于改进多目标蛙跳算法的协同空战火力分配方法 |
CN105243458A (zh) * | 2015-11-10 | 2016-01-13 | 河海大学 | 一种基于多目标混合蛙跳差分算法的水库调度方法 |
CN106445881A (zh) * | 2016-09-21 | 2017-02-22 | 景德镇陶瓷大学 | 一种基于多层次信息反馈的多目标混合蛙跳算法 |
CN107102552A (zh) * | 2017-06-12 | 2017-08-29 | 合肥工业大学 | 基于混合集合蛙跳与变邻域算法的平行机调度方法及*** |
CN108776845A (zh) * | 2018-05-07 | 2018-11-09 | 北京创源微致软件有限公司 | 一种基于双目标作业车间调度的混合果蝇算法 |
CN109886589A (zh) * | 2019-02-28 | 2019-06-14 | 长安大学 | 一种基于改进鲸鱼优化算法求解低碳车间调度的方法 |
CN110417061A (zh) * | 2019-07-24 | 2019-11-05 | 东北大学 | 一种基于改进蛙跳算法的电热联合***调度方法 |
-
2019
- 2019-11-29 CN CN201911200006.2A patent/CN111026051B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6288727B1 (en) * | 1997-03-31 | 2001-09-11 | Peter Akemann | Method and apparatus for providing realistic movement of human characters on a video display |
CN103425840A (zh) * | 2013-08-14 | 2013-12-04 | 西北工业大学 | 一种基于改进多目标蛙跳算法的协同空战火力分配方法 |
CN105243458A (zh) * | 2015-11-10 | 2016-01-13 | 河海大学 | 一种基于多目标混合蛙跳差分算法的水库调度方法 |
CN106445881A (zh) * | 2016-09-21 | 2017-02-22 | 景德镇陶瓷大学 | 一种基于多层次信息反馈的多目标混合蛙跳算法 |
CN107102552A (zh) * | 2017-06-12 | 2017-08-29 | 合肥工业大学 | 基于混合集合蛙跳与变邻域算法的平行机调度方法及*** |
CN108776845A (zh) * | 2018-05-07 | 2018-11-09 | 北京创源微致软件有限公司 | 一种基于双目标作业车间调度的混合果蝇算法 |
CN109886589A (zh) * | 2019-02-28 | 2019-06-14 | 长安大学 | 一种基于改进鲸鱼优化算法求解低碳车间调度的方法 |
CN110417061A (zh) * | 2019-07-24 | 2019-11-05 | 东北大学 | 一种基于改进蛙跳算法的电热联合***调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111026051A (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111026051B (zh) | 一种基于改进蛙跳智能算法的柔性离散制造流水车间低碳调度方法 | |
CN109190857B (zh) | 一种基于多目标资源受限项目调度模型的优化算法 | |
Zhou et al. | A multi-objective differential evolution algorithm for parallel batch processing machine scheduling considering electricity consumption cost | |
Zolfaghari et al. | Comparative study of simulated annealing, genetic algorithms and tabu search for solving binary and comprehensive machine-grouping problems | |
Cao et al. | An efficient scheduling approach for an iron-steel plant equipped with self-generation equipment under time-of-use electricity tariffs | |
CN108805403A (zh) | 一种基于改进遗传算法的单件车间调度方法 | |
CN112947319B (zh) | 一种钢铁企业冷轧区多产线的集批排程优化方法及*** | |
CN113805545B (zh) | 一种考虑批处理的柔性流水车间组合调度规则生成方法 | |
CN112668901B (zh) | 一种考虑能耗的钢厂生产调度方法及*** | |
Li et al. | A Production Planning Model for Make‐to‐Order Foundry Flow Shop with Capacity Constraint | |
CN115145235A (zh) | 一种铸造全流程的多目标智能调度方法 | |
CN110135752A (zh) | 一种具有切换时间的成套订单的调度方法 | |
CN113139720B (zh) | 一种具有学习效应的铸造车间生产调度优化方法 | |
CN114298567B (zh) | 连铸机浇次计划排程及开浇时间动态决策方法及*** | |
CN110008531A (zh) | 一种含间隙原子的无序固溶材料原子结构建模方法 | |
CN112836353A (zh) | 一种面向批次生产的锻造调度方法 | |
CN111445050B (zh) | 一种特种合金铝锭生产的组炉方法 | |
Duda et al. | Optimization methods for lot-sizing problem in an automated foundry | |
Zhang et al. | Multi-objective job shop scheduling using a modified artificial bee colony algorithm | |
CN109213094A (zh) | 一种基于多工厂间炼钢-连铸生产钢坯过程的优化调度方法 | |
CN115685920B (zh) | 基于改进布谷鸟算法的铸造混合流水车间调度方法及*** | |
CN115356992A (zh) | 基于多智能体的铸轧界面物流仿真***、方法、计算机设备 | |
Xu et al. | Optimization Model and Algorithm for Ingredient Optimization of Copper Strip Production | |
Liu et al. | Discrete differential evolution algorithm for the job shop scheduling problem | |
He et al. | Research on modeling and solving algorithm of steelmaking continuous casting production scheduling |
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 |