高频感应淬火仿真

感应淬火的物理与理论基础

趋肤效应决定一切

感应淬火是利用感应线圈的交变磁场在工件表层激起的涡流所产生的焦耳热，只把表面快速加热，紧接着急冷形成马氏体硬化层的工艺。支配加热空间分布的正是趋肤效应：电流密度自表面起按透入深度 \( \delta \) 指数衰减。

$$ \delta = \sqrt{\frac{2\rho}{\mu \omega}} = \sqrt{\frac{\rho}{\pi \mu f}} $$

对室温下的钢，\( f = 10 \) kHz时 \( \delta \) 为0.1 mm量级，\( f = 1 \) kHz时略小于1 mm——频率的选择本身就是硬化层深度的设计变量。这里还要加上感应加热特有的一层曲折：温度升高时电阻率 \( \rho \) 增大，而在居里点（钢约770 °C）铁磁性消失、\( \mu \) 骤降，透入深度会一下子跳到数倍。也就是说，加热过程中"被加热的位置"本身在动态变化——这种强非线性，正是感应淬火仿真区别于单纯热分析的核心。

电磁、热、组织的三物理耦合

由于电源频率与加热时间的尺度相差6个数量级，标准做法是把电磁场用频域（时谐）分析求出单周期平均的发热密度，再把它作为热分析的源项，两者分离求解。整体则构成弱耦合迭代：温度每变化一次就更新物性（ρ、μ）并重新求解电磁场。

淬火侧的物理——何谓"硬化"

硬化的条件有两级。①加热要超过奥氏体化温度（Ac3，随升温速度上移）；②冷却速度要高于临界冷却速度，从而发生马氏体相变（Ms点以下按Koistinen-Marburger定律）。因此仿真的交付物不是温度云图，而是"完成奥氏体化的区域∩能够急冷的区域"＝硬化层形貌与硬度分布。冷却由喷淋或浸淬的沸腾传热（温度依赖极强的沸腾曲线）支配，这一侧同样是强非线性。

数值分析方法

电磁-热耦合的实现

1. 电磁场（时谐分析）——给定线圈电流（或电源功率），计算涡流与发热密度分布。网格的铁律是趋肤深度内至少布置2～3个单元（以居里点之前那个较薄的δ为准）
2. 热分析（瞬态）——以发热密度为源项推进温度
3. 物性更新与再求解——随温度升高更新ρ与μ（B-H曲线），重新求解电磁场。更新间隔按温度变化控制，尤其是通过居里点时
4. 冷却与相变分析——用沸腾传热急冷，并基于CCT／TTT数据的相变模型（扩散型用JMAK系，马氏体用KM式）计算组织分数与硬度

收敛与精度的要害

• 通过居里点——μ骤降会使发热分布剧变。这一区间若不细化物性更新（温度步长10～20 °C），加热历程就会算偏
• B-H非线性——近表面区域工作在磁饱和状态。线性μ近似会高估发热量。实测的B-H曲线（含温度依赖）是精度的地基
• 移动（扫描淬火）——线圈或工件运动时，可用逐位置的时谐分析组合成移动热源，或在随动坐标系中按定常问题处理
• 电源的建模——实机采用恒功率／恒流控制加匹配电路。若"线圈电流恒定"的假设与实机不符，就要把电源控制纳入迭代回路

冷却换热的标定

喷淋冷却的换热系数依赖表面温度（膜态沸腾→过渡沸腾→核态沸腾），变化超过一个数量级，文献值的离散也很大——冷却侧的h是应当与加热侧独立标定的参数。实务上是用试件内埋设热电偶测得的冷却曲线，反向辨识h的温度依赖曲线（与回流焊辨识h的思路相同，只是非线性更强）。

实务应用指南

标准工作流程

1. 材料数据整备——温度相关的ρ、B-H曲线与热物性，CCT曲线，以及（若要看变形）高温力学物性。这里是最大的准备工作量
2. 加热的标定——按实机条件加热试件，把电磁-热模型拟合到红外测温仪或热电偶的温度历程上（辨识有效电流与效率）
3. 冷却的标定——由冷却曲线辨识沸腾h曲线
4. 硬化形貌的预测与验证——与截面宏观金相、硬度分布（维氏硬度扫描）对比。硬化层深度±10～15%左右，是标定过的模型的常见水平
5. 条件寻优——通过频率、功率、加热时间、线圈形状、冷却条件的参数研究，锁定既得到目标硬化层又变形最小的条件

仿真发挥作用的典型场景

• 感应线圈设计——对阶梯轴、齿轮这类形状，权衡棱边过热（角部磁通集中，有淬裂风险）与齿根未硬化，确定线圈形状与耦合间隙。减少试制线圈的效果很大
• 新产品的条件确定——形状变更后无法沿用同类件条件时，用于缩小初始条件的范围
• 变形与开裂对策——预测相变膨胀＋热应力造成的变形与残余应力，识别淬裂（拉伸残余应力×急冷）的高风险部位

与实测对照的做法

验证的基准数据，最有力的是"硬度分布的截面图"。它比温度历程的点测量信息量更大，能同时检验加热（奥氏体化区域）与冷却（达到临界冷却速度的区域）两侧的合理性。要看的不只是硬化层深度，还有形貌的形状（棱边处的绕入、齿根处的中断）是否吻合——若深度对得上而形貌不同，就提示电磁侧（线圈模型、耦合间隙）有误。

主要工具的支持情况

工具对比

选型的思路

按三点决定：①电磁侧的精度要求（若要做到线圈优化，就需要电磁专业级的B-H处理能力），②是否要一路贯通到组织、硬度与变形（是否需要DANTE级的材料模型），③使用者是谁（专职分析人员，还是线圈设计现场）。有一条共通的提醒：工具自带的材料数据库往往与自家钢种对不上，至少B-H曲线与CCT数据要靠自家整备，这是精度的前提。

前沿动向

感应线圈形状的自动优化

以目标硬化形貌为约束，对线圈形状、匝数、耦合间隙做参数化或拓扑式优化的研究与实现正在推进。由于电磁-热耦合的单个算例已能压缩到分钟量级，具有实用迭代次数的优化成为可能。再结合3D打印制造复杂线圈，"把计算设计出的线圈原样做出来"正在成为现实。

与数字孪生、工艺控制的融合

把电源的实测波形与红外测温实时同化进模型，以补偿波动（材料批次、线圈磨损、冷却介质劣化）的控制研究正在进行。物理模型被缩约为降阶模型或代理模型以实现毫秒级响应，进而通向全数质量保证（把每个工件的推定硬化层深度记录下来）——这一构想正是产业侧的驱动力。

组织与特性预测的深化

处理快速短时加热特有的非平衡奥氏体化（Ac3随升温速度变化、碳化物未完全溶解）的动力学模型，原奥氏体晶粒尺寸与硬度、疲劳强度的关系，以及包含残余应力在内直至疲劳寿命预测的一贯耦合，都是研究前沿。预测的目标正从"硬化层深度"向"零件性能（疲劳极限）"再深入一层。

故障排查

按症状的原因与对策

算不准时，先怀疑哪个物理

相关文章：工业热过程文章一览、电弧焊仿真（同样是移动热源与相变的近亲）、回流焊炉仿真。