实际钢件淬冷过程仿真与应用.doc

7 实际钢件淬冷与应用 7.1 问题的提出 除轴类、板类和盘类等简单形状的零部件外，实际工程中还常常涉及到复杂异形的工件材质选择及其热处理等问题。复杂异形钢件的热处理就不能采用前面所介绍的模型算法及应用平台处理。 近十年来，有限元或有限差分法淬冷温度场[191~193]计算技术的发展为复杂异形钢件的热处理模拟提供了强有力的工具。本章节以淬冷温度场有限元计算为手段，通过几个工程实例来揭示在材料设计或热处理工艺制订过程中的指导作用和意义。 7.2钢件淬冷温度场有限元计算法 7.2.1淬冷过程热传导方程 根据文献[192]，对于热物性指标呈各向同性的三维复杂形状淬冷钢件而言，其热传导平衡方程在直角坐标下的形式为： (7-1) 式中，ρ、c、λ和Q分别代表钢件材料的密度、比热容、热传导系数和相变潜热，其中，过冷奥氏体转变为铁素体—珠光体、贝氏体和马氏体的相变潜热可分别按305MJ/m3 、440MJ/m3和648MJ/m3选取[193]；钢件材料的密度、比热容c和热传导系数λ可按文献[194,195]选择；T和τ分别代表钢件温度和淬冷时间。就淬火过程而言，c和λ与钢件的具体温度有关，因此，式(1)为非线性导热方程。 钢件淬火时的初始条件为淬火温度，即 (7-2) 钢件淬火时，边界条件属于对流换热形式： (7-3) 其中h 、s和n分别表示表面换热系数、换热边界和换热边界的法向方向，Tw为工件表面温度，Tq为淬火介质温度。钢件水冷和油冷淬火时，其对流换热系数与工件表面温度之间关系可分别按文献[102,191]选取。除与空气的对流换热外，工件空冷淬火时还应有辐射换热存在。根据文献[196]，其综合换热系数按下式计算： (7-4) 式中，Tw和Tq单位均为绝对温标，K。 7.2.2求解步骤 1) 简化对象：基于钢件的对称性，将其剖分成N等份，选其中之一作为研究对象以减少计算量。 2)对象结构的离散化：将对象划分为e个有限单元，每个单元含p个节点。 3)选择插值函数：每一单元内对温度设定一个合适的插值形式Ni(x,y,z)，将单元e中任一位置的温度用各个结点的温度来表示，即 (7-5) 4)单元矩阵和载荷向量的推导应满足对流边界条件(见式7-3)的导热微分方程在每个单元上的积分形式，整理后所有单元的热平衡方程可集合成系统总的平衡方程组，具体可表示为下列矩阵形式[192]： (7-6) 式中，[K]为传导矩阵，包含导热系数、对流系数及辐射率等；[C]、{T}、和{q}分别为比热容矩阵、节点温度向量、温度对时间的导数和节点热流率向量(含相变潜热)。 采用向后差分将式(7-6)离散化 (7-7) 应用式(7-7)可以在每一时间段求出所有结点上的温度数值，再由插值函数计算出各单元的温度分布，从而实现钢件三维淬冷温度场的仿真。 7.3应用实例 7.3.1X38CrMo16模具钢模坯淬冷温度场仿真与参数选择 图7-1. X38CrMo16塑模钢模坯淬火裂纹(淬火温度1040℃) Fig.7-1 The crack occurred in X38CrMo16 plastic mould steel bloom quenched at 1040℃ 因具有高的耐腐蚀性能和良好的抛光性能，马氏体型不锈钢X38CrMo16(材料号为DIN1.2316)通过热处理调质(1020~1050℃油淬+600~650℃)预硬化处理至32~35HRC后主要用于制作要求耐腐蚀的塑料模具，如PVC、再生塑料和橡胶制品挤压模。然而，由于合金化程度高和导热性能较差，大规格X38CrMo16模坯淬冷过程中因油冷时间掌握不当，淬裂(油冷时间过长，见图7-1中所示)或模坯心部未淬硬(油冷时间过短)现象时有发生。故有必要对大型X38CrMo16模坯的淬冷参数进行研究、优化和选择。本实例以X38CrMo16塑料模具钢的连续冷却转变曲线为准则，以瞬态热传导平衡方程为计算工具，在对尺寸为1700mm╳280mm╳650mm的模坯在油冷介质条件下的淬冷温度场仿真基础上，分析并确定了其最小淬冷时间和油淬空冷转换时机以及最大许可尺寸，为获得无淬火缺陷的模坯提供了理论依据。 1) 模坯淬冷温度场计算过程 a) 简化对象：基于长方体模坯的对称性，将其剖分成8等份，见图7-2(a)所示，选其中之一作为研究对象，坐标原点定在模坯的心部中心处。 b) 对象结构的离散化：见图7-2(b)中所示。 c)按式7-7计算：X38CrMo16模具钢CCT曲线见图7-3中所示，密度ρ=7.870g/cm2[197]、定压比热容c和热传导系数λ与温度的关系可近似按3Cr13