CN114672753A

CN114672753A - 一种热镀锌沉没辊系转动失效预报方法

Info

Publication number: CN114672753A
Application number: CN202210357913.3A
Authority: CN
Inventors: 白振华; 钱胜; 王孝建; 徐烨明
Original assignee: Yanshan University
Current assignee: Yanshan University
Priority date: 2022-04-06
Filing date: 2022-04-06
Publication date: 2022-06-28
Anticipated expiration: 2042-04-06
Also published as: CN114672753B

Abstract

本发明提供了一种热镀锌沉没辊系转动失效预报方法，涉及热镀锌技术领域，该方法基于条元法对与辊面接触带钢划分单元，基于辊面粗糙度测量值计算辊面摩擦系数，再根据能量平衡方程、达朗贝尔原理建立带钢条元动力学方程求出条元摩擦力及接触压力，并依据打滑临界条件对任意条元判断是否打滑，求出带钢对辊面总滚动摩擦力和总接触压力，将总接触压力与轴端间隙所承受最大承载比较来判断轴端间隙接触关系，依次分别求出带钢对辊面摩擦驱动力矩、刮刀对辊面阻力矩、轴端摩擦阻力矩，依据辊系转动失效临界条件判断实时在线预测辊的运行状态。本发明为后续的沉没辊系各辊工艺参数优化的实现提供数据支撑，从而实现对沉没辊系受力状态精确动态调控。

Description

一种热镀锌沉没辊系转动失效预报方法

技术领域

本发明涉及热镀锌技术领域，特别涉及一种热镀锌沉没辊系转动失效预报方法。

背景技术

热镀锌是板带产品的最后一道工序，板带表面质量，如镀锌厚度均匀性、上下表面镀层厚度偏差量是热镀板带表面质量的主要衡量指标。板带在热镀锌工艺时，首先经过连续退火工艺过程，而后进入锌锅内的沉没辊系(沉没辊、稳定辊及矫正辊构成)，随后出锌锅进入气刀处，刮去板带表面的多余的锌层，最后冷却成形。

由于锌锅内含有很多锌渣，如浮渣和沉渣，同时锌锅是处于高温、锌液场等多物理场作用下沉没辊系各辊受力比较复杂，导致沉没辊系各辊容易出现转动失效，如轴端的轴瓦与轴套间锌渣沉积过量而出现卡死现象，另外，带钢作用辊面的驱动摩擦力小于辊受到阻力矩而出现打滑现象等。为研究实现沉没辊系在锌锅内处于最佳工作状态，需要掌握辊系受力状态，避免辊系发生转动失效，实现对其精确化控制。

目前现场难以实现沉没辊系各辊的转动失效预报，也就无法实现相关工艺参数优化。因此，为了保证热镀锌机组的正常生产，提高带钢产品镀锌质量，减小沉没辊系卡死或打滑发生率，亟需一种沉没辊系转动失效预报新方法。

发明内容

有鉴于此，本发明提供了一种沉没辊系转动失效预报方法，根据锌锅内浸没在锌液内沉没辊系工作状态，实现对沉没辊、稳定辊、矫正辊的转动失效预报，以保证热镀锌机组的正常生产，提高带钢产品镀锌质量。

为此，本发明提供了以下技术方案：

本发明提供了一种热镀锌沉没辊系转动失效预报方法，该方法包括以下步骤：

步骤A、设定沉没辊系的各辊的辊号k的初值k＝1；

步骤B、获取镀锌机组沉没辊系各辊的设备参数以及带钢生产工艺参数；

步骤C、对包角范围内接触带钢划分单元，单元的编号初值i＝1，j＝1，打滑单元数n_h＝0，单元(1，1)的旋入段张力的初值

辊的驱动力矩的初始值M_s＝0；其中，T₀为带钢生产工艺参数，表示带钢初始设定张力；

步骤D、利用镀锌机组沉没辊系辊号为k的辊设备参数，根据现场辊面摩擦系数经验模型确定辊面最大静摩擦系数τ_max和平均摩擦系数τ；

步骤E、利用镀锌机组沉没辊系辊号为k的辊设备参数以及带钢生产工艺参数，根据能量守恒原理和达朗贝尔原理确定带钢单元与辊面间接触压力N_i,j、滚动摩擦力f_i,j及总接触压力N分量；

步骤F、判断带钢单元打滑临界条件是否成立：

若带钢单元打滑临界条件成立，则单元(i,j)不发生打滑；若带钢单元打滑临界条件不成立，则单元(i,j)发生打滑，并令n_h＝n_h+1，判断i≤m且j≤n是否成立；若i≤m且j≤n成立，则令i＝i+1，j＝j+1，T_i,j＝T_i,′_j，其中，T_i,′_j为带钢单元旋出端的张力，并返回步骤F；若i≤m且j≤n不成立，则判断辊打滑条件

是否成立；若辊打滑条件成立，则认定辊发生打滑；若辊打滑条件不成立，则认定辊未发生打滑；式中，ξ_n为带钢子单元发生打滑数量的最大占比，根据机组实际工况及生产效率预先设定；

若单元(i,j)不发生打滑，则计算带钢有效驱动力矩M_s和刮刀力对辊的总阻力矩M_z，计算公式如下：M_s＝f_i,jR+M_s；M_z＝M_d+M_g；其中，M_g为刮刀力对辊的阻力矩；M_d为轴端摩擦阻力矩；

步骤G、判断辊发生转动条件δ≥[δ]是否成立，其中，转动动力因数

[δ]为转动动力因数临界值；若辊发生转动条件成立，则辊未发生转动失效；若辊发生转动条件不成立，则辊发生转动失效。

进一步地，沉没辊系包括：沉没辊、稳定辊及矫正辊。

进一步地，所述方法还包括：

步骤I、判断k≤3是否成立；其中，沉没辊、稳定辊及矫正辊的辊号为1、2、3；若k≤3成立，则k＝k+1，并转入步骤B判断辊是否发生转动失效。

进一步地，包角范围内接触带钢划分单元，包括：

将辊接触区带钢沿辊周向划分m个横向条元；再沿任意辊接触区带钢横向条元的横向划分n个子单元；

计算与沉没辊接触带钢的任意(i,j)子单元坐标值(θ_i，z_j)：

其中，L为辊的设备参数，表示辊身长，B、θ、θ_-为带钢生产工艺参数，B表示带材宽度，θ表示带钢与辊包角初始值，θ_-表示包角的负角。

进一步地，确定辊面最大静摩擦系数τ_max和平均摩擦系数τ，包括：

在辊面沿辊身方向在带钢接触辊的部分均匀选取m_s个点，并测量这些点所在区域锌渣厚度h_k，其中k＝1,2,...,m_s，根据测量得到的h_k计算沉没辊、稳定辊及矫正辊的辊面平均锌渣沉积厚度

计算辊面子单元所在辊的全周向区域摩擦系数τ：

τ₀为带钢生产工艺参数，表示带钢与辊的辊面间的初始摩擦系数；

计算辊面的最大静摩擦系数τ_max：

其中，h_max为h_k中的最大值，

进一步地，确定带钢单元与辊面间接触压力N_i,j、滚动摩擦力f_i,j及总接触压力N分量，包括：

计算带钢单元受到重力与浮力的合力Δg_i,j：

Δg_i,j＝(ρ_Fe-ρ_Zn)g·Δθ·R·Δz·H；其中，

依据能量平衡建立方程：

式中，I为带钢单元的惯性矩，

E为高温下带弹性模量；

沿带钢单元沿辊的径向和切向建立动力学方程：

式中，R为辊的设备参数，表示辊半径；

H、ρ_Fe、ρ_Zn、v为带钢生产工艺参数，分别表示带材的厚度、带钢密度、锌锅内高温下锌液密度、带钢运行速度；

联立上述方程求出接触压力N_i,j、摩擦力f_i,j；

计算总接触压力X向分力N_x、Y向分力N_y：

进一步地，M_g的计算方式如下：

M_g＝τF_gRcosγ；其中，F_g、γ为带钢生产工艺参数，分别表示刮刀力、刮刀力与辊径向的夹角。

进一步地，M_d的计算方式如下：

计算辊的最大承载量[P]、辊受到的外载荷P：

式中，β为轴颈在轴瓦间隙内转动时最小油膜厚度至出口间的夹角；χ表示辊的轴端偏心率；ζ^*表示辊系轴端安全承载系数；ΔG为沉没辊系重力与浮力的合力，ΔG＝G-ρ_ZngV；其中，g为重力加速常量；V、G、η、L_d、d、D为辊的设备参数，分别表示辊的有效体积、辊的重力、锌锅内锌液的粘度、轴端的轴瓦长度、轴颈半径、轴瓦内径；

判断辊的轴端的轴瓦与轴颈是否接触：判断P≤[P]是否成立；若成立，则计算轴端摩擦阻力矩：

若不成立，则转入计算轴端摩擦阻力矩：

式中，

ζ₁为速度影响系数；ζ₂为镀锌量影响系数，Q_s为带钢生产工艺参数，表示带钢表面单位面积上镀锌量，μ₀为辊的设备参数，表示辊的轴套与轴瓦间初始摩擦系数。

本发明的优点和积极效果：本发明结合热镀锌机组的设备特点，充分考虑沉没辊系特性及各辊结构及受力特点，首先基于条元法对与辊面接触带钢划分单元，随后基于辊面粗糙度测量值计算辊面摩擦系数，再根据能量平衡方程、达朗贝尔原理建立带钢条元动力学方程求出条元摩擦力及接触压力，并依据打滑临界条件对任意条元判断是否打滑，从而求出带钢对辊面总滚动摩擦力和总接触压力，将总接触压力与轴端间隙所承受最大承载比较来判断轴端间隙接触关系，依次分别求出带钢对辊面摩擦驱动力矩、刮刀对辊面阻力矩、轴端摩擦阻力矩，最后依据辊系转动失效临界条件判断实时在线预测辊的运行状态。建立热镀锌机组沉没辊系工作过程中转动失效预报模型(本发明中热度锌沉没辊系转动失效预报方法的整个流程)，实时在各种工况下预测沉没辊系转动失效预测，并为辊系发生卡死或打滑等现象提供预警，为后续的沉没辊系各辊工艺参数优化的实现提供数据支撑，从而最终实现对沉没辊系受力状态精确动态调控。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例中沉没辊系组成示意图；

图2(a)为本发明实施例中对与辊接触带钢划分横向条元；

图2(b)为本发明实施例中对带钢横向条元划分子单元；

图3为本发明实施例中沉没辊系各辊转动失效预报流程图。

具体实施方式

为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

以如图1所示的某镀锌机组沉没辊系各辊(包括沉没辊、稳定辊及矫正辊)为例，结合图3示出的沉没辊系各辊转动失效预报流程图，对本发明提供的一种热镀锌沉没辊系转动失效预报方法进行详细说明。

为计算及表述方便，以沉没辊系某辊辊身长度中点处为坐标原点，以过辊身中点的截面为XOY坐标平面，以辊轴线方向为Z轴方向，建立如图2所示的空间直角坐标系O-XYZ。同时在沉没辊系的端面上建立角坐标系，以沉没辊的Y轴负向为角坐标起始边，顺时针为正角区，逆时针方向为负角区，如图2(a)所示。并以围绕在沉没辊系表面上包角范围内带钢为对象，将其分别沿带钢的纵向(带钢运行方向)划分横向条元，再在任意横向条元上沿横向(辊系轴向)划分子单元，如图2(b)所示。

实施例1：

本发明实施例中提供的一种热镀锌沉没辊系转动失效预报方法包括以下步骤：

步骤A1、设定沉没辊系各辊的辊号k的初值k＝1，对沉没辊取k＝1，稳定辊k＝2，矫正辊k＝3；

