CN107999547A - 一种层流冷却的自学习方法及装置 - Google Patents
一种层流冷却的自学习方法及装置 Download PDFInfo
- Publication number
- CN107999547A CN107999547A CN201810040503.XA CN201810040503A CN107999547A CN 107999547 A CN107999547 A CN 107999547A CN 201810040503 A CN201810040503 A CN 201810040503A CN 107999547 A CN107999547 A CN 107999547A
- Authority
- CN
- China
- Prior art keywords
- mrow
- pyrometer
- cooling
- msup
- laminar
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001816 cooling Methods 0.000 title claims abstract description 184
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000017525 heat dissipation Effects 0.000 claims abstract description 49
- 230000005855 radiation Effects 0.000 claims abstract description 41
- 229910000831 Steel Inorganic materials 0.000 claims description 52
- 239000010959 steel Substances 0.000 claims description 52
- 230000004907 flux Effects 0.000 claims description 29
- 238000004364 calculation method Methods 0.000 claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 13
- 238000012546 transfer Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000004801 process automation Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 14
- 238000005098 hot rolling Methods 0.000 abstract description 3
- 238000002474 experimental method Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000005096 rolling process Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000007921 spray Substances 0.000 description 2
- 238000005507 spraying Methods 0.000 description 2
- 239000002436 steel type Substances 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 239000000498 cooling water Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B37/00—Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
- B21B37/74—Temperature control, e.g. by cooling or heating the rolls or the product
- B21B37/76—Cooling control on the run-out table
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B38/00—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product
- B21B38/006—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product for measuring temperature
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B2261/00—Product parameters
- B21B2261/20—Temperature
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Control Of Heat Treatment Processes (AREA)
- Heat Treatment Of Strip Materials And Filament Materials (AREA)
Abstract
本发明涉及一种层流冷却的自学习方法及装置,属于热轧带钢生产领域。该装置包括层流冷却区高温计I、层流冷却区高温计II、基础自动化系统和层流冷却控制系统;层流冷却区高温计I位于层流冷却区下游出口处,层流冷却区高温计II位于层流冷却区下游出口处至卷取机之间;在两个高温计之间构建一个空冷区域,基础自动化系统将采集到的层流冷却区高温计I和层流冷却区高温计II的实测值上传到层流冷却控制系统,通过数值计算求解等效的辐射散热系数,并用于层流冷却区的空冷散热量的计算。本发明仅需要在传统热轧生产线上增加一处测温点,对带钢生产过程中的空冷散热能力进行实时准确求解,用于提升层流冷却模型在后续计算时的预报精度。
Description
技术领域
本发明属于热轧带钢生产领域,涉及一种层流冷却的自学习方法及装置。
背景技术
在热轧带钢生产过程中,层流冷却控制系统的控制精度直接影响着成品带钢的产品性能。在精轧机出口到卷取机之间的输出辊道上,通常设置若干个冷却段,每个冷却段由设置在辊道上下两侧的若干喷水集管构成。从精轧机轧出的带钢在经过输出辊道时,通过开启一定数量的集管,使带钢长度方向上的各点在运行至卷取入口高温计时,达到目标的卷取温度。集管的开启数量由层流冷却模型计算确定。
为了掌握带钢在层流冷却区入口和出口的实际温度,通常在层流冷却区入口和出口设置有高温计,一些生产线还会在层流冷却区中间位置另外设置高温计,用于检测带钢在冷却区中间位置的实际温度。冷却区入口的实测温度还作为层流冷却模型的输入数据,用于层流冷却模型的设定计算。
带钢在输出辊道上的冷却过程可分为两类:水冷温降和空冷温降。同时还包含有带钢和层冷辊道之间的传导散热,但由于难以准确计算,通常将其等效至水冷温降和空冷温降之中。水冷温降主要由带钢与冷却水之间的换热引起,其换热能力的大小由传热系数表示。空冷温降主要由辐射散热和对流散热引起,由辐射散热系数表示。
层流冷却模型具备自学习功能,通过收集带钢上的样本点在两个高温计之间的冷却历程和实测温度,通过一定的自学习逻辑,对设定计算时使用的水冷换热参数进行修正,用于提高后续计算的预报精度。基于已知条件,在两个高温计之间仅能构建一个非线性方程,难以同时求解水冷换热系数和辐射散热系数,因此在解方程时,辐射散热系数设置为固定值,仅对水冷换热系数进行求解。
在模型调试中,辐射散热系数通常采用空冷实验的方法进行分析确定,随后各钢种在生产中即采用此实验分析值。通过对生产数据的分析来看,采用实验分析值并不够精确,带来较大的误差。
采用空冷实验分析辐射散热系数的方法,将影响实验带钢的产品性能,成本高;对生产的各个钢种都进行空冷实验分析辐射散热系数,成本高,不具有可操作性;空冷实验分析出的辐射散热系数具备时效性,随着外部环境的改变,辐射散热系数将发生改变,从而影响模型预测精度。
发明内容
有鉴于此,本发明的目的在于提供一种层流冷却的自学习方法及装置,仅需要在传统热轧生产线上增加一处测温点,对带钢生产过程中的空冷散热能力进行实时准确求解,用于提升层流冷却模型在后续计算时的预报精度。
为达到上述目的,本发明提供如下技术方案:
一种层流冷却的自学习装置,该装置包括层流冷却区高温计I、层流冷却区高温计II、基础自动化系统和层流冷却控制系统;
所述层流冷却区高温计I位于层流冷却区下游出口处,层流冷却区高温计II位于层流冷却区下游出口处至卷取机之间;在两个高温计之间构建一个空冷区域,基础自动化系统将采集到的层流冷却区高温计I和层流冷却区高温计II的实测值上传到层流冷却控制系统,通过数值计算求解等效的辐射散热系数,并用于层流冷却区的空冷散热量的计算。
进一步,所述层流冷却区高温计I与卷取机之间还包括若干个高温计。
一种层流冷却的自学习方法,该方法包括以下步骤:
S1:在层流冷却区下游出口处的层流冷却区高温计I,与其相距若干米的位置上,安装层流冷却区高温计II,并将两个高温计测量数据接入基础自动化系统;
S2:过程自动化系统的网关进程负责接收新增高温计的实测数据,并传给层流冷却控制系统进行处理,将处理后的数据将转发给层流冷却模型;
S3:层流冷却模型在收到数据后,进行自学习处理,学习结果为等效辐射散热系数;
S4:层流冷却模型在执行设定计算时,将经过学习修正的辐射散热系数,作为输入参数用于设定计算,提升准确度。
进一步,在步骤S1中,所述安装层流冷却区高温计II的数量为若干,适当调整各个高温计的间距,并将层流冷却区下游出口处的与卷取机之间所有的高温计的测量数据接入基础自动化系统。
进一步,所述步骤S3具体为:
S31:通过一维非稳态热传导方程来拟合带钢内部的热传导,
其中,T=T(t,l,y)表示在t时刻带钢在长度方向l处厚度方向y处的温度分布,ρ表示带钢密度,c表示带钢的比热容,λ(T)表示带钢的热传导系数;
S32:根据初始条件,T(t0,0,y)=T0,分析边界条件,
其中,表示在冷却区x处,通过上表面的热流密度,表示在冷却区x处,通过下表面的热流密度,t0表示初始时刻,T0表示初始温度,
即,
其中,x表示带钢所处的冷却区位置,d表示带钢厚度的一半,L表示带钢的总长度,φtop,l(x,Tl(x,d))表示水冷换热时带钢上表面的热流密度,φbot,l(x,Tl(x,-d))表示水冷换热时带钢下表面的热流密度,xs表示热传导方程计算开始位置,xE表示热传导方程计算结束位置;
S33:将带钢表面的热流密度分为空气热流密度和水对流热流密度两部分,
其中,φi,l(x,T)表示带钢在冷却区x处,通过上下表面的热流密度,表示带钢在冷却区x处,通过上下表面的空气热流密度,表示带钢在冷却区x处,通过上下表面的水对流热流密度,表示带钢上下表面辐射散热的热流密度
其中,ε表示等效辐射系数,σ表示玻尔兹曼常数,Tamb表示环境温度;
S34:对一维非稳态热传导方程进行化简为,
其中,v为带钢运行速度;
S35:根据化简后的一维非稳态热传导方程以及步骤S33,并根据层流冷却区高温计I和层流冷却区高温计II的实测温度,采用数值方法求解出带钢上下表面辐射散热的热流密度,进而计算出等效热辐射散热系数。
本发明的有益效果在于:本发明仅需要在传统热轧生产线上增加一处测温点,对带钢生产过程中的空冷散热能力进行实时准确求解,用于提升层流冷却模型在后续计算时的预报精度。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为传统层流冷却系统高温计布置图;
图2为本发明的层流冷却系统高温计布置图。
图中标记:1-精轧机末机架;2-层流冷却下部冷却段;3-层流冷却上部冷却段;4-冷却中的带钢;5-卷取机;6-精轧出口高温计;7-中间高温计;8-层流冷却区下游出口处的层流冷却区高温计I;9-层流冷却区高温计I至卷取机之间的层流冷却区高温计II。
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述。
1)传统层流冷却控制系统的高温计布置,如图1所示:1为精轧机末机架;2为层流冷却下部冷却段;3为层流冷却上部冷却段;4为冷却中的带钢;5为卷取机;6为精轧出口高温计;7为中间高温计;8为层流冷却区下游出口处的层流冷却区高温计I。
在任意两个高温计之间,基于入口和出口高温计的测量数据,以及带钢样本点在两个高温计之间的冷却历程,可建立一个自学习区间。每一个自学习区间的学习结果,通常仅用于本区间的设定计算。
每一个自学习区间内同时存在水冷换热和空冷散热,在执行自学习计算时,基于已知的数据,任一个自学习区间仅能构建一个非线性方程,水冷传热系数和辐射散热系数之间,仅能求解二者中的一个。水冷散热的效率远高于空冷换热,水冷散热为主要影响因素,因此在自学习计算时,辐射散热系数固定采用空冷实验的分析值,仅对水冷传热系数进行求解。
通过空冷实验分析的辐射散热系数,存在局限性:①采用空冷实验分析辐射散热系数,需要在生产中对某一卷带钢全长或部分长度上,不进行喷水冷却,用于分析带钢在此过程中的等效辐射散热系数,将会对实验带钢的产品性能带来不利影响。②带钢在空冷区散热时所适用的辐射散热系数,实际上是包含了对流传热、辊道接触传热的等效辐射散热系数,不同钢种之间的辐射散热系数可能存在较大差异。如果对各类钢种都进行空冷实验分析其辐射散热系数,成本较高,不具备可操作性。③通常在投产初期进行空冷实验,对辐射散热系数进行分析。随着外部条件的改变,等效的辐射散热系数将会发生改变,投产初期经空冷实验分析而来的辐射散热系数可能不再适用。
2)本发明中的层流冷却控制系统的高温计布置,如图2所示:9为层流冷却区高温计I至卷取机之间的层流冷却区高温计II。
在卷取层流冷却区高温计I的下游,增设一个新的层流冷却区高温计II,在层流冷却区高温计I和层流冷却区高温计II之间构建了一个纯空冷区域。高温计层流冷却区高温计I、层流冷却区高温计II的实测温度经由基础自动化系统采集后,上传给层流冷却控制系统;带钢样本点从层流冷却区高温计I运行至层流冷却区高温计II的时间,可根据带钢的实际速度和层流冷却区高温计I-层流冷却区高温计II之间的距离计算而来;
带钢在层流冷却区高温计I和层流冷却区高温计II之间的温度场变化,可由传热微分方程进行求解,边界条件,即等效辐射散热系数为未知数,在层流冷却区高温计I-层流冷却区高温计II之间可构建一个非线性方程,通过数值计算求解出等效辐射散热系数,步骤具体为:
由于,带钢内部热传导可以用一维非稳态热传导方程来描述,其中T=T(t,l,y),是关于时刻t、带钢与头部距离的位置点l和带钢厚度方向的位置y的函数:
其中,T=T(t,l,y)表示在t时刻带钢在长度方向l处厚度方向y处的温度分布,t=0时刻,带钢头部通过x=0位置,ρ表示带钢密度,c表示带钢的比热容,λ(T)表示带钢的热传导系数;
初始条件:
T(t0,0,y)=T0(T0为两个高温计之间,起始点的带钢表面温度,可由高温计实测而来);
边界条件(傅立叶定律):
top:
bot:
其中,表示在冷却区x处,通过上表面的热流密度,表示在冷却区x处,通过下表面的热流密度,t0表示初始时刻,T0表示初始温度,
即:
其中,x表示带钢所处的冷却区位置,d表示带钢厚度的一半,L表示带钢总长度,φtop,l(x,Tl(x,d))表示水冷换热时带钢上表面的热流密度,φbot,l(x,Tl(x,-d))表示水冷换热时带钢下表面的热流密度,xs表示热传导方程计算开始位置,xE表示热传导方程计算结束位置;
带钢表面的热流密度φi,l(T)(i∈{top,bot})由两部分构成:
热辐射产生的热流密度与空气/水对流的热流密度即
其中,φi,l(x,T)表示带钢在冷却区x处,通过上下表面的热流密度,表示带钢在冷却区x处,通过上下表面的空气热流密度,表示带钢在冷却区x处,通过上下表面的水对流热流密度,表示带钢上下表面辐射散热的热流密度。
其中辐射散热的热流密度为:
其中,ε表示等效辐射系数,σ表示玻尔兹曼常数,Tamb表示环境温度。
由得到v为带钢的运行速度;
且则热传导方程(1)变形为:
该方程描述在冷却区位置x处,带钢长度方向位置l处,厚度方向上y处的温度分布T。
带钢表面上某一样本点在层流冷却区高温计I 8和层流冷却区高温计II 9处的表面温度为T0和Tl,均由高温计测量而来,样本点在层流冷却区高温计I 8和层流冷却区高温计II 9之间的运动时的初始速度v0和加速度a来自基础自动化系统或精轧设定模型。
联立方程(3)和(4),已知T0、Tl、v0、a,采用数值计算方法求解出未知数ε,即等效辐射散热系数。
层流冷却模型在执行设定计算时,经过学习修正的辐射散热系数,将作为输入参数用于设定计算,提升设定计算的准确度;
层流冷却区的每一个冷却段,均由设置在辊道上下两侧的若干喷水集管构成,由电磁阀门控制每个集管的开闭状态。
设定计算的目的,是根据给定的开阀方式,计算出所需的阀门开启个数,以使得带钢上的每一点在运动到卷取入口高温计8时,达到所需的目标温度。
已知精轧出口高温计6处的带钢表面温度(由高温计实测,或者由其他模型预测而来)、带钢点运动到高温计8处的表面温度目标、带钢点运动到高温计6处的起始速度、带钢点在6~8之间运行时的每个加速段的加速度和加速时间,已知带钢在开阀区进行喷水冷却时的热流密度,在关阀区域为空冷散热,其热流密度由自学习确定的等效辐射散热系数ε求得,开阀的个数p为未知数。采用数值方法求解可得到所需要的开阀个数p。
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。
Claims (5)
1.一种层流冷却的自学习装置,其特征在于:该装置包括层流冷却区高温计I、层流冷却区高温计II、基础自动化系统和层流冷却控制系统;
所述层流冷却区高温计I位于层流冷却区下游出口处,层流冷却区高温计II位于层流冷却区下游出口处至卷取机之间;在两个高温计之间构建一个空冷区域,基础自动化系统将采集到的层流冷却区高温计I和层流冷却区高温计II的实测值上传到层流冷却控制系统,通过数值计算求解等效的辐射散热系数,并用于层流冷却区的空冷散热量的计算。
2.根据权利要求1所述的一种层流冷却的自学习装置,其特征在于:所述层流冷却区高温计I与卷取机之间还包括若干个高温计。
3.一种层流冷却的自学习方法,其特征在于:该方法包括以下步骤:
S1:在层流冷却区下游出口处的层流冷却区高温计I,与其相距若干米的位置上,安装层流冷却区高温计II,并将两个高温计测量数据接入基础自动化系统;
S2:过程自动化系统的网关进程负责接收新增高温计的实测数据,并传给层流冷却控制系统进行处理,将处理后的数据将转发给层流冷却模型;
S3:层流冷却模型在收到数据后,进行自学习处理,学习结果为等效辐射散热系数;
S4:层流冷却模型在执行设定计算时,将经过学习修正的辐射散热系数,作为输入参数用于设定计算,提升准确度。
4.根据权利要求3所述的一种层流冷却的自学习方法,其特征在于:在步骤S1中,所述安装层流冷却区高温计II的数量为若干,适当调整各个高温计的间距,并将层流冷却区下游出口处的与卷取机之间所有的高温计的测量数据接入基础自动化系统。
5.根据权利要求3所述的一种层流冷却的自学习方法,其特征在于:所述步骤S3具体为:
S31:通过一维非稳态热传导方程来拟合带钢内部的热传导,
<mrow>
<mi>&rho;</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>&lambda;</mi>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中,T=T(t,l,y)表示在t时刻带钢在长度方向l处厚度方向y处的温度分布,ρ表示带钢密度,c表示带钢的比热容,λ(T)表示带钢的热传导系数;
S32:根据初始条件,T(t0,0,y)=T0,分析边界条件,
其中,表示在冷却区x处,通过上表面的热流密度,表示在冷却区x处,通过下表面的热流密度,t0表示初始时刻,T0表示初始温度,
即,
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>T</mi>
<mi>l</mi>
</msup>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>d</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>T</mi>
<mi>l</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>&phi;</mi>
<mrow>
<mi>t</mi>
<mi>o</mi>
<mi>p</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<msup>
<mi>T</mi>
<mi>l</mi>
</msup>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>d</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>l</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mi>L</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>x</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mi>S</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>E</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>T</mi>
<mi>l</mi>
</msup>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mo>-</mo>
<mi>d</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>T</mi>
<mi>l</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mo>-</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>&phi;</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>t</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<msup>
<mi>T</mi>
<mi>l</mi>
</msup>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mo>-</mo>
<mi>d</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>l</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mi>L</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>x</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mi>S</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>E</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,x表示带钢所处的冷却区位置,d表示带钢厚度的一半,L表示带钢的总长度,φtop,l(x,Tl(x,d))表示水冷换热时带钢上表面的热流密度,φbot,l(x,Tl(x,-d))表示水冷换热时带钢下表面的热流密度,xs表示热传导方程计算开始位置,xE表示热传导方程计算结束位置;
S33:将带钢表面的热流密度分为空气热流密度和水对流热流密度两部分,
<mrow>
<msup>
<mi>&phi;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mtable>
<mtr>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
<mi>r</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>q</mi>
<mi>u</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mi>t</mi>
<mi>o</mi>
<mi>p</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>t</mi>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mrow>
其中,φi,l(x,T)表示带钢在冷却区x处,通过上下表面的热流密度,表示带钢在冷却区x处,通过上下表面的空气热流密度,表示带钢在冷却区x处,通过上下表面的水对流热流密度,表示带钢上下表面辐射散热的热流密度
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&epsiv;</mi>
<mi>&sigma;</mi>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>T</mi>
<mo>+</mo>
<mn>273</mn>
<mo>,</mo>
<mn>15</mn>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>273</mn>
<mo>,</mo>
<mn>15</mn>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mi>t</mi>
<mi>o</mi>
<mi>p</mi>
<mo>,</mo>
<mi>b</mi>
<mi>o</mi>
<mi>t</mi>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,ε表示等效辐射系数,σ表示玻尔兹曼常数,Tamb表示环境温度;
S34:对一维非稳态热传导方程进行化简为,
<mrow>
<mi>&rho;</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mi>v</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>&lambda;</mi>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中,v为带钢运行速度;
S35:根据化简后的一维非稳态热传导方程以及步骤S33,并根据层流冷却区高温计I和层流冷却区高温计II的实测温度,采用数值方法求解出带钢上下表面辐射散热的热流密度,进而计算出等效热辐射散热系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810040503.XA CN107999547B (zh) | 2018-01-16 | 2018-01-16 | 一种层流冷却的自学习方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810040503.XA CN107999547B (zh) | 2018-01-16 | 2018-01-16 | 一种层流冷却的自学习方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107999547A true CN107999547A (zh) | 2018-05-08 |
CN107999547B CN107999547B (zh) | 2023-10-13 |
Family
ID=62050760
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810040503.XA Active CN107999547B (zh) | 2018-01-16 | 2018-01-16 | 一种层流冷却的自学习方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107999547B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112122361A (zh) * | 2019-06-25 | 2020-12-25 | 上海梅山钢铁股份有限公司 | 一种防止中高碳钢开裂的层流冷却控制方法 |
CN113083913A (zh) * | 2021-03-31 | 2021-07-09 | 北京首钢股份有限公司 | 卷取带钢温度控制方法、装置及层冷系统 |
CN114682632A (zh) * | 2020-12-29 | 2022-07-01 | 唐山学院 | 汽车用热轧双相钢的控冷工艺设计方法 |
CN115931152A (zh) * | 2022-12-08 | 2023-04-07 | 鞍钢股份有限公司 | 一种基于层流控制的高温计信号优化方法、系统及存储介质 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5886922A (ja) * | 1981-11-18 | 1983-05-24 | Mitsubishi Heavy Ind Ltd | 鋼板の急速冷却装置 |
JP2004160532A (ja) * | 2002-11-15 | 2004-06-10 | Nippon Steel Corp | 熱延鋼板の冷却制御方法 |
CN101745551A (zh) * | 2008-12-11 | 2010-06-23 | 宝山钢铁股份有限公司 | 一种热轧带钢自由冷却方法 |
CN102284511A (zh) * | 2011-07-28 | 2011-12-21 | 山西太钢不锈钢股份有限公司 | 一种带钢层流冷却温度自适应方法 |
CN102392119A (zh) * | 2011-10-28 | 2012-03-28 | 重庆赛迪工业炉有限公司 | 热镀锌连续退火炉的在线综合控制方法 |
KR20140005461A (ko) * | 2012-07-04 | 2014-01-15 | 주식회사 포스코 | 열간압연장치 및 열간압연재의 제조방법 |
CN103559334A (zh) * | 2013-10-15 | 2014-02-05 | 华中科技大学 | 一种层流冷却温度场的建模方法及系统 |
CN103878186A (zh) * | 2014-03-24 | 2014-06-25 | 华中科技大学 | 一种确定热轧带钢层流冷却温度的方法 |
CN104226699A (zh) * | 2014-08-21 | 2014-12-24 | 南京钢铁股份有限公司 | 一种层流冷却正、负反馈的闭环控制方法 |
CN105327949A (zh) * | 2015-11-25 | 2016-02-17 | 北京金自天正智能控制股份有限公司 | 一种热轧带钢卷取温度的流量控制方法 |
CN106862283A (zh) * | 2017-02-10 | 2017-06-20 | 武汉钢铁工程技术集团自动化有限责任公司 | 热轧超快冷过程控制系统 |
-
2018
- 2018-01-16 CN CN201810040503.XA patent/CN107999547B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5886922A (ja) * | 1981-11-18 | 1983-05-24 | Mitsubishi Heavy Ind Ltd | 鋼板の急速冷却装置 |
JP2004160532A (ja) * | 2002-11-15 | 2004-06-10 | Nippon Steel Corp | 熱延鋼板の冷却制御方法 |
CN101745551A (zh) * | 2008-12-11 | 2010-06-23 | 宝山钢铁股份有限公司 | 一种热轧带钢自由冷却方法 |
CN102284511A (zh) * | 2011-07-28 | 2011-12-21 | 山西太钢不锈钢股份有限公司 | 一种带钢层流冷却温度自适应方法 |
CN102392119A (zh) * | 2011-10-28 | 2012-03-28 | 重庆赛迪工业炉有限公司 | 热镀锌连续退火炉的在线综合控制方法 |
KR20140005461A (ko) * | 2012-07-04 | 2014-01-15 | 주식회사 포스코 | 열간압연장치 및 열간압연재의 제조방법 |
CN103559334A (zh) * | 2013-10-15 | 2014-02-05 | 华中科技大学 | 一种层流冷却温度场的建模方法及系统 |
CN103878186A (zh) * | 2014-03-24 | 2014-06-25 | 华中科技大学 | 一种确定热轧带钢层流冷却温度的方法 |
CN104226699A (zh) * | 2014-08-21 | 2014-12-24 | 南京钢铁股份有限公司 | 一种层流冷却正、负反馈的闭环控制方法 |
CN105327949A (zh) * | 2015-11-25 | 2016-02-17 | 北京金自天正智能控制股份有限公司 | 一种热轧带钢卷取温度的流量控制方法 |
CN106862283A (zh) * | 2017-02-10 | 2017-06-20 | 武汉钢铁工程技术集团自动化有限责任公司 | 热轧超快冷过程控制系统 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112122361A (zh) * | 2019-06-25 | 2020-12-25 | 上海梅山钢铁股份有限公司 | 一种防止中高碳钢开裂的层流冷却控制方法 |
CN114682632A (zh) * | 2020-12-29 | 2022-07-01 | 唐山学院 | 汽车用热轧双相钢的控冷工艺设计方法 |
CN113083913A (zh) * | 2021-03-31 | 2021-07-09 | 北京首钢股份有限公司 | 卷取带钢温度控制方法、装置及层冷系统 |
CN115931152A (zh) * | 2022-12-08 | 2023-04-07 | 鞍钢股份有限公司 | 一种基于层流控制的高温计信号优化方法、系统及存储介质 |
CN115931152B (zh) * | 2022-12-08 | 2023-10-20 | 鞍钢股份有限公司 | 一种基于层流控制的高温计信号优化方法、系统及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107999547B (zh) | 2023-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107999547A (zh) | 一种层流冷却的自学习方法及装置 | |
CN107066737B (zh) | 一种预测热轧过程板带温度场的二维交替差分方法 | |
RU2404000C2 (ru) | Способ управления охлаждением, устройство управления охлаждением и устройство расчета количества охлаждающей воды | |
CN101664793B (zh) | 基于红外热成像的连铸坯实时温度场在线预测方法 | |
CN111079275B (zh) | 一种板带热轧生产线轧件温度获取方法及装置 | |
CN103878186B (zh) | 一种确定热轧带钢层流冷却温度的方法 | |
JPWO2008012881A1 (ja) | 圧延ラインの材質予測および材質制御装置 | |
Wang et al. | Inverse problem-coupled heat transfer model for steel continuous casting | |
CN109013717B (zh) | 一种热连轧中间坯心部温度计算方法 | |
Muhin et al. | Simulation of accelerated strip cooling on the hot rolling mill run-out roller table | |
CN103028615A (zh) | 一种预测带钢热连轧过程温度演变的方法 | |
CN112036056A (zh) | 一种热轧带钢层流冷却有限元建模方法 | |
CN115121626B (zh) | 一种基于误差补偿的热轧带钢瞬态热辊型预报方法 | |
CN103559334B (zh) | 一种层流冷却温度场的建模方法及系统 | |
CN105414205B (zh) | 一种基于plc的钢板温度在线预测方法 | |
Pian et al. | A hybrid soft sensor for measuring hot-rolled strip temperature in the laminar cooling process | |
CN103878185B (zh) | 一种用于热轧层流冷却的动态段冷却控制方法 | |
CN103990653B (zh) | 精轧入口温度命中精度确保方法 | |
CN208162301U (zh) | 一种层流冷却的自学习装置 | |
Zhi-Kun et al. | Application of weighted multiple models adaptive controller in the plate cooling process | |
CN105149361B (zh) | 一种热连轧机带材换热系数自学习方法 | |
CN105537284A (zh) | 一种热连轧精轧入口温度预报方法 | |
CN116144888B (zh) | 基于横纵温差的双相钢板带均质化挂卷及控冷调质方法 | |
CN101633004B (zh) | 厚板轧后控制冷却广义观测器设计方法 | |
CN103586290A (zh) | 一种层流冷却过程带钢温度软测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |