淬火作为改变和提高材料性能的重要手段,在汽车、航空、运输、建筑及其他制造业领域有着广泛而重要的应用。但由于淬火过程是一个复杂的多场耦合过程,传统热处理技术已经不能适应机械工业现代化的要求,应用计算机技术改造热处理已成为当务之急。在这种背景下热处理数学模型与数值模拟技术日益受到重视,已成为当今各国研究开发的热点。上世纪七八十年代以来,计算机技术的迅猛发展,使计算机作为一个强有力的工具在各工程领域获得了越来越广泛的应用,淬火领域也不例外。淬火过程的数值模拟作为一种计算机辅助工艺设计方法,可以优化设计工艺、降低废品率,从而大大提高工作效率。淬火过程计算机模拟技术近年来虽然取得了一些进展,并显示出巨大的优越性,但它仍处于发展阶段,许多重要问题如换热系数测算、温度组织应力耦合分析模型建立等有待于进一步研究,它的潜力也没有真正发挥出来。本文针对上述淬火模拟存在的问题,建立了淬火冷却过程三维温度场、组织转变和应力/应变场的耦合分析模型,提出了基于有限元和最优化方法计算淬火换热系数的方法,把遗传算法引入气体淬火工艺参数优化过程,得到了可以提高淬火质量的淬火工艺参数。具体如下;针对淬火冷却过程温度和组织变化的特点,综合考虑了边界换热系数、相变潜热以及热物性参数等非线性因素的影响,采用C++语言编写了求解淬火冷却过程三维温度场的有限元程序。该程序大大减小了通用有限元软件处理淬火过程非线性问题时只能以常数代替所引起的模拟误差。提出了相变潜热的计算与处理方法,把相变潜热作为每个单元的内热源进行处理,并迭代求解至收敛,提高了温度场的计算精度。编写了模拟相变过程和组织转变的程序。研究了网格细化、动态时间步长和集中热容矩阵对求解温度场时出现的数值振荡问题的影响,计算结果表明;在时间步长选择合适的情况下,集中热容矩阵能够很好地解决数值振荡问题;如果不进行集中热容矩阵,无论是细化网格还是减小时间步长,其解决数值振荡问题的效果都不好。但是集中热容矩阵造成边界节点的温度曲线与解析解的误差较大。自主设计了端部淬火实验工装,对P20试样进行端部淬火实验,并对试样相应部分进行金相分析及硬度检验,得到了准确的实验数值。建立了P20钢端部淬火的有限元模型,用本文编写的程序对P20钢端部淬火实验进行了有限元模拟,得到了相应实验位置的各新生成相的含量及相应位置的硬度计算值,将这些数据与实验数据进行比较,对有限元模拟程序进行了检验,检验表明本文的有限元程序是准确的、可靠的。结合非线性有限元方法和最优化方法,对淬火冷却过程的反向热传导问题进行了研究。利用文献中的换热系数数据和热物性参数,计算得到淬火零件相应位置的温度变化曲线,然后根据温度变化曲线,用本文开发的基于有限元和最优化方法的反传热分析程序求解随温度变化的换热系数。利用三维非线性有限元方法,以进退法确定初始搜索区间,结合黄金分割法、逐次二次插值法及其混合法等优化方法,根据已知温度曲线,对淬冷过程反传热中的换热系数的求解时间问题进行了研究。在求解换热系数时,对端淬和全淬的情况进行了比较。由于在两种情况下首先计算换热系数的点选取不同,在端淬情况下所得换热系数与参考值吻合程度好于全淬的情况。自主设计了采集试样内部相应位置处冷却曲线的实验装置及基于ISA总线的高速数据采集系统,利用该系统采集了P20钢在20℃和60℃水中冷却时的冷却曲线,根据实验测量的温度曲线利用本文中的有限元模拟程序和最优化程序得到了20℃和60℃水随淬火件温度变化的换热系数曲线。推导和建立了淬火过程材料处于弹性区、塑性区、过渡区的应力/应变关系,并编写了相应的三维有限元模拟程序。根据淬火过程的特点,提出了一种新的计算弹塑性比例方法,并提出了加快弹塑性比例系数收敛速度的方法。计算结果表明;该方法可以明显地加快迭代时弹塑性比例系数的收敛速度,避免产生不收敛现象,而且可以较好的保证计算精度。采用C++语言编写了淬火过程的三维应力应变场的模拟程序。选择1080钢水冷淬火和板的热应力问题两个算例对该程序进行了验证。结果表明,模拟结果与实验结果比较接近,用本文的有限元模拟程序,可以得到比较准确的应力应变值。建立了温度场、组织转变和应力场的三维耦合分析模型。选择轴对称空心圆柱体问题、偏心圆柱筒为研究对象,建立相应的有限元分析模型,研究了其温度场、组织转变和应力场在淬火过程中的变化情况,并将其组织分布情况与实验结果进行了比较。研究了三维偏心圆柱筒由温度引起的初应变、相变引起的初应变及相变塑性对淬火残余应力的影响。把遗传算法引入气体淬火工艺参数优化过程。提出了阶段性换热系数模型,并根据该模型以淬火变形为目标对P20钢砖型结构件气体淬火工艺进行了优化,以淬火硬度、残余应力等作为评估函数,对优化前后钢件的淬火结果进行了比较,得到了可以提高淬火质量的淬火工艺参数。