CN109858780A - 一种炼钢-连铸生产调度优化方法 - Google Patents
一种炼钢-连铸生产调度优化方法 Download PDFInfo
- Publication number
- CN109858780A CN109858780A CN201910030357.7A CN201910030357A CN109858780A CN 109858780 A CN109858780 A CN 109858780A CN 201910030357 A CN201910030357 A CN 201910030357A CN 109858780 A CN109858780 A CN 109858780A
- Authority
- CN
- China
- Prior art keywords
- heat
- time
- continuous casting
- model
- lagrange multiplier
- 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
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/30—Computing systems specially adapted for manufacturing
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- General Factory Administration (AREA)
Abstract
本发明涉及工业过程计划与调度优化技术领域,提供一种炼钢‑连铸生产调度优化方法,首先获取炼钢‑连铸生产数据,然后建立炼钢‑连铸生产调度优化模型,并建立拉格朗日乘子松弛模型、每个炉次的带有拉格朗日乘子的计划调度优化子问题模型及对偶问题模型,接着设置拉格朗日乘子初始值,采用后向动态规划法,求解子问题模型,判断最优决策变量值是否满足代理最优条件,然后设置拉格朗日乘子更新的梯度方向,并判断相邻两次迭代的梯度方向所成的角度是否为锐角,再计算更新步长,更新拉格朗日乘子并在其满足停止准则时构造可行的炼钢‑连铸生产调度优化方案。本发明能够缩小调度规模、降低求解难度、提高求解质量及求解速度。
Description
技术领域
本发明涉及工业过程计划与调度优化技术领域,特别是涉及一种炼钢-连铸生产调度优化方法。
背景技术
在工业生产中,炼钢-连铸生产调度,就是在炉次的工艺路径即每个炉次加工工序总数及每道工序的所选设备种类已知的情况下,在保证连铸机连续浇铸过程中每个浇次内的炉次之间及同一设备上两相邻炉次之间不产生作业冲突的前提下,以浇次在连铸机上的开浇时间与该浇次的理想开浇时间之差及相邻工序之间的炉次等待时间为节能降耗因素,以炉次在任意加工时间内的设备唯一性和炉次在任意加工设备内的时间连续性为约束条件,确定每个炉次在已知工艺路径下的每个工序上对应的加工处理设备及该设备对应的加工开始时间,形成炼钢-连铸生产作业时间表,即钢厂的炼钢-连铸计划调度。炼钢-连铸生产调度对炼钢-连铸生产乃至企业的生产效益具有重要意义。
炼钢-连铸生产调度优化问题属于NP(Non-deterministic Polynomial)即非确定性多项式难题,其目标解的搜索涉及解空间的组合爆炸。现有的炼钢-连铸生产调度优化方法主要有运筹学方法、启发式方法和人工智能方法。经典的运筹学方法受限于生产过程的规模而难以在短时间得到高质量的调度方案。启发式方法虽然具有描述简单且适用于实际调度环境的特点,但是对于大规模的复杂多目标调度优化问题的求解效果不理想。人工智能方法很难完全由计算机处理实际调度生产过程结构不良的问题,主要表现在系统中的知识不完备及不一致,从而导致调度优化速度慢。从而,上述现有方法难以保证炼钢-连铸生产调度优化问题的求解质量,而且在短时间内难以获得满足实际生产要求的近似优化可行解。
发明内容
针对现有技术存在的问题,本发明提供一种炼钢-连铸生产调度优化方法,能够缩小调度规模、降低调度问题求解难度、提高求解质量及求解速度,从而提高企业的生产效益。
本发明的技术方案为:
一种炼钢-连铸生产调度优化方法,其特征在于,包括下述步骤:
步骤1:获取钢厂的炼钢-连铸生产数据,所述炼钢-连铸生产数据包括工序总数J、第t时刻第j个工序的可用设备数Mjt、第m台连铸机上已指派的浇次有序集合Bm、第k个浇次中的炉次有序集合Ωk、第i个炉次在第j个工序的加工开始时间ti,j、第i个炉次在第j个工序的作业时间Pi,j、炉次在第j个工序和第j+1个工序之间的标准传输时间Tj,j+1、相邻浇次间更换结晶器所需最小间隔时间Su、炉次在连铸工序的等待时间惩罚系数W1、炉次在精炼工序的等待时间惩罚系数W2、实际开浇时间晚于理想开浇时间时的惩罚系数W3、实际开浇时间早于理想开浇时间时的惩罚系数W4;
其中,j=1,2,3,…,J,第1个工序为炼钢工序,第2至J-1个工序为精炼工序的子工序,第J个工序为连铸工序,Mjt≥1;第J个工序的设备为连铸机,第t时刻第J个工序的设备数即连铸机数为MJt=MJ,MJ为常数;t∈{0,1,2,…,T},ti,j∈{0,1,2,…,T},T为第J个工序的截止时间;
Bm={b(m-1)+1,b(m-1)+2,…,b(m)},b(m)为第m台连铸机上的末浇次序号,b(m)=b(m-1)+|Bm|,|Bm|为第m台连铸机上的浇次个数,b(0)=0,m=1,2,3,…,MJ,b(MJ)=K,K为浇次总数;
Ωk={s(k-1)+1,s(k-1)+2,…,s(k)},s(k)为第k个浇次中的末炉次序号,s(k)=s(k-1)+|Ωk|,|Ωk|为第k个浇次中的炉次个数,s(0)=0,k=1,2,3,…,K,s(K)=|Ω|,|Ω|为炉次总数,Ω为所有炉次的集合,Ω=Ω1∪Ω2∪…∪Ωk∪…∪ΩK, 且k1≠k2;i∈Ω;
步骤2:建立炼钢-连铸生产调度优化模型,具体包括下述步骤:
步骤2.1:建立加权所有炉次在连铸工序的等待时间惩罚之和模型:
步骤2.2:建立加权所有炉次在精炼工序的等待时间惩罚之和模型:
步骤2.3:建立加权所有浇次在连铸工序的开浇时间惩罚之和模型:
其中,当ti,J>di时,x1i=ti,J-di,x2i=0;当ti,J<di时,x1i=0,x2i=di-ti,J;当ti,J=di时,x1i=0,x2i=0;ti,J为第i个炉次在连铸工序的加工开始时间也即实际开浇时间,di为第i个炉次在连铸工序的理想开浇时间;
步骤2.4:建立目标函数模型:
min Z=Z1+Z2+Z3
步骤2.5:建立相邻工序-炉次加工顺序约束模型,即同一炉次在前一工序加工完毕并运达到下一工序后,才能开始加工:
ti,j+1-Tj,j+1-Pi,j-ti,j≥0 i∈Ω,j=1,2…,J-1
步骤2.6:建立浇次加工顺序约束模型,即相邻浇次间需要时间来更换结晶器:
ti+1,J-Pi,J-ti,J≥Su
其中,i=s(k),k∈{b(m-1)+1,b(m-1)+2,…,b(m)-1},m=1,2,3,…,MJ;
步骤2.7:建立连铸工序-炉次加工顺序约束模型,即同一浇次的相邻炉次有加工前后顺序:
ti+1,J-Pi,J-ti,J≥0
其中,i∈{s(k-1)+1,s(k-1)+2,…,s(k)-1},k∈{b(m-1)+1,b(m-1)+2,…,b(m)},m=1,2,3,…,MJ;
步骤2.8:建立设备能力约束模型,即在同一设备上不能同时加工多个炉次:
∑i∈Ωδijt≤Mjt j=1,2…,J-1,t=1,2,…,T
其中,当第t时刻第i个炉次在第j个工序被加工时,δijt=1;当第t时刻第i个炉次在第j个工序没有被加工时,δijt=0;
步骤2.9:上述步骤2.4-步骤2.8建立的所有模型共同构成炼钢-连铸生产调度优化模型;
步骤3:建立拉格朗日乘子松弛“不断浇、不冲突”约束即“加工顺序约束、设备能力约束”模型;
其中,ui为第i个炉次的拉格朗日乘子,vjt为第t时刻第j个工序的拉格朗日乘子;
步骤4:建立第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型:
其中,
步骤5:建立对偶问题模型;
步骤6:用n表示迭代次数,n的初始值设为0,第n次迭代中拉格朗日乘子ui、vjt的值分别为设置拉格朗日乘子ui、vjt的初始值分别为其中,i∈Ω,j=1,2…,J,t=0,1,…,T;
步骤7:令i的初始值为1,采用后向动态规划法,求解步骤4中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型,得到第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的最优决策变量值
步骤8:判断步骤7中获得的第n次迭代中第i个炉次的最优决策变量值是否满足代理最优条件若满足,则进入步骤9;若不满足,则令i=i+1,返回步骤7,继续求解剩余炉次的带有拉格朗日乘子的计划调度优化子问题模型,直到满足代理最优条件;
步骤9:设置第n次迭代中拉格朗日乘子在第j个工序中更新所需的梯度方向分别为
步骤10:判断相邻两次迭代的梯度方向所成的角度,如果为锐角,则进入步骤11;如果为钝角,则引入拉格朗日乘子的偏移次梯度分别为 使相邻两次迭代的梯度方向为锐角;其中,
步骤11:计算第n次迭代中拉格朗日乘子的更新步长为
其中,
步骤12:基于梯度方向和更新步长,来更新拉格朗日乘子:
步骤13:检验拉格朗日乘子是否满足停止准则||Hn+1-Hn||<ε1,如果满足停止准则,则停止拉格朗日乘子更新,得到最优决策即单炉次调度优化方案,进入步骤14;如果不满足停止准则,则令n=n+1,返回步骤7,进入下一次迭代计算;其中, ε1为常数,ε1>0;
步骤14:根据单炉次调度优化方案,获取炼钢-连铸生产调度优化方案,并判断炼钢-连铸生产调度优化方案是否有违反所述步骤3中“加工顺序约束、设备能力约束”的情况出现,如果有,则按照改变加工开始时间代价规则进一步优化炼钢-连铸生产调度优化方案;如果没有,则所获取的炼钢-连铸生产调度优化方案为可行方案,具体包括下述步骤:
步骤14.1:从时间节点0开始,判断所有炉次是否都满足“加工顺序约束、设备能力约束”,如果都满足,则进入步骤14.4;如果不都满足,则进入步骤14.2;
步骤14.2:在不满足“加工顺序约束、设备能力约束”的炉次中,分别计算每个炉次的加工开始时间节点推迟一个单位的代价,对代价最小的炉次的加工开始时间推迟一个单位;
步骤14.3:重复步骤14.1-14.2,直到该时间节点所有炉次都满足“加工顺序约束、设备能力约束”;
步骤14.4:判断下一时间节点是否达到最后时间,如果达到最后时间,则得到可行的炼钢-连铸生产调度优化方案;如果未达到最后时间,则转入步骤14.1,在下一时间节点对所有炉次进行“加工顺序约束、设备能力约束”的判断。
所述步骤7包括下述步骤:
步骤7.1:令i的初始值为1、j的初始值为J,将第n次迭代中决策变量ti,j、δijt的值分别表示为将第n次迭代中第j个工序的惩罚值函数Fij(tij,δijt)表示为在第j个工序中,在第n次迭代中决策变量的可行值中找出第j个工序的最优决策变量值使得第n次迭代中第j个工序的惩罚值函数最小,得到第n次迭代中第j个工序的最优惩罚值为其中,第j个工序的最优决策变量值和最优惩罚值可通过下式得到:Fij*(tij,δijt)=min{Sij(tij,δijt)+ρj*Fi,j+1*(ti,j+1,δi,j+1,t)};Sij(tij,δijt)为在决策变量tij和δijt下第j个工序的当前惩罚值, 其中,当j=J时,ρj=0、αj=1,否则ρj=1、αj=0;
步骤7.2:令j=j-1,反向迭代计算第j个工序的最优决策变量值和最优惩罚值,直到第1个工序;第1个工序的最优惩罚值即为第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的目标函数最优值,最优惩罚值对应的最优决策变量值即为第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的最优解。
本发明的有益效果为:
本发明基于高效拉氏松弛迭代算法进行炼钢-连铸生产调度的优化,能够缩小调度规模、降低调度问题求解难度、提高求解质量及求解速度,尤其适用于大规模炼钢-连铸生产过程的多工序、多耦合、多约束、间歇与连续作业方式相混杂的流程式特点。本发明在寻找最优拉格朗日乘子的过程中,采用无需预估最优值的基于梯度方向可控的迭代优化策略,避免了传统方法在可行域内部搜索时出现锯齿震荡的问题,并在保证求解质量的前提下,提高了收敛速度,能够有效提高炼钢-连铸生产过程调度编制的求解质量以及响应速度,提高企业的生产效益。
附图说明
图1为本发明的炼钢-连铸生产调度优化方法的流程图。
具体实施方式
下面将结合附图和具体实施方式,对本发明作进一步描述。
本发明的目的是提供一种炼钢-连铸生产调度优化方法,能够缩小调度规模、降低调度问题求解难度、提高求解质量及求解速度,从而提高企业的生产效益。
如图1所示,为本发明的炼钢-连铸生产调度优化方法的流程图。
本发明的炼钢-连铸生产调度优化方法,其特征在于,包括下述步骤:
步骤1:获取钢厂的炼钢-连铸生产数据,所述炼钢-连铸生产数据包括工序总数J、第t时刻第j个工序的可用设备数Mjt、第m台连铸机上已指派的浇次有序集合Bm、第k个浇次中的炉次有序集合Ωk、第i个炉次在第j个工序的加工开始时间ti,j、第i个炉次在第j个工序的作业时间Pi,j、炉次在第j个工序和第j+1个工序之间的标准传输时间Tj,j+1、相邻浇次间更换结晶器所需最小间隔时间Su、炉次在连铸工序的等待时间惩罚系数W1、炉次在精炼工序的等待时间惩罚系数W2、实际开浇时间晚于理想开浇时间时的惩罚系数W3、实际开浇时间早于理想开浇时间时的惩罚系数W4;
其中,j=1,2,3,…,J,第1个工序为炼钢工序,第2至J-1个工序为精炼工序的子工序,第J个工序为连铸工序,Mjt≥1;第J个工序的设备为连铸机,第t时刻第J个工序的设备数即连铸机数为MJt=MJ,MJ为常数;t∈{0,1,2,…,T},ti,j∈{0,1,2,…,T},T为第J个工序的截止时间;
Bm={b(m-1)+1,b(m-1)+2,…,b(m)},b(m)为第m台连铸机上的末浇次序号,b(m)=b(m-1)+|Bm|,|Bm|为第m台连铸机上的浇次个数,b(0)=0,m=1,2,3,…,MJ,b(MJ)=K,K为浇次总数;
Ωk={s(k-1)+1,s(k-1)+2,…,s(k)},s(k)为第k个浇次中的末炉次序号,s(k)=s(k-1)+|Ωk|,|Ωk|为第k个浇次中的炉次个数,s(0)=0,k=1,2,3,…,K,s(K)=|Ω|,|Ω|为炉次总数,Ω为所有炉次的集合,Ω=Ω1∪Ω2∪…∪Ωk∪…∪ΩK, 且k1≠k2;i∈Ω。
本实施例中,第i个炉次在第j个工序的作业时间Pi,j为在区间[2,6]上满足均匀分布的整数;炉次在第j个工序和第j+1个工序之间的标准传输时间Tj,j+1为在区间[1,5]上分布的整数;相邻浇次间更换结晶器所需最小间隔时间Su=5;惩罚系数W1=10,W2=10,W3=15,W4=15。
步骤2:建立炼钢-连铸生产调度优化模型,具体包括下述步骤:
步骤2.1:建立加权所有炉次在连铸工序的等待时间惩罚之和模型:
步骤2.2:建立加权所有炉次在精炼工序的等待时间惩罚之和模型:
步骤2.3:建立加权所有浇次在连铸工序的开浇时间惩罚之和模型:
其中,当ti,J>di时,x1i=ti,J-di,x2i=0;当ti,J<di时,x1i=0,x2i=di-ti,J;当ti,J=di时,x1i=0,x2i=0;ti,J为第i个炉次在连铸工序的加工开始时间也即实际开浇时间,di为第i个炉次在连铸工序的理想开浇时间;
步骤2.4:建立目标函数模型:
min Z=Z1+Z2+Z3
步骤2.5:建立相邻工序-炉次加工顺序约束模型,即同一炉次在前一工序加工完毕并运达到下一工序后,才能开始加工:
ti,j+1-Tj,j+1-Pi,j-ti,j≥0 i∈Ω,j=1,2…,J-1
步骤2.6:建立浇次加工顺序约束模型,即相邻浇次间需要时间来更换结晶器:
ti+1,J-Pi,J-ti,J≥Su
其中,i=s(k),k∈{b(m-1)+1,b(m-1)+2,…,b(m)-1},m=1,2,3,…,MJ;
步骤2.7:建立连铸工序-炉次加工顺序约束模型,即同一浇次的相邻炉次有加工前后顺序:
ti+1,J-Pi,J-ti,J≥0
其中,i∈{s(k-1)+1,s(k-1)+2,…,s(k)-1},k∈{b(m-1)+1,b(m-1)+2,…,b(m)},m=1,2,3,…,MJ;
步骤2.8:建立设备能力约束模型,即在同一设备上不能同时加工多个炉次:
∑i∈Ωδijt≤Mjt j=1,2…,J-1,t=1,2,…,T
其中,当第t时刻第i个炉次在第j个工序被加工时,δijt=1;当第t时刻第i个炉次在第j个工序没有被加工时,δijt=0;
步骤2.9:上述步骤2.4-步骤2.8建立的所有模型共同构成炼钢-连铸生产调度优化模型。
步骤3:建立拉格朗日乘子松弛“不断浇、不冲突”约束即“加工顺序约束、设备能力约束”模型;
其中,ui为第i个炉次的拉格朗日乘子,vjt为第t时刻第j个工序的拉格朗日乘子;
步骤4:建立第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型:
其中,
步骤5:建立对偶问题模型;
至此,模型建立完成,从步骤6开始,通过采用拉氏框架下无需预估最优值的基于梯度方向可控的迭代优化策略的乘子寻优方法,对炼钢-连铸生产调度优化问题进行求解,具体如下:
步骤6:用n表示迭代次数,n的初始值设为0,第n次迭代中拉格朗日乘子ui、vjt的值分别为设置拉格朗日乘子ui、vjt的初始值分别为其中,i∈Ω,j=1,2…,J,t=0,1,…,T;
步骤7:令i的初始值为1,采用后向动态规划法,求解步骤4中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型,得到第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的最优决策变量值具体步骤如下:
步骤7.1:令i的初始值为1、j的初始值为J,将第n次迭代中决策变量ti,j、δijt的值分别表示为将第n次迭代中第j个工序的惩罚值函数Fij(tij,δijt)表示为在第j个工序中,在第n次迭代中决策变量的可行值中找出第j个工序的最优决策变量值使得第n次迭代中第j个工序的惩罚值函数最小,得到第n次迭代中第j个工序的最优惩罚值为其中,第j个工序的最优决策变量值和最优惩罚值可通过下式得到:Fij*(tij,δijt)=min{Sij(tij,δijt)+ρj*Fi,j+1*(ti,j+1,δi,j+1,t)};Sij(tij,δijt)为在决策变量tij和δijt下第j个工序的当前惩罚值,αj;其中,当j=J时,ρj=0、αj=1,否则ρj=1、αj=0;
步骤7.2:令j=j-1,反向迭代计算第j个工序的最优决策变量值和最优惩罚值,直到第1个工序;第1个工序的最优惩罚值即为第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的目标函数最优值,最优惩罚值对应的最优决策变量值即为第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的最优解;
步骤8:判断步骤7中获得的第n次迭代中第i个炉次的最优决策变量值是否满足代理最优条件若满足,则进入步骤9;若不满足,则令i=i+1,返回步骤7,继续求解剩余炉次的带有拉格朗日乘子的计划调度优化子问题模型,直到满足代理最优条件;
步骤9:设置第n次迭代中拉格朗日乘子在第j个工序中更新所需的梯度方向分别为
步骤10:判断相邻两次迭代的梯度方向所成的角度,如果为锐角,则进入步骤11;如果为钝角,则引入拉格朗日乘子的偏移次梯度分别为 使相邻两次迭代的梯度方向为锐角;其中,
本实施例中,β=1.2;
步骤11:计算第n次迭代中拉格朗日乘子的更新步长为
其中,本实施例中,Y=1.5,p=0.8;
步骤12:基于梯度方向和更新步长,来更新拉格朗日乘子:
步骤13:检验拉格朗日乘子是否满足停止准则||Hn+1-Hn||<ε1,如果满足停止准则,则停止拉格朗日乘子更新,得到最优决策即单炉次调度优化方案,进入步骤14;如果不满足停止准则,则令n=n+1,返回步骤7,进入下一次迭代计算;其中, ε1为常数,ε1>0;本实施例中,ε1=1e-4;
步骤14:根据单炉次调度优化方案,获取炼钢-连铸生产调度优化方案,并判断炼钢-连铸生产调度优化方案是否有违反所述步骤3中“加工顺序约束、设备能力约束”的情况出现,如果有,则按照改变加工开始时间代价规则进一步优化炼钢-连铸生产调度优化方案;如果没有,则所获取的炼钢-连铸生产调度优化方案为可行方案,具体包括下述步骤:
步骤14.1:从时间节点0开始,判断所有炉次是否都满足“加工顺序约束、设备能力约束”,如果都满足,则进入步骤14.4;如果不都满足,则进入步骤14.2;
步骤14.2:在不满足“加工顺序约束、设备能力约束”的炉次中,分别计算每个炉次的加工开始时间节点推迟一个单位的代价,对代价最小的炉次的加工开始时间推迟一个单位;
步骤14.3:重复步骤14.1-14.2,直到该时间节点所有炉次都满足“加工顺序约束、设备能力约束”;
步骤14.4:判断下一时间节点是否达到最后时间,如果达到最后时间,则得到可行的炼钢-连铸生产调度优化方案;如果未达到最后时间,则转入步骤14.1,在下一时间节点对所有炉次进行“加工顺序约束、设备能力约束”的判断。
本实施例中,通过c#语言编程,在Intel Core i5 5200 CPU 4GB Windows 10/64位操作系统的PC上执行传统代理梯度优化方法(LR-1)以及本发明的炼钢-连铸生产调度优化方法(LR-2),来同时求解不同规模(不同工序总数和炉次总数)的炼钢-连铸生产调度优化问题,用对偶间隙gap和运行时间作为评价两种方法求解质量的依据,具体结果如表1所示。其中,对偶间隙UB为历史上得到的最好上界,LB为代理对偶目标值。
由表1可以看出,本发明的炼钢-连铸生产调度优化方法与传统代理梯度优化方法相比,求解质量相差不多,求解效率优于传统代理梯度优化方法。特别地,对于规模大的炼钢-连铸生产调度优化问题,本发明的方法具有更大的优势。
显然,上述实施例仅仅是本发明的一部分实施例,而不是全部的实施例。上述实施例仅用于解释本发明,并不构成对本发明保护范围的限定。基于上述实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,也即凡在本申请的精神和原理之内所作的所有修改、等同替换和改进等,均落在本发明要求的保护范围内。
表1
Claims (2)
1.一种炼钢-连铸生产调度优化方法,其特征在于,包括下述步骤:
步骤1:获取钢厂的炼钢-连铸生产数据,所述炼钢-连铸生产数据包括工序总数J、第t时刻第j个工序的可用设备数Mjt、第m台连铸机上已指派的浇次有序集合Bm、第k个浇次中的炉次有序集合Ωk、第i个炉次在第j个工序的加工开始时间ti,j、第i个炉次在第j个工序的作业时间Pi,j、炉次在第j个工序和第j+1个工序之间的标准传输时间Tj,j+1、相邻浇次间更换结晶器所需最小间隔时间Su、炉次在连铸工序的等待时间惩罚系数W1、炉次在精炼工序的等待时间惩罚系数W2、实际开浇时间晚于理想开浇时间时的惩罚系数W3、实际开浇时间早于理想开浇时间时的惩罚系数W4;
其中,j=1,2,3,…,J,第1个工序为炼钢工序,第2至J-1个工序为精炼工序的子工序,第J个工序为连铸工序,Mjt≥1;第J个工序的设备为连铸机,第t时刻第J个工序的设备数即连铸机数为MJt=MJ,MJ为常数;t∈{0,1,2,…,T},ti,j∈{0,1,2,…,T},T为第J个工序的截止时间;
Bm={b(m-1)+1,b(m-1)+2,…,b(m)},b(m)为第m台连铸机上的末浇次序号,b(m)=b(m-1)+|Bm|,|Bm|为第m台连铸机上的浇次个数,b(0)=0,m=1,2,3,…,MJ,b(MJ)=K,K为浇次总数;
Ωk={s(k-1)+1,s(k-1)+2,…,s(k)},s(k)为第k个浇次中的末炉次序号,s(k)=s(k-1)+|Ωk|,|Ωk|为第k个浇次中的炉次个数,s(0)=0,k=1,2,3,…,K,s(K)=|Ω|,|Ω|为炉次总数,Ω为所有炉次的集合,Ω=Ω1∪Ω2∪…∪Ωk∪…∪ΩK, 且k1≠k2;i∈Ω;
步骤2:建立炼钢-连铸生产调度优化模型,具体包括下述步骤:
步骤2.1:建立加权所有炉次在连铸工序的等待时间惩罚之和模型:
步骤2.2:建立加权所有炉次在精炼工序的等待时间惩罚之和模型:
步骤2.3:建立加权所有浇次在连铸工序的开浇时间惩罚之和模型:
其中,当ti,J>di时,x1i=ti,J-di,x2i=0;当ti,J<di时,x1i=0,x2i=di-ti,J;当ti,J=di时,x1i=0,x2i=0;ti,J为第i个炉次在连铸工序的加工开始时间也即实际开浇时间,di为第i个炉次在连铸工序的理想开浇时间;
步骤2.4:建立目标函数模型:
min Z=Z1+Z2+Z3
步骤2.5:建立相邻工序-炉次加工顺序约束模型,即同一炉次在前一工序加工完毕并运达到下一工序后,才能开始加工:
ti,j+1-Tj,j+1-Pi,j-ti,j≥0 i∈Ω,j=1,2…,J-1
步骤2.6:建立浇次加工顺序约束模型,即相邻浇次间需要时间来更换结晶器:
ti+1,J-Pi,J-ti,J≥Su
其中,i=s(k),k∈{b(m-1)+1,b(m-1)+2,…,b(m)-1},m=1,2,3,…,MJ;
步骤2.7:建立连铸工序-炉次加工顺序约束模型,即同一浇次的相邻炉次有加工前后顺序:
ti+1,J-Pi,J-ti,J≥0
其中,i∈{s(k-1)+1,s(k-1)+2,…,s(k)-1},k∈{b(m-1)+1,b(m-1)+2,…,b(m)},m=1,2,3,…,MJ;
步骤2.8:建立设备能力约束模型,即在同一设备上不能同时加工多个炉次:
∑i∈Ωδijt≤Mjt j=1,2…,J-1,t=1,2,…,T
其中,当第t时刻第i个炉次在第j个工序被加工时,δijt=1;当第t时刻第i个炉次在第j个工序没有被加工时,δijt=0;
步骤2.9:上述步骤2.4-步骤2.8建立的所有模型共同构成炼钢-连铸生产调度优化模型;
步骤3:建立拉格朗日乘子松弛“不断浇、不冲突”约束即“加工顺序约束、设备能力约束”模型;
其中,ui为第i个炉次的拉格朗日乘子,vjt为第t时刻第j个工序的拉格朗日乘子;
步骤4:建立第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型:
其中,
步骤5:建立对偶问题模型;
步骤6:用n表示迭代次数,n的初始值设为0,第n次迭代中拉格朗日乘子ui、vjt的值分别为设置拉格朗日乘子ui、vjt的初始值分别为其中,i∈Ω,j=1,2…,J,t=0,1,…,T;
步骤7:令i的初始值为1,采用后向动态规划法,求解步骤4中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型,得到第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的最优决策变量值
步骤8:判断步骤7中获得的第n次迭代中第i个炉次的最优决策变量值是否满足代理最优条件若满足,则进入步骤9;若不满足,则令i=i+1,返回步骤7,继续求解剩余炉次的带有拉格朗日乘子的计划调度优化子问题模型,直到满足代理最优条件;
步骤9:设置第n次迭代中拉格朗日乘子在第j个工序中更新所需的梯度方向分别为
步骤10:判断相邻两次迭代的梯度方向所成的角度,如果为锐角,则进入步骤11;如果为钝角,则引入拉格朗日乘子的偏移次梯度分别为 使相邻两次迭代的梯度方向为锐角;其中,
步骤11:计算第n次迭代中拉格朗日乘子的更新步长为
其中,
步骤12:基于梯度方向和更新步长,来更新拉格朗日乘子:
步骤13:检验拉格朗日乘子是否满足停止准则||Hn+1-Hn||<ε1,如果满足停止准则,则停止拉格朗日乘子更新,得到最优决策即单炉次调度优化方案,进入步骤14;如果不满足停止准则,则令n=n+1,返回步骤7,进入下一次迭代计算;其中, ε1为常数,ε1>0;
步骤14:根据单炉次调度优化方案,获取炼钢-连铸生产调度优化方案,并判断炼钢-连铸生产调度优化方案是否有违反所述步骤3中“加工顺序约束、设备能力约束”的情况出现,如果有,则按照改变加工开始时间代价规则进一步优化炼钢-连铸生产调度优化方案;如果没有,则所获取的炼钢-连铸生产调度优化方案为可行方案,具体包括下述步骤:
步骤14.1:从时间节点0开始,判断所有炉次是否都满足“加工顺序约束、设备能力约束”,如果都满足,则进入步骤14.4;如果不都满足,则进入步骤14.2;
步骤14.2:在不满足“加工顺序约束、设备能力约束”的炉次中,分别计算每个炉次的加工开始时间节点推迟一个单位的代价,对代价最小的炉次的加工开始时间推迟一个单位;
步骤14.3:重复步骤14.1-14.2,直到该时间节点所有炉次都满足“加工顺序约束、设备能力约束”;
步骤14.4:判断下一时间节点是否达到最后时间,如果达到最后时间,则得到可行的炼钢-连铸生产调度优化方案;如果未达到最后时间,则转入步骤14.1,在下一时间节点对所有炉次进行“加工顺序约束、设备能力约束”的判断。
2.根据权利要求1所述的炼钢-连铸生产调度优化方法,其特征在于,所述步骤7包括下述步骤:
步骤7.1:令i的初始值为1、j的初始值为J,将第n次迭代中决策变量ti,j、δijt的值分别表示为将第n次迭代中第j个工序的惩罚值函数Fij(tij,δijt)表示为在第j个工序中,在第n次迭代中决策变量的可行值中找出第j个工序的最优决策变量值使得第n次迭代中第j个工序的惩罚值函数最小,得到第n次迭代中第j个工序的最优惩罚值为其中,第j个工序的最优决策变量值和最优惩罚值可通过下式得到:Fij*(tij,δijt)=min{Sij(tij,δijt)+ρj*Fi,j+1*(ti,j+1,δi,j+1,t)};Sij(tij,δijt)为在决策变量tij和δijt下第j个工序的当前惩罚值, 其中,当j=J时,ρj=0、αj=1,否则ρj=1、αj=0;
步骤7.2:令j=j-1,反向迭代计算第j个工序的最优决策变量值和最优惩罚值,直到第1个工序;第1个工序的最优惩罚值即为第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的目标函数最优值,最优惩罚值对应的最优决策变量值即为第n次迭代中第i个炉次的带有拉格朗日乘子的计划调度优化子问题模型的最优解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910030357.7A CN109858780B (zh) | 2019-01-14 | 2019-01-14 | 一种炼钢-连铸生产调度优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910030357.7A CN109858780B (zh) | 2019-01-14 | 2019-01-14 | 一种炼钢-连铸生产调度优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109858780A true CN109858780A (zh) | 2019-06-07 |
CN109858780B CN109858780B (zh) | 2023-01-31 |
Family
ID=66894562
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910030357.7A Active CN109858780B (zh) | 2019-01-14 | 2019-01-14 | 一种炼钢-连铸生产调度优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109858780B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111769555A (zh) * | 2020-07-08 | 2020-10-13 | 山东大学 | 计及工序限制的短流程钢铁企业负荷优化调控方法及系统 |
CN112091192A (zh) * | 2020-09-02 | 2020-12-18 | 沈阳建筑大学 | 一种基于炼钢连铸生产优化调度装置 |
CN112488391A (zh) * | 2020-11-30 | 2021-03-12 | 合肥工业大学 | 一种基于拉格朗日松弛的工业烟草物流调度方法 |
CN115729193A (zh) * | 2022-11-24 | 2023-03-03 | 沈阳建筑大学 | 一种炼钢-连铸中间包批量计划优化方法 |
CN116088457A (zh) * | 2023-04-10 | 2023-05-09 | 清华大学 | 分布式鲁棒联合机会约束模型的炼钢连铸调度方法和装置 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10235540A (ja) * | 1997-02-25 | 1998-09-08 | Nippon Steel Corp | 製造工程のスケジューリングシステム |
CN101339634A (zh) * | 2008-06-18 | 2009-01-07 | 苏志雄 | 考虑钢水驻留时间期限的炼钢连铸生产调度方法 |
US20100314069A1 (en) * | 2009-06-16 | 2010-12-16 | Nucor Corporation | High efficiency plant for making steel |
CN104375478A (zh) * | 2014-09-04 | 2015-02-25 | 太极计算机股份有限公司 | 一种轧钢生产过程产品质量在线预测与优化的方法及装置 |
CN105678332A (zh) * | 2016-01-08 | 2016-06-15 | 昆明理工大学 | 火焰图像cnn识别建模的转炉炼钢终点判断方法及系统 |
CN106126842A (zh) * | 2016-06-29 | 2016-11-16 | 中冶赛迪工程技术股份有限公司 | 一种炼钢连铸车间动态生产运行的仿真方法 |
CN106611221A (zh) * | 2016-12-21 | 2017-05-03 | 重庆大学 | 一种用于解决连铸机故障的炼钢‑连铸重调度方法 |
CN106651003A (zh) * | 2016-11-18 | 2017-05-10 | 沈阳工程学院 | 一种条件偏转近似次梯度的多项式动态规划方法 |
WO2017088674A1 (zh) * | 2015-11-23 | 2017-06-01 | 东北大学 | 一种面向全流程生产的炼钢组批与排产方法 |
CN108376315A (zh) * | 2018-02-11 | 2018-08-07 | 重庆大学 | 一种基于图模型的不确定炼钢-连铸柔性优化调度方法及系统 |
CN108676955A (zh) * | 2018-05-02 | 2018-10-19 | 中南大学 | 一种转炉炼钢终点碳含量和温度控制方法 |
-
2019
- 2019-01-14 CN CN201910030357.7A patent/CN109858780B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10235540A (ja) * | 1997-02-25 | 1998-09-08 | Nippon Steel Corp | 製造工程のスケジューリングシステム |
CN101339634A (zh) * | 2008-06-18 | 2009-01-07 | 苏志雄 | 考虑钢水驻留时间期限的炼钢连铸生产调度方法 |
US20100314069A1 (en) * | 2009-06-16 | 2010-12-16 | Nucor Corporation | High efficiency plant for making steel |
CN104375478A (zh) * | 2014-09-04 | 2015-02-25 | 太极计算机股份有限公司 | 一种轧钢生产过程产品质量在线预测与优化的方法及装置 |
WO2017088674A1 (zh) * | 2015-11-23 | 2017-06-01 | 东北大学 | 一种面向全流程生产的炼钢组批与排产方法 |
CN105678332A (zh) * | 2016-01-08 | 2016-06-15 | 昆明理工大学 | 火焰图像cnn识别建模的转炉炼钢终点判断方法及系统 |
CN106126842A (zh) * | 2016-06-29 | 2016-11-16 | 中冶赛迪工程技术股份有限公司 | 一种炼钢连铸车间动态生产运行的仿真方法 |
CN106651003A (zh) * | 2016-11-18 | 2017-05-10 | 沈阳工程学院 | 一种条件偏转近似次梯度的多项式动态规划方法 |
CN106611221A (zh) * | 2016-12-21 | 2017-05-03 | 重庆大学 | 一种用于解决连铸机故障的炼钢‑连铸重调度方法 |
CN108376315A (zh) * | 2018-02-11 | 2018-08-07 | 重庆大学 | 一种基于图模型的不确定炼钢-连铸柔性优化调度方法及系统 |
CN108676955A (zh) * | 2018-05-02 | 2018-10-19 | 中南大学 | 一种转炉炼钢终点碳含量和温度控制方法 |
Non-Patent Citations (4)
Title |
---|
KUN MAO ET AL.: "Lagrangian heuristic for scheduling a steelmaking-continuous casting process", 《 2013 IEEE SYMPOSIUM ON COMPUTATIONAL INTELLIGENCE IN SCHEDULING (CISCHED)》 * |
KUNKUNPENG ET AL.: "An improved artificial bee colony algorithm for steelmaking–refining–continuous casting scheduling problem", 《CHINESE JOURNAL OF CHEMICAL ENGINEERING》 * |
LIANGLIANG SUN ET AL.: "Research on Scheduling of Iron and Steel Scrap Steelmaking and Continuous Casting Process Aiming at Power Saving and Carbon Emissions Reducing", 《IEEE ROBOTICS AND AUTOMATION LETTERS》 * |
韩大勇等: "基于拉格朗日下界求解的炼钢", 《武汉科技大学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111769555A (zh) * | 2020-07-08 | 2020-10-13 | 山东大学 | 计及工序限制的短流程钢铁企业负荷优化调控方法及系统 |
CN111769555B (zh) * | 2020-07-08 | 2021-04-30 | 山东大学 | 计及工序限制的短流程钢铁企业负荷优化调控方法及系统 |
CN112091192A (zh) * | 2020-09-02 | 2020-12-18 | 沈阳建筑大学 | 一种基于炼钢连铸生产优化调度装置 |
CN112488391A (zh) * | 2020-11-30 | 2021-03-12 | 合肥工业大学 | 一种基于拉格朗日松弛的工业烟草物流调度方法 |
CN112488391B (zh) * | 2020-11-30 | 2022-10-04 | 合肥工业大学 | 一种基于拉格朗日松弛的工业烟草物流调度方法 |
CN115729193A (zh) * | 2022-11-24 | 2023-03-03 | 沈阳建筑大学 | 一种炼钢-连铸中间包批量计划优化方法 |
CN116088457A (zh) * | 2023-04-10 | 2023-05-09 | 清华大学 | 分布式鲁棒联合机会约束模型的炼钢连铸调度方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109858780B (zh) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109858780A (zh) | 一种炼钢-连铸生产调度优化方法 | |
CN111966050B (zh) | 基于ammas-ga嵌套算法的双资源模具作业车间调度优化方法 | |
CN108376315B (zh) | 一种基于图模型的不确定炼钢-连铸柔性优化调度方法及系统 | |
CN113705866B (zh) | 基于资源受限项目调度问题模型的排程优化方法和系统 | |
Kim et al. | Synchronized production planning and scheduling in semiconductor fabrication | |
CN111259314B (zh) | 一种分布式混凝土预制构件流水车间生产调度方法 | |
Wang et al. | Synchronizing production scheduling with resources allocation for precast components in a multi-agent system environment | |
Pan et al. | Steelmaking-continuous casting scheduling problem with multi-position refining furnaces under time-of-use tariffs | |
CN115204497A (zh) | 一种基于强化学习的预制构件生产调度优化方法及系统 | |
CN101344781B (zh) | 一种冷轧热镀锌机组板卷自动调度方法及系统 | |
Geng et al. | Research on multi-objective hybrid flow shop scheduling problem with dual resource constraints using improved memetic algorithm | |
Dai et al. | Resource levelling in repetitive construction projects with interruptions: an integrated approach | |
Liu et al. | An integrated scheduling approach considering dispatching strategy and conflict-free route of AMRs in flexible job shop | |
Chen et al. | Digital twin-oriented collaborative optimization of fuzzy flexible job shop scheduling under multiple uncertainties | |
Shin et al. | A Study on Method of Vertical Zoning of Construction Lift for High-rise Building based on Lift Planning & Operation History Database | |
Gou et al. | Improved shuffled frog leaping algorithm for multi-objection flexible job-shop scheduling problem | |
Shidong et al. | Multi-objectives optimization on flights landing sequence at busy airport | |
Li | Multi-ant colony-based sequencing method for semiconductor wafer fabrication facilities with multi-bottleneck | |
Sun et al. | Research on steelmaking-continuous casting production scheduling problem with uncertain processing time based on Lagrangian relaxation framework | |
Zhang et al. | Mathematical modelling and a discrete cuckoo search particle swarm optimization algorithm for mixed model sequencing problem with interval task times | |
Liang et al. | Optimization of precast production scheduling for rail track slab | |
CN113642763B (zh) | 基于预算约束的高端装备研制资源配置与优化调度方法 | |
CN108985605B (zh) | 用于施工计划编制的进度管控方法 | |
CN114819321B (zh) | 一种面向分布式机器学习的参数传输通信优化方法 | |
CN114827142B (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 |