发明名称 圆坯连铸结晶器电磁搅拌参数的制定方法
摘要 本发明涉及圆坯连铸结晶器电磁搅拌参数的制定方法。为克服目前对搅拌参数的确定研究较少或研究结论相差较大的问题,发明一种结合钢种、铸坯尺寸、铸机和电磁搅拌设备,通过数值仿真分析和试验验证,确定电磁搅拌的电流强度和频率等参数的制定方法。该方法是在电磁搅拌磁场分布模拟结果的基础之上,对电磁搅拌过程结晶器内钢液的流场、温度场及夹杂物运动轨迹进行耦合计算,针对钢液的流动、凝固传热及夹杂物轨迹等现象做考察判断电磁搅拌参数的合理性,并再将计算结果与现场取样试验进行反馈对比分析,以确定电磁搅拌参数。用本发明的方法,可快速、准确有效地确定不同条件下电磁搅拌的电流强度和频率,明显减少铸坯表面和皮下夹杂物及气泡,扩大等轴晶率,改善芯部质量,提高组织致密性和均匀性。
申请公布号 CN100566888C 申请公布日期 2009.12.09
申请号 CN200710060348.X 申请日期 2007.12.19
申请人 天津钢铁有限公司;东北大学 发明人 曾小平;许克亮;朱苗勇;凌遵峰;吴波;郑淑国
分类号 B22D11/115(2006.01)I 主分类号 B22D11/115(2006.01)I
代理机构 天津才智专利商标代理有限公司 代理人 吕志英
主权项 1、圆坯连铸结晶器电磁搅拌参数的制定方法,其特征在于:(1)采用的磁场偏微分方程<maths num="0001"><![CDATA[<math><mrow><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>A</mi><mo>-</mo><mi>&mu;&epsiv;</mi><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>A</mi></mrow><mrow><mo>&PartialD;</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>=</mo><mo>-</mo><mi>&mu;j</mi><mo>,</mo></mrow></math>]]></maths>和电场偏微分方程<img file="C2007100603480002C2.GIF" wi="432" he="143" />确定电磁搅拌电、磁场分布,式中,A为矢量磁势;j为电流密度;μ为磁导率;ρ为电荷密度;ε为介质电容率;<img file="C2007100603480002C3.GIF" wi="30" he="37" />为标量电势;t为时间;(2)结晶器电磁搅拌下,铸坯内流场和温度场的控制方程包括:连续方程:<maths num="0002"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;</mi><msub><mi>u</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>动量方程:<maths num="0003"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;</mi><msub><mi>u</mi><mi>j</mi></msub><msub><mi>u</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><mo>-</mo><mfrac><mrow><mo>&PartialD;</mo><mi>p</mi></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>i</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>[</mo><mrow><mo>(</mo><mi>&mu;</mi><mo>+</mo><msub><mi>&mu;</mi><mi>t</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><mrow><mo>&PartialD;</mo><msub><mi>u</mi><mi>i</mi></msub></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>&PartialD;</mo><msub><mi>u</mi><mi>j</mi></msub></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>i</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msub><mi>F</mi><mi>i</mi></msub><mo>,</mo></mrow></math>]]></maths>能量方程:<maths num="0004"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;</mi><msub><mi>u</mi><mi>j</mi></msub><mi>T</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>[</mo><mrow><mo>(</mo><mfrac><mi>&mu;</mi><mi>Pr</mi></mfrac><mo>+</mo><mfrac><msub><mi>&mu;</mi><mi>t</mi></msub><msub><mi>&sigma;</mi><mi>T</mi></msub></mfrac><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><mi>T</mi></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>]</mo><mo>-</mo><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;</mi><msub><mi>u</mi><mi>j</mi></msub><mi>&Delta;H</mi><mo>)</mo></mrow></mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>j</mi></msub></mfrac><mo>,</mo></mrow></math>]]></maths>以及k方程:<maths num="0005"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;</mi><msub><mi>u</mi><mi>j</mi></msub><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>[</mo><mrow><mo>(</mo><mi>&mu;</mi><mo>+</mo><mfrac><msub><mi>&mu;</mi><mi>t</mi></msub><msub><mi>&sigma;</mi><mi>k</mi></msub></mfrac><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><mi>k</mi></mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>j</mi></msub></mfrac><mo>]</mo><mo>+</mo><mi>G</mi><mo>-</mo><mi>&rho;&epsiv;</mi><mo>-</mo><msub><mi>D</mi><mi>k</mi></msub><mo>,</mo></mrow></math>]]></maths>ε方程:<maths num="0006"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;</mi><msub><mi>u</mi><mi>j</mi></msub><mi>&epsiv;</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>[</mo><mrow><mo>(</mo><mi>&mu;</mi><mo>+</mo><mfrac><msub><mi>&mu;</mi><mi>t</mi></msub><msub><mi>&sigma;</mi><mi>&epsiv;</mi></msub></mfrac><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><mi>&epsiv;</mi></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>]</mo><mo>+</mo><msub><mi>f</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub><mi>G</mi><mfrac><mi>&epsiv;</mi><mi>k</mi></mfrac><mo>-</mo><msub><mi>f</mi><mn>2</mn></msub><msub><mi>C</mi><mn>2</mn></msub><mi>&rho;</mi><mfrac><msup><mi>&epsiv;</mi><mn>2</mn></msup><mi>k</mi></mfrac><mo>+</mo><msub><mi>E</mi><mi>&epsiv;</mi></msub></mrow></math>]]></maths>以上各式中,ρ为钢液密度;u<sub>i</sub>和u<sub>j</sub>(i,j=1,2,3)为不同坐标方向下的速度;F<sub>i</sub>为电磁力;μ为动力黏度;p为压力;μ<sub>t</sub>为湍流黏度;k为湍流动能;ε为湍流耗散率;T为温度;Pr为普朗特数;C<sub>p</sub>为等压比热容;λ为热导率;σ<sub>T</sub>=1.0,C<sub>1</sub>=1.44,C<sub>2</sub>=1.92,σ<sub>k</sub>=1.0,σ<sub>ε</sub>=1.3,f<sub>1</sub>=1.0;<maths num="0007"><![CDATA[<math><mrow><mi>G</mi><mo>=</mo><msub><mi>&mu;</mi><mi>t</mi></msub><mrow><mo>(</mo><mfrac><mrow><mo>&PartialD;</mo><msub><mi>u</mi><mi>i</mi></msub></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>+</mo><mfrac><msub><mrow><mo>&PartialD;</mo><mi>u</mi></mrow><mi>j</mi></msub><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>i</mi></msub></mrow></mfrac><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><msub><mi>u</mi><mi>i</mi></msub></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>;</mo></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>D</mi><mi>k</mi></msub><mo>=</mo><mn>2</mn><mi>&mu;</mi><msup><mrow><mo>(</mo><mfrac><mrow><mo>&PartialD;</mo><msqrt><mi>k</mi></msqrt></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>i</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msub><mi>E</mi><mi>&epsiv;</mi></msub><mo>=</mo><mn>2</mn><mfrac><mrow><msub><mi>&mu;</mi><mi>t</mi></msub><mi>&mu;</mi></mrow><mi>&rho;</mi></mfrac><msup><mrow><mo>(</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><msub><mi>u</mi><mi>i</mi></msub></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>j</mi></msub><mo>&PartialD;</mo><msub><mi>x</mi><mi>k</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths>f<sub>2</sub>=1-0.3exp(-Re<sub>t</sub><sup>2</sup>),流场计算采用有限体积法离散流场微分方程,用交错网格存储速度分量,采用Simplec算法,当连续方程的质量源小于10<sup>-4</sup>时,即认为迭代收敛,采用Fortran语言编程,将ANSYS计算的电磁力加入动量源项,结晶器计算长度为1.8m;(3)在上述流场计算结果的基础上,加入夹杂物运动轨迹离散模型控制方程<maths num="0010"><![CDATA[<math><mrow><mi>m</mi><mfrac><msub><mi>dV</mi><mi>p</mi></msub><mi>dt</mi></mfrac><mo>=</mo><msub><mi>F</mi><mi>P</mi></msub><mo>+</mo><msub><mi>F</mi><mi>B</mi></msub><mo>+</mo><msub><mi>F</mi><mi>A</mi></msub><mo>+</mo><msub><mi>F</mi><mi>D</mi></msub><mo>,</mo></mrow></math>]]></maths>式中,m为粒子质量;V<sub>p</sub>为粒子速度;t表示时间;F<sub>p</sub>为钢液压力梯度引起的力;F<sub>B</sub>为钢液对粒子的浮力;F<sub>A</sub>为钢液附加质量力;F<sub>D</sub>为粒子所受粘性阻力,设置夹杂物粒度5~50μm,密度3500kg/m<sup>3</sup>,质量流率4×10<sup>-5</sup>kg/s参数及离散边界条件,继续对电磁搅拌条件下的流场、温度场及夹杂物轨迹进行耦合计算;(4)根据模拟计算出的圆坯结晶器电磁搅拌的磁场、流场、温度场以及夹杂物运动轨迹的数据,并针对结晶器内有效搅拌区的钢水流动速度在0.2-0.25m/s之间,可产生足够大的离心力和剪切力;可使保护渣熔融良好;可造成热顶效应;可使夹杂物在结晶器内停留的时间足够长,使其尽量上浮做考察判断,以修正模拟计算出的数据;(5)在采用经过模拟计算、考察、判断及修正后所得到的励磁电流强度和频率条件下,检测铸坯的皮下裂纹、芯部裂纹、缩孔、中心疏松、皮下气泡及等轴晶率,最后确定电磁搅拌参数。
地址 300301天津市东丽区津塘公路398号