CN102728624B - 一种精轧带钢负荷分配设定方法 - Google Patents

Abstract

一种精轧带钢负荷分配设定方法，涉及精轧带钢生产控制领域，解决现有压下分配模式中轧制状况改变时各机架轧制力比例波动大，轧制力分配模式中轧制力模型偏差大时各机架厚度分配波动大的问题，本方法实现了压下分配及轧制力分配模式的结合，用压下分配确定的分布范围来约束轧制力分配结果，在减小各机架厚度分配波动的同时尽量满足各机架轧制力按比例下降的要求，在轧制力模型偏差不大时，通过压下率学习表的学习使两种分配模式结果越来越一致，迭代次数越来越少；在轧制力模型偏差较大时，通过压下分配的限幅有效控制各机架压下率的波动范围，本方法综合两种模式的优点，同时克服两者各自的缺点，有效提高精轧带钢轧制规程设定精度和轧制稳定性。

Description

一种精轧带钢负荷分配设定方法

技术领域

本发明涉及精轧带钢生产控制领域，尤其涉及一种精轧带钢负荷分配设定方法。

背景技术

带钢热连轧机组中，轧制规程的制定对于提高产品质量起着重要的作用，其合理与否，对产品质量的高低、轧制设备调整的难易、机组运转的可靠性和生产过程的稳定性均有直接影响。轧制规程一般包括压下制度，速度制度，张力制度，温度制度等。其中压下制度指各机架（或各道次）压下量的分配制度，业内通常称之为负荷分配或厚度分配。

负荷分配是轧制规程的核心，它直接影响到板形、板厚精度等产品质量，负荷分配还对轧制能耗、辊耗、生产过程的稳定性和作业率等项指标有重要影响。带钢热连轧机组的负荷分配方法从以前的经验表格法开始，还经历了能耗法、动态负荷分配法、负荷分配系数法、人工智能方法等几个阶段。

目前，现代化热连轧机组负荷分配较普遍采用的方法是负荷分配系数法，它按一定的负荷分配系数给定各机架的压下量和轧件厚度。负荷分配系数可以参照品种和规格来确定，也可以采用离线优化的方法来确定。负荷分配系数法主要包括压下量、压下率、轧制力和功率等类型。其中轧制力、功率类型需要迭代计算，为方便分析都归为轧制力分配模式；压下量、压下率类型则可直接获得各机架的压下量和厚度，都归为压下分配模式。这两种模式的优缺点：

压下分配的优点在于各机架厚度分配波动小，计算过程简单；缺点在于轧制状况改变（如换辊）时各机架轧制力比例波动大，不能消除轧制力倒挂，不能保证产品板形。

轧制力分配的优点在于各机架轧制力下降比例基本恒定，有利于产品板形和轧制稳定性；缺点在于当轧制力模型偏差大或工艺条件变化大时各机架厚度分配波动大，计算过程复杂。

目前，有关精轧带钢负荷分配的设定方法也出现一些公开的专利文件，例，名称为《一种热轧板带连轧轧制规程动态设定的方法》（CN101733289A）的文件中记载的方法：将压下分配系数作为寻优起始点，计算此时相应的负荷，并与目标负荷分配系数相比较，根据差值调整各道次压下量；经过多次比较与调整，达到设定负荷分配逼近目标负荷分配的目的。其局限在于：①每次只调整两个机架的分配，机架间协调性不好，调整效率较低；②仅将压下分配系数作为寻优起始点，本质上还是一种轧制力分配，不能克服轧制力分配的缺点。

发明内容

本发明的目的是为了解决现有技术存在的问题，提出一种精轧带钢负荷分配设定方法，从而提高精轧带钢轧制规程设定精度和轧制稳定性。

本发明的目的是通过以下技术方案实现的：

一种精轧带钢负荷分配设定方法，包括以下步骤：

S1、根据压下率分配系数表和压下率学习表，按压下分配模式计算精轧各机架的压下分配结果，所述压下分配结果指压下分配的压下率和压下分配的出口厚度S2、将上述得到的各机架压下分配的出口厚度作为初始值，按轧制力分配模式迭代计算各机架的绝对压下量，直到轧制力比例满足目标要求，得到各机架的轧制力分配结果，所述轧制力分配结果包括轧制力分配的出口厚度及轧制力分配的压下率S3、对所述各机架的轧制力分配结果进行限幅处理，使其偏离所述压下分配结果在一定范围内，从而使各机架的最终压下率在由压下分配确定的一定范围内，只要没超出范围就取轧制力分配结果，超出范围则将其限制在边界处；S4，带钢头部穿带完成并收集到实际轧制力后，判断各机架轧制力设定值与实际轧制力的偏差；S5，判定各机架偏差值是否都在θ以内，θ为偏差控制参数，取5%-10%，如果不满足各机架轧制力偏差都在θ以内的条件，则不进行压下率学习，此时，所述压下率学习表中相应层别的数据不进行更新，流程进入结束步骤；S6，如果各机架偏差值在θ以内，则对本次设定的最终压下率进行学习，将学习值更新到压下率学习表格相应层别，经学习过程，该层别轧制力会逐步达到目标轧制力比例的要求；S7，负荷分配设定结束。

所述步骤S1包括以下步骤：

S11，取表，从压下率分配系数表和压下率学习表取到各机架的压下率和压下率学习值，两者进行相加，得到各机架的压下率分配表格值

S12，进行相对化计算，得到各机架的压下分配结果：根据精轧各机架的压下量之和必须等于精轧的总压下量，得到压下率的比例因子xrk为：

$\mathrm{xrk}=\frac{-\mathrm{kb}\·{\mathrm{xs}}_1+\sqrt{{(\mathrm{kb}\·{\mathrm{xs}}_1)}^2+4\·\mathrm{ka}\·{\mathrm{xs}}_2\·1\mathrm{nx}}}{2\·\mathrm{ka}\·{\mathrm{xs}}_2}$      公式1

其中，h_n为终轧厚度，h₀为中间坯厚度；

n为精轧道次数，i为机架号；

kb=0.906501，ka=0.959597。

各机架的压下率分配表格值相对化后的压下率为压下分配的压下率

${\mathrm{\ϵ}}_{\mathrm{red}}^i=\mathrm{xrk}\·{\mathrm{\ϵ}}_0^i$        公式2

根据压下分配的压下率得到前面n-1个机架各自的压下分配的出口厚度 $h_{\mathrm{red}}^i:$

$\begin{array}{cc}{h}_{\mathrm{red}}^i=h_0\·\underset{i=1}{\overset{i}{\mathrm{\Π}}}(1-{\mathrm{\ϵ}}_{\mathrm{red}}^i)& i=\mathrm{1,2},...,n-1.\end{array}$       公式3

