1.一种热镀锌沉没辊系转动失效预报方法，其特征在于，所述方法包括以下步骤：步骤A、设定沉没辊系的各辊的辊号k，初值k＝1；

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

步骤C、将辊接触区带钢沿辊周向划分m个横向条元；再沿任意辊接触区带钢横向条元的横向划分n个子单元，并计算与沉没辊接触带钢的任意(i,j)子单元坐标值(θi，zj)；单元的编号初值i＝1，j＝1，打滑单元数nh＝0，单元(1，1)的旋入段张力的初值 辊的驱动力矩的初始值Ms＝0；其中，T0为带钢生产工艺参数，表示带钢初始设定张力；

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

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

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

若带钢单元打滑临界条件成立，则单元(i,j)不发生打滑；若带钢单元打滑临界条件不成立，则单元(i,j)发生打滑，并令nh＝nh+1，判断i≤m且j≤n是否成立；若i≤m且j≤n成立，则令i＝i+1，j＝j+1，Ti,j＝T′i,j，其中，T′i,j为带钢单元旋出端的张力，并返回步骤F；若i≤m且j≤n不成立，则判断辊打滑条件 是否成立；若辊打滑条件成立，则认定辊发生打滑；若辊打滑条件不成立，则认定辊未发生打滑；式中，ξn为带钢子单元发生打滑数量的最大占比，根据机组实际工况及生产效率预先设定；

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

步骤G、判断辊发生转动条件δ≥[δ]是否成立，其中，转动动力因数 [δ]为转动动力因数临界值；若辊发生转动条件成立，则辊未发生转动失效；若辊发生转动条件不成立，则辊发生转动失效。

2.根据权利要求1所述的一种热镀锌沉没辊系转动失效预报方法，其特征在于，沉没辊系包括包：沉没辊、稳定辊及矫正辊。

3.根据权利要求2所述的一种热镀锌沉没辊系转动失效预报方法，其特征在于，所述方法还包括：步骤I、判断k≤3是否成立；其中，沉没辊、稳定辊及矫正辊的辊号为1、2、3；若k≤3成立，则k＝k+1，并转入步骤B判断辊是否发生转动失效。

4.根据权利要求1所述的一种热镀锌沉没辊系转动失效预报方法，其特征在于，计算与沉没辊接触带钢的任意(i,j)子单元坐标值(θi，zj)：其中，L为辊的设备参数，表示辊身长，B、θ、θ‑为带钢生产工艺参数，B表示带材宽度，θ表示带钢与辊包角初始值，θ‑表示包角的负角。

5.根据权利要求4所述的一种热镀锌沉没辊系转动失效预报方法，其特征在于，确定辊面最大静摩擦系数τmax和平均摩擦系数τ，包括：在辊面沿辊身方向在带钢接触辊的部分均匀选取ms个点，并通过测量传感器测量这些点所在区域锌渣厚度hk，其中k＝1,2,...,ms，ms≤n根据测量得到的hk计算各辊的辊面平均锌渣沉积厚度计算辊面子单元所在辊的全周向区域摩擦系数τ： τ0为带钢生产工艺参数，表示带钢与辊的辊面间的初始摩擦系数；

计算辊面的最大静摩擦系数τmax： 其中，hmax为hk中的最大值，

6.根据权利要求4所述的一种热镀锌沉没辊系转动失效预报方法，其特征在于，确定带钢单元与辊面间接触压力Ni,j、滚动摩擦力fi,j及总接触压力N分量，包括：计算带钢单元受到重力与浮力的合力Δgi,j：

Δgi,j＝(ρFe‑ρZn)g·Δθ·R·Δz·H；其中，依据能量平衡建立方程：

式中，I为带钢单元的惯性矩， E为高温下带弹性模量；

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

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

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

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

计算总接触压力X向分力Nx、Y向分力Ny：

7.根据权利要求6所述的一种热镀锌沉没辊系转动失效预报方法，其特征在于，Mg的计算方式如下：Mg＝τFgRcosγ；其中，Fg、γ为带钢生产工艺参数，分别表示刮刀力、刮刀力与辊径向的夹角。

8.根据权利要求6所述的一种热镀锌沉没辊系转动失效预报方法，其特征在于，Md的计算方式如下：计算辊的最大承载量[P]、辊受到的外载荷P：

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

判断辊的轴端的轴瓦与轴颈是否接触：判断P≤[P]是否成立；若成立，则计算轴端摩擦阻力矩： 若不成立，则转入计算轴端摩擦阻力矩： 式中，ζ1为速度影响系数；ζ2为镀锌量影响系数，Qs为带钢生产工艺参数，表示带钢表面单位面积上镀锌量，μ0为辊的设备参数，表示辊的轴套与轴瓦间初始摩擦系数。