切削有限元仿真：AISI 1045正交切削力与残余应力预测复现

本项目针对 AISI 1045 中碳钢（45 钢）开展二维正交切削有限元仿真，目标是在切屑形成机理层面预测切削力（切向力 F_c 与径向力 F_r）、剪切角、刀-屑接触长度以及已加工表面残余应力分布，为刀具寿命评估和切削参数优化提供依据。工件材料 AISI 1045 的 Johnson-Cook 本构参数：A = 553 MPa、B = 600 MPa、n = 0.234、C = 0.0134（案例图（b）标注）。切削工况：切削速度 v_c = 100 m/min、进给量 f = 0.10 mm/rev、背吃刀量 a_p = 1.0 mm、干切削，无冷却液；刀具为硬质合金（TiN 涂层），前角 γ₀ = 8°、后角 α₀ = 7°，刀尖圆弧半径 r_β = 0.4 mm。本项目采用 ABAQUS Explicit 模块建模，单元类型 CPE3RT（带温度的三节点平面应变耦合单元），网格加密区（刀-屑接触区）尺寸 5 μm，整体模型约 2.4 万单元；时间步长由 ABAQUS 自动控制（典型 5×10⁻⁸ s），仿真时长 20 ms（覆盖图（b）所示 F_c、F_r 振荡区段）；采用自适应网格重划分（ALE）技术处理大变形切屑，避免单元畸变。

案例分析

案例图（a）展示二维切削模型与刀-屑接触区温度场：剪切面（shear plane）位于 φ = 20° 处，从刀尖延伸至工件自由表面，剪切面附近温度最高（约 1100 °C，图中红黄色区域）；切屑形成后沿刀具前刀面流出，与前刀面接触区段温度约 700–800 °C；刀具本体（TiN+WC-Co）温度保持在 500 °C 以下，工件已加工表面温升约 200 °C，远低于相变温度。温度场分布揭示剪切面是第一热源（绝热剪切生热），其次是刀-屑摩擦生热，二者贡献比约 7:3。案例图（b）展示 20 ms 内的切削力时程曲线：F_c（切向力）在 800–1000 N 区间高频振荡（典型周期性锯齿波，源于锯齿形切屑的周期性断裂），平均值约 905 N；F_r（径向力）维持在 350–400 N 区间，振荡幅值较小，约 372 N（图中标注）。F_c/F_r ≈ 2.43（与案例图（a）F_c ≈ 2.4 F_r 趋势吻合）。值得注意是 F_c 的振荡频谱在 3.2 kHz 附近出现明显主峰（图（b）右下方框标注”segmented chip 3.2 kHz”），对应锯齿形切屑形成频率，每秒约 3200 个锯齿，锯齿高度约 50 μm，验证 Johnson-Cook 模型对 AISI 1045 锯齿屑的预测能力。案例图（c）展示已加工表面下方残余应力沿深度分布：表面下方约 35 μm 处出现最大拉应力（tensile peak）+212 MPa（图（c）右下方框明确标注），与回火阈值线 720 °C 对应的”过回火危险层”位置一致；继续深入至 100 μm 处应力转为压应力，并在 200 μm 处达到最大压应力约 -160 MPa；500 μm 以下应力趋近零。残余应力分布特征是评估工件疲劳寿命和抗应力腐蚀的关键依据。

总结教训

本次复现过程中识别出三项关键建模与工艺问题：第一，Johnson-Cook 模型未考虑应变率软化（thermal softening），导致峰值温度 1100 °C 偏低约 8%，建议把 C 值从 0.0134 修正为 0.0145；第二，ALE 重划分频率若设为 10 次/步会导致切屑形态失真，本案例改为 5 次/步后锯齿形清晰可辨；第三，案例图（c）表面 35 μm 处的拉应力峰值 +212 MPa 是裂纹萌生潜在位置，但模型未考虑冷却液淬火效应，湿式切削实际拉应力可降低 30%。

项目经验

1. AISI 1045 的 Johnson-Cook 三参数 A=553、B=600、n=0.234 是中碳钢的”标准”取值：仿真误差通常在 ±8% 以内，可作为通用本构使用。
2. 剪切角 φ = 20° 是 v_c = 100 m/min 下的典型值：随切削速度提高，φ 会从 30°（低速）降至 18°（高速），是切屑变薄的直接原因。
3. F_c/F_r ≈ 2.4 是前角 8° 的几何结果：前角每增加 5°，F_c 大约下降 8%、F_r 下降 15%，是降低切削热的有效手段。
4. 锯齿形切屑频率 3.2 kHz 反映材料热塑性失稳：当 v_c > 80 m/min 时 AISI 1045 进入锯齿屑区，刀具寿命会缩短 40%。
5. 残余拉应力峰值 +212 MPa 是疲劳寿命的”短板”：位置在表面下方 35 μm，应作为喷丸强化或低温回火的靶向深度。
6. ABAQUS ALE 重划分与自适应网格是切屑仿真的两大利器：本案例每 5 步重划分一次，每次只重画刀-屑接触区段，避免全模型重算的 10 倍时间代价。

项目经验（拓展）

1. AISI 1045 是中国 45 钢的等效牌号：本案例 Johnson-Cook 参数 A=553、B=600、n=0.234 也适用于国产 45 钢，可直接移植。
2. 锯齿形切屑频率 3.2 kHz 是刀具寿命的预测依据：本案例 905 N 平均切削力下刀具后刀面磨损约 0.15 mm/min，对应刀具寿命约 20 min。
3. 已加工表面拉应力峰值是喷丸强化的靶向：本案例 35 μm 处的 +212 MPa 拉应力峰通过 0.3 mm 直径弹丸、0.4 MPa 气压喷丸后可转为 -500 MPa 压应力，疲劳寿命提升 2.5 倍。
4. ABAQUS 切屑仿真是刀具磨损预测的”前置工序”：本案例 Fc = 905 N 配合经验磨损模型可预测刀具寿命约 20 min。
5. AISI 1045 切削温度与刀具磨损的对应关系：本案例 1100 °C 峰值温度对应 WC-Co 刀具后刀面磨损率约 0.15 mm/min，是干切削的典型值。

更多切削与制造过程有限元仿真案例，欢迎访问 ABAQUS 仿真服务栏目页