所述步骤S11中各机架的压下率分配表格值取得过程为：根据当前带钢的钢种大类、厚度等级和宽度等级信息，从压下分配系数表相应层别取到精轧各机架的压下率其中i为机架号；根据当前带钢的出钢记号、厚度等级和宽度等级信息，从压下率学习表相应层别取到精轧各机架的压下率学习值其中i为机架号，所述压下率学习表是由程序自动学习得到的数据表，表中数据通过学习不断更新，初始值取0；所述各机架的压下率分配表格值为：

${\mathrm{\ϵ}}_0^i={\mathrm{\ϵ}}_{\mathrm{ori}}^i+{\mathrm{\ϵ}}_{\mathrm{lrn}}^i.$

所述步骤S2包括以下步骤：

S21，根据当前带钢的钢种大类和厚度等级信息，从轧制力分配系数表相应层别获取到各机架初始轧制力分配系数α_0i；

S22，进行操作修正处理，获得各机架修正后的轧制力分配系数α_1i：从操作画面HMI收集各机架的操作修正因子ζ_i，各机架修正后的轧制力分配系数α_1i为：

${\mathrm{\α}}_{1i}={\mathrm{\α}}_{0i}\·(1+{\mathrm{\ζ}}_i)$      公式4

其中，ζ_i取值范围为-20%～20%；

进行归一化，归一化后的轧制力分配系数αi为：

S23，计算在当前厚度分配下各机架的轧制力和轧制力对绝对压下量的导数：计算轧制力设：i为机架号，j为迭代次数，为第i机架第j次迭代的轧制力，第i机架第j次迭代的出口厚度，第i-1机架第j次迭代的出口厚度；第一次迭代计算前第i机架出口厚度初始值为上述公式(3)得到的压下分配的出口厚度 在其它工艺参数固定的情况下，轧制力是入口厚度和出口厚度的函数，存在：

计算轧制力对绝对压下量的导数

其中，为绝对压下量正扰动后的轧制力；为绝对压下量负扰动后的轧制力，为第i机架第j次迭代的绝对压下量；存在以下函数关系： $p_{i1}^j=p(h_{i-1}^j,h_{i-1}^j-{\mathrm{\Δh}}_{i1}^j),p_{i2}^j=p(h_{i-1}^j,h_{i-1}^j-{\mathrm{\Δh}}_{i2}^j),$ 其中，为正扰动后的绝对压下量，为负扰动后的绝对压下量，

${\mathrm{\Δh}}_{i1}^j=(1+\frac{\mathrm{dt}}{2})\·{\mathrm{\Δh}}_i^j,{\mathrm{\Δh}}_{i2}^j=(1-\frac{\mathrm{dt}}{2})\·{\mathrm{\Δh}}_i^j,$ 其中dt为扰动，取1%；

S24，根据上述各机架的轧制力和轧制力对绝对压下量导数的计算结果，计算各机架绝对压下量的修正值：

轧制分配迭代计算公式：

$\mathrm{\δ}\left({\mathrm{\Δh}}_i^j\right)=({\mathrm{\Σp}}_i^j\·\frac{{\mathrm{\α}}_i}{\mathrm{\Σ}{\mathrm{\α}}_i}-p_i^j)/\frac{{\∂p}_i^j}{\∂\left(\mathrm{\Δ}{h}_i^j\right)}-\mathrm{\Σ}((\mathrm{\Σ}{p}_i^j\·\frac{{\mathrm{\α}}_i}{\mathrm{\Σ}{\mathrm{\α}}_i}-p_i^j)/\frac{{\∂p}_i^j}{\∂\left({\mathrm{\Δh}}_i^j\right)})\·\frac{{\mathrm{\Δh}}_i^j}{\mathrm{\Σ\Δ}{h}_i^j}$     公式5

其中：为第i机架第j次迭代的绝对压下量的修正值；

S25，各机架绝对压下量的修正值得到后，计算各机架绝对压下量的更新值，该更新值为第j+1次迭代的绝对压下量计算公式为：

${\mathrm{\Δh}}_i^{j+1}={\mathrm{\Δh}}_i^j+{\mathrm{damp}}_j\·\mathrm{\δ}\left({\mathrm{\Δh}}_i^j\right)$      公式6

其中damp_j为阻尼系数，damp_j=damp_mpy·β+(1-β)·(1-e^-j)，

damp_mpy取1.0，β取0.6；

S26，各机架绝对压下量的更新值得到后，计算各机架出口厚度的更新值，该更新值为第j+1次迭代的出口厚度

以及中间坯厚度h₀为已知量，根据仅第一个机架其余都是前一个机架的输出为后一个机架的输入，从第一个机架计算到最后一个机架，得到各机架出口厚度的更新值，这个值将用于下一次迭代计算轧制力用，其中指第i-1机架第j+1次迭代的出口厚度；

S27，判断轧制力比例是否满足收敛条件：

轧制力分配迭代计算的收敛条件为：

$|\frac{p_i^j}{\mathrm{\Σ}{p}_i^j}-\frac{{\mathrm{\α}}_i}{\mathrm{\Σ}{\mathrm{\α}}_i}|\≤\frac{\mathrm{\τ}}{{\mathrm{\Σ\α}}_i}$        公式7

其中，τ为正数，取0.01；当式(7)成立，或者迭代次数j超过设定次数时结束迭代计算，否则继续执行所述S23步骤，同时迭代次数j累加1；

S28，上述步骤S27中如果轧制力比例满足收敛条件，则判断某机架的轧制力是否超过设备允许的最大轧制力：

如果存在轧制力超限：

当第i机架的轧制力p_i超过设备允许的最大轧制力时，进入S281步骤，调整该机架的轧制力分配系数，其余机架的分配系数不作调整，存在：

时，则执行 ${\mathrm{\α}}_i^{\′}={\mathrm{\α}}_i\·\frac{1}{(p_i/p_i^{\max }-1)\·0.8+1}$    公式8

其中，α_i为所述归一化后的轧制力分配系数α_i，α_i’为超限调整后的轧制力分配系数；

α_i’得到后回到步骤S23，用调整后的的轧制力分配系数α_i’重新进行绝对压下量迭代计算，同时迭代次数j赋值0，取步骤S27最后一次迭代的绝对压下量和出口厚度作为绝对压下量和出口厚度的初始值；

如果不存在轧制力超限，则直接执行下一步骤；

S29，绝对压下量迭代计算和超限检查结束，得到轧制力分配结果，取步骤S27最后一次迭代的出口厚度作为各机架轧制力分配的出口厚度同时得到各机架轧制力分配的压下率为：式中为第i机架轧制力分配的出口厚度，为第i-1机架轧制力分配的出口厚度。

所述步骤3包括以下步骤：

S31，对轧制力分配结果进行限幅处理：

${\mathrm{\&epsiv;}}_{\mathrm{xf}}^i=\left\{\begin{array}{cc}{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})& {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i>{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})\\ {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i& {\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})\&le;{\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i\&le;{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})\\ {\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})& {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i<{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})\end{array}\right.$       公式9

