发明名称 DC轧机轧制压力、辊间压力的预报方法
摘要 本发明公开一种DC轧机轧制压力、辊间压力的预报方法。该方法包括以下步骤:(a)收集实际DC轧机的设备参数及工艺参数;(b)将辊系及轧件沿辊身方向离散化;(c)假设初始辊缝形状;(d)设定轧制压力;(e)设定辊间压力;(f)计算辊缝形状h1(y),并判断是否收敛,不收敛转到步骤(d);(g)验证工作辊的力(力矩)平衡条件。本发明根据DC轧机偏移交叉轧制特性,充分考虑滚动摩擦等因素的影响,合理设定DC轧机的轧制压力、辊间压力,耦合金属模型和辊系变形模型预报DC轧机的轧制压力、辊间压力。该发明是一种精度很高的预报方法,不仅有利于提高辊系变形模型计算的准确程度,而且有利于提高DC轧机板厚板形控制精度,从而提高产品的质量。
申请公布号 CN101716604B 申请公布日期 2012.06.13
申请号 CN200910227864.6 申请日期 2009.12.17
申请人 燕山大学 发明人 彭艳;刘宏民;陆小武
分类号 B21B37/00(2006.01)I;B21B37/48(2006.01)I 主分类号 B21B37/00(2006.01)I
代理机构 秦皇岛市维信专利事务所 13102 代理人 鄂长林
主权项 1.一种DC轧机轧制压力、辊间压力的预报方法,其特征是:所述方法包括以下步骤:(a)收集实际DC轧机的设备参数及工艺参数:包括支承辊的辊身长度L<sub>b</sub>、压下支点与辊身端部距离L<sub>bz</sub>、水平力作用点与辊身端部距离L<sub>bx</sub>、辊身直径D<sub>b</sub>或半径R<sub>b</sub>、辊颈直径D<sub>b1</sub>、弹性模量E<sub>b</sub>、泊松比v<sub>b</sub>,工作辊辊身长度L<sub>w</sub>、弯辊力作用点与辊身端部距离L<sub>wz</sub>、水平力作用点与辊身端部距离L<sub>wx</sub>、辊身直径D<sub>w</sub>或半径R<sub>w</sub>、辊颈直径D<sub>w1</sub>、弹性模量E<sub>w</sub>、泊松比v<sub>w</sub>、交叉角θ、偏移量e<sub>c</sub>、摩擦圆半径ρ<sub>w</sub>,工作辊和支承辊之间的辊间摩擦系数μ<sub>R</sub>、滚动摩擦力臂M,轧件来料规格尺寸B×h×l、入口厚度横向分布拟合曲线、弹性模量E<sub>s</sub>、泊松比v<sub>s</sub>,前后总张力T<sub>1</sub>、T<sub>0</sub>;(b)将辊系及轧件沿辊身方向离散化:在工作辊和支承辊辊身长度L范围内,沿辊身方向,将辊系划分为m个单元,将轧件划分为n个单元,其中m=n+2d,m、n均为奇数,d为辊身端部到同侧轧件边部的辊系单元划分个数,以左压下支点为坐标原点,单元宽度为Δy<sub>i</sub>,i=1,2,…,m,各单元的中点坐标为y<sub>i</sub>,i=1,2,…,m;将作用在轧辊上的轧制压力、辊间压力及辊系变形也按相同单元离散化;(c)假设初始辊缝形状:假设初始辊缝形状为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>h</mi><mn>1</mn><mo>*</mo></msubsup><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>b</mi><mn>0</mn></msub><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>y</mi></mrow><mi>B</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>b</mi><mn>4</mn></msub><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>y</mi></mrow><mi>B</mi></mfrac><mo>)</mo></mrow><mn>4</mn></msup></mrow></math>]]></maths>式中:b<sub>0</sub>-轧后厚度横向分布零次回归系数;b<sub>2</sub>-轧后厚度横向分布二次回归系数;b<sub>4</sub>-轧后厚度横向分布四次回归系数;(d)设定轧制压力:设定轧制压力方向α<sub>i</sub>为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mo>-</mo><mi>arcsin</mi><mfrac><mrow><mn>2</mn><msub><mi>D</mi><mi>w</mi></msub><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><mi>L</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mi>&theta;</mi></mrow><mrow><msubsup><mi>D</mi><mi>w</mi><mn>2</mn></msubsup><mo>+</mo><mn>2</mn><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><mi>L</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup><msup><mi>&theta;</mi><mn>2</mn></msup></mrow></mfrac><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mi>d</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>d</mi><mo>+</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>d</mi><mo>+</mo><mi>n</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>(e)设定辊间压力,包括以下有计算机系统执行的步骤:e1)设定辊间压力方向初始假设值β<sub>i0</sub>为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>=</mo><mi>arctan</mi><msub><mi>&mu;</mi><mi>R</mi></msub><mo>-</mo><mfrac><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><mi>L</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mi>&theta;</mi><msub><mrow><mo>+</mo><mi>e</mi></mrow><mi>c</mi></msub></mrow><mrow><msub><mi>R</mi><mi>b</mi></msub><mo>+</mo><msub><mi>R</mi><mi>w</mi></msub></mrow></mfrac><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>;</mo></mrow></math>]]></maths>e2)基于影响函数法,建立DC轧机辊系变形模型和金属模型,计算DC轧机辊系所受的未知力、各单元的工作辊水平方向挠度X<sub>wi</sub>和支承辊水平方向挠度X<sub>bi</sub>;e3)重新设定辊间压力方向β<sub>i</sub>为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>=</mo><mi>arctan</mi><msub><mi>&mu;</mi><mi>R</mi></msub><mo>-</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><msub><mi>X</mi><mi>wi</mi></msub><mo>-</mo><msub><mi>X</mi><mi>bi</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>D</mi><mi>w</mi></msub><mo>+</mo><msub><mi>D</mi><mi>b</mi></msub></mrow></mfrac><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>;</mo></mrow></math>]]></maths>e4)以β<sub>i</sub>、β<sub>i0</sub>的均方差是否在0~1.0×10<sup>-5</sup>度范围之内为收敛判据,比较β<sub>i</sub>、β<sub>i0</sub>,判断是否收敛,若收敛进入步骤(f),否则,用松弛因子法修改β<sub>i0</sub>,转入步骤e2);(f)计算辊缝形状h<sub>1</sub>(y),并以出口厚度横向分布变化量的最大值是否在0~1.0×10<sup>-5</sup>mm范围之内为收敛判据,比较h<sub>1</sub>(y)、<img file="FSB00000594886200023.GIF" wi="140" he="59" />判断是否收敛,若收敛,输出各单元单位宽度轧制压力p<sub>i</sub>、辊间压力q<sub>i</sub>和工作辊水平力F<sub>wxl</sub>、F<sub>wxr</sub>,进入步骤(g),否则,用松弛因子法修改<img file="FSB00000594886200024.GIF" wi="140" he="58" />转入步骤(d);(g)根据步骤(f)输出结果,验证工作辊的力或力矩平衡条件,包括以下由计算机系统执行的步骤:g1)计算各单元轧制压力对工作辊轴心的力臂a<sub>i</sub>为:a<sub>i</sub>=R<sub>w</sub> sin(α<sub>i</sub>+φ<sub>i</sub>)         i=1,2,…,m式中:R<sub>w</sub>-工作辊半径;φ<sub>i</sub>-轧制压力作用点所对应的轧辊中心角,设定为咬入角的一半;g2)计算各单元辊间压力对工作辊轴心的力臂b<sub>i</sub>为:b<sub>i</sub>=R<sub>w</sub>sin(β<sub>i</sub>+γ<sub>i</sub>)-Mcosβ<sub>i</sub> i=1,2,…,m;g3)受力模型中已经满足工作辊水平和铅垂方向的力平衡条件,只需验证工作辊在轧制压力、辊间压力、前后张力和水平支承力作用下的力矩平衡条件,验证条件为:<maths num="0005"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>d</mi><mo>+</mo><mn>1</mn></mrow><mrow><mi>d</mi><mo>+</mo><mi>n</mi></mrow></munderover><msub><mi>a</mi><mi>i</mi></msub><msub><mi>p</mi><mi>i</mi></msub><mi>&Delta;</mi><msub><mi>y</mi><mi>i</mi></msub><msub><mi>R</mi><mi>w</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>F</mi><mi>wxl</mi></msub><mo>+</mo><msub><mi>F</mi><mi>wxr</mi></msub><mo>)</mo></mrow><msub><mi>&rho;</mi><mi>w</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>b</mi><mi>i</mi></msub><msub><mi>q</mi><mi>i</mi></msub><mi>&Delta;</mi><msub><mi>y</mi><mi>i</mi></msub><mo>+</mo><msub><mi>R</mi><mi>w</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><mn>2</mn><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi></mrow></math>]]></maths>式中:F<sub>wxl</sub>-工作辊左水平支承力;F<sub>wxr</sub>-工作辊右水平支承力;辊系受力是否满足工作辊平衡条件,可以用误差ε来判断;误差ε的具体含义是:上式两边之差的绝对值与两边最小值之比,则误差ε可表示为:<maths num="0006"><![CDATA[<math><mrow><mi>&epsiv;</mi><mo>=</mo><mo>|</mo><mfrac><mrow><mi>A</mi><mo>-</mo><mi>B</mi></mrow><mrow><mi>MIN</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>)</mo></mrow></mrow></mfrac><mo>|</mo></mrow></math>]]></maths>其中:<maths num="0007"><![CDATA[<math><mrow><mi>A</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>d</mi><mo>+</mo><mn>1</mn></mrow><mrow><mi>d</mi><mo>+</mo><mi>n</mi></mrow></munderover><msub><mi>a</mi><mi>i</mi></msub><msub><mi>p</mi><mi>i</mi></msub><mi>&Delta;</mi><msub><mi>y</mi><mi>i</mi></msub><msub><mi>R</mi><mi>w</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>F</mi><mi>wxl</mi></msub><mo>+</mo><msub><mi>F</mi><mi>wxr</mi></msub><mo>)</mo></mrow><msub><mi>&rho;</mi><mi>w</mi></msub><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mi>B</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>b</mi><mi>i</mi></msub><msub><mi>q</mi><mi>i</mi></msub><mi>&Delta;</mi><msub><mi>y</mi><mi>i</mi></msub><mo>+</mo><msub><mi>R</mi><mi>w</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><mn>2</mn><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi></mrow></math>]]></maths>MIN(A,B)表示取A和B中的最小值;若误差ε满足<img file="FSB00000594886200034.GIF" wi="485" he="140" />则可认为辊系受力满足工作辊平衡条件,说明辊系受力模型准确,轧制压力、辊间压力预报精度高。
地址 066004 河北省秦皇岛市海港区河北大街438号