发明名称 一种连铸机钢水浇铸顺序的控制方法
摘要 本发明提供一种连铸机钢水浇铸顺序的控制方法,从炼钢连铸生产车间的三级制造执行系统获取生产数据;建立数学模型,描述生产过程的工艺和钢水浇铸顺序的控制方案;以基本工艺约束为依据,制定钢水浇铸顺序的预控制方案;优化制定的钢水浇铸顺序的预控制方案,将其转化为生产指令下发给二级过程控制系统,二级过程控制系统执行生产。本发明通过数学模型描述连铸生产过程中的工艺特征,采用智能优化算法决策钢水浇铸顺序,使得连铸生产系统处于最优的工作状态,从而保证机组稳定运行,均衡生产,提高产品质量和钢水收得率,提高设备利用率,缩短生产周期和提高生产效率。
申请公布号 CN102266927B 申请公布日期 2013.02.13
申请号 CN201110218618.1 申请日期 2011.08.01
申请人 东北大学 发明人 唐立新;汪恭书
分类号 B22D11/18(2006.01)I 主分类号 B22D11/18(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 梁焱
主权项 1.一种连铸机钢水浇铸顺序的控制方法,其特征在于:包括如下步骤:步骤1:获取生产过程工艺数据;步骤2:建立生产过程的工艺特征和钢水浇铸顺序联系的描述方式;步骤3:以基本工艺约束为依据,制定钢水浇铸顺序的预控制方案;步骤4:以描述方式为依据,优化连铸生产过程各项工艺指标,改进预控制方案,设定钢水浇铸顺序控制方案;所述步骤1获取生成过程工艺数据,包括:从炼钢连铸生产车间的三级制造执行系统MES获取以下数据项,每炉钢水的炉次号、钢级、首宽、尾宽、铸坯厚度、钢水精炼方式、钢水重量,首宽为第一块铸坯宽度,尾宽为最后一块铸坯宽度;从炼钢连铸生产车间的二级过程控制系统获取工艺数据,包括连铸机数目,连铸机结晶器厚度范围,连铸机结晶器宽度范围及在线调整规则,中间包寿命和结晶器寿命,钢级拉速度;所述步骤2建立生产过程的工艺特征和钢水浇铸顺序联系的描述方式,以数学模型描述生产过程工艺特征与钢水浇铸顺序的联系,具体包括:1)将生产数据和工艺数据以及数据间的关联关系映射为数学模型的参数,其具体过程如下:1.1)对于任意炉次i,计算其浇铸时间t<sub>i</sub>;<img file="FDA00001884616600011.GIF" wi="1879" he="124" />其中,Q<sub>i</sub>为炉次i的钢水重量,v<sub>o</sub>为奇流浇铸速度,w<sub>o</sub>为奇流结晶器宽度,h<sub>o</sub>为奇流结晶器厚度,v<sub>e</sub>为偶流浇铸速度,w<sub>e</sub>为偶流结晶器宽度,h<sub>e</sub>为偶流结晶器厚度,ρ为钢水密度取7.6吨/m<sup>3</sup>;1.2)对任意两个需要插铁板工艺的炉次i和j,计算插铁板工艺对钢水浇铸的影响时间τ<sub>ij</sub>:<img file="FDA00001884616600012.GIF" wi="1878" he="65" />其中,v<sub>i</sub>为炉次i的浇铸拉速,v<sub>j</sub>为炉次j的浇铸拉速,<img file="FDA00001884616600013.GIF" wi="313" he="59" />对任意两个不需要插铁板工艺的炉次i和j,τ<sub>ij</sub>=0;1.3)计算中间包更换过程对钢水浇铸的影响时间:<img file="FDA00001884616600014.GIF" wi="1878" he="64" />其中,v<sub>i</sub>为换前旧中间包内最后一炉钢水i的浇铸拉速,v<sub>j</sub>为换后新中间包内第一炉钢水j的浇铸拉速,<img file="FDA00001884616600015.GIF" wi="313" he="59" />1.4)对在浇次内的首炉i,计算开浇过程对钢水浇铸的影响时间ζ<sub>i</sub>;ζ<sub>i</sub>=(2v<sub>i</sub>-0.3-v<sub>0</sub>)/v<sub>i</sub>+(v<sub>i</sub>-v<sub>0</sub>)<sup>2</sup>/(2av<sub>i</sub>)=4v<sub>i</sub>-0.06/v<sub>i</sub>-1.2,        (4)开浇过程对钢水浇铸影响时间的计算依据如下:开浇升速过程为塞棒打开钢水注入,钢水上升到中间包液面时,以一定速度v<sub>0</sub>开浇并保持一定时间,再按一定的加速度a将速度上升到首炉标准拉速;1.5)对在浇次内的尾炉j,计算终浇过程对钢水浇铸的影响时间υ<sub>j</sub>:υ<sub>j</sub>=1-v<sub>0</sub>/v<sub>j</sub>+(v<sub>j</sub>-v<sub>0</sub>)<sup>2</sup>/(2av<sub>j</sub>)=4v<sub>j</sub>+0.24/v<sub>j</sub>-2.2                    (5);1.6)按式(6)建立任意两个炉次i和j之间的钢级连浇映射关系α<sub>ij</sub>;<img file="FDA00001884616600021.GIF" wi="1893" he="302" />1.7)按式(7)建立任意两个炉次i和j之间的铸坯宽度在线调节映射关系β<sub>ij</sub>;<img file="FDA00001884616600022.GIF" wi="1878" he="240" />其中B<sub>j</sub>和E<sub>i</sub>分别为j炉第一块和炉i最后一块铸坯的宽度,σ<sub>1</sub>为结晶器宽度在线正向最大允许调整幅度,σ<sub>2</sub>为结晶器宽度在线逆向最大允许调整幅度,k<sub>3</sub>为实际浇铸过程中i炉和j炉钢水在线调宽时铸出的梯形板坯导致的钢水收得率损失;1.8)按式(8)建立任意一个炉次i和任意一台连铸机c之间的浇铸映射关系d<sub>ic</sub>;<img file="FDA00001884616600023.GIF" wi="1894" he="141" />2)将连铸生产过程中钢水浇铸顺序方案映射为数学模型决策变量,其具体过程如下:2.1)设定炉次间浇铸顺序的0-1决策变量x<sub>ijbc</sub>,当决定炉次i和j都在铸机c的第b个浇次内,且炉次j紧接着炉次i之后浇铸时,x<sub>ijbc</sub>取1;否则x<sub>ijbc</sub>取0;2.2)设定浇次的开浇炉次0-1决策变量x<sub>0jbc</sub>,当决定炉次j为铸机c的第b个浇次的开浇炉次时,x<sub>0jbc</sub>取1;否则x<sub>0jbc</sub>取0;2.3)设定浇次的终浇炉次0-1决策变量x<sub>j0bc</sub>,当决定炉次j为铸机c的第b个浇次的终浇炉次时,x<sub>j0bc</sub>取1;否则x<sub>j0bc</sub>取0;2.4)设定中间包更换方案的0-1决策变量z<sub>ilbc</sub>,浇铸完炉次i后决定更换铸机c上第b个浇次内 的第l中间包,z<sub>ilbc</sub>取1;否则z<sub>ilbc</sub>取0;2.5)设定浇次更换方案0-1决策变量y<sub>bc</sub>,铸机c上第b个浇次内所有炉次浇铸完成后,需要更换浇次,y<sub>bc</sub>取1;否则,铸机c上第b个浇次为空浇次不需更换,y<sub>bc</sub>取0;2.6)设定附加变量s<sub>j</sub>,e<sub>j</sub>,分别表示炉次j的开浇时间和终浇时间,设定附加变量S<sub>bc</sub>,E<sub>bc</sub>,分别表示铸机c上第b个浇次的开浇时间和终浇时间;3)将连铸生产过程的工艺限制客观映射为数学模型约束条件,其具体过程如下3.1)建立连铸生产过程中控制钢水浇铸顺序约束条件:<img file="FDA00001884616600031.GIF" wi="1879" he="104" /><img file="FDA00001884616600032.GIF" wi="1904" he="103" /><img file="FDA00001884616600033.GIF" wi="1902" he="101" />其中,N为所有待生产的炉次集合,C为铸机集合,B<sub>c</sub>为某一铸机c上的预设浇次集合;3.2)建立连铸生产过程中更换浇次的约束条件:<img file="FDA00001884616600034.GIF" wi="1904" he="101" />其中,M为一个足够大的正数,y<sub>bc</sub>≥y<sub>b+1c</sub>,c∈C,b∈B<sub>c</sub>/{|B<sub>c</sub>|},                    (13);3.3)建立连铸生产过程中更换中间包的约束条件:<img file="FDA00001884616600035.GIF" wi="1904" he="100" /><img file="FDA00001884616600036.GIF" wi="1904" he="93" />其中,L<sub>b</sub>为浇次b内预设的中间包集合;3.4)建立连铸生产过程中与生产设备和工艺相关的约束条件:<img file="FDA00001884616600037.GIF" wi="1916" he="104" /><img file="FDA00001884616600038.GIF" wi="1904" he="104" /><img file="FDA00001884616600039.GIF" wi="1916" he="104" /><img file="FDA000018846166000310.GIF" wi="1903" he="104" /><img file="FDA000018846166000311.GIF" wi="1903" he="104" /><img file="FDA000018846166000312.GIF" wi="1903" he="104" /><img file="FDA00001884616600041.GIF" wi="1916" he="93" />其中,T<sub>tudish</sub>为中间包使用寿命;<img file="FDA00001884616600042.GIF" wi="1916" he="101" />其中,T<sub>cast</sub>为工艺允许的一个浇次最大时长;<img file="FDA00001884616600043.GIF" wi="1916" he="101" />其中,T<sub>setup</sub>为重设新浇次所需的机器清理时间;4)将连铸生产过程中优化的工艺指标映射为数学模型目标函数,其具体过程如下:4.1)将钢水收得率指标映射为式(25)的目标函数,即最小化异钢级连浇产生的交接板坯和结晶器在线调宽形成的梯形板坯所带来的钢水收得率损失:<img file="FDA00001884616600044.GIF" wi="1911" he="108" />4.2)将质量成本指标映射为式(26)的目标函数,即最小化异钢级连浇产生的交接板坯带来以优充次成本损失:<img file="FDA00001884616600045.GIF" wi="1911" he="108" />4.3)将中间包使用效率指标映射为式(27)的目标函数,即最小化中间包更换次数:<img file="FDA00001884616600046.GIF" wi="1911" he="104" />4.4)将铸机作业率指标映射为式(28)的目标函数,即最小化浇次更换次数:<img file="FDA00001884616600047.GIF" wi="1911" he="104" />4.5)将整个连铸系统生产效率指标映射为式(29)的目标函数,即最小化完成所有炉次钢水浇铸所需的时间:Min max<sub>c∈C,b∈B</sub>(E<sub>bc</sub>);                                        (29)所述步骤3以基本工艺约束为依据,制定钢水浇铸顺序的预控制方案的过程如下:步骤3.1:组炉批:对所有炉次N按照钢级、铸坯宽度和厚度进行分组,使每一组中的炉次具有相同宽度、厚度和钢级,对于每一组同宽度、厚度和钢级的炉次N<sub>b</sub>,将它们划分为<img file="FDA00001884616600048.GIF" wi="319" he="94" />个炉批;步骤3.2:炉批组浇:对所有的炉批,按照其对应钢级的钢种类别和铸坯厚度分组多个炉批组,对每组炉批,按照精炼类别,异钢种连浇规则及调宽规则进行合并,形成一个或多个浇次; 步骤3.3:浇次分配:对所有浇次,按照浇次中所包含炉次的属性与铸机属性之间的对应关系,并考虑铸机负荷平衡因素,将浇次分配到炉次;步骤3.4:浇次排程:对每个铸机上的分配的浇次,枚举所有可能的序列,选择其中一个使得精炼生产均衡化的序列为铸机上的浇次排程;步骤3.5:确定变量取值:由浇次排程结果获得所有与顺序相关的决策变量的取值,包括x<sub>ijbc</sub>,x<sub>0jbc</sub>,x<sub>j0bc</sub>,y<sub>bc</sub>的取值,计算各个炉次、浇次的开浇和终浇时间,即s<sub>j</sub>,e<sub>j</sub>和S<sub>bc</sub>,E<sub>bc</sub>,按照逻辑关系获得变量z<sub>ilbc</sub>的取值;所述步骤4采用禁忌搜索算法或变邻域搜索算法优化连铸生产过程各项工艺指标,改进预控制方案,设定钢水浇铸顺序控制方案;所述禁忌搜索算法,具体包括如下步骤:第1步:初始化算法参数,包括最大运行代数MaxIter、最大连续未改进代数UnimproveIter、禁忌表长度TabuListLength、禁忌表T=φ、历史最优目标值HisBestObj=∞、历史最优排序π<sup>HBST</sup>=φ、初始排序π=φ、当前迭代次数iter=0、连续未改进代数unIter=0、炉次信息;第2步:通过步骤三中所述的确定钢水浇铸顺序预控制方案的过程,获得所有决策变量x<sub>ijbc</sub>,x<sub>0jbc</sub>,x<sub>j0bc</sub>,y<sub>bc</sub>,z<sub>ilbc</sub>,s<sub>j</sub>,e<sub>j</sub>,S<sub>bc</sub>和E<sub>bc</sub>的取值,将变量取值,计算各个目标函数分量的值,并对目标函数分量值进行加权计算,获得总目标值f(π<sup>0</sup>),其中π<sup>0</sup>=(x<sub>ijbc</sub>,x<sub>0jbc</sub>,x<sub>j0bc</sub>,y<sub>bc</sub>,z<sub>ilbc</sub>,s<sub>j</sub>,e<sub>j</sub>,S<sub>bc</sub>,E<sub>bc</sub>),令HisBestObj=f(π<sup>0</sup>),π<sup>HBST</sup>=π<sup>0</sup>;第3步:根据邻域搜索策略求出邻域最优移动<img file="FDA00001884616600051.GIF" wi="116" he="60" />以及邻域最优排列<img file="FDA00001884616600052.GIF" wi="119" he="61" />如果<img file="FDA00001884616600053.GIF" wi="474" he="60" />则<img file="FDA00001884616600054.GIF" wi="275" he="61" /><img file="FDA00001884616600055.GIF" wi="502" he="60" />unIter=0,转到第5步;否则,unIter=unIter+1;第4步:如果unIter=UnimproveIter,则停止;第5步:更新禁忌表,将<img file="FDA00001884616600056.GIF" wi="90" he="60" />所表征的相关信息v<sup>lbest</sup>加入到禁忌表中,如果禁忌表未满,则依次添加;若禁忌表已满,禁忌表中最后一个元素解禁,其他的元素依次下移一个位置,然后将新的元素加入到第一个位置;第6步:iter=iter+1,如果iter=MaxIter,则停止;否则转到第3步。
地址 110819 辽宁省沈阳市和平区文化路3号巷11号