发明名称 一种带钢轧制在线控制摩擦系数模型优化方法
摘要 本发明涉及一种带钢轧制在线控制摩擦系数模型优化方法,将带钢轧制过程的实测轧制力代入轧制力理论计算模型中。由于轧制力是摩擦系数的函数,因此通过变式处理得到含有以摩擦系数为自变量的非线性方程。除了将实测轧制力值代入轧制力理论计算模型外,还要将轧制过程中带钢的实际入口厚度、实际出口厚度、实际张力等数据同时代入轧制力计算模型。通过求解该非线性方程获得摩擦系数的实际计算值,以实际摩擦系数为基础,对摩擦系数理论计算模型进行自适应优化。通过该方法可以提高摩擦系数在线设定精度。
申请公布号 CN103191919A 申请公布日期 2013.07.10
申请号 CN201210002100.9 申请日期 2012.01.05
申请人 鞍钢股份有限公司 发明人 王军生;候永钢;宋蕾;张岩;刘宝权;秦大伟;宋君;费静;吴萌
分类号 B21B1/28(2006.01)I;B21B38/08(2006.01)I;B21B38/06(2006.01)I;B21B38/02(2006.01)I 主分类号 B21B1/28(2006.01)I
代理机构 鞍山嘉讯科技专利事务所 21224 代理人 张群
主权项 一种带钢轧制在线控制摩擦系数模型优化方法,其特征在于,该方法包括以下步骤:1)在线检测带钢轧制过程的实测轧制力、实测带钢轧制入口厚度及出口厚度、实测带钢轧制过程的前后张力,将上述工艺参数代入轧制力计算模型中,构成以摩擦系数为未知量的非线性方程,求解该方程得到轧制过程摩擦系数的后计算值,将该摩擦系数的后计算值作为实际摩擦系数值;所述的轧制力计算模型为: <mrow> <msup> <msub> <mi>&mu;</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <mi>Fi</mi> </msub> <mo>-</mo> <mn>1.08</mn> <mo>+</mo> <mn>1.02</mn> <mo>&CenterDot;</mo> <msup> <msub> <mi>r</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> </mrow> <mrow> <mn>1.79</mn> <mo>&CenterDot;</mo> <msup> <msub> <mi>r</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <msqrt> <mn>1</mn> <mo>-</mo> <msup> <msub> <mi>r</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> </msqrt> <msqrt> <mfrac> <msubsup> <mi>R</mi> <mi>i</mi> <mrow> <mo>&prime;</mo> <mi>mdl</mi> </mrow> </msubsup> <msup> <msub> <mi>h</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> </mfrac> </msqrt> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>Q</mi> <mi>Fi</mi> </msub> <mo>=</mo> <mfrac> <msup> <msub> <mi>F</mi> <mi>i</mi> </msub> <mi>measure</mi> </msup> <mrow> <mo>{</mo> <msup> <msub> <mi>k</mi> <mi>mi</mi> </msub> <mi>mdl</mi> </msup> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> <msup> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mi>back</mi> </msup> <mo>}</mo> <mi>W</mi> <msqrt> <msup> <msub> <mi>R</mi> <mi>i</mi> </msub> <mrow> <mo>&prime;</mo> <mi>mdl</mi> </mrow> </msup> <mo>&CenterDot;</mo> <mi>&Delta;</mi> <msup> <msub> <mi>h</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> </msqrt> </mrow> </mfrac> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <msqrt> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>v</mi> <mn>2</mn> </msup> </mrow> <mi>E</mi> </mfrac> <msup> <msub> <mi>k</mi> <mi>mi</mi> </msub> <mi>mdl</mi> </msup> <mfrac> <msup> <msub> <mi>h</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mrow> <mi>&Delta;</mi> <msup> <msub> <mi>h</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> </mrow> </mfrac> </msqrt> </mrow>式中:μiback摩擦系数后计算值实际值;riback压下率后计算实际值;R′mdl轧辊压扁半径模型计算值;hiback带钢厚度后计算实际值;Fimeasure轧制力实测值;kmimdl带钢变形抗力模型计算值;tiback带钢张力后计算实际值;Δhiback带钢压下量后计算实际值;i机架号;2)将上述的摩擦系数后计算值作为实际摩擦系数值,与摩擦系数模型计算值进行比较,计算出用于补偿摩擦系数模型计算值与实际值偏差的模型自适应学习系数,通过该模型自适应学习系数可以提高轧制过程摩擦系数模型的计算精度;所述的摩擦系数模型为:μ=μL·μv+μmin                         (5)(1)轧制长度影响项 <mrow> <msub> <mi>&mu;</mi> <mi>L</mi> </msub> <mo>=</mo> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>&CenterDot;</mo> <msub> <mi>a</mi> <mi>L</mi> </msub> <mo>&CenterDot;</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>b</mi> <mi>L</mi> </msub> <mo>&CenterDot;</mo> <mi>L</mi> <mo>+</mo> <msub> <mi>c</mi> <mi>L</mi> </msub> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>d</mi> <mi>L</mi> </msub> <mo>&CenterDot;</mo> <msup> <mi>L</mi> <mn>3</mn> </msup> <mo>+</mo> <msub> <mi>e</mi> <mi>L</mi> </msub> <mo>&CenterDot;</mo> <msup> <mi>L</mi> <mn>4</mn> </msup> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>(2)速度影响项 <mrow> <msub> <mi>&mu;</mi> <mi>v</mi> </msub> <mo>=</mo> <msub> <mi>a</mi> <mi>v</mi> </msub> <mo>&CenterDot;</mo> <mo>[</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>C</mi> <mn>2</mn> </msub> <mo>{</mo> <msub> <mi>b</mi> <mi>v</mi> </msub> <mrow> <mo>(</mo> <mi>v</mi> <mo>-</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mi>v</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>v</mi> <mo>-</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>d</mi> <mi>v</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>v</mi> <mo>-</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>3</mn> </msup> <mo>+</mo> <msub> <mi>e</mi> <mi>v</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>v</mi> <mo>-</mo> <msub> <mi>v</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>4</mn> </msup> <mo>}</mo> <mo>|</mo> </mrow> </msup> <mo>]</mo> <mo>+</mo> <msub> <mi>f</mi> <mi>v</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>μ     摩擦系数L      工作辊轧制长度v      出口带钢速度μmin  摩擦系数基准值aL~eL 轧制长度影响系数av~fv 速度影响系数v0     标准速度C1     轧制长度影响项的自适应学习系数C2     速度影响项的自适应学习系数通过公式(13)得出的摩擦系数的后计算值即实际摩擦系数值,再结合摩擦系数模型公式(5)、公式(6)、公式(7),可以得出轧制长度影响项的自适应学习系数C1及速度影响项的自适应学习系数C2;C1i、C2i为第i对机架的自适应学习系数; <mrow> <msup> <msub> <mi>C</mi> <mrow> <mn>1</mn> <mi>i</mi> </mrow> </msub> <mi>current</mi> </msup> <mo>=</mo> <mfrac> <mrow> <msup> <msub> <mi>&mu;</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>&mu;</mi> <mrow> <mi>min</mi> <mi>i</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&mu;</mi> <mi>vi</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>b</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>L</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>c</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <msub> <mi>L</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>d</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <msub> <mi>L</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <mo>+</mo> <msub> <mi>e</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <msub> <mi>L</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> </mrow> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&mu;</mi> <mi>vi</mi> </msub> <mo>=</mo> <msub> <mi>a</mi> <mi>vi</mi> </msub> <mo>&CenterDot;</mo> <mo>{</mo> <msup> <mi>e</mi> <mrow> <msup> <msub> <mi>c</mi> <mrow> <mn>2</mn> <mi>i</mi> </mrow> </msub> <mi>before</mi> </msup> <mo>[</mo> <msub> <mi>b</mi> <mi>vi</mi> </msub> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mrow> <mi>vi</mi> <mo>&CenterDot;</mo> </mrow> </msub> <msup> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>d</mi> <mi>vi</mi> </msub> <mo>&CenterDot;</mo> <msup> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>3</mn> </msup> <mo>+</mo> <msub> <mi>e</mi> <mi>vi</mi> </msub> <mo>&CenterDot;</mo> <msup> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>4</mn> </msup> <mo>]</mo> </mrow> </msup> <mo>}</mo> <mo>+</mo> <msub> <mi>f</mi> <mi>vi</mi> </msub> </mrow>极限检查:   C1ilower≤C1icurrent≤C1iupper                        (15)学习系数更新 C1inext=(1‑αC1i)·C1icurrent+αC1i·C1ibefore       (16) <mrow> <msup> <msub> <mi>C</mi> <mrow> <mn>2</mn> <mi>i</mi> </mrow> </msub> <mi>current</mi> </msup> <mo>=</mo> <mfrac> <mrow> <mi>ln</mi> <mo>{</mo> <mfrac> <mrow> <msup> <msub> <mi>&mu;</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>&mu;</mi> <mrow> <mi>min</mi> <mi>i</mi> </mrow> </msub> </mrow> <msub> <mi>&mu;</mi> <mi>Li</mi> </msub> </mfrac> <mo>-</mo> <msub> <mi>f</mi> <mi>vi</mi> </msub> <mo>}</mo> <mo>-</mo> <mi>ln</mi> <msub> <mi>a</mi> <mi>vi</mi> </msub> </mrow> <mrow> <mo>{</mo> <msub> <mi>b</mi> <mi>vi</mi> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mi>vi</mi> </msub> <msup> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>d</mi> <mi>vi</mi> </msub> <msup> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>3</mn> </msup> <mo>+</mo> <mi>e</mi> <msup> <mrow> <mo>(</mo> <msup> <msub> <mi>v</mi> <mi>i</mi> </msub> <mi>back</mi> </msup> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>4</mn> </msup> <mo>}</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>式中 <mrow> <msub> <mi>&mu;</mi> <mi>Li</mi> </msub> <mo>=</mo> <msup> <msub> <mi>C</mi> <mrow> <mn>1</mn> <mi>i</mi> </mrow> </msub> <mi>before</mi> </msup> <msub> <mi>a</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>b</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>L</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>c</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <msub> <mi>L</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>d</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <msub> <mi>L</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <mo>+</mo> <msub> <mi>e</mi> <mi>Li</mi> </msub> <mo>&CenterDot;</mo> <msup> <msub> <mi>L</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> </mrow> </msup> </mrow>极限检查:    C2ilower≤C2icurrent≤C2iupper                       (18)学习系数更新:C2inext=(1‑αC2i)·C2icurrent +αC2i·C2ibefore     (19)式中C1inew,C2inew‑下一次计算摩擦系数模型自适应学习系数    C1icurrent,C2icurrent‑当前卷摩擦系数模型自适应学习系数计算值    C1ibefore,C2ibefore‑上一次计算摩擦系数模型自适应学习系数    αC1i,αC2i‑摩擦系数模型自适应学习系数平滑因子i‑机架号(1、2、3、4、5......)3)将计算出的轧制长度影响项的自适应学习系数C1i及速度影响项的自适应学习系数C2i,代入摩擦系数模型公式(6)、公式(7),根据公式(5)计算出摩擦系数计算值。
地址 114000 辽宁省鞍山市铁西区环钢路1号
您可能感兴趣的专利