辙叉淬火过程多相流数值模拟及温度分析
 
 摘要
 
辙叉是铁路轨道结构中承受最大动载荷的关键部件，其淬火质量直接影响行车安全与服役寿命。高锰钢导热性能较差，导致辙叉在淬火时心部冷却速度较慢，容易在心部产生超标碳化物，是制约辙叉产品质量提升的核心瓶颈。生物化工以高锰钢辙叉水韧淬火过程为研究对象，系统分析了辙叉淬火多相流数值模拟的理论基础与技术方法，从多相流控制方程、流固耦合传热机制及相变动力学三个层面构建了辙叉淬火过程的数值模拟框架。在此基础上，探讨了多相流模拟在辙叉淬火设备优化中的工程应用，揭示了多相流场对辙叉温度分布及冷却均匀性的影响规律。研究表明，多相流数值模拟技术能够有效预测辙叉淬火过程中的温度场演变规律，为淬火工艺参数优化与冷却设备改进提供重要的理论依据与技术支撑。
 
关键词：辙叉；淬火；多相流；数值模拟；温度场
 
 一、引言
 
辙叉是铁路轨道结构中实现车轮从一股钢轨转移到另一股钢轨的关键部件，服役过程中承受着列车轮载的反复冲击与剧烈摩擦，工作条件极为严酷。高锰钢因其优异的加工硬化特性与良好的韧性，成为辙叉制造的主导材料。然而，高锰钢辙叉的质量在很大程度上取决于其热处理工艺，尤其是水韧淬火过程的控制。高锰钢导热性能较差，导致辙叉在淬火时心部冷却速度较慢，容易在心部产生超标碳化物，严重削弱了辙叉的力学性能与服役寿命。
 
在淬火过程中，高温辙叉浸入水中后，工件表面与冷却介质之间发生剧烈的热交换，涉及液体介质的流动、汽化相变以及气-液-固三相之间的复杂相互作用，构成典型的多相流动与传热问题。传统“试错法”依赖经验调整工艺参数，开发周期长、成本高，且难以实现对淬火过程温度场的精准控制。由于淬火试验存在高温作业、在线测量困难、成本高等问题，淬火过程的数值模拟成为目前淬火技术研究的重要手段。
 
数值模拟技术能够突破物理试验的时空限制，以较低的成本获取工件内部温度场、组织场和应力场的完整演变信息。近年来，随着计算流体动力学（CFD）和有限元方法的快速发展，多相流数值模拟已逐渐成为淬火工艺优化与设备改进的核心技术手段。本文旨在系统阐述辙叉淬火过程多相流数值模拟的理论基础、技术方法及其工程应用，以期为辙叉淬火工艺的科学化设计与精准化控制提供理论支撑。
 
 二、辙叉淬火工艺与多相流物理特征
 
 2.1 辙叉淬火工艺特点
 
高锰钢辙叉的水韧处理是将铸态辙叉加热至奥氏体化温度（通常为1000～1100℃）并保温足够时间后，迅速浸入水中快速冷却的工艺过程。水韧处理的目的是使碳化物充分固溶到奥氏体中，获得单一的奥氏体组织，从而保证辙叉在使用过程中具有良好的加工硬化能力和韧性。
 
高锰钢导热性能较差，这一材料特性给辙叉淬火带来了根本性的技术难题。大型辙叉在淬火时，表面与心部的冷却速度存在显著差异：表面与冷却水直接接触，冷却速度较快；心部则因材料导热系数低而冷却迟缓。当心部冷却速度不足时，过饱和的碳原子将从奥氏体中析出，形成粗大的碳化物，导致心部组织中出现超标碳化物。热模拟试验结果表明，随着冷却速度的增加，高锰钢析出碳化物逐渐减少；在5℃/s的冷却速度下，析出碳化物级别可达到X2级。
 
 2.2 淬火过程的多相流物理特征
 
辙叉浸入水中的淬火过程涉及极为复杂的物理现象。高温辙叉（约1000℃）浸入常温水中后，江西化工工件表面附近的水瞬间被加热至沸腾，形成大量的水蒸气气泡。这些气泡脱离工件表面向上运动，在工件周围形成气-液两相混合的流动区域。随着冷却的进行，工件表面温度逐渐降低，沸腾强度减弱，气泡生成速率下降，流动逐渐从强烈的核态沸腾过渡到温和的对流换热阶段。
 
从多相流的角度来看，辙叉淬火过程呈现出以下典型特征：
 
气-液两相共存与演化。淬火初期，工件表面强烈的核态沸腾产生大量气泡，形成气-液两相流。气泡的运动、聚并、破裂及其与工件表面的相互作用，直接决定了热交换的强度与均匀性。
 
