热轧带钢负荷分配目标厚度函数法.docVIP

第27卷第7期2006年7月东北大学学报(自然科学版)JournalofNortheasternUniversity(NaturalScience)Vol127,No.7Jul.2006 文章编号:100523026(2006)0720768203 热轧带钢负荷分配目标厚度函数法 江　潇1,胡贤磊1,刘相华1,郭立合2 (1.东北大学轧制技术及连轧自动化国家重点实验室,辽宁沈阳　110004; 2.本溪钢铁集团热连轧厂,辽宁本溪　117000) 摘　　　要:在分析了传统的解非线性方程组负荷分配方法缺陷的基础上提出目标厚度函数法? 该方法首先建立总负荷与目标出口厚度的单调函数关系,而后给定目标厚度并对n个一维非线性方程求根得出各道次压下规程,避免了初始分配率限定问题,同时也对迭代过程中可能出现的轧制道次出口厚度大于入口厚度的异常情况做了处理?使用该方法并采用实际生产数据进行优化计算,结果表明该方法运算精确快捷,可以满足实时在线控制的要求?关　键　词:热轧;带钢;负荷分配;非线性方程组;目标厚度函数法中图分类号:TG33515　　　文献标识码:A 板带轧制过程控制的目的是根据轧件、轧机 厚度和负荷分配,据,、品种都准确预测能耗与压下量间定量关系;负荷分配系数方法首先确定各机架/道次轧制力或功率的比例系数[2～7],而后由此计算出各道次的厚度,文献[8,9]中通过建立非线性方程组求解,本文在分析其方法缺陷的基础上,提出目标厚度函数法,实现负荷成比例的优化目标? i j=1 6? j n (2) Pαi - =0, αi+1 P(3) 式中,Pi为第i机架负荷(轧制力或电机功率);αi为第i机架目标负荷分配比? 由式(3)可建立方程组: α1 -- =0, α2=0,α3 (4) α2 1　解非线形方程组方法 设来料厚度h0,目标厚度ht,机架数n,目标 负荷分配比αi,需计算各机架出口带钢厚度hi(i=1,2,…,n-1)? 　　… αn-1 - αn=0式(4)中共包括n-1个方程?根据轧钢工艺理论可知任意一个机架i(i=1,2,…,n)的负荷Pi是该机架入口厚度hi-1和出口厚度hi的函数,式 (4)中的未知量即为机架间带钢厚度(h1～hn-1),共n-1个?方程数与未知量数相等,因此 实际应用中,通常综合考虑层别表中对应该规格产品的标准负荷分配比和人工干预量来确定目标负荷分配比 3 α(1)i=αi?ci,式中,αi为第i机架目标负荷分配比;αi为第i机架标准负荷分配比(层别数据);ci为第i机架负荷分配比修正系数? 3 可通过解此方程组实现目标负荷分配比? 文献[8]使用Newton2Raphson 法求解,需先给定一个严格的迭代初始值以确保收敛,并且当 迭代过程中出现上游机架轧件厚度小于下游机架 收稿日期:2005209209 基金项目:国家自然科学基金资助项目?作者简介:江　潇(1980-),男,浙江衢州人,东北大学博士研究生;刘相华(1953-),男,黑龙江双鸭山人,东北大学教授,博士生 导师? ? 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved. 第7期　　　　　　　　江　潇等: 热轧带钢负荷分配目标厚度函数法轧件厚度(hi-1hi)时将无法计算轧制负荷Pi;文献[9]使用平衡迭代法求解,不易给定修正系数mi,亦无法完全避免出现hi-1hi的情况?为此本文提出了目标厚度函数法,以避免求解非线形方程组? 769 求解过程亦可采用直线搜索方法?当求解得到的Pt超出轧机负荷能力上限或终轧温度超限时,需调整穿带速度,重新计算? 上述整个求解过程即为目标厚度函数法,其流程如图1所示 ? 2　目标厚度函数法 对于某单个机架而言,在入口厚度hi-1固定的情况下,轧制负荷P是出口厚度hi的单调减函数,即有P=P(hi)?在迭代计算中若变形程度取值过大,屈服极限随变形程度加大有可能降低,从而破坏轧制力对出口厚度的单调性?在本算法中将变形程度取值限定在一定范围内,既符合现场实际,也保持了函数的局部单调性?在函数单调性成立的前提下,就可以通过给定某个负荷值P,使用直线搜索方法求出满足min[P-P(h0)]的 h0? 对整个轧制过程而言,首先给定某个总轧制负荷值Pt,按照目标负荷分配比αi分配各机架负荷,即 ii) 图1　算法流程图 Fig.1　Flowchartofalgorithm ,,P1P1(h1)得 h1,P2P2(h2)得h2,如此一直到最后求解Pn=Pn(hn)得末机架出口厚度hn,显然hn是Pt的函数,称其为目标厚度函数: hn=f(Pt)? (6) 3　计算实例