步骤B1、收集镀锌机组沉没辊系设定辊的设备参数，包括辊半径R＝0.3m，辊的重力G＝9.7kN，辊的有效体积V＝0.124m³，辊的辊身长为L＝1.6m，轴颈半径d＝80mm，轴端的轴瓦长度L_d＝175mm，轴瓦内径D＝110mm，辊的轴套与轴瓦间初始摩擦系数μ₀＝0.01，锌锅内锌液的粘度η＝2.75pa.s。

步骤C1、收集带钢生产工艺参数，包括带钢运行速度v＝120m/min，带钢初始设定张力T₀＝15kN，带材的宽度B＝1000mm及厚度H＝1.2mm，带钢密度ρ_Fe＝7850kg/m3，带钢与辊包角初始值θ＝130^°，包角的负角θ_-＝87^°，锌锅内高温下锌液密度ρ_Zn＝6510kg/m³，带钢与辊的辊面间的初始摩擦系数τ₀＝0.02，刮刀力的数值F_g＝875N，刮刀力与辊面径向的角度γ＝8°，带钢表面单位面积上镀锌量Q_s＝45g/m²。

步骤D1、计算沉没辊系重力与浮力的合力ΔG：ΔG＝G-ρ_ZngV＝994.7N；

式中，g为重力加速常量，取9.8N/kg。

步骤E1、包角范围内接触带钢划分单元；

将辊接触区带钢沿辊周向划分横向(辊轴向)条元，其数目分别为m＝10；再沿任意辊接触区带钢横向条元的横向划分n＝20个子单元：

计算与辊接触带钢的任意(i,j)子单元坐标值：

表1

(i,j)	(1,1)	(1,2)	(2,3)	(2,4)	(3,5)	(3,6)	(4,7)	(4,9)	(5,1)
										θ<sub>i</sub>/°	-93.5	-83	-76	-62.5	-51.5	-34	-20.5	16	22.5
z<sub>i</sub>/mm	300	350	400	450	500	550	600	700	300

步骤F1、确定辊面最大静摩擦系数τ_max、平均摩擦系数τ；

在辊面沿辊身方向在带钢接触辊的部分均匀选取m_s＝20点，并通过测量传感器测量这些点所在区域锌渣厚度h_k，其中k＝1,2,...,n^*；n^*≤n；

运行过程线上测量，测量传感器沿辊轴向均匀移动，等间隔在辊面采样出若干点，每个采样点对应带钢子单元；

计算沉没辊、稳定辊及矫正辊的辊面平均锌渣沉积厚度

计算辊面子单元所在辊的全周向区域摩擦系数τ：

计算辊面的最大静摩擦系数τ_max：

步骤G1、设定带钢相关参数的初值，单元的编号初值i＝1，j＝1，打滑单元数n_h＝0，单元(1，1)的旋入段张力的初值

辊的驱动力矩的初始值M_s＝0。

步骤H1、计算带钢单元(1,1)与辊面间接触压力N_1,1、滚动摩擦力f_1,1及总接触压力N分量；具体地：

计算带钢单元受到重力与浮力的合力Δg_1,1：

Δg_1,1＝(ρ_Fe-ρ_Zn)g·Δθ·R·Δz·H＝0.753kg；

依据能量平衡建立方程：

式中，I为带钢单元的惯性矩，

根据欧拉公式T_i,′_j＝T_i,je^τΔθ，E为高温下带弹性模量。

带钢单元沿辊的径向和切向建立动力学方程：

联立上述方程，求出接触压力N_1,1＝28.5N、摩擦力f_1,1＝27.9N。

计算总接触压力X向分力N_x、Y向分力N_y：

步骤I1、判断带钢单元打滑临界条件是否成立：

若条件成立，则单元(i,j)不发生打滑，并转入步骤J1；若条件不成立，则单元(i,j)发生打滑，并令n_h＝n_h+1且转入步骤K1。

步骤J1、计算带钢有效驱动力矩M_s＝f_i,jR+M_s＝15.3N.m。

步骤K1、判断i≤m且j≤n是否成立；若条件成立，则令i＝i+1，j＝j+1，令T_i,j＝T_i,′_j并转入步骤H1；若条件不成立，则转入步骤L1。

