[发明专利]一种四辊轧机负载辊缝形状预报方法

权利要求书

1.一种四辊轧机负载辊缝形状预报方法，其特征在于：它包括以下由计算机执行的步骤：

a、以左压下支点为坐标原点，给定四辊轧机设备参数及工艺参数，具体包括：工作辊辊身长度L_w、支撑辊辊身长度L_b、工作辊弯辊缸间距L_fw、压下油缸距离L_s、工作辊直径横向分布函数D_w(y)、支撑辊直径横向分布函数D_b(y)、支撑辊辊颈直径d_b、工作辊弹性模量E_w、工作辊泊松比υ_w、支撑辊弹性模量E_b、支撑辊泊松比υ_b、工作弯辊力F_w、工作辊窜辊量S_w、辊间压力左边界坐标y_bl、辊间压力右边界坐标y_br、辊间压力多项式最高次数n、计算收敛精度ε；

b、给定轧制压力横向分布多项式，具体包括：轧制压力左边界坐标y_wl，轧制压力右边界坐标y_wr，轧制压力多项式的最高次数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、将辊间压力多项式各次项系数a₀～a_n、工作辊辊身左端部刚性位移C₁与工作辊辊身右端部刚性位移C₂作为未知数，根据辊间变形协调条件及工作辊力与力矩平衡条件形成线性方程组，具体包括：

d1、令r＝0；

d2、计算yr=ybl+r(ybr-ybl)n;]]>

d3、令i＝0；

d4、根据莫尔积分法计算辊间压力多项式中i次项在y_r点引起的支撑辊挠度影响系数α_b(i,y_r)；

d5、根据莫尔积分法计算辊间压力多项式中i次项在y_r点引起的工作辊挠度影响系数

d6、判断i＝n是否成立，若成立，则进入d7；若不成立，令i＝i+1，返回d4；

d7、令i＝0；

d8、根据莫尔积分法计算轧制压力多项式中i次项在y_r点引起的工作辊挠度影响系数

d9、判断i＝m是否成立，若成立，则进入d10；若不成立，令i＝i+1，返回d8；

d10、根据材料力学计算工作辊弯辊力在y_r引起的工作辊挠曲位移

d11、计算工作辊与支撑辊间空载间隙△D_wb(y_r)＝[D_w(L_s/2)-D_w(y_r)+D_b(L_s/2)-D_b(y_r)]/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次项在y_r点引起的辊间弹性压扁影响系数λ_wb(i,y_r)；

d19、判断i＝n是否成立，若成立，则进入d20；若不成立，令i＝i+1，返回d18；

d20、根据辊间变形协调条件形成y_r点处的线性方程：

Σ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个方程，得到未知数a₀～a_n、C₁与C₂共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、给定需要计算负载辊缝形状的坐标y_j；

f3、根据莫尔积分法计算辊间压力多项式、轧制压力多项式与工作辊弯辊力在y_j点共同引起的工作辊挠度f_w(y_j)；

f4、根据弹性半空间理论计算轧制压力在y_j点引起的工作辊弹性压扁量δ_ws(y_j)；

f5、计算轧制压力y_j点的工作辊直径D_w(y_j)；

f6、判断j＝num是否成立，若成立，则进入f7；若不成立，令j＝j+1，返回f2；

f7、令j＝1；

f8、计算y_j处负载辊缝形状：

s(y_j)＝2f_w(L_s/2)-[f_w(y_j)+f_w(L_s-y_j)]+2δ_ws(L_s/2)-

；

[δ_ws(y_j)+δ_ws(L_s-y_j)]+[D_w(y)+D_w(L_s-y_j)]/2-D_w(L_s/2)

f9、判断j＝num是否成立，若成立，计算结束；若不成立，令j＝j+1，返回f8。