CN105033214A - 一种宽厚板坯连铸机基础辊缝制定方法 - Google Patents
一种宽厚板坯连铸机基础辊缝制定方法 Download PDFInfo
- Publication number
- CN105033214A CN105033214A CN201510541805.1A CN201510541805A CN105033214A CN 105033214 A CN105033214 A CN 105033214A CN 201510541805 A CN201510541805 A CN 201510541805A CN 105033214 A CN105033214 A CN 105033214A
- Authority
- CN
- China
- Prior art keywords
- casting
- dimensional
- finite element
- element model
- mrow
- 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
- 238000005266 casting Methods 0.000 title claims abstract description 274
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000007711 solidification Methods 0.000 claims abstract description 74
- 230000008023 solidification Effects 0.000 claims abstract description 74
- 238000012546 transfer Methods 0.000 claims abstract description 56
- 238000001816 cooling Methods 0.000 claims abstract description 54
- 230000008602 contraction Effects 0.000 claims abstract description 25
- 238000009749 continuous casting Methods 0.000 claims description 39
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 30
- 230000009467 reduction Effects 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 14
- 229910000831 Steel Inorganic materials 0.000 claims description 11
- 239000010959 steel Substances 0.000 claims description 11
- 239000007790 solid phase Substances 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 230000005855 radiation Effects 0.000 claims description 5
- 230000017525 heat dissipation Effects 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 13
- 230000005499 meniscus Effects 0.000 description 8
- 239000007787 solid Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 230000004907 flux Effects 0.000 description 4
- 238000005204 segregation Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000000498 cooling water Substances 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000005452 bending Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000011946 reduction process Methods 0.000 description 2
- 238000005507 spraying Methods 0.000 description 2
- 239000002352 surface water Substances 0.000 description 2
- 239000010963 304 stainless steel Substances 0.000 description 1
- 229910000975 Carbon steel Inorganic materials 0.000 description 1
- 229910000589 SAE 304 stainless steel Inorganic materials 0.000 description 1
- 241000860832 Yoda Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000010962 carbon steel Substances 0.000 description 1
- 238000005097 cold rolling Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000005098 hot rolling Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
Landscapes
- Continuous Casting (AREA)
Abstract
本发明涉及一种宽厚板坯连铸机基础辊缝制定方法,该方法为采集连铸机的浇铸工况参数,建立板坯的凝固传热三维有限元模型,根据连铸机的浇铸工况参数计算得到浇铸过程中的三维温度场;建立板坯的三维热收缩有限元模型,以浇铸过程中的三维温度场作为温度载荷求解该三维热收缩有限元模型,得到三维热收缩规律,根据三维热收缩规律及板坯的初始厚度制定相应工况条件下的连铸机基础辊缝;本发明以二冷区内的实测水流密度作为边界条件,保证了凝固传热及热收缩规律计算的准确性;根据板坯在凝固传热及不同厚度方向热收缩存在非均匀性特点,结合铸坯的初始厚度求得各扇形段基础辊缝,能够更加准确的计算铸坯不同位置及不同浇铸时刻的热收缩情况。
Description
技术领域
本发明属于板坯连铸领域,具体涉及一种宽厚板坯连铸机基础辊缝制定方法。
背景技术
板坯连铸机主要由结晶器及后面由不同类型密排分节辊构成的扇形段组成。板坯浇铸过程中,钢水首先通过结晶器一次冷却形成一定坯壳厚度的铸坯,再经过后续扇形段中的喷淋二次冷却及弯曲矫直过程,最终完全冷却凝固形成板坯。板坯经过后续的热轧或冷轧过程成为最终的板材,被广泛应用于国民生产的各个领域,例如,船舶制造、桥梁工程、电站锅炉、重型机械、海洋石油及军事领域。
板坯连铸机生产的板坯质量直接影响最终轧材质量。在板坯浇铸过程中,经常出现的铸坯质量问题包括中心偏析、疏松及各种类型的裂纹等。通过一定的方法制定合理的基础辊缝,并在此基础上,实施凝固末端压下是改善上述质量问题的有效途径之一。
板坯自然热收缩是指板坯浇铸过程中,由于铸坯逐渐受冷凝固而产生的自然收缩现象。基于板坯浇铸过程中板坯的自然热收缩行为制定合理的基础辊缝对避免铸坯鼓肚、降低矫直应力、防止滞坯等具有重要意义,同时也是凝固末端压下工艺实施的重要前提保障。
目前针对连铸机辊缝设计及铸坯热收缩计算方法已有所报道。中国专利CN101362196公开了一种基于大方坯自然热收缩辊缝工艺的控制方法,通过建立大方坯的二维在线模型,实时计算铸坯的二维温度场,并在此基础上,通过热-力直接耦合的方式来计算大方坯的二维热收缩规律,最终根据热收缩规律计算出铸坯实时控制辊缝值;曹学欠(连铸大方坯热收缩行为的有限元分析,连铸,2010,6:1-4),林启勇(连铸板坯自然热收缩行为有限元模拟,东北大学学报(自然科学版),27(S2):8-10)通过建立二维切片有限元模型,分别研究了大方坯及板坯浇铸过程中热收缩规律;陈洪智(不锈钢板坯连铸自由线收缩与辊缝研究,中国冶金,22(2):25-30)通过二维差分方法,研究了304不锈钢板坯浇铸过程中的线收缩规律。上述计算方法均采用了二维切片模型进行建模,忽略了沿拉坯方向上的凝固传热及热收缩过程中铸坯变形方面的相互影响。然而,在铸流比较靠前的位置,铸坯受到的冷却强度大,沿拉坯方向上温度降低速度比较快,在此阶段,铸坯沿拉坯方向上的凝固传热及铸坯热收缩变形相互影响会比较显著,因此,采用二维建模势必会影响模拟结果的准确性。此外,上述方法在二冷区各区设置均匀的等效换热系数作为边界条件,然而,在板坯浇铸过程中,二冷区采用多喷嘴的交错布置方式,铸坯表面受到的水流密度分布趋势比较复杂,如果单纯采用均匀等效换热系数作为边界条件,势必引起计算结果的偏差,无法真实反应连铸坯宽向的非均匀凝固过程及其对辊缝制定的影响作用。
历英(连铸机辊列辊缝的优化控制,中国钢铁年会论文集,2009:542-545)及高永龙(钢水凝固收缩与板坯连铸机开口度参数设计,山东冶金,2012,34(1):24-26)结合由物理实验得到的钢液凝固过程各阶段的收缩系数,估算了浇铸过程中铸坯的热收缩量,并据此制定了连铸机静态压下辊缝。然而此方法存在如下不足:将连铸过程分为几个凝固阶段,只能粗略的计算整个浇铸过程中铸坯总收缩量。而有限元或有限差分的方法将整个连铸过程进行空间离散及时间离散,能够更加准确的计算铸坯不同位置及不同浇铸时刻的热收缩情况,因此,该方法计算的准确性不及有限元或有限差分等数值模拟法。
发明内容
针对现有技术的不足,本发明提出一种宽厚板坯连铸机基础辊缝制定方法,通过三维建模方法建立板坯浇铸过程中的凝固传热三维有限元模型和三维热收缩有限元模型,模拟计算出板坯的三维温度场及三维热收缩规律,在此基础上,综合考虑板坯的凝固传热及表面收缩的不均匀性,最终制定板坯浇铸过程中的基础辊缝。
一种宽厚板坯连铸机基础辊缝制定方法,包括以下步骤:
步骤1:采集板坯连铸机的浇铸工况参数,所述浇铸工况参数包括:浇铸钢种的断面尺寸、浇铸钢种成分、浇铸温度、工作拉速、结晶器水流量及回水温差、二冷区划分参数及各二冷区内的水量;
步骤2:建立板坯的凝固传热三维有限元模型,根据连铸机的浇铸工况参数计算得到浇铸过程中的三维温度场;
步骤2.1:取板坯浇铸过程中一定高度的1/4铸坯,建立其凝固传热三维有限元模型;
步骤2.2:根据连铸机的浇铸工况参数确定凝固传热三维有限元模型的不同铸流位置处的边界条件,求解凝固传热三维有限元模型,得到浇铸过程中的三维温度场;
步骤2.3:通过实测板坯的表面温度及坯壳厚度,对凝固传热三维有限元模型进行修正,从而对板坯浇铸过程中的三维温度场进行修正;
步骤3:建立板坯的三维热收缩有限元模型,以浇铸过程中的三维温度场作为温度载荷求解该三维热收缩有限元模型,得到浇铸过程中的三维热收缩规律;
步骤3.1:取板坯浇铸过程中一定高度的1/4铸坯,建立其三维热收缩有限元模型;
步骤3.2:将浇铸过程中的三维温度场作为温度载荷,通过刚体对称面完成对该有限元模型宽度方向和厚度方向的位移约束,求解该三维热收缩有限元模型,得到浇铸过程中的三维热收缩规律;
步骤4:根据获取的浇铸过程中的三维热收缩规律及板坯的初始厚度制定相应工况条件下的连铸机基础辊缝;
步骤4.1:根据浇铸过程中的三维温度场和实施轻压下时压下区间的起始固相率确定铸坯压下起始铸流位置和铸坯凝固终点铸流位置;
步骤4.2:根据铸坯压下起始铸流位置和铸坯凝固终点铸流位置确定当前铸流位置的域加权值;
步骤4.3:根据浇铸过程中的三维热收缩规律、板坯的初始厚度、以及当前铸流所处浇铸过程中的不同位置的域加权值,制定相应工况条件下的连铸机基础辊缝。
所述的制定相应工况条件下的连铸机基础辊缝的计算公式如下:
其中,Gap为基础辊缝值,InGap为铸坯初始厚度,Sb为已凝固区域起始节点,Se为已凝固区域结束节点,Mb为未凝固区域起始节点,Me为未凝固区域结束节点,i为凝固传热三维有限元模型中已凝固区域的节点,j为凝固传热三维有限元模型中未凝固区域的节点,Si为已凝固区域内的节点i的厚度方向收缩量,Sj为未凝固区域内的节点j的厚度方向收缩量,为已凝固区域内节点i的热收缩加权值,N1为凝固传热三维有限元模型中已凝固区域的节点个数,为未凝固区域内的节点j的热收缩加权值,N2为凝固传热三维有限元模型中未凝固区域的节点个数,γ1为已凝固区域的域加权值,γ2为未凝固区域的域加权值。
所述的凝固传热三维有限元模型的不同铸流位置处的边界条件分别为:以铸坯在结晶器中的热流密度作为结晶器的边界条件,以实测水流密度确定的等效换热系数作为二冷区边界条件,以空冷区内的辐射散热作为空冷区的边界条件。
所述的根据铸坯压下起始铸流位置和铸坯凝固终点铸流位置确定当前铸流位置的域加权值具体方法为:
(1)当前铸流位置小于铸坯压下起始铸流位置时,已凝固区域的域加权值取值为0≤γ1≤0.5,未凝固区域的域加权值取值为0.5≤γ2≤1.0;
(2)当前铸流位置介于铸坯压下起始铸流位置和铸坯凝固终点铸流位置之间时,已凝固区域的域加权值取值为0.3≤γ1≤0.8,未凝固区域的域加权值取值为0.2≤γ2≤0.7;
(3)当前铸流位置大于铸坯凝固终点铸流位置时,已凝固区域的域加权值取值为γ1=1.0,未凝固区域消失,则γ2=0。
本发明的有益效果:
本发明提出一种宽厚板坯连铸机基础辊缝制定方法,与现有技术方法相比,本发明具有如下特点:1)以板坯浇铸过程中的一段铸坯为对象,建立了凝固传热三维有限元模型和三维热收缩有限元模型,并以二冷区内的实测水流密度作为边界条件,保证了凝固传热及热收缩规律计算的准确性;2)根据板坯在凝固传热及不同宽面位置处的厚度方向热收缩存在非均匀性特点,首先准确描述板连铸坯非均匀凝固的形貌特点,然后通过加权的方式综合考虑已凝固区域与未凝固区域对铸坯热收缩的整体影响作用,求得铸坯整体的热收缩值,最后结合铸坯的初始厚度最终求得各扇形段基础辊缝,能够更加准确的计算铸坯不同位置及不同浇铸时刻的热收缩情况。
附图说明
图1为本发明具体实施方式中宽厚板坯连铸机基础辊缝制定方法的流程图;
图2为本发明具体实施方式中水压、气压均为0.1MPa条件下的三喷嘴组合水流密度分布趋势;
图3为本发明具体实施方式中计算浇铸过程中的三维温度场的流程图;
图4为本发明具体实施方式中板坯三维有限元模型示意图;
其中,I为厚度方向的刚体对称面,II为宽度方向的刚体对称面,III为铸坯窄面,IV为铸坯宽面,V为弯月面位置;
图5为本发明具体实施方式中板坯的三维凝固传热模型的计算结果与实测值之间的对比图;
图6为本发明具体实施方式中计算浇铸过程中的三维热收缩规律的流程图;
图7为本发明具体实施方式中三维热收缩模型在二冷8区末端位置处的宽面不同位置处的热收缩值;
图8为本发明具体实施方式中凝固传热三维有限元模型中已凝固区域和未凝固区域的划分示意图;
其中,a为未完全凝固区域,b为完全凝固区域。
具体实施方式
下面结合附图对本发明具体实施方式加以详细的说明。
本发明具体实施方式中结合现场板坯连铸机的实际浇铸工况参数,以大型商用有限元软件MSC.Marc为计算工具,通过建立板坯浇铸过程中的凝固传热三维有限元模型和三维热收缩有限元模型,最终获取了相应工况条件下板坯三维热收缩规律,基于上述规律制定了相应工况条件下的板坯连铸机基础辊缝。
一种宽厚板坯连铸机基础辊缝制定方法,如图1所示,包括以下步骤:
步骤1:采集板坯连铸机的浇铸工况参数,所述浇铸工况参数包括:浇铸钢种的断面尺寸、浇铸钢种成分、浇铸温度、工作拉速、结晶器水流量及回水温差、二冷区划分参数及各二冷区内的水量。
本实施方式中,板坯连铸机参数:浇铸温度为1540℃,结晶器高度900mm,有效高度800mm,结晶器宽面水流量为5029L/min,回水温差为7.0℃,结晶器窄面水流量为584L/min,回水温差为7.8℃,铸机包含8个二冷区和2个空冷区,其中二冷区长度20.57m,总冷却长度34.725m,二冷区划分参数及各二冷区内的水量如表1所示;板坯连铸机包含15个扇形段,能够浇铸铸坯的宽面尺寸为1600mm、1800mm及2100mm,厚度尺寸为180mm、230mm及250mm,本实施方式中的铸坯尺寸为2100mm×250mm。
表1二冷区划分参数及各二冷区内的水量
浇铸钢种AH36,钢种主要成分百分含量为:0.15C-0.25Si-1.36Mn-0.025P-0.010S;浇铸钢种的断面尺寸为2100mm×250mm;现场工作拉速0.9m/min;二冷区比水量0.86L/Kg。
板坯连铸机的浇铸工况参数采集的准确性直接影响后续的凝固传热模型及热收缩模型边界条件施加的准确性,并最终关系到模拟结果的准确性及制定的基础辊缝的合理性。
本实施方式中,由于在实际生产过程中,仅能采集到二冷各区的冷却回路总水量、气量及各区的喷嘴布置方式及每个喷嘴的喷淋参数,而不能直接掌握浇铸过程中铸坯表面不同位置处的受水密度。因此,为了掌握宽厚板坯浇铸过程中,板坯宽面水流密度分布趋势,还需结合二冷区参数,通过冷态实验方法实测获取板坯表面的二冷区水流密度分布。当二冷区喷嘴采用三喷嘴错列布置且气压水压均为0.1MPa时,测得的三喷嘴组合水流密度分布趋势如图2所示。
步骤2:建立板坯的凝固传热三维有限元模型,根据连铸机的浇铸工况参数计算得到浇铸过程中的三维温度场,如图3所示。
步骤2.1:取板坯浇铸过程中一定高度的1/4铸坯,建立其凝固传热三维有限元模型。
本实施方式中,板坯三维有限元模型示意图如图4所示,其中,x,y,z为有限元模型所在坐标系中的三个坐标轴,I为厚度方向的刚体对称面,用来约束位于该平面上的所有节点的x方向位移,II为宽度方向的刚体对称面,用来约束位于该平面上的所有节点的z方向位移,III为铸坯窄面;IV为铸坯宽面;V为弯月面位置。
根据板坯浇铸过程中板坯冷却具有的对称性(宽度方向及厚度方向对称),取板坯浇铸过程中高度为800m的1/4铸坯,板坯凝固传热三维有限元模型如式(1)所示:
其中,T为铸坯温度,ρ为铸坯密度,C为铸坯比热容,keff为导热系数,So为内热源项,v为工作拉速,x,y,z为铸坯的坐标,t为浇铸时间。
凝固传热三维有限元模型初始温度为现场浇铸温度T0=1540℃。
步骤2.2:根据连铸机的浇铸工况参数确定凝固传热三维有限元模型的不同铸流位置处的边界条件,求解凝固传热三维有限元模型,得到浇铸过程中的三维温度场。
凝固传热三维有限元模型的不同铸流位置处的边界条件分别为:
(一)以铸坯在结晶器中的热流密度作为结晶器的边界条件。
本实施方式中,采用Savage、Pritchard的方法测得的铸坯在结晶器中的热流密度作为边界条件,如式(2)所示:
其中,q(y)为凝固传热三维有限元模型距弯月面距离为y位置处的结晶器热流密度,MW/m2;l为凝固传热三维有限元模型每个单元距弯月面距离,m,v为工作拉速,m/s,A一般取经验值2.64,B为根据现场浇铸时结晶器水流量及回水温差求得,公式如式(3)和式(4)所示:
其中,为结晶器内平均热流密度,MW/m2,ρw为冷却水密度,kg/m3,Vw为冷却水流量,m3/s;Cw为冷却水热容,4200J/(kg·℃),Δt为结晶器内的回水温差,As为结晶器传热面的面积,m2,tmold为结晶器内的冷却时间,可根据结晶器的有效高度h及工作拉速v计算得到,tmold=h/v,q(t)为铸坯进入结晶器内不同时刻的热流密度,MW/m2。
通过联立上述两式,最终求得B的表达式如式(5)所示:
本实施方式中,结合结晶器相关工况参数及式(5)最终制定的A、B为:
结晶器宽面:A=2.688;B=0.2628,
结晶器窄面:A=2.688;B=0.2288。
(二)以实测二冷区水流密度确定的等效换热系数作为二冷区边界条件。
本实施方式中,采用菲格洛及岛田测定的二冷区等效换热系数作为二冷区内的边界条件,如式(6)所示:
hi=αiwi (z) 0.55(1-0.0075Tw)(6)
其中,hi为第i个二冷区的等效换热系数,W/(m2·℃);w(z)为实测第i个二冷区宽面位置z处的水流密度,L/(m2·s),该水流密度根据相应二冷区内的喷嘴及水流量经过实际测量得出;αi-第i个二冷区内的修正系数。
本实施方式中,得到的二冷区的等效换热系数如表2所示:
表2二冷区的等效换热系数
(三)以空冷区内的辐射散热作为空冷区的边界条件。
本实施方式中,空冷区内的辐射散热如式(7)所示:
qB=σε((T+273)4-(Tamb+273)4)(7)
式中,qB为空冷区铸坯表面的辐热散热热流密度,w/m2;σ=5.67×10-8W/(m2·K4)为斯蒂芬-波尔兹曼常数,ε=0.8为辐射系数;Tsur为铸坯表面温度,℃;Tamb为环境温度,℃。
本实施方式中,为了避免已进入结晶器部分的凝固传热三维有限元模型在计算过程中对未进入结晶器部分模型在凝固传热及热收缩变形方面的影响,通过MSC.Marc中的二次开发子程序Uactive对上述过程进行处理,具体过程为:通过Uactive子程序判断铸坯有限元模型中的每个单元与弯月面的位置关系,当单元处于弯月面以上时,说明单元还未进入结晶器进行冷却降温及变形过程,因此,维持该单元的初始状态;当有限元模型中的单元处于结晶器弯月面以下时,该单元正常参与冷却降温及热收缩变形计算过程。此外,借助MSC.Marc中的二次开发子程序Uflux及Ufilm完成凝固传热计算过程中边界条件的施加,具体方法:以铸坯宽面及窄面单元外表面上的每个积分点为单位,实时计算每个积分点到弯月面的距离,根据该距离及步骤1中采集的现场铸机二冷区划分参数,对每个积分点施加相应铸流位置处的冷却边界条件。
步骤2.3:通过实测板坯的表面温度及坯壳厚度,对凝固传热三维有限元模型进行修正,从而对板坯浇铸过程中的三维温度场进行修正。
本实施方式中,为了验证与校正板坯的三维凝固传热模型的准确性,在一定的铸流位置采用红外热成像仪(ThermalCAMTMResearcher)测定板坯的表面温度,通过射钉实验测量板坯相应铸流位置处的坯壳厚度,通过实测数据对板坯的三维凝固传热模型进行校正,最终得出的板坯的三维凝固传热模型的计算结果与实测值之间的对比如图5所示,根据附图说明了板坯的三维凝固传热模型的准确性。
步骤3:建立板坯的三维热收缩有限元模型,以浇铸过程中的三维温度场作为温度载荷求解该三维热收缩有限元模型,得到浇铸过程中的三维热收缩规律,如图6所示。
步骤3.1:取板坯浇铸过程中一定高度的1/4铸坯,建立其三维热收缩有限元模型。
本实施方式中,同样取板浇铸过程中高度为800m的1/4铸坯,建立其三维热收缩有限元模型如式(8)和式(9)所示:
d{ε}=d{ε}e+d{ε}p+d{ε}T=d{ε}e+d{ε}p+{α}dT(8)
d{σ}=[D]ep(d{ε}-d{ε0})+d{σ0}(9)
其中,dT为三维热收缩有限元模型的温度增量,d{ε}为三维热收缩有限元模型的总应变增量;d{ε}e为三维热收缩有限元模型的弹性应变增量;d{ε}p为三维热收缩有限元模型的塑性应变增量;d{ε}T为三维热收缩有限元模型的热应变增量;{α}为三维热收缩有限元模型的热收缩系数;d{σ}为三维热收缩有限元模型的总应力增量;[D]ep为三维热收缩有限元模型的刚度矩阵;d{ε0}为三维热收缩有限元模型的初始应变增量;d{σ0}为三维热收缩有限元模型的初始应力增量。
三维热收缩有限元模型的初始温度为现场浇铸温度T0=1540℃;初始应力应变均为0。
步骤3.2:将浇铸过程中的三维温度场作为温度载荷,通过刚体对称面完成对该有限元模型宽度方向和厚度方向的位移约束,求解该三维热收缩有限元模型,得到浇铸过程中的三维热收缩规律。
三维热收缩模型在二冷8区末端位置处的宽面不同位置处的热收缩值如图7所示。
步骤4:根据获取的浇铸过程中的三维热收缩规律及板坯的初始厚度制定相应工况条件下的连铸机基础辊缝。
步骤4.1:根据铸过程中的三维温度场和实施轻压下时压下区间的起始固相率确定铸坯压下起始铸流位置Pstart和铸坯凝固终点铸流位置Psolid。
本实施方式中,通过浇铸过程中的三维温度场确定铸坯在不同铸流位置的固相率,当计算得到的固相率与现场实施压下工艺的起始固相率相等时,此时的铸流位置即为铸坯压下起始铸流位置Pstart,当计算得到的固相率为1时所对应的铸流位置即为铸坯凝固终点铸流位置Psolid。
步骤4.2:根据铸坯压下起始铸流位置和铸坯凝固终点铸流位置确定当前铸流位置的域加权值。
根据铸坯压下起始铸流位置Pstart和铸坯凝固终点铸流位置Psolid确定当前铸流位置的域加权值具体方法为:
(1)当前铸流位置小于铸坯压下起始铸流位置即P<Pstart时,已凝固区域的域加权值取值为0≤γ1≤0.5,未凝固区域的域加权值取值为0.5≤γ2≤1.0。
本实施方式中,将当前铸流位置小于铸坯压下起始铸流位置即P<Pstart定义为第一阶段,在该阶段,铸坯表面温度较高,坯壳比较薄,受铸坯中心未凝固区域钢水静压力作用,铸坯易在相应区域产生鼓肚缺陷,且在铸坯弯曲矫直过程中,铸坯温度较低的角部区域易发生裂纹,基于上述因素,此阶段的已凝固区域权值较小,未凝固区域适用较大权值,结合现场实验,已凝固区域的域加权值取值为0≤γ1≤0.5,未凝固区域的域加权值取值为0.5≤γ2≤1.0。
(2)当前铸流位置介于铸坯压下起始铸流位置和铸坯凝固终点铸流位置之间即Pstart<P<Psolid时,已凝固区域的域加权值取值为0.3≤γ1≤0.8,未凝固区域的域加权值取值为0.2≤γ2≤0.7。
本实施方式中,将当前铸流位置介于铸坯压下起始铸流位置和铸坯凝固终点铸流位置之间即Pstart<P<Psolid定义为第二阶段,此阶段通过实施轻压下,对铸坯中心未凝固区域通过挤压作用改善铸坯的中心疏松、偏析等缺陷,如果基础辊缝制定的偏大,则势必削弱压下作用效果,因此,在此阶段需要适当加大收缩值比较大的已凝固区域的权值,结合现场实验,已凝固区域的域加权值取值为0.3≤γ1≤0.8,未凝固区域的域加权值取值为0.2≤γ2≤0.7。
(3)当前铸流位置大于铸坯凝固终点铸流位置即P>Psolid时,已凝固区域的域加权值取值为γ1=1.0,未凝固区域消失,则γ2=0。
本实施方式中,将当前铸流位置大于铸坯凝固终点铸流位置即P>Psolid定义为第三阶段,该阶段铸坯已完全凝固,铸坯未凝固区域消失,已凝固区域的域加权值取值为γ1=1.0,未凝固区域消失,则γ2=0。
凝固传热三维有限元模型中已凝固区域和未凝固区域的划分示意图如图8所示,铸坯表面的“已凝固区域”、“未凝固区域”是根据铸坯的中心部位是否已完全凝固进行划分的。其中已凝固区域对应的铸坯中心完全凝固区域如图8中b所示,未凝固区域对应的铸坯中心未完全凝固区域如图8中a所示。
浇铸过程中三个阶段已凝固区域的域加权值和未凝固区域的域加权值如表3所示:
表3浇铸过程中三个阶段已凝固区域的域加权值和未凝固区域的域加权值
浇铸阶段 | 第一阶段 | 第二阶段 | 第三阶段 |
已凝固区域权值 | 0.25 | 0.45 | 1.00 |
未凝固区域权值 | 0.75 | 0.55 | 0.00 |
步骤4.3:根据浇铸过程中的三维热收缩规律、板坯的初始厚度、以及当前铸流所处浇铸过程中的不同位置的域加权值,制定相应工况条件下的连铸机基础辊缝。
确定板坯连铸机各扇形段基础辊缝值的计算公式如(10)所示:
其中,Gap为基础辊缝值,InGap为铸坯初始厚度,Sb为已凝固区域起始节点,Se为已凝固区域结束节点,Mb为未凝固区域起始节点,Me为未凝固区域结束节点,i为凝固传热三维有限元模型中已凝固区域的节点,j为凝固传热三维有限元模型中未凝固区域的节点,Si为已凝固区域内的节点i的厚度方向收缩量,Sj为未凝固区域内的节点j的厚度方向收缩量,为已凝固区域内节点i的热收缩加权值,N1为凝固传热三维有限元模型中已凝固区域的节点个数,为未凝固区域内的节点j的热收缩加权值,N2为凝固传热三维有限元模型中未凝固区域的节点个数,γ1为已凝固区域的域加权值,γ2为未凝固区域的域加权值。
本实施方式中,确定的板坯连铸机各扇形段基础辊缝值如表4所示
表4板坯连铸机各扇形段基础辊缝值
扇形段 | 入口位置,m | 出口位置,m | 入口辊缝,mm | 出口辊缝,mm |
0段 | 1.040 | 4.185 | 261.70 | 260.15 |
1段 | 4.446 | 6.078 | 260.07 | 259.72 |
2段 | 6.368 | 8.000 | 259.67 | 259.35 |
3段 | 8.290 | 9.922 | 259.30 | 258.94 |
4段 | 10.212 | 11.844 | 258.89 | 258.60 |
5段 | 12.134 | 13.766 | 258.55 | 258.26 |
6段 | 14.056 | 15.688 | 258.21 | 257.93 |
7段 | 16.028 | 18.044 | 257.88 | 257.50 |
8段 | 18.409 | 20.445 | 257.43 | 256.90 |
9段 | 20.785 | 22.825 | 256.75 | 256.08 |
10段 | 23.165 | 25.205 | 255.98 | 255.42 |
11段 | 25.545 | 27.585 | 255.33 | 254.82 |
12段 | 27.925 | 29.965 | 254.74 | 254.29 |
13段 | 30.305 | 32.345 | 254.21 | 253.79 |
14段 | 32.685 | 34.725 | 253.72 | 253.32 |
本发明对现场铸坯质量改善作用明显,在使用本发明方法之前,本实例所涉及的铸坯中心偏析≤1.0比例为62.56%,中心疏松≤1.0比例为89.62%;实施本发明方法之后,经过连续4个月的生产统计,铸坯中心偏析≤1.0比例为89.67%,中心疏松≤1.0比例为100%。
Claims (7)
1.一种宽厚板坯连铸机基础辊缝制定方法,其特征在于,包括以下步骤:
步骤1:采集板坯连铸机的浇铸工况参数;
步骤2:建立板坯的凝固传热三维有限元模型,根据连铸机的浇铸工况参数计算得到浇铸过程中的三维温度场;
步骤3:建立板坯的三维热收缩有限元模型,以浇铸过程中的三维温度场作为温度载荷求解该三维热收缩有限元模型,得到浇铸过程中的三维热收缩规律;
步骤4:根据获取的浇铸过程中的三维热收缩规律及板坯的初始厚度制定相应工况条件下的连铸机基础辊缝;
步骤4.1:根据浇铸过程中的三维温度场和实施轻压下时压下区间的起始固相率确定铸坯压下起始铸流位置和铸坯凝固终点铸流位置;
步骤4.2:根据铸坯压下起始铸流位置和铸坯凝固终点铸流位置确定当前铸流位置的域加权值;
步骤4.3:根据浇铸过程中的三维热收缩规律、板坯的初始厚度、以及当前铸流所处浇铸过程中的不同位置的域加权值,制定相应工况条件下的连铸机基础辊缝。
2.根据权利要求1所述的宽厚板坯连铸机基础辊缝制定方法,其特征在于,所述浇铸工况参数包括:浇铸钢种的断面尺寸、浇铸钢种成分、浇铸温度、工作拉速、结晶器水流量及回水温差、二冷区划分参数及各二冷区内的水量。
3.根据权利要求1所述的宽厚板坯连铸机基础辊缝制定方法,其特征在于,所述的步骤2包括以下步骤:
步骤2.1:取板坯浇铸过程中一定高度的1/4铸坯,建立其凝固传热三维有限元模型;
步骤2.2:根据连铸机的浇铸工况参数确定凝固传热三维有限元模型的不同铸流位置处的边界条件,求解凝固传热三维有限元模型,得到浇铸过程中的三维温度场;
步骤2.3:通过实测板坯的表面温度及坯壳厚度,对凝固传热三维有限元模型进行修正,从而对板坯浇铸过程中的三维温度场进行修正。
4.根据权利要求1所述的宽厚板坯连铸机基础辊缝制定方法,其特征在于,所述的步骤3包括以下步骤:
步骤3.1:取板坯浇铸过程中一定高度的1/4铸坯,建立其三维热收缩有限元模型;
步骤3.2:将浇铸过程中的三维温度场作为温度载荷,通过刚体对称面完成对该有限元模型宽度方向和厚度方向的位移约束,求解该三维热收缩有限元模型,得到浇铸过程中的三维热收缩规律。
5.根据权利要求1所述的宽厚板坯连铸机基础辊缝制定方法,其特征在于,所述的制定相应工况条件下的连铸机基础辊缝的计算公式如下:
其中,Gap为基础辊缝值,InGap为铸坯初始厚度,Sb为已凝固区域起始节点,Se为已凝固区域结束节点,Mb为未凝固区域起始节点,Me为未凝固区域结束节点,i为凝固传热三维有限元模型中已凝固区域的节点,j为凝固传热三维有限元模型中未凝固区域的节点,Si为已凝固区域内的节点i的厚度方向收缩量,Sj为未凝固区域内的节点j的厚度方向收缩量,为已凝固区域内节点i的热收缩加权值,N1为凝固传热三维有限元模型中已凝固区域的节点个数,为未凝固区域内的节点j的热收缩加权值,N2为凝固传热三维有限元模型中未凝固区域的节点个数,γ1为已凝固区域的域加权值,γ2为未凝固区域的域加权值。
6.根据权利要求1所述的宽厚板坯连铸机基础辊缝制定方法,其特征在于,所述的根据铸坯压下起始铸流位置和铸坯凝固终点铸流位置确定当前铸流位置的域加权值具体方法为:
(1)当前铸流位置小于铸坯压下起始铸流位置时,已凝固区域的域加权值取值为0≤γ1≤0.5,未凝固区域的域加权值取值为0.5≤γ2≤1.0;
(2)当前铸流位置介于铸坯压下起始铸流位置和铸坯凝固终点铸流位置之间时,已凝固区域的域加权值取值为0.3≤γ1≤0.8,未凝固区域的域加权值取值为0.2≤γ2≤0.7;
(3)当前铸流位置大于铸坯凝固终点铸流位置时,已凝固区域的域加权值取值为γ1=1.0,未凝固区域消失,则γ2=0。
7.根据权利要求3所述的宽厚板坯连铸机基础辊缝制定方法,其特征在于,所述的凝固传热三维有限元模型的不同铸流位置处的边界条件分别为:以铸坯在结晶器中的热流密度作为结晶器的边界条件,以实测水流密度确定的等效换热系数作为二冷区边界条件,以空冷区内的辐射散热作为空冷区的边界条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510541805.1A CN105033214B (zh) | 2015-08-28 | 2015-08-28 | 一种宽厚板坯连铸机基础辊缝制定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510541805.1A CN105033214B (zh) | 2015-08-28 | 2015-08-28 | 一种宽厚板坯连铸机基础辊缝制定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105033214A true CN105033214A (zh) | 2015-11-11 |
CN105033214B CN105033214B (zh) | 2017-03-22 |
Family
ID=54440510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510541805.1A Active CN105033214B (zh) | 2015-08-28 | 2015-08-28 | 一种宽厚板坯连铸机基础辊缝制定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105033214B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105798253A (zh) * | 2016-06-06 | 2016-07-27 | 重庆大学 | 一种异型坯连铸二冷辊缝调整方法 |
CN106001478A (zh) * | 2016-07-08 | 2016-10-12 | 首钢总公司 | 一种板坯连铸机基础辊缝工艺的制定方法 |
CN108052757A (zh) * | 2017-12-22 | 2018-05-18 | 江苏理工学院 | 一种确定板带钢热轧过程工作辊温度场的有限元方法 |
CN109492317A (zh) * | 2018-11-20 | 2019-03-19 | 中冶赛迪工程技术股份有限公司 | 基于连铸机二维温度场仿真方法及监控模型的运行方法 |
CN109783994A (zh) * | 2019-03-08 | 2019-05-21 | 东北大学 | 一种基于坯壳鼓肚和机械压下的厚板坯宏观偏析计算方法 |
CN110568010A (zh) * | 2019-09-05 | 2019-12-13 | 中冶赛迪工程技术股份有限公司 | 一种板坯连铸内部裂纹在线预测及定位的方法 |
CN111299534A (zh) * | 2020-03-23 | 2020-06-19 | 中冶南方连铸技术工程有限责任公司 | 连铸轻重压下基准辊缝确定方法及装置 |
CN111914360A (zh) * | 2020-07-10 | 2020-11-10 | 中国三冶集团有限公司 | 一种基于有限元分析连铸机扇形段导辊在线周期的方法 |
CN112475251A (zh) * | 2019-09-12 | 2021-03-12 | 上海梅山钢铁股份有限公司 | 一种板坯连铸动态轻压下压下区间的确定方法 |
CN113000804A (zh) * | 2021-02-24 | 2021-06-22 | 中冶赛迪工程技术股份有限公司 | 一种板坯连铸动态轻压下压下区间最佳起始点和结束点的选择方法 |
CN114239354A (zh) * | 2021-12-14 | 2022-03-25 | 中冶赛迪重庆信息技术有限公司 | 连铸工艺过程中铸坯鼓肚量的预测方法 |
CN114833320A (zh) * | 2022-05-19 | 2022-08-02 | 中冶南方工程技术有限公司 | 一种支持多坯型的连铸三维温度场可视化系统 |
CN115401178A (zh) * | 2021-05-28 | 2022-11-29 | 宝山钢铁股份有限公司 | 一种改善齿轮钢内部质量的压下工艺确定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04231158A (ja) * | 1990-12-28 | 1992-08-20 | Kobe Steel Ltd | 連続鋳造機のシミュレーション装置 |
CN101362196A (zh) * | 2008-09-12 | 2009-02-11 | 攀钢集团研究院有限公司 | 一种基于大方坯自然热收缩辊缝工艺的控制方法 |
CN101362195A (zh) * | 2008-09-12 | 2009-02-11 | 攀钢集团研究院有限公司 | 一种大方坯连铸动态轻压下压下量的控制方法 |
CN101554650A (zh) * | 2009-05-22 | 2009-10-14 | 东北大学 | 一种大方坯连铸机轻压下拉矫辊辊缝在线标定的控制方法 |
CN101704079A (zh) * | 2009-09-28 | 2010-05-12 | 田陆 | 用于连铸坯浇铸的控制方法 |
-
2015
- 2015-08-28 CN CN201510541805.1A patent/CN105033214B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04231158A (ja) * | 1990-12-28 | 1992-08-20 | Kobe Steel Ltd | 連続鋳造機のシミュレーション装置 |
CN101362196A (zh) * | 2008-09-12 | 2009-02-11 | 攀钢集团研究院有限公司 | 一种基于大方坯自然热收缩辊缝工艺的控制方法 |
CN101362195A (zh) * | 2008-09-12 | 2009-02-11 | 攀钢集团研究院有限公司 | 一种大方坯连铸动态轻压下压下量的控制方法 |
CN101554650A (zh) * | 2009-05-22 | 2009-10-14 | 东北大学 | 一种大方坯连铸机轻压下拉矫辊辊缝在线标定的控制方法 |
CN101704079A (zh) * | 2009-09-28 | 2010-05-12 | 田陆 | 用于连铸坯浇铸的控制方法 |
Non-Patent Citations (1)
Title |
---|
陈洪智等: "不锈钢板坯连铸自由线收缩与辊缝研究", 《中国冶金》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105798253B (zh) * | 2016-06-06 | 2017-08-04 | 重庆大学 | 一种异型坯连铸二冷辊缝调整方法 |
CN105798253A (zh) * | 2016-06-06 | 2016-07-27 | 重庆大学 | 一种异型坯连铸二冷辊缝调整方法 |
CN106001478A (zh) * | 2016-07-08 | 2016-10-12 | 首钢总公司 | 一种板坯连铸机基础辊缝工艺的制定方法 |
CN106001478B (zh) * | 2016-07-08 | 2018-01-23 | 首钢总公司 | 一种板坯连铸机基础辊缝工艺的制定方法 |
CN108052757A (zh) * | 2017-12-22 | 2018-05-18 | 江苏理工学院 | 一种确定板带钢热轧过程工作辊温度场的有限元方法 |
CN109492317A (zh) * | 2018-11-20 | 2019-03-19 | 中冶赛迪工程技术股份有限公司 | 基于连铸机二维温度场仿真方法及监控模型的运行方法 |
CN109492317B (zh) * | 2018-11-20 | 2023-06-20 | 中冶赛迪工程技术股份有限公司 | 基于连铸机二维温度场仿真方法及监控模型的运行方法 |
CN109783994A (zh) * | 2019-03-08 | 2019-05-21 | 东北大学 | 一种基于坯壳鼓肚和机械压下的厚板坯宏观偏析计算方法 |
CN109783994B (zh) * | 2019-03-08 | 2023-08-25 | 东北大学 | 一种基于坯壳鼓肚和机械压下的厚板坯宏观偏析计算方法 |
CN110568010B (zh) * | 2019-09-05 | 2022-03-11 | 中冶赛迪工程技术股份有限公司 | 一种板坯连铸内部裂纹在线预测及定位的方法 |
CN110568010A (zh) * | 2019-09-05 | 2019-12-13 | 中冶赛迪工程技术股份有限公司 | 一种板坯连铸内部裂纹在线预测及定位的方法 |
CN112475251A (zh) * | 2019-09-12 | 2021-03-12 | 上海梅山钢铁股份有限公司 | 一种板坯连铸动态轻压下压下区间的确定方法 |
CN112475251B (zh) * | 2019-09-12 | 2022-03-15 | 上海梅山钢铁股份有限公司 | 一种板坯连铸动态轻压下压下区间的确定方法 |
CN111299534B (zh) * | 2020-03-23 | 2021-12-17 | 中冶南方连铸技术工程有限责任公司 | 连铸轻重压下基准辊缝确定方法及装置 |
CN111299534A (zh) * | 2020-03-23 | 2020-06-19 | 中冶南方连铸技术工程有限责任公司 | 连铸轻重压下基准辊缝确定方法及装置 |
CN111914360A (zh) * | 2020-07-10 | 2020-11-10 | 中国三冶集团有限公司 | 一种基于有限元分析连铸机扇形段导辊在线周期的方法 |
CN113000804A (zh) * | 2021-02-24 | 2021-06-22 | 中冶赛迪工程技术股份有限公司 | 一种板坯连铸动态轻压下压下区间最佳起始点和结束点的选择方法 |
CN115401178A (zh) * | 2021-05-28 | 2022-11-29 | 宝山钢铁股份有限公司 | 一种改善齿轮钢内部质量的压下工艺确定方法 |
CN115401178B (zh) * | 2021-05-28 | 2023-07-07 | 宝山钢铁股份有限公司 | 一种改善齿轮钢内部质量的压下工艺确定方法 |
CN114239354A (zh) * | 2021-12-14 | 2022-03-25 | 中冶赛迪重庆信息技术有限公司 | 连铸工艺过程中铸坯鼓肚量的预测方法 |
CN114239354B (zh) * | 2021-12-14 | 2024-06-04 | 中冶赛迪信息技术(重庆)有限公司 | 连铸工艺过程中铸坯鼓肚量的预测方法 |
CN114833320A (zh) * | 2022-05-19 | 2022-08-02 | 中冶南方工程技术有限公司 | 一种支持多坯型的连铸三维温度场可视化系统 |
CN114833320B (zh) * | 2022-05-19 | 2023-10-27 | 中冶南方工程技术有限公司 | 一种支持多坯型的连铸三维温度场可视化系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105033214B (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105033214B (zh) | 一种宽厚板坯连铸机基础辊缝制定方法 | |
CN104493121B (zh) | 一种大方坯连铸生产过程的凝固末端位置在线检测方法 | |
CN112475251B (zh) | 一种板坯连铸动态轻压下压下区间的确定方法 | |
CN103920859B (zh) | 一种连铸板坯内部裂纹在线预测方法 | |
CN103028615B (zh) | 一种预测带钢热连轧过程温度演变的方法 | |
CN103273272B (zh) | 一种综合提升宽厚板坯内外质量的锻、轧复合成形方法 | |
CN112570675B (zh) | 宽厚板连铸板坯轻压下过程最小理论压下量的确定方法 | |
CN103302104A (zh) | 热轧硅钢的制造方法 | |
CN105798253B (zh) | 一种异型坯连铸二冷辊缝调整方法 | |
CN109013717B (zh) | 一种热连轧中间坯心部温度计算方法 | |
CN107243611B (zh) | 一种大方坯连铸凝固末端单辊压下位置确定方法 | |
CN105930594B (zh) | 一种立轧轧件狗骨形状预测方法 | |
CN104874754A (zh) | 一种板坯窄面内凸型曲面结晶器及其设计方法 | |
CN104331540A (zh) | 一种连铸二冷区冷却水对流换热系数的优化方法 | |
Zhang et al. | Evolution of transverse weld during porthole extrusion of AA7N01 hollow profile | |
CN101362196A (zh) | 一种基于大方坯自然热收缩辊缝工艺的控制方法 | |
CN102083573B (zh) | 金属连铸坯的连铸方法 | |
WO2020181575A1 (zh) | 一种大方坯凸形辊连铸重压下偏析行为预测方法 | |
CN104399917B (zh) | 一种具有增强水冷结构的渐进弧型倒角结晶器及设计方法 | |
CN116738518B (zh) | 一种连铸轻压下裂纹位置的数值模拟验证方法及内部质量控制方法 | |
Kim et al. | Continuous Casting and Rolling for Aluminum Alloy Wire and Rod | |
CN116776668B (zh) | 一种小方坯连铸结晶器内坯壳凝固收缩计算方法 | |
CN106001478B (zh) | 一种板坯连铸机基础辊缝工艺的制定方法 | |
CN115344957A (zh) | 一种渐变曲率曲型辊的设计方法 | |
CN107096898B (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 |