CN101109948A - 一种非均匀有理b样条曲线插补与位置跟踪方法 - Google Patents
一种非均匀有理b样条曲线插补与位置跟踪方法 Download PDFInfo
- Publication number
- CN101109948A CN101109948A CNA2007100295154A CN200710029515A CN101109948A CN 101109948 A CN101109948 A CN 101109948A CN A2007100295154 A CNA2007100295154 A CN A2007100295154A CN 200710029515 A CN200710029515 A CN 200710029515A CN 101109948 A CN101109948 A CN 101109948A
- Authority
- CN
- China
- Prior art keywords
- interpolation
- nurbs
- calculation
- buffer zone
- plugging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Numerical Control (AREA)
Abstract
本发明涉及一种非均匀有理B样条曲线(Non-Uniform Rational B-Spline,简称NURBS)插补与位置跟踪方法,其特点是NURBS插补计算采用预估与迭代算法,以提高其插补计算的速度;采用缓冲PID位置跟踪方法消除上述插补计算不均匀的影响;采用本发明方法可实现NURBS插补的快速计算,且计算过程稳定,速度控制精度高,可以用于实现中高档数控***的NURBS曲线插补。
Description
技术领域
本发明涉及计算机控制领域,特别是数控技术领域的一种非均匀有理B样条曲线插补与位置跟踪方法。
背景技术
目前,NURBS参数曲线是广为认可的CAD/CAM(Computer Aided Design andComputer aided Manufacture,计算机辅助设计与计算机辅助制造)***的产品形状表达曲线,NURBS插补功能一方面可以减小复杂曲线的程序代码量,另一方面可以保证加工表面光滑,但是,市面上现有的数控***中具有NURBS曲线插补功能的数控***产品不多,现行公开的NURBS插补方法多为理论计算方法,缺乏实用性。
发明内容
本发明的目的在于提供一种具有实用性的非均匀有理B样条曲线(NURBS)插补与位置跟踪方法,应用到数控***中,其控制器可以采用嵌入式***,也可以采用PC(Personal Computer,个人计算机)或IPC(Industrial Personal Computer,工业个人计算机),本发明方法均允许在上述硬件***上实现。
本发明的NURBS曲线插补与位置跟踪方案包括三个部分,一阶预估迭代插补方法;插补、位置跟踪、***通信与监控多任务调度方案;缓冲PID位置跟踪方法。
本发明采用一阶预估迭代方法计算曲线上的坐标点进行NURBS曲线插补;
采用缓冲比例-积分-微分(Proportional-Integral-Derivative,简称PID)位置跟踪方法消除上述插补算法的计算时延不均匀性,并跟踪指令位置。
将NURBS插补计算任务放在后台程序中执行,在控制器中开设缓冲区队列,NURBS插补计算的结果暂存于缓冲区中队列中;
位置跟踪PID算法放在前台程序中,它从缓冲区中读取插补结果,用于跟踪指令位置;
缓冲区队列用于平衡NURBS插补计算的时延不均匀。
技术原理:
(1)一阶预估插补方法
NURBS参数曲线的计算十分复杂,公认的计算方法为为德布尔方法,已知参数ui来计算NURBS曲线上的坐标点。但ui的计算方法有两类,一类是直接计算方法,如一阶泰勒方法、二阶泰勒方法、ADAMS方法等,上述算法计算耗时、特殊情况下不能求解,计算稳定性差。另一类是预估与迭代方法。它先用预估公式计算ui,然后再用德布尔方法进行坐标点验算,不满足设定的速度精度条件时,修改ui的值,再进行坐标点验算,直到满足条件为止;其计算流程如附图1所示。现行的预估公式为二阶预估公式:
当参数ui的进给步长过小时,ui+1<ui,这就造成插补运动方向折反。本发明采用一阶预估公式:
这就消除了上述插补运动方向折反的现象。
(2)前后台多任务调度
NURBS曲线插补与跟踪算法的任务划分为三个子任务,NURBS插补计算、位置跟踪PID控制、通信及***监控。其多任务调度方案如图2所示。
常规的数控***任务调度方案是前后台工作方式,NURBS插补计算、位置跟踪PID计算、通信及***监控三个任务均为强实时任务,他们被安排在前台程序中实时处理。由于本发明的NURBS插补算法采用预估与迭代方法,每一次的计算时延不均匀,该任务被安排到后台程序处理,所有的***空闲时间都被用于NURBS插补计算;位置跟踪PID计算任务和通信及***监控任务安排在前台程序中,位置跟踪PID计算任务每个采样周期调用一次,通信及***监控任务根据情况,若干个采样周期调用一次。这种多任务调度方案是本发明的特色之一。
(3)缓冲位置跟踪PID方法
传统的跟踪方法为计算出NURBS曲线坐标点后,即对相应的坐标轴进行位置跟踪PID计算。由于本发明的NURBS插补算法采用预估与迭代方法,计算时延不均匀,插补算法的中间结果预存在一个缓冲队列中,位置跟踪PID算法响应实时时钟Ts的信号,从缓冲区队列中取得插补位置指令,进行PID跟踪控制。本方法的计算任务配置如图3所示。
缓冲区队列的管理选用了队列上、下限两个指针,缓冲队列的长度大于等于队列上限指针时,停止插补计算,当队列长度指针小于队列下限指针时,开始插补计算。这就避免了插补计算任务的频繁切换。
本发明与现有技术相比具有如下优点和效果:
(1)提高了程序计算的稳定性
当NURBS曲线的运动控制为速度控制时,本发明方法不需要计算曲线曲率等,算法的稳定性高。
(2)提高了计算速度
一阶泰勒方法、二阶泰勒方法、ADAMS方法计算耗时,采用一阶预估与迭代方法平均计算时延比上述方法减少20%-30%。
附图说明
图1是本发明的插补计算流程框图;
图2是本发明方法的多任务调度流程图;
图3是本发明的缓冲位置跟踪流程图。
具体实施方式
如附图3所示,插补计算任务和位置跟踪控制任务分开进行,中间选用一个缓冲队列存放插补数据;插补计算任务提前进行,以保证缓冲区中有一定量的数据供实时的位置跟踪PID算法使用。
如附图2所示,插补任务放在背景程序中运行,位置跟踪PID程序和通信及***监控程序放在前台程序中运行,位置跟踪PID程序响应***的采样周期时钟(Ts),完成位置跟踪任务后即退出运行;通信及***监控程序在n次采样周期Ts后运动一次。***的空闲时间由插补程序计算曲线上的参数点坐标;中间的缓冲队列平衡插补计算的时延不均匀性,以防止位置跟踪任务中断。
如附图1所示,NURBS插补计算采用一阶预估公式,首先估计参数ui,然后计算速度v,如果速度不满足精度要求,按流程图中的公式修改估计值,再计算速度,判断速度精度,直到速度精度达到精度要求为止。
Claims (1)
1.一种非均匀有理B样条曲线(NURBS)插补与位置跟踪方法,其特征在于采用一阶预估迭代方法计算曲线上的坐标点,以进行NURBS曲线插补;
采用缓冲比例-积分-微分(Proportional-Integral-Derivative,简称PID)位置跟踪方法消除上述插补算法的计算时延不均匀影响,并跟踪指令位置;
将NURBS插补计算任务放在后台程序中执行,在控制器中开设缓冲区队列,NURBS插补计算的结果暂存于缓冲区中队列中;
位置跟踪PID算法放在前台程序中,它从缓冲区中读取插补结果,用于跟踪指令位置;
缓冲区队列用于平衡NURBS插补计算的时延不均匀。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007100295154A CN101109948A (zh) | 2007-08-01 | 2007-08-01 | 一种非均匀有理b样条曲线插补与位置跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007100295154A CN101109948A (zh) | 2007-08-01 | 2007-08-01 | 一种非均匀有理b样条曲线插补与位置跟踪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101109948A true CN101109948A (zh) | 2008-01-23 |
Family
ID=39042047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007100295154A Pending CN101109948A (zh) | 2007-08-01 | 2007-08-01 | 一种非均匀有理b样条曲线插补与位置跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101109948A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298358A (zh) * | 2011-04-30 | 2011-12-28 | 上海交通大学 | 五轴数控加工双nurbs刀具轨迹速度规划方法 |
CN102591257A (zh) * | 2012-02-27 | 2012-07-18 | 山东理工大学 | 面向参数曲线刀具轨迹的数控***轮廓误差控制方法 |
CN102794767A (zh) * | 2012-08-31 | 2012-11-28 | 江南大学 | 视觉引导的机器人关节空间b样条轨迹规划方法 |
CN103345201A (zh) * | 2013-07-01 | 2013-10-09 | 浙江工业大学 | Nurbs曲线降阶插补的运动控制*** |
CN104133422A (zh) * | 2014-07-16 | 2014-11-05 | 北京航空航天大学 | 一种非均匀有理b样条实时插补过程中的快速求值求导方法 |
CN104238458A (zh) * | 2013-06-08 | 2014-12-24 | 沈阳高精数控技术有限公司 | 一种nurbs曲线高速前瞻插补方法 |
CN104133422B (zh) * | 2014-07-16 | 2017-01-04 | 北京航空航天大学 | 一种非均匀有理b样条实时插补过程中的快速求值求导方法 |
WO2018205248A1 (zh) * | 2017-05-12 | 2018-11-15 | 深圳配天智能技术研究院有限公司 | 样条路径的插补方法及相关设备 |
CN114296398A (zh) * | 2021-11-16 | 2022-04-08 | 中南大学 | 一种用于激光切割的高速高精度插补方法 |
-
2007
- 2007-08-01 CN CNA2007100295154A patent/CN101109948A/zh active Pending
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298358B (zh) * | 2011-04-30 | 2013-05-01 | 上海交通大学 | 五轴数控加工双nurbs刀具轨迹速度规划方法 |
CN102298358A (zh) * | 2011-04-30 | 2011-12-28 | 上海交通大学 | 五轴数控加工双nurbs刀具轨迹速度规划方法 |
CN102591257A (zh) * | 2012-02-27 | 2012-07-18 | 山东理工大学 | 面向参数曲线刀具轨迹的数控***轮廓误差控制方法 |
CN102591257B (zh) * | 2012-02-27 | 2013-09-18 | 山东理工大学 | 面向参数曲线刀具轨迹的数控***轮廓误差控制方法 |
CN102794767A (zh) * | 2012-08-31 | 2012-11-28 | 江南大学 | 视觉引导的机器人关节空间b样条轨迹规划方法 |
CN104238458A (zh) * | 2013-06-08 | 2014-12-24 | 沈阳高精数控技术有限公司 | 一种nurbs曲线高速前瞻插补方法 |
CN103345201A (zh) * | 2013-07-01 | 2013-10-09 | 浙江工业大学 | Nurbs曲线降阶插补的运动控制*** |
CN104133422A (zh) * | 2014-07-16 | 2014-11-05 | 北京航空航天大学 | 一种非均匀有理b样条实时插补过程中的快速求值求导方法 |
CN104133422B (zh) * | 2014-07-16 | 2017-01-04 | 北京航空航天大学 | 一种非均匀有理b样条实时插补过程中的快速求值求导方法 |
WO2018205248A1 (zh) * | 2017-05-12 | 2018-11-15 | 深圳配天智能技术研究院有限公司 | 样条路径的插补方法及相关设备 |
CN109564417A (zh) * | 2017-05-12 | 2019-04-02 | 深圳配天智能技术研究院有限公司 | 样条路径的插补方法及相关设备 |
CN109564417B (zh) * | 2017-05-12 | 2021-07-09 | 深圳配天智能技术研究院有限公司 | 样条路径的插补方法及相关设备 |
CN114296398A (zh) * | 2021-11-16 | 2022-04-08 | 中南大学 | 一种用于激光切割的高速高精度插补方法 |
CN114296398B (zh) * | 2021-11-16 | 2024-04-05 | 中南大学 | 一种用于激光切割的高速高精度插补方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yau et al. | Fast Bezier interpolator with real-time lookahead function for high-accuracy machining | |
CN101109948A (zh) | 一种非均匀有理b样条曲线插补与位置跟踪方法 | |
CN101493687B (zh) | 实时前瞻全程加减速控制的nurbs曲线自适应分段插补方法 | |
Erkorkmaz et al. | Feedrate optimization for freeform milling considering constraints from the feed drive system and process mechanics | |
CN102681487B (zh) | 数控***中操作设备的轨迹平滑方法和装置及数控机床 | |
CN109976262B (zh) | 一种针对微线段加工的全局曲率连续光顺方法 | |
Hu et al. | Five-axis tool path generation based on machine-dependent potential field | |
WO2019119268A1 (zh) | 数控***、具有存储功能的装置及nurbs曲线插补方法 | |
Zhang et al. | Curve fitting and optimal interpolation on CNC machines based on quadratic B-splines | |
Du et al. | An error-bounded B-spline curve approximation scheme using dominant points for CNC interpolation of micro-line toolpath | |
CN105843174B (zh) | 一种样条曲线插补参数计算方法 | |
CN101791801A (zh) | 工业机器人运动规划与性能测试***及其实现方法 | |
CN102722140A (zh) | 基于s曲线加减速控制的多周期拐角小直线段插补方法 | |
CN103809521A (zh) | 基于弦截法的样条曲线插补方法 | |
CN110244658B (zh) | 一种基于改进bp神经网络提高裁床插补位置精度的方法 | |
CN201856249U (zh) | 工业机器人运动规划与性能测试*** | |
CN108170101A (zh) | 面向多项式样条曲线的插补方法及*** | |
CN108062071B (zh) | 参数曲线轨迹伺服轮廓误差的实时测定方法 | |
CN114925462A (zh) | 一种基于切削力与刚度关联演变的薄壁件加工变形预测方法 | |
CN104865898A (zh) | 数控机床中处理微小路径段的方法及设备 | |
WO2024124794A1 (zh) | 一种五轴联动同步刀具路径插补方法和*** | |
CN107486860A (zh) | 一种工程机械机械臂的控制方法 | |
CN110045681B (zh) | 一种基于位置预测的数控机床位置相关型误差的外置补偿方法 | |
CN104020718A (zh) | 一种nurbs曲线参考模型自适应实时插补方法 | |
CN110154020A (zh) | 一种自适应动态规划的柔性机械臂控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |