基于切比雪夫分割法的铣削稳定性分析-李鹤.pdf

《基于切比雪夫分割法的铣削稳定性分析-李鹤.pdf》由会员分享，可在线阅读，更多相关《基于切比雪夫分割法的铣削稳定性分析-李鹤.pdf（5页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、收稿日期: 2015 - 02 - 16基金项目:国家自然科学基金资助项目( 51175071) ;中央高校基本科研业务费专项资金资助项目( N120203001, N130803001) .作者简介:李 鹤( 1975 - ) ,男,河南方城人,东北大学教授,博士生导师.第37卷第4期20 1 6年4月东北大学学报(自然科学版)Journal ofNortheastern University(Natural Science)Vol. 37,No. 4Apr. 2 0 1 6doi: 10. 3969/j. issn. 1005 - 3026. 2016. 04. 018基于切比雪夫分割法的

2、铣削稳定性分析李 鹤,韩 萍,张玲利(东北大学机械工程与自动化学院,辽宁沈阳 110819)摘 要:利用切比雪夫分割法离散延时微分方程中的时间周期项,然后采用二阶Newton差商方法求出每个时间元端点处铣刀的运动向量,进而得到铣削系统的传递矩阵,并用Eloquet理论判断铣削系统的稳定性.利用以上方法,研究了一个两自由度铣削系统稳定性问题,得到了系统稳定性叶瓣图.结果表明:在同等计算精度的要求下,切比雪夫分割法比平均分割法所取的点数要少,计算效率较高;将时间周期分割成相同的份数,切比雪夫分割法计算精度高.关 键 词:铣削;切比雪夫分割法;传递矩阵;二阶Newton差商法;稳定性叶瓣图中图分类号

3、: TG 502. 14 文献标志码: A 文章编号: 1005 - 3026( 2016) 04 - 0538 - 05M illing Stability Analysis Based on Chebyshev SegmentationLI He, HAN Ping, ZHANG Ling- li( School of M echanical Engineering & Automation, Northeastern University, Shenyang 1 10819, China.Corresponding author: LI He, professor, E-mail: hl

4、i mail. neu. edu. cn)Abstract: Using the Chebyshev segmentation method to discrete the time period contained indelay differential equation, the Newton second- order difference quotient method was used tocalculate the motion vector of cutter at each time endpoint, then the transfer matrix of millings

5、ystem was obtained, and the system stability was determined by the Eloquet theory. Using theabove methods, a two- degrees of freedom milling system stability issues was investigated to getsystem stability lobe diagrams. The results showed that under the same calculation accuracy, thepoints needed to

6、 cut the time period by the Chebyshev segmentation are less than those by theaverage segmentation, and the computational efficiency ofthe Chebyshev segmentation is higher.W hen the time period is divided into the same parts, the stability lobe diagrams gotten by theChebyshev segmentation method are

7、more accurate than that ofthe other.Key words: milling; Chebyshev segmentation; transfer matrix; Newton second- order differencequotient method; stability lobe diagram在模具、电子、航空、航天等领域,常常要求高效、高精度地从毛坯中铣削去除大量材料.实践表明,颤振是影响高效、高精度铣削的重要因素之一 1 .颤振使得刀具与工件的相对位置发生改变,降低铣削效率和加工表面的精度,还会加剧刀具的磨损甚至破坏,缩短机床和刀具的使用寿命.因此,

8、铣削过程稳定性分析和铣削颤振的控制理论研究已经成为提高铣削加工效率的重要手段之一,具有重要的研究价值.铣削稳定性分析是合理选取主轴转速和切深等加工参数以达到优化工艺参数,避免发生颤振,提高生产效率的重要途径.铣削过程中多齿刀具旋转而引起的切削力周期变化和切削过程间断等因素,使得铣削稳定性分析比连续正交切削稳定性分析困难得多.目前国内外学者针对这类铣削稳定性问题,已经提出许多近似计算方法预测由临界加工参数构成的铣削过程稳定性边界曲线即稳定性叶瓣图( stability lobe diagram, SLD) 2 .文献 3提出了用于判断延时系统稳定性的半离散法( semi -discretizat

9、ion method, SDM ) .文献 4针对单自由度断续切削的情况,提出了时间有限元分析法AcroPDF - A Quality PDF Writer and PDF Converter to create PDF files. To remove the line, buy a license. ( temporal finite element analysis, TEEA) .文献 5基于铣削过程动态响应的直接积分格式提出了适用于大(小)径向切深、大(小)轴向切深等多种工况半解析的全离散法( full - discretizationmethod, EDM ) .文献 6证明,全离

10、散法与零阶半离散法具有同阶局部截断误差.以上的铣削稳定性预测方法,不论是半离散法、时间有限元法还是全离散法都是将铣削延时微分方程的时间周期平均分割.本文提出的切比雪夫分割方法是将时间周期T不均匀分割成m个时间元,刀具任意时刻的状态项和时间周期项用时间元端点的状态项线性表示,延时项采用二阶Newton差商的方法用时间元端点的状态项插值表示.首先求出铣削系统在每个时间元上的传递矩阵,再求铣削系统在一个时间周期内的传递矩阵,再用Eloquet理论 7来判断铣削系统的稳定性.通过时域内仿真结果可以看出:在相同精度要求下,切比雪夫不平均分割时间周期方法的计算效率比平均分割时间周期方法的计算效率高.1 切

11、比雪夫分割法逼近法近似求函数解析表达式是在函数定义域区间段上取n + 1个互异的节点,给定这些节点的函数值,采用Lagrange或者Newton插值方法 8构造出一个唯一的n次插值多项式来逼近原函数.采用均匀取点法构造出来的n + 1个节点的插值型积分公式一般具有n次代数精度.如果适当地选取节点来构造求积公式,就可以得到代数精度更高的求积公式.切比雪夫分割法 9就是这样一种取点分割方法,它的取点分割过程如下.如图1所示,将半圆弧按照角分割,圆弧上的分割点在x轴上的投影点即是该区间段 - 1,1上的分割点. = 2N0 + 1 - 2k2( N0 + 1) ( k = 0, 1, , N0) ,

12、 ( 1)x k = cos = cos 2N0 + 1 - 2k2( N0 + 1) , ( 2)在区间 a, b上, x轴上取点的位置表达式为xk = b - a2 x k + b + a2 =b - a2 cos2N0 + 1 - 2k2( N0 + 1) +b + a2 ,k = 0, 1, , N0. ( 3)图1 切比雪夫分割图( N0 = 4)Fig. 1 Image of Chebyshev segmentation2 全离散法只考虑再生颤振的铣削延时微分方程为觶X( t) = AX( t) + B( t) X( t) - X( t - T) .( 4)式中: A为铣削系统特性

13、的常数矩阵; B为再生效应动态铣削力决定的周期矩阵B( t) = B( t + T) ;T为时间周期等于延时量.切比雪夫分割方法将时间周期T不均匀分割为m个时间元:T= m- 1k= 0 tk( k = 0, 1, , m - 1) ; ( 5) tk = tk + 1 - tk = T2 ( tk + 1 - tk) ,( k = 0, 1, , m - 1) ; ( 6)tk = cos 2m + 1 - 2k2( m + 1) ; ( 7)tk + 1 = cos 2m + 1 - 2( k + 1)2( m + 1) . ( 8)在第一个时间元内tk t tk + 1 ( k = 0,

14、 1, ,m - 1) ,式( 4)的响应可以写成直接积分格式:X( t) = eA( t- tk) X( tk) +乙 ttk eA( t- ) B( ) X( ) d -乙 ttk eA( t- ) B( ) X( - T) d . ( 9)X( tk)表示在tk时刻刀具的状态向量 6 ,则在tk + 1时刻刀具的状态向量X( tk + 1 )为X( tk+ 1 ) = eA( tk) X( tk) +乙 tk0 eA B( tk+ 1 - ) X( tk+ 1 - ) d -乙 tk0 eA B( tk+ 1 - ) X( tk+ 1 - - T) d .( 10)此时式( 10)中的积

15、分时间区间为 0, tk .对周期系数矩阵B( tk + 1 - ) , X( tk + 1 - )在时间元tk t tk + 1上做线性逼近 10 12 ;对于延时项X935第4期 李 鹤等:基于切比雪夫分割法的铣削稳定性分析AcroPDF - A Quality PDF Writer and PDF Converter to create PDF files. To remove the line, buy a license. ( tk + 1 - - T)可以用时间元 tk - T, tk + 1 - T两端点时刻铣刀的状态向量X( tk + 1 - T)和X( tk - T)采用二阶

16、Newton差商的方法线性表示,见式( 11) ;将式( 11)带入式( 10)得式( 12) .X( tk + 1 - ) = 2 - tk tk( tk + tk - 1 ) X( tk - 1 ) + ( tk + tk - 1 - ) tk tk - 1 X( tk) +( tk - 1 + tk - ) ( tk - ) tk( tk + tk - 1 ) X( tk + 1 ) ,( 1 1)X( tk + 1 ) = Fk + 1Xk + 1 + ( F0 + F0, k) Xk +Fk - 1Xk - 1 - Fm - 1Xk + 1 - m - FmXk - m, ( 12)

17、 0 = eA tk, 1 = 乙 tk0eA d , 2 = 乙 tk0eA d , 3 = 乙 tk0eA 2d , 4 = 乙 tk0eA 3d ,式( 12)中的变量Fk + 1 , Fk - 1 , F0, k分别为Fk + 1 = ( 1 - 2 tk- 2 tk + tk - 1+ 3 tk( tk + tk - 1 ) ) B( k)0 + ( 2 - 3 tk - 3 tk + tk - 1 + 4 tk( tk + tk - 1 ) ) B( k)1 ,Fk - 1 = ( 3 tk( tk + tk - 1 )- 2 tk + tk - 1 ) B( k)0 + ( 4

18、tk( tk + tk - 1 ) - 3 tk + tk - 1 ) B( k)1 ,F0, k = ( 2 tk - 1+ 2 tk- 3 tk tk - 1) B( k)0 +( 3 tk - 1+ 3 tk- 4 tk tk - 1) B( k)1 .设 I - Fk + 1 - 1存在,则铣削系统每个时间元上的传递矩阵 k可以表示为 k = I - Fk + 1 - 1 ( F0 + F0, k) I - Fk + 1 - 1Fk - 1 0 0 - I - Fk + 1 - 1Fm - 1 - I - Fk + 1 - 1FmI 0 0 0 0 00 I 0 0 0 0左 左 左

19、左 左 左0 0 0 0 0 00 0 0 I 0 00 0 0 0 I 0.铣削系统在一个时间周期T上的传递矩阵可以用这个周期内时间元的传递矩阵序列 k( k = 0, 1, , m - 1)表示,即 = m - 1 m - 2 .最后,根据Eloquet理论,铣削系统的稳定性可由系统传递矩阵的特征值决定:若传递矩阵的所有特征值的模长均小于1,则系统稳定;否则系统不稳定.3 切比雪夫分割法下的稳定性分析3. 1 铣削系统参数两自由度铣削系统如下 12 :两个刀齿的槽式直铣刀;铣刀齿数N = 2;质量mtx = 0. 039 93 kg,mty = 0. 039 93 kg;阻尼比 = 0.

20、011;系统的固有频率fn = n/ ( 2 ) = 922 Hz;铣刀转速变化范围为0. 5 104 2. 5 104 r/ min;切深w的变化范围为0 0. 01 m; Kt = 6 108 N/ m2; Kn =2 108 N/ m2;铣削系统特性的常数矩阵A:A =- 63. 72 0 25. 04 00 - 63. 72 0 25. 04- 1. 34 106 0 - 63. 72 00 - 1. 34 106 0 - 63. 72.再生效应动态铣削力决定的周期矩阵B:B( t) =0 0 0 00 0 0 0- whxx( t) - whxy( t) 0 0- whyx( t)

21、- whyy( t) 0 0.Kt是线性化切向切削力系数, Kn是线性化法向切削力系数, a/ D表示铣刀的径向侵入比.3. 2 稳定性分析结果将转速-轴向切深平面平均划分为200 100网格点,即转速步长100 r/ min,轴向切深步长为0. 1 mm.分别计算径向侵入比等于0. 1和0. 2两种工况下的铣削稳定性极限.根据前面所述的理论基础和算法原理,编写计算程序,利用计算机绘制稳定性叶瓣图,如表1所示.045东北大学学报(自然科学版) 第37卷AcroPDF - A Quality PDF Writer and PDF Converter to create PDF files. To

22、 remove the line, buy a license. 表1 两自由度铣削模型的算法比较TabIe 1 AIgorithm comparison of two degrees of freedom miIIing modeI切比雪夫分割法下的全离散( m = 30)平均分割法下的全离散法( m = 40)a/ D = 0. 1计算时间/ s 242. 819 501. 195a/ D = 0. 2计算时间/ s 251. 978 509. 060由表1可以看出在用叶瓣图来预测铣削系统的稳定性时,切比雪夫分割法取的点数比平均分割法取的点数少的情况下,叶瓣图仍能达到相同的精度.由两种分割

23、方法的计算时间的对比可知,获得相同精度叶瓣图时,切比雪夫分割法的计算效率比平均分割法计算效率高.图2表示a/ D = 0. 1, m = 40时切比雪夫分割法和平均分割法得到的稳定性叶瓣对比图,这两种分割方法得到的铣削系统的稳定性范围有所不同.在图2中取4个工作点:点1,点2,点3,点4,这4个工作点对应的铣削参数和两种分割方法预测的4个工作点的稳定性如表2所示.对表2中所列工作点进行时域仿真,仿真结果如图3所示.工作点的仿真结果如图3所示.由图3可知,工作点1,点2,点3对应的铣削过程是稳定的,点4对应的铣削过程是不稳定的.图3表明4个工作点的时域仿真结果与表2中切比雪夫分割法预测的4个工作

24、点的稳定性完全吻合,而与平均分割法预测的4个工作点的稳定性相差较多.因此,取相同分割点时,切比雪夫分割法对铣削系统的稳定性预测结果比平均分割法对铣削系统的稳定性预测结果更准确.图2 两种分割方法下的稳定性叶瓣图对比(顺铣N= 2, a/ D= 0. 1 )Fig. 2 Comparison of two different segmentation sstabiIity Iobes diagram表2 两种取点分割法得到的部分工作点的稳定性TabIe 2 StabiIity of sectionaI working points for twoway to take segmentation工

25、作点r min- 1 wmm切比雪夫分割法平均分割法点1 15 000 1. 5稳定不稳定点2 20 000 1稳定稳定点3 20 000 2稳定不稳定点4 20 000 3不稳定不稳定145第4期 李 鹤等:基于切比雪夫分割法的铣削稳定性分析AcroPDF - A Quality PDF Writer and PDF Converter to create PDF files. To remove the line, buy a license. 图3 部分工作点的时域仿真结果图Fig. 3 Time domain simuIation s resuIts for partof workin

26、g points( a) 点1; ( b) 点2; ( c) 点3; ( d) 点4.4 结 论本文介绍了切比雪夫分割方法和Newton二阶差商法对离散项进行插值积分的方法.分析两自由度铣削系统示例,由绘制的系统稳定性叶瓣图可知:在达到相同预测精度的条件下,切比雪夫不平均分割法需要取的点数比平均分割法所需的点数少,并且切比雪夫不平均分割法的计算效率比平均分割法计算效率高;在分割点数相同的情况下,切比雪夫不平均分割法预测的系统稳定性比平均分割法预测的系统稳定性更加准确.参考文献: 1 Tlusty J, Polacek M . The stability of machine tools aga

27、instself- excited vibrations in machining J . InternationalResearch in Production Engineering, ASM E, 1963, 1 ( 8 ) :465 - 474. 2 M errit H E. Theory ofself- excited machine tool chatter J .ASM E Journal of Engineering for Industry, 1965, 8 ( 7 ) :447 - 454. 3 Insperger T, St p n G. Semi- discretiza

28、tion method fordelayed systems J . International Journal for NumericalM ethods in Engineering, 2002, 55( 5) : 503 - 518. 4 Bayly P V, Halley J E, M ann B P, et al. Stability ofinterrupted cutting by temporal finite element analysis J .Journal of M anufacturing Science and Engineering,Transactions of

29、the ASM E, 2003, 125( 2) : 220 - 225. 5 Ding Y, Zhu L M , Zhang X J, et al. A full-discretizationmethod for prediction of milling stability J . InternationalJournal of M achine Tools and M anufacture, 2010, 50 ( 5 ) :502 - 509. 6 Insperger T. Eull-discretization and semi- discretization formilling s

30、tability prediction: some comments J .International Journal of M achine Tools and M anufacture,2010, 50( 7) : 658 - 662. 7 Earkas M . Periodic motions M . New York: Springer-Verlag, 1994. 8 张铁,阎家斌.数值分析 M .北京:冶金工业出版社,2007: 141 - 187.( Zhang tie, Yan Jia-bin. Numerical analysis M . Beijing:M etallurgi

31、cal Industry Press, 2007: 141 - 187. ) 9 Yang W Y, Cao W W , Chung T S, et al. Applied numericalmethods using M ATLAB M . New York : John W iley &Sons, 2005: 132 - 135. 10 Tan S J, Zhong W X. Precise integration method for duhamelterms arising from non-homogenous dynamic systems J .Chinese Journal o

32、f Theoretical and Applied M echanics,2007, 39( 3) : 374 - 381. 1 1 Zhong W X, W illiams E W . A precise time step integrationmethod J . Proceedings of the Institution of M echanicalEngineers. Part C: M echanical Engineering Science, 1994,208( 6) : 427 - 430. 12 Ding Y, Zhu L M , Zhang X J. Second- order full-discrimination method for milling stability prediction J .International Journal ofM achine Tools & M anufacture, 2010,50( 10) : 926 - 932.245东北大学学报(自然科学版) 第37卷AcroPDF - A Quality PDF Writer and PDF Converter to create PDF files. To remove the line, buy a license.