作者简介:戴玉琦(1995—),深圳大学硕士研究生.研究方向:超精密加工技术.E-mail:daiyuqi2017@email.szu.edu.cn
中文责编:方 圆; 英文责编:溯 心
College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen 518060, Guangdong Province, P.R.China
cutting theory; diamond turning; tool height error; center cone; tool interference; cutting force
DOI: 10.3724/SP.J.1249.2020.05543
为探究工件中心刀具干涉的形成机理,实现刀高误差的在线辨识,围绕刀高误差存在情况下,刀具干涉的产生及加工过程参数的影响等方面展开研究.通过分别建立车削行程中工件中心圆台及刀具后刀面的三维形貌方程,基于两者间的位置关系推导刀具干涉区域半径的数学模型,探讨加工过程参数对刀具干涉区域三维形貌的影响.根据所推导的数学模型,利用刀具与工件挤压面积的投影来表征切削力的变化,建立刀高误差的预测模型.理论与试验研究结果表明,刀高误差、切削深度、刀尖半径及刀具后角直接影响刀具干涉区域的范围与形貌,工件刀高误差可通过监测切削行程中的切削力进行在线辨识.
In order to explore formation mechanism of the tool interference at the workpiece center and realize the online identification of tool height error, tool interference generation and the influence of the processing parameters were investigated under the condition of the tool height error. By establishing the equations of the three-dimensional(3D)topography of round platform and the tool clearance face respectively during the turning process, a mathematical model for calculating the radius of tool interference zone was derived based on the relationship between the two positions. Meanwhile, the influence of processing parameters on the 3D topography of the tool interference zone was discussed. According to the derived mathematical model, the cutting force was characterized by the projection of the extrusion area between tool and workpiece, whereby a prediction model for the tool height error was established. The theoretical and experimental results show that the tool height error, cutting depth, tool nose radius and tool clearance angle directly affect the range and shape of the tool interference zone. The tool height error can be identified online by monitoring cutting force during the turning process.
单点金刚石车削是一种利用单晶金刚石刀具进行超精密加工的技术手段[1].由于单晶金刚石具有纳米级的边缘锋利度和极高的硬度,可直接用于加工具有纳米级表面粗糙度和亚微米级面形精度的超光滑表面[2],被广泛应用于光学工程、航空航天及军事装备等相关领域[3- 4].在单点金刚石车削过程中,机床误差的存在直接影响加工表面的面形精度及表面粗糙度.单点金刚石车削中的机床误差通常包括主轴振动误差、导轨行进误差、刀具振动误差、机床热误差及刀具偏移误差等[5-7].随着对光学产品表面加工精度的要求越发严格,单点金刚石车削中表面生成与误差控制研究日趋重要[8].
近年来,国内外诸多学者对单点金刚石车削中误差来源及对加工表面粗糙度的影响进行深入研究.ZHANG等[9]基于轴承流体动力学建立主轴电机磁力分析模型,研究主轴不平衡引起强迫振动的动态特性及其对加工表面的影响.KHIM等[10]利用传递函数和轨道形状误差之间的几何关系,探究线性导轨形状对运动误差的影响.WANG等[11]通过对切削力信号的监测提出一种具有刀尖振动效应的几何模型,指出刀尖振动和过程阻尼效应是影响工件表面粗糙度的主要因素.ABDULSHAHED等[12]提出一种基于自适应神经模糊推理系统(Anfis)的热误差建模法,使机床精度提升80%.然而,有关单点金刚石车削过程中刀具中心误差的研究却鲜有报道.
在单点金刚石车削中,加工表面成形于刀具几何形貌沿其轨迹上的复印,故刀具与工件中心偏离所产生的刀具中心误差,将直接影响加工表面质量.相关研究表明,因刀具在水平方向上没有对中所产生的中心偏移误差可分为刀具欠中心偏差和刀具过中心偏差,其根据加工表面形状的凹凸不同将导致W形误差或M形误差[13].刀具在竖直方向上没有对中所产生的误差可分为刀高误差与刀低误差,前者将导致在切削行程末端,因刀具后刀面与工件之间的干涉在工件表面中心形成中心圆锥,而后者将导致工件表面中心圆柱的形成.与以往研究不同,本研究发现当存在刀高误差时,加工表面中心圆锥周围将形成圆台状的特征凸起.该圆台的高度高于工件表面,且相较于中心圆锥具有更大区域,但在中心圆柱的形成中却未见相似特征.为此,本研究根据刀具与工件之间的切削运动学关系建立刀具干涉区域计算几何模型,探讨刀具干涉的产生机理及演变规律.通过对切削过程中的切削力分析建模,进一步阐明干涉过程中刀具-工件的干涉面积与切削力之间的作用关系,建立基于干涉作用力的刀高误差预测模型.
1.1 中心圆台侧面方程在车削过程中,主轴转速恒定,工件通过真空吸盘被吸附在主轴上,金刚石刀具以恒定的进给速率完成对工件的端面切削.当存在刀高误差时,为描述车削过程中刀具与工件的空间位置,建立如下坐标系:以工件中心为原点o建立工件坐标系o-xyz, 其中, x轴正向与刀具进给方向相反, z轴正方向垂直于工件表面向外; 以刀尖为原点o'建立刀具坐标系o'-x'y'z', 其各坐标轴方向与坐标系o-xyz相同,如图1.
本研究所采用金刚石刀具的切削刃为圆弧切削刃,其几何形状可在坐标系o'-x'y'z'中用圆弧方程进行描述,如图1(b).经过坐标系平移变换,切削刃在工件坐标系o-xyz中可表示为
{(x1-s)2+(z1-R)2=R2
y1=h
0≤z1≤ap(1)
其中, 点(x1, y1, z1)为切削刃上任意点坐标; R为刀尖圆弧半径; ap为切削深度; h为刀高误差; s为切削刃尖端与原点o之间的水平距离.
因中心圆台由刀具圆弧切削刃切削工件而形成,故其侧面方程为刀具圆弧切削刃方程绕工件坐标系z轴旋转而形成的回转方程.刀具坐标系o'-x'y'z'相对于工件坐标系o-xyz在x和y方向分别具有偏移量s和h. 图1(b)中切削刃上任意一点的回转半径R'为
R'=(x2+y2)1/2=(x21+h2)1/2(2)
结合式(1)和式(2)可得中心圆台的侧面方程为
{((x2+y2-h2)1/2-s)2+(z-R)2=R2
0≤z≤ap(3)
在单点金刚石车削中,所用金刚石刀具后刀面为圆锥形后刀面,其几何示意图如图2.
根据图2的刀具几何形貌,考虑实际加工中切削深度的限制,后刀面在工件坐标系o-xyz下的几何方程可通过坐标变换表示为
{Z(x-s)2+(z-R)2=((y-h+Rcot δ)/(cot δ))2
-(2Rap-a2p)1/2≤x-s≤(2Rap-a2p)1/2
(-apcot δ+h)≤y≤h(4)
其中,δ为刀具后角.
在实际加工中,为达到所需的表面微纳结构,通常需要对工件进行多层车削.对于首层表面车削,中心圆锥完全由中心圆台侧面与工具后刀面之间的接触干涉生成.而对于多层车削,上层车削所形成的中心圆锥可能首先会接触到刀具后刀面,从而影响本层切削中刀具干涉状态,进而改变刀具干涉的初始位置,如图3.可见,当刀具存在刀高误差时,会导致车削表面始终残留中心圆锥微结构特征,该特征在o-xyz坐标系中可表示为
图3 刀具与上层车削残留中心圆锥示意图
Fig.3 Schematic of the relation between tool and residual center cone generated from upper layer turning
{x2+y2=(h-(z-ap)cot δ)2
ap≤z≤htan δ+ap(5)
给定切削参量为h=120 μm、 ap=10 μm、 δ =10°及R =0.281 mm时,计算得刀具后刀面与上层车削残留中心圆锥之间初始干涉点位置坐标为x=86.76 μm、 y=82.90 μm,而刀具后刀面与中心圆台侧面初始干涉位置坐标为x=179.90 μm、 y=120.00 μm.比较可得后者干涉初始位置离工件中心更远,故刀具后刀面与中心圆台侧面先发生干涉,再与上层车削所形成的中心圆锥发生2次干涉.
随着刀具以恒定速率从工件边缘进给到工件中心,中心圆台的直径逐渐变小.刀具后刀面会在刀具进给的某一位置与工件中心圆台侧面发生干涉.对于刀具进给到任意处形成的工件中心圆台,其截面圆半径随着中心圆台高度的增加逐渐减小,而曲率逐渐增大,即当ap3<ap2<ap1时,工件中心圆台截面圆弦切角α3<α2<α1, 如图4.因此,最大弦切角出现在中心圆台的顶圆平面内,弦切角在顶圆平面内顶点坐标为(s-(2Rap-a2p)1/2, h, ap). 顶圆平面内的弦切角可通过求其正切值进行计算,即
(∂ y)/(∂ x)=-(s-(2Rap-a2p)1/2)/h(6)
对于金刚石刀具而言,随着切深的增加,切深平面与刀具后刀面交线的曲率也相应变化.如图4所示,随着切深的增加(ap3<ap2<ap1), 刀具后刀面与中心圆台在切深平面交点处的弦切角β会越小(β3>β2>β1), 刀具后刀面最大弦切角出现在中心圆台顶面内.当切深平面与中心圆台顶面重合时,刀具后刀面与中心圆台顶面交点(s-(2Rap-a2p)1/2, h, ap)处交线弦切角β的正切值为
(∂ y)/(∂ x)=-((R2-(ap-R)2)1/2)/(Rtan β)(7)
图4 中心圆台不同截面处刀具与工件干涉位置
Fig.4 Interference location diagram of tool and workpiece at different sections of center circular table
由于β与α成负相关关系,故可判断α1=β1为刀具干涉发生的临界点,刀具干涉最先出现在中心圆台顶圆平面内.此外, α随着刀具从工件边缘进给到工件中心越来越大,而β在此过程中保持不变,则刀具干涉临界点发生在α取值最大时.
联立式(6)和式(7),便可求得刀具发生临界干涉时所处的水平位置为
is=(h/(Rtan δ)+1)(2Rap-a2p)1/2(8)
当刀具在车削过程中存在刀高误差h时,刀具干涉区域半径为
r=((2Rap-a2p)(h/(Rtan δ)+1)2+h2)1/2(9)
由式(9)可见,刀尖半径、刀具后角、切削深度及刀高误差将直接影响刀具干涉区域的范围.
本研究中工件的端面切削试验由Nanotech 450UPL单点金刚石车床完成,试验搭建如图5.加工工件通过夹具夹持后吸附在真空吸盘上,金刚石刀具安装在3通道Kistler 9017C测力计上,测力计通过夹具固定在可调刀座上.金刚石刀具在机床y轴方向的刀高误差通过机床附带可调刀座与千分表配合进行精确调节,直接得到期望的刀高误差.为避免刀具水平方向误差对刀高误差的分析产生影响,需在对刀过程中将其控制在0.1 μm以内.
1为真空吸盘; 2为工件夹具; 3为工件; 4为金刚石刀具;
5为测力计; 6为喷嘴; 7为切削力采集装置;
8为可调刀座; 9为金相显微镜
为探究加工过程参数与刀具干涉区域拓扑形貌之间的交互关系,车削试验在不同加工过程参数及刀高误差下进行.表1为开展车削试验的加工过程参数,采用金刚石刀具主要从刀具后角、刀尖圆弧半径、切削深度及进给速率对刀具干涉的影响展开.在完成对工件端面的切削后,分别利用光学显微镜(Leica DM 2700M)和白光干涉仪(Bruker GT-X)测量工件表面,观察工件表面因刀具干涉所引起的几何特征创成.
在单点金刚石车削过程中,当刀具进给轨迹低于工件回转中心(h<0)时,由于切削刃无法完成对工件中心部分的切削将在工件表面残留中心圆柱,如图6(a); 刀具进给轨迹高于工件回转中心(h>0)时,将在工件表面残留中心圆锥,如图6(b).由白光干涉仪的测量结果可知,当工件表面形成中心圆锥时,在其周围会形成圆台状凸起的干涉区域,而中心圆柱周围却无类似的干涉区域.
图7为不同切削深度与进给速度下的刀具干涉区域示意图.由图7(a)可见,从凸状干涉圆台边缘到中心圆锥底圆,其截面曲线可划分为4个曲线段,记为S1、 S2、 S3和S4. 其中, S1和S3较为平缓,而S2和S4曲线斜率较大,大斜率截面曲线的形成是由于金刚石刀具在干涉作用力下,在该段切削区域经历被动抬升,因此, S2和S4曲线可以反映刀具干涉作用力幅值与方向的变化情况.
从S1曲线刀具干涉开始到S2刀具被动抬升,干涉作用力有一个明显增加的过程,可认为刀具的被动抬升作用力有一个阈值.当干涉作用力垂直加工表面的分力大于该阈值时,刀具被动抬升从而形成干涉圆台; 反之,即使刀具发生了干涉,因干涉作用力小于刀具抬升阈值,也无法使金刚石刀具抬升而形成干涉圆台.当改变切削深度与刀尖半径时,中心圆锥周围的干涉圆台形貌变成凹状干涉圆台,如图7(b).从凹状干涉圆台边缘到中心圆锥底圆,凹状干涉圆台的截面曲线只划分为S1和S2段,并无平滑过渡的干涉曲线段,证明凸状或凹状干涉圆台的形成与加工过程参数有关.
刀具干涉区域是在刀具后刀面与中心圆台之间干涉力作用下,由刀具的被动抬升所形成.试验研究发现,刀具干涉作用力与刀具干涉区域拓扑形貌之间存在密切联系.因此,可通过监测切削力实现刀高误差的在线辨识.
图9为中心圆柱和中心圆锥形貌及与之对应切削力形态.由图9(a)可见,当工件表面形成中心圆柱时,切削力在切削行程中无明显突变; 而当工件表面形成中心圆锥时,切削力在切削行程末端有显著突变,呈力脉冲形式.由图9(b)和(c)可见,当h从29.84 μm变为49.91 μm时,切削力的3个分量力脉冲幅值均有较大增幅.因此,在中心圆锥形成过程中,刀具后刀面对工件材料挤压而产生刀具干涉,从而在切削行程末端产生脉冲形式的干涉作用力.干涉作用力与刀高误差关系密切,刀高误差越大,干涉作用力幅值也越大.
图 10为不同刀高误差下刀具干涉过程中的干涉作用力形态.由图 10(a)可见,在中心圆柱形成过程中,切削行程末端未发现脉冲状切削力; 而在中心圆锥形成时,随着刀高误差的增加,其3个方向上的切削力分量均有很大增加,其中, Fx的幅值改变最明显,代表了干涉作用力的敏感方向.表3比较了h=29.84 μm及h=49.91 μm时的切削力幅度, h=49.91 μm处的干涉作用力分量Fx, Fy和Fz均大于h=29.84 μm处的干涉作用力分量,且这两处刀高误差下各干涉作用力分量的比值(F49.91/F29.84)分别为2.15、2.05及2.26,稳定在2.1左右.证明干涉作用力于各个方向上的分量的变化趋势也基本相同.
图 11为刀具干涉区域投影图.由图 11(a)可见,刀具与工件的干涉位置最先发生在中心圆台顶面与刀具切削刃相交的点处,随着刀具的持续进给,干涉边界一方面会在中心圆台顶面内沿刀具后刀面方向迁移,另一方面会沿着切削刃向圆锥底部迁移.由于在干涉过程中切削力分量Fx的变化最明显,抗误差干扰能力最强,故利用Fx与干涉区域形貌之间的映射关系表征刀高误差.由图 11可见,干涉区域在xoy平面的投影为封闭曲面ABC. 因封闭曲面ABC的BC边长度较短,可近似为直线,故可将封闭曲面ABC简化为三角形ABC.
点A为中心圆台顶面与刀具前刀面相交点,其位于切削刃上,则其x轴坐标为
xA={s-(2Rap-a2p)1/2,(2Rap-a2p)1/2≤s≤is
0, 0≤s≤(2Rap-a2p)1/2(10)
由图 11(b)可见,点B为干涉点向下迁移的边界点.对于xoy截面与刀具后刀面交线的斜率为
(∂ y)/(∂ x)=((x-s)cot2δ)/(y-h+Rtan δ)(11)
因点B在切削刃上,其y值恒为h. 同时,对于刀具后刀面和中心圆台,两者干涉点处斜率相等,代入式(11)可得
-x/h=(x-s)/(Rtan δ)(12)
求解式(12)可得到点B的x轴坐标为xB=(sh)/(Rtan δ+h). 点C为中心圆台顶面圆与刀具后刀面的交点,其坐标点(xc,yc,zc)为
{(xc-s)2+(zc-R)2=((yc-h+Rtan δ)/(cot δ))2
(s-xc)/((yc-h+Rtan δ)tan δ2)=(xc)/(yc)
zc=ap(13)
将三角形ABC的面积带入下式,便可计算干涉作用力幅值为
Fx={K((-sRtan δ)/(Rtan δ+h)+c)((h-yc)/2), c≤s≤is
(Ksh)/(Rtan δ+h)((h-yc)/2), 0≤s≤c(14)
其中, c=(2Rap-a2p)1/2; K通过标定获得, K=9.2×109 N/m2.
用时间t替换式(14)中的水平距离s, 可得到干涉作用力随时间变化的表达式为
Fx=
{K[((tfr-is)Rtanδ)/(Rtanδ+h)+c]((h-yc)/2), c/(fr)≤t≤(is)/(fr)
(Kh(is-tfr))/(Rtanδ+h)((h-yc)/2), 0≤t≤c/(fr)(15)
其中, fr为刀具的进给速率.
图 12为在相同切削参量下(刀尖半径为0.5 mm、刀具后角为15°、切削深度为10 μm及刀高误差为49.91 μm),Fx随时间变化的模拟值与试验值对比结果.可见,随着刀具的持续进给,刀具干涉作用力先增加后减小.
图 12 干涉作用力分量Fx的模拟值与试验值对比图
Fig.12 Comparison of simulated and experimental values of interference force Fx
图 13 不同刀高误差下切削力分量Fx幅值
Fig.13 The amplitude of the cutting force Fx under different tool height errors
对比不同刀高误差下(刀尖半径为0.5 mm、刀具后角为15°及切削深度为10 μm)刀具干涉作用力分量Fx的模拟值与试验值,如图 13.随着刀高误差的增加, Fx分量的模拟值与试验值都相应地趋于线性增加,且两值在增加过程中吻合较好,证明利用干涉作用力进行刀高误差表征方法可行.因此,加工中可通过监测车削行程中的切削力,无需拆卸工件即可将刀高误差值控制在极小范围内.
本研究根据工件与刀具之间的切削运动学关系建立刀具干涉区域的几何计算模型,系统讨论刀具干涉的演化生成机理,利用刀具与工件之间挤压面积的变化推导刀具干涉行程中力学特征模型.结合试验分析比较发现,因刀具抬升阈值的存在,刀具干涉区域理论值均略大于试验值,且该几何计算模型的最大相对误差为6.71%; 通过对刀具干涉行程中切削力的监测可实现刀高误差的在线辨识,将极大提升单点金刚石车削加工的效率,也为实现智能化加工提供理论基础.
深圳大学学报理工版
JOURNAL OF SHENZHEN UNIVERSITY SCIENCE AND ENGINEERING
(1984年创刊 双月刊)
主 管 深圳大学
主 办 深圳大学
编辑出版 深圳大学学报理工版编辑部
主 编 李清泉
国内发行 深圳市邮电局
国外发行 中国国际图书贸易集团有限公司(北京399信箱)
地 址 北京东黄城根北街16号
邮 编 100717
电 话 0755-26732266
0755-26538306
Email journal@szu.edu.cn
标准刊号 ISSN 1000-2618
CN 44-1401/N