基于物理的16MnCr5低压渗碳模型：渗层深度对比与参数可辨识性

《Journal of Manufacturing and Materials Processing》：A Physics-Based Low-Pressure Carburizing Model for 16MnCr5 with Case-Depth Comparisons and Parameter Identifiability

　　

低压渗碳(LPC)模型可支持工艺配方筛选，但数值一致性必须与实验验证和校准区分开来。研究人员采用包含浓度依赖性?gren扩散系数、预设增强期碳通量、零通量扩散窗口、碳化物储库和非等温尾段的一维模型，与16MnCr5钢在940–950 °C条件下的四次工业试验进行了比较。数值验证包括解析解、网格细化以及独立的COMSOL实现。在未针对本批次拟合任何参数的情况下，两个自由表面试验在0.36 wt.% C处的深度分别为894和866 μm，与在HV0.3下测得的基于硬度的CHD₅₅₀值1036和907 μm相比，点偏差为?13.7%和?4.6%。联合通量拟合的残差为?4.9%和+5.2%；两点留一法传递结果则为+10.6%和?9.6%。这些均属条件性偏差，而非包含不确定性的精度界限。经验硬度桥接使未校准深度改变约11%，且现有硬度曲线无法支撑统计置信区间。主要的实验局限在于缺少实测碳浓度分布：比较仅局限于深度判据，预测的表面碳浓度仍未得到验证。对于这两种工艺制度，通量与扩散系数的灵敏度近乎共线。单次内孔观测值比另一次试验的自由表面观测值浅36%，其中载荷、淬火及实验室效应尚未厘清。标称供气计算表明流量缩减实验存在空间，但与气瓶称重结果不一致，且无法确立传递效率。实际应用限于所研究条件下的初步配方筛选及后续信息性测量的选择；更广泛的精度、浓度分布验证和几何修正仍需进一步实验。

低压渗碳(LPC)以乙炔为渗碳介质，通过短时增强脉冲在低压下引入碳，随后在真空扩散阶段使碳重新分布，这种增强–扩散循环已成为脉冲渗碳研究的成熟描述。LPC配方计算需将表面传递与碳扩散耦合，但气体可获性、表面浓度演化及反应产物会使实际碳吸收具有时间依赖性。工业应用的有效性还取决于模型评估方式：试验前预测、继承参数的非校准计算与针对同次试验的拟合回答的是不同问题，而既有研究常将数值一致性与实验验证混为一谈。为此，研究人员基于16MnCr5钢四次工业LPC试验(940–950 °C)，明确区分了试验前记录的设计输出、未校准通量模型比较与后续单参数拟合三个证据层次，围绕未拟合参数模型的渗层深度误差、历史配方预测与受控边界比较的差异、渗层深度数据能否区分候选参数、以及哪些观测超出模型适用范围四个问题展开研究。

研究人员建立了一维半隐式有限体积模型，采用浓度依赖的?gren扩散系数、预设增强碳通量、零通量扩散窗口、碳化物储库和非等温尾段，未针对本批次拟合任何参数。与两个自由表面硬度法CHD₅₅₀(HV0.3)值比较，0.36 wt.% C阈值深度的点偏差为?13.7%(Test 1)和?4.6%(Test 4)；联合通量拟合残差为?4.9%和+5.2%，两点留一法传递误差为+10.6%和?9.6%。通量与扩散系数灵敏度近乎共线，仅凭渗层深度无法可靠区分二者。研究建立了从数值验证、非校准比较到校准传递的三级证据层次，为LPC配方筛选提供了条件性定量参考，并指明后续实验应补充的观测量。

关键技术方法：一维半隐式有限体积格式；?gren浓度依赖碳扩散系数乘以Gunnarson型合金因子k_alloy（本炉次为0.868）；预设恒定碳通量增强边界与零通量扩散边界；碳化物储库与过饱和回溶算法；非等温降温段积分。数值验证包括解析互补误差函数基准、网格收敛研究（Δx=8至1 μm，深度变化0.025%）及COMSOL Multiphysics 6.4独立有限元交叉验证。试样为同一炉次(2323340V, Deutsche Edelstahlwerke)的16MnCr5+H圆柱，硬度测试由TüRKAK认可的Teknotest实验室(Tests 1、2、4自由表面)和Yeditepe大学(Test 3内孔)完成。

