提出了一种考虑刀杆柔性的铣削过程模型。该模型将“刀具”系统简化为分布动态铣削力激励下在其静平衡位置附近发生双向横振动的柔性悬臂梁。基于该模型对刀杆的动态变形进行数值仿真，在这之前，需要获得刀杆的各阶固有频率值。然而，由于螺旋棒铣刀的截面形状不规则（见图1），无法用公式直接计算截面惯性矩。若用公式：
图1 圆柱螺旋铣刀横截面示意图
(1)
式中：Cq——当量系数；
D——铣刀直径。
求得铣刀截面的惯性矩，则当量系数的选取存在很大的主观性。运用GREEN公式将区域积分化成边界积分，试图通过数值方法确定刀杆的截面惯性矩，然而该方法的精度对所划分的边界单元的个数十分敏感，而且单元的划分没有一定的规则可循，人为因素显得过分突出。况且，即便能够获得准确的截面惯性矩，而材料的弹性模量也只能在一个大致范围内选取，因此，不可能借助公式(2)获得准确的刀杆各阶固有频率。
（2）
式中：ωn——铣刀实际工况下的第n阶固有频率；
l——铣刀工作时的实际悬臂长度；
knl——频率方程
ρ——铣刀材料单位体积质量；
A——铣刀有效横截面积；
E——铣刀材料的弹性模量；
I——铣刀截面的有效惯性矩。
表1 频率方程的前5个根
为了解决这一问题，本文采用实验的方法先测得刀杆的有效截面系数
2 有效截面系数测定原理
图2所示为测量装置示意图。测量时，用橡皮槌敲击紧固在三爪卡盘上的被测铣刀，响应信号通过加速度传感器，经电荷放大器放大，再经A/D转换，进入计算机。通过模态分析软件包作FFT变换，即可获得响应频谱。若测得铣刀的基本固有频率为ω1，测得的悬臂长度为l，那么其有效截面系数：
1.被测铣刀；2.加速度传感器；3.有机玻璃套；4.激振用橡皮槌。
图2 铣刀有效截面系数测定原理示意
(3)
实际测定时，可取其多次测量计算的平均值。
有效截面系数ζ测定之后，将其代入公式(2)，便可获得铣刀比较精确的各阶固有频率值。
3 应用实例
3.1 有效截面系数的测定
为了表明该方法的有效性，现以8的直柄平头圆柱螺旋铣刀为例，进行实测计算。该铣刀的相关参数见表2。
表2 铣刀及相关参数
表3 有效截面系数
(mm)
率实测值
(Hz)
系数ζ
频率计算
值(Hz)
频率实测
值(Hz)
3.2 刀杆变形仿真
现以实际切削实验数据为依据，提出的铣削过程模型，对立铣过程中刀杆产生的总体变形量进行仿真。已知：仿真时间间隔Δt＝0.02 ms，仿真时间t=130 ms，铣削过程参数见表4。
表4 铣削过程参数
图3 刀杆总体位移仿真曲线
从仿真结果可以看出，刀杆的动力响应是在静态变形之上迭加一个高频动态分量而成。这说明铣刀的总体变形由动、静两部分组成，且在一般稳态切削条件下，静变形占主要地位，从宏观上将导致一定量的几何形状误差和波纹度。
4 结论
本文通过引入刀杆有效截面系数的概念，提出了一种在不同悬臂长度下确定刀具各阶固有频率的新方法。对于同一把铣刀，经一次性测量计算得到其有效截面系数之后，就可用它确定该铣刀在不同悬臂长度时比较精确的各阶固有频率值，无需再去关心刀具材料特性参数和几何参数。
