基于机床刀具加工变形研究的铣削工艺参数优化方法
2015-2-6 来源：数控机床市场网 作者：上海理工大学机械工程学院迟玉伦李郝林
摘要: 针对机床刀具变形影响加工精度问题，据铣床刀具系统动刚度测量计算及铣削力模型提出基于遗传算法的铣削工艺参数优化方法。通过阶跃响应实验获得机床刀具系统动刚度及用正交实验法可准确铣削力模型，对实验机床铣削工件工艺参数进行数学优化;开发设计在线测量刀具变形及铣削力实验装置并分别对优化前后工艺参数实验验证，结果已证明该方法的有效性、实用性。
关键词: 铣削;刀具变形;工艺优化;实验设计
数控铣削加工过程中在切削力作用下机床刀具变形引起的刀具与工件间相对位置变化成为影响机床加工精度的重要因素。已有对铣削刀具变形建模补偿进行的诸多研究。张臣等以球头铣刀为研究对象，建立球头铣刀刀具变形引起的加工误差仿真预测模型。Denkena 等提出适应主轴系统铣削加工过程中对刀具静、动态变形进行补偿并取得较好效果。
Habibi 等通过对铣削刀具切削变形建模，经数控系统程序对刀具切削路径进行补偿，并实验证明该方法的有效性;但铣削刀具变形补偿研究中建模控制过程复杂，难以实际应用。
因此，本文提出通过优化铣削工艺参数减小刀具变形量、提高加工精度方法。
本文据铣床刀具系统动刚度测量计算及铣削力模型，提出基于遗传算法的铣削工艺参数优化方法。通过阶跃响应实验获得机床刀具系统动刚度，并利用正交实验法准确获得铣削力模型;对实验机床铣削工件主轴转速、进给速度、切削宽度及切削深度等工艺参数进行数学优化;开发设计在线测量刀具变形及铣削力实验装置，并分别对优化前后工艺参数进行实验验证，结果证明该方法有效、实用。
1 机床铣削工艺参数优化模型
1. 1 机床刀具系统动刚度模型
图 1 的机床－刀具系统可简化为 X，Y 方向相互垂直的二自由度振动系统，即
式中:Mx，My分别为 X、Y 方向机床－刀具系统质量;Cx，Cy分别为 X、Y 方向机床－刀具系统阻尼;Ksx，Ksy分别为 X、Y 方向机床－刀具系统刚度;Fx，Fy分别为 X、Y方向刀具切削力分力。
经拉普拉斯转换可得
式中:Kdx，Kdy分别为 X、Y 方向机床－刀具系统动刚度，可表示为
式中: Ksx， Ksy分别为 X、Y 方向机床 － 刀具系统静刚度;ξx，y分别为 X、Y 方向机床－刀具系统阻尼比;λx，λy分别为 X、Y 方向机床－刀具系统交变频率 ω 与固有频率 ωn之比，ω =jvs60 ，j 为铣刀齿数，vs为主轴转速。
1. 2 铣削力模型
铣削磨削X、Y、Z 三方向平均铣削力 Fx、Fy、Fz计算式为
式中:vs为刀具主轴转速;fw为工件进给速度;ae为刀具切削深度;N 为刀具齿数;Ktc为刀具切向切削力系数;Krc为刀具径向切削力系数;Kte为刀具切向刃口力系数;Kre为刀具径向刃口力系数;Kac为刀具轴向切削力系数;Kae为刀具轴向刃口力系数;ex为切出角;st为切入角。顺铣时，切入、切出角可表示为
逆铣时，切入、切出角可表示为
式中:ap为刀具切削宽度;Ｒ为切削铣刀半径。
1． 3 铣削工艺参数优化方法
由加工工艺参数可确定铣刀 X，Y 方向铣削力大小。结合式(3)、(4)可求出在一定切削工艺参数下铣削刀具 X，Y 向的动态变形量为
结合铣削加工效率及尺寸要求，即可分析在铣削工艺参数 vs，ap，fw，ae下铣削刀具最大切削变形量。本文利用遗传算法对铣削工艺参数进行优化。设主轴转速为 vs，切削宽度为 ap，工件进给速度为 fw，刀具切削深度为 ae;每齿切削量 ft= fw/(jvs) 的约束条件为
正规化后，各铣削工艺参数作为个体的表现形式为
设世代数为 n，个体数为 k，变量个体位长为 l，交叉率为 m1，突然变异率为 m 。将所得遗传算法优化后主轴转速 vs，切削宽度 ap，工件进给速度 fw，刀具切削深度 ae作为机床加工工艺参数。
2 实验研究
2. 1 铣床刀具系统静刚度测量计算
本实验在沈阳机床厂VMC1165B 立式加工中心，工件材料 45 钢，铣刀型号 MG－2EA35M－D10－M。实验过程见图 2 (a)，压力传感器型号 KISTLEＲ9323AA,量程设为 1 kN，灵敏度－9．6 pc/N;电涡流位移传感器型号为 MICＲO － EPSILON eddyNCDT 3010，量程 1 mm，分辨率小于0．05 μm。通过力传感器对刀具刀尖位置施加 X 方向力 F1，电涡流位移传感器测出刀具刀尖变形量 Δx，求出刀具系统在 X 方向静刚度为Ksx=F1Δx=3．46 ×106N/m;同样得刀具系统在 Y 方向静刚度为 Ksy=3． 32 ×106N/m。
2. 2 铣床刀具系统固有频率及阻尼比测量计算
实验过程见图 3，压力传感器型号为 KISTLEＲ9301B，量程设为 2kN，灵敏度－4 pc/N。在刀具 X 方向施加拉力 F2= 100 N 作为系统输入信号;用剪刀将细绳瞬间剪断，利用电涡流位移传感器获得刀具位移振动输出信号。为减小拉绳弹性影响应尽量缩短细绳长度或其它部分改用刚性较大钢丝绳(图 3(b)中白色钢丝绳)。利用阶跃响应频响函数求出机床－刀具系统 X 方向固有频率 ωnx= 890 Hz，阻 尼 比 ξx=0.00227;同样可得刀具系统 Y 方向固有频率 ωny=886Hz，阻尼比 ξy=0.00342
据刀具系统各方向静刚度、固有频率及阻尼比测量结果，可得刀具系统 X，Y 方向动刚度关系为
2. 3 铣床力模型系数测量计算
用正交实验与偏最小二乘回归法确定铣削力模型中各系数值以获得准确的铣削力模型。实验见图 4。用 KISTLEＲ 9257B 三向测力仪对铣削加工进行测试，选 L9(34)正交试验方案安排各种测试。
选主轴转速、进给速度、铣削深度及铣削宽度等工艺参数不同数值作为该正交实验因素水平，见表 1。
在每组正交设计实验中用三向测力仪测得 X，Y 方向铣削力见表 2。
利用偏最小二乘回归法将铣削力模型与实验测量数据拟合求得铣削力模型各系数值为刀具切向切削力系数 Ktc=1 865．1，刀具径向切削力系数 Krc= 736．2，刀具切向刃口力系数 Kte=113．5，刀具径向刃口力系数Kre=158．7，由此获得较可靠的铣削力数学模型。
2．4 铣削工艺参数优化结果
设铣削各工艺参数约束条件为主轴转速 1 000
选世代数30，个体数100，变量个体位长 l =9，交叉率 0．6，突然变异率 0．001。得铣削工艺参数优化结果为主轴转速 vs=7 500 r/min，切削深度 ae=3．2 mm，切削宽度 ap=7．5 mm，工件进给速度 fw=301 mm/min。
2．5 实验验证
为验证铣削工艺参数优化方法的有效性、实用性，本文设计的在线测量铣削刀具变形及铣削力装置可同时测量铣削刀具变形及铣削力。采用 2 个电涡流位移传感器与 KISTLEＲ 9257B 三向测量仪在线测量铣削刀具加工过程中变形量及铣削力，实验方案见图 5。
通过专用夹具将电涡流位移传感器 1、电涡流位移传感器 2 固定在刀具处， 利用电涡流非接触式原理在线测量刀具位置变形量。电涡流传感器 1 测量刀具 X向变形，电涡流传感器 2 测量刀具 Y 向变形。两电涡流位移传感器测量信号、铣削力信号通过数据采集卡由计算机获得，采样频率设为 2 000 Hz。据两电涡流位移传感器测量结果，利用悬臂梁挠度计算式计算刀具刀尖变形量。将优化前后两组参数分别进行实验验证，结果见表 3。由表 3 已证明工艺参数优化方法的有效性。
3 结论
(1) 针对铣削加工过程中刀具变形问题本文据机床刀具系统动态特性及铣削力模型提出基于遗传算法的铣削工艺参数优化方法。
(2) 本文设计的机床刀具系统静、动刚度求解阶跃响应实验方法，据实验测量数据确定铣削力模型系数，基于实验测量结果对该机床铣削工艺参数进行优化。
(3) 为验证铣削工艺参数优化方法的有效性，提出在线测量铣削刀具变形及铣削力实验方法，分别对优化前、后工艺参数在线测量刀具变形及铣削力，实验结果已验证优化方法的有效性、实用性。
(4) 由于本文优化数学模型未考虑颤振因素，故该方法仅适用非颤振铣削开环系统的工艺参数优化。
投稿箱：
如果您有机床行业、企业相关新闻稿件发表，或进行资讯合作，欢迎联系本网编辑部， 邮箱：skjcsc@vip.sina.com
如果您有机床行业、企业相关新闻稿件发表，或进行资讯合作，欢迎联系本网编辑部， 邮箱：skjcsc@vip.sina.com
更多相关信息
