考虑相变过程高碳钢线材温度场计算.docVIP

考虑相变过程的高碳钢线材温度场计算 蒲玉梅　蔡庆伍　米振莉　唐荻 摘　要：根据高碳钢线材的冷却条件，考虑了相变过程的相变潜热，计算出线材搭接点与非搭接点在斯太尔摩冷却线上的温度分布，计算结果与实测结果吻合较好。关键词：线材；控制冷却；温度场；计算 Simulation of temperature field considering phase transformation of high carbon steel wire PU Yu-mei1,CAI Qing-wu2,MI Zhen-li2,TANG Di2(1.Maanshan Iron Steel Co.,Ltd.,Maanshan 243000,China;2.University of Science Technology Beijing) Abstract:Based on cooling condition,considering potential heat of phase transformation,the temperature field of lap joint and non-lap joint of high carbon steel wire on STELMOR cooling line are calculated.Thecalcu～lated results are good agreement with practical test.Key words:wire;control cooling;temperature field;calculation 　　目前，控制冷却逐步应用于钢材生产中，通过采用不同的冷却制度来控制成品的组织和性能。　　作为预应力钢丝、钢铰线原料的大规格高碳钢线材，要求严格的组织分布，工艺上主要通过选定合理的吐丝温度及在其后的斯太尔摩冷却曲线实现。在斯太尔摩冷却线上，由于线材是成圈平辅的，造成搭接点与非搭接点处冷却不均匀。因此，准确计算高碳钢线材在斯太尔摩冷却线上的温度场，可为控制线材均匀冷却提供依据。温度场的数值模拟广泛采用有限元法，同时，ANSYS热分析使有限元法具有准确处理稳态、瞬态温度场的功能。 1　数学模型建立及初始、边界条件　　工件冷却过程中，其热传导符合Fourier导热方程：
（1） 式中，q为内热源的热流密度/W.m-2;T为温度/℃；t为过程进行的时间/s;ρ为材料密度/kg.m-3;Cp为材料定压比热/J.(kg.℃)-1;k为导热系数/W.(m.k)-1。　　初始条件：t=0,0≤r≤R0,T=Tin　　边界条件：线材中心：t0　R=0　-k=0　　线材表面：式中：Tin为吐丝温度/℃；TR0为线材表面温度/℃；Ta为环境温度/℃；R0为线材半径/m;h为热交换系数/W.m-2.K-1;r为线材中心到表面的任意位置．　　依照文献［1～2］介绍的方法，分别选取搭接点与非搭接点的热交换系数作为温度场计算的边界条件．　　依据文献［3］，可确定（1）式中的热物性参数． 2　组织转变潜热的考虑　　高碳钢由奥氏体向低温冷却时发生珠光体转变，此为放热效应，因此应考虑其对温度场的影响，即用有内热源的导热微分方程计算这种过程．在冷却过程中，计算出组织转变量，按转变量可计算内热源的热流密度．2.1　转变过程的确定　　对于连续冷却条件下的转变过程有很多模型表示，在此选取转变量与温度成线性关系的模型： V=Vo(Ts-T)/Ts-Tf　(2) 式中：V为组织转变量；Vo为该组织在该冷却速度下的最大转变量；Ts、Tf分别为转变开始、终了温度，T为连续冷却所达到的温度。2.2　平均冷却速度的处理　　按文献［4］介绍的方法，将临界温度作为起始点，计算平均冷却速度。当温度进入转变区后，用内插法计算各冷却速度对应的转变开始点和转变终了点，按（2）式计算转变量。2.3　模拟计算步骤　　（1）确定温度达到临界点Ac1的时间t0,作为冷却开始时间；　　（2）随着温度T降低，不断计算冷却速度CR,CR=(Ac1-T)/(t-t0)　　（3）根据CCT曲线，求出该冷却速度下的开始转变温度Tt;　　（4）如果T=Tt,且为扩散型转变，记下达到转变温度的时间t1;　　（5）随温度T降低，不断计算冷速CR,CR=(Tt-T)/(t-t1);　　（6）根据CCT曲线，求出该冷却速度下的转变终了温度和最大转变量；　　（7）如果T=Tt为马氏体转变，则按V=1-exp（-α(Ms-T）），根据当时温度计算马氏体转变量；　　（8）计算直到温度降到规定值或转变量达到规定值为止。 3　斯太尔摩冷却模拟结果3.1　搭接与非搭接点的温度场模拟　　对12.5mm线材搭接点与非搭接点在斯太尔