流固耦合传热。工件内部的热量通过导热方式传递至表面，再由表面对流传热传递给冷却水，同时伴随沸腾传热。工件的温度分布既受材料热物性参数的影响，也受表面换热条件的制约，形成典型的热-流耦合问题。
 
剧烈的温度梯度。由于高锰钢导热系数低且淬火冷却速度快，辙叉截面内将形成极大的温度梯度，这种温度梯度是导致组织不均匀和残余应力的根本原因。
 
蒸汽膜的形成与破坏。在淬火初期，若工件表面形成的蒸汽膜不能及时被水流冲散，将严重阻碍热交换，造成局部冷却不足。多相流耦合搅拌系统正是通过增强水流湍动来破坏辙叉底部空腔的蒸汽膜，从而保证均匀的冷却效果。
 
 三、多相流数值模拟的理论与方法
 
 3.1 多相流控制方程
 
辙叉淬火过程中的多相流动与传热问题，需要基于两相流体动力学及流固耦合传热建立完整的数学描述。由于淬火过程中涉及气-液两相的剧烈相互作用，通常采用欧拉-欧拉双流体模型或欧拉-拉格朗日混合模型进行描述。
 
连续性方程描述各相的质量守恒。对于不可压缩的液相，其连续性方程为：
 
$$\frac{\partial}{\partial t}(\alpha_l \rho_l) + \nabla \cdot (\alpha_l \rho_l \mathbf{u}_l) = 0$$
 
其中，$\alpha_l$为液相的体积分数，$\rho_l$为液相密度，$\mathbf{u}_l$为液相速度矢量。对于气相，由于气泡的体积可随压力变化，其连续性方程需考虑可压缩性效应。
 
动量方程描述各相的动量守恒。液相动量方程可写为：
 
$$\frac{\partial}{\partial t}(\alpha_l \rho_l \mathbf{u}_l) + \nabla \cdot (\alpha_l \rho_l \mathbf{u}_l \mathbf{u}_l) = -\alpha_l \nabla p + \nabla \cdot (\alpha_l \boldsymbol{\tau}_l) + \alpha_l \rho_l \mathbf{g} + \mathbf{F}_{lg}$$
 
其中，$p$为压力，$\boldsymbol{\tau}_l$为液相粘性应力张量，$\mathbf{g}$为重力加速度，$\mathbf{F}_{lg}$为气-液相间作用力（包括曳力、升力、虚拟质量力等）。
 
能量方程描述各相的能量守恒。对于液相，能量方程为：
 
$$\frac{\partial}{\partial t}(\alpha_l \rho_l h_l) + \nabla \cdot (\alpha_l \rho_l \mathbf{u}_l h_l) = \nabla \cdot (\alpha_l k_l \nabla T_l) + Q_{lg}$$
 
其中，$h_l$为液相焓值，$k_l$为液相导热系数，$T_l$为液相温度，$Q_{lg}$为气-液相间传热。
 
 3.2 流固耦合传热模型
 
辙叉工件内部的传热由导热控制，其能量方程为：
 
$$\rho_s c_p \frac{\partial T_s}{\partial t} = \nabla \cdot (k_s \nabla T_s) + \dot{q}$$
 
其中，$\rho_s$为材料密度，$c_p$为比热容，$T_s$为工件温度，$k_s$为导热系数，$\dot{q}$为相变潜热释放率。该方程与流体域的能量方程在工件表面通过换热边界条件耦合：
 
$$-k_s \frac{\partial T_s}{\partial n} = h(T_s - T_l)$$
 
其中，$h$为表面换热系数，它不仅是温度的函数，还受沸腾状态、气泡行为和流体流动状态的强烈影响。精确确定换热系数是多相流数值模拟的关键难点。
 
 3.3 数值求解方法
 
辙叉淬火多相流数值模拟通常采用有限体积法（FVM）或有限元法（FEM）。对于流体域，FVM因其良好的守恒特性而成为主流选择，常用软件包括ANSYS Fluent、STAR-CCM+等。对于固体域，FEM在应力分析方面具有优势，常用软件包括ABAQUS、ANSYS Mechanical等。
 
在求解策略上，通常采用分区耦合或整体耦合方法。分区耦合方法分别求解流体域和固体域，通过数据交换实现耦合；整体耦合方法则将流体域和固体域的控制方程联立求解。由于淬火问题的强非线性和多时间尺度特征，瞬态求解是必要的，时间步长需兼顾计算精度和计算效率。
 
对于多相流模拟中的气泡行为，可采用群体平衡模型（PBM） 描述气泡的成核、生长、聚并和破碎过程。通过耦合PBM与双流体模型，可以更准确地预测淬火过程中的气含率分布和界面面积密度，从而提高传热预测的精度。
 
 四、多相流模拟在辙叉淬火中的应用
 
 4.1 温度场分布的预测与分析
 
