发明名称 一种冷轧薄带钢宽度方向起浪形式的预测方法
摘要 本发明提供了一种冷轧薄带钢宽度方向起浪形式的预测方法。该方法在得到冷轧带钢残余应力值及其坐标值的基础上,经如下步骤实现:(1)计算左辅助向量和右辅助向量;(2)计算包括端点值的残余应力正极大值坐标向量;(3)计算包括端点值的残余应力的负极小值向量及其坐标向量;(4)计算宽度向量;(5)计算幅值向量;(6)确定宽度方向起浪形式函数。利用该方法可以预测冷轧薄带钢宽度方向的起浪区域和起浪形式,提高冷轧生产中对带钢质量的评估能力。
申请公布号 CN104923573A 申请公布日期 2015.09.23
申请号 CN201510254770.3 申请日期 2015.05.18
申请人 北京科技大学 发明人 李洪波;李秾;张杰;曹建国
分类号 B21B37/28(2006.01)I 主分类号 B21B37/28(2006.01)I
代理机构 北京金智普华知识产权代理有限公司 11401 代理人 皋吉甫
主权项 一种冷轧薄带钢宽度方向起浪形式的预测方法,其特征在于:生产现场板形仪测得的冷轧薄带钢某一横截面上残余应力检测值对应向量σ=(σ<sub>1</sub>,σ<sub>2</sub>,σ<sub>3</sub>,…,σ<sub>n‑2</sub>,σ<sub>n‑1</sub>,σ<sub>n</sub>),采用公式1得到左辅助向量σ<sup>left</sup>和右辅助向量σ<sup>right</sup><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msup><mi>&sigma;</mi><mi>left</mi></msup><mo>=</mo><mrow><mo>(</mo><msub><mi>&sigma;</mi><mn>2</mn></msub><mo>,</mo><msub><mi>&sigma;</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&sigma;</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&sigma;</mi><mrow><mi>n</mi><mo>-</mo><mn>3</mn></mrow></msub><mo>,</mo><msub><mi>&sigma;</mi><mrow><mi>n</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>,</mo><msub><mi>&sigma;</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msup><mi>&sigma;</mi><mi>right</mi></msup><mo>=</mo><mrow><mo>(</mo><msub><mi>&sigma;</mi><mn>2</mn></msub><mo>,</mo><msub><mi>&sigma;</mi><mn>3</mn></msub><mo>,</mo><msub><mi>&sigma;</mi><mn>4</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&sigma;</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&sigma;</mi><mi>n</mi></msub><mo>,</mo><msub><mi>&sigma;</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mn>1</mn></mrow>]]></math><img file="FDA0000719450570000011.GIF" wi="1520" he="170" /></maths>式中:σ<sub>1</sub>,σ<sub>2</sub>,σ<sub>3</sub>,…,σ<sub>n‑2</sub>,σ<sub>n‑1</sub>,σ<sub>n</sub>为检测得到的带钢残余应力在同一横截面上不同点的检测值,也就是向量σ的分量;σ<sup>left</sup>为左辅助向量;σ<sup>right</sup>为右辅助向量;n为检测点个数,采用公式2得到包括端点值的残余应力正极大值坐标向量c<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>u</mi><mo>=</mo><mi>N</mi><mo>{</mo><mi>U</mi><mo>[</mo><mi>U</mi><mrow><mo>(</mo><mi>&sigma;</mi><mo>-</mo><msup><mi>&sigma;</mi><mi>left</mi></msup><mo>)</mo></mrow><mo>.</mo><mo>*</mo><mi>U</mi><mrow><mo>(</mo><mi>&sigma;</mi><mo>-</mo><msup><mi>&sigma;</mi><mi>right</mi></msup><mo>)</mo></mrow><mo>.</mo><mo>*</mo><mi>&sigma;</mi><mo>]</mo><mo>.</mo><mo>*</mo><mrow><mo>(</mo><mi>r</mi><mo>+</mo><msub><mi>o</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>}</mo></mtd></mtr><mtr><mtd><mi>c</mi><mo>=</mo><mi>u</mi><mo>-</mo><msub><mi>o</mi><mi>u</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mn>2</mn></mrow>]]></math><img file="FDA0000719450570000012.GIF" wi="1523" he="164" /></maths>式中:运算U表示将一个向量中大于0的分量置1,小于等于0的分量置0,进而得到一个新的向量;运算.*表示将两个相同维数的向量对应元素相乘,进而得到一个新的向量,其运算优先度低于算子U;r为冷轧薄带钢某一横截面上残余应力检测点坐标对应的向量;o<sub>r</sub>为与r维数相同且分量全为1的向量;运算N表示将一个向量中零分量全部去掉后,剩余分量构成一个新的向量;u为计算过程中的一个中间向量;c为包括端点值的残余应力正极大值坐标向量;o<sub>u</sub>为与u维数相同且分量全为1的向量,以相同的方法,采用公式3得到包括端点值的残余应力负极小值向量s及其坐标向量b<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>s</mi><mo>=</mo><mi>N</mi><mo>{</mo><mi>U</mi><mo>[</mo><mi>U</mi><mrow><mo>(</mo><mo>-</mo><mi>&sigma;</mi><mo>+</mo><msup><mi>&sigma;</mi><mi>left</mi></msup><mo>)</mo></mrow><mo>.</mo><mo>*</mo><mi>U</mi><mrow><mo>(</mo><mo>-</mo><mi>&sigma;</mi><mo>+</mo><msup><mi>&sigma;</mi><mi>right</mi></msup><mo>)</mo></mrow><mo>.</mo><mo>*</mo><mrow><mo>(</mo><mo>-</mo><mi>&sigma;</mi><mo>)</mo></mrow><mo>]</mo><mo>.</mo><mo>*</mo><mi>&sigma;</mi><mo>}</mo></mtd></mtr><mtr><mtd><mi>w</mi><mo>=</mo><mi>N</mi><mo>{</mo><mi>U</mi><mo>[</mo><mi>U</mi><mrow><mo>(</mo><mo>-</mo><mi>&sigma;</mi><mo>+</mo><msup><mi>&sigma;</mi><mi>left</mi></msup><mo>)</mo></mrow><mo>.</mo><mo>*</mo><mi>U</mi><mrow><mo>(</mo><mo>-</mo><mi>&sigma;</mi><mo>+</mo><msup><mi>&sigma;</mi><mi>right</mi></msup><mo>)</mo></mrow><mo>.</mo><mo>*</mo><mrow><mo>(</mo><mo>-</mo><mi>&sigma;</mi><mo>)</mo></mrow><mo>]</mo><mo>.</mo><mo>*</mo><mrow><mo>(</mo><mi>r</mi><mo>+</mo><msub><mi>o</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>}</mo></mtd></mtr><mtr><mtd><mi>b</mi><mo>=</mo><mi>w</mi><mo>-</mo><msub><mi>o</mi><mi>w</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mn>3</mn></mrow>]]></math><img file="FDA0000719450570000021.GIF" wi="1538" he="261" /></maths>式中:s为包括端点值的残余应力负极小值向量;w为计算过程中的一个中间向量;b为包括端点值的残余应力负极小值坐标向量;o<sub>w</sub>为与w维数相同且分量全为1的向量,采用公式4得到宽度向量a<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mi>min</mi><mo>[</mo><mi>abs</mi><mrow><mo>(</mo><msubsup><mi>o</mi><mi>c</mi><mi>T</mi></msubsup><mi>b</mi><mo>-</mo><msup><mi>c</mi><mi>T</mi></msup><msub><mi>o</mi><mi>b</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mn>4</mn></mrow>]]></math><img file="FDA0000719450570000022.GIF" wi="1548" he="80" /></maths>式中:min表示将矩阵每一列的最小值取出后,按原有的列顺序排列构成的行向量;abs表示将矩阵每一元素取绝对值;a为宽度向量;o<sub>c</sub>为与c维数相同且分量全为1的向量;<img file="FDA0000719450570000024.GIF" wi="65" he="73" />为o<sub>c</sub>的转置;c<sup>T</sup>为c的转置;o<sub>b</sub>为与b维数相同且分量全为1的向量;<img file="FDA0000719450570000025.GIF" wi="65" he="83" />与b之间做矩阵乘法;c<sup>T</sup>与o<sub>b</sub>之间做矩阵乘法,采用公式5得到幅值向量ε<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&epsiv;</mi><mo>=</mo><mo>-</mo><mfrac><mrow><mn>3</mn><mi>&pi;</mi></mrow><mrow><mn>4</mn><mi>E</mi></mrow></mfrac><mi>a</mi><mo>.</mo><mo>*</mo><mi>s</mi><mo>-</mo><mo>-</mo><mo>-</mo><mn>5</mn></mrow>]]></math><img file="FDA0000719450570000023.GIF" wi="1556" he="130" /></maths>式中:π为圆周率;E为钢的弹性模量2.1×10<sup>5</sup>MPa;ε为幅值向量,采用公式6确定宽度方向起浪形式函数δ(x)<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>&delta;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><msub><mi>&epsiv;</mi><mi>j</mi></msub><msubsup><mi>a</mi><mi>j</mi><mn>4</mn></msubsup></mfrac><mtext>H</mtext><mrow><mo>(</mo><msub><mi>a</mi><mi>j</mi></msub><mo>-</mo><msub><mi>b</mi><mi>j</mi></msub><mo>+</mo><mi>x</mi><mo>)</mo></mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>j</mi></msub><mo>+</mo><msub><mi>b</mi><mi>j</mi></msub><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mo>[</mo><msubsup><mi>a</mi><mi>j</mi><mn>2</mn></msubsup><mo>-</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>b</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><mo>]</mo><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mn>6</mn></mrow>]]></math><img file="FDA0000719450570000031.GIF" wi="1539" he="179" /></maths>式中:m为向量a、b、ε的维数;ε<sub>j</sub>为向量ε的分量;a<sub>j</sub>为向量a的分量;b<sub>j</sub>为向量b的分量;x为沿带钢宽度方向的自变量,范围从零到带钢宽度;δ(x)为带钢宽度方向起浪形式函数;H(x)为Heaviside函数,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><mi>x</mi><mo>></mo><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo></mtd><mtd><mi>x</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><mi>x</mi><mo>&lt;</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>.</mo></mrow>]]></math><img file="FDA0000719450570000032.GIF" wi="476" he="294" /></maths>
地址 100083 北京市海淀区学院路30号