胡 腾，王 虎，阳 红，郑华林，刘 雁

(1.西南石油大学 机电工程学院，四川 成都 610500；2.中国工程物理研究院 机械制造工艺研究所，四川 绵阳 621999；3.四川普什宁江机床有限公司，四川 都江堰 611830)

目前，颤振已成为高效铣削领域亟待解决的关键共性问题[1]。准确预报铣削稳定性并将其作为约束条件，可为铣削工艺参数优化提供坚实的理论基础和数据支撑，有利于实现高效、无颤振铣削加工，亦有利于延长刀具及机床寿命，具有较重要的科学意义与工程价值。

近年来，国内外学者针对铣削稳定性预报方法开展了较深入的研究，取得了不少成果：Ramirez等[2]提出了基于Lyapunov-Krasovskii分析的铣削稳定性预报方法，揭示了稳定性与切削参数之间的复杂联系；Sekar等[3]联合自适应步长R-K方法与误差控制算法，提出了铣削稳定域高精度仿真方法，该方法亦可推广至多时滞动力学特征值问题的求解；Li等[4]提出了一种直接利用加工系统模态参数或频响函数快速绘制稳定性叶瓣图的方法，较显著地提升了切削稳定性预报效率；Ahmadi等[5]分别基于多频率法和半离散法对铣削稳定性进行了预报，并借助等效粘性模型探明了过程阻尼对系统低转速区切削颤振的抑制作用；Li等[6]建立了考虑模态耦合和过程阻尼的铣削动力学模型，并利用二阶半离散法对铣削稳定性进行了预报；Ezugwu等[7]提出了基于最小二乘近似的全离散法，并以此分别对再生型极限轴向切深和极限径向切宽进行了分析与预报。王民等[8]从再生型颤振机理出发，基于零阶频域法构建了3维稳定性叶瓣图，系统探讨了主轴转速、轴向切深与径向切宽之间复杂的非线性映射关系；……