D0I:10.13374/i.is8n1001-053x.1992.01.028 第14卷第1期 北京科技大学学报 Vo1.14No.1 1992年1月 Journal of University of Science and Technology Beijing Jan。1992 棒材轧后控冷计算机控制系统 方安祥·蔡庆伍·宋佩莼… 摘要:钢材轧后控冷中引入计算机控制系统,可以根据终轧温度调节冷却水量以达到 预定自回火蓝度,保证产品质量的意定。唐钢小型厂棒材轧后控冷计算构控制系统的分析、 硬件与软件系统的构成、数学模型的推导和自学习的实现,并给出了控制结果。 关健词:控冷,计算机控制系统,棒材 Computer Control System for the Controlled Cooling of Rod after Rolling Fang Anxiang Cai Qingwu.Shong Peichun. ABSTRACT:To assure the stableness of quality of the products,the computer control system adjusts the quantity of water,depending on the temperature of rod after rolling,in order to reach the predetermined temperature of rod's self-temper.This paper introduced the analysis of this computer control system which has been used in Tangshan Iron and steel mill,constitution of the ha- rdware and software of this system,the inferring of mathematical model and the self-learning of model.The results of controlling were also given. KEY WORDS:controlled cooling,computer control system,rod 1991-03-02收到初稿,1991-09-20稍到修改稿 ,计算机系(Department of Computer Science). ··金属压力加工系(Department of Metals Forming) 70
第 卷第 期 年 月 北 京 科 技 大 学 学 报 。 。 。 棒材轧后控冷计算机控制系统 方安祥 蔡庆伍” 宋佩药“ 摘 要 舰材轧后控冷中引人计算机控制系统 , 可以根据 终轧 温 度调节冷却水量以达到 预定 自回火沮度 , 保证产品质量的稳定 。 唐钢小型厂棒材轧后控冷计算构控制系统 的分析 、 硬件与软件系统 的构成 、 数学模型 的推导和 自学习 的实现 , 并给出 了控制结果 。 关健饲 控冷 , 计算机控制系统 , 棒材 侧 ” 劣 ‘ 夕 ” 口“ , 夕 二 , , , , 一 扭 。 五 , 五 , 一 。 , , ,, 一 一 收到初稿 , ,, 一 一 稿到修改稿 计算机系 几 八 金属压力加工系 丹子 DOI :10.13374/j .issn1001-053x.1992.01.028
棒材轧后控冷是提高产品质量、性能,节省能源的新工艺技术。具体方法是依据轧制规 格、实测终轧温度和目标自回火福温度,进行水量设定,对轧后棒材冷却,以达到规定自回火 温度范围。国外从60年代开始研究并在生产中普遍采用。国内从1972年开始研究,应用尚不 普遍。已采用该工艺的系统仅为人工操作,不能及时依据终轧温度调节水量,造成自回火温 度波动较大,产品质量不稳定。 唐钢小型厂建立的棒材轧后控冷生产线,在国内首次引入计算机控制系统,保证了实时 准确地控制,并通过自学习克服生产过程中存在的一些环境、特性等变化因素的影响,使系 统工作处在最佳状态。 1系统的分析与构成 1.1系统的分析 水量调节阀及流量计反应时间使得水流量调节变化较为迟缓。经现场实测,总水流量每 秒平均调节能力8一10t/h,精轧出口轧件速度为8m/s。轧件测量终轧温度后至进入水冷器的 时间约为0.8s。因此,依据精轧出口实测轧件温度进行水量模型计算并调节水量,在时间上是 不允许的。为此,在精轧入口设置测温点,测量轧件温度,经终轧温度预报模型计算,预报 该轧件在精轧出口温度,并以此进行水量的计算与设定。在精轧出口对该轧件实测终轧温 度,进行终轧温度预报模型的自学习,以不断提高预报精度。 轧件在穿水冷却后的终冷表面温度T。与自回火温度T:的关系及变化情况示于图1。考 虑到唐钢小型厂在冷床设置测温点的困难,测量不易淮确等问题,把水冷后测温点设置在冷 床入口前端。以该点实测温度做为终冷温度,建立该温度与自回火温度的关系模型,并以终 冷温度代替自回火温度值。利用实测终冷温度对水量设定模型进行自学习,提高控制精度。 Tst Water-pipe Operation quide PID Computer 00 Regulator Full-auto Time/s 图1自回火示意图 图2系统结构示意 Fig.1 Diagrammatic sketch of self- Fig.2 Diagrammatic sketch of system temper of steel constitution 1.2系统结构与工作方式 系统结构如图2,采用2种工作方式: (1)计算机全自动控制方式计算机依据输入数据,经模型计算输出水量设定值,由流 量调节仪表进行水量控制。计算机利用自学习,保证系统工作在最佳状态”。 71
棒 材轧 后控冷是 提高产品质量 、 性能 , 节省能源的新工艺技 术 。 具 体方法是依据轧制规 格 、 实测终轧温度和 目标 自回火温度 , 进行水量设定 , 对轧后棒 材冷却 , 以达到规定 自回火 温度范 围 。 国外 从 年代开始研究并在生 产中普遍采 用 。 国内从 年开始研究 , 应 用尚不 普遍 。 已采用该工艺的 系统仅为 人 工操作 , 不能及时依 据终轧温度调 节水量 , 造 成 自回火温 度波动较大 , 产品质量不稳定 。 唐钢 小型厂建 立的棒材轧后控冷生 产线 , 在 国内首 次 引入 计算机控制 系统 , 保证了实时 准确地控制 , 并通过 自学 习克服生 产过程 中存在 的 一些环 境 、 特性等变化 因素的影响 , 使系 统 工作处在 最佳状态 。 系统的分析与构成 系统 的分析 水量调 节阀及 流 量计反应时间使得 水流 量调 节变化较为迟 缓 。 经 现场 实测 , 总 水流量每 秒 平均 调 节能 力 一 , 精轧 出 口 轧 件速 度为 。 轧 件测量终轧温 度后至进入 水冷器的 时 间约为 。 ‘ 。 因此 ,依据 精轧 出 口 实测轧件 温度进行 水量模型 计算并 调 节水量 , 在 时 间上是 不 允许 的 。 为此 , 在 精轧 入 口 设置 测温点 , 测量轧件温度 , 经 终轧温 度预报模型计算 , 预报 该轧件在 精轧 出 口 温度 , 并以 此进行 水量的 计算与 设定 。 在 精轧 出 口 对 该轧件 实 测终轧温 度 , 进行终轧温度预报模型的 自学 习 , 以不 断 提高预报精度 。 轧件在 穿水冷却后的终冷表 面温度 。 。 与 自回火 温度 , 的关 系及变化情况示于图 。 考 虑到 唐钢 小型厂在 冷床设置测温点的 困难 , 测量不 易准确等问题 , 把水冷后测温点设置在 冷 床入 口前端 。 以该点实测温度做为 终冷温度 , 建 立该温度与 自回火温度 的关系模型 , 并以终 冷温度代替 自回火温度值 。 利 用实测 终冷温度对 水量设定模型进行 自学习 , 提高控制精度 。 一 弓 尸 。 口 。 ‘洲尸。 一 一一一 一 一 占一 一 〔〕 下三 厂 图 自回火示 意图 七 尺 年。门自 口卜。﹄,门田任。 图 系统 结构示 意 妞 系统结构 与工作 方式 系统结构如图 , 采用 种工作 方式 计算机全 自动控制方式 计算机依 据输人数据 , 经模型 计算输 出水量设定值 , 由流 量调 节仪表进行水量控制 。 计算机利 用 自学习 , 保证 系统 工作在最佳状态 〔 ‘ ’
(2)计算机操作指导计算机显示模型计 算水量设定值,由人工通过调节仪表系统进行 水量给定。 erminal Printer 1.3硬件系统的构成 由主机4MAC-5000,D/A输出扩展板 Ts1 -Qs1 uMAC-4030,开关量I/O板uMAC-4020及 Ts2” :4AC-5000 AC-4030 Qs2 终端和打印机构成,见图3〔23) -Qs3 Q1- 计算机系统的I/O信号有: Qs Q2- 模拟量输入: Q 一HD1 Q- #MAC-4020 -HMD2 Ts 精轧入口实测轧件温度; SW Ts2 精轧出口实测轧件温度影 Tss 终冷实测温度, 图3硬件系统 Twi 水冷器入水口实测温度: Fig.3 Hardware system Q1-Q4 实测4个分水管流量。 开关量输入: HMD 精轧人口轧件位置检测信号: HMD2 精轧出口轧件位置检测信号; SW 轧件进入冷床前端位置检测信号。 模拟量输出: Q;1-Qg4计算机对流量PID调节仪表的设定值。 1,4软件系统 由主程序、HMD2中断处理程序和S形中断处理程序构成,见图4。 HMD2 interrupt Start Self-learn for the predicting model of Ta Return 2's predicted calculating: Output Qat-Qs2 SW interrupt HMD OFF?> Self-learn for the setting model of quantity of water; Record the data of steel 图4程序泷程图 - Fig.4 The flowchart of program Return 72
计算机操作指导 计算机显示模型计 算水量设定值 , 由人工通过调节仪表系统进行 水量给定 。 。 硬件系统的构成 由 主 机 产 一 。 , 输 出扩展板 产 一 , 开 关 量 板 娜 一 及 终端和打印机构成 , 见 图 ‘ “ ’ ’ 。 计算机系统的 信号有 模拟量输人 、 精轧人 口实测轧件温度 , 精轧 出 口 实测轧件温度 , 终冷实测温度 , 水冷器人水 口实侧温度 一 ‘ 实测 个分水管流量 。 开关量输人 、 引 一 一 丁幻一 司。 叫 ,二 八 一 口 曰 ” ” 。一 图 硬件系统 卜】 时 叽 下 环 精轧人 口轧件位置检测信号 精轧 出 口轧件位置检测信号 轧件进入冷床前端位置检测信号 。 模拟量输 出 一 计算机对流量 调节仪表的设定值 。 软件系统 由主程序 、 中断处理程序和 牙 中断处理程序构成 , 见 图 。 ’ 、, 口。 厂 。 图 程序流程图 厂 一 日 门 七 主 印 、,吕 卜 助 ﹃ 兮声厅
系统启动后进入主程序。主程序对中断、通道及程序初始化和水量的初始设定。主程序 监视查询HMD:状态。当HMD,ON时,测量Ts,并进行Ts2预报,以此进行模型计算,输出 Qs1-Qs4设定值。 HMD2ON时请求中断。HMD2中断处理程序检测Ts2并进行终轧温度预报模型自学习。 SWON时请求中断。SW中断处理程序检测Ts3并进行水量设定模型自学习。 考虑到在轧制和冷却线上存在多轧件情况。只有前次终轧温度预报自学习完成后,进入精 轧机组的软件才参加终轧温度预报模型自学习。同样,只有前次水量自学习完成后进入水冷 区的轧件才参加水量设定模型自学习。为此,设立自学习指针对自学习跟踪。 对同一轧件在HMD:和HMD2间超过最大允许时间时,判断为事故状态。此时,该轧件 需从跟踪队列中取消,进行自学习跟踪调整。 温度、流量等信号超出合理值范围时,判断为系统异常,需进行自学习跟踪调整。 系统中使用WATCHDOG监视计算机系统。异常时发出声光报警。 人机联系通过终端显示器和键盘实现。显示器显示正常冷却参数和异常情况。系统初始 化中含有人工键入品种、规格等参数功能。在正常控制过程中,允许人工干预选择打印数 据,停止或继续系统运行。 2 数学模型 2.1终轧温度预报摸型 利用精轧机入口实测轧件温度Ts1,对该轧件的终轧温度(即开冷温度)进行预报。轧 件在通过精轧机组时,温度约1000℃左右,主要是辐射散热<4?。根据斯蒂芬一波尔茨曼定 律,热量损失dQ为: T 4 dQ=F.8. 100 (1) 式中: F一轧件的热辐射面积,m2; ε一热辐射系数, 0一斯蒂芬一波尔茨曼常数: T一一轧件的绝对温度,K; T一热辐射时间,S。 由于该热量损失造成的轧件温度降dt与dQ的关系式为: dQ=-G.cdi (2) 式中G一轧件质量,kg; t一一轧件温度,℃; c一钢材的比热kJ/kg·℃。 由(1),(2)式可得热平衡方程: -G.c.dt=F.e.o (3) 73
系统 启动后进入 主程序 。 主程序对 中断 、 通道及程序初始化和 水量的初始设定 。 主程序 监视查询 状态 。 当 时 , 测量 。 ,并进行 预报 , 以此进行模型计算 ,输 出 ,一 ‘ 设定值 。 时 请求中断 。 中断处理程序检测 并进行终轧温度 预报模型 自学 习 。 平 时 请求中断 。 班 中断处理程序检测 。 并进行 水量设定模型 自学习 。 考虑 到 在轧制和 冷却线上 存在 多轧 件情况 。 只有前次终轧温度预报 自学 习完成后 ,进人精 轧机组的 软件才参加终轧温度预报 模型 自学 习 。 同样 , 只有前次水量 自学 习完成后进人水冷 区 的轧 件才参加 水量设定模型 自学 习 。 为 此 , 设立 自学 习指针对 自学习跟踪 。 对 同一轧 件在 ,和 间超过最 大允许时 间时 , 判断为事故状态 。 此时 , 该轧件 需从跟踪队 列中取 消 , 迸行 自学习跟踪调整 。 温度 、 流 量等信号超 出合理值范 围时 , 判断 为系统异常 , 需进行 自学习跟踪调整 。 系 统 中使用 监视 计算机系统 。 异常时发 出声光报警 。 人机 联 系通过终端显示器 和键盘实现 。 显 示 器显示正常冷却参数和 异常情况 。 系统初始 化中含有人 工键人品 种 、 规格等参数功能 。 在正常控制过程 中 , 允许人 工干预 选 择 打 印数 据 , 停止或继续系统运行 。 数学模型 。 终 轧温 度预 报模型 利 用精轧机入 口 实测轧件温度 、 , 对该轧件的终轧温度 即开冷温度 进行预报 。 轧 件在通过精轧机组时 , 温度约 ℃ 左右 , 主要是辐射散热 〔 难 ’ 。 根据斯蒂芬一波 尔 茨曼定 律 , 热 量损失 为 二 。 £ 。 式 中 - 轧件 的热辐射 面积 , “ “ -热辐射系数 ‘ - 斯蒂芬一波 尔茨曼 常数 , - 轧 件的 绝对温度 , 丁 - 热辐射时 间 , 。 由于该热量损失造 成的轧件温度降 与 的 关系式为 一 一 式中 - 轧件质量 , 。 - 钢 材的 比热 吨 · ℃ 。 由 , 式可得热平衡方程 才一 一 轧件温度 , ℃ 一 ‘ 一 、 , 一 行 了 尸 『 侣 - , 丁
对(3)式整理、积分后得到: 个5:=100 0+(Tt23)÷-279 100G.c 100 (4) 式中: 个,2一预报开冷温度,℃,Ts1一实测精轧人口温度,℃; 下一轧件由精轧机入口至水冷器入口之间的行走时间,s。 当品种规格及轧制速度确定时,(4)式中的、0、c、G、F及T也为定值。对(4)式泰 勒展开,并舍去高次项得到: T52=a+b.Ts1+△ (5) 式中·、b为常数,△为误差项。使用现场实测数据进行回归,确定有关常数后得到实际使用 的终轧(开冷)温度预报模型: fs2=A+B.Ts1 (6) 2,2终轧温度预报模型自学习 该自学习是对(6)式中A进行修正。在各参数后加()表示第条钢的参数,则有: A(n)=Ts2(n)-B.Ts(n) (7) A(n+1)=A(n)+at4(n)-A(n)] (8) 式中a为加权系数,实验确定为0.3~0.5。 2.3水量设定模型 轧后控冷以自回火温度为目标值,根据开冷温度计算冷却水量(5)。由图1可知,在自 回火过程中,终冷温度与自回火温度基本为线性关系: Tss=a+b.T. (9) 式中:T一自回火温度;a、b一常数,由现场实测数据回归确定。 轧件在水冷器中冷却的传热过程为: ,dT=aAs(T-Tm) -Vpcp dT (10) 式中:p—一密度; Cp一一比热影 T一轧件温度; Tm一水温度: V一轧件体积影 A,一轧件在水冷器中的表面积; 2一水的传热系数。 对(10)式整理并积分得到: Ts3-Tw =e perr Ts2-Tw (11) 74
对 式整理 、 积分后得到 九 ” 〔瓷歌于 十 一 孺丝 一 〕 一 了 一 式中 -预报开冷温度 , ℃ - 实侧精轧 人 口温度 , ℃ 了 - 轧 件 由精轧机 人 口 至水冷器人 口 之 间的行走时间 , 。 当品种 规格及轧制速度确定时 , 式 中的 ‘ 、 、 、 、 及, 也为 定值 。 对 式 泰 勒展开 , 并舍去高次项得到 · 乙 式 中 、 为 常数 , 」为误差项 。 使 用现场实测数据进行 回归 , 确定有关常数后得到 实际使用 的 终轧 开冷 温度预报模型 二 一 。 终轧沮 度预 报模型 自学 习 该 自学习是对 式 中 进行修正 。 在 各参数后加 表 示第 条钢 的参数 , 则 有 称 , 一 王 〔 ” 一 〕 式 中 为加权系数 , 实验确定为 。 一 。 。 。 水级设 定模型 轧 后控冷以 自回火温度为 目标值 , 根据开冷温度计算冷却 水量 〔 “ ’ 。 由图 可 知 , 在 自 回火过程 中 , 终冷温度与 自回火温度基本为线性关系 一 式中 , - 自回火温度 、 - 常数 , 轧件在 水冷器 中冷却的传热过程 为 。 月 , 。 。 、 一 犷 。 一其土 几 、 一 二 。 一 尸 下 式中 - 密度 - 轧件温度 犷- 轧件体积 几-水的 传热系数 。 对 式整理并积分得到 由现场实侧数据回归确定 。 , - 比热 二 -水温度 , - 车 件在 水冷器中的表面积 一 二 一 二 , 尹 犷
式中: 下,一轧件在水冷器中滞留时间。对(11)式求对数得到: 2=-bV.cnCla(Ts:-Tv)-In(Ts2-T)] (12) As.T1 入与水流量及冷却器特性等有关。 轧件在水冷器中的冷却可认为是绕过潜体强制对流的换热过程r)。可采用McAdams:关 系式: NuPr-03=0.26Re08 (13) 式中:N。一努塞耳谁数,N:=,D为水冷器中藏冷却的长圆柱体的直径。 Pr一普兰特准数,Pr=C”,y为流体的粘度。 R Re一雷诺谁数,Re=V·D·巳,V为流体的流速。 由(12)和(13)式可得: w=-a026R…()ca(T-T)- P.cpV.D ln(T2-Tm)门o.6.” D.P (14) 当水冷器截面积为S时,流量Q为: Q=VS (15) 由(9)、(14)和(15)式得到: Q=ci-diCln(a+6Ts2-Tw)-In(Ts2-Tw)]0.6 (16) 当钢种和规格一定时,c,和d:为常数。 对(16)式泰勒展开并使用现场实测系数后得到水量设定计算模型为: Qs =c+d.Ts2+e-Tss+f.Tw (17) 式中Ts3为目标终冷温度值。 2,4水量设定自学习 采用指数平滑法,并确定(17)式中c为自学习项,则有: c(n)=Q(n)-[d.(Ts2(n)+e.Tss(n)+f.Tw(n)] (18) c(n+1)=c(n)+〔c(n)-c'(n)门 (19) 其中a为加权系数,取0.3~0.5。 75
式 中 了 , -轧件在 水冷器中滞 留时伺 。 对 式求对数得到 、 一 全竺竺〔 。 一 二 一 一 , 〕 月 下 凡与 水流量及 冷却器特性等有 关 。 轧件在水冷器中的冷却可 认为是绕过潜体强制对流的换热过程 〔 。 ’ 。 系式 一 ’ 。 可采 用 关 , 、 、 , 、 , 凡 式中 , 了 “ - 劳 签耳雌戮 , “ 、 了 , 为水冷 器中被冷却的长圆 柱体的直径 。 尸 -普兰特准数 , 尸 尸 峥 , ,为流体的粘度 。 -雷诺准数 , 。 二 由 和 式可得 犷尸 刀 。 竺 , 犷二 为流体的流速 。 罕 一 , 。 犷 。 , 一 丁 一 。 · 旱 一 。 ” 〔‘ · ‘ 一 甲 , 一︺ 一 一 二 〕 当 水冷器截 面积为 时 , 流 量 为 甲 由 、 和 式得到 · 。 一 二 一 、 一 二 〕 万 当钢种和规格一定时 , 、 和 为常数 。 对 式泰勒展开井使用现场实测系数后得到水量设定计算模型 为 , · 。 · 二 式 中 。 为 目标终冷温度值 。 。 水 设定 自学习 采用指数平滑法 , 并确定 式 中 。 为 自学习项 , 则有 二 一 〔小 一 , 一 , 〕 作 一 〔 一 , 〕 其中“ 为加权系数 , 取。 。
3控制结果 以下给出了批量生产中中14C20控冷结果。冷床实测自回火温度满足不大于士15℃的范 围要求。 炉批号 批号 冷床实测温度(℃) 水量(t/h) 01189-3 1981 550~560 160+150 D1-3067 单 560~580 130+130 D1-3707 单 570~580 120+120 D1-2436 单 560~580 125+125 0792-3 1987 550~580 105+105 D1-3035 1909 560~580 115+120 4结 束语 唐钢小型厂轧后控冷计算机控制系统投产以来工作稳定可靠,操作简单灵活。控制精度 满足终冷温祖度±15℃的工艺要求。应用该系统提高了产品等级和表面质量,减少二次氧化, 提高了成材率,减少了冷床弯废和设备事故,改善了劳动环境。综合经济效益为年500万 元。 参考文献 1刘植桢等。计算机控制,北京:清华大学出版社,1981 2 Analog Devices Inc.,uMAC-5000 Operation,1983 3 Analog Devices Inc.,uMAC-4030 User's Manual,1982 4俞昌铭。热传导及其数值分析,北京:清华大学出版社,1981 5石皎。相似理论在棒材轧后控冷中应用,北京科技大学硕士论文,1988 6G·H·盖格,D·R·波伊里尔;俞景禄,魏季和译,治金中的传热传质现象,北 京:治金工业出版社,1981 76
控 制 结 果 以下给出 了批量生 产中价 。 控冷结果 。 冷床实测 自回火温度满足不大于 士 ℃ 的范 围要求 。 炉批号 一 一 一 一 一 一 批号 冷床实测 温度 ℃ , 水量 单 结 束 语 唐钢小型厂轧后控冷计算机控制 系统投产以来工作稳定可靠 , 操作 简单灵活 。 控制精度 满足终冷温度 土 ℃ 的工艺要求 。 应 用该系统 提高了产品等级和表 面质量 , 减少二次氧化 , 提高了成材率 , 减少 了冷床弯废和设备事故 , 改善了劳动环境 。 综合经济 效 益 为 年 万 元 。 勺 参 考 文 献 刘植祯等 。 计算机控制 , 北京 清华大学 出版社 , , 产 一 , 。 , 产 一 ’ , 俞 昌 铭 热传导及其数值分析 , 北京 清华大学出 版社 , 石 皎 。 相似理论在棒材轧后控冷中应 用 , 北京科技 大学硕士论文 , · · 盖 格 , · · 波伊里尔 俞景禄 , 魏季和译 , · 冶 金中的传热传质 现 象 , 北 京 冶 金工业 出版社 ,