其中为限幅处理后压下率，η为限幅控制参数，取5%～10%；

S32，对限幅处理后压下率进行相对化计算：

上述步骤S31中压下率如果被修改，计算终轧厚度h′_n与目标终轧厚度h_n间会存在一定偏差，还要进行一轮相对化计算：

计算终轧厚度h′_n：

${h_n}^{\&prime;}=h_0\&CenterDot;\underset{i=1}{\overset{n}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i)$

计算终轧相对化因子(1-xc)，其中

$\mathrm{xc}=\frac{1n\frac{h_n}{{h_n}^{\&prime;}}}{\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}\frac{{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i}{1-{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i}}$

终轧相对化后的压下率称为最终分配压下率：

${\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i=(1-\mathrm{xc})\&CenterDot;{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i$

取作为各机架的最终输出的压下率，同时得到各机架的最终出口厚度

$\begin{array}{cc}{h}_{\mathrm{xdh}}^i:h_{\mathrm{xdh}}^i=h_0\&CenterDot;\underset{i=1}{\overset{i}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i)& i=\mathrm{1,2},...,n-1,\end{array}$

根据最终分配压下率和最终出口厚度进行精轧带钢的轧制规程设定计算。

所述步骤S6中，最终压下率进行学习的过程为：

求各机架最终分配与压下分配的压下率之差

${\mathrm{\&epsiv;}}_{\mathrm{dif}}^i=\frac{{\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i-{\mathrm{\&epsiv;}}_{\mathrm{red}}^i}{\mathrm{xrk}}$      公式10

其中，为最终分配的压下率，为压下分配的压下率，xrk为公式（1）中的压下率的比例因子；

得到轧制规程设定后压下率学习表相应层别更新的学习值

${\mathrm{\&epsiv;}}_{\mathrm{lrn}}^{i,\mathrm{new}}={\mathrm{\&epsiv;}}_{\mathrm{lrn}}^i\&CenterDot;\mathrm{\&mu;}+{\mathrm{\&epsiv;}}_{\mathrm{dif}}^i\&CenterDot;(1-\mathrm{\&mu;})$      公式11

其中，μ为平滑系数，取0.7；为轧制规程设定前压下率学习表相应层别原来的学习值。

本发明的有益效果：

本发明的精轧带钢负荷分配设定方法，实现“压下分配+轧制力分配”的综合模式，用压下分配确定的分布范围来约束轧制力分配结果，综合了两种模式的优点，同时克服两者各自的缺点，在减小各机架厚度分配波动的同时尽量满足各机架轧制力按比例下降的要求，达到提高精轧带钢轧制规程设定精度和轧制稳定性的目的。

为进一步说明本发明的上述目的、结构特点和效果，以下将结合附图对本发明进行详细说明。

附图说明

图1为本发明的精轧带钢负荷分配设定方法流程图；

图2为图1流程中计算各机架绝对压下量的具体流程图；

图3为轧制力分配和最终分配由轧制力模型偏差导致的压下率波动。

具体实施方式

下面结合实施例的附图对本发明的具体实施方式进行详细说明。

参见图1，图1为本发明的精轧带钢负荷分配设定方法流程图，本发明的精轧带钢负荷分配设定方法，包括以下步骤：

S1、根据压下率分配系数表和压下率学习表，按压下分配模式计算精轧各机架的压下分配结果，所述压下分配结果指压下分配的压下率和压下分配的出口厚度其包括以下步骤：

S11，取表，从压下率分配系数表和压下率学习表取到各机架的压下率和压下率学习值，两者进行相加，得到各机架的压下率分配表格值，具体过程为：

根据当前带钢的钢种大类、厚度等级和宽度等级信息，从压下分配系数表相应层别取到精轧各机架（例：F1～F7）的压下率（i为机架号），所述压下率分配系数表是业内常用表，表中数据在生产前确定，一般为考虑设备能力和产品特性的经验数据，可收集实际生产中的优秀案例；

根据当前带钢的出钢记号、厚度等级和宽度等级信息，从压下率学习表相应层别取到精轧各机架（F1～F7）的压下率学习值（i为机架号），所述压下率学习表是由程序自动学习得到的数据表，表中数据通过学习不断更新，初始值取0；

将两者相加得到精轧各机架的压下率分配表格值即（i为机架号）。

两层表格设计的好处在于通过压下率学习表可以反映同一钢种大类下不同出钢记号的差异(一个钢种大类包括若干出钢记号)。

S12，进行相对化计算，得到各机架的压下分配结果，即压下分配的压下率和压下分配的出口厚度具体过程为：

为了得到各种中间坯厚度和终轧厚度下的具体压下率数据，需要对上述压下率分配表格值进行相对化计算，根据精轧各机架的压下量之和必须等于精轧的总压下量，得到压下率的比例因子xrk为：

$\mathrm{xrk}=\frac{-\mathrm{kb}\&CenterDot;{\mathrm{xs}}_1+\sqrt{{(\mathrm{kb}\&CenterDot;{\mathrm{xs}}_1)}^2+4\&CenterDot;\mathrm{ka}\&CenterDot;{\mathrm{xs}}_2\&CenterDot;1\mathrm{nx}}}{2\&CenterDot;\mathrm{ka}\&CenterDot;{\mathrm{xs}}_2}---\left(1\right)$

其中，h_n为终轧厚度，h₀为中间坯厚度（h_n与h₀是精轧轧制规程设定计算前已知的值）；

n为精轧道次数，i为机架号；

kb=0.906501，ka=0.959597。

各机架的压下率分配表格值相对化后的压下率（即压下分配的压下率）为：

${\mathrm{\&epsiv;}}_{\mathrm{red}}^i=\mathrm{xrk}\&CenterDot;{\mathrm{\&epsiv;}}_0^i---\left(2\right)$

根据压下分配的压下率可得到前面n-1个机架各自的压下分配的出口厚度 $h_{\mathrm{red}}^i:$

$\begin{array}{cc}{h}_{\mathrm{red}}^i=h_0\&CenterDot;\underset{i=1}{\overset{i}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{red}}^i)& i=\mathrm{1,2},...,n-1.\end{array}---\left(3\right)$

式（3）是根据各机架入口厚度、出口厚度与压下率之间的关系得到，是业内使用的基本公式。

S2、以上述S12中得到的各机架压下分配的出口厚度为初始值，按轧制力分配模式迭代计算各机架的绝对压下量，直到轧制力比例满足目标要求，得到各机架的轧制力分配结果，所述轧制力分配结果包括轧制力分配的出口厚度及轧制力分配的压下率具体流程参见图2，其中包括绝对压下量迭代计算过程（S21-S27）及超限检查过程（S28）：

S21，根据当前带钢的钢种大类和厚度等级信息，从轧制力分配系数表相应层别获取到各机架初始轧制力分配系数α_0i。

上述分配系数表中数据在生产前确定，主要考虑带钢轧制稳定性和板形指标等因素，可以在经验数据基础上不断收集优秀轧制案例进行优化。上述步骤S1中的压下率分配系数表为压下分配模式所必须，本步骤中的轧制力分配系数表为轧制力分配模式所必须，也是业内所公知的数据表，两者数据是独立的，都是工艺经验数据，本发明实现了该两种分配的综合。

S22，进行操作修正处理，获得各机架修正后的轧制力分配系数α_1i。

从操作画面HMI收集各机架的操作修正因子ζ_i，各机架修正后的轧制力分配系数α_1i为：

α_1i=α_0i·(1+ζ_i)   (4)

其中，ζ_i取值范围为-20%～20%。

进行归一化，归一化后的轧制力分配系数α_i为：

这里的归一化指所有机架（例F1～F7共7个机架）中最大的α_1i值归为1，其余项目进行相应处理。

S23，计算在当前厚度分配下各机架的轧制力和轧制力对绝对压下量的导数：

(1)计算轧制力

下面公式中设置的变量说明：i为机架号，j为迭代次数，为第i机架第j次迭代的轧制力，第i机架第j次迭代的出口厚度，第i-1机架第j次迭代的出口厚度。

第一次迭代计算前第i机架出口厚度初始值为上述公式(3)得到的压下分配的出口厚度即

在其它工艺参数固定的情况下，轧制力是入口厚度和出口厚度的函数，即存在：（该函数关系为业内公知，在此不作说明），轧制力的大小受到轧制力模型精度的影响；要注意的是，各机架厚度分配的变化会导致各机架温度的变化，在计算各机架轧制力前，应重新计算一次各机架的温度（温度与轧制力互相影响，其互相影响的关系为业内所知，在此不作说明）；

（2）计算轧制力对绝对压下量的导数

计算公式为：该计算过程的公式原理为轧制力变化量除以绝对压下量变化量，即导数的定义，是业内公知计算过程，在此作简要说明。

公式中，为绝对压下量正扰动后的轧制力；为绝对压下量负扰动后的轧制力，为第i机架第j次迭代的绝对压下量；

存在以下函数关系： $p_{i1}^j=p(h_{i-1}^j,h_{i-1}^j-{\mathrm{\&Delta;h}}_{i1}^j),p_{i2}^j=p(h_{i-1}^j,h_{i-1}^j-{\mathrm{\&Delta;h}}_{i2}^j)$

其中，为正扰动后的绝对压下量，为负扰动后的绝对压下量，

${\mathrm{\&Delta;h}}_{i1}^j=(1+\frac{\mathrm{dt}}{2})\&CenterDot;{\mathrm{\&Delta;h}}_i^j,{\mathrm{\&Delta;h}}_{i2}^j=(1-\frac{\mathrm{dt}}{2}){\&CenterDot;\mathrm{\&Delta;h}}_i^j$

其中dt为很小的扰动，可取1%。

S24，根据上述各机架的轧制力和轧制力对绝对压下量导数的计算结果，计算各机架绝对压下量的修正值。

轧制分配迭代计算的核心公式：

$\mathrm{\&delta;}\left({\mathrm{\&Delta;h}}_i^j\right)=({\mathrm{\&Sigma;p}}_i^j\&CenterDot;\frac{{\mathrm{\&alpha;}}_i}{\mathrm{\&Sigma;}{\mathrm{\&alpha;}}_i}-p_i^j)/\frac{{\&PartialD;p}_i^j}{\&PartialD;\left(\mathrm{\&Delta;}{h}_i^j\right)}-\mathrm{\&Sigma;}((\mathrm{\&Sigma;}{p}_i^j\&CenterDot;\frac{{\mathrm{\&alpha;}}_i}{\mathrm{\&Sigma;}{\mathrm{\&alpha;}}_i}-p_i^j)/\frac{{\&PartialD;p}_i^j}{\&PartialD;\left({\mathrm{\&Delta;h}}_i^j\right)})\&CenterDot;\frac{{\mathrm{\&Delta;h}}_i^j}{\mathrm{\&Sigma;\&Delta;}{h}_i^j}---\left(5\right)$

其中：为第i机架第j次迭代的绝对压下量的修正值。

上述公式（5）是根据一种改进Newton法导出的结果，选用各道次的绝对压下量作为变量，将n维的非线性方程变为一维的n个非线性方程，并对每个非线性方程用Newton迭代法进行求解。它借鉴了Newton-Raphson法有良好收敛性的特点，同时又可以避免求Jacobian矩阵及其逆矩阵，使算法得到有效简化，计算速度较快，精度较高。

S25，各机架绝对压下量的修正值得到后，计算各机架绝对压下量的更新值，也就是第j+1次迭代的绝对压下量

计算公式为：

$\mathrm{\&Delta;}{h}_i^{j+1}=\mathrm{\&Delta;}{h}_i^j+{\mathrm{damp}}_j\&CenterDot;\mathrm{\&delta;}\left(\mathrm{\&Delta;}{h}_i^j\right)---\left(6\right)$

其中damp_j为阻尼系数，damp_j=damp_mpy·β+(1-β)·(1-e^-j)，

damp_mpy可取1.0，β可取0.6。

每次迭代时仅吸收绝对压下量修正值的一定百分比，随着迭代次数增加，吸收百分比也增加。

S26，各机架绝对压下量的更新值得到后，计算各机架出口厚度的更新值，也就是第j+1次迭代的出口厚度

以及中间坯厚度h₀为已知量，根据仅第一个机架其余都是前一个机架的输出为后一个机架的输入，从第一个机架计算到最后一个机架，得到各机架出口厚度的更新值，这个值将用于下一次迭代计算轧制力用，其中指第i-1机架第j+1次迭代的出口厚度；

例对F1机架，i=1,无论哪次迭代，都有不变；

对F2机架，i=2，F1输出作为F2输入；

对F3机架，i=3,依次递推.......。

所以归纳为仅当F1时，其余都是前一个机架的输出为后一个机架的输入。

S27，判断轧制力比例是否满足收敛条件。

轧制力分配迭代计算的收敛条件为：

$|\frac{p_i^j}{\mathrm{\&Sigma;}{p}_i^j}-\frac{{\mathrm{\&alpha;}}_i}{\mathrm{\&Sigma;}{\mathrm{\&alpha;}}_i}|\&le;\frac{\mathrm{\&tau;}}{{\mathrm{\&Sigma;\&alpha;}}_i}---\left(7\right)$

其中，τ为很小的正数，可取0.01。

当式(7)成立，或者迭代次数j超过设定次数时结束迭代计算，否则继续执行S23步骤，同时迭代次数j累加1。举例，迭代次数j的设定次数为6次，即j的取值范围为0≤j≤5，第1次迭代计算j取0，当式(7)成立或者迭代次数j超过6次时迭代计算结束；如果第一次迭代计算还没达到收敛条件，继续执行S23步骤，同时迭代次数j累加1。

迭代计算结束后，可得到各机架的绝对压下量和出口厚度

S28，上述步骤S27中如果轧制力比例满足收敛条件，则判断某机架的轧制力是否超过设备允许的最大轧制力：

如果存在轧制力超限：

当第i机架的轧制力pi超过设备允许的最大轧制力时，进入S281步骤，调整该机架的轧制力分配系数，其余机架的分配系数不作调整，存在：

时，则执行 ${\mathrm{\&alpha;}}_i^{\&prime;}={\mathrm{\&alpha;}}_i\&CenterDot;\frac{1}{(p_i/p_i^{\max }-1)\&CenterDot;0.8+1}---\left(8\right)$

其中，α_i为所述归一化后的轧制力分配系数α_i，α_i’为超限调整后的轧制力分配系数；

α_i’得到后回到步骤S23，用调整后的的轧制力分配系数α_i’重新进行绝对压下量迭代计算。同时迭代次数j赋值0，取S27最后一次迭代的绝对压下量和出口厚度作为绝对压下量和出口厚度的初始值；

如果不存在轧制力超限，则直接执行下一步骤S29。

S29，绝对压下量迭代计算和超限检查结束，得到轧制力分配结果，取步骤S27最后一次迭代的出口厚度作为各机架轧制力分配的出口厚度同时得到各机架轧制力分配的压下率为：式中为第i机架轧制力分配的出口厚度，为第i-1机架轧制力分配的出口厚度。

在S2步骤中以压下分配结果为轧制力分配的初值可以减少迭代次数，轧制力分配可以控制各机架轧制力下降的比例，进而改善精轧带钢的轧制稳定性及板形指标。

S3、对上述步骤S2中各机架的轧制力分配结果进行限幅处理，使其偏离上述步骤S1压下分配结果在一定范围内，从而使各机架的最终压下率在由压下分配确定的一定范围内，只要没超出范围就取轧制力分配结果，超出范围则将其限制在边界处。

S31，对轧制力分配结果进行限幅处理：

${\mathrm{\&epsiv;}}_{\mathrm{xf}}^i=\left\{\begin{array}{cc}{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})& {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i>{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})\\ {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i& {\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})\&le;{\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i\&le;{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})\\ {\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})& {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i<{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})\end{array}\right.---\left(9\right)$

其中为限幅处理后压下率，η为限幅控制参数，取5%～10%；

S32，对限幅处理后压下率进行相对化计算：

上述步骤S31中压下率如果被修改，计算终轧厚度与目标终轧厚度h_n间会存在一定偏差，还要进行一轮相对化计算：

计算终轧厚度

${h_n}^{\&prime;}=h_0\&CenterDot;\underset{i=1}{\overset{n}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i)$

计算终轧相对化因子(1-xc)，其中

$\mathrm{xc}=\frac{1n\frac{h_n}{{h_n}^{\&prime;}}}{\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}\frac{{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i}{1-{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i}}$

终轧相对化后的压下率称为最终分配压下率：

${\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i=(1-\mathrm{xc})\&CenterDot;{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i$

取作为各机架的最终输出的压下率，同时得到各机架的最终出口厚度

$\begin{array}{cc}{h}_{\mathrm{xdh}}^i:h_{\mathrm{xdh}}^i=h_0\&CenterDot;\underset{i=1}{\overset{i}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i)& i=\mathrm{1,2},...,n-1,\end{array}$

根据最终分配压下率和最终出口厚度进行精轧带钢的轧制规程设定计算。

S4，带钢头部穿带完成并收集到实际轧制力后，判定轧制力设定值（该轧制力设定值由轧制规程计算时的轧制力模型计算得到，是公有知识，与上述S23中计算轧制力一样）与实际轧制力的偏差。

S5，判断各机架偏差值是否都在θ以内，θ为偏差控制参数，可取5%-10%，如果不满足各机架轧制力偏差都在θ以内的条件，则不进行压下率学习，此时，压下率学习表（下述表格2）中相应层别的数据不进行更新，流程进入结束步骤。

S6，如果各机架偏差值在θ以内，则对本次设定的最终压下率进行学习，学习值更新到压下率学习表格相应层别，为：

${\mathrm{\&epsiv;}}_{\mathrm{dif}}^i=\frac{{\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i-{\mathrm{\&epsiv;}}_{\mathrm{red}}^i}{\mathrm{xrk}}---\left(10\right)$

其中为第i机架最终分配与压下分配的压下率之差，xrk为公式（1）中的压下率的比例因子。

得到轧制规程设定后压下率学习表相应层别更新的学习值为：

${\mathrm{\&epsiv;}}_{\mathrm{lrn}}^{i,\mathrm{new}}={\mathrm{\&epsiv;}}_{\mathrm{lrn}}^i\&CenterDot;\mathrm{\&mu;}+{\mathrm{\&epsiv;}}_{\mathrm{dif}}^i\&CenterDot;(1-\mathrm{\&mu;})---\left(11\right)$

其中，μ为平滑系数，可取0.7；为压下率学习表（下述表格2）相应层别原来的学习值（轧制规程设定前）。

S7，计算结束，在将更新到压下率学习表（表格2）相应层别后，由于学习的作用，该层别轧制力会逐步达到目标轧制力比例的要求。

下面举实施例说明。

实施例1：

着重说明整个计算流程。最终结果为经过压下分配限幅后的轧制力分配，克服在单用轧制力分配时机架压下率波动大的弊端（此种情况在轧制力模型误差大和工艺条件变化大时发生，轧制力分配要通过调节压下率来达到各机架轧制力比例的要求，所以压下率波动大；而使用压下分配时各机架的压下率几乎保持不变），有利于提高提高精轧带钢轧制规程设定精度和轧制稳定性。

某带钢信息如下：中间坯厚度40.69mm，终轧厚度2.0mm，带钢宽度1233.6mm，精轧入口温度FET为1000℃，终轧目标温度FDT为880℃。精轧轧制道次数n=7。层别信息：钢种大类为11002，厚度等级为6，宽度等级为3，出钢记号为GV4924E1。从3个表格取到相应层别的数据如下：

表1：压下率分配系数表

表2：压下率学习表：

表3：轧制力分配系数表

上述步骤S1：根据压下分配系数表，按压下分配模式计算各机架的厚度分配将压下分配系数表和压下率学习表的数据相加，得到下表4：压下分配表格值

为了得到在中间坯厚度为40.69mm和终轧厚度为2.0mm下的具体的压下率数据，需要对进行相对化计算。

根据公式（1），计算压下率的比例因子xrk

$x=\frac{h_0}{h_7}=20.346,{\mathrm{xs}}_1=\underset{i=1}{\overset{7}{\mathrm{\&Sigma;}}}{\mathrm{\&epsiv;}}_0^i=2.36,{\mathrm{xs}}_2=\underset{i=1}{\overset{7}{\mathrm{\&Sigma;}}}{\left({\mathrm{\&epsiv;}}_0^i\right)}^2=0.9768$

$\mathrm{xrk}=\frac{-\mathrm{kb}\&CenterDot;{\mathrm{xs}}_1+\sqrt{{(\mathrm{kb}\&CenterDot;{\mathrm{xs}}_1)}^2+4\&CenterDot;\mathrm{ka}\&CenterDot;{\mathrm{xs}}_2\&CenterDot;1\mathrm{nx}}}{2\&CenterDot;\mathrm{ka}\&CenterDot;{\mathrm{xs}}_2}=0.984$

得到下表5：各机架相对化后的压下率及出口厚度数据

上述步骤S2：以压下分配结果为初始值，得到各机架的轧制力分配结果。

先执行绝对压下量迭代计算过程（S23～S27），得到下表6：各次迭代的轧制力比例变化

迭代次数	F1	F2	F3	F4	F5	F6	F7
								1	1	0.90617	0.85469	0.75396	0.55891	0.45492	0.31708
2	1	0.93258	0.86329	0.75989	0.5853	0.45067	0.33803
								3	1	0.95718	0.87182	0.76879	0.60745	0.44919	0.35346
4	1	0.97063	0.87697	0.7759	0.61707	0.44984	0.35862
								目标	1	0.98	0.88	0.78	0.62	0.45	0.36

再执行超限检查过程（S28）,各机架轧制力都没有超出其允许的最大轧制力，故不需要重新回到S23。

下表7：轧制力分配的压下率和出口厚度

上述步骤S3：对轧制力分配结果做限幅处理，使其偏离压下分配结果在一定范围内。

1、对轧制力分配结果做限幅处理。取η=6%，发现所有机架轧制力分配结果都在压下分配结果的±6%以内，故厚度分配不变。

2、因为上述1中数据没有变化，故无需相对化计算。

上述步骤S4－S6：带钢头部穿带完成并收集到实际轧制力后，判断轧制力设定值与实测值的偏差，如果各机架偏差都在θ以内，则对本次设定的压下率进行学习。

带钢头部穿带完成并收集到实际轧制力后，判断轧制力设定值与实测值的偏差，发现各机架偏差都在5%以内（取θ=5%），得到下表8：进行压下率学习

将更新到压下率学习表中的相应层别。由于学习的作用，该层别压下分配工艺表与轧制力分配工艺表数据会越来越吻合，即压下分配结果对应的轧制力会更加接近轧制力分配目标比例的要求。

这种情况最为理想，不但轧制力分配达到目标比例要求，而且偏离压下分配结果在一定范围内，无需限幅处理，最终分配取轧制力分配的结果。

下面假设对同样的实施例，输入数据不变，仅轧制力模型偏差增大为下表9：

相同的计算过程进行一遍。

步骤S1结果不变。由于轧制力模型的偏差不一样，步骤S2计算结果不同。迭代5次才满足终止条件，各次迭代的轧制力比例变化如下表10：

迭代次数	F1	F2	F3	F4	F5	F6	F7
								1	1	0.86302	0.7326	0.64625	0.5323	0.47658	0.27178
2	1	0.89134	0.77522	0.69524	0.57158	0.46032	0.3122
								3	1	0.92722	0.82495	0.74563	0.60658	0.45129	0.34476
4	1	0.95498	0.85953	0.77391	0.62107	0.45078	0.35705
								5	1	0.97054	0.87596	0.78306	0.62332	0.45095	0.3599
目标	1	0.98	0.88	0.78	0.62	0.45	0.36

轧制力分配的压下率和出口厚度如下表11：

上述步骤S3，对轧制力分配结果进行限幅处理，发现F1、F3、F4、F6、F7机架都超出压下分配结果的±6%（η=6%），限幅处理结果如下表12：

限幅处理后要做相对化计算。执行1次满足结束条件，xc=0.0079，最终输出的压下率和出口厚度如下表13：

最后，将前后两种情况的轧制力分配与最终分配的压下率波动做对比如下表14：

图3为轧制力分配和最终分配由轧制力模型偏差导致的压下率波动。可见，如果直接将轧制力分配作为最终分配，那么后者相对于前者的波动很大，F6、F7将近15%，即在其它条件都不变仅因轧制力模型偏差就可以导致这么大的波动，对轧制规程设定计算和轧制稳定性都不利；但限幅处理后的最终分配，后者相对于前者的波动并不大，基本控制在5%以内。这就是本专利提出要用压下分配确定的分布范围来约束轧制力分配结果的原因，以往技术要么是单纯的“压下分配”，要么是单纯的“轧制力分配”，本专利实现了“压下分配+轧制力分配”,它的优势在于：在轧制力模型偏差不大时，通过表格2的学习使得两种分配模式结果越来越一致，迭代次数越来越少；在轧制力模型偏差较大或工艺条件变化较大时，通过压下分配的限幅作用能有效控制各机架压下率的波动范围，有利于提高精轧规程设定计算精度和轧制稳定性。

实施例2：

机架设备极限检验。

带钢信息：中间坯厚度42.308mm，终轧厚度4.597mm，道次数n=7。步骤S2中要求的轧制力分配系数如下表15：

SGF	THICK	α1	α2	α3	α4	α5	α6	α7
									11002	16	1	0.98	0.88	0.78	0.62	0.45	0.36

轧制力超限检查

为便于比较，分别计算两种情况下的负荷分配：一种是对F2的最大轧制力没有限制；一种是对F2的最大轧制力进行限制，假设计算结果比较如下表16。

	F1	F2	F3	F4	F5	F6	F7
								αi	1	0.98	0.88	0.78	0.62	0.45	0.36
αi修改后	1	0.94	0.88	0.78	0.62	0.45	0.36

可见，在设置最大轧制力限制时，F2机架轧制力满足

综上所述，本发明针对现有技术中压下分配和轧制力分配的不足，提出一种新的精轧带钢负荷分配方法，实现“压下分配+轧制力分配”，用压下分配确定的分布范围来约束轧制力分配结果，综合两种模式的优点，同时克服两者各自的缺点，在减小各机架厚度分配波动的同时尽量满足各机架轧制力按比例下降的要求，达到提高精轧带钢轧制规程设定精度和轧制稳定性的目的。主要创新点说明：从单个步骤看，主要为步骤S2、3；从总体看，实现了两种分配模式的结合，优势：在轧制力模型偏差不大时，通过表格2的学习使得两种分配模式结果越来越一致，迭代次数越来越少；在轧制力模型偏差较大或工艺条件变化较大时，通过压下分配的限幅作用能有效控制各机架压下率的波动范围，有利于提高精轧规程设定计算精度和轧制稳定性。

本技术领域中的普通技术人员应当认识到，以上的实施例仅是用来说明本发明的目的，而并非用作对本发明的限定，只要在本发明的实质范围内，对以上所述实施例的变化、变型都将落在本发明的权利要求的范围内。

Claims (6)

1.一种精轧带钢负荷分配设定方法，其特征在于包括以下步骤：

S1、根据压下率分配系数表和压下率学习表，按压下分配模式计算精轧各机架的压下分配结果，所述压下分配结果指压下分配的压下率和压下分配的出口厚度

S2、将上述得到的各机架压下分配的出口厚度作为初始值，按轧制力分配模式迭代计算各机架的绝对压下量，直到轧制力比例满足目标要求，得到各机架的轧制力分配结果，所述轧制力分配结果包括轧制力分配的出口厚度及轧制力分配的压下率

S3、对所述各机架的轧制力分配结果进行限幅处理，使其偏离所述压下分配结果在一定范围内，从而使各机架的最终压下率在由压下分配确定的一定范围内，只要没超出范围就取轧制力分配结果，超出范围则将其限制在边界处；

S4，带钢头部穿带完成并收集到实际轧制力后，判断各机架轧制力设定值与实际轧制力的偏差；

S5，判定各机架偏差值是否都在θ以内，θ为偏差控制参数，取5%-10%，如果不满足各机架轧制力偏差都在θ以内的条件，则不进行压下率学习，此时，所述压下率学习表中相应层别的数据不进行更新，流程进入结束步骤；

S6，如果各机架偏差值在θ以内，则对本次设定的最终压下率进行学习，将学习值更新到压下率学习表格相应层别，经学习过程，该层别轧制力会逐步达到目标轧制力比例的要求；

S7，负荷分配设定结束。

2.如权利要求1所述的精轧带钢负荷分配设定方法，其特征在于：

所述步骤S1包括以下步骤：

S11，取表，从压下率分配系数表和压下率学习表取到各机架的压下率和压下率学习值，两者进行相加，得到各机架的压下率分配表格值

S12，进行相对化计算，得到各机架的压下分配结果：

根据精轧各机架的压下量之和必须等于精轧的总压下量，得到压下率的比例因子xrk为：

$\mathrm{xrk}=\frac{-\mathrm{kb}\&CenterDot;{\mathrm{xs}}_1+\sqrt{{(\mathrm{kb}\&CenterDot;{\mathrm{xs}}_1)}^2+4\&CenterDot;\mathrm{ka}\&CenterDot;{\mathrm{xs}}_2\&CenterDot;1\mathrm{nx}}}{2\&CenterDot;\mathrm{ka}\&CenterDot;{\mathrm{xs}}_2}$       公式1

其中，h_n为终轧厚度，h₀为中间坯厚度；

n为精轧道次数，i为机架号；

kb=0.906501，ka=0.959597。

各机架的压下率分配表格值相对化后的压下率为压下分配的压下率

${\mathrm{\&epsiv;}}_{\mathrm{red}}^i=\mathrm{xrk}\&CenterDot;{\mathrm{\&epsiv;}}_0^i$          公式2

根据压下分配的压下率得到前面n-1个机架各自的压下分配的出口厚度

$\begin{array}{cc}{h}_{\mathrm{red}}^i=h_0\&CenterDot;\underset{i=1}{\overset{i}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{red}}^i)& i=\mathrm{1,2},...,n-1.\end{array}$       公式3

3.如权利要求2所述的精轧带钢负荷分配设定方法，其特征在于：

所述步骤S11中各机架的压下率分配表格值取得过程为：

根据当前带钢的钢种大类、厚度等级和宽度等级信息，从压下分配系数表相应层别取到精轧各机架的压下率其中i为机架号；

根据当前带钢的出钢记号、厚度等级和宽度等级信息，从压下率学习表相应层别取到精轧各机架的压下率学习值其中i为机架号，所述压下率学习表是由程序自动学习得到的数据表，表中数据通过学习不断更新，初始值取0；

所述各机架的压下率分配表格值为：

4.如权利要求1所述的精轧带钢负荷分配设定方法，其特征在于：

所述步骤S2包括以下步骤：

S21，根据当前带钢的钢种大类和厚度等级信息，从轧制力分配系数表相应层别获取到各机架初始轧制力分配系数α_0i；

S22，进行操作修正处理，获得各机架修正后的轧制力分配系数α_1i：

从操作画面HMI收集各机架的操作修正因子ζ_i，各机架修正后的轧制力分配系数α_1i为：

α_1i=α_0i·(1+ζ_i)        公式4

其中，ζ_i取值范围为-20%～20%；

进行归一化，归一化后的轧制力分配系数α_i为：

S23，计算在当前厚度分配下各机架的轧制力和轧制力对绝对压下量的导数：

计算轧制力

设：i为机架号，j为迭代次数，为第i机架第j次迭代的轧制力，第i机架第j次迭代的出口厚度，第i-1机架第j次迭代的出口厚度；

第一次迭代计算前第i机架出口厚度初始值为上述公式3得到的压下分配的出口厚度

在其它工艺参数固定的情况下，轧制力是入口厚度和出口厚度的函数，存在： $p_i^j=p(h_{i-1}^j,h_i^j);$

计算轧制力对绝对压下量的导数

其中，为绝对压下量正扰动后的轧制力；为绝对压下量负扰动后的轧制力，为第i机架第j次迭代的绝对压下量；

存在以下函数关系： $p_{i1}^j=p(h_{i-1}^j,h_{i-1}^j-{\mathrm{\&Delta;h}}_{i1}^j),p_{i2}^j=p(h_{i-1}^j,h_{i-1}^j-{\mathrm{\&Delta;h}}_{i2}^j)$

其中，为正扰动后的绝对压下量，为负扰动后的绝对压下量，

${\mathrm{\&Delta;h}}_{i1}^j=(1+\frac{\mathrm{dt}}{2})\&CenterDot;{\mathrm{\&Delta;h}}_i^j,{\mathrm{\&Delta;h}}_{i2}^j=(1-\frac{\mathrm{dt}}{2}){\&CenterDot;\mathrm{\&Delta;h}}_i^j$

其中dt为扰动，取1%；

S24，根据上述各机架的轧制力和轧制力对绝对压下量导数的计算结果，计算各机架绝对压下量的修正值：

轧制分配迭代计算公式：

$\mathrm{\&delta;}\left({\mathrm{\&Delta;h}}_i^j\right)=({\mathrm{\&Sigma;p}}_i^j\&CenterDot;\frac{{\mathrm{\&alpha;}}_i}{\mathrm{\&Sigma;}{\mathrm{\&alpha;}}_i}-p_i^j)/\frac{{\&PartialD;p}_i^j}{\&PartialD;\left(\mathrm{\&Delta;}{h}_i^j\right)}-\mathrm{\&Sigma;}((\mathrm{\&Sigma;}{p}_i^j\&CenterDot;\frac{{\mathrm{\&alpha;}}_i}{\mathrm{\&Sigma;}{\mathrm{\&alpha;}}_i}-p_i^j)/\frac{{\&PartialD;p}_i^j}{\&PartialD;\left({\mathrm{\&Delta;h}}_i^j\right)})\&CenterDot;\frac{{\mathrm{\&Delta;h}}_i^j}{\mathrm{\&Sigma;\&Delta;}{h}_i^j}$       公式5

其中：为第i机架第j次迭代的绝对压下量的修正值；

S25，各机架绝对压下量的修正值得到后，计算各机架绝对压下量的更新值，该更新值为第j+1次迭代的绝对压下量

计算公式为：

$\mathrm{\&Delta;}{h}_i^{j+1}=\mathrm{\&Delta;}{h}_i^j+{\mathrm{damp}}_j\&CenterDot;\mathrm{\&delta;}\left(\mathrm{\&Delta;}{h}_i^j\right)$           公式6

其中damp_j为阻尼系数，damp_j=damp_mpy·β+(1-β)·(1-e^-j)，

damp_mpy取1.0，β取0.6；

S26，各机架绝对压下量的更新值得到后，计算各机架出口厚度的更新值，该更新值为第j+1次迭代的出口厚度

以及中间坯厚度h₀为已知量，根据仅第一个机架其余都是前一个机架的输出为后一个机架的输入，从第一个机架计算到最后一个机架，得到各机架出口厚度的更新值，这个值将用于下一次迭代计算轧制力用，其中指第i-1机架第j+1次迭代的出口厚度；

S27，判断轧制力比例是否满足收敛条件：

轧制力分配迭代计算的收敛条件为：

$|\frac{p_i^j}{\mathrm{\&Sigma;}{p}_i^j}-\frac{{\mathrm{\&alpha;}}_i}{\mathrm{\&Sigma;}{\mathrm{\&alpha;}}_i}|\&le;\frac{\mathrm{\&tau;}}{{\mathrm{\&Sigma;\&alpha;}}_i}$         公式7

其中，τ为正数，取0.01；

当式(7)成立，或者迭代次数j超过设定次数时结束迭代计算，否则继续执行所述S23步骤，同时迭代次数j累加1；

S28，上述步骤S27中如果轧制力比例满足收敛条件，则判断某机架的轧制力是否超过设备允许的最大轧制力：

如果存在轧制力超限：

当第i机架的轧制力p_i超过设备允许的最大轧制力时，进入S281步骤，调整该机架的轧制力分配系数，其余机架的分配系数不作调整，存在：

时，则执行 ${\mathrm{\&alpha;}}_i^{\&prime;}={\mathrm{\&alpha;}}_i\&CenterDot;\frac{1}{(p_i/p_i^{\max }-1)\&CenterDot;0.8+1}$     公式8

其中，α_i为所述归一化后的轧制力分配系数α_i，α_i’为超限调整后的轧制力分配系数；

α_i’得到后回到步骤S23，用调整后的的轧制力分配系数α_i’重新进行绝对压下量迭代计算，同时迭代次数j赋值0，取步骤S27最后一次迭代的绝对压下量和出口厚度作为绝对压下量和出口厚度的初始值；

如果不存在轧制力超限，则直接执行下一步骤；

S29，绝对压下量迭代计算和超限检查结束，得到轧制力分配结果，取步骤S27最后一次迭代的出口厚度作为各机架轧制力分配的出口厚度同时得到各机架轧制力分配的压下率为：式中为第i机架轧制力分配的出口厚度，为第i-1机架轧制力分配的出口厚度。

5.如权利要求1所述的精轧带钢负荷分配设定方法，其特征在于：

所述步骤3包括以下步骤：

S31，对轧制力分配结果进行限幅处理：

${\mathrm{\&epsiv;}}_{\mathrm{xf}}^i=\left\{\begin{array}{cc}{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})& {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i>{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})\\ {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i& {\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})\&le;{\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i\&le;{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1+\mathrm{\&eta;})\\ {\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})& {\mathrm{\&epsiv;}}_{\mathrm{rfd}}^i<{\mathrm{\&epsiv;}}_{\mathrm{red}}^i\&CenterDot;(1-\mathrm{\&eta;})\end{array}\right.$      公式9

其中为限幅处理后压下率，η为限幅控制参数，取5%～10%；

S32，对限幅处理后压下率进行相对化计算：

上述步骤S31中压下率如果被修改，计算终轧厚度h′_n与目标终轧厚度hn间会存在一定偏差，还要进行一轮相对化计算：

计算终轧厚度h′_n：

${h_n}^{\&prime;}=h_0\&CenterDot;\underset{i=1}{\overset{n}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i)$

计算终轧相对化因子(1-xc)，其中

$\mathrm{xc}=\frac{1n\frac{h_n}{{h_n}^{\&prime;}}}{\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}\frac{{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i}{1-{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i}}$

终轧相对化后的压下率称为最终分配压下率：

${\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i=(1-\mathrm{xc})\&CenterDot;{\mathrm{\&epsiv;}}_{\mathrm{xf}}^i$

取作为各机架的最终输出的压下率，同时得到各机架的最终出口厚度

$\begin{array}{cc}{h}_{\mathrm{xdh}}^i:h_{\mathrm{xdh}}^i=h_0\&CenterDot;\underset{i=1}{\overset{i}{\mathrm{\&Pi;}}}(1-{\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i)& i=\mathrm{1,2},...,n-1,\end{array}$

根据最终分配压下率和最终出口厚度进行精轧带钢的轧制规程设定计算。

6.如权利要求1所述的精轧带钢负荷分配设定方法，其特征在于：

所述步骤S6中，最终压下率进行学习的过程为：

求各机架最终分配与压下分配的压下率之差

${\mathrm{\&epsiv;}}_{\mathrm{dif}}^i=\frac{{\mathrm{\&epsiv;}}_{\mathrm{xdh}}^i-{\mathrm{\&epsiv;}}_{\mathrm{red}}^i}{\mathrm{xrk}}$         公式10

其中，为最终分配的压下率，为压下分配的压下率，xrk为公式（1）中的压下率的比例因子；

得到轧制规程设定后压下率学习表相应层别更新的学习值

${\mathrm{\&epsiv;}}_{\mathrm{lrn}}^{i,\mathrm{new}}={\mathrm{\&epsiv;}}_{\mathrm{lrn}}^i\&CenterDot;\mathrm{\&mu;}+{\mathrm{\&epsiv;}}_{\mathrm{dif}}^i\&CenterDot;(1-\mathrm{\&mu;})$         公式11

其中，μ为平滑系数，取0.7；为轧制规程设定前压下率学习表相应层别原来的学习值。