铣削受力图与受力分析：从原理到代码的切削力计算指南

简介：这份资源围绕铣削加工中工件受力这一核心问题展开，面向金属切削工艺学习者、机械加工技术人员及机械专业学生，帮助理解铣削力、切削力与摩擦力的合成关系，并对比逆铣与顺铣两种方式的受力差异与适用场景。压缩包内共1个docx文档，约3.51MB，以图文结合的方式呈现受力图与受力分析，便于对照阅读与整理笔记。资源中配有铣削时工件受力图，明确标注合力FR、竖直分力FV与水平分力FH的方向关系，并说明逆铣时FV向上、振动较大，顺铣时FV向下、切削平稳，一般粗加工用逆铣、精加工用顺铣。目前已有3501人学习下载，适合需要快速掌握铣削受力规律、为工艺选择与参数优化提供理论依据的读者参考。

1. 铣削时工件受力图及其受力分析：从一张图看懂切削力怎么算、怎么用

车间里经常出现一种争论：机床功率够不够、夹具夹得紧不紧、刀具为什么崩刃，大家各说各话，最后往往靠“手感”拍板。其实这些问题背后都指向同一张图——铣削时工件受力图。它把刀齿作用在工件上的力，按空间方向拆成切向力、径向力、轴向力，再合成到进给、法向、轴向三个方向，配合受力分析就能算出大小和方向。搞懂这张图，你能判断夹具该往哪顶、工艺系统刚性够不够、参数该往哪调；做工艺、夹具、刀具、仿真的从业者都用得上。下面按“图怎么画、力怎么算、坑在哪”一步步拆开讲。

2. 铣削受力图怎么画：坐标系、分力与合成路径

画受力图不是随手画几个箭头，核心是先定坐标系，再定各分力的正方向，最后做矢量合成。坐标系选错，后面所有计算都是白费。

2.1 先定坐标系：工件坐标系与刀具坐标系的区别

铣削受力分析里有两套坐标系，混用是新手最常见的翻车点。

工件坐标系（通常记作 Xw-Yw-Zw）固定在工件上，Xw 沿进给方向，Yw 沿横向（法向），Zw 沿工件表面法线（一般朝上）。这套坐标系的好处是：夹具设计、工件变形、振动分析都在这套系里做，因为夹具和工件是“死”的。

刀具坐标系（记作 Xc-Yc-Zc）跟着刀齿走，原点在刀尖接触点，Xc 沿切削速度方向（切向），Yc 沿径向（指向刀轴），Zc 沿刀轴方向。切削力的原始定义就在这套系里：切向力 Ft、径向力 Fr、轴向力 Fa。

两套坐标系之间差一个旋转矩阵，旋转角就是刀齿的瞬时位置角 φ。对圆柱铣刀（周铣），φ 从切入角扫到切出角；对面铣刀，φ 是刀齿在端面上的方位角。常见做法是：先在刀具坐标系里算 Ft、Fr、Fa，再通过旋转矩阵投影到工件坐标系，得到 Fx、Fy、Fz。这样夹具和变形分析就能直接用。

提示：如果你只做刀具寿命评估，用刀具坐标系就够了；只要涉及夹具、工件变形、颤振，必须转到工件坐标系。

2.2 三个分力的物理含义与正方向约定

铣削中每个刀齿受到的力可以拆成三个正交分量：

正方向约定必须统一：一般取 Ft 与切削速度同向为正，Fr 指向刀轴为正，Fa 沿刀轴指向工件为正（或按右手定则统一）。约定不统一，后面合成时符号会乱，算出来的合力方向可能完全反。

对螺旋刃立铣刀，螺旋角 β 会把一部分切向力转到轴向，所以 Fa 不可忽略。β 越大，Fa 越大，Ft 相对减小。这就是为什么大螺旋角刀具切削轻快但轴向力大，主轴轴承和夹具有时要为此买单。

2.3 从单齿到多齿：受力图的合成逻辑

实际铣削是多齿同时参与，受力图要按“单齿 → 多齿叠加 → 整周期”三步走。

第一步，算单齿在位置角 φ 处的瞬时力。用切削力系数模型：