研究结果：
3.1 数值验证与内部灵敏度：模型通过解析基准验证，最大浓度误差6.8×10^?5 wt.%；COMSOL与有限体积代码的深度差小于0.05 μm；碳化物储库容量0.1–3.0 wt.%仅改变深度约0.5 μm；温度设定值每±10 °C引起6.1%渗层深度变化，非等温尾段对Tests 3/4贡献15 μm(1.8%)。
3.2 实测渗层深度：Test 1的CHD₅₅₀为1036 μm(表面942 HV0.3)，Test 4为907 μm(表面802 HV0.3)，芯部293–314 HV0.3；Test 3内孔为579 μm。硬度压痕采样间距限制了插值分辨率，但无法追溯标准不确定度。
3.3 非校准比较与边界条件控制：未校准模型偏差?13.7%/?4.6%；历史配方设计器输出偏差+25.4%。碳判据下通量边界优于Dirichlet浓度边界，但采用经验硬度桥接后Test 4排序反转，说明排名依赖比较判据。
3.4 校准、传递与不可分辨参数：联合通量拟合φ=1.213，残差?4.9%/+5.2%；留一法传递+10.6%/?9.6%；扩散系数拟合λ=1.576给出几乎相同深度。局部灵敏度矩阵奇异值比约3–4×103，表明这些工艺制度下参数实际病态，需实测碳浓度分布改善可辨识性。
3.5 探索性内孔观测：相同配方下Test 3内孔579 μm比Test 4自由表面907 μm浅36%，但混杂载荷、淬火气体(N? vs H?)和实验室差异；内孔单独拟合φ_bore=0.552，仅为自由表面的约45%。需同试样配对测量才能归因于几何效应。
3.6 供气状态：标称供碳量与模型吸收速率之比为49–181，但未测量实际到达表面的碳通量；气瓶称重与设定点估算存在约1.7–10倍差异，无法建立质量平衡或传递效率。
3.7 硬度：经验硬度桥接(Andrews M_s关系、Koistinen–Marburger转变及马氏体硬度曲线)在近表面低估(728 vs 802/942 HV)、芯部高估(435 vs 293–314 HV)；桥接深度d_HV550为995.5和963.7 μm，与d_C0.36相差约11%，但该转换未经验证。
3.8 欠渗碳试验：Test 2总增强时间短且含138.5 min冷却尾段，模型预测表面碳仅0.286 wt.%、无0.36 wt.%交叉；实测硬度近表面非单调(472→621→502 HV0.3)但存在约210 μm下降交叉，碳判据无法定量复现，模型仍可将其标记为渗碳不足候选配方。

讨论部分界定了四层结论边界：已建立的是数值验证模型在特定条件下的条件性非校准比较，可支持初步配方筛选；未建立的是绝对碳浓度分布验证、表面微观组织解释、芯部性能预测与几何修正。参数可辨识性应作为未来实验设计的约束，需补充实测碳浓度分布、独立吸收速率测量及不同灵敏度方向的工艺制度。工艺启示包括分步流量缩减试验、装炉面积作为缩放变量以及复杂特征的配对测量与传输模型。研究结论为：有限体积格式通过解析解、网格细化和COMSOL交叉验证，但数值一致不验证材料规律；两个自由表面未校准偏差为?13.7%/?4.6%，缺少碳浓度分布是主要实验局限，不推断通用精度区间；历史设计器偏差+25.4%因数值来源不全无法全部归因于边界条件；联合通量拟合残差?4.9%/+5.2%，两点传递+10.6%/?9.6%，经验硬度桥接使未校准深度移动约11%；通量与扩散系数难以区分，需更广泛灵敏度数据；标称供气比值为条件性模型估计，与气瓶称重不一致；36%内孔差异受多重混杂影响，φ_bore非有效几何修正；表面与芯部硬度并非经过验证的模型输出，模型对绝对碳浓度分布和复杂内腔几何仍未校验。