CN104164537B - 一种炼钢过程中真空精炼物理模拟试验方法及装置 - Google Patents
一种炼钢过程中真空精炼物理模拟试验方法及装置 Download PDFInfo
- Publication number
- CN104164537B CN104164537B CN201410426265.8A CN201410426265A CN104164537B CN 104164537 B CN104164537 B CN 104164537B CN 201410426265 A CN201410426265 A CN 201410426265A CN 104164537 B CN104164537 B CN 104164537B
- Authority
- CN
- China
- Prior art keywords
- section
- comer
- vertical cross
- laser
- fluid velocity
- 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.)
- Active
Links
Landscapes
- Laser Beam Processing (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明提供一种炼钢过程中真空精炼物理模拟试验方法及装置,方法包括:根据待测RH真空精炼装置建立模拟装置;运行模拟装置,并在模拟装置中添加示踪粒子;利用激光器照射下降管垂直截面;利用相机对激光器照射的垂直截面进行连续拍摄;根据测量图像中示踪粒子的移动距离和示踪粒子的拍摄时间生成垂直截面的流体速度分布;对不同垂直截面的流体速度进行拟合和积分,得到RH的循环流量。通过在模拟装置中添加示踪粒子,并利用激光照射垂直截面,利用相机进行连续拍照,获取示踪粒子的移动距离和时间生成垂直截面的流体速度分布,经过拟合和积分得到RH的循环流量,不仅提高了截面的流体速度的测量精度,而且降低了对测量模型内流场的影响。
Description
技术领域
本发明涉及炼钢领域中的真空精炼物理模拟试验领域,尤其涉及炼钢过程中真空精炼物理模拟测量方法及装置。
背景技术
目前,RH物理模拟循环流量测量的主要方法分为直接法和间接法,直接法主要有节流法和溢流法两种,间接法主要有皮托管、示踪法、流速仪测速法、超声波流量计等几种方法。节流法和溢流法可以很准确地测量RH的循环流量,但这两种方法的缺点是改变了钢包内钢液的流场,不适合同时考虑模型内流动情况又测量循环流量的试验。皮托管法是通过测量流体动压和静压的压力差来测量流速的,原理简单、成本较低,测量时对模型内流场影响很小,在RH物理模拟循环流量测量中使用较多。但是使用皮托管测量RH循环流量时,试验结果的稳定性不好、误差比较大。示踪法测量RH循环流量的误差较大些。流速仪测速法是将智能流速仪和流速旋浆传感器放到水中,因此也会对流场产生一定影响。超声波流量计是一种不改变模型流场的测量方法,但这种方法对测量小管径,特别是小流量时的误差较大。
粒子成像测速(Particle Image Velocimetry,简称PIV)是一种瞬态、多点、无接触式的激光流体力学测速方法。PIV超出了单点测速技术的局限性,能在同一瞬态记录下大量空间点上的速度分布信息,并可提供丰富的流场空间结构以及流动特性。PIV是通过测量示踪粒子在已知很短时间间隔内的位移来间接地测量流场的瞬态速度分布,所有测量装置并不介入流场,具有较高的测量精度。
PIV这一在上世纪七十年代发展起来的流动显示技术,如今已经成为人们对涉及流体力学的诸多学科的得力科研工具,如在水洞、风洞、发动机、管道、河流等各种水和空气中已进行了大量试验研究。在炼钢领域的水模拟试验中也有了较多的应用,但目前研究主要集中在结晶器、中间包的流场测量,利用PIV测量RH的循环流量并无先例。
发明内容
为了克服现有RH物理模拟循环流量测量方法测量精度差或对流场产生干扰的不足,本发明提供一种炼钢过程中真空精炼物理模拟测量方法及装置,具有不干扰被测流场、测量精度高的特点。
本发明提供一种炼钢过程中真空精炼物理模拟试验方法,所述方法包括:
根据待测RH真空精炼装置建立模拟装置;
运行模拟装置,并在模拟装置中添加示踪粒子;
在模拟装置的测量位置选取垂直截面,利用激光器照射所述垂直截面;
利用相机对激光器照射的垂直截面进行连续拍摄,生成测量图像;
识别测量图像中的示踪粒子,根据测量图像中示踪粒子的移动距离和示踪粒子的拍摄时间生成垂直截面的流体速度分布;
对不同垂直截面中心高度处的流体速度进行拟合和积分,得到RH的循环流量。
优选的,所述根据待测RH真空精炼装置建立模拟装置,包括:
采用有机玻璃按照模型与原型1∶N的比例制作RH物理模型,其中2≤N≤10,并在真空室下降管外部及钢包外部加装水盒子;
所述运行模拟装置,并在模拟装置中添加示踪粒子,包括:
将适量PIV专用颗粒在水中搅拌、润湿后加入钢包,调节供气装置与真空泵以控制吹气量与抽气量,使RH物理模型循环达到稳定;
所述在模拟装置的测量位置选取垂直截面,利用激光器照射所述垂直截面,包括:
在钢包中加入预定量的水,调节真空泵使真空室内压强达到原型对应值,水盒子内充满水,启动PIV,待激光器电源预热完毕,调节片光源镜头组,使片状激光照射在RH下降管某一垂直截面,将标定板或直尺放置在该垂直面内,即片状激光照射区域,开启摄像模式,调节CCD相机的焦距与光圈,使画面中标定板或直尺刻度显示达到最清晰,即使CCD相机对焦在片状激光所照射的垂直截面上;
所述利用相机对激光器照射的垂直截面进行连续拍摄,生成测量图像,包括;
由计算机内置的可编程时间控制器调控激光器电源与CCD相机,连续拍摄照片,生成测量图像;
所述识别测量图像中的示踪粒子,根据测量图像中示踪粒子的移动距离和示踪粒子的拍摄时间生成垂直截面的流体速度分布,包括:
利用PIV图像数据处理***,计算出所拍图像瞬时速度的平均值,即为该下降管垂直截面的流体速度分布;
所述对不同垂直截面中心高度处的流体速度进行拟合和积分,得到RH的循环流量,包括:
根据下降管不同垂直截面中心高度处的流体速度,得到下降管中心水平截面内若干点的流体速度,利用Matlab数学软件,采用多项式拟合流体速度曲面并得到拟合曲面函数,用该软件编写程序,计算下降管内流体速度的积分,进一步转换为RH的循环流量。
优选的,所述水盒子为有机玻璃制作的上端开口的长方体水槽,加装于CCD相机所拍摄流体区域的圆柱体外侧。
优选的,所述激光为532nm绿光,该激光峰值能量400mJ/Pulse。
优选的,所述PIV专用颗粒是直径50μm的空心玻璃球。
优选的,所述CCD相机是透过下降管外的水盒子进行拍摄。
优选的,所述下降管不同垂直截面是
从所选取的垂直截面按照一定间距向两侧选取下降管截面直到下降管边界。
本发明提供一种炼钢过程中真空精炼物理模拟试验模拟装置,所述装置包括模拟装置和粒子成像测速装置;所述模拟装置包括钢包,真空室,供气装置,真空泵,U型压力计,钢包水盒子,下降管,下降管水盒子;所述钢包顶部与所述真空室的底部连接,所述真空室顶部与所述真空泵链接,所述供气装置与所述真空室的侧壁连接,所述钢包水盒子设置在所述钢包的侧壁外侧,所述上升管和下降管穿过所述真空室底部,所述上升管管壁设置有通气孔,所述上升管通过所述通气孔与所述供气装置的输出端连接,所述下降管水盒子设置在所述下降管的管壁外侧;
所述粒子成像测速装置包括:激光器电源,双腔激光器,导光臂,计算机内置可编程时间控制器,PIV图像数据处理***,CCD相机,片光源镜头组,PIV专用颗粒;所述激光器电源与所述激光器连接,所述激光器的输出端与所述导光臂连接,所述导光臂的输出端与所述片光源镜头组连接,所述PIV图像数据处理***与所述计算机内置可编程时间控制器和所述CCD相机连接。
本发明的上述技术方案的有益效果如下:
上述方案中,通过在模拟装置中添加示踪粒子,并利用激光照射垂直截面,利用相机进行连续拍照,获取示踪粒子的移动距离和时间生成垂直截面的流体速度分布,经过拟合和积分得到RH的循环流量,不仅提高了截面的流体速度的测量精度,而且降低了对测量模型内流场的影响。
附图说明
图1为本发明的RH真空精炼循环流体模拟测量方法流程图;
图2为本发明的RH真空精炼循环流体模拟测量装置的模拟装置结构图;
图3为本发明的RH真空精炼循环流体模拟测量装置的粒子成像测速装置结构图;
图4为实施例下降管中心水平截面内流体速度取样点示意图;
图5为实施例1运用Matlab数学软件拟合效果图;
图6为实施例2运用Matlab数学软件拟合效果图;
图7为实施例3运用Matlab数学软件拟合效果图;
图8为实施例4运用Matlab数学软件拟合效果图。
[主要元件符号说明]
1、钢包;
2、真空室;
3、供气装置;
4、真空泵;
5、钢包水盒子;
6、上升管;
7、下降管;
8、下降管水盒子;
9、气体流量计;
10、分气盘;
11、U型压力计;
12、激光器电源;
13、激光器;
14、导光臂;
15、片光源镜头组;
16、时间控制器;
17、数据处理器;
18、CCD相机。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
如图1所示,本发明的实施例的RH真空精炼循环流体模拟测量方法流程图,所述方法包括:
步骤101:根据待测RH真空精炼装置建立模拟装置。
其中,采用有机玻璃按照模型与原型1∶N的比例制作RH物理模型,其中2≤N≤10,并在真空室下降管外部及钢包外部加装水盒子。
步骤102:运行模拟装置,并在模拟装置中添加示踪粒子。
其中,将适量PIV专用颗粒在水中搅拌、润湿后加入钢包,调节供气装置与真空泵以控制吹气量与抽气量,使RH物理模型循环达到稳定。
步骤103:在模拟装置的测量位置选取垂直截面,利用激光器照射所述垂直截面。
其中,在钢包中加入预定量的水,调节真空泵使真空室内压强达到原型对应值,水盒子内充满水,启动PIV,待激光器电源预热完毕,调节片光源镜头组,使片状激光照射在RH下降管某一垂直截面,将标定板或直尺放置在该垂直面内,即片状激光照射区域,开启摄像模式,调节CCD相机的焦距与光圈,使画面中标定板或直尺刻度显示达到最清晰,即使CCD相机对焦在片状激光所照射的垂直截面上。
步骤104:利用相机对激光器照射的垂直截面进行连续拍摄,生成测量图像。
其中,由计算机内置的可编程时间控制器调控激光器电源与CCD相机,连续拍摄照片, 生成测量图像。
步骤105:识别测量图像中的示踪粒子,根据测量图像中示踪粒子的移动距离和示踪粒子的拍摄时间生成垂直截面的流体速度分布。
其中,利用PIV图像数据处理***,计算出所拍图像瞬时速度的平均值,即为该下降管垂直截面的流体速度分布。
步骤106:对不同垂直截面中心高度处的流体速度进行拟合和积分,得到RH的循环流量。
其中,根据下降管不同垂直截面中心高度处的流体速度,得到下降管中心水平截面内若干点的流体速度,利用Matlab数学软件,采用多项式拟合流体速度曲面并得到拟合曲面函数,用该软件编写程序,计算下降管内流体速度的积分,进一步转换为RH的循环流量。
优选的,所述水盒子为有机玻璃制作的上端开口的长方体水槽,加装于CCD相机所拍摄流体区域的圆柱体外侧。
优选的,所述激光为532nm绿光,该激光峰值能量400mJ/Pulse。
优选的,所述PIV专用颗粒是直径50μm表面镀银的空心玻璃球。
优选的,所述CCD相机是指透过下降管外的水盒子进行拍摄。
优选的,所述下降管不同垂直截面是
从所选取的垂直截面按照一定间距向两侧选取下降管截面直到下降管边界。
通过在模拟装置中添加示踪粒子,并利用激光照射垂直截面,利用相机进行连续拍照,获取示踪粒子的移动距离和时间生成垂直截面的流体速度分布,经过拟合和积分得到RH的循环流量,不仅提高了截面的流体速度的测量精度,而且降低了对测量模型内流场的影响。
如图2所示,本发明的RH真空精炼循环流体模拟测量装置的模拟装置结构图,待测装置与模拟装置的模拟比例可以为1∶N,2≤N≤10。所述装置包括模拟装置和粒子成像测速装置;
所述模拟装置包括钢包1、真空室2、供气装置3、真空泵4、钢包水盒子5、上升管6、下降管7和下降管水盒子8;所述钢包1顶部与所述真空室2的底部连接,所述真空室2顶部与所述真空泵4链接,所述供气装置3与所述真空室2的侧壁连接,所述钢包水盒子5设置在所述钢包1的侧壁外侧,所述上升管6和下降管7穿过所述真空室2底部,所述上升管7管壁设置有通气孔,所述上升管6通过所述通气孔与所述供气装置3的输出端连接,所述下降管水盒子8设置在所述下降管7的管壁外侧。
优选的,所述模拟装置还包括:气体流量计9、分气盘10、U型压力计11;所述气体流量计9和分气盘10设置在所述上升管6与供气装置3之间,所述U型压力计11与所述真空室2连接。
所述粒子成像测速装置包括:激光器电源12、激光器13、导光臂14、片光源镜头组15、时间控制器16、数据处理器17和CCD相机18;所述激光器电源12与所述激光器13连接,所述激光器13的输出端与所述导光臂14连接,所述导光臂14的输出端与所述片光源镜头组15连接,所述数据处理器17与所述时间控制器16和CCD相机18连接。
优选的,所述数据处理器用于识别测量图像中的示踪粒子,根据测量图像中示踪粒子的移动距离和示踪粒子的拍摄时间生成测量点的瞬时速度的平均值,所述瞬时速度的平均值即为测量点的流体速度。
优选的,所述数据处理器还用于获取不同测量的流体速度,对不同垂直截面中心高度处的流体速度进行拟合和积分,得到测量位置的循环流量。
在使用时,将激光器对准模拟装置的下降管进行照射,CCD相机进行拍摄,并将生成的测量数据发送至数据处理器进行处理。
RH真空精炼循环流体模拟测量装置的优点为;
1.激光照射流体中的示踪粒子同时,CCD相机成像测速,在测量RH循环流量的过程中,避免了对RH装置内的流场产生干扰。
2.高分辨率CCD相机记录短时间内两次脉冲激光图像,具有很高的测量精度;圆柱形下降管外部加装水盒子,消除光折射对测速的影响,进一步提高测量精度。
3.对下降管中心横截面内的大量点进行速度,运用软件编程积分得到循环流量,处理结果可靠。
4.本发明对RH循环流量物理模拟的各项参数进行了优化,使得模拟过程更精准,测量精度高,误差小。
实施例1
(1)制作物理模型
钢包材质为有机玻璃,模型与原型按照1∶4的比例制作。钢包外为有机玻璃制作的钢包水盒子,测量时注满水,用以消除钢包曲面对光的折射。真空室为压力容器,其材质为有机玻璃,模型与原型的比例为1∶4。真空室的上部连接有真空泵,真空泵型号的选择应根据吹气量匹配。真空室的侧壁上用软管连接U型压力计。真空室的下降管外设有下降管水盒子,上升管圆周均匀错位排布两层共12个直径0.75mm的吹气孔。各个吹气孔由软管连接至分气盘,通过气体流量计与供气装置联通。
待测装置(原型)与模拟装置(模型)主要尺寸可以如表1所示。
表1 原型与模型尺寸(mm)对照表
激光器电源控制双腔激光器产生双脉冲Nd:YAG激光,该激光峰值能量400mJ/Pulse,为532nm绿光。经导光臂11传入片光源镜头组,脉冲激光束变换为片状光束。其中光束成型器的聚焦范围在0.3m—2m之间连续可调,在PIV试验过程中,可旋动该器件旋钮进行焦距的调节。下降管内PIV专用颗粒(直径50μm表面镀银的空心玻璃球)紧随流体运动,被激光照射后由双帧高分辨率CCD相机记录下来。ImagerPro11M型CCD相机,是逐行扫描相机,分辨率4008×2672像素,像素尺寸9.0μm×9.0μm,帧频5帧/秒。CCD相机配有10纳米带宽滤光片,用于接收示踪粒子的散射的532nm激光,滤掉杂散光,相机镜头焦距为50毫米,数值孔径为1.8。激光器电源与CCD相机由计算机内置的可编程时间控制器调控,时间分辨率可以达到10ns,控制信号的时序延迟抖动小于1ns。
(2)利用PIV标定RH真空室下降管某一垂直截面
在钢包内注入825mm水,启动真空泵使真空室内液面高度达到200mm,关闭真空泵,逆止阀使真空室内液面维持不变。启动PIV,待激光器电源预热完毕,调节片光源镜头组,使片状激光照射在下降管中心垂直截面。将直尺放置在下降管中心垂直面内,即片状激光照射区域。开启摄像模式,调节CCD相机的焦距与光圈,使画面最清晰,即使CCD相机对焦在所照射的中心垂直截面上。
(3)启动并调节RH物理模型使循环达到稳定
启动供气装置并调节吹气量,使气体流量计示数达到2.6m3/h;同时启动真空泵并调节抽气量,使U型压力计达到3615Pa,模型循环达到稳定。
(4)利用PIV测量下降管垂直截面的流体速度
用天平称取2.0gPIV专用颗粒,在水中搅拌、润湿后加入钢包内。由计算机内置的可编程时间控制器调控激光器电源与CCD相机,连续拍摄照片。利用PIV图像数据处理***,计算出所拍图像瞬时速度的平均值,即为该下降管垂直截面的流体速度分布。
(5)重复步骤(2)-(4),得到下降管不同垂直截面的流体速度分布
在下降管中心垂直截面两侧每隔2cm取一个截面,重复步骤(2)-(4)中相关操作,得到下降管不同垂直截面的流体速度分布。
(6)通过积分计算将所得数据转换为循环流量
根据下降管不同垂直截面中心高度处的流体速度,得到下降管中心水平截面内90个点的流体速度,如图4所示。利用Matlab数学软件中的sftool工具,采用5阶多项式拟合流体速度曲面,如图5所示,并得到拟合公式f=0.3013+4.522×10-17x+0.01848y-0.1397x2+5.853×10-18xy-0.103y2-1.068×10-16x3+0.1834x2y+2.893×10-16xy2-0.06864y3+0.002311x4-1.563×10-17x3y+0.01313x2y2-1.084×10-17xy3+0.01224y4+4.94×10-17x5-0.05735x4y-6.904×10-17x3y2-0.006353x2y3-1.042×10-16xy4+0.009138y5。
运用Matlab数学软件编写程序:
clearall
symsxy
f=0.3013+4.522×10-17x+0.01848y-0.1397x2+5.853×10-18xy-0.103y2-1.068×10-16x3+0.1834x2y+2.893×10-16xy2-0.06864y3+0.002311x4-1.563×10-17x3y+0.01313x2y2-1.084×10-17xy3+0.01224y4+4.94×10-17x5-0.05735x4y-6.904×10-17x3y2-0.006353x2y3-1.042×10-16xy4+0.009138y5;
vpa(int(int(f,y,-sqrt(0.0066015625-x^2),sqrt(0.0066015625-x^2)),x,-0.08125,0.08125))
得到下降管内流体速度的积分0.006241m3/s,转换为RH的循环流量为83.88t/min(原型)。
实施例2
采用本发明技术方案对某钢厂210tRH在浸入深度720mm、吹气量1600NL/min(原型参数)工况下的循环流量进行测定,步骤为:
(1)制作物理模型
模型与原型按照1:5的比例制作,其它试验装置与实施例1相同。
原型与模型主要尺寸如表2所示。
表2 原型与模型尺寸(mm)对照表
(2)利用PIV标定RH真空室下降管某一垂直截面
在钢包内注入660mm水,启动真空泵使真空室内液面高度达到160mm,关闭真空泵,逆止阀使真空室内液面维持不变。启动PIV,待激光器电源预热完毕,调节片光源镜头组,使片状激光照射在RH下降管中心垂直截面。将直尺放置在下降管中心垂直面内,即片状激光照射区域。开启摄像模式,调节CCD相机的焦距与光圈,使画面最清晰,即使CCD相机对焦在所照亮的中心垂直截面上。
(3)启动并调节RH物理模型使循环达到稳定
启动空气压缩机并调节吹气量,使气体流量计示数达到3.2m3/h;同时启动真空泵并调节抽气量,使U型压力计达到2895Pa,模型循环达到稳定。
(4)利用PIV测量下降管垂直截面的流体速度
用天平称取2.0gPIV专用颗粒,在水中搅拌、润湿后加入钢包内。由计算机内置的可编程时间控制器调控激光器电源与CCD相机,连续拍摄照片。利用PIV图像数据处理***,计算出所拍图像瞬时速度的平均值,即为该下降管垂直截面的流体速度分布。
(5)重复步骤(2)-(4),得到下降管不同垂直截面的流体速度分布
在下降管中心垂直截面两侧每隔2cm取一个截面,重复步骤(2)-(4)中相关操作,得到下降管不同垂直截面的流体速度分布。
(6)通过积分计算将所得数据转换为循环流量
根据下降管不同垂直截面中心高度处的流体速度,得到下降管中心水平截面内90个点的流体速度,如图4所示。利用Matlab数学软件中的sftool工具,采用5阶多项式拟合流体速度曲面,如图6所示,并得到拟合公式
f=0.4544-0.04378x+0.001745y-0.1221x2+0.002757xy-0.1121y2+0.03902x3+0.2238x2y+0.04369xy2-0.06796y3-0.008904x4+0.002749x3y+0.01419x2y2-0.002897xy3+0.01356y4-0.007962x5-0.07099x4y-0.01276x3y2-0.004875x2y3-0.00903xy4+0.009442y5。
运用Matlab数学软件编写程序:
clearall
symsxy
f=0.4544-0.04378x+0.001745y-0.1221x2+0.002757xy-0.1121y2+0.03902x3+0.2238x2y+0.04369xy2-0.06796y3-0.008904x4+0.002749x3y+0.01419x2y2-0.002897xy3+0.01356y4-0.007962x5-0.07099x4y-0.01276x3y2-0.004875x2y3-0.00903xy4+0.009442y5;
vpa(int(int(f,y,-sqrt(0.004225-x^2),sqrt(0.004225-x^2)),x,-0.065,0.065))
得到下降管内流体速度的积分0.006028m3/s,转换为RH的循环流量为141.53t/min(原型)。
实施例3
采用本发明技术方案研究在浸入深度400mm、吹气量1200NL/min(原型参数)工况下浸渍管内径对循环流量的影响,步骤为:
(1)制作物理模型
以某钢厂210tRH真空精炼装置为原型,根据相似原理模型与原型按照1:4建立水模型,其中下降管内径增大至170mm、外径增大至331mm,其它试验装置与实施例1相同。
原型与模型主要尺寸如表3所示。
表3 原型与模型尺寸(mm)对照表
(2)利用PIV标定RH真空室下降管某一垂直截面
在钢包内注入825mm水,启动真空泵使真空室内液面高度达到199.2mm,关闭真空泵,逆止阀使真空室内液面维持不变。启动PIV,待激光器电源预热完毕,调节片光源镜头组,使片状激光照射在下降管中心垂直截面。将直尺放置在下降管中心垂直面内,即片状激光照射区域。开启摄像模式,调节CCD相机的焦距与光圈,使画面最清晰,即使CCD相机对焦在所照射的中心垂直截面上。
(3)启动并调节RH物理模型使循环达到稳定
启动供气装置并调节吹气量,使气体流量计示数达到2.6m3/h;同时启动真空泵并调节抽气量,使U型压力计达到3615Pa,模型循环达到稳定。
(4)利用PIV测量下降管垂直截面的流体速度
用天平称取2.0gPIV专用颗粒,在水中搅拌、润湿后加入钢包内。由计算机内置的可编程时间控制器调控激光器电源与CCD相机,连续拍摄照片。利用PIV图像数据处理***,计算出所拍图像瞬时速度的平均值,即为该下降管垂直截面的流体速度分布。
(5)重复步骤(2)-(4),得到下降管不同垂直截面的流体速度分布
在下降管中心垂直截面两侧每隔2cm取一个截面,重复步骤(2)-(4)中相关操作,得到下降管不同垂直截面的流体速度分布。
(6)通过积分计算将所得数据转换为循环流量
根据下降管不同垂直截面中心高度处的流体速度,得到下降管中心水平截面内90个点的流体速度,如图4所示。利用Matlab数学软件中的sftool工具,采用5阶多项式拟合流体速度曲面,如图7所示,并得到拟合公式
f=0.3237+0.01414x+0.01628y-0.1307x2+0.002286xy-0.09918y2-0.005805x3+0.1759x2y-0.01543xy2-0.0658y3+0.001446x4-0.0009229x3y+0.01281x2y2-0.0004443xy3+0.0122yx-0.000301x5-0.05515x4y+0.004821x3y2-0.006364x2y3+0.003248xy4+0.008848y5。
运用Matlab数学软件编写程序:
clearall
symsxy
f=0.3237+0.01414x+0.01628y-0.1307x2+0.002286xy-0.09918y2-0.005805x3+0.1759x2y-0.01543xy2-0.0658y3+0.001446x4-0.0009229x3y+0.01281x2y2-0.0004443xy3+0.0122y4-0.000301x5-0.05515x4y+0.004821x3y2-0.006364x2y3+0.003248xy4+0.008848y5;
vpa(int(int(f,y,-sqrt(0.007225-x^2),sqrt(0.007225-x^2)),x,-0.085,0.085))
得到下降管内流体速度的积分0.00734m3/s,转换为RH的循环流量为98.65t/min(原型)。对比实施例1可知,增大浸渍管内径能提高循环流量。
实施例4
采用本发明技术方案研究在浸入深度400mm、吹气量1200NL/min(原型参数)工况下浸渍管内径对循环流量的影响,步骤为:
(1)制作物理模型
以某钢厂210tRH真空精炼装置为原型,根据相似原理模型与原型按照1:4建立水模型,其中下降管内径减小至150mm、外径减小至311mm,其它试验装置与实施例1相同。
原型与模型主要尺寸如表4所示。
表4原型与模型尺寸(mm)对照表
(2)利用PIV标定RH真空室下降管某一垂直截面
在钢包内注入825mm水,启动真空泵使真空室内液面高度达到201.4mm,关闭真空泵,逆止阀使真空室内液面维持不变。启动PIV,待激光器电源预热完毕,调节片光源镜头组,使片状激光照射在下降管中心垂直截面。将直尺放置在下降管中心垂直面内,即片状激光照射区域。开启摄像模式,调节CCD相机的焦距与光圈,使画面最清晰,即使CCD相机对焦在所照射的中心垂直截面上。
(3)启动并调节RH物理模型使循环达到稳定
启动供气装置并调节吹气量,使气体流量计示数达到2.6m3/h;同时启动真空泵并调节抽气量,使U型压力计达到3615Pa,模型循环达到稳定。
(4)利用PIV测量下降管垂直截面的流体速度
用天平称取2.0gPIV专用颗粒,在水中搅拌、润湿后加入钢包内。由计算机内置的可编程时间控制器调控激光器电源与CCD相机,连续拍摄照片。利用PIV图像数据处理***,计算出所拍图像瞬时速度的平均值,即为该下降管垂直截面的流体速度分布。
(5)重复步骤(2)-(4),得到下降管不同垂直截面的流体速度分布
在下降管中心垂直截面两侧每隔2cm取一个截面,重复步骤(2)-(4)中相关操作,得到下降管不同垂直截面的流体速度分布。
(6)通过积分计算将所得数据转换为循环流量
根据下降管不同垂直截面中心高度处的流体速度,得到下降管中心水平截面内90个点的流体速度,如图4所示。利用Matlab数学软件中的sftool工具,采用5阶多项式拟合流体速度曲面,如图8所示,并得到拟合公式
f=0.5038+0.3759x+0.3748y-74.42x2-1.323xy-72.79y2-109.3x3+3043x2y-234.8xy2-799.2y3-3156x4+198.6x3y+4136x2y2+517.1xy3+6221y4-3971x5-6.603×105x4y+5.239×104x3y2-6.276×104x2y3+2.906×104xy4+7.186×104y5。
运用Matlab数学软件编写程序:
clearall
symsxy
f=0.5038+0.3759x+0.3748y-74.42x2-1.323xy-72.79y2-109.3x3+3043x2y-234.8xy2-799.2y3-3156x4+198.6x3y+4136x2y2+517.1xy3+6221y4-3971x5-6.603×105x4y+5.239×104x3y2-6.276×104x2y3+2.906×104xy4+7.186×104y5
vpa(int(int(f,y,-sqrt(0.005625-x2),sqrt(0.005625-x2)),x,-0.075,0.075))
得到下降管内流体速度的积分0.005555m3/s,转换为RH的循环流量为74.66t/min(原型)。对比实施例1可知,减小浸渍管内径使循环流量随之减小。
上述实施例仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。
对比例1
将实施例1中的PIV专用颗粒16由直径50μm换为直径10μm,其他试验条件及操作与实施例1一致。
利用Matlab数学软件中的sftool工具,采用5阶多项式拟合流体速度曲面,并得到拟合公式f=0.4216-6.181×10-15x-0.3522y-80.05x2+1.704×10-14xy-75.72y2+5.388×10-12x3+3172x2y+2.422×10-12xy2-837.8y3-2881x4-5.508×10-12x3y+4270x2y2-1.674×10- 11xy3+6303y4-9.16×10-10x5-6.866×105x4y-5.685×10-10x3y2-6.265×104x2y3-4.7×10- 10xy4+7.422×104y5。
运用Matlab数学软件编写程序:
clearall
symsxy
f=0.4216-6.181×10-15x-0.3522y-80.05x2+1.704×10-14xy-75.72y2+5.388×10-12x3+3172x2y+2.422×10-12xy2-837.8y3-2881x4-5.508×10-12x3y+4270x2y2-1.674×10- 11xy3+6303y4-9.16×10-10x5-6.866×105x4y-5.685×10-10x3y2-6.265×104x2y3-4.7×10- 10xy4+7.422×104y5
vpa(int(int(f,y,-sqrt(0.0066015625-x2),sqrt(0.0066015625-x2)),x,-0.08125,0.08125))
得到下降管内流体速度的积分0.003959m3/s,转换为RH的循环流量为53.21t/min(原型)。
由实施例1与对比例1发现,当PIV专用颗粒直径由50μm减小到10μm后,所测循环流量也随之减小。因为PIV颗粒越小,越难捕捉其运动状态,所测量的流体速度偏小,导致循环流量变小。所以采用直径50μm的PIV专用颗粒试验能获得更好的试验结果。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种炼钢过程中真空精炼物理模拟试验方法,其特征在于,包括:
步骤101:根据待测RH真空精炼装置建立模拟装置;
其中,采用有机玻璃按照模型与原型1:N的比例制作RH物理模型,其中2≤N≤10,并在真空室下降管外部及钢包外部加装水盒子;
步骤102:运行模拟装置,并在模拟装置中添加示踪粒子;
其中,将适量PIV专用颗粒在水中搅拌、润湿后加入钢包,调节供气装置与真空泵以控制吹气量与抽气量,使RH物理模型循环达到稳定;
步骤103:在模拟装置的测量位置选取垂直截面,利用激光器照射所述垂直截面;
其中,在钢包中加入预定量的水,调节真空泵使真空室内压强达到原型对应值,水盒子内充满水,启动PIV,待激光器电源预热完毕,调节片光源镜头组,使片状激光照射在RH下降管某一垂直截面,将标定板或直尺放置在该垂直面内,即片状激光照射区域,开启摄像模式,调节CCD相机的焦距与光圈,使画面中标定板或直尺刻度显示达到最清晰,即使CCD相机对焦在片状激光所照射的垂直截面上;
步骤104:利用相机对激光器照射的垂直截面进行连续拍摄,生成测量图像;
其中,由计算机内置的可编程时间控制器调控激光器电源与CCD相机,连续拍摄照片,生成测量图像;
步骤105:识别测量图像中的示踪粒子,根据测量图像中示踪粒子的移动距离和示踪粒子的拍摄时间生成垂直截面的流体速度分布;
其中,利用PIV图像数据处理***,计算出所拍图像瞬时速度的平均值,即为该下降管垂直截面的流体速度分布;
步骤106:对不同垂直截面中心高度处的流体速度进行拟合和积分,得到RH的循环流量;
其中,根据下降管不同垂直截面中心高度处的流体速度,得到下降管中心水平截面内若干点的流体速度,利用Matlab数学软件,采用多项式拟合流体速度曲面并得到拟合曲面函数,用该软件编写程序,计算下降管内流体速度的积分,进一步转换为RH的循环流量;
其中,所述水盒子为有机玻璃制作的上端开口的长方体水槽,加装于CCD相机所拍摄流体区域的圆柱体外侧;
其中,所述激光为532nm绿光,该激光峰值能量400mJ/Pulse;
其中,所述PIV专用颗粒是直径50μm的空心玻璃球;
其中,所述CCD相机是透过下降管外的水盒子进行拍摄;
其中,所述下降管不同垂直截面是从所选取的垂直截面按照一定间距向两侧选取下降管截面直到下降管边界;
其中,所述方法采用了如下物理模拟试验模拟装置,包括模拟装置和粒子成像测速装置;所述模拟装置包括钢包,真空室,供气装置,真空泵,U型压力计,钢包水盒子,上升管,下降管,下降管水盒子;所述钢包顶部与所述真空室的底部连接,所述真空室顶部与所述真空泵连接,所述供气装置与所述真空室的侧壁连接,所述钢包水盒子设置在所述钢包的侧壁外侧,所述上升管和下降管穿过所述真空室底部,所述上升管管壁设置有通气孔,所述上升管通过所述通气孔与所述供气装置的输出端连接,所述下降管水盒子设置在所述下降管的管壁外侧;
所述粒子成像测速装置包括:激光器电源,双腔激光器,导光臂,计算机内置可编程时间控制器,PIV图像数据处理***,CCD相机,片光源镜头组,PIV专用颗粒;所述激光器电源与所述激光器连接,所述激光器的输出端与所述导光臂连接,所述导光臂的输出端与所述片光源镜头组连接,所述PIV图像数据处理***与所述计算机内置可编程时间控制器和所述CCD相机连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410426265.8A CN104164537B (zh) | 2014-08-26 | 2014-08-26 | 一种炼钢过程中真空精炼物理模拟试验方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410426265.8A CN104164537B (zh) | 2014-08-26 | 2014-08-26 | 一种炼钢过程中真空精炼物理模拟试验方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104164537A CN104164537A (zh) | 2014-11-26 |
CN104164537B true CN104164537B (zh) | 2016-05-25 |
Family
ID=51908518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410426265.8A Active CN104164537B (zh) | 2014-08-26 | 2014-08-26 | 一种炼钢过程中真空精炼物理模拟试验方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104164537B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104458522B (zh) * | 2014-11-27 | 2017-05-03 | 中盟科创(深圳)科技发展有限公司 | 基于移动终端的大气微粒检测***与检测方法 |
CN104777164A (zh) * | 2015-03-30 | 2015-07-15 | 天津大学 | 基于piv的机舱内空气流动大尺寸测量实验装置及方法 |
CN106153294A (zh) * | 2016-07-21 | 2016-11-23 | 西安建筑科技大学 | 一种大型钢锭多炉合浇的水力学模拟装置 |
EP3550035A4 (en) * | 2017-02-27 | 2020-04-22 | Nippon Steel Corporation | DIAGNOSTIC ASSISTANCE DEVICE FOR VACUUM DEGASSING TANK, DIAGNOSTIC ASSISTANCE METHOD, DIAGNOSTIC METHOD, AND REPAIR METHOD |
CN108315529B (zh) * | 2018-04-21 | 2023-12-19 | 江西理工大学 | 一种rh精炼炉上升管吹气孔布置结构 |
CN108765578B (zh) * | 2018-05-15 | 2021-02-12 | 华中科技大学 | 离心铸造物理模拟的充型流体形貌三维拍摄***及方法 |
CN111289564B (zh) * | 2020-04-07 | 2022-05-06 | 东北石油大学 | 储油罐热油喷洒加热过程传热和流动规律模拟装置及方法 |
CN111289563B (zh) * | 2020-04-07 | 2022-05-27 | 东北石油大学 | 研究储油罐管式加热过程传热与流动规律的装置及方法 |
CN111521841B (zh) * | 2020-05-22 | 2021-06-15 | 北京科技大学 | 用于结晶器水口模型内壁流场粒子图像测速的装置及方法 |
CN114002457B (zh) * | 2021-11-17 | 2023-05-09 | 北京航空航天大学 | 粒子图像测速装置 |
CN114134288B (zh) * | 2021-12-03 | 2023-04-07 | 北京科技大学 | 一种用于rh浸渍管的热状态及破损监测***及方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102965475A (zh) * | 2012-11-29 | 2013-03-13 | 重庆大学 | Rh精炼装置物理模型制作方法 |
-
2014
- 2014-08-26 CN CN201410426265.8A patent/CN104164537B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102965475A (zh) * | 2012-11-29 | 2013-03-13 | 重庆大学 | Rh精炼装置物理模型制作方法 |
Non-Patent Citations (5)
Title |
---|
Katoh T,OkamotoT.Mixing of molten steel in a ladle with RH reactor by water model experiment.《Electric Furnace Steel》.1979,第50卷(第2期),第128-137页. * |
朱德平等.真空循环精炼过程中钢液的流动和混合特性.《包头钢铁学院学报》.2001,第20卷(第1期),第12-18页. * |
肖亚克等.PIV技术中的示踪粒子发生和布撤.《气动研究与实验》.2005,第22卷(第4期),第8-15页. * |
舒宏富等.RH-MFB真空精炼过程中循环流量的物理模拟研究.《材料与冶金学报》.2004,第3卷(第2期),第107-112页. * |
郁能文等.RH精炼过程中钢液流动特性的水模型.《东北大学学报(自然科学版)》.1998,第19卷第118-120页. * |
Also Published As
Publication number | Publication date |
---|---|
CN104164537A (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104164537B (zh) | 一种炼钢过程中真空精炼物理模拟试验方法及装置 | |
CN101852814B (zh) | 一种滴灌灌水器迷宫流道内流动的全场测试方法 | |
CN105738648B (zh) | 一种多相体系中颗粒速度的在线测量方法 | |
CN202599852U (zh) | 气液两相泡状流中气泡的识别及测量装置 | |
CN113324727B (zh) | 一种针对压缩拐角超声速流场结构的纹影图像处理方法 | |
CN214472543U (zh) | 一种液体粘度检测*** | |
CN103439230A (zh) | 一种气泡参数测量方法及测量装置 | |
CN101393107B (zh) | 一种微气泡动态显微测试装置 | |
CN110044278B (zh) | 一种非接触式流体薄膜厚度测量装置及方法 | |
CN105928841B (zh) | 一种浸入式在线多相测量仪及测量方法 | |
CN102680076B (zh) | 一种模拟人体声带振动装置及其实现方法 | |
CN109372034A (zh) | 上拔过程中吸力桶基础内部破坏机理的试验装置及方法 | |
CN114061964A (zh) | 一种多功能雾化试验*** | |
CN111220514B (zh) | 一种水下微气泡及其絮体测量装置 | |
CN107588885A (zh) | 一种仿生鱼摆尾的压力场测量装置与方法 | |
CN102680203A (zh) | 微小通道气液两相流空隙率测量装置及方法 | |
CN109141814A (zh) | 一种基于piv技术测量超重力旋转床流场特性的装置和方法 | |
CN108007669A (zh) | 一种实时追踪上浮气泡动态特性和流场特征的装置 | |
CN102564905B (zh) | 超导磁体模拟微重力环境下接触角测量装置及其方法 | |
CN108844959B (zh) | 一种圆管内气液两相环状流截面相含率的测量及修正方法 | |
CN109323845A (zh) | 基于高速成像技术的雾滴发生装置 | |
CN106984785B (zh) | 一种用3d摄像或照相来监测结晶器内液面波动的方法 | |
Wenyuan et al. | Study on the flow field around two parallel moving bubbles and interaction between bubbles rising in CMC solutions by PIV | |
CN112611680A (zh) | 一种液体粘度检测装置 | |
CN113624446B (zh) | 一种非接触式分层流场测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |