关于2.25Cr1Mo钢大型锻件热处理工艺数值模拟研究的学术报告
本研究报告基于中国原子能科学研究院龚雪婷、武志广、李鑫以及辽宁科技大学贾红帅共同完成，并于2019年3月发表于《Heat Treatment of Metals》（金属热处理）第44卷第3期的一篇学术论文。该研究题为“数值模拟2.25Cr1Mo钢大型锻件热处理工艺数值模拟”（Numerical simulation on heat treatment process of 2.25Cr1Mo steel large forgings）。
一、 研究的学术背景
本研究属于材料科学与工程领域，具体聚焦于大型金属锻件热处理工艺的数值模拟与优化。随着重型装备制造业（如核电、石化）的发展，对大型锻件的性能、组织均匀性和稳定性提出了极高要求。热处理是决定大型锻件最终性能的关键工序。然而，大型锻件尺寸巨大、结构复杂、壁厚不均的特点，使得其热处理过程温度场和组织场极不均匀，传统上往往需要通过成本高昂的整体解剖试验来摸索工艺参数，效率低且浪费严重。
在此背景下，有限元数值模拟技术为预测和优化大型锻件热处理工艺提供了强有力的工具。通过计算机模拟，可以在实际生产前预测锻件在加热和冷却过程中的温度变化与组织演变，从而优化工艺参数，提高产品质量并降低成本。
本研究的具体研究对象是2.25Cr1Mo钢，这是一种广泛应用于核电、火电及石化领域高温高压部件的铁素体-珠光体型耐热钢。其性能高度依赖于热处理后获得的显微组织。在不同的冷却条件下，该钢可能形成贝氏体、珠光体、铁素体或其混合组织。对于壁厚差异大的大型锻件（如研究中涉及的“深碗状”管板锻件），若不精确控制冷却条件，极易导致心部与表面组织不一致，从而影响产品整体性能的均匀性。
因此，本研究的主要目标在于：利用有限元模拟软件，结合实验测得的连续冷却转变（Continuous Cooling Transformation, CCT）曲线，建立2.25Cr1Mo大型锻件热处理过程的温度场与组织场模型。通过模拟水冷、油冷和空冷三种典型冷却工艺，探究锻件内部温度分布与组织演变规律，特别是心部“难温降”区域的冷却行为，进而评估并优化热处理工艺，以期为获得以理想贝氏体组织为主的均匀组织提供理论指导和工艺参考。
二、 详细研究流程
本研究的工作流程是一个典型的“建模-模拟-分析-验证-优化”循环，具体步骤如下：
第一步：研究对象与材料基础数据准备 研究以某大型2.25Cr1Mo管板锻件为物理背景，其毛坯呈“深碗状”，最大壁厚达480mm，内外直径差750mm，这种几何特征导致冷却时内外温差大，组织均匀性控制困难。研究使用的材料化学成分如表1所示。最关键的基础数据是研究团队通过膨胀法与金相-硬度法相结合实测得到的2.25Cr1Mo钢的CCT曲线（图3）。该曲线明确了不同冷却速率下相变的产物：冷却速度小于0.083℃/s时，转变为铁素体-珠光体；在0.083～3.33℃/s之间时，以贝氏体组织为主；大于3.33℃/s时，转变为马氏体。图4展示了不同冷速下的实际显微组织，为后续组织模拟提供了真实参照。
第二步：有限元模型的建立 1. 几何建模与网格划分：基于锻件图纸（图1）建立三维几何模型。为了精确捕捉温度梯度和组织变化，研究采用了六面体单元进行网格划分，基本单元尺寸为80mm x 80mm x 80mm，将整个锻件划分为9440个网格单元（图2）。这种细致的划分旨在获得更详细的微观组织模拟结果。 2. 温度场模型与边界条件定义：温度场模拟的核心是定义锻件表面与冷却介质之间的换热系数。研究区分了不同冷却方式的换热模型： * 水冷及油冷：考虑了淬火过程中可能存在的膜沸腾、核沸腾和对流换热三个阶段，指出大锻件淬火时介质的搅动会显著影响换热系数规律，但未给出具体拟合公式，暗示了实际工艺的复杂性。 * 空冷：将换热系数作为变量处理，由辐射换热和对流换热两部分构成。研究采用了匡琦、潘健生等人提出的经验公式（公式3）来计算空冷换热系数，并给出了在15℃室温下，界面温度从20℃到900℃变化时对应的综合换热系数表（表2）。这为模拟提供了量化的边界条件。 3. 组织场模型建立：将实验测得的CCT曲线（图3）导入有限元模拟软件中（转化为软件可识别的格式，如图5）。软件将根据模拟计算出的每个单元在冷却过程中经历的温度-时间曲线（即冷却曲线），与CCT曲线进行比对，从而预测该单元最终转变的组织类型（铁素体-珠光体、贝氏体或马氏体）及比例。这是连接物理场（温度）与最终性能（组织）的关键桥梁。
第三步：不同冷却工艺的数值模拟与结果分析 研究对水冷、油冷和空冷三种工艺进行了系统的模拟。 1. 温度场模拟：模拟结果显示，在冷却过程中，锻件表面温度下降快，心部温度下降慢，形成了从表面到心部逐渐升高的温度梯度（图6）。特别是锻件1/4壁厚区域（模拟图中黑框标示）及“碗底”壁厚中心，由于热量传递路径长，成为“难温降”区域。 2. 组织场模拟： * 水冷：冷却速度最快。模拟结果显示（图7a），锻件内外表面形成了大量的贝氏体组织，仅在1/4处难温降区域及壁厚中心出现少量铁素体-珠光体组织。 * 油冷：换热系数小于水冷。模拟结果（图7b）显示，仅在内外表层接触介质处形成少量贝氏体，距表面约220mm以内的区域开始出现铁素体-珠光体，冷却结束后锻件主体以铁素体-珠光体组织为主。 * 空冷：冷却速度最慢。模拟结果（图7c）显示，整个锻件降温缓慢，最终全部转变为铁素体-珠光体组织。 3. 关键部位冷却行为深入分析——质点追踪法：为了更精确地对比锻件不同部位的冷却历程，研究采用了“质点追踪”方法。在空冷模拟中，从表面到心部1/4处选取了5个追踪点（图8a），绘制其温度-时间曲线（图8b）。曲线清晰表明，靠近表面的点（1、2点）冷却快，而位于心部区域的点（3、4、5点）冷却缓慢。计算关键点（2、3、5点）在不同冷却方式下的冷却速率（表3）。数据显示，在距外表面约210mm（2、3点）的位置，冷却速率随介质不同变化显著（从空冷0.054℃/s到水冷0.330℃/s），说明该厚度范围内的材料对冷却介质敏感；而在更深的300-360mm处（5点附近），冷却速率受介质影响小，且在所有冷却方式下均低于形成贝氏体所需的最低临界冷却速率（约0.078℃/s，来自CCT曲线），这解释了为何在心部难以获得纯贝氏体组织。
第四步：工艺优化探索 基于上述发现，研究指出，理想的均匀贝氏体组织介于马氏体（过快冷却，韧性差）和铁素体-珠光体（过慢冷却，易辐照脆化）之间。传统单一水冷虽能在外层获得贝氏体，但难以确保心部完全转变。因此，研究团队利用有限元模拟进行了反向优化探索：通过反复调整和验证模拟中的换热系数，发现当将综合换热系数提升至35000-45000 W/(m³·K)范围时，模拟结果显示（图9），锻件绝大部分区域可转变为贝氏体，仅在最难冷却的1/4区域残留极少铁素体-珠光体。这启示了实际工艺改进方向：需要采取强制冷却措施（如锻件倒置、内孔喷水、加强水循环和流速）来打破淬火时表面的蒸汽膜，显著提高有效换热系数，从而促进心部冷却，实现整体组织的贝氏体化。
三、 主要研究结果
四、 研究结论与价值
本研究得出以下主要结论： 1. 利用有限元数值模拟技术，可以有效建立2.25Cr1Mo钢大型复杂锻件热处理过程的温度场与组织场预测模型。 2. 对于所研究的大型管板锻件，传统水冷是三种冷却方式中最优的，能使其表面和大部分区域获得所需的贝氏体组织，但仍无法保证心部完全转变。 3. 为实现锻件整体均匀的贝氏体组织，必须在传统水冷基础上进一步强化冷却过程，例如通过强制循环、喷淋等方式提高换热系数，以加速心部冷却，使其冷速进入贝氏体转变区。
本研究的价值体现在： * 科学价值：将CCT曲线这一基础相变数据与有限元数值模拟深度结合，实现了从冷却条件到微观组织的定量化、可视化预测，深化了对大型锻件热处理过程中传热与相变耦合规律的理解。 * 应用价值：为2.25Cr1Mo钢大型锻件的热处理工艺制定提供了重要的预研工具和优化方向。通过“模拟先行”，可以减少甚至替代昂贵的实物解剖试验，显著降低研发成本和周期，提高产品质量的可靠性和一致性。提出的提高换热系数的具体思路对实际生产具有直接的指导意义。
五、 研究亮点
六、 其他有价值内容
研究在讨论部分坦率指出了当前模型的局限性：基于理论或经验公式的水冷换热系数与实际生产（尤其是带有搅动、循环的淬火过程）存在偏差。作者提出将在后续工作中通过实际生产数据来修正换热系数模型，以提高模拟的精确度。这种实事求是的科学态度和对模型持续改进的追求，是数值模拟研究走向工程实用化的关键一环。此外，研究引用的空冷换热系数经验公式及对换热机理的简要讨论，也为相关领域的研究者提供了有价值的参考。