1.一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，所述方法包括以下步骤：步骤A、设定沉没辊系中辊的辊号λ，初值λ＝1；

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

步骤C、根据几何法推导的公式，并结合实际工况的矫正辊插入量确定沉没辊系各辊在工作状态下的包角θ；

步骤D、将辊接触区带钢沿辊周向划分m个横向条元；再沿任意辊接触区带钢横向条元的横向划分n个子单元，并计算与沉没辊接触带钢的任意(i,j)子单元坐标值(θi，zj)；

步骤E、计算辊号为λ的辊对带钢单元的接触压力Ni,j和滚动摩擦力fi,j，计算公式为：fi,j＝τk·Ni,j；

其中，Ti,j为带钢单元旋入端的张力，Ti,j＝T0；T0为带钢初始设定张力；T′i,j为带钢单元旋出端的张力， Δgi,j为带钢单元受到重力与浮力的合力，H、ρFe、ρZn、v、B为带钢生产工艺参数，分别表示带材的厚度、带钢密度、锌锅内高温下锌液密度、带钢运行速度、带材宽度；R为辊的设备参数，表示辊半径；τk为辊面子单元所在辊的全周向区域摩擦系数；

步骤F、计算带钢对各辊的辊面滚动摩擦力矩M以及各辊受到的总接触压力N：步骤G、计算各辊的轴端支撑力F以及摩擦阻力矩Md，计算公式为：其中，

Md＝μFd/2；

其中，ΔG为沉没辊系重力与浮力的合力，ΔG＝G‑ρZngV；其中，g为重力加速常量；V、G、L、d、μ为辊的设备参数，分别表示辊的有效体积、辊的重力、辊的辊身长、辊轴颈直径、辊的轴颈与轴瓦间初始摩擦系数；B为带钢生产工艺参数，表示带材宽度；

步骤H、判断λ≤3是否成立，若λ≤3成立，则令λ＝λ+1，并返回执行步骤B，进行下一个辊的受力预报；若λ≤3不成立，则输出沉没辊系各辊受到的总接触压力N、带钢的摩擦驱动力矩M、轴端支持力F及其摩擦阻力矩Md。

2.根据权利要求1所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，所述沉没辊系包括：沉没辊、稳定辊及矫正辊。

3.根据权利要求2所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，确定沉没辊系各辊在工作状态下的包角，包括：判断λ＝1是否成立，若λ＝1成立，则：

* *

并令R＝R，θ＝θ， 辊包角θ的正向角θ+、负向角θ‑为：若不成立，则

判断λ＝2是否成立，若λ＝2成立，则

若λ＝2不成立，则

* *

其中，θ0为带钢与辊包角初始值，Δs为矫正辊插入量，R、θ、 是过程变量，初值均为0。

4.根据权利要求3所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，分别计算与辊接触带钢的任意(i,j)子单元坐标值(θi,zj)包括：其中，θ为辊的包角，θ‑为辊包角的负向角。

5.根据权利要求1所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，还包括：计算带钢对各辊的接触压力N的方向角φ ，计算公式为：

6.根据权利要求1所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，辊面子单元所在辊的全周向区域摩擦系数τk的计算方式如下：*

在各辊沿其辊面轴向均匀选取若干个子单元n ，并通过测量传感器测量得出在各辊表* *面上对应这些点的锌渣沉积厚度hk，其中k＝1,2,...,n，n≤n；

计算辊面子单元所在辊的全周向区域摩擦系数τk，计算公式为： 其中，为辊的辊面平均锌渣沉积厚度， τ为带钢生产工艺参数，表示带钢与辊的辊面间的初始摩擦系数。