发明名称 一种热连轧机带材换热系数自学方法
摘要 一种热连轧机带材换热系数自学方法,它主要包括以下由计算机执行的步骤:1、收集轧机参数、轧制工艺参数以及换热系数当前值;2、对除鳞区对流换热系数、机架间水冷对流换热系数、带材两侧边水冷对流换热系数、空冷换热系数、轧件与轧辊之间接触热传导换热系数进行寻优。本发明通过对五个换热系数在较大区间进行整体寻优,再固定各个换热系数对其它换热系数在较小区间进行寻优的方法,维持了换热系数之间相对的定性关系,保证了自学之后的换热系数在不失去物理意义的前提下更加准确的预报带材温度。
申请公布号 CN105149361B 申请公布日期 2017.04.05
申请号 CN201510398450.5 申请日期 2015.07.08
申请人 燕山大学 发明人 王东城;张威;吴燕林;刘宏民
分类号 B21B37/00(2006.01)I 主分类号 B21B37/00(2006.01)I
代理机构 秦皇岛一诚知识产权事务所(普通合伙) 13116 代理人 崔凤英
主权项 一种热连轧机带材换热系数自学习方法,其特征在于:它主要包括以下由计算机执行的步骤:a、收集轧机参数、轧制工艺参数以及换热系数当前值,具体包括:a1、收集轧机参数包括:热连轧机机架数n、精轧入口测温仪距离第一机架间距L<sub>M01</sub>、机架间距L<sub>M12</sub>~L<sub>M(n‑1)n</sub>、末机架距离精轧出口测温仪间距L<sub>Mn(n+1)</sub>、除鳞箱中心至精轧入口测温仪间距L<sub>T0</sub>、各机架后水冷节点至该机架中心间距L<sub>T1</sub>~L<sub>Tn</sub>、除鳞箱喷洒覆盖区间长度L<sub>C0</sub>、各机架后水冷节点喷洒覆盖区间长度L<sub>C1</sub>~L<sub>Cn</sub>、各机架后水冷节点是否喷水的状态变量Bu<sub>1</sub>~Bu<sub>n</sub>,取值0或1,取0表示不存在或者不喷水,取1表示正常喷水,各机架工作辊半径D<sub>1</sub>~D<sub>n</sub>;a2、收集轧制工艺参数包括:当前卷带钢各机架轧制速度v<sub>1</sub>~v<sub>n</sub>、当前卷带钢来料宽度B、当前卷带钢来料平均厚度与各机架出口平均厚度<img file="FDA0001182584660000011.GIF" wi="163" he="70" />当前卷带钢各机架摩擦系数μ<sub>1</sub>~μ<sub>n</sub>、当前卷带钢各机架轧件轧制压力P<sub>1</sub>~P<sub>n</sub>、当前卷带钢精轧入口测温仪实测温度T<sub>in</sub>、当前卷带钢精轧出口测温仪实测温度T<sub>out</sub>、冷却水温度T<sub>c</sub>、空气温度T<sub>a</sub>;a3、收集换热系数当前值包括:除鳞区对流换热系数<img file="FDA0001182584660000012.GIF" wi="115" he="55" />机架间水冷对流换热系数<img file="FDA0001182584660000013.GIF" wi="110" he="55" />带材两侧边水冷对流换热系数<img file="FDA0001182584660000014.GIF" wi="111" he="63" />空冷换热系数<img file="FDA0001182584660000015.GIF" wi="115" he="59" />轧件与轧辊之间接触热传导换热系数<img file="FDA0001182584660000016.GIF" wi="113" he="55" />b、采用行向量<img file="FDA0001182584660000017.GIF" wi="44" he="55" />表示轧机参数,即:<maths num="0001"><math><![CDATA[<mrow><mover><mi>U</mi><mo>&OverBar;</mo></mover><mo>=</mo><mrow><mo>(</mo><mi>n</mi><mo>,</mo><msub><mi>L</mi><mrow><mi>M</mi><mn>01</mn></mrow></msub><mo>,</mo><msub><mi>L</mi><mrow><mi>M</mi><mn>12</mn></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>L</mi><mrow><mi>M</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>L</mi><mrow><mi>M</mi><mi>n</mi><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>,</mo><msub><mi>L</mi><mrow><mi>T</mi><mn>0</mn></mrow></msub><mo>,</mo><msub><mi>L</mi><mrow><mi>T</mi><mn>1</mn></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>L</mi><mrow><mi>T</mi><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>L</mi><mrow><mi>C</mi><mn>0</mn></mrow></msub><mo>,</mo><msub><mi>L</mi><mrow><mi>C</mi><mn>1</mn></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>L</mi><mrow><mi>C</mi><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>Bu</mi><mn>1</mn></msub><mo>...</mo><mo>,</mo><msub><mi>Bu</mi><mi>n</mi></msub><mo>,</mo><msub><mi>D</mi><mn>1</mn></msub><mo>...</mo><mo>,</mo><msub><mi>D</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001182584660000018.GIF" wi="1854" he="94" /></maths>采用行向量<img file="FDA0001182584660000019.GIF" wi="39" he="55" />表示轧制工艺参数,即:<maths num="0002"><math><![CDATA[<mrow><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mo>=</mo><mrow><mo>(</mo><msub><mi>v</mi><mn>1</mn></msub><mo>...</mo><mo>,</mo><msub><mi>v</mi><mi>n</mi></msub><mo>,</mo><mi>B</mi><mo>,</mo><mover><msub><mi>h</mi><mn>0</mn></msub><mo>&OverBar;</mo></mover><mo>,</mo><mo>...</mo><mo>,</mo><mover><msub><mi>h</mi><mi>n</mi></msub><mo>&OverBar;</mo></mover><mo>,</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>&mu;</mi><mi>n</mi></msub><mo>,</mo><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>T</mi><mrow><mi>i</mi><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>T</mi><mi>c</mi></msub><mo>,</mo><msub><mi>T</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00011825846600000110.GIF" wi="1110" he="86" /></maths> 采用有限差分法计算带材温度场,表示为:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>T</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow><mrow><mi>c</mi><mi>a</mi><mi>l</mi></mrow></msubsup><mo>=</mo><msubsup><mi>T</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow><mrow><mi>c</mi><mi>a</mi><mi>l</mi></mrow></msubsup><mrow><mo>(</mo><mover><mi>U</mi><mo>&RightArrow;</mo></mover><mo>,</mo><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mo>,</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>,</mo><msub><mi>&alpha;</mi><mn>3</mn></msub><mo>,</mo><msub><mi>&alpha;</mi><mn>4</mn></msub><mo>,</mo><msub><mi>&alpha;</mi><mn>5</mn></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA00011825846600000111.GIF" wi="701" he="93" /></maths>式中:α<sub>1</sub>,α<sub>2</sub>,α<sub>3</sub>,α<sub>4</sub>,α<sub>5</sub>为给定的除鳞区对流换热系数、机架间水冷对流换热系数、带材两侧边水冷对流换热系数、空冷换热系数、轧件与轧辊之间接触热传导换热系数;根据收集的参数对换热系数进行自学习,具体包括:b1、对换热系数进行大区间整体寻优,具体包括:b1(1)、令mine=1.0e10,i=0,<img file="FDA00011825846600000112.GIF" wi="1181" he="63" />b1(2)、令<img file="FDA00011825846600000113.GIF" wi="1114" he="71" /><img file="DA00011825846628365.GIF" wi="556" he="70" /><img file="FDA0001182584660000021.GIF" wi="1169" he="70" />采用有限差分法计算带材出口温度<img file="FDA0001182584660000022.GIF" wi="667" he="91" />计算精轧出口测温仪实测温度与计算温度之差<img file="FDA0001182584660000023.GIF" wi="307" he="86" />b1(3)、判断e&lt;mine是否成立,如果e&lt;mine成立,令mine=e,记录当前最佳的换热系数值<img file="FDA0001182584660000024.GIF" wi="1074" he="55" />转入b1(4),否则直接转入b1(4);b1(4)、判断i=20是否成立,如果i=20成立,转入b2,否则令i=i+1,转入b1(2);b2、固定除鳞区对流换热系数,对其它换热系数进行小区间寻优,具体包括:b2(1)、令<img file="FDA0001182584660000025.GIF" wi="179" he="55" />mine=1.0e10,i=0;b2(2)、令<img file="FDA0001182584660000026.GIF" wi="1119" he="71" /><img file="DA00011825846628390.GIF" wi="557" he="71" /><img file="FDA0001182584660000027.GIF" wi="553" he="71" />采用有限差分法计算带材出口温度<img file="FDA0001182584660000028.GIF" wi="670" he="87" />计算精轧出口测温仪实测温度与计算温度之差<img file="FDA0001182584660000029.GIF" wi="302" he="89" />b2(3)、判断e&lt;mine是否成立,如果e&lt;mine成立,令mine=e,记录当前最佳的换热系数值<img file="FDA00011825846600000210.GIF" wi="1070" he="55" />转入b2(4),否则直接转入b2(4);b2(4)、判断i=20是否成立,如果i=20成立,转入b3,否则令i=i+1,转入b2(2);b3、固定机架间水冷对流换热系数,对其它换热系数进行小区间寻优,具体包括:b3(1)、令<img file="FDA00011825846600000211.GIF" wi="187" he="55" />mine=1.0e10,i=0;b3(2)、令<img file="FDA00011825846600000212.GIF" wi="1107" he="67" /><img file="DA00011825846628420.GIF" wi="556" he="67" /><img file="FDA00011825846600000213.GIF" wi="556" he="67" />采用有限差分法计算带材出口温度<img file="FDA00011825846600000214.GIF" wi="675" he="91" />计算精轧出口测温仪实测温度与计算温度之差<img file="FDA00011825846600000215.GIF" wi="302" he="83" />b3(3)、判断e&lt;mine是否成立,如果e&lt;mine成立,令mine=e,记录当前最佳的换热系数值<img file="FDA00011825846600000216.GIF" wi="1072" he="63" />转入b3(4),否则直接转入b3(4);b3(4)、判断i=20是否成立,如果i=20成立,转入b4,否则令i=i+1,转入b3(2);b4、固定带材两侧边水冷对流换热系数,对其它换热系数进行小区间寻优,具体包括:b4(1)、令<img file="FDA0001182584660000031.GIF" wi="182" he="55" />mine=1.0e10,i=0;b4(2)、令<img file="FDA0001182584660000032.GIF" wi="1103" he="71" /><img file="DA00011825846628445.GIF" wi="572" he="70" /><img file="FDA0001182584660000033.GIF" wi="556" he="70" />采用有限差分法计算带材出口温度<img file="FDA0001182584660000034.GIF" wi="675" he="91" />计算精轧出口测温仪实测温度与计算温度之差<img file="FDA0001182584660000035.GIF" wi="302" he="85" />b4(3)、判断e&lt;mine是否成立,如果e&lt;mine成立,令mine=e,记录当前最佳的换热系数值<img file="FDA0001182584660000036.GIF" wi="1075" he="55" />转入b4(4),否则直接转入b4(4);b4(4)、判断i=20是否成立,如果i=20成立,转入b5,否则令i=i+1,转入b4(2);b5、固定空冷换热系数,对其它换热系数进行小区间寻优,具体包括:b5(1)、令<img file="FDA0001182584660000037.GIF" wi="162" he="50" />mine=1.0e10,i=0;b5(2)、令<img file="FDA0001182584660000038.GIF" wi="1107" he="69" /><img file="DA00011825846628473.GIF" wi="553" he="70" /><img file="FDA0001182584660000039.GIF" wi="554" he="70" />采用有限差分法计算带材出口温度<img file="FDA00011825846600000310.GIF" wi="669" he="93" />计算精轧出口测温仪实测温度与计算温度之差<img file="FDA00011825846600000311.GIF" wi="302" he="86" />b5(3)、判断e&lt;mine是否成立,如果e&lt;mine成立,令mine=e,记录当前最佳的换热系数值<img file="FDA00011825846600000312.GIF" wi="1073" he="63" />转入b5(4),否则直接转入b5(4);b5(4)、判断i=20是否成立,如果i=20成立,转入b6,否则令i=i+1,转入b5(2);b6、固定轧件与轧辊之间接触热传导换热系数,对其它换热系数进行小区间寻优,具体包括:b6(1)、令<img file="FDA00011825846600000317.GIF" wi="181" he="62" />mine=1.0e10,i=0;b6(2)、令<img file="FDA00011825846600000315.GIF" wi="1109" he="71" /><img file="DA00011825846628497.GIF" wi="554" he="65" /><img file="FDA00011825846600000316.GIF" wi="562" he="71" />采用有限差分法计算带材出口温度<img file="FDA0001182584660000041.GIF" wi="670" he="91" />计算精轧出口测温仪实测温度与计算温度之差<img file="FDA0001182584660000042.GIF" wi="302" he="91" />b6(3)、判断e&lt;mine是否成立,如果e&lt;mine成立,令mine=e,记录当前最佳的换热系数值<img file="FDA0001182584660000043.GIF" wi="1073" he="63" />转入b6(4),否则直接转入b6(4);b6(4)、判断i=20是否成立,如果i=20成立,转入b7,否则令i=i+1,转入b6(2);b7、通过指数平滑法对换热系数进行更新,令<img file="FDA0001182584660000044.GIF" wi="526" he="70" /><img file="DA00011825846628521.GIF" wi="519" he="71" /><img file="FDA0001182584660000045.GIF" wi="1038" he="71" /><img file="FDA0001182584660000046.GIF" wi="518" he="71" />其中β为指数平滑因子,自学习结束。
地址 066004 河北省秦皇岛市海港区河北大街西段438号