多相流数值模拟的核心目标之一是获取辙叉在淬火过程中任意时刻、任意位置的温度分布。通过求解前述耦合方程，可以获得辙叉截面内温度随淬火时间的演变规律。
 
模拟结果通常揭示以下关键信息：心部温度下降速率远低于表面，这种温差在淬火初期最为显著；辙叉底部和凹槽等复杂几何区域因水流受阻，冷却速度往往低于其他部位，形成“冷却死区”，这是导致局部碳化物超标的主要原因之一。通过对温度场的分析，可以精确定位辙叉内部的“慢冷区”，为冷却设备的结构优化提供依据。
 
 4.2 多相流场对冷却均匀性的影响
 
多相流场的分布直接决定了辙叉表面的换热条件。数值模拟可以揭示以下重要规律：
 
气泡行为对换热的影响。淬火初期，气泡在辙叉表面成核、生长和脱离，显著增强了局部扰动，使换热系数大幅提高。但随着气泡在复杂几何区域的聚集，可能形成阻碍换热的“气膜层”，导致局部换热恶化。
 
水流速度与方向的影响。水流速度的提高能够增强对流换热，同时有助于冲刷表面气泡，破坏蒸汽膜。但过高的流速可能导致局部冷却过快，增加开裂风险。模拟可以优化水流速度和方向，实现“快而不裂”的理想淬火效果。
 
多相流耦合搅拌的工程价值。最新工程实践表明，在淬火水槽中引入多相流耦合搅拌系统，通过文丘里结构利用负压效应增强水流湍动，能够有效破坏辙叉底部空腔的蒸汽膜。60米长的超大水韧池在设计上首次采用多相流模拟技术，保证最优的水冷却效果。
 
 4.3 基于模拟结果的工艺优化
 
数值模拟的核心价值在于指导工艺优化。通过多相流数值模拟，可实现以下优化目标：
 
优化冷却设备结构。基于模拟确定的“慢冷区”位置，调整冷却水入口和溢流口的布置，增加局部水流速度。改进后高锰钢心部冷却速度可达到5℃/s，心部析出碳化物级别由X3级提升到X2级。
 
优化淬火工艺参数。通过模拟不同入水速度、水温、水流速度下的温度场演变，确定最优工艺参数组合，兼顾冷却速度与温度均匀性。
 
预测组织与性能。将温度场模拟结果与TTT/CCT曲线耦合，可以预测不同位置的最终组织，从而预判力学性能。
 
 五、结论
 
辙叉淬火过程涉及高温工件与冷却介质之间的剧烈热交换，伴随沸腾相变、气液两相流动以及流固耦合传热等复杂的多相流物理现象。高锰钢导热性能较差的材料特性与大型辙叉复杂的几何结构叠加，使心部冷却速度不足成为制约辙叉产品质量的核心瓶颈。
 
本文从多相流控制方程、流固耦合传热机制及数值求解方法三个层面，系统阐述了辙叉淬火多相流数值模拟的理论框架与技术路径。研究表明，多相流数值模拟能够有效揭示辙叉淬火过程中温度场的时空演变规律，精确定位心部“慢冷区”，为冷却设备的结构改进与淬火工艺参数的优化提供科学依据。通过基于模拟结果的设备优化，可使高锰钢心部冷却速度达到5℃/s，心部析出碳化物级别由X3级提升到X2级。最新工程实践表明，采用多相流模拟技术指导设计的超大水韧池与多相流耦合搅拌系统，能够有效破坏蒸汽膜，实现均匀、高效的冷却效果。
 
随着计算流体动力学与有限元方法的持续进步以及高性能计算能力的不断提升，多相流数值模拟将在辙叉淬火工艺的精细化设计与智能化控制中发挥日益重要的作用，为高锰钢辙叉产品质量的持续提升提供坚实的技术支撑。
 
 
 参考文献
 
[1] 张晓宇, 袁勇堂, 吴添奇, 等. 辙叉淬火过程多相流数值模拟及温度分析[J]. 锻压装备与制造技术, 2026(3).
 
[2] 柳亚楠, 刘洪武, 李学通. 高锰钢辙叉淬火冷却过程模拟及淬火设备优化[J]. 金属加工（热加工）, 2024(3): 83-86,90.
 
[3] 曹瑞, 孙会. 淬火过程数值模拟技术的研究进展[J]. 材料导报, 2015(5).
 
[4] TB/T 447-2004, 高锰钢辙叉技术条件[S].
 
[5] GB/T 5680-1998, 奥氏体锰钢铸件[S].