基于时间分割法的两种抛物线插补算法的比较研究及n分步长算法
摘要 对基于时间分割法的两种抛物线插补算法在插补精度和插补运算时间两方面进行了比较分析，通过实验证明了数学分析的正确性，并提出了n分步长插补算法，该算法适用于多种曲线，可以有效地提高插补精度。
关键词：插补算法 时间分割法 n分步长1 引言
对于任何一个复杂的零件，其轮廓线都可以由一些简单的几何因素组成，如直线、圆弧、抛物线等。即使一些自由型曲线、曲面体也可以由样条曲线进行拟合，zui终还是要用直线、圆弧或抛物线进行二次插值或逼近，实现其数控加工，即完成插补运算的工作。插补的算法很多，但共同特点是必须在保证加工精度要求的前提下确定zui大的插补间距，并计算出每一插补段的起止点在被加工零件轮廓线上的坐标值。时间分割法是根据系统设定的采样周期及设定的运动轨迹切向速度，确定每步所走距离，再分别计算每个单元时间两轴进给量，时间分割法广泛应用于现代数控机床中级以上的cnc系统［1，2］。目前，基于时间分割法的直线和圆弧的插补算法已经被广泛研究，但是对于象抛物线一类的较复杂曲线的插补问题还少有人研究，致使缺乏有效的算法，造成编程复杂，插补精度低等问题，影响数控系统的实用性。32位机处理速度快，运算能力强，便于人们进一步研究复杂曲线的插补问题。
本文通过对基于时间分割法的两种抛物线插补算法的比较分析，从插补精度及插补运算时间两方面进一步研究了抛物线插补问题；并在此基础上，提出了n分步长插补算法，该算法充分发挥32位计算机的性能，同时可以应用于许多曲线插补，有效地提高了插补精度。
2 插补算法
下面首先给出基于时间分割法的两种抛物线插补算法。
2.1 插补算法一
对于抛物线方程y＝－ax2＋bx，（a＞0，b＞0）。f＝ft，f为进给量，f为进给速度，t为采样时间，a（xi，yi）为本次插补点，b（xi＋1，yi＋1）为下次插补点，δxi，δyi为x、y轴进给距离。
ki＝－2axi＋bi
xi＋1＝xi＋δxi
yi＋1＝－ax2i＋1＋bxi＋1
δyi＝yi＋1－yi (1)
2.2 插补算法二［1］
抛物线方程为y＝－ax2＋bx（a＞0，b＞0），a（xi－1，yi－1）、b（xi，yi）为两已知插补点，c（xi＋1，yi＋1）为下一插补点，f为进给量，f为进给速度，t为插补时间。在抛物线顶点左侧：
δyi-1＝δyi－1－2aδx2i－1
xi＋1＝xi＋δxi
yi＋1＝－ax2i＋1＋bxi＋1
δyi＝yi＋1－yi （2）
在抛物线顶点右侧：
δyi-1＝δyi－1－2aδx2i－1
xi＋1＝xi＋δxi
yi＋1＝－ax2i＋1＋bxi＋1
δyi＝yi＋1－yi （3）
3 两种插补算法比较分析
对于时间分割法插补曲线，关键是弦长代替弧长，而弧长应为合成插补进给量f＝f.t。每次插补实际所走弦长与f差值越小，则相应的插补精度就越高。对上面提出的两种算法作如下比较。设抛物线方程为y＝－ax2＋bx（a＞0，b＞0），考虑插补点在顶点左侧情况，见图1。
图1 两种插补算法比较分析图
图中a（xi－1，yi－1）、b（xi，yi）为两已知插补点，c为下一理想插补点，ae为过a点所作抛物线切线，其斜率k＝tgα1，α1为ae与x轴向夹角，连接ab并延长，交以b为圆心、f为半径的圆于c1点，a、c1均为圆上的点。过c1作y轴平行线，交抛物线于c′点，延长c1c′，交过b点所作x轴平行线于d1点。c2为过b点作抛物线切线交于圆上的点，β1为bc2于x轴向夹角。延长bc′交圆于c3点，β2为bc3与x轴向夹角，过c2、c3、c分别作y轴平行线，交bd1延长线于d2、d3、d点。
可以看出c2点为*种算法的插补点，其插补量δxi＝｜bd2｜，δyi′＝｜c2d2｜，c3为第二种算法的插补点，其插补量δyi″＝｜c3d3｜，而c为理想插补点，其插补量δxi＝｜bd｜，δyi＝｜cd｜。
下面求bc2与bc3的斜率。先看a点，从假设方程及a点相对顶点位置可知ae斜率必大于ab斜率，即kae＝tgα1＞kab＝tgα2。因a、b、c1三点共线，所以bc1斜率kbc1＝tgα2，分别考虑a、b点，对a点有：y′i－1＝－2axi－1＋b；对b点有：
yi′＝－2axi＋b所以有
yi′－y′i－1＝－2a（xi－xi－1）＝－2aδx′i－1
得： yi′＝y′i－1－2aδx′i－1＝kab－2aδx′i－1
而对于b点，有｜c′d1｜＝y″i－1－2ax″i－12
对上一次插补，可以假设，δy′i－1＝δy″i－1，δx′i－1＝δx″i－1，于是有：
tgβ1＝kbc2＝kae－2aδx′i－1
tgβ2＝kbc3＝kab－2aδx′i－1
由于kae＞kab，所以kbc2＞kbc3
现在考虑b点局部坐标系，则有b(0，0)，圆方程为x′2＋y′2＝f2。对于过b点且斜率为k的直线交于圆的交点坐标，可求得
则有x′c2＜x′c3。而在该范围内k＝－2ax＋b＞0，所以y′c2＜y′c3，也就是
由上述推导可看出，算法二的合成插补实际进给量比算法一的更逼近于合成插补理想进给量，也就是说，算法二的插补精度理论上比算法一的插补精度高。但通过比较两种算法基本思想，可以看出算法二比算法一复杂，即算法二实际运行应较慢。在66mhzpc-486机型上进行实验，两种算法在不同的理想合成插补距离下对每次插补的实际距离作了统计平均。表1列出了两种算法在不同条件下实际合成插补量统计平均值。
表1 两种抛物线插补算法实际合成插补量平均值
0.02 0.04 0.06 0.08
a=0.01,b=0 0.019956 0.039669 0.059607 0.079293
0.019955 0.039665 0.059589 0.079272
a=0.1,b=0 0.019987 0.039682 0.059744 0.079345
0.019984 0.039680 0.059742 0.079339
a=1,b=0 0.019992 0.039710 0.059851 0.079452
0.019986 0.039709 0.059748 0.079447
在实验中测出两种算法每次插补运算时间，算法一t1＝0．000013s；算法二t2＝0．000021s。由以上实验数据可以看出，实验结果同理论分析一致，即算法一的运算速度比算法二的运算速度快。但算法二的匀速性及插补精度比算法一要好。 4 n分步长的改进算法
随着32位微机的广泛使用，提供了高速的运算能力，便得cnc插补算法可以适当增加复杂程度，而达到提高插补精度的目的。为此提出了下面的n分步长的改进算法。
以二分步长说明该算法的基本思想，以上面介绍的抛物线插补算法一为例，见图2。
图2 二分步长改进算法示意图
设抛物线方程为y＝－ax2＋bx（a＞0，b＞0），a（xi，yi）为已知插补点，f为理论合成插补量，ac′为过a点的抛物线切线，a，c′，b′1三点共线，作ab′1＝f，c′，c，d′三点共线，该线平行于y轴，a、d′、d三点共线，平行于x轴。｜ac′｜＝f／2。c点为c′d′与抛物线的交点，cb″为过c点抛物线切线，b″d与抛物线交于b点，｜cb″｜＝f／2。二分步长基本思想与算法一相同，但二分步长每步算两次，每次以f／2为理想合成插补量。两次后算出实际插补点b。
依据前述对算法一、二的比较分析相同方法，可以得出结论，采用二分步长法比整步长法进行插补的精度要高。同样，二分步长可推广至三分步长、四分步长直至n分步长的n分数，相应的插补精度都会提高。但同时也可看出，从n分步长到n+1分步长，插补时间都会增加，通过在66mhz486机型上对算法一和算法二进行n-1到n-10的n分步长算法实验，得出以下的每步实际平均插补时间，见表2。
表2 n分步长插补时间 (ms)
n 1 2 3 4 5 6 7 8 9 10
算法一 0.013 0.026 0.037 0.046 0.058 0.069 0.079 0.091 0.100 0.112
算法二 0.021 0.037 0.051 0.068 0.086 0.101 0.116 0.134 0.149 0.164
为了提高插补精度，可参考表2中的数据，根据实际需要的插补运算时间，选择合适的n值。其中的插补算法可适用于其它曲线的时间分割法基本插补算法。图3给出n分步长插补算法流程图。
图3 n分步长插补算法流程图
——可采用直线、圆弧、抛物线插补算法公式
5 结论
通过本文对两种抛物线插补算法的比较分析，可知算法二比算法一的插补精度要高，但算法一的运算速度较快；同时，根据对n分步长插补算法的研究，可看出通过选择适当n，可以有效地提高插补精度。
NDIS电桥接头701955&701956横河YOKOGAWA操作使用
根据受力情况采用不同结构尺寸的压铆螺钉
透平油滤油机选用合适的胶粘剂要考虑4个因素
阿托斯柱塞泵
深井潜水泵安装与使用要求
基于时间分割法的两种抛物线插补算法的比较研究及n分步长算法
氢型阳离子交换树脂的离子交换容量与性能
邯郸钢厂高炉煤气干法脉冲袋式除尘器工艺特点
北京不锈钢金库门解决方案
上银导轨的模组化设计
CLFS-1A便携式氯、氟离子水质检测仪
美国posi-flate不锈钢蝶阀操作手册
空调木托使用产品种类
促使我国智能导热油流量计产业走在前列
贺德克温度传感器该如何安装使用
浪涌保护器的工作类型及其应用
菏泽轻集料混凝土在建筑垫层中普遍使用了
冠层分析仪促进了农业的发展
如何选择实验室测试台
cems伴热采样管烟气防凝取样
关键词：插补算法 时间分割法 n分步长1 引言
对于任何一个复杂的零件，其轮廓线都可以由一些简单的几何因素组成，如直线、圆弧、抛物线等。即使一些自由型曲线、曲面体也可以由样条曲线进行拟合，zui终还是要用直线、圆弧或抛物线进行二次插值或逼近，实现其数控加工，即完成插补运算的工作。插补的算法很多，但共同特点是必须在保证加工精度要求的前提下确定zui大的插补间距，并计算出每一插补段的起止点在被加工零件轮廓线上的坐标值。时间分割法是根据系统设定的采样周期及设定的运动轨迹切向速度，确定每步所走距离，再分别计算每个单元时间两轴进给量，时间分割法广泛应用于现代数控机床中级以上的cnc系统［1，2］。目前，基于时间分割法的直线和圆弧的插补算法已经被广泛研究，但是对于象抛物线一类的较复杂曲线的插补问题还少有人研究，致使缺乏有效的算法，造成编程复杂，插补精度低等问题，影响数控系统的实用性。32位机处理速度快，运算能力强，便于人们进一步研究复杂曲线的插补问题。
本文通过对基于时间分割法的两种抛物线插补算法的比较分析，从插补精度及插补运算时间两方面进一步研究了抛物线插补问题；并在此基础上，提出了n分步长插补算法，该算法充分发挥32位计算机的性能，同时可以应用于许多曲线插补，有效地提高了插补精度。
2 插补算法
下面首先给出基于时间分割法的两种抛物线插补算法。
2.1 插补算法一
对于抛物线方程y＝－ax2＋bx，（a＞0，b＞0）。f＝ft，f为进给量，f为进给速度，t为采样时间，a（xi，yi）为本次插补点，b（xi＋1，yi＋1）为下次插补点，δxi，δyi为x、y轴进给距离。
ki＝－2axi＋bi
xi＋1＝xi＋δxi
yi＋1＝－ax2i＋1＋bxi＋1
δyi＝yi＋1－yi (1)
2.2 插补算法二［1］
抛物线方程为y＝－ax2＋bx（a＞0，b＞0），a（xi－1，yi－1）、b（xi，yi）为两已知插补点，c（xi＋1，yi＋1）为下一插补点，f为进给量，f为进给速度，t为插补时间。在抛物线顶点左侧：
δyi-1＝δyi－1－2aδx2i－1
xi＋1＝xi＋δxi
yi＋1＝－ax2i＋1＋bxi＋1
δyi＝yi＋1－yi （2）
在抛物线顶点右侧：
δyi-1＝δyi－1－2aδx2i－1
xi＋1＝xi＋δxi
yi＋1＝－ax2i＋1＋bxi＋1
δyi＝yi＋1－yi （3）
3 两种插补算法比较分析
对于时间分割法插补曲线，关键是弦长代替弧长，而弧长应为合成插补进给量f＝f.t。每次插补实际所走弦长与f差值越小，则相应的插补精度就越高。对上面提出的两种算法作如下比较。设抛物线方程为y＝－ax2＋bx（a＞0，b＞0），考虑插补点在顶点左侧情况，见图1。
图1 两种插补算法比较分析图
图中a（xi－1，yi－1）、b（xi，yi）为两已知插补点，c为下一理想插补点，ae为过a点所作抛物线切线，其斜率k＝tgα1，α1为ae与x轴向夹角，连接ab并延长，交以b为圆心、f为半径的圆于c1点，a、c1均为圆上的点。过c1作y轴平行线，交抛物线于c′点，延长c1c′，交过b点所作x轴平行线于d1点。c2为过b点作抛物线切线交于圆上的点，β1为bc2于x轴向夹角。延长bc′交圆于c3点，β2为bc3与x轴向夹角，过c2、c3、c分别作y轴平行线，交bd1延长线于d2、d3、d点。
可以看出c2点为*种算法的插补点，其插补量δxi＝｜bd2｜，δyi′＝｜c2d2｜，c3为第二种算法的插补点，其插补量δyi″＝｜c3d3｜，而c为理想插补点，其插补量δxi＝｜bd｜，δyi＝｜cd｜。
下面求bc2与bc3的斜率。先看a点，从假设方程及a点相对顶点位置可知ae斜率必大于ab斜率，即kae＝tgα1＞kab＝tgα2。因a、b、c1三点共线，所以bc1斜率kbc1＝tgα2，分别考虑a、b点，对a点有：y′i－1＝－2axi－1＋b；对b点有：
yi′＝－2axi＋b所以有
yi′－y′i－1＝－2a（xi－xi－1）＝－2aδx′i－1
得： yi′＝y′i－1－2aδx′i－1＝kab－2aδx′i－1
而对于b点，有｜c′d1｜＝y″i－1－2ax″i－12
对上一次插补，可以假设，δy′i－1＝δy″i－1，δx′i－1＝δx″i－1，于是有：
tgβ1＝kbc2＝kae－2aδx′i－1
tgβ2＝kbc3＝kab－2aδx′i－1
由于kae＞kab，所以kbc2＞kbc3
现在考虑b点局部坐标系，则有b(0，0)，圆方程为x′2＋y′2＝f2。对于过b点且斜率为k的直线交于圆的交点坐标，可求得
则有x′c2＜x′c3。而在该范围内k＝－2ax＋b＞0，所以y′c2＜y′c3，也就是
由上述推导可看出，算法二的合成插补实际进给量比算法一的更逼近于合成插补理想进给量，也就是说，算法二的插补精度理论上比算法一的插补精度高。但通过比较两种算法基本思想，可以看出算法二比算法一复杂，即算法二实际运行应较慢。在66mhzpc-486机型上进行实验，两种算法在不同的理想合成插补距离下对每次插补的实际距离作了统计平均。表1列出了两种算法在不同条件下实际合成插补量统计平均值。
表1 两种抛物线插补算法实际合成插补量平均值
0.02 0.04 0.06 0.08
a=0.01,b=0 0.019956 0.039669 0.059607 0.079293
0.019955 0.039665 0.059589 0.079272
a=0.1,b=0 0.019987 0.039682 0.059744 0.079345
0.019984 0.039680 0.059742 0.079339
a=1,b=0 0.019992 0.039710 0.059851 0.079452
0.019986 0.039709 0.059748 0.079447
在实验中测出两种算法每次插补运算时间，算法一t1＝0．000013s；算法二t2＝0．000021s。由以上实验数据可以看出，实验结果同理论分析一致，即算法一的运算速度比算法二的运算速度快。但算法二的匀速性及插补精度比算法一要好。 4 n分步长的改进算法
随着32位微机的广泛使用，提供了高速的运算能力，便得cnc插补算法可以适当增加复杂程度，而达到提高插补精度的目的。为此提出了下面的n分步长的改进算法。
以二分步长说明该算法的基本思想，以上面介绍的抛物线插补算法一为例，见图2。
图2 二分步长改进算法示意图
设抛物线方程为y＝－ax2＋bx（a＞0，b＞0），a（xi，yi）为已知插补点，f为理论合成插补量，ac′为过a点的抛物线切线，a，c′，b′1三点共线，作ab′1＝f，c′，c，d′三点共线，该线平行于y轴，a、d′、d三点共线，平行于x轴。｜ac′｜＝f／2。c点为c′d′与抛物线的交点，cb″为过c点抛物线切线，b″d与抛物线交于b点，｜cb″｜＝f／2。二分步长基本思想与算法一相同，但二分步长每步算两次，每次以f／2为理想合成插补量。两次后算出实际插补点b。
依据前述对算法一、二的比较分析相同方法，可以得出结论，采用二分步长法比整步长法进行插补的精度要高。同样，二分步长可推广至三分步长、四分步长直至n分步长的n分数，相应的插补精度都会提高。但同时也可看出，从n分步长到n+1分步长，插补时间都会增加，通过在66mhz486机型上对算法一和算法二进行n-1到n-10的n分步长算法实验，得出以下的每步实际平均插补时间，见表2。
表2 n分步长插补时间 (ms)
n 1 2 3 4 5 6 7 8 9 10
算法一 0.013 0.026 0.037 0.046 0.058 0.069 0.079 0.091 0.100 0.112
算法二 0.021 0.037 0.051 0.068 0.086 0.101 0.116 0.134 0.149 0.164
为了提高插补精度，可参考表2中的数据，根据实际需要的插补运算时间，选择合适的n值。其中的插补算法可适用于其它曲线的时间分割法基本插补算法。图3给出n分步长插补算法流程图。
图3 n分步长插补算法流程图
——可采用直线、圆弧、抛物线插补算法公式
5 结论
通过本文对两种抛物线插补算法的比较分析，可知算法二比算法一的插补精度要高，但算法一的运算速度较快；同时，根据对n分步长插补算法的研究，可看出通过选择适当n，可以有效地提高插补精度。
NDIS电桥接头701955&701956横河YOKOGAWA操作使用
根据受力情况采用不同结构尺寸的压铆螺钉
透平油滤油机选用合适的胶粘剂要考虑4个因素
阿托斯柱塞泵
深井潜水泵安装与使用要求
基于时间分割法的两种抛物线插补算法的比较研究及n分步长算法
氢型阳离子交换树脂的离子交换容量与性能
邯郸钢厂高炉煤气干法脉冲袋式除尘器工艺特点
北京不锈钢金库门解决方案
上银导轨的模组化设计
CLFS-1A便携式氯、氟离子水质检测仪
美国posi-flate不锈钢蝶阀操作手册
空调木托使用产品种类
促使我国智能导热油流量计产业走在前列
贺德克温度传感器该如何安装使用
浪涌保护器的工作类型及其应用
菏泽轻集料混凝土在建筑垫层中普遍使用了
冠层分析仪促进了农业的发展
如何选择实验室测试台
cems伴热采样管烟气防凝取样