import numpy as np # 切削力系数（单位：N/mm^2），需按工件-刀具组合标定 Kt = 1800.0 # 切向系数 Kr = 0.35 # 径向与切向之比 Ka = 0.25 # 轴向与切向之比 def single_tooth_force(phi, ap, fz, Kt=Kt, Kr=Kr, Ka=Ka): """ phi: 刀齿瞬时位置角（弧度），0 表示切入 ap: 轴向切深（mm） fz: 每齿进给（mm/z） 返回：单齿在刀具坐标系下的 Ft, Fr, Fa """ # 瞬时切屑厚度 h = fz * sin(phi)，仅当 sin(phi) > 0 时切削 h = fz * np.sin(phi) if h <= 0: return 0.0, 0.0, 0.0 # 微元切削力 = 系数 × 切屑厚度 × 切深 Ft = Kt * h * ap Fr = Kr * Ft Fa = Ka * Ft return Ft, Fr, Fa

这段代码的关键在h = fz * sin(phi)：切屑厚度随位置角变化，切入时为零，中间最大，切出时又归零。Kt、Kr、Ka不是常数，跟材料、刀具涂层、刃口钝圆都有关，必须用实验标定，不能照抄手册。ap是轴向切深，对周铣是侧刃参与长度，对面铣是端刃切深。

第二步，把单齿力投影到工件坐标系：

def tooth_to_workpiece(Ft, Fr, Fa, phi): """把刀具坐标系的力投影到工件坐标系（周铣简化模型）""" # 切向沿速度方向，径向指向刀轴，投影到 Xw-Yw 平面 Fx = Ft * np.cos(phi) + Fr * np.sin(phi) Fy = -Ft * np.sin(phi) + Fr * np.cos(phi) Fz = Fa return Fx, Fy, Fz

第三步，对同时参与的所有刀齿求和，再扫一个刀齿周期，得到 Fx、Fy、Fz 随时间的波形。这个波形就是受力图的“动态版”，平均值用于功率和夹具设计，峰值用于刀具强度和颤振判断。

2.4 用一段最小代码把受力图画出来

把上面两步串起来，扫一个周期就能出图：

import numpy as np import matplotlib.pyplot as plt N = 4 # 刀齿数 phi_p = 2 * np.pi / N # 齿间角 ap, fz = 2.0, 0.1 # 切深 2mm，每齿进给 0.1mm phis = np.linspace(0, 2 * np.pi, 720) Fx_t, Fy_t, Fz_t = [], [], [] for phi in phis: Fx_sum = Fy_sum = Fz_sum = 0.0 for k in range(N): phi_k = (phi + k * phi_p) % (2 * np.pi) Ft, Fr, Fa = single_tooth_force(phi_k, ap, fz) Fx, Fy, Fz = tooth_to_workpiece(Ft, Fr, Fa, phi_k) Fx_sum += Fx; Fy_sum += Fy; Fz_sum += Fz Fx_t.append(Fx_sum); Fy_t.append(Fy_sum); Fz_t.append(Fz_sum) plt.plot(phis, Fx_t, label='Fx') plt.plot(phis, Fy_t, label='Fy') plt.plot(phis, Fz_t, label='Fz') plt.xlabel('刀具转角 (rad)'); plt.ylabel('力 (N)') plt.legend(); plt.title('铣削力周期波形') plt.show()

逻辑说明：外层扫刀具转角，内层对每个刀齿判断是否在切削区，再叠加。参数说明：N改刀齿数，ap改切深，fz改每齿进给，Kt/Kr/Ka改材料组合。跑出来的波形能直接看出峰值和均值，均值乘线速度就是功率，峰值用于判断冲击。

3. 受力分析怎么算：从经验公式到系数标定

画完图只是第一步，真正要落地得把力算准。算不准，夹具设计就是拍脑袋。

3.1 经验公式的适用边界与误差来源

车间里最常用的是经验公式，比如 F = Kc × ap × fz × z / (π × D) 这类形式。它的好处是快，坏处是误差大，通常 ±30% 起步。误差来源主要有四个：

