热处理数值模拟进展之——扩展求解域热处理数值模拟.pdfVIP

热处理数值模拟进展之 扩展求解域热处理数值模拟 潘健生，王 婧，顾剑锋，郝晓伟 (上海市激光制造与材料改性重点实验室，上海交通大学材料科学与工程学院，上海200240) 摘要：金属热处理是复杂的多场耦合非线性物理过程。气-固、液-固界面的换热是热处理数值模拟必须首先解决的难题，也 是热处理工艺设计、热处理设备设计的基础。仅以热处理工件为研究对象、通过设定界面换热系数的来确定界面换热的传统方 法不叮避免会带来较大的误差，显著影响模拟精度。流体力学和传热学领域中有关“耦合换热”的一系列研究成果是破解现有 热处理模型中所遇到难题的有效途径。本文首先简要归纳前人有关“耦合换热”的研究成果，在此基础上首次提出了扩展求解 域的热处理模型。它将温度场的求解域扩展到加热或冷却的热处理工件及其周围环境。例如在加热过程中，考虑炉体、炉膛、 发热体、所有的炉内构件，如工装夹具，作为一个整体进行计算。炉膛空间内进行流场与温度场耦合计算，用壁面函数法处理 流场与固体之间换热，多物体多表面之间的辐射热交换用离散坐标法处理。固体区域的内部进行温度场与相变的耦合计算，各 物体(子域)界面上的热流密度和温度是耦合换热的计算结果而不是预设定的边界条件，避免了现有的热处理模型中预先设定 工件表面的边界条件所遇到的困难和小确定性，有助于更接近于实际生产的情况。最后，以新型罩式渗氮炉的设计为例，运用 扩展求解域的热处理数值模拟方法成功地进行了设备虚拟设计，证明了该模型具有的巨大的优越性。作为热处理数值模拟的重 要进展之一，扩展求解域热处理模型还很不成熟，尚有许多地方需完善。 关键词：热处理数值模拟数学模型 早在上世纪七十年代至八十年代初，Inoue【l】作为工件与环境之间的边界条件，这就需要预先给 Denis【21等分别建立了温度场、相变、应力场相互耦出换热系数。它包括对流换热系数和辐射换热系数 合的模型(又称为热一机械一冶金模型)为热处理的 两个部分，前者与边界上流体的流动状况密切相 数学建模和计算机模拟的研究奠定了基础，引起了 关，然而工件不同位置的表面附近的流速各不相 高度重视，此后本领域的研究迅速发展。各国学者 同，因炉内其他物体的存在而改变。在辐射换热系 先后发表了大量的研究成果，热处理模型不断得到 数中包含着环境温度与工件温度的三次方项，而工 改进并逐步走上实用化。 件表面的位置不同，外法线方向不同，与它辐射热 另一方面，几种商品化有限元分析软件的推 交换的“环境”也就各不相同。所以，预先给定一 出，为推动热处理计算模拟技术的研究与生产应用 个准确反映整个工件表面热交换状况的“换热系 创造了条件。在2000年第一届国际热加工建模与计数”常常是很困难的，甚至是不可能的。例如对图 算机会议上，T．Inoue【3】，Meel(ishoH和K．Funatani【5】1所示的钟罩式渗氮炉中的工件或是如图2所示的在 等对20世纪后期国际上热处理建模与计算机模拟的 一个真空加热高压气淬炉内的工件都不可能预先准 发展情况作了归纳和总结，进入新世纪以来，热处 确给出换热系数。 理模拟的研究继续深化之外，热处理计算机模拟在 另一方面在现有模型中一般将“炉膛”视作 实际生产中得到愈来愈多的应用【6。1u并有热处理模 “环境”。然而如图1和图2所示，在实际生产中 拟的专用软件推向市场【9l。 “炉膛”的状况(几何形状、尺寸、温度场、流 热处理是一个十分复杂的物理、化学过程。热 场)是很复杂的而且是多变的。因此将求解域局限 处理建模比其他的热加工过程(铸、锻、焊等)更 于工件本身的建模方法，就将面对一个复杂的带有 为复杂和困难，以致无论是模型的成熟程度或是计 不确定性的“环境”，以致于其边界条件也就将难 算机模拟的生产应用，热处理都明显落后于其他的 以确定。 热加工领域。在热处理模型的研究中，还有不少难 流体力学和传热学领域中有关“耦合换热”的 题有待解决，其中之一是现有热处理模型将求解域 一