本发明提供一种确定单机架可逆冷轧机压下分配的方法，根据钢卷的原始数据、轧机的设备参数和工艺要求和给定的负荷分配比例系数，然后根据轧制工艺数学模型对轧制工艺参数进行计算，在计算过程中对设备能力进行极限校核，如果超限则对超限的工艺参数进行修正计算，对于不同的工艺条件分别采用不同的压下分配调整策略，进行迭代计算，直到满足算法的收敛条件为止，如果超过系统最大迭代次数或不满足收敛条件，则在画面上显示出错信息提示操作工，修改总道次数或者张力分配后重新进行计算，直到获得的满意的压下分配为止，最后将此压下分配保存到数据库中，为设定计算模型准备数据。
1.一种确定单机架可逆冷轧机压下分配的方法，其特征在于：它包括以下步骤：步骤1、输入带钢原始数据、设备参数、设备能力参数、模型系数和工艺要求；
步骤2、输入初始轧制总道次数：依据经验给定初始的轧制总道次数；
步骤3、给定各道次带钢出口厚度初值：将总压下率在各道次进行平均分配，计算得到的厚度分配作为迭代计算的初始值；
上式中，h[0,j]为第0次迭代时各道次轧制带钢出口厚度；j为道次号，j=1,2,3，…，n为总轧制道次数；hn为成品带钢厚度；h0为原料带钢厚度；
步骤4、确定带钢单位张力：查带钢单位张力表确定各道次带钢的单位前张力和单位后张力，如果带钢单位张力表中没有，则采用插值或人工输入确定；
步骤5、给定各道次轧制速度初值：确定各道次的轧制速度，取各道次轧制速度的最大值作为其轧制速度初值；
步骤6、检查开卷、卷取功率是否超限：对开卷、卷取功率进行校核，如果超限根据开卷、卷取额定功率的δ1反算轧制速度，并返回到步骤4，修改带钢单位张力后继续计算，直至开卷、卷取功率没有超限；δ1为比例系数，为经验值，取90%-98%；
步骤7、利用工艺数学模型计算轧制工艺参数，工艺参数包括材料变形抗力、各道次摩擦系数、前滑、轧制力、压扁半径；
步骤8、确定第i次迭代的最大轧制力和第i次迭代的最小轧制力所处的道次号；
上式中，Fmax[i]为第i次迭代的最大轧制力；Fmin[i]为第i次迭代的最小轧制力；
Fm[i]为第i次迭代的平均轧制力；ε为对数应变系数，由人为设定；
若不满足轧制力平衡条件则执行步骤10；若满足轧制力平衡条件，则执行步骤11；
步骤10、判断i是否大于最大迭代次数：如果i大于系统设定的最大迭代次数则结束，并输出本次计算的厚度分配；否则根据不同情况分别采用10a或10b调整压下分配：
10a、如果本次迭代的最小轧制力处于最大轧制力的右侧，则：将Fmin[i]所处道次的压下率增加压下率调节量dr作为下次迭代的压下率，其余道次的压下率保持不变；
采用倒推的方法计算各道次的入口厚度，即根据出口厚度和压下率计算入口厚度：h[i+1,k]=h[i+1,k+1]/(1-r[i+1,k+1])，上式中k为Fmax与Fmin之间的道次，h[i+1,k]为第i+1次迭代时，第k道次的入口厚度，h[i+1,k+1]为第i+1次迭代时，第k道次的出口厚度，r[i+1,k+1]为第i+1次迭代时，第k+1道次的压下率；
10b、如果本次迭代的最小轧制力处于最大轧制力的左侧，则：将Fmin[i]所处道次的压下率增加压下率调节量dr作为下次迭代的压下率，其余道次的压下率保持不变；
采用顺推的方法计算各道次的出口厚度，即根据入口厚度和压下率计算出口厚度：h[i+1,k+1]=h[i+1,k]*(1-r[i+1,k+1]);上式中k为Fmax与Fmin之间的道次，h[i+1,k+1]为第i+1次迭代时，第k道次的出口厚度，h[i+1,k]为第i+1次迭代时，第k道次的入口厚度，r[i+1,k+1]为第i+1次迭代时，第k+1道次的压下率；
然后根据厚度，重新计算道次压下率；从而计算相关工艺参数；
步骤11、利用工艺数学模型计算轧制力矩、电机转矩、电机功率；对主电机功率进行校核，检查主电机功率是否超限：
11a、若主电机功率超限，则采用工艺数学模型，用主电机额定功率的δ2，反算各道次的轧制速度，并返回步骤5；δ2为比例系数，小于δ1；
11b、若主电机功率未超限，则输出压下分配及轧制工艺参数，将计算得到的厚度分配保存到数据库中；
所述的步骤11采用如下工艺数学模型计算轧制力矩、电机转矩、电机功率：
上式中：ε为对数应变系数，k0、m、n为模型系数，kp为带钢的平均变形抗力，ε＝ln(H1/hm)，其中H1为热轧后带钢来料厚度，hm为平均带钢厚度，hm=(1-β)×H+β×h，其中H为入口带钢厚度，h为出口带钢厚度，β为加权系数；
上式中：μ为摩擦系数，vr为轧制速度，m0、m1、m2为模型系数；
上式中：r为压下率( )，μ为摩擦系数，R为工作辊半径，H为入口带钢厚度，h为出口带钢厚度，b为带钢宽度，kp为平均变形抗力，nt为张力影响系数，CH为压扁系数，A1,A2,A3,A4,A5,A6为中间变量，F为轧制力；
其中，tb为单位后张力，tf为单位前张力，kp为平均变形抗力；
上式中：R’为工作辊压扁半径，R为工作辊半径，CH为压扁系数，F为轧制力，b为带钢宽度，H为入口带钢厚度，h为出口带钢厚度；
上式中：R’为工作辊压扁半径，r为压下率，H为入口带钢厚度，h为出口带钢厚度，tb为单位后张力，tf为单位前张力，kp为平均变形抗力，μ为摩擦系数，α为模型系数；
上式中，R为工作辊半径，R为工作辊压扁半径，r为压下率，H为入口带钢厚度，h为出口带钢厚度，b为带钢宽度，kp为平均变形抗力，μ为摩擦系数；
上式中R为工作辊半径，tb为单位后张力，tf为单位前张力，H为入口带钢厚度，h为出口带钢厚度，b为带钢宽度；
上式中，KL为模型系数，F为轧制力，V为轧制速度；
HP=0.2192×10 ·(v/R)·GM/1.34上式中，v为轧制速度，R为工作辊半径，GM为电机转矩。
2.根据权利要求1所述的确定单机架可逆冷轧机压下分配的方法，其特征在于：所述的带钢原始数据包括钢种、来料厚度、成品厚度和宽度；所述的设备参数包括工作辊直径和杨氏弹性模量；设备能力参数包括最大轧制力、最大轧制速度、主电机额定功率、开卷机额定功率和卷取机额定功率；工艺要求包括各道次轧制力平衡负荷分配比例系数。
3.根据权利要求1所述的确定单机架可逆冷轧机压下分配的方法，其特征在于：所述的对数应变系数ε取0.05或者0.1。
4.根据权利要求1所述的确定单机架可逆冷轧机压下分配的方法，其特征在于：所述的压下率调节量dr取0.001。
一种确定单机架可逆冷轧机压下分配的方法
技术领域
[0001] 本发明涉及冶金自动化领域，尤其涉及一种确定单机架可逆冷轧机压下分配的方法。
背景技术
[0002] 压下分配是冷轧轧制中轧机设定计算的重要组成部分，它是根据轧件特性、轧机特性和工艺要求计算出特定意义下最优的厚度分配。合理的压下分配既可以有效地利用现有设备，充分发挥设备潜力，提高轧机的产品质量，同时还能降低轧制能耗。因此，寻找一种合理的压下分配计算方法对于实际生产有着重要的现实意义。
[0003] 压下分配最主要功能是在钢种、轧制道次、张力、入口厚度、出口厚度已知的情况下，确定各道次的压下率，同时得到相关的轧制工艺参数。单机架可逆冷轧机压下分配计算方法可以分为三个发展时期：以能耗曲线为代表的经验分配法；以轧制负荷比例关系为代表的轧制理论法；以目标优化函数为代表的优化计算方法。
[0004] 早期生产按照经验进行压下分配，各种规格带钢的压下分配按经验表格形式存储，表格的经验通常以能耗曲线为基础，轧制负荷分配以最小能耗消耗为目标。根据待轧钢卷的钢种、规格从表格中获取经验数据。这种方法没有考虑带钢材料特性的差异和轧制生产条件的变化对控制参数的影响，同时，由于本方法需要以前类似轧机的生产经验数据，它不适应现代冷轧机生产的要求。
[0005] 随着计算机控制技术的发展，以轧制工艺参数计算模型为基础的轧制理论用于在线控制。以各道次轧制负荷相关关系为目标，通过轧制模型在线计算的方法获得各道次的压下率。
[0006] 随着冷轧机自动控制技术的发展，提出以单一或多个工艺参数为优化目标，并建立起相应的约束条件，在寻优过程中确定轧机各道次的带钢厚度。计算过程中不同的负荷分配方法使用了动态规划、人工智能等优化计算方法。使用这种方法目标函数和约束条件的结构复杂，并且单纯从理论上构建的目标优化函数和约束条件与轧机实际工况有较大差异。
[0007] 另外，对于在线过程控制而言，轧制压下分配目标函数的建立可以由生产操作工进行手动干预，因此要求目标函数和约束条件的建立及其物理含义尽可能简单明了，操作工在轧制过程中可以根据实际生产情况对各道次的压下量的分配进行快速干预，但这种方法的使用依赖于操作工的操作水平和经验。
发明内容
[0008] 本发明所要解决的技术问题是：提供一种确定单机架可逆冷轧机压下分配的方法，能够满足不同钢种、不同规格的产品的轧制工艺要求。
[0009] 本发明为解决上述技术问题所采取的技术方案为：一种确定单机架可逆冷轧机压下分配的方法，其特征在于：它包括以下步骤：
[0010] 步骤1、输入带钢原始数据、设备参数、设备能力参数、模型系数和工艺要求；
[0011] 步骤2、输入初始轧制总道次数：依据经验给定初始的轧制总道次数；
[0012] 步骤3、给定各道次带钢出口厚度初值：将总压下率在各道次进行平均分配，计算得到的厚度分配作为迭代计算的初始值；
[0013] 步骤4、确定带钢单位张力：查带钢单位张力表确定各道次带钢的单位前张力和单位后张力，如果此表中没有，则采用插值或人工输入确定；
[0014] 步骤5、给定各道次轧制速度初值：确定各道次的轧制速度，取各道次轧制速度的最大值作为其轧制速度初值；
[0015] 步骤6、检查开卷、卷取功率是否超限：对开卷、卷取功率进行校核，如果超限根据开卷、卷取额定功率的δ1反算轧制速度，并返回到步骤4，修改带钢单位张力后继续计算，直至开卷、卷取功率没有超限；δ1为比例系数，为经验值，取90%-98%；
[0016] 步骤7、利用工艺数学模型计算轧制工艺参数，工艺参数包括材料变形抗力、各道次摩擦系数、前滑、轧制力、压扁半径；
[0017] 步骤8、确定第i次迭代的最大轧制力和第i次迭代的最小轧制力所处的道次号；
[0018] 步骤9、检查是否满足轧制力平衡条件：
[0019]
[0020] 上式中，Fmax[i]为第i次迭代的最大轧制力；Fmin[i]为第i次迭代的最小轧制力；Fm[i]为第i次迭代的平均轧制力；ε为对数应变系数，由人为设定；
[0021] 若不满足轧制力平衡条件则执行步骤10；若满足轧制力平衡条件，则执行步骤11；
[0022] 步骤10、判断i是否大于最大迭代次数：如果i大于系统设定的最大迭代次数则结束，并输出本次计算的厚度分配；否则根据不同情况分别采用10a或10b调整压下分配：
[0023] 10a、如果本次迭代的最小轧制力处于最大轧制力的右侧，则：
[0024] 将Fmin[i]所处道次的压下率增加压下率调节量dr作为下次迭代的压下率，其余道次的压下率保持不变；
[0025] 采用倒推的方法计算各道次的入口厚度，即根据出口厚度和压下率计算入口厚度；然后根据厚度，重新计算道次压下率；
[0026] 返回步骤4继续计算；
[0027] 10b、如果本次迭代的最小轧制力处于最大轧制力的左侧，则：
[0028] 将Fmin[i]所处道次的压下率增加压下率调节量dr作为下次迭代的压下率，其余道次的压下率保持不变；
[0029] 采用顺推的方法计算各道次的出口厚度，即根据入口厚度和压下率计算出口厚度；然后根据厚度，重新计算道次压下率；从而计算相关工艺参数；
[0030] 返回步骤4继续计算；
[0031] 步骤11、利用工艺数学模型计算轧制力矩、电机转矩、电机功率；对主电机功率进行校核，检查主电机功率是否超限：
[0032] 11a、若主电机功率超限，则采用工艺数学模型，用主电机额定功率的δ2，反算各道次的轧制速度，并返回步骤5；δ2为比例系数，小于δ1；
[0033] 11b、若主电机功率未超限，则输出压下分配及轧制工艺参数，将计算得到的厚度分配保存到数据库中。
[0034] 按上述方案，所述的带钢原始数据包括钢种、来料厚度、成品厚度和宽度；所述的设备参数包括工作辊直径和杨氏弹性模量；设备能力参数包括最大轧制力、最大轧制速度、主电机额定功率、开卷机额定功率和卷取机额定功率；工艺要求包括各道次轧制力平衡负荷分配比例系数。
[0035] 按上述方案，所述的对数应变系数ε取0.05或者0.1。
[0036] 按上述方案，所述的压下率调节量dr取0.001。
[0037] 本发明的工作原理为：根据钢卷的原始数据、轧机的设备参数和工艺要求和给定的负荷分配比例系数，然后根据轧制工艺数学模型对轧制工艺参数进行计算，在计算过程中对设备能力进行极限校核，如果超限则对超限的工艺参数进行修正计算，对于不同的工艺条件分别采用不同的压下分配调整策略，进行迭代计算，直到满足算法的收敛条件为止，如果超过系统最大迭代次数或不满足收敛条件，则在画面上显示出错信息提示操作工，修改总道次数或者张力分配后重新进行计算，直到获得的满意的压下分配为止，最后将此压下分配保存到数据库中，为设定计算模型准备数据。
[0038] 本发明的有益效果为：
[0039] 1、本方法对不同的工艺条件分别采用不同的压下分配调整策略，进行迭代计算，原理简单、计算可靠，避免了传统技术中复杂的函数处理；本发明操作简便、实现容易，无须依赖操作工的操作水平和经验，适用于生产现场在线控制。
[0040] 2、本发明具有推广应用价值，可推广应用于单机架两辊、四辊可逆冷轧机以及多机架冷连轧机的压下分配及过程控制中。
附图说明
[0041] 图1为本发明的工作流程图。
具体实施方式
[0042] 基于本发明的确定单机架可逆冷轧机压下分配的方法可适用于两辊、四辊、六辊单机架冷轧机。本实施例公开的是某单机架六辊可逆冷轧机过程控制系统中的确定单机架可逆冷轧机压下分配的方法。该六辊轧机可轧制的产品包括普碳钢、高强钢、部分不锈钢和硅钢等。本实施例轧制的是中高牌号硅钢，机型为UCM轧机。
[0043] 本实施例提供的确定单机架可逆冷轧机压下分配的方法如图1所示，包括以下步骤：
[0044] 步骤1、输入带钢原始数据、设备参数、设备能力参数、模型系数和工艺要求。
[0045] 带钢原始数据包括钢种、来料厚度、成品厚度和宽度；所述的设备参数包括工作辊直径和杨氏弹性模量；设备能力参数包括最大轧制力、最大轧制速度、主电机额定功率、开卷机额定功率和卷取机额定功率；工艺要求包括各道次轧制力平衡负荷分配比例系数。
[0046] 带钢原始数据：以钢种HNGO（高牌号无取向硅钢）为例，来料厚度h0=2.3mm，成品厚度hn=0.5，宽度b=1270mm；设备参数：工作辊泊松系数=0.3，工作辊杨氏弹性2
模量=20600kg/mm ；设备能力参数：轧机最大轧制力Fmax=18000kN，主电机额定功率Nmax=5500kW，开卷机额定功率Nuncoilermax=650kW，卷取机额定功率Ncoilermax=2640kW，工作辊直径wr_diam=340mm，轧机最大出口速度vmax=900m/min，开卷机最大速度vmaxuncoiler=500m/min，卷取机最大速度vmaxcoiler=950m/min，轧制力平衡负荷分配比例系数=1:1:1:1:1。
[0047] 步骤2、输入初始轧制总道次数：依据经验给定初始的轧制总道次数，PASSCOUNT＝5。
[0048] 步骤3、给定各道次带钢出口厚度初值：将总压下率在各道次进行平均分配，计算得到的厚度分配作为迭代计算的初始值，计算公式如下：
[0049]
[0050] 上式中，h[0，j]为第0次迭代时各道次轧制带钢出口厚度；j为道次号，j=1,2,3，…，n为PASSCOUNT总道次数，PASSCOUNT=5（以轧制道次=5为例）；hn为成品带钢厚度；h0为原料带钢厚度。
[0051] 步骤4、确定带钢单位张力：查带钢单位张力表确定各道次带钢的单位前张力和单位后张力，如果表中没有，则采用插值或人工输入确定。
[0052] 在本实例中查带钢单位张力表得到的单位张力如下表所示：
[0053]
[0054] 步骤5、给定各道次轧制速度初值：确定各道次的轧制速度，取各道次轧制速度的最大值Vmax[j]作为其轧制速度初值Vr[0，j]，j为道次号，j=1～5。
[0055] 步骤6、检查开卷、卷取功率是否超限：对开卷、卷取功率进行校核，如果超限根据开卷、卷取额定功率的96%反算轧制速度，并返回到步骤4，修改带钢单位张力后继续计算，直至开卷、卷取功率没有超限。
[0056] 步骤7、利用工艺数学模型计算轧制工艺参数，工艺参数包括材料变形抗力、各道次摩擦系数、前滑、轧制力、压扁半径。
[0057] 步骤8、确定第i次迭代的最大轧制力和第i次迭代的最小轧制力所处的道次号。
[0058] 步骤9、检查是否满足轧制力平衡条件：
[0059]
[0060] 上式中，Fmax[i]为第i次迭代的最大轧制力；Fmin[i]为第i次迭代的最小轧制力；Fm[i]为第i次迭代的平均轧制力；ε为对数应变系数，由人为设定，一般取0.05或者0.1；
[0061] 若不满足轧制力平衡条件则执行步骤10；若满足轧制力平衡条件，则执行步骤11。
[0062] 步骤10、判断i是否大于最大迭代次数：如果i大于系统设定的最大迭代次数则结束，并输出本次计算的厚度分配；否则根据不同情况分别采用10a或10b调整压下分配：
[0063] 10a、如果本次迭代的最小轧制力处于最大轧制力的右侧，则：
[0064] 将Fmin[i]所处道次的压下率增加压下率调节量dr作为下次迭代的压下率，其余道次的压下率保持不变；压下率调节量dr一般取0.001。
[0065] 采用倒推的方法计算各道次的入口厚度，即根据出口厚度和压下率计算入口厚度：
[0066] h[i+1,k]=h[i+1,k+1]/(1-r[i+1,k+1])，
[0067] 上式中k为Fmax与Fmin之间的道次，h[i+1,k]为第i+1次迭代时，第k道次的入口厚度，h[i+1,k+1]为第i+1次迭代时，第k道次的出口厚度，r[i+1,k+1]为第i+1次迭代时，第k+1道次的压下率；
[0068] 然后根据厚度，重新计算道次压下率；
[0069] 返回步骤4继续计算；
[0070] 10b、如果本次迭代的最小轧制力处于最大轧制力的左侧，则：
[0071] 将Fmin[i]所处道次的压下率增加压下率调节量dr作为下次迭代的压下率，其余道次的压下率保持不变；压下率调节量dr一般取0.001。
[0072] 采用顺推的方法计算各道次的出口厚度，即根据入口厚度和压下率计算出口厚度：
[0073] h[i+1,k+1]=h[i+1,k]*(1-r[i+1,k+1]);
[0074] 上式中k为Fmax与Fmin之间的道次，h[i+1,k+1]为第i+1次迭代时，第k道次的出口厚度，h[i+1,k]为第i+1次迭代时，第k道次的入口厚度，r[i+1,k+1]为第i+1次迭代时，第k+1道次的压下率；
[0075] 然后根据厚度，重新计算道次压下率；从而计算相关工艺参数；
[0076] 返回步骤4继续计算。
[0077] 步骤11、利用工艺数学模型计算轧制力矩、电机转矩、电机功率；对主电机功率进行校核，检查主电机功率是否超限：
[0078] 11a、若主电机功率超限，则采用工艺数学模型，用主电机额定功率的93%，反算各道次的轧制速度，并返回步骤5；
[0079] 11b、若主电机功率未超限，则输出压下分配及轧制工艺参数，将计算得到的厚度分配保存到数据库中。
[0080] 下表为本实施例的计算结果，由下表可知，采用本方法能够得到在满足轧制力平衡的条件下，同时满足在设备能力允许范围内的满意的压下分配。
[0081]道次号 出口厚度(mm) 轧制力(kN) 电机功率(kW) 压扁半径(mm)
1 1.813 12125.4 3021.8 241.15
2 1.335 12717.7 5452.8 246.40
3 0.977 12412.2 5461.3 269.40
4 0.697 12085.0 5450.0 293.83
5 0.5 12088.1 5459.3 346.37
[0082] 本发明所采用的主要工艺数学模型：
[0083] 1．变形抗力模型
[0084] kp＝kO×(ε+m)n
[0085] 上式中：ε为对数应变系数，k0、m、n为模型系数，kp为带钢的平均变形抗力。
[0086] ε＝ln(H1/hm)，其中H1为热轧后带钢来料厚度，hm为平均带钢厚度。
[0087] hm＝(1-β)×H+β×h，其中H为入口带钢厚度，h为出口带钢厚度，β为加权系数=0.75。
[0088] 2．摩擦模型
[0089]
[0090] 上式中：μ为摩擦系数，vr为轧制速度，m0、m1、m2为模型系数。
[0091] 3.轧制力模型
[0092]
[0093]
[0094] A3＝1.08-1.02×γ
[0095]
[0096]
[0097]
[0098] F＝b×(A4+A6)
[0099] 上式中：r为压下率
[0100] μ为摩擦系数，R为工作辊半径，H为入口带钢厚度，h
[0101] 为出口带钢厚度，b为带钢宽度，kp为平均变形抗力，nt为张力影响系数，CH为压扁系数，A1，A2,A3,A4，A5，A6为中间变量，F为轧制力。
[0102]
[0103] 其中，tb为单位后张力，tf为单位前张力，kp为平均变形抗力。
[0104] 4.工作辊压扁半径模型
[0105]
[0106] 上式中：R’为工作辊压扁半径，R为工作辊半径，CH为压扁系数，F为轧制力，b为带钢宽度，H为入口带钢厚度，h为出口带钢厚度
[0107] 5.中性角模型
[0108]
[0109] 上式中：R’为工作辊压扁半径，r为压下率，H为入口带钢厚度，h为出口带钢厚度，tb为单位后张力，tf为单位前张力，kp为平均变形抗力，μ为摩擦系数，α为模型系数。
[0110] 6.前滑模型
[0111]
[0112] f=0:φ<0
[0113] Ф为中性角。
[0114] 7.轧制力矩、张力力矩、损失力矩模型
[0115] 轧制力矩GR模型：
[0116] GR＝b×kp×R×(H-h)×DG
[0117]
[0118] 上式中，R为工作辊半径，R为工作辊压扁半径，r为压下率，H为入口带钢厚度，h为出口带钢厚度，b为带钢宽度，kp为平均变形抗力，μ为摩擦系数。
[0119] 张力力矩GT模型
[0120] GT＝R×b×(tb×H-tf×h)
[0121] 上式中R为工作辊半径，tb为单位后张力，tf为单位前张力，H为入口带钢厚度，h为出口带钢厚度，b为带钢宽度。
[0122] 损失力矩GL模型
[0123]
[0124] 上式中，KL为模型系数，F为轧制力，V为轧制速度。
[0125] 8.电机转矩模型
[0126] GM＝GR+GT+GL
[0127] 9.电机功率模型
[0128] HP=0.2192×10-3·(v/R)·GM/1.34
[0129] 上式中，v为轧制速度，R为工作辊半径，GM为电机转矩。
[0130] 以上实施例仅用于说明本发明的计算思想和特点，其目的在于使本领域内的技术人员能够了解本发明的内容并据以实施，本发明的保护范围不限于上述实施例，凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰，均在本发明的保护范围之内。
