等误差直线逼近非圆曲线的点计算及刀具轨迹模拟.docVIP

等误差直线逼近非圆曲线的节点计算及刀具轨迹模拟/content/2009-06-26/56746.html 1 引言 由于大部分数控机床不具备对非圆曲线刀尖轨迹的插补指令，因此在编制此类曲线刀尖轨迹的数控程序时通常用直线段或圆弧段予以替代。由于直线替代法 简单、直观，因此使用较多。用直线段替代非圆曲线的方法如图1所示。在满足精度要求的条件下，可用折线段替代非圆曲线。图中a、b、c、d等称为节点，实 现刀尖轨迹数控编程的关键就是确定这些节点。为简化计算，常采用等间距法和等步长法来确定节点。等间距法是在理论曲线与直线的最大偏差小于允许偏差(?最大≤?允)的条件下，令各节点在x轴上的投影的间距?x相等。等步长法是在理论曲线与直线的最大偏差小于允许偏差(?最大≤?允)的条件下，令各节点间的直线长度?L相等。它们的共同特点是计算较为简单。但当各节点之间曲线的曲率变化较大时，由于?x和?L为定值，因此会造成被加工零件的表面粗糙度变化较大，从而影响工件的表面加工质量；同时，曲线曲率的变化也使工件的加工误差?发生变化。另一方面，等间距法的间距和等步长法的步长均是根据加工精度由非圆曲线的最小曲率半径确定的，因此这两种方法在整个非圆曲线内会产生很多节点，使计算和编程相当繁琐。如采用等误差直线逼近法则可有效避免上述问题。 2 等误差直线逼近的理论计算 如图1所示，等误差直线逼近法的特点是令各节点间非圆曲线与直线的误差?相等。其具体求解步骤如下： 以起点a( xa，ya)为圆心、?为半径作圆，确定允许误差的圆方程为 (x-xa)2+(y-ya)2=?2 (1) 圆与曲线的公切线PT的斜率为 K= yT-yp xT-xp (2) 式中的xT、yT、xp、yp需通过求解下列联立方程获得： { yT-yp=f1( xp)( xT-xp) yp=f1( xp) yT-yp=f2(xT)(xT-xp) yT=f2(xT) (3) 式中：f1(x)——误差圆函数 f2(x)——加工曲线函数 可知与PT 平行的弦ab的斜率为K，则弦ab的直线方程为 y-ya=K(x-xa) (4) 联立曲线方程和弦ab方程，可求得b点坐标为 { y=f2(x) y-ya=k(x-xa) (5) 重复上述步骤即可顺次求得c、d、e 等各点坐标。 3 等误差直线逼近法的节点计算 在数控加工中，刀尖轨迹的理论曲线一般取为抛物线y=ax2( a0，x0)，则有y=2ax。根据允差圆方程式(1)可得 { y=ya-[?2-(x-xa)2]? y=- x-xa y-ya (6) 因此方程组(3)可改写为 { yT-yp=- xp-xT (xT-xp) yp-yT yp=ya-[?2-(xp-xT)2]? yT-yp=2axT(xT-xp) yT=axT2 (7) 联立方程组(7)可得 4au3-4au2ya-t3+4aut3+4autxa=0 (8) 式中：t=xp-xa u=(?2-t2)? 图2 程序流程图 由已知条件可知：0≤t≤?。为求解t值，可在计算机上采用逐步有哪些信誉好的足球投注网站法或二分法等数值解法求根(本文选用逐步有哪些信誉好的足球投注网站法)。由式(8)解出t 值后即可求解方程。 根据式(5)可得 { y=ax2 y-ya=k(x-xa) (9) 式中k=(yT-yp)/(xT-xp) 求解方程组(9)，可得第一个节点坐标为 x= k+del 2a (10) 式中del=[k2-4a( kxa-ya)]? 然后以该节点作为新的起点，反复求解，即可顺次求得所有节点。 4 计算程序的编制 在AutoCAD内嵌的VBA(Visual Basic for Applications)开发环境下开发计算程序，程序流程如图2所示。 数控加工刀尖轨迹的理想曲线为抛物线。按?=0.05求出的节点绘制的折线段模拟曲线轨迹如图3所示。计算出的节点数据保存于“c:/data.txt”文件中，见图4。 图3 模拟曲线轨迹(?=0.05) 图4 保存节点坐标数据的文件(?=0.01) 求解方程的子程序如下： //求解方程 For tana=delta To 0 Step -0.0001 u=Sq(r delta ^ 2-tana ^ 2) texpr=4 * a * u ^ 3 + 4 * a * u * tana ^ 2 - 4 * a * u ^ 2 * ya-tana ^ 2 + 4 * a * u * tana * xa If texpr 0.00001 And texpr 0 Then Exit