基于CALPHAD与析出相表征的Inconel 718合金物理本构模型及其在钻削切屑形成有限元模拟中的应用

《JOURNAL OF MATERIALS PROCESSING TECHNOLOGY》：A physics-based flow stress model for chip formation simulation when machining Alloy 718

　　

摘要：研究人员开发并验证了基于微观组织强化机制（位错强化、Hall–Petch晶界强化、固溶强化及γ′/γ″析出相强化）的物理基流动应力(Flow Stress)模型，用于描述时效态Inconel 718合金在宽温度与应变率范围内的塑性变形行为，并将其嵌入三维钻削有限元(FEM)仿真中以预测热机械载荷与切屑形态。该模型通过CALPHAD（Computer Coupling of Phase Diagrams and Thermochemistry）方法计算γ′与γ″相的尺寸及体积分数，并结合Split-Hopkinson Pressure Bar(SHPB)高温高应变率实验进行参数标定。二维正交切削与全三维钻削仿真结果表明，相较于经验性的Johnson–Cook(J–C)本构方程，所提出的物理基本构模型能更准确地预测不同晶粒尺寸(d = 9 μm、15 μm、33 μm)材料条件下的钻削扭矩与刀–屑接触长度，且在考虑微观组织演变（晶粒尺寸、析出相分数）时无需重新进行逆标定。摩擦条件通过类钻削工况下的专用摩擦实验确定Coulomb摩擦系数μ = 0.4。研究表明，物理基模型因引入动态回复(Dynamic Recovery)与温度相关的析出相溶解修正，能更合理地描述切削区的高应变率热软化行为，为镍基高温合金加工过程数值仿真提供了更具通用性的本构描述框架。

论文解读——基于CALPHAD与微观组织强化的Inconel 718物理基本构模型及其在钻削FEM仿真中的应用

---

一、研究背景与意义

Inconel 718（Alloy 718）等镍基高温合金因优异的高温强度广泛用于航空航天件，其切削加工伴随极高应力与复杂切屑形态，数值仿真成工艺优化的重要工具。传统广泛使用的Johnson–Cook(J–C)本构模型属唯象(Phenomenological)模型，无法显式纳入晶粒尺寸、γ′（Ni₃(Al,Ti)）与γ″（Ni₃Nb）析出相尺寸/体积分数等微观组织参数，且需针对不同热处理/晶粒条件反复逆标定，在高温与大应变下易过预测硬化。为克服上述局限，研究人员提出一种基于位错密度演化与多种强化机制叠加的物理基(Physics-based)流动应力模型，借助CALPHAD计算时效析出相信息，经SHPB实验标定后分别用于二维正交切削与三维钻削仿真，并与J–C模型对照验证。该文发表于《Journal of Materials Processing Technology》。

---

二、主要关键技术方法

① 材料表征：取三种显微组织状态(MC1–MC3)的Inconel 718（锻/热轧，AMS 5663双级时效），测定平均晶粒尺寸d（9 μm、15 μm、33 μm）及硬度；

② CALPHAD热力学/动力学计算：采用ThermoCalc?配TCNI8与MOBNI4数据库及TC-PRISMA?模块（KWN模型——Kampmann–Wagner Numerical model），估算固溶处理后基体溶质浓度及γ′（球状）、γ″（盘状，纵横比6）相在720℃×8h+620℃×8h双级时效后的半径r、半厚t与体积分数f_v，并按Semiatin伪平衡修正高温回溶；

③ 本构建模：流动应力σ = σ₀（点阵阻力Peierls–Nabarro）+ σ_s（多组元固溶强化Labusch–Gypen模型）+ σ_H–P（Hall–Petch晶界强化）+ σ_P（γ′/γ″共格应变强化Gerold–Haberkorn/Oblak模型）+ σ_G（Taylor位错林强化），位错密度ρ_i由存储率（与d、位错胞尺寸、析出间距相关）与热激活动态回复率控制，采用径向返回算法(Radial Return Algorithm)与Newton–Raphson迭代更新；待拟合参数由SHPB多温-多应变率流变曲线约束最小二乘标定；

④ SHPB（Split-Hopkinson Pressure Bar）实验：感应加热至20~750℃，冲击气压1.5~2.5 bar，获真应力–真应变数据；

⑤ 摩擦标定：Berger PFS 5558-1正交摩擦实验台模拟刀–屑二次剪切区工况（TiN/TiAlN涂层对摩Inconel 718，干摩擦，相对速度≈15 m/min），测定Coulomb摩擦系数μ ≈ 0.4；

⑥ 有限元仿真：DEFORM? 2D正交切削（vc = 30、60 m/min）与3D钻削（vc = 25 m/min，f = 0.1 mm/rev，实心硬质合金TiN/TiAlN麻花钻?8 mm），J–C模型取自Malakizadi等文献作参照，物理基模型按d = 9 μm与d = 33 μm生成列表化流动应力数据，热接触HTC = 1×10⁵kW/(m²K)；

⑦ 实验验证：KISTLER 9170A动态扭矩仪测钻削扭矩，高速摄像记录透明钻孔通道内切屑形成，金相测切屑厚度。

---

三、研究结果

3.1 材料表征与CALPHAD结果