步骤L1、判断辊是否发生打滑：判断条件

是否成立；若条件成立，则认定辊发生打滑，并转入步骤P1；若条件不成立，则认定辊未发生打滑；

式中，ξ_n为带钢子单元发生打滑数量的最大占比。

步骤M1、计算总接触压力N及载荷方向角φ：

步骤N1、计算轴端摩擦阻力矩M_d；

具体地：

计算辊的最大承载量[P]、辊受到的外载荷P；

式中，β—轴颈在轴瓦间隙内转动时最小油膜厚度至出口间的夹角；

判断辊的轴端的轴瓦与轴颈是否接触：判断P≤[P]是否成立；若成立，则计算轴端摩擦阻力：

若不成立，则计算轴端摩擦阻力：

式中，ζ₁为速度影响系数；ζ₂为镀锌量影响系数。

步骤O1、判断辊转动是否失效；

计算刮刀力对辊的阻力矩M_g、总阻力矩M_z：

判断条件δ≥[δ]是否成立；若条件成立，则辊未发生转动失效；若条件不成立，则辊发生转动失效。

步骤P1、判断k≤3是否成立；若条件成立，则k＝k+1，并转入步骤B1。

对于稳定辊及矫正辊，继续对k循环，依次按照上面顺序计算和判断。

步骤Q1、输出辊转动的判断结果，程序结束。

最后输出镀锌机组沉没辊系的摩擦驱动力矩、轴端阻力矩、刮刀力对辊的阻力矩、转动动力因数，如表2所示。

表2

从表2可以看出沉没辊、矫正辊转动动力因数大于其转动动力因数临界值，故沉没辊、矫正辊未出现转动失效；而稳定辊转动动力因数小于其转动动力因数临界值，稳定辊出现转动失效。通过现场观测发现与实际工况相符。

实施例2：

步骤A2、设定沉没辊系各辊的辊号k的初值k＝1，对沉没辊取k＝1，稳定辊k＝2，矫正辊k＝3；

步骤B2、收集镀锌机组沉没辊系设定辊的设备参数，包括辊半径R＝0.4m，辊的重力G＝11.7kN，辊的有效体积V＝0.157m³，辊的辊身长为L＝1.6m，轴颈半径d＝90mm，轴端的轴瓦长度L_d＝175mm，轴瓦内径D＝117mm，辊的轴套与轴瓦间初始摩擦系数μ₀＝0.02，锌锅内锌液的粘度η＝2.75pa.s。

步骤C2、收集带钢生产工艺参数，包括带钢运行速度v＝125m/min，带钢初始设定张力T₀＝18kN，带材的宽度B＝1200mm及厚度H＝0.5mm，带钢密度ρ_Fe＝7850kg/m3，带钢与辊包角初始值θ＝135^°，包角的负角θ_-＝90^°，锌锅内高温下锌液密度ρ_Zn＝6510kg/m³，带钢与辊的辊面间的初始摩擦系数τ₀＝0.03，刮刀力的数值F_g＝950N，刮刀力与辊面径向的角度γ＝7.1°，带钢表面单位面积上镀锌量Q_s＝40g/m²。

步骤D2、计算沉没辊系重力与浮力的合力ΔG：ΔG＝G-ρ_ZngV＝1683N；

式中，g为重力加速常量，取9.8N/kg。

步骤E2、包角范围内接触带钢划分单元；

将辊接触区带钢沿辊周向划分横向(辊轴向)条元，其数目分别为m＝15；再沿任意辊接触区带钢横向条元的横向划分n＝20个子单元：

计算与辊接触带钢的任意(i,j)子单元坐标值：

表3

(i,j)	(1,1)	(1,2)	(2,3)	(2,4)	(3,5)	(3,6)	(4,7)	(4,9)	(5,1)
										θ<sub>i</sub>/°	-90.5	-82	-70	-58.5	-46.5	-30	-18.5	10	20.5
