本发明属于冶金
技术领域：
，特别涉及一种带钢平整生产工艺方法。
背景技术：
：近年来，随着汽车、航天、家电、建筑、包装等行业的迅猛发展，高等级高品质冷轧板带产品的需求也不断扩大，六辊轧机因其刚度大、板形控制能力强而越来越得到钢铁企业的青睐。特别是部分中小型钢铁企业为了提高轧机的板形控制能力与满足轧薄的需要，新建了许多工作辊“长径比(即工作辊辊身长与直径之比)”较大且直径绝对值较小的小型六辊轧机。与正常工作辊辊径的轧机不一样，“长径比”较大且工作辊辊径绝对值较小的轧机在轧制过程中工作辊除了产生垂直方向的挠曲变形之外，水平方向的挠曲也不可忽略。与此同时，由于工作辊水平挠曲的存在，使得工作辊与带材之间、工作辊与中间辊之间的接触状态以及受力状态发生了较大的变化，导致出口板形预报不准确，板形精度达不到要求，影响成品带材质量。但是纵观国内外相关文献[1‐5]，关于六辊轧机的板形问题的研究成果都是建立在常规的板形控制基础上，研究的焦点集中于轧辊的垂直挠曲，都没有考虑到工作辊的水平挠曲问题，这样，如何定量分析水平挠曲对板形的影响，并提出相应的板形控制措施依然是现场技术攻关的焦点。(参考文献：[1]连家创，刘宏民.板厚板形控制[M].兵器工业出版社，1995.[2]刘宏民，郑振忠，彭艳.六辊CVC宽带轧机板形控制特性的计算机模拟[J].钢铁研究学报.2001,13(1):15～18.[3]刘宏民.三维轧制理论及其应用.北京:科学出版社,1999.[4]杜凤山，薛涛，孙静娜.六辊CVC轧机轧辊弯曲和压扁变形的有限元分析[J].燕山大学学报.2011,35(5):396～401.[5]白振华.平整轧制工艺模型[M].北京：冶金工业出版社，2010.)技术实现要素：本发明的目的是针对长径比大且直径绝对值小的六辊轧机因工作辊产生水平挠曲，而导致成品带材板形预报不准确的问题，提供了一种工作辊长径比大且直径绝对值小的六辊轧机板形预报方法。该方法建立了一套适合长径比大的六辊轧机板形模型，能够定量预报工作辊发生水平挠曲时出口板形值的分布，为成品带材的板形控制提供了依据。本发明包括以可由计算机执行的步骤：(a)收集待预报的六辊轧机的主要设备参数，主要包括：工作辊弯辊力Fw；中间辊弯辊力Fm；中间辊窜辊量δc；支撑辊压下螺丝与轧制中心线的距离lbl、lbr；工作辊弯辊缸与轧制中心线的距离lwl、lwr；中间辊弯辊缸与轧制中心线的距离lml、lmr；工作辊、中间辊、支撑辊的辊身长度Lw、Lm、Lb；工作辊、中间辊、支撑辊的直径Dw、Dm、Db；工作辊、中间辊、支撑辊的辊型Dwi、Dmi、Dbi；带材跑偏量为δp；倾辊量为η；(b)收集待轧制带材的特征参数和轧制工艺参数，主要包括：带材的变形抗力σs；带材的来料宽度B；带材来料厚度平均值h0；带材来料厚度横向分布值h0i；带材的弹性模量E；带材的泊松比v；轧制速度V；压下率ε；前、后张力平均值T0、T1；(c)定义预报过程中所涉及的过程变量，主要包括：垂直方向上工作辊相对支撑辊的刚性转角为βvw；中间辊相对支撑辊的刚性转角βm；带材的出口板形分布shapei；前、后张力横向分布σ1i、σ0i；带材出口厚度横向分布h1i,h1i′；支撑辊沿辊身分段数N；支撑辊各段宽度Δx；带材沿宽度方向分段数M；上、下支撑辊单元划分过程参数n；带材单元划分过程参数m；过程变量i,j；带材跑偏所占单元数np；中间辊窜动引起的支承辊与中间辊、工作辊与中间辊辊间压力分布区间变化单元数nc1、nc2；支撑辊j段载荷引起i段挠度的影响系数Gbij；支撑辊的支撑力对支撑辊i段挠度的影响系数GFbi；中间辊j段载荷引起i段挠度的影响系数Gmij；中间辊的弯辊力对中间辊i段挠度的影响系数GFmi；工作辊j段载荷引起i段挠度的影响系数Gwij；工作辊的弯辊力对工作辊i段挠度的影响系数GFwi；工作辊的水平支撑力对工作辊i段挠度的影响系数GFlwi；工作辊与带材接触在第j单元内垂直方向的接触压力qvj；工作辊与带材接触在第j单元内水平方向的接触压力qlj；中间辊、支撑辊辊间压力分布qmbj；中间辊、工作辊辊间压力分布qmwj；工作辊水平方向挠度分布fwli；工作辊垂直方向挠度分布fwvi；工作辊横向凸度值ΔDwi；中间辊横向凸度值ΔDmi；支撑辊横向凸度值ΔDbi；过程变量Bli,Bri,Wvli,Wvri；中间辊和工作辊在第i单元内辊凸度增量Δi；工作辊在第i单元内偏移角αi；工作辊左右两侧的水平支撑力Flwz、Flwy；(d)将轧辊和带材进行单元划分并求解影响系数，主要包括以下步骤：d1)将支撑辊沿辊身长度方向划分为N等分，并计算出支撑辊各段宽度d2)计算待轧制带材沿宽度方向分段数M，并令d3)计算上、下支撑辊单元划分过程参数n；带材单元划分过程参数m，并令d4)计算带材跑偏所占单元数np、上下中间辊窜动引起辊系压力分布区间变化单元数nc1、nc2，并令d5)计算工作辊挠度影响系数Gwij、工作辊弯辊力影响系数GFwi、工作辊水平支撑力影响系数GFlwi、中间辊挠度影响系数Gmij、中间辊弯辊力影响系数GFmi、支撑辊挠度影响系数Gbij、支撑辊支撑力影响系数GFbi；(e)预报六辊轧机工作辊水平挠曲时带材前张力横向分布值，主要包括以下步骤：e1)给定带材出口厚度横向分布初始值h1i′；e2)将给定的厚度分布值带入到金属变形模型中计算前张力横向分布值σ1i、后张力横向分布值σ0i；e3)根据前后张力值计算工作辊与带材在各单元内水平接触压力1≤i＝j≤2n+1；e4)根据轧制压力模型，计算出当前带材出口厚度横向分布、前后张力分布下的轧制压力横向分布qvi；e5)根据工作辊水平方向上力和力矩平衡方程，计算出工作辊水平支撑力Flwz、Flwy，如下所示：其中xi为第i单元到轧制中心线的距离；e6)由材料力学的知识可得到工作辊水平方向上挠度可由下式表示：e7)根据工作辊水平方向上的挠度，得到工作辊在各单元内的偏移角e8)根据中间辊与工作辊的受力图，可以得到中间辊和工作辊在各单元内产生的附加凸度值e9)引入过程参数ρ1i、ρ2i、bmwzi、bmwyi、kmwi，其中ρ1i、ρ2i分别表示工作辊水平挠曲后在轧制中线两侧曲率半径，bmwzi、bmwyi分别表示工作辊第i段挠曲后左右部分距轧制中线的宽度，为中间辊与工作辊第i段单位辊身长度上的触压力，kmwi表示中间辊和工作辊的压扁系数；e10)考虑到工作辊的水平挠曲，中间辊与工作辊不对称接触，建立工作辊水平挠曲时压扁系数模型，计算出辊间压扁宽度，如下所示：式中，E1、ν1分别表示工作辊的弹性模量和泊松比，E2、ν2分别表示中间辊的弹性模量和泊松比；e11)根据辊间压扁宽度，计算出工作辊水平挠曲时各单元内辊间压扁系数，如下所示Kmwi=2[(1-v12)πE1(1-fwli2Rw2)(bmwzibmwzi+bmwyiln2Rwbmwzi+bmwyibmwzi+bmwyiln2Rwbmwyi+0.407)+(1-v22)πE2(1-fwli2Rm2)(bmwzibmwzi+bmwyiln2Rmbmwzi+bmwyibmwzi+bmwyiln2Rmbmwyi+0.036)]·cosai]]>e12)根据辊系的变形协调关系，可以给出工作辊弯辊力Fw，带材跑偏量为δp，中间辊窜辊量为δc；中间辊弯辊力Fm，支撑辊支撑力Fbl、Fbr，支撑辊倾辊量η，工作辊辊型Dwi，中间辊辊型Dmi，支撑辊辊型Dbi等相关工艺与设备参数间的关系：Σj=1nqmwj(Gwij+Gmij)+Kmwiqmwi-Kmwiqmw(n+1)-Σj=1nqmbjGmij-Σj=1nqvjGwij-(βvw-βm)x(i)=Wvli,1≤i≤n]]>Σj=n+22n+1qmwj(Gwij+Gmij)+Kmwiqmwi-Kmwiqmw(n+1)-Σj=n+22n+1qmbjGmij-Σn+22n+1qvjGwij-(βvw-βm)x(i)=Wvri,n+2≤i≤2n+1]]>Σj=1nqmbj(Gmij+Gbij)+Kmbiqmbi-Kmbiqmb(n+1)-Σj=1nqmwjGmij-βmx(i)=Bli,1≤i≤n]]>Σj=n+22n+1qmbj(Gmij+Gbij)+Kmbiqmbi-Kmbiqmb(n+1)-Σj=n+22n+1qmwjGmij-βmx(i)=Bri,n+2≤i≤2n+1]]>其中：Wvli=FwGFwi-FmGFmi-ΔDwi+ΔDmi2-Δi,1≤i≤n]]>Wvri=FwGFwi-FmGFmi-ΔDwi+ΔDmi2-Δi,n+2≤i≤2n+1]]>Bli=FmGFmi+FblGFbi-ΔDmi+ΔDbi2+ηLbξxi,1≤i≤n]]>Bri=FmGFmi+FbrGFbi-ΔDmi+ΔDbi2+ηLbξxi,n+2≤i≤2n+1]]>式中：ξ为考虑轧辊弹性变形时支撑辊倾辊量影响系数；Kmbi为支承辊与中间辊在第i单元内辊间压扁系数，Kmwi为工作辊与中间辊在第i单元内辊间压扁系数；e13)根据支撑辊的受力以及力矩平衡，给出相应的平衡方程，如下所示：Σi=12n+1qmbi=Fbr+Fbl]]>Σi=12n+1qmbixi=Fbrlbr+Fbllbl]]>e14)根据中间辊的受力以及力矩平衡，给出相应的平衡方程，如下所示：Σi=12n+1qmwi=Fbl+Fbr-2Fm]]>Σi=12n+1qmwixi=Fbllbl-Fbrlbr+Fmlml-Fmlmr]]>e15)根据e12)-e14)即可求得支撑辊、中间辊辊间压力分布qmbi；工作辊、中间辊辊间压力分布qmwi；垂直方向上工作辊相对支撑辊的刚性转角为βvw；中间辊相对支撑辊的刚性转角βm；e16)根据中间辊、工作辊辊间压力分布qmwi；工作辊相对支撑辊的刚性转角为βvw；计算出垂直方向上工作辊的挠度分布，计算模型如下所示：fwvi=Σj=1n(qmwj-qvj)Gwij-FwGFwi+βvwxi1≤i≤nfwvi=Σj=n+22n+1(qmwj-qvj)Gwij-FwGFwi+βvwxin+2≤i≤2n+1]]>e17)根据垂直方向上工作辊的挠度分布fwvi计算出口厚度横向分布h1i；e18)判断不等式是否成立？如果不等式成立，转入步骤e19)；如果不等式不成立，则令h1i′＝h1i，转入步骤e2)重新计算；e19)输出六辊轧机轧制时工作辊发生水平挠曲现象带材前张力横向分布值σ1i；(f)根据带材的前张力横向分布值预报出六辊轧机工作辊水平挠曲时板形值分布本发明与现有技术相比具有如下优点：可以定量预报出工作辊产生水平挠曲时对轧机成品板形的影响，为成品带材的板形控制提供了依据。附图说明图1是本发明的总计算流程图；图2是本发明将轧辊和带材进行单元划分并求解影响系数计算流程图；图3是本发明六辊轧机工作辊水平挠曲时带材前张力横向分布预报计算流程图；图4是本发明实施例1中工作辊辊型曲线图；图5是本发明实施例1中中间辊辊型曲线图；图6是本发明实施例1中支撑辊辊型曲线图；图7是本发明实施例1中来料厚度横向分布图；图8是本发明实施例1中轧制力横向分布图；图9是本发明实施例1中工作辊水平方向挠曲分布图；图10是本发明实施例1中支撑辊与中间辊辊间压力横向分布图；图11是本发明实施例1中工作辊与中间辊辊间压力横向分布图；图12是本发明实施例1中工作辊的挠度曲分布图；图13是本发明实施例1中根据挠度曲线计算出带材出口厚度横向分布图；图14是本发明实施例1中前张力横向分布图；图15是本发明实施例1中板形值横向分布图；图16是本发明实施例2中工作辊辊型曲线示意图；图17是本发明实施例2中中间辊辊型曲线示意图；图18是本发明实施例2中支撑辊辊型曲线示意图；图19是本发明实施例2中来料厚度横向分布图；图20是本发明实施例2中轧制力横向分布图；图21是本发明实施例2中工作辊水平方向挠曲分布图；图22是本发明实施例2中支撑辊与中间辊辊间压力横向分布图；图23是本发明实施例2中工作辊与中间辊辊间压力横向分布图；图24是本发明实施例2中工作辊的挠度曲分布图；图25是本发明实施例2中根据挠度曲线计算出带材出口厚度横向分布图；图26是本发明实施例2中前张力横向分布图；图27是本发明实施例2中板形值横向分布图。具体实施方式实施例1工作辊长径比大且直径绝对值小的六辊轧机板形预报方法，其计算流程如图1所示：首先，在步骤1中，收集六辊轧机的主要设备参数，主要包括：工作辊弯辊力Fw＝49.5t；中间辊弯辊力Fm＝88.6t；中间辊窜辊量δc＝10.6mm；支撑辊压下螺丝与轧制中心线的距离lbl＝760mm、lbr＝760mm；工作辊弯辊缸与轧制中心线的距离lwl＝760mm、lwr＝760mm；中间辊弯辊缸与轧制中心线的距离lml＝760mm、lmr＝760mm；工作辊、中间辊、支撑辊的辊身长度Lw＝450mm、Lm＝490mm、Lb＝450mm；工作辊、中间辊、支撑辊的直径Dw＝80mm、Dm＝175mm、Db＝350mm；工作辊、中间辊、支撑辊的辊型Dwi、Dmi、Dbi分布曲线示意图如图4、图5和图6所示；带材跑偏量为δp＝0mm；倾辊量为η＝0.01mm；随后，在步骤2中，收集待轧制带材的特征参数和轧制工艺参数，主要包括：带材的变形抗力σs＝806MPa；带材的来料宽度B＝310mm；带材来料厚度平均值h0＝0.2832mm；带材来料厚度横向分布值h0i，如图7所示；带材的弹性模量E＝210GPa；带材的泊松比v＝0.3；轧制速度V＝1360m/min；压下率ε＝0.290；前、后张力平均值T0＝120MPa、T1＝60MPa；随后，在步骤3中，定义预报过程中所涉及的过程变量，主要包括：垂直方向上工作辊相对支撑辊的刚性转角为βvw；中间辊相对支撑辊的刚性转角βm；带材的出口板形分布shapei；前、后张力横向分布σ1i、σ0i；带材出口厚度横向分布h1i,h1i′；支撑辊沿辊身分段数N；支撑辊各段宽度Δx；带材沿宽度方向分段数M；上、下支撑辊单元划分过程参数n；带材单元划分过程参数m；过程变量i,j；带材跑偏所占单元数np；中间辊窜动引起的支承辊与中间辊、工作辊与中间辊辊间压力分布区间变化单元数nc1、nc2；支撑辊j段载荷引起i段挠度的影响系数Gbij；支撑辊的支撑力对支撑辊i段挠度的影响系数GFbi；中间辊j段载荷引起i段挠度的影响系数Gmij；中间辊的弯辊力对中间辊i段挠度的影响系数GFmi；工作辊j段载荷引起i段挠度的影响系数Gwij；工作辊的弯辊力对工作辊i段挠度的影响系数GFwi；工作辊的水平支撑力对工作辊i段挠度的影响系数GFlwi；工作辊与带材接触在第j单元内垂直方向的接触压力qvj；工作辊与带材接触在第j单元内水平方向的接触压力qlj；中间辊、支撑辊辊间压力分布qmbj；中间辊、工作辊辊间压力分布qmwj；工作辊水平方向挠度分布fwli；工作辊垂直方向挠度分布fwvi；工作辊横向凸度值ΔDwi；中间辊横向凸度值ΔDmi；支撑辊横向凸度值ΔDbi；过程变量Bli,Bri,Wvli,Wvri；中间辊和工作辊在第i单元内辊凸度增量Δi；工作辊在第i单元内偏移角αi；工作辊左右两侧的水平支撑力Flwz、Flwy；随后，如图2所示，在步骤4中，将轧辊和带材进行单元划分并求解影响系数，包括以下计算步骤：随后，在步骤4-1中，将支撑辊沿辊身长度方向划分为N＝31等分，并计算出支撑辊各段宽度随后，在步骤4-2中，计算待轧制带材沿宽度方向分段数随后，在步骤4-3中，计算上、下支撑辊单元划分过程参数带材单元划分过程参数随后，在步骤4-4中，计算带材跑偏所占单元数辊系辊间压力分布区间变化单元数随后，在步骤4‐5中，计算工作辊挠度影响系数Gwij、工作辊弯辊力影响系数GFwi、工作辊水平支撑力影响系数GFlwi、中间辊挠度影响系数Gmij、中间辊弯辊力影响系数GFmi、支撑辊挠度影响系数Gbij、支撑辊支撑力影响系数GFbi；随后，如图3所示，在步骤5中，预报六辊轧机工作辊水平挠曲时带材前张力横向分布值，包括以下计算步骤：随后，在步骤5‐1中，给定带材出口厚度横向分布初始值h1i′；随后，在步骤5‐2中，将给定的厚度分布值带入到金属变形模型中计算前张力横向分布值σ1i、后张力横向分布值σ0i；随后，在步骤5-3中，根据前后张力值计算工作辊与带材在各单元内水平接触压力1≤i＝j≤31；随后，在步骤5-4中，根据轧制压力模型，计算出当前带材出口厚度横向分布、前后张力分布下的轧制压力横向分布qvi，其分布曲线如图8所示；随后，在步骤5-5中，根据工作辊水平方向上力和力矩平衡方程，计算出工作辊水平支撑力Flwz、Flwy，如下所示：其中xi为第i单元到轧制中心线的距离；随后，在步骤5-6中，由材料力学的知识可得到工作辊水平方向上挠度可由下式表示：随后，在步骤5-7中，根据工作辊水平方向上的挠度，分布曲线如图9所示，得到工作辊在各单元内的偏移角随后，在步骤5-8中，根据中间辊与工作辊的受力情况，可以得到中间辊和工作辊在各单元内产生的附加凸度值随后，在步骤5-9中，引入过程参数ρ1i、ρ2i、bmwzi、bmwyi、kmwi，其中ρ1i、ρ2i分别表示工作辊水平挠曲后在轧制中线两侧曲率半径，bmwzi、bmwyi分别表示工作辊第i段挠曲后左右部分距轧制中线的宽度，为中间辊与工作辊第i段单位辊身长度上的触压力，kmwi表示中间辊和工作辊的压扁系数；随后，在步骤5-10中，考虑到工作辊的水平挠曲，中间辊与工作辊不对称接触，建立工作辊水平挠曲时压扁系数模型，计算出辊间压扁宽度，如下所示：式中，E1、ν1分别表示工作辊的弹性模量和泊松比，E2、ν2分别表示中间辊的弹性模量和泊松比；随后，在步骤5-11中，根据辊间压扁宽度，计算出工作辊水平挠曲时各单元内辊间压扁系数，如下所示：Kmwi=2[(1-v12)πE1(1-fwli2Rw2)(bmwzibmwzi+bmwyiln2Rwbmwzi+bmwyibmwzi+bmwyiln2Rwbmwyi+0.407)+(1-v22)πE2(1-fwli2Rm2)(bmwzibmwzi+bmwyiln2Rmbmwzi+bmwyibmwzi+bmwyiln2Rmbmwyi+0.036)]·cosai]]>随后，在步骤5-12中，根据辊系的变形协调关系，可以给出工作辊弯辊力Fw，带材跑偏量为δp，中间辊窜辊量为δc；中间辊弯辊力Fm，支撑辊支撑力Fbl、Fbr，支撑辊倾辊量η，工作辊辊型Dwi，中间辊辊型Dmi，支撑辊辊型Dbi等相关工艺与设备参数间的关系：Σj=115qmwj(Gwij+Gmij)+Kmwiqmwi-Kmwiqmw(n+1)-Σj=115qmbjGmij-Σj=115qvjGwij-(βvw-βm)x(i)=Wvli,1≤i≤15]]>Σj=1731qmwj(Gwij+Gmij)+Kmwiqmwi-Kmwiqmw(n+1)-Σj=1731qmbjGmij-Σj=1731qvjGwij-(βvw-βm)x(i)=Wvri,17≤i≤31]]>Σj=115qmbj(Gmij+Gbij)+Kmbiqmbi-Kmbiqmb(n+1)-Σj=115qmwjGmij-βmx(i)=Bli,1≤i≤15]]>Σj=1731qmbj(Gmij+Gbij)+Kmbiqmbi-Kmbiqmb(n+1)-Σj=1731qmwjGmij-βmx(i)=Bri,17≤i≤31]]>其中：Wvli=FwGFwi-FmGFmi-ΔDwi+ΔDmi2-Δi,1≤i≤15]]>Wvri=FwGFwi-FmGFmi-ΔDwi+ΔDmi2-Δi,17≤i≤31]]>Bli=FmGFmi+FblGFbi-ΔDmi+ΔDbi2+ηLbξxi,1≤i≤15]]>Bri=FmGFmi+FbrGFbi-ΔDmi+ΔDbi2+ηLbξxi,17≤i≤31]]>式中：ξ为考虑轧辊弹性变形时支撑辊倾辊量影响系数；Kmbi为支承辊与中间辊在第i单元内辊间压扁系数，Kmwi为工作辊与中间辊在第i单元内辊间压扁系数；随后，在步骤5-13中，根据支撑辊的受力以及力矩平衡，给出相应的平衡方程，如下所示：Σi=131qmbi=Fbr+Fbl]]>Σi=131qmbixi=Fbrlbr+Fbllbl]]>随后，在步骤5-14中，根据中间辊的受力以及力矩平衡，给出相应的平衡方程，如下所示：Σi=131qmwi=Fbl+Fbr-2Fm]]>Σi=131qmwixi=Fbllbl-Fbrlbr+Fmlml-Fmlmr]]>随后，在步骤5-15中，综合步骤5-12、5-13、5-14，即可求得支撑辊、中间辊辊间压力分布qmbi，分布曲线如图10所示；工作辊、中间辊辊间压力分布qmwi，分布曲线如图11所示；垂直方向上工作辊相对支撑辊的刚性转角为βvw；中间辊相对支撑辊的刚性转角βm；随后，在步骤5-16中，根据中间辊、工作辊辊间压力分布qmwi；工作辊相对支撑辊的刚性转角为βvw；计算出垂直方向上工作辊的挠度分布，分布曲线如图12所示，计算模型如下所示：fwvi=Σj=115(qmwj-qvj)Gwij-FwGFwi+βvwxi1≤i≤15fwvi=Σj=1731(qmwj-qvj)Gwij-FwGFwi+βvwxi17≤i≤31]]>随后，在步骤5-17中，根据垂直方向上工作辊的挠度分布fwvi计算出口厚度横向分布h1i，分布曲线如图13所示；随后，在步骤5-18中，判断不等式是否成立？显然不等式不成立，则令h1i′＝h1i，转入步骤5-12重新计算，循环直至不等式成立，转入步骤5-19；随后，在步骤5-19中，输出六辊轧机轧制时工作辊发生水平挠曲现象带材前张力横向分布值σ1i，其分布曲线如图14所示；最后，在步骤6中，根据带材的前张力横向分布值预报出长径比大且工作辊直径绝对值小的六辊轧机轧制时的板形值分布其分布曲线如图15所示。实施例2首先，在步骤1中，收集六辊轧机的主要设备参数，主要包括：工作辊弯辊力Fw＝55.4t；中间辊弯辊力Fm＝78.2t；中间辊窜辊量δc＝8.3mm；支撑辊压下螺丝与轧制中心线的距离lbl＝760mm、lbr＝760mm；工作辊弯辊缸与轧制中心线的距离lwl＝760mm、lwr＝760mm；中间辊弯辊缸与轧制中心线的距离lml＝760mm、lmr＝760mm；工作辊、中间辊、支撑辊的辊身长度Lw＝450mm、Lm＝490mm、Lb＝450mm；工作辊、中间辊、支撑辊的直径Dw＝80mm、Dm＝175mm、Db＝350mm；工作辊、中间辊、支撑辊的辊型Dwi、Dmi、Dbi分布曲线如图16、图17和图18所示；带材跑偏量为δp＝0mm；倾辊量为η＝0.03mm；随后，在步骤2中，收集待轧制带材的特征参数和轧制工艺参数，主要包括：带材的变形抗力σs＝806MPa；带材的来料宽度B＝330mm；带材来料厚度平均值h0＝0.3126mm；带材来料厚度横向分布值h0i，如图19所示；带材的弹性模量E＝210GPa；带材的泊松比v＝0.3；轧制速度V＝1360m/min；压下率ε＝0.201；前、后张力平均值T0＝110MPa、T1＝80MPa；随后，在步骤3中，定义预报过程中所涉及的过程变量，主要包括：垂直方向上工作辊相对支撑辊的刚性转角为βvw；中间辊相对支撑辊的刚性转角βm；带材的出口板形分布shapei；前、后张力横向分布σ1i、σ0i；带材出口厚度横向分布h1i,h1i′；支撑辊沿辊身分段数N；支撑辊各段宽度Δx；带材沿宽度方向分段数M；上、下支撑辊单元划分过程参数n；带材单元划分过程参数m；过程变量i,j；带材跑偏所占单元数np；中间辊窜动引起的支承辊与中间辊、工作辊与中间辊辊间压力分布区间变化单元数nc1、nc2；支撑辊j段载荷引起i段挠度的影响系数Gbij；支撑辊的支撑力对支撑辊i段挠度的影响系数GFbi；中间辊j段载荷引起i段挠度的影响系数Gmij；中间辊的弯辊力对中间辊i段挠度的影响系数GFmi；工作辊j段载荷引起i段挠度的影响系数Gwij；工作辊的弯辊力对工作辊i段挠度的影响系数GFwi；工作辊的水平支撑力对工作辊i段挠度的影响系数GFlwi；工作辊与带材接触在第j单元内垂直方向的接触压力qvj；工作辊与带材接触在第j单元内水平方向的接触压力qlj；中间辊、支撑辊辊间压力分布qmbj；中间辊、工作辊辊间压力分布qmwj；工作辊水平方向挠度分布fwli；工作辊垂直方向挠度分布fwvi；工作辊横向凸度值ΔDwi；中间辊横向凸度值ΔDmi；支撑辊横向凸度值ΔDbi；过程变量Bli,Bri,Wvli,Wvri；中间辊和工作辊在第i单元内辊凸度增量Δi；工作辊在第i单元内偏移角αi；工作辊左右两侧的水平支撑力Flwz、Flwy；随后，在步骤4中，将轧辊和带材进行单元划分并求解影响系数，包括以下计算步骤：随后，在步骤4-1中，将支撑辊沿辊身长度方向划分为N＝31等分，并计算出支撑辊各段宽度随后，在步骤4-2中，计算待轧制带材沿宽度方向分段数随后，在步骤4-3中，计算上、下支撑辊单元划分过程参数带材单元划分过程参数随后，在步骤4-4中，计算带材跑偏所占单元数辊系辊间压力分布区间变化单元数随后，在步骤4‐5中，计算工作辊挠度影响系数Gwij、工作辊弯辊力影响系数GFwi、工作辊水平支撑力影响系数GFlwi、中间辊挠度影响系数Gmij、中间辊弯辊力影响系数GFmi、支撑辊挠度影响系数Gbij、支撑辊支撑力影响系数GFbi；随后，在步骤5中，预报六辊轧机工作辊水平挠曲时带材前张力横向分布值，包括以下计算步骤：随后，在步骤5‐1中，给定带材出口厚度横向分布初始值h1i′；随后，在步骤5‐2中，将给定的厚度分布值带入到金属变形模型中计算前张力横向分布值σ1i、后张力横向分布值σ0i；随后，在步骤5-3中，根据前后张力值计算工作辊与带材在各单元内水平接触压力1≤i＝j≤31；随后，在步骤5-4中，根据轧制压力模型，计算出当前带材出口厚度横向分布、前后张力分布下的轧制压力横向分布qvi，其分布曲线如图20所示；随后，在步骤5-5中，根据工作辊水平方向上力和力矩平衡方程，计算出工作辊水平支撑力Flwz、Flwy，如下所示：其中xi为第i单元到轧制中心线的距离；随后，在步骤5-6，由材料力学的知识可得到工作辊水平方向上挠度可由下式表示：随后，在步骤5-7中，根据工作辊水平方向上的挠度，分布曲线如图21所示，得到工作辊在各单元内的偏移角随后，在步骤5-8中，根据中间辊与工作辊的受力情况，可以得到中间辊和工作辊在各单元内产生的附加凸度值随后，在步骤5-9中，引入过程参数ρ1i、ρ2i、bmwzi、bmwyi、kmwi，其中ρ1i、ρ2i分别表示工作辊水平挠曲后在轧制中线两侧曲率半径，bmwzi、bmwyi分别表示工作辊第i段挠曲后左右部分距轧制中线的宽度，为中间辊与工作辊第i段单位辊身长度上的触压力，kmwi表示中间辊和工作辊的压扁系数；随后，在步骤5-10中，考虑到工作辊的水平挠曲，中间辊与工作辊不对称接触，建立工作辊水平挠曲时压扁系数模型，计算出辊间压扁宽度，如下所示：式中，E1、ν1分别表示工作辊的弹性模量和泊松比，E2、ν2分别表示中间辊的弹性模量和泊松比；随后，在步骤5-11中，根据辊间压扁宽度，计算出工作辊水平挠曲时各单元内辊间压扁系数，如下所示Kmwi=2[(1-v12)πE1(1-fwli2Rw2)(bmwzibmwzi+bmwyiln2Rwbmwzi+bmwyibmwzi+bmwyiln2Rwbmwyi+0.407)+(1-v22)πE2(1-fwli2Rm2)(bmwzibmwzi+bmwyiln2Rmbmwzi+bmwyibmwzi+bmwyiln2Rmbmwyi+0.036)]·cosai]]>随后，在步骤5-12中，根据辊系的变形协调关系，可以给出工作辊弯辊力Fw，带材跑偏量为δp，中间辊窜辊量为δc；中间辊弯辊力Fm，支撑辊支撑力Fbl、Fbr，支撑辊倾辊量η，工作辊辊型Dwi，中间辊辊型Dmi，支撑辊辊型Dbi等相关工艺与设备参数间的关系：Σj=115qmwj(Gwij+Gmij)+Kmwiqmwi-Kmwiqmw(n+1)-Σj=115qmbjGmij-Σj=115qvjGwij-(βvw-βm)x(i)=Wvli,1≤i≤15]]>Σj=1731qmbj(Gmij+Gbij)+Kmbiqmbi-Kmbiqmb(n+1)-Σj=1731qmwjGmij-Σj=1731qvjGwij-(βvw-βm)x(i)=Wvri,17≤i≤31]]>Σj=115qmbj(Gmij+Gbij)+Kmbiqmbi-Kmbiqmb(n+1)-Σj=115qmwjGmij-βmx(i)=Bli,1≤i≤15]]>Σj=1731qmbj(Gmij+Gbij)+Kmbiqmbi-Kmbiqmb(n+1)-Σj=1731qmwjGmij-βmx(i)=Bri,17≤i≤31]]>其中：Wvli=FwGFwi-FmGFmi-ΔDwi+ΔDmi2-Δi,1≤i≤15]]>Wvri=FwGFwi-FmGFmi-ΔDwi+ΔDmi2-Δi,17≤i≤31]]>Bli=FmGFmi+FblGFbi-ΔDmi+ΔDbi2+ηLbξxi,1≤i≤15]]>Bri=FmGFmi+FbrGFbi-ΔDmi+ΔDbi2+ηLbξxi,17≤i≤31]]>式中：ξ为考虑轧辊弹性变形时支撑辊倾辊量影响系数；Kmbi为支承辊与中间辊在第i单元内辊间压扁系数，Kmwi为工作辊与中间辊在第i单元内辊间压扁系数；随后，在步骤5-13中，根据支撑辊的受力以及力矩平衡，给出相应的平衡方程，如下所示：Σi=131qmbi=Fbr+Fbl]]>Σi=131qmbixi=Fbrlbr+Fbllbl]]>随后，在步骤5-14中，根据中间辊的受力以及力矩平衡，给出相应的平衡方程，如下所示：Σi=131qmwi=Fbl+Fbr-2Fm]]>Σi=131qmwixi=Fbllbl-Fbrlbr+Fmlml-Fmlmr]]>随后，在步骤5-15中，综合步骤5-12、5-13、5-14，即可求得支撑辊、中间辊辊间压力分布qmbi，分布曲线如图22所示；工作辊、中间辊辊间压力分布qmwi，分布曲线如图23所示；垂直方向上工作辊相对支撑辊的刚性转角为βvw；中间辊相对支撑辊的刚性转角βm；随后，在步骤5-16中，根据中间辊、工作辊辊间压力分布qmwi；工作辊相对支撑辊的刚性转角为βvw；计算出垂直方向上工作辊的挠度分布，分布曲线如图24所示，计算模型如下所示：fwvi=Σj=115(qmwj-qvj)Gwij-FwGFwi+βvwxi1≤i≤15fwvi=Σj=1731(qmwj-qvj)Gwij-FwGFwi+βvwxi17≤i≤31]]>随后，在步骤5-17中，根据垂直方向上工作辊的挠度分布fwvi计算出口厚度横向分布h1i，分布曲线如图25所示；随后，在步骤5-18中，判断不等式是否成立？显然不等式不成立，则令h1i′＝h1i，转入步骤5-2重新计算，循环直至不等式成立，转入步骤5-19；随后，在步骤5-19中，输出六辊轧机轧制时工作辊发生水平挠曲现象带材前张力横向分布值σ1i，其分布曲线附图26所示；最后，在步骤6中，根据带材的前张力横向分布值预报出长径比大且工作辊直径绝对值小的六辊轧机轧制时的板形值分布其分布曲线如图27所示。当前第1页1 2 3