SEM确认MC1(d≈33 μm)、MC2(d≈15 μm)、MC3(d≈9 μm，热轧)。EDS鉴别晶界NbC型碳化物与针状δ(Ni₃Nb)相。固溶后γ基体中残余溶质主要为Cr、Fe、Mo；经TC-PRISMA?计算γ′半径r_γ′≈6.1 nm、f_v,γ′≈4%，γ″等效半径r_γ″≈21.4 nm、f_v,γ″≈14.6%，与文献峰值时效相符；高温下按Solvus温度(T_s,γ′≈1118 K、T_s,γ″≈1200 K)及Semiatin缩放因子修正体积分数与尺寸。

3.2 物理基本构模型标定

SHPB流变曲线拟合得标定参数（表5）：τ₀=6.5×10^?3、Δf₀=0.43、q=2.60、p=0.37、Z=5×10^?3、Ω=2.0、Φ=1.5、k_grain=0.001、k_cell=190、k_prec=2×10^?4、ω=6.5。模型较好复现不同T、?下的屈服与加工硬化/动态回复平台。

3.3 摩擦特性实验

干摩擦实验得平均法向力F_n,avg=1709 N、摩擦力F_r,avg=676 N，对应σ_N=907 MPa、τ_f=359 MPa，Coulomb摩擦系数μ_avg=0.4（v_r≈14.94 m/min），适用于后续仿真。

3.4 二维正交切削仿真结果(v_c= 30、60 m/min)

J–C模型(d≈16 μm)预测剪切角Φ较大(30 m/min时Φ_JC=26.9°，60 m/min时Φ_JC=32.2°)，导致接触长度偏短(误差Δ_exp,JC16=?44.6%@30 m/min)、最大等效应力偏高(σ_eff,max≈1990 MPa)；物理基模型(d=9/33 μm)给出较小剪切角(Φ_PB≈21.7°~24.4°)，接触长度预测误差降至?28.9%(30 m/min)与?21.5%(60 m/min)，最大等效应力降低至1780~1850 MPa。切削力预测：J–C偏差?4.4%(30 m/min)、?3.6%(60 m/min)；物理基偏差+3.9%~+5.5%(30 m/min)、+5.8%~+7.4%(60 m/min)，均具较高精度。切屑厚度J–C预测较准(误差8.5%@30 m/min)，物理基因较小剪切角致切屑厚度被高估(≈50%)，此差异源于二维简化与摩擦模型近似。Coulomb μ=0.4与压力相关剪切摩擦模型对比显示二者对主切削力预测相当，后者略改善切向力预测。

3.5 三维钻削仿真结果与实验对比

• •
  钻削扭矩：实验MC1(d=33 μm)：M_z,exp=3.65 Nm；MC3(d=9 μm)：M_z,exp=3.69 Nm。J–C仿真偏差ΔM_z,JC16=+5.38%；物理基模型偏差ΔM_z,PB33=?1.86%(d=33 μm)、ΔM_z,PB9=?0.88%(d=9 μm)，显著优于J–C模型。
• •
  切屑流速验证：二次剪切区平均切屑流速v_ssz,JC≈13.86 m/min、v_ssz,PB9≈14.07 m/min，与摩擦实验设定v_r=14.94 m/min吻合，佐证μ=0.4合理性。
• •
  切屑形态：高速摄像显示典型向内弯曲→贴前刀面卷曲→碰容屑槽折断过程；仿真定性再现卷曲形态，但数值曲率角略大于实测（实验κ_exp=54°，J–C κ=62°，物理基 κ=60°），表明物理基仍轻微高估材料流动抗力，但不影响扭矩与接触长度精度优势。

---

四、讨论与结论总结

讨论指出：物理基模型因显式含热软化、应变率敏感性与动态回复项，流变应力在大应变趋于稳态，较J–C无界硬化更贴近切削实际；其可随CALPHAD输入的析出相信息与实测晶粒尺寸自适应调整，无需如J–C反复逆标定，对不同热处理/工艺状态具更好迁移性。当前局限含SHPB应变/应变率低于真实切削区、二维切屑厚度预测偏差及未引入绝热剪切带(DASB)与损伤模型，未来建议发展SHPB数据+正交切削逆标定混合方法，并细化速度/温度依赖摩擦模型。

结论翻译：

① 建立了含晶界(Hall–Petch)、固溶(Cr/Fe/Mo)、γ′/γ″共格应变强化及位错存储–动态回复机制的Inconel 718物理基流动应力模型，结合CALPHAD预测特定成分与热处理后析出相尺寸/分数，可在宽T–?范围内描述流动应力且随微观组织自适应。

② 物理基模型只需已知化学成分与热处理制度即可建立，免去唯象模型对不同显微组织的重复逆标定，降低实验成本。

③ 实验标定的恒定视在Coulomb摩擦系数(μ=0.4)在钻削运动学范围内足以合理描述刀–屑界面复杂摩擦，与压力相关模型相比未显著降低预测质量。

④ 2D正交切削中物理基模型显著改善接触长度预测(误差由?44.6%降至?28.9%)，切削力误差相当。

⑤ 3D钻削扭矩预测：J–C偏差5.38%，物理基偏差最低至?0.88%(d=9 μm)与?1.86%(d=33 μm)，提升明显。

⑥ 该物理基本构框架为不同热处理态镍基高温合金切削仿真提供可靠基础，可支撑航空制造中刀具热机械载荷与工件微观组织演变预测。