一种四辊轧机负载辊缝形状预报方法技术

一种四辊轧机负载辊缝形状预报方法，它主要包括以下由计算机执行的步骤：1、给定四辊轧机设备参数及工艺参数；2、给定轧制压力横向分布多项式；3、计算相关参数，初始化辊间压力；4、根据辊间变形协调条件及工作辊力与力矩平衡条件形成线性方程组；5、求解线性方程组；6、计算四辊轧机负载辊缝形状。本发明专利技术采用高阶多项式描述辊间压力与轧制压力的横向分布，采用莫尔积分法计算出多项式分布力作用下轧辊挠度的解析函数，最终使辊间压力求解方程组阶数较常规影响函数法显著减小，使计算速度的明显减少，计算精度显著提高，是一种实用高效的负载辊缝形状预报方法。

【技术实现步骤摘要】

本专利技术属于冶金轧制金属领域，特别涉及一种辊缝形状的预报方法。
技术介绍
板带轧制过程中，乳件的出口断面形状就是负载辊缝的形状，板形控制手段都是 通过改变负载辊缝形状实现板形控制的。因此，实现负载辊缝形状的高效预报是进行高精 度板形控制的先决条件。 目前研宄辊系弹性变形的方法中，影响函数法是一种比较实用的计算方法。其基 本思路是将辊系沿轴向划分一定数量的单元，假设每个单元上承受均布压力，即通过分段 均布压力逼近真实的曲线分布力，因此其计算精度与单元划分数直接相关。用于在线计算 时，计算速度与计算精度存在冲突，当单元划分较多时，虽然计算精度缓慢提高，但计算时 间显著增加，无法满足在线计算要求，当单元划分较少时，虽然计算时间显著降低，计算精 度又明显下降，无法得到真实最优的计算结果。
技术实现思路
本专利技术的目的在于提供一种既能减少计算时间，又能提高计算精度的四辊轧机负 载辊缝形状预报方法。本专利技术主要是采用高阶多项式描述辊间压力与轧制压力的横向分 布，应用莫尔积分法计算多项式分布力作用下轧辊挠度函数，最终使辊间压力求解方程组 阶数较常规影响函数法显著减小，在保证计算精度的前提下，实现计算时间的显著降低。 本专利技术四辊轧机负载辊缝形状预报方法包括以下由计算机执行的步骤： a、以左压下支点为坐标原点，给定四辊轧机设备参数及工艺参数，具体包括：工作 辊辊身长度L w、支撑辊辊身长度Lb、工作辊弯辊缸间距Lfw、压下油缸距离Ls、工作辊直径横 向分布函数D w (y)、支撑辊直径横向分布函数Db (y)、支撑辊辊颈直径db、工作辊弹性模量Ew、 工作辊泊松比uw、支撑辊弹性模量E b、支撑辊泊松比u b、工作弯辊力Fw、工作辊窜辊量Sw、 棍间压力左边界坐标y bl、间压力右边界坐标ybr、间压力多项式最高次数n、计算收敛精 度e ; b、给定轧制压力横向分布多项式，具体包括：乳制压力左边界坐标ywl，乳制压力 c、计算相关参数，初始化辊间压力，具体如下： cl、计算支撑辊辊身惯性矩<、支撑辊辊颈惯性矩/【、支撑辊辊身横截面积、 支撑辊辊颈横截面积4、工作辊辊身惯性矩、工作辊辊身横截面积;c3、计算辊间压力多项式各次项系数初始值< - a+，令 al=(P+ 2Fw)/{ybr-ybl),O1=O2=......= 〇:=〇,令临时变量％ = a0% C4、令i=〇; c5、令辊间压力多项式的i次项系数a,.= < ; c6、判断i = n是否成立，若成立，则进入d ;若不成立，令i = i+1，返回c5 ; d、将辊间压力多项式各次项系数a(l~an、工作辊辊身左端部刚性位移(^与工作辊 辊身右端部刚性位移C 2作为未知数，根据辊间变形协调条件及工作辊力与力矩平衡条件形 成线性方程组，具体包括： dl、令 r = 0 ; d3、令i=0; d4、根据莫尔积分法计算辊间压力多项式中i次项在y,点引起的支撑辊挠度影响 系数 a b(i，yJ ; d5、根据莫尔积分法计算辊间压力多项式中i次项在^点引起的工作辊挠度影响 系数；C 6 (/，兄.）； d6、判断i = n是否成立，若成立，则进入d7 ;若不成立，令i = i+1，返回d4 ; d7、令 i = 0 ; d8、根据莫尔积分法计算轧制压力多项式中i次项在^点引起的工作辊挠度影响 系数 d9、判断i = m是否成立，若成立，贝丨」进入dlO ;若不成立，令i = i+1，返回d8 ; dlO、根据材料力学计算工作棍弯棍力在yr引起的工作棍烧曲位移乂f (.Vr); dll、计算工作辊与支撑辊间空载间隙 ADwb(yr) = /2 ； dl2、判断r = n是否成立，若成立，则进入dl3 ;若不成立，令r = r+1，返回d2 ; dl3、根据工作辊力平衡条件建立线性方程如下： dl8、根据弹性半平面理论计算辊间压力多项式中i次项在y,点引起的辊间弹性 压扁影响系数x wb(i, yj ; dl9、判断i = n是否成立，若成立，贝丨」进入d20;若不成立，令i = i+l，返回dl8; d20、根据辊间变形协调条件形成h点处的线性方程： d21、判断r = n是否成立，若成立，则进入e;若不成立，令r = r+1，返回dl6; e、求解线性方程组，判断是否收敛，具体如下： el、利用LU分解法求解步骤dl3、步骤dl4与步骤d20形成的n+3个方程，得到未 知数a。~a n、(^与C 2共n+3个未知数； e2、令 r = 0，q = 0; e6、令 i = 0; el、令O1 =Cii; e8、判断i = n是否成立，若成立，进入dl5。若不成立，令i = i+1，返回e7; f、计算四辊轧机负载辊缝形状，具体包括： H、给定需要计算负载辊缝形状的点数num，令j=1 ; f2、给定需要计算负载辊缝形状的坐标yj; f3、根据莫尔积分法计算辊间压力多项式、乳制压力多项式与工作辊弯辊力在 7」 点共同引起的工作辊挠度fw(yP ; f4、根据弹性半空间理论计算轧制压力在7」点引起的工作辊弹性压扁量 5Ws(Yj)； f5、计算乳制压力y」点的工作棍直径Dw(y」）； f6、判断j = num是否成立，若成立，则进入f7 ;若不成立，令j = j+1，返回f2 ;f7、令 j = I; f8、计算y」处负载辊缝形状： s(yj)= 2fw (Ls/2) - +2 8 ws(Ls/2)- + /2_DW (Ls/2); f9、判断j = num是否成立，若成立，计算结束；若不成立，令j = j+1，返回f8。 本专利技术与现有技术相比具有如下优点： 采用高阶多项式描述辊间压力与轧制压力的横向分布，采用莫尔积分法计算出多 项式分布力作用下轧辊挠度的解析函数，最终使辊间压力求解方程组阶数较常规影响函数 法显著减小，使计算时间明显减少，计算速度显著提高，是一种实用高效的负载辊缝形状预 报方法。【附图说明】 图1为本专利技术的计算总流程图。 图2为本专利技术步骤d的计算流程图。 图3为本专利技术步骤f的计算流程图。图4为本专利技术实施例1中计算得到的辊间压力分布与普通影响函数法计算得到的 辊间压力分布对比图。 图5为实施例1中计算得到的负载辊缝形状与普通影响函数法计算得到的负载辊 缝形状对比图。 图6为实施例2中计算得到的辊间压力分布与普通影响函数法计算得到的辊间压 力分布对比图。 图7为实施例2中计算得到的负载辊缝形状与普通影响函数法计算得到的负载辊 缝形状对比图。【具体实施方式】 实施例1 对某实际四辊轧机的负载辊缝形状进行预报，其计算流程如图1所示，首先在步 骤a中，以左压下支点为坐标原点，给定四辊轧机设备参数及工艺参数，具体包括：工作辊 棍身长度L w= 2000mm、支撑棍棍身长度Lb= 1750mm、工作棍弯棍缸间距Lfw= 3000mm、压 下油缸距离Ls= 2900mm、工作棍直径横向分布函数Dw (y) = 700mm、支撑棍直径横向分布函 数Db(y) = 1500mm、支撑辊辊颈直径db= 955mm、工作辊弹性模量Ew= 210000MPa、工作辊 泊松比Uw= 0. 3、支撑辊弹性模量Eb= 210000MPa、支本文档来自技高网...

【技术保护点】
一种四辊轧机负载辊缝形状预报方法，其特征在于：它包括以下由计算机执行的步骤：a、以左压下支点为坐标原点，给定四辊轧机设备参数及工艺参数，具体包括：工作辊辊身长度Lw、支撑辊辊身长度Lb、工作辊弯辊缸间距Lfw、压下油缸距离Ls、工作辊直径横向分布函数Dw(y)、支撑辊直径横向分布函数Db(y)、支撑辊辊颈直径db、工作辊弹性模量Ew、工作辊泊松比υw、支撑辊弹性模量Eb、支撑辊泊松比υb、工作弯辊力Fw、工作辊窜辊量Sw、辊间压力左边界坐标ybl、辊间压力右边界坐标ybr、辊间压力多项式最高次数n、计算收敛精度ε；b、给定轧制压力横向分布多项式，具体包括：轧制压力左边界坐标ywl，轧制压力右边界坐标ywr，轧制压力多项式的最高次数m，轧制压力多项式p(y)=Σj=0mbjyj;]]>c、计算相关参数，初始化辊间压力，具体如下：c1、计算支撑辊辊身惯性矩支撑辊辊颈惯性矩支撑辊辊身横截面积支撑辊辊颈横截面积工作辊辊身惯性矩工作辊辊身横截面积c2、计算总轧制力P=Σj=0mbj(ywrj+1-ywlj+1)j+1;]]>c3、计算辊间压力多项式各次项系数初始值令a1*=a2*=......=an*=0,]]>令临时变量q0=a0*;]]>c4、令i＝0；c5、令辊间压力多项式的i次项系数c6、判断i＝n是否成立，若成立，则进入d；若不成立，令i＝i+1，返回c5；d、将辊间压力多项式各次项系数a0～an、工作辊辊身左端部刚性位移C1与工作辊辊身右端部刚性位移C2作为未知数，根据辊间变形协调条件及工作辊力与力矩平衡条件形成线性方程组，具体包括：d1、令r＝0；d2、计算yr=ybl+r(ybr-ybl)n;]]>d3、令i＝0；d4、根据莫尔积分法计算辊间压力多项式中i次项在yr点引起的支撑辊挠度影响系数αb(i,yr)；d5、根据莫尔积分法计算辊间压力多项式中i次项在yr点引起的工作辊挠度影响系数d6、判断i＝n是否成立，若成立，则进入d7；若不成立，令i＝i+1，返回d4；d7、令i＝0；d8、根据莫尔积分法计算轧制压力多项式中i次项在yr点引起的工作辊挠度影响系数d9、判断i＝m是否成立，若成立，则进入d10；若不成立，令i＝i+1，返回d8；d10、根据材料力学计算工作辊弯辊力在yr引起的工作辊挠曲位移d11、计算工作辊与支撑辊间空载间隙△Dwb(yr)＝[Dw(Ls/2)‑Dw(yr)+Db(Ls/2)‑Db(yr)]/2；d12、判断r＝n是否成立，若成立，则进入d13；若不成立，令r＝r+1，返回d2；d13、根据工作辊力平衡条件建立线性方程如下：Σi=0nai(ybri+1-ybli+1)i+1=Σj=0nbj(ywri+1-ywli+1)j+1+2Fw;]]>d14、根据工作辊力矩平衡条件建立线性方程如下：Σi=0nai(ybri+2-ybli+2)i+2=Σj=0mbj(ywrj+2-ywlj+2)j+2+FwLs;]]>d15、令r＝0；d16、计算yr=ybl+r(ybr-ybl)n;]]>d17、令i＝0；d18、根据弹性半平面理论计算辊间压力多项式中i次项在yr点引起的辊间弹性压扁影响系数λwb(i,yr)；d19、判断i＝n是否成立，若成立，则进入d20；若不成立，令i＝i+1，返回d18；d20、根据辊间变形协调条件形成yr点处的线性方程：Σi=0nai(αb(i,yr)+γwwb(i,yr)+λwb(i,yr))+yr-Lw-DLwC1+D-yrLwC2=Σj=0mbjγwp(i,yr)-fwF(yr)-ΔDwb(yr);]]>d21、判断r＝n是否成立，若成立，则进入e；若不成立，令r＝r+1，返回d16；e、求解线性方程组，判断是否收敛，具体如下：e1、利用LU分解法求解步骤d13、步骤d14与步骤d20形成的n+3个方程，得到未知数a0～an、C1与C2共n+3个未知数；e2、令r＝0，q＝0；e3、计算q=q+(Σi=0naiyri-Σi=0nai*yri)2;]]>e4、判断r＝n是否成立，若成立，进入e5；若不成立，令r＝r+1，返回e3；e5、判断是否成立，若成立，进入f；若不成立，则进入e6；e6、令i＝0；e7、令ai*=ai;]]>e8、判断i＝n是否成立，若成立，进入d15，若不成立，令i＝i+1，返回e7；f、计算四辊轧机负载辊缝形状，具体包括：f1、给定需要计算负载辊缝形状的点数num，令j＝1；f2、给定需要计算负载辊缝形状的坐标yj；f3、根据莫尔积分法计算辊间压力多项式、轧制压力多项式与工作...

【技术特征摘要】

【专利技术属性】
技术研发人员：王东城，吴燕林，张龙，刘宏民，
申请(专利权)人：燕山大学，
类型：发明
国别省市：河北;13