1045钢淬火时非线性表面换热系数和连续冷却曲线的数值模拟材料学院沙军.ppt

1045钢淬火时非线性表面换热系数和连续冷却曲线的数值模拟 材料学院 沙军威 2011208100 李 响 2011208083 陈利霞 2011208039 Contents 钢淬火时连续冷却的数学模型 4 背景介绍 1 实验过程 2 非线性导热方程和有限差分格式 3 表面换热系数和冷却曲线的计算及结论 5 背景介绍 热处理是改善金属材料的机械性能和工件加工性能的重要手段。热处理过程会极大影响钢铁材料内部的显微组织，而组织的不同会导致钢铁材料具有不同的力学性能。 常用的热处理工艺 退 火 淬 火 回 火 正 火 在淬火的过程中，控制淬火工件的热应力和热应变是比较困难的。如果不能控制好冷却速度，就会造成工件的变形和开裂，导致工件报废。淬火时温度场的计算成为了研究的重点。目前，数值模拟方法被广泛应用于温度场的计算，而表面换热系数和冷却曲线的数值模拟又是温度场数值模拟的基础。 实验过程 本文采用圆柱体1045钢作为工件，应用热电偶测量钢淬火时的冷却曲线，利用显式有限差分和非线性估计法求解淬火时非线性导热方程的逆问题，给出一种求解钢淬火时表面换热系数的方法。基于1045钢淬火时的连续冷却曲线和相成分体积含量的变化，求解了1045钢淬火时的非线性表面换热系数。 工件为Φ45mm*120mm的1045钢，垂直于圆柱体端面钻孔，并将3支直径为0.02mm镍铬-镍铝热电偶埋入距端面60mm的位置。热电偶1位于距圆柱表面1mm的位置，热电偶2位于距圆柱体轴线11.25mm的位置，热电偶3位于圆柱体的轴线上。将试件加热至860℃开始淬火，淬入20℃的静水中，淬火时试件温度的变化由热电偶测量，并通过计算机自动采样系统实时采样记录，采样频率为12次/s。热电偶2和3所测量的温度主要用于检验计算结果。而热电偶1所测量的温度值是求解非线性热传导逆问题的关键参数。 非线性导热方程和有限差分格式 淬火时试件中的导热问题是典型的轴对称问题，且试件在淬火之前温度均匀一致，因此试件中部的热传导方程为 式中， 为钢的密度，cp是定压比热，θ为温度，r为圆柱体的半径，t为时间， 为导热系数，与温度和相成分有关。因此，式（1）为非线性导热方程。 (1) 非线性导热方程和有限差分格式 非线性导热方程的有限差分格式为： (2) (3) (4) (5) (6) 其中： 非线性导热方程和有限差分格式 将860℃的试件放入20℃的水中时，淬火换热边界条件为 对应的差分格式为： 是淬火介质的温度， 是试件表面温度，h是表面换热系数，该系数是表面温度和相变体积量的函数。 钢淬火时连续冷却的数学模型 基于1045钢的等温转变曲线，应用数学转化可以模拟1045钢淬火时的连续转变情况。铁素体/珠光体和贝氏体的等温转变动力学可以应用如下的JMA定律： (9) 系数Cjk和Djk可以由1045钢的TTT曲线和最小二乘法确定。 b和n是相变时的材料参数。这里，k=1定义为铁素体/珠光体，k=2定义为贝氏体。 材料参数b和n可以由1045钢等温转变曲线和具体的实验数据确定。参数b和n可以由多项式来拟合得到： 钢淬火时连续冷却的数学模型 为了应用式（9）来计算钢淬火过程连续冷却时的相成分，可采用Ivan Tzitzelkov数学变换法，计算时采用虚拟时间步长τ*作为迭代计算时的时间步，该时间步长为： (12) 马氏体相变体积分数可由下式来计算： tm+1为m+1时刻， 为m时刻珠光体或贝氏体的体积分数。 将式（12）代入式（9），可得： 表面换热系数和冷却曲线的计算 1、应用非线性估计法，由淬火时换热边界条件可得钢淬火时表面换热系数的增量： (12) 表面换热系数的增量△h已经导出，其中 为非线性系数。在表面换热系数的计算中就可以用表面换热系数的增量△h来进行迭代计算。 表面换热系数和冷却曲线的计算 2、计算结果 下图为计算得到的1045钢淬火时表面换热系数与试件表面温度的关系。 下图为应用有限差分法计算得到的1045钢试件淬火时的冷却曲线。 淬火时试件中截面上距表面1.5mm处的相变情况 淬火后试件中各主要部分相成分的计算值和实测值: 结论 钢淬火时的导热、换热情况是非线性的。在钢淬火时表面换热系数测试技术和换热机理不太完善的情况下，应用有限差分法、非线性估计法和实测的冷却曲线来确定钢淬火时表面换热系数和温度场等非线性热问题是非常简便有效的。它不仅利用了温度测试技术较简单成熟的优点，而且适用于各种淬火液，计算速度快，因此便于工程应用。 谢谢大家！