发明名称 一种热轧带钢开卷张力综合优化计算方法
摘要 一种热轧带钢开卷张力综合优化计算方法,它主要包括以下可由计算机执行的步骤:1、收集开卷机的设备工艺参数;2、收集带钢的品种规格及工艺参数;3、收集过程参数;4、给定张力的优化步长和优化过程;5、依次求开卷第k层时钢卷内任意第i层第s条的半径r<sub>si</sub>、第s条第i层与第i-1层间钢卷内部径向应力p<sub>si</sub>、周向应力q<sub>si</sub>、层间摩擦力τ<sub>si</sub>,开卷松卷及层间滑移的控制目标函数<img file="DDA0000570635770000011.GIF" wi="70" he="52" />6、输出最终确定的最优张力值。本发明能够预报张力综合优化后的开卷周向应力、开卷径向应力、开卷层间摩擦应力分布情况,有效防止带钢在开卷过程中出现松卷、层间滑移划伤带钢现象的发生。
申请公布号 CN104275368B 申请公布日期 2015.12.09
申请号 CN201410467747.8 申请日期 2014.09.15
申请人 燕山大学 发明人 白振华;王瑞;朱厚鸿;陈双玉;陈麒麟;陈博文
分类号 B21C51/00(2006.01)I;B21C47/18(2006.01)I 主分类号 B21C51/00(2006.01)I
代理机构 石家庄一诚知识产权事务所 13116 代理人 李合印
主权项 一种热轧带钢开卷张力综合优化计算方法,其特征在于:它主要包括以下由计算机执行的步骤:(a)收集开卷机的设备工艺参数,主要包括:开卷最大线速度v<sub>max</sub>,最大开卷张力T<sub>max</sub>,最小开卷张力T<sub>min</sub>,卷筒内半径r<sub>a</sub>,卷筒外半径r<sub>b</sub>,卷筒最大涨径r<sub>max</sub>,卷筒最小涨径r<sub>min</sub>,卷筒宽度B<sub>jt</sub>,允许钢卷最大重量G<sub>max</sub>,卷筒的弹性模量E<sub>1</sub>,卷筒的泊松比ν<sub>1</sub>;(b)收集带钢的品种规格范围及工艺参数,主要包括:带钢钢种,带钢平均厚度h,带钢厚度分布h<sub>s</sub>,带钢宽度B,钢卷卷重G,带钢表面粗糙度Ra,带钢摩擦系数μ,带钢的泊松比v<sub>2</sub>,带钢的弹性模量E<sub>2</sub>,带钢的密度ρ,钢卷外半径和内半径R<sub>1</sub>、R<sub>2</sub>,带材的周向弹性模量E,带材的径向弹性模量E<sub>r</sub>,带钢的板凸度Δh,来料钢卷的温度t,来料钢卷的层数C;(c)收集过程参数,主要包括:带材表面不平度的最大值δ<sub>max</sub>,考虑理论与实际存在差别的修正系数k<sub>0</sub>,加权系数α,带钢横向条元数m,钢卷层数n;(d)给定张力的优化步长ΔT,张力优化过程参数j<sub>1</sub>、j<sub>2</sub>、j<sub>1y</sub>,开卷松卷及层间滑移的控制目标函数初始值<img file="FDA0000570635740000011.GIF" wi="195" he="73" />计算过程的收敛极限ε;(e)令j<sub>1</sub>=0,开卷张力T=T<sub>min</sub>+j<sub>1</sub>·ΔT;(f)定义开卷时钢卷最外层层数为k;(g)令j<sub>2</sub>=0;(h)令k=n‑j<sub>2</sub>;(i)依次求开卷第k层时钢卷内任意第i层第s条的半径r<sub>si</sub>、第s条第i层与第i‑1层间钢卷内部径向应力p<sub>si</sub>、周向应力q<sub>si</sub>、层间摩擦力τ<sub>si</sub>,开卷松卷及层间滑移的控制目标函数<img file="FDA0000570635740000012.GIF" wi="76" he="63" />主要包括以下步骤:i1)计算开卷第k层时钢卷内任意第i层第s条的钢卷半径r<sub>si</sub>=r<sub>b</sub>+ih<sub>s</sub>;i2)计算开卷第k层时带钢在第s条、第i层与第i‑1层间钢卷内部径向应力p<sub>si</sub>,其中有相关的关系式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>p</mi><mi>si</mi></msub><mo>[</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><msub><mi>&nu;</mi><mn>2</mn></msub><mn>2</mn></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>r</mi><mi>si</mi></msub><mo>+</mo><mfrac><msub><mi>&nu;</mi><mn>2</mn></msub><mn>2</mn></mfrac><msub><mi>r</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>]</mo><mo>-</mo><msub><mi>p</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>[</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><msub><mi>&nu;</mi><mn>2</mn></msub><mn>2</mn></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>r</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>+</mo><mfrac><msub><mi>&nu;</mi><mn>2</mn></msub><mn>2</mn></mfrac><msub><mi>r</mi><mi>si</mi></msub><mo>]</mo><mo>=</mo><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mi>si</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><mrow><msub><mi>r</mi><mi>si</mi></msub><mo>-</mo><msub><mi>r</mi><mrow><mn>01</mn><mi>si</mi></mrow></msub></mrow><msub><mi>r</mi><mrow><mn>01</mn><mi>si</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000570635740000021.GIF" wi="1630" he="160" /></maths>式中:r<sub>si+1</sub>开卷当前层时第i+1层带钢在第s条处的内半径,p<sub>si+1</sub>开卷当前层时带钢在第s条、第i+1层与第i层间钢卷内部径向应力,r<sub>01si</sub>开卷当前层的前一层时第i层带钢在第s条处的内半径;i3)计算开卷第k层时的卷筒位移<img file="FDA0000570635740000022.GIF" wi="696" he="165" />i4)计算开卷第k层考虑到层间变形时任意第i层第s条的钢卷半径<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><msub><mi>r</mi><mi>si</mi></msub><mo>&prime;</mo></msup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>[</mo><msub><mi>E</mi><mn>2</mn></msub><msub><mi>h</mi><mi>s</mi></msub><mo>+</mo><mfrac><mrow><msub><mi>&nu;</mi><mn>2</mn></msub><msub><mi>p</mi><mi>si</mi></msub><msub><mi>r</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mn>2</mn></mfrac><mo>+</mo><mrow><mo>(</mo><mfrac><msub><mi>&nu;</mi><mn>2</mn></msub><mn>2</mn></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>p</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub><msub><mi>r</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><mn>0.5</mn><mi>&rho;g</mi><msub><mi>h</mi><mi>s</mi></msub><mi>sin</mi><mi>&theta;</mi><mo>]</mo></mrow><mrow><mo>[</mo><mfrac><mrow><msub><mi>E</mi><mn>2</mn></msub><msub><mi>h</mi><mi>s</mi></msub></mrow><msub><mi>r</mi><mrow><mn>01</mn><mi>si</mi></mrow></msub></mfrac><mo>+</mo><mfrac><mrow><msub><mi>&nu;</mi><mn>2</mn></msub><msub><mi>p</mi><mrow><mi>si</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mn>2</mn></mfrac><mo>+</mo><mrow><mo>(</mo><mfrac><msub><mi>&nu;</mi><mn>2</mn></msub><mn>2</mn></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>p</mi><mi>si</mi></msub><mo>+</mo><mn>0.5</mn><mi>&rho;g</mi><msub><mi>h</mi><mi>s</mi></msub><mi>sin</mi><mi>&theta;</mi><mo>]</mo></mrow></mfrac></mtd><mtd><mi>i</mi><mo>></mo><mn>1</mn></mtd></mtr><mtr><mtd><msub><mi>r</mi><mi>b</mi></msub><mo>-</mo><msub><mi>u</mi><mi>s</mi></msub></mtd><mtd><mi>i</mi><mo>=</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000570635740000023.GIF" wi="1231" he="358" /></maths>θ为极角,g为重力常数;i5)判断不等式||r<sub>si</sub>‑r'<sub>si</sub>||≤ε是否成立,如果不等式不成立则令r<sub>si</sub>=r'<sub>si</sub>,转入步骤i2),如果不等式成立,则直接转入步骤i6);i6)计算开卷第k层时带钢在第s条、第i层与第i‑1层间钢卷周向应力q<sub>si</sub>=q<sub>si+1</sub>‑π(r<sub>si</sub>τ<sub>si</sub>‑r<sub>si+1</sub>τ<sub>si+1</sub>)/(r<sub>si+1</sub>‑r<sub>si</sub>),式中:q<sub>si+1</sub>开卷第k层时带钢在第s条、第i+1层与第i层间钢卷周向应力,τ<sub>si+1</sub>开卷当前层时第i+1层带钢与第i层带钢在第s条处的摩擦应力;i7)计算开卷第k层时第i层带钢与第i‑1层带钢在第s条处的紧密系数<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>m</mi><mi>si</mi></msub><mo>=</mo><mn>1</mn><mo>+</mo><mn>1.27</mn><msub><mi>k</mi><mn>0</mn></msub><mfrac><msub><mi>&delta;</mi><mi>max</mi></msub><msub><mi>h</mi><mi>s</mi></msub></mfrac><msup><mrow><mo>(</mo><mfrac><msub><mi>E</mi><mn>2</mn></msub><msub><mi>p</mi><mi>si</mi></msub></mfrac><mo>)</mo></mrow><mn>0.8</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0000570635740000024.GIF" wi="605" he="168" /></maths>i8)计算开卷第k层时带钢在第s条、第i层与第i‑1层间钢卷内部层与层间摩擦应力的τ<sub>si</sub>=r<sub>si+1</sub>τ<sub>si+1</sub>/r<sub>si</sub>+(q<sub>si+1</sub>‑q<sub>si</sub>)(r<sub>si+1</sub>‑r<sub>si</sub>)/πr<sub>si</sub>;i9)求解开卷第k层时的最小摩擦平衡裕量K<sub>l</sub>=min{μp<sub>si</sub>‑τ<sub>si</sub>};i10)求解开卷第k层时的开卷松卷及层间滑移的控制目标函数<img file="FDA0000570635740000031.GIF" wi="894" he="264" />(j)判断不等式<img file="FDA0000570635740000032.GIF" wi="136" he="61" />是否成立?如果不等式成立则令<img file="FDA0000570635740000033.GIF" wi="164" he="59" />j<sub>1y</sub>=j<sub>1</sub>,转入步骤(k);如果不等式不成立,则直接转入步骤(k);(k)判断不等式k<2是否成立?如果不等式成立则令j<sub>1</sub>=j<sub>1</sub>+1,T=T<sub>min</sub>+j<sub>1</sub>·ΔT,转入步骤(l);如果不等式不成立则令j<sub>2</sub>=j<sub>2</sub>+1,转入步骤(h);(l)判断不等式T≤T<sub>max</sub>是否成立?如果不等式成立则转入步骤(f);如果不等式不成立则转入步骤(m);(m)输出最优开卷张力值T<sub>y</sub>=T<sub>min</sub>+j<sub>1y</sub>·ΔT。
地址 066004 河北省秦皇岛市海港区河北大街西段438号