基于时间分割法的两种抛物线插补算法的比较研究及n分步长算法.docVIP

基于时间分割法的两种抛物线插补算法的比较研究及n分步长算法 　　摘要　对基于时间分割法的两种抛物线插补算法在插补精度和插补运算时间两方面进行了比较分析，通过实验证明了数学分析的正确性，并提出了n分步长插补算法，该算法适用于多种曲线，可以有效地提高插补精度。 　　关键词：插补算法　时间分割法　n分步长 1　引言 　　对于任何一个复杂的零件，其轮廓线都可以由一些简单的几何因素组成，如直线、圆弧、抛物线等。即使一些自由型曲线、曲面体也可以由样条曲线进行拟合，最终还是要用直线、圆弧或抛物线进行二次插值或逼近，实现其数控加工，即完成插补运算的工作。插补的算法很多，但共同特点是必须在保证加工精度要求的前提下确定最大的插补间距，并计算出每一插补段的起止点在被加工零件轮廓线上的坐标值。时间分割法是根据系统设定的采样周期及设定的运动轨迹切向速度，确定每步所走距离，再分别计算每个单元时间两轴进给量，时间分割法广泛应用于现代数控机床中级以上的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