一是切削力系数 Kc 随材料硬度、热处理状态波动，同一牌号不同批次能差 15%。二是刃口钝圆半径，新刀和磨损刀的 Kc 能差 20% 以上。三是切削液，干切和湿切的摩擦分量不同。四是刀具跳动，实际每齿进给不均匀，公式假设的是理想均匀。

所以经验公式只适合粗估，比如判断机床功率够不够、夹具大概要多大夹紧力。要做精细分析，必须标定。

3.2 用测力仪标定 Kt、Kr、Ka 的最小流程

标定不复杂，但步骤要严。常见做法是：用测力仪（三向 dynamometer）测 Fx、Fy、Fz，同时记录刀具转角或主轴编码器信号，然后反算系数。

import numpy as np from scipy.optimize import least_squares # 实测数据：phi 为转角，Fx_m/Fy_m/Fz_m 为测得的力 phi = np.loadtxt('phi.txt') Fx_m = np.loadtxt('Fx.txt') Fy_m = np.loadtxt('Fy.txt') Fz_m = np.loadtxt('Fz.txt') def model(params, phi, ap, fz, N): Kt, Kr, Ka = params Fx = np.zeros_like(phi); Fy = np.zeros_like(phi); Fz = np.zeros_like(phi) phi_p = 2 * np.pi / N for i, p in enumerate(phi): for k in range(N): pk = (p + k * phi_p) % (2 * np.pi) h = fz * np.sin(pk) if h <= 0: continue Ft = Kt * h * ap Fr = Kr * Ft; Fa = Ka * Ft Fx[i] += Ft * np.cos(pk) + Fr * np.sin(pk) Fy[i] += -Ft * np.sin(pk) + Fr * np.cos(pk) Fz[i] += Fa return np.concatenate([Fx, Fy, Fz]) def residual(params): pred = model(params, phi, ap=2.0, fz=0.1, N=4) meas = np.concatenate([Fx_m, Fy_m, Fz_m]) return pred - meas res = least_squares(residual, x0=[1800, 0.35, 0.25], bounds=([500, 0.1, 0.0], [4000, 0.8, 0.6])) print('标定结果 Kt, Kr, Ka =', res.x)

逻辑说明：用最小二乘拟合，让模型预测的力和实测力残差最小。参数说明：x0是初值，bounds是物理合理范围，防止拟合出负系数。ap、fz、N必须和实验条件一致，否则标定无效。

标定时要注意：至少做三组不同 fz，覆盖实际加工范围；每组重复三次取平均；测力仪采样率要高于刀齿通过频率的 5 倍以上，否则峰值会被削掉。

3.3 把受力分析接到夹具与功率校核上

算出力之后，两个最直接的用途是夹具夹紧力校核和机床功率校核。

夹具校核：把 Fx、Fy、Fz 的峰值作为外载荷，检查工件是否滑动或翻转。滑动判据是 F_friction ≥ sqrt(Fx² + Fy²)，其中 F_friction = μ × N_clamp，μ 一般取 0.1～0.15（钢-钢干摩擦）。翻转判据是力矩平衡，要算合力作用点到支撑点的力臂。

功率校核：切削功率 Pc = Ft_avg × vc / 60000（kW），vc 是切削速度（m/min）。再除以机床效率 η（一般 0.7～0.85），得到主轴所需功率。如果超过机床额定功率的 80%，就要降参数。

注意：峰值力用于夹具和刀具强度，平均力用于功率和热。两者不能混用，混用是常见错误。

4. 避坑与排查：受力分析里最容易翻车的五件事

受力分析看起来是套公式，实际翻车点很多。下面五条是血泪经验，每条按现象、原因、解决写。

4.1 现象：算出来的力方向反了，夹具设计全错

原因：坐标系正方向约定不统一，或者旋转矩阵符号写反。最常见的是把 Fr 的投影符号搞错，导致 Fy 整体反号。

解决：先画一张示意图，标清每个分力的正方向，再用一个已知工况验证。比如单齿切入瞬间，Ft 应该和切削速度同向，Fr 应该指向刀轴。验证通过再往下算。

