发明名称 一种冷轧连续退火机组钢卷优化排序方法及其系统
摘要 一种冷轧连续退火机组钢卷优化排序方法及其系统,属于金属材料加工信息技术领域,优化方法包括以下步骤:一、将候选钢卷按退火温度分别从高到低和从低到高排序,形成两种初始排序方案,对每种初始排序方案采取宽度优先排序或厚度优先排序方法进行优化,获得多组初始可行排序方案;二、从初始钢卷排序方案中选择优化目标值最小的排序方案作为初始可行生产计划;三、对初始可行生产计划,以优化排序模型目标值最小为目标,用交换邻域禁忌搜索和交替路径变换邻域搜索进行调整。基于此方法提出相应系统,本发明使得生产计划执行过程中切换减少,过渡平滑,提高了产品质量,挖掘了产能。
申请公布号 CN101334660B 申请公布日期 2011.01.26
申请号 CN200810012090.0 申请日期 2008.06.30
申请人 东北大学 发明人 唐立新;杨阳;庞博
分类号 G05B19/418(2006.01)I 主分类号 G05B19/418(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 梁焱
主权项 1.一种冷轧连续退火机组钢卷排序方法,其特征在于:基于加工工艺制定连续退火机组初始排序方案,将连续退火机组工艺约束描述为以最小化切换成本和波动附加成本为目标的数学模型,依据此目标对排序方案进行改进,减少退火温度、宽度和厚度的切换和波动,在连续退火机组钢卷排序方法的基础上,提供了一套连续退火计划排产系统,该系统衔接钢铁企业的内部信息管理平台和生产现场自动控制系统,依据信息管理平台提供的信息通过连续退火优化调度系统中的自动钢卷优化排序模块进行排序并由现场自动控制系统予以实施,该排序方法包括以下步骤:步骤一:建立一系列连续退火机组初始可行钢卷排序方案;将候选钢卷按照退火温度分别从高到低和从低到高进行排序,形成两种初始钢卷排序方案,对每种初始钢卷排序方案采取宽度优先排序或者厚度优先排序进行预排;所述宽度优先排序是指在每种初始钢卷排序方案中,对于具有相同退火温度的钢卷再根据钢卷的宽度,分别采用从宽到窄和从窄到宽两种不同的排序规则排序,获得四种钢卷排序方案;再对具有相同温度、相同宽度的钢卷根据钢卷的厚度从厚到薄和从薄到厚两种不同的排序规则排序,获得八个初始可行钢卷排序方案;所述厚度优先排序是指在每种初始钢卷排序方案中,对于具有相同退火温度的钢卷再根据钢卷的厚度,分别采用从厚到薄和从薄到厚两种不同的排序规则排序,获得四种钢卷排序方案;再对具有相同温度、相同厚度的钢卷根据钢卷的宽度从宽到窄和从窄到宽两种不同的排序规则排序,获得八个初始可行钢卷排序方案;步骤二:比较连续退火机组初始可行钢卷排序方案;从已获得的一系列钢卷排序方案中,选择优化目标值最小的排序方案作为初始可行钢卷排序方案;所述优化目标值最小,指切换成本和波动附加成本目标值最小,所述的波动附加成本包括退火温度、宽度、厚度变化趋势发生波动所造成的波动附加成本,所述的切换成本可根据如下公式进行计算,相邻钢卷切换成本为<img file="FSB00000311127200011.GIF" wi="271" he="97" />其中:<img file="FSB00000311127200012.GIF" wi="1043" he="173" /><img file="FSB00000311127200013.GIF" wi="46" he="67" />为相邻钢卷退火温度切换成本,v<sub>e</sub>为炉温变化速度,ρ<sub>steel</sub>为钢的密度,T<sub>i</sub>表示钢卷i要求的退火温度,炉温从T<sub>i</sub>切换到T<sub>j</sub>的过程中,利益损失为r元,δ<sup>T</sup>为允许最大的退火温度跳跃,λ<sup>T</sup>为相邻钢卷退火温度切换成本权重;<img file="FSB00000311127200021.GIF" wi="1218" he="225" /><img file="FSB00000311127200022.GIF" wi="46" he="59" />为相邻钢卷规格跳跃造成的切换成本,s<sub>i</sub>为钢卷i的横截面积,即s<sub>i</sub>=w<sub>i</sub>×t<sub>i</sub>,δ<sup>s</sup>为最大横截面积切换比例;<img file="FSB00000311127200023.GIF" wi="963" he="243" /><img file="FSB00000311127200024.GIF" wi="48" he="60" />为相邻钢卷宽度切换成本,w<sub>i</sub>表示钢卷i要求的宽度,δ<sup>w</sup>为允许最大的宽度跳跃,λ<sup>w</sup>为相邻钢卷宽度切换成本权重;<img file="FSB00000311127200025.GIF" wi="837" he="276" /><img file="FSB00000311127200026.GIF" wi="42" he="63" />是相邻钢卷厚度切换成本,t<sub>i</sub>表示钢卷i要求的厚度,δ<sup>t</sup>为允许最大的厚度跳跃,λ<sup>t</sup>为相邻钢卷厚度切换成本权重;步骤三:确定满足冷轧连续退火工艺约束的钢卷排序方案并下发生产机组;以优化排序模型中给定的优化目标值最小为目的,满足优化排序模型约束条件,利用交换邻域禁忌搜索和交替路径变换邻域搜索两种策略调整初始可行钢卷排序方案,确定满足冷轧连续退火工艺约束的钢卷排序方案;所述优化排序模型可描述为如下数学模型:<maths num="0001"><![CDATA[<math><mrow><mi>min</mi><mi>imize</mi><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msubsup><mi>c</mi><mi>ij</mi><mi>T</mi></msubsup><mo>+</mo><msubsup><mi>c</mi><mi>ij</mi><mi>g</mi></msubsup><mo>)</mo></mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msubsup><mi>p</mi><mi>i</mi><mi>T</mi></msubsup><msubsup><mi>y</mi><mi>i</mi><mi>T</mi></msubsup><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mi>w</mi></msubsup><msubsup><mi>y</mi><mi>i</mi><mi>w</mi></msubsup><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mi>t</mi></msubsup><msubsup><mi>y</mi><mi>i</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>约束条件:<maths num="0002"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>ij</mi></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>ij</mi></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>ik</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>kj</mi></msub><mo>=</mo><mn>0</mn><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><mi>S</mi></mrow></munder><msub><mi>x</mi><mi>ij</mi></msub><mo>&le;</mo><mo>|</mo><mi>S</mi><mo>|</mo><mo>-</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>S</mi><mo>&SubsetEqual;</mo><mo>{</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>T</mi><mi>i</mi></msub><mo>-</mo><msub><mi>T</mi><mi>j</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>]</mo><mo>&times;</mo><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>T</mi><mi>j</mi></msub><mo>-</mo><msub><mi>T</mi><mi>k</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>jk</mi></msub><mo>]</mo><mo>&times;</mo><msubsup><mi>y</mi><mi>i</mi><mi>T</mi></msubsup><mo>-</mo><mi>M</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>y</mi><mi>i</mi><mi>T</mi></msubsup><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mo>,</mo><mo>&ForAll;</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>-</mo><msub><mi>w</mi><mi>j</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>]</mo><mo>&times;</mo><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>w</mi><mi>j</mi></msub><mo>-</mo><msub><mi>w</mi><mi>k</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>jk</mi></msub><mo>]</mo><mo>&times;</mo><msubsup><mi>y</mi><mi>i</mi><mi>w</mi></msubsup><mo>-</mo><mi>M</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>y</mi><mi>i</mi><mi>w</mi></msubsup><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mo>,</mo><mo>&ForAll;</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mi>j</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>]</mo><mo>&times;</mo><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>t</mi><mi>j</mi></msub><mo>-</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>jk</mi></msub><mo>]</mo><mo>&times;</mo><msubsup><mi>y</mi><mi>i</mi><mi>t</mi></msubsup><mo>-</mo><mi>M</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>y</mi><mi>i</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mo>,</mo><mo>&ForAll;</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>,</mo><msubsup><mi>y</mi><mi>i</mi><mi>T</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mi>i</mi><mi>w</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mi>i</mi><mi>t</mi></msubsup><mo>&Element;</mo><mo>{</mo><mn>0,1</mn><mo>}</mo><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,x<sub>ij</sub>表示连续退火钢卷排序方案中,当钢卷i在钢卷j紧前进行退火,x<sub>ij</sub>取值为1,否则x<sub>ij</sub>取值为0;<img file="FSB00000311127200036.GIF" wi="54" he="60" />表示连续退火钢卷排序方案中,退火温度的变化趋势在加工钢卷i处是否发生调整,如果发生调整,<img file="FSB00000311127200037.GIF" wi="53" he="61" />取值为1,否则<img file="FSB00000311127200038.GIF" wi="53" he="61" />取值为0;<img file="FSB00000311127200039.GIF" wi="41" he="56" />表示连续退火钢卷排序方案中,厚度的变化趋势在钢卷i处是否发生调整,如果发生调整,<img file="FSB000003111272000310.GIF" wi="41" he="55" />取值为1,否则<img file="FSB000003111272000311.GIF" wi="40" he="55" />取值为0;<img file="FSB000003111272000312.GIF" wi="47" he="47" />表示连续退火钢卷排序方案中,宽度的变化趋势在钢卷i处是否发生调整,如果发生调整,<img file="FSB000003111272000313.GIF" wi="56" he="54" />取值为1,否则<img file="FSB000003111272000314.GIF" wi="48" he="45" />取值为0,<img file="FSB000003111272000315.GIF" wi="45" he="67" />为由于相邻钢卷i和j的退火温度差异引起的设备切换成本,<img file="FSB000003111272000316.GIF" wi="46" he="60" />为由于相邻钢卷i和j的规格引起的设备切换成本,<img file="FSB000003111272000317.GIF" wi="50" he="50" />是由于在钢卷排序方案中钢卷i处退火温度的变化趋势发生调整引起的设备调整成本,<img file="FSB000003111272000318.GIF" wi="42" he="47" />是由于在钢卷排序方案中钢卷i处厚度的变化趋势发生调整引起的设备调整成本,<img file="FSB000003111272000319.GIF" wi="52" he="47" />是由于在钢卷排序方案中钢卷i处宽度的变化趋势发生调整引起的设备调整成本,T<sub>i</sub>表示钢卷i要求的退火温度,w<sub>i</sub>表示钢卷i的宽度,t<sub>i</sub>表示钢卷i的厚度;λ<sub>1</sub>和λ<sub>2</sub>分别是切换成本与波动附加成本在目标中所占的权重;所述的交换邻域禁忌搜索策略,指交换初始可行钢卷排序方案中任意两个钢卷的退火位置,寻找优化排序模型目标值最小的交换改进方案及相应的钢卷排序方案;所述的交替路径变换邻域搜索,指用如下动态规划迭代式进行,设(S,L)表示辅助改进图中从I<sub>0</sub>到I<sub>L</sub>的改进路径,其中I<sub>S</sub>是在I<sub>L</sub>紧前被访问;f(S,L)是改进路径(S,L)对应的模型中优化目标的最好改进量;辅助改进图中的每条弧(I<sub>j</sub>,I<sub>k</sub>)对应的改进移动所引起的对优化目标值中切换成本的改进量定义为D[j,k],可按照如下公式进行计算:当1≤j≤n时,D[0,j]=-c[I<sub>n</sub>,I<sub>0</sub>]+c[I<sub>0</sub>,I<sub>j</sub>],当1≤j<k≤n时,D[j,]=-c[I<sub>j-1</sub>,I<sub>J</sub>]+c[I<sub>j-1</sub>,I<sub>k</sub>].P[j]表示在辅助改进图的改进路径末端插入节点I<sub>j</sub>所对应的对初始可行钢卷排序方案的改变所带来的对优化目标中变化趋势波动成本的改进量,其动态规划迭代式给定如下:<img file="FSB00000311127200041.GIF" wi="1363" he="176" />此迭代式可被初始化为f(0,L)=D[0,L].式中,f(S′,s)为从I<sub>0</sub>到I<sub>S</sub>并且I<sub>S′</sub>在I<sub>S</sub>紧前被访问的改进路径对应的最好改进量,即优化排序模型目标值最小,其中S′<S<L,如上迭代式试图找到一条从I<sub>0</sub>到I<sub>L</sub>并且I<sub>S</sub>在I<sub>L</sub>紧前被访问的,对初始可行钢卷排序方案目标值改进最好的改进路径,当完成寻找一条从I<sub>0</sub>到I<sub>n</sub>并且I<sub>S</sub>在I<sub>n</sub>紧前被访问的具有最好改进量的改进路径(S,n)的任务时,迭代过程完毕,搜索交替路径变换邻域所得的最优钢卷排序方案可通过求解<img file="FSB00000311127200042.GIF" wi="498" he="105" />获得。
地址 110004 辽宁省沈阳市和平区文化路3号巷11号