z<sub>i</sub>/mm	350	400	450	500	550	600	650	700	300

步骤F2、确定辊面最大静摩擦系数τ_max、平均摩擦系数τ；

在辊面沿辊身方向在带钢接触辊的部分均匀选取m_s＝20点，并测量这些点所在区域锌渣厚度，计算沉没辊、稳定辊及矫正辊的辊面平均锌渣沉积厚度

计算辊面子单元所在辊的全周向区域摩擦系数τ：

计算辊面的最大静摩擦系数τ_max：

步骤G2、设定带钢相关参数的初值，单元的编号初值i＝1，j＝1，打滑单元数n_h＝0，单元(1，1)的旋入段张力的初值

辊的驱动力矩的初始值M_s＝0。

步骤H2、计算带钢单元(1,1)与辊面间接触压力N_1,1、滚动摩擦力f_1,1及总接触压力N分量；具体地：

计算带钢单元受到重力与浮力的合力Δg_1,1：

Δg_1,1＝(ρ_Fe-ρ_Zn)g·Δθ·R·Δz·H＝0.689kg；

依据能量平衡建立方程：

式中，I为带钢单元的惯性矩，

E为高温下带弹性模量。

带钢单元沿辊的径向和切向建立动力学方程：

计算总接触压力X向分力N_x、Y向分力N_y：

步骤I2、判断带钢单元打滑临界条件是否成立：

步骤J1、计算带钢有效驱动力矩M_s＝f_i,jR+M_s＝17.2N.m。

步骤L2、判断辊是否发生打滑：判断条件

式中，ξ_n为带钢子单元发生打滑数量的最大占比。

步骤M2、计算总接触压力N及载荷方向角

步骤N2、计算轴端摩擦阻力矩M_d；

具体地：

计算辊的最大承载量[P]、辊受到的外载荷P；

若不成立，则计算轴端摩擦阻力：

式中，ζ₁为速度影响系数；ζ₂为镀锌量影响系数。

步骤O2、判断辊转动是否失效；

计算刮刀力对辊的阻力矩M_g、总阻力矩M_z：

步骤P2、判断k≤3是否成立；若条件成立，则k＝k+1，并转入步骤B2。

步骤Q2、输出辊转动的判断结果，程序结束。

最后输出镀锌机组沉没辊系的摩擦驱动力矩、轴端阻力矩、刮刀力对辊的阻力矩、转动动力因数，如表4所示。

表4

从表4可以看出沉没辊、矫正辊、稳定辊转动动力因数大于其转动动力因数临界值，故沉没辊、矫正辊、稳定辊未出现转动失效。通过现场观测发现与实际工况相符。

实施例1、2分别针对两种规格带钢，实施例1针对生产厚度大于1mm带钢的条件下沉没辊系转动失效的实时监测，实施例2针对厚度小于1mm带钢的条件下沉没辊系转动失效的实时监测，这两种实施例1、2具有代表性，可作为机组生产相关规格带钢的参考。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims

1.一种热镀锌沉没辊系转动失效预报方法，其特征在于，所述方法包括以下步骤：

步骤A、设定沉没辊系的各辊的辊号k，初值k＝1；

步骤C、将辊接触区带钢沿辊周向划分m个横向条元；再沿任意辊接触区带钢横向条元的横向划分n个子单元，并计算与沉没辊接触带钢的任意(i,j)子单元坐标值(θ_i，z_j)；单元的编号初值i＝1，j＝1，打滑单元数n_h＝0，单元(1，1)的旋入段张力的初值

步骤F、判断带钢单元打滑临界条件是否成立：

若带钢单元打滑临界条件成立，则单元(i,j)不发生打滑；若带钢单元打滑临界条件不成立，则单元(i,j)发生打滑，并令n_h＝n_h+1，判断i≤m且j≤n是否成立；若i≤m且j≤n成立，则令i＝i+1，j＝j+1，T_i,j＝T′_i,j，其中，T′_i,j为带钢单元旋出端的张力，并返回步骤F；若i≤m且j≤n不成立，则判断辊打滑条件