4.2 现象：仿真波形和实测对不上，峰值差一倍

原因：刀具跳动没考虑。实际刀具跳动会让每齿进给不均匀，有的齿切得多，有的齿切得少，峰值可能差 50% 以上。

解决：在模型里加跳动量。常见做法是给每个刀齿一个径向偏移 δk，修正切屑厚度 h = fz × sin(φ) + δk。δk 用百分表测，一般 5～20 μm。

4.3 现象：标定出来的 Kr 是负数，物理上不可能

原因：实测数据里混入了非切削力，比如测力仪零点漂移、工件自重、装夹预紧力。这些力被拟合进去，系数就乱了。

解决：标定前先空转采一段，做零点扣除；装夹后采一段静态，扣除静态分量。拟合时加 bounds 约束，防止出负值。

4.4 现象：功率校核通过，实际加工却闷车

原因：只算了平均力，没算峰值。峰值力可能是平均力的 2～3 倍，瞬时功率超限。

解决：用峰值力算瞬时功率，或者用平均力但留 1.5 倍安全系数。更稳的做法是看主轴电流波形，峰值电流对应的功率才是真实需求。

4.5 现象：换了一批材料，原来的系数完全不能用

原因：切削力系数对材料状态极敏感。同一牌号，退火和调质的 Kt 能差 30%。

解决：每换一批材料，至少做一组快速标定。如果不想每次都标，就按材料硬度做插值，但误差会到 20% 左右，只适合粗估。

5. 进阶技巧：用受力图反推刀具与参数优化方向

受力图不只是算力，它还能反过来指导刀具选型和参数优化。下面三个技巧是我平时用得最多的。

5.1 用 Fx/Fy 比值判断刀具是否偏磨

正常切削时，Fx 和 Fy 的比值在一个稳定区间。如果某个刀齿的 Fx/Fy 突然变大，说明该齿径向力异常，通常是偏磨或崩刃的前兆。做法很简单：把每齿的 Fx、Fy 单独提取出来，算比值，画趋势图。比值漂移超过 15% 就停机检查。

# 按齿提取 Fx/Fy 比值 ratio = [] for k in range(N): Fx_k = [fx for i, fx in enumerate(Fx_t) if i % N == k] Fy_k = [fy for i, fy in enumerate(Fy_t) if i % N == k] ratio.append(np.mean(np.abs(Fx_k)) / np.mean(np.abs(Fy_k))) print('各齿 Fx/Fy 比值:', ratio)

参数说明：N是刀齿数，Fx_t、Fy_t是整周期波形。比值异常的那个齿就是要检查的对象。

5.2 用轴向力 Fa 判断螺旋角是否匹配

Fa 太大，主轴轴承负载高，工件容易轴向窜动；Fa 太小，排屑可能不畅。一般立铣刀螺旋角 30°～45°，Fa 约占 Ft 的 20%～35%。如果实测 Fa 超过 Ft 的 50%，说明螺旋角偏大或轴向切深过大，可以考虑换小螺旋角刀具或降 ap。

5.3 用受力波形判断颤振风险

颤振的根源是切削力波形和工艺系统固有频率耦合。如果波形里出现明显的高频分量，且频率接近系统固有频率，就要警惕。做法是对 Fx 波形做 FFT，看主频和系统模态频率是否接近。

from numpy.fft import fft, fftfreq Fx_arr = np.array(Fx_t) n = len(Fx_arr) yf = np.abs(fft(Fx_arr))[:n // 2] xf = fftfreq(n, d=1.0 / 720)[:n // 2] # 采样间隔按实际转角步长换算 peak_idx = np.argmax(yf[1:]) + 1 print('主频对应转角频率:', xf[peak_idx])

参数说明：d是采样间隔，要按实际主轴转速换算成时间。主频接近固有频率时，优先降 ap 或改齿数，改转速有时反而更糟。

我自己的习惯是：每上一个新工况，先跑一遍受力图，看峰值、均值、波形主频三个指标，再决定参数。这套流程帮我省过好几次崩刃和闷车。希望帮到你。

本文还有配套的精品资源，点击获取