主权项 |
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>α</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>θ</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>θ</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>β</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>=</mo><mi>arctan</mi><msub><mi>μ</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>θ</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>β</mi><mi>i</mi></msub><mo>=</mo><mi>arctan</mi><msub><mi>μ</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>Σ</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>Δ</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>ρ</mi><mi>w</mi></msub><mo>=</mo><munderover><mi>Σ</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>Δ</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>ϵ</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>Σ</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>Δ</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>ρ</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>Σ</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>Δ</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" />则可认为辊系受力满足工作辊平衡条件,说明辊系受力模型准确,轧制压力、辊间压力预报精度高。 |