CN102073301A - 具有刀具长度补偿功能的五轴样条插补器 - Google Patents

具有刀具长度补偿功能的五轴样条插补器 Download PDF

Info

Publication number
CN102073301A
CN102073301A CN 201110027578 CN201110027578A CN102073301A CN 102073301 A CN102073301 A CN 102073301A CN 201110027578 CN201110027578 CN 201110027578 CN 201110027578 A CN201110027578 A CN 201110027578A CN 102073301 A CN102073301 A CN 102073301A
Authority
CN
China
Prior art keywords
cutter
coordinate
module
batten
interpolation cycle
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.)
Granted
Application number
CN 201110027578
Other languages
English (en)
Other versions
CN102073301B (zh
Inventor
陈良骥
李慧莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN2011100275782A priority Critical patent/CN102073301B/zh
Publication of CN102073301A publication Critical patent/CN102073301A/zh
Application granted granted Critical
Publication of CN102073301B publication Critical patent/CN102073301B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

具有刀具长度补偿功能的五轴样条插补器涉及样条技术在五轴数控机床及其数控系统中的具体应用,以解决现行技术普遍采用的刀具轨迹曲线被离散成直线段的“以直代曲”逼近加工模式存在的精度问题。本发明的语法词法分析模块负责读取数控程序、语法词法检查、存储控制点坐标等数据信息;刀具长度补偿模块负责计算刀具长度改变后新的控制点坐标、更新存储的控制点坐标数据;加减速处理模块负责按照梯形规律完成插补过程的加速、减速控制;刀位生成模块负责计算每个插补周期刀心点的坐标和刀轴单位矢量;刀位后处理模块根据不同的机床结构形式将实时生成的刀位数据转换为五个伺服电机的运动控制指令。

Description

具有刀具长度补偿功能的五轴样条插补器
技术领域
本发明涉及一种五轴数控机床数字控制系统的具有样条直接生成功能的插补器,具体涉及的技术领域是非均匀有理B样条(Non-Uniform Rational B Spline,NURBS)曲线在五轴数控机床及其数控系统中的实现。
背景技术
利用五轴数控机床加工复杂曲面零部件的传统技术是:先借助计算机辅助设计(Computer-AidedDesig,CAD)技术绘制零部件加工表面,再借助计算机辅助制造(Computer-Aided Manufacturing,CAM)技术生成加工刀具走刀轨迹,并将这些刀具轨迹曲线在加工误差范围内离散成直线段,最后再由后置处理模块将这些直线段转换成五轴线性插补指令格式,由机床数控装置插补出各坐标轴的运动指令驱动伺服电机来完成零件曲面的加工。这种传统数控加工技术存在诸多不足之处,进而严重制约加工表面质量的提高。随着计算机等相关技术的发展和市场对加工质量要求的不断提高,客观上为一种新的数字控制器的产生提供了有力条件。这种新的技术发展趋势是:摒弃以直线段逼近复杂曲线加工轨迹的方式,不断增强数字控制系统的复杂曲线轨迹直接生成能力,使刀具在很长一段加工轨迹上可以较高的进给速度切削加工表面而且有较高的加工表面质量和加工效率。
发明内容
本发明针对传统五轴逼近式加工技术的不足之处,从数控系统的插补技术角度出发,为五轴数控机床的数控系统提供一种具有刀具长度补偿的样条插补器,现将本发明所包括的内容陈述如下。
插补器由语法词法分析模块1、刀具长度补偿模块2、加减速处理模块3、刀位生成模块4以及刀位后处理模块5五个部分组成,各部分具体包括如下内容,
语法词法分析模块1,读取含有五轴样条指令的数控加工程序、语法检查、词法检查、存储布尔变量TLC、控制点坐标、长度补偿值、节点矢量、样条阶次以及切削进给速度F等数据信息;发明的五轴样条指令格式定义如下:
N_SPLINE
P_F_TLC_H_
X_Y_Z_TX_TY_TZ_R_K_......
X_Y_Z_TX_TY_TZ_R_K_
K_......
K_EOL
(注释:N表示程序段段号,SPLINE表示样条插补开始,P表示样条阶次,F表示切削进给速度,TLC表示刀具长度补偿,H表示长度补偿寄存器号,X、Y、Z表示刀心点样条的控制点坐标,TX、TY、TZ表示方位样条的控制点坐标,R表示权重,K表示节点值,EOL为程序段结束符)
刀具长度补偿模块2,如果语法词法分析模块1中的布尔变量TLC=TRUE则按照以下过程进行刀具长度补偿:将语法词法分析模块1中寄存器H存储的刀具长度值或者长度补偿值、形成刀心点样条和方位样条的所有控制点坐标等数据信息经数学计算后得出补偿后新的控制点坐标、更新语法词法分析模块1中原来存储的控制点坐标数据;
加减速处理模块3,采用梯形加减速控制规律、根据用户事先设定的加速度值a和插补周期Ts应用数值分析方法预测出减速点的参数值ud和减速时插补周期个数kd、将ud和kd写入语法词法分析模块1的数据区、进给速度V关于刀位点参数u对应的插补周期序号k的函数关系:
Figure BSA00000426739300021
刀位生成模块4,结合加速度a、切削进给速度、插补周期Ts、长度补偿后的控制点坐标、样条阶次、节点矢量等,利用一阶泰勒方法将刀心点样条C0(u)展开后得出下一个插补周期刀位点参数值uk+1与当前插补周期刀位点参数值uk的递推表达式
u k + 1 = u k + T s V ( k ) | | dC 0 ( u k ) du | | ,
进而得出插补后的刀心点坐标C0(uk+1)和刀轴上一点的坐标C1(uk+1),计算刀轴方位矢量并进行单位化
T k + 1 = C 1 ( u k + 1 ) - C 0 ( u k + 1 ) | | C 1 ( u k + 1 ) - C 0 ( u k + 1 ) | | ;
刀位后处理模块5,结合刀位生成模块4得到的刀心点坐标C0(uk+1)和刀轴方位单位矢量Tk+1,将刀轴方位单位矢量Tk+1转化为五轴数控机床的两个转动坐标、利用克莱姆法则解出五轴数控机床的三个平动坐标。
附图说明
图1是本发明的整体结构示意图;图2是刀具长度补偿示意图;图3是切削进给速度示意图;图4是插补数据输出示意图。
具体实施方式
本实施方式实现的五轴样条插补软件包是基于微软公司的可视化编程工具Visual C++ 6.0开发的,首先定义如下数据结构:
struct CPt{double x;double y;doublez;};//空间点坐标
typedef CClassArray<CPt,CPt&>CPtArray;//CClassArray为自定义类
typelef CClassArray<double,double&>CDoubleArray;//双精度浮点数数组
class Spline{public://样条类
      int m_P,m_H;//样条阶次和长度补偿寄存器号
      double m_F;//切削进给速度
      bool m_TLC;//刀具长度补偿布尔变量
      CPtArray m_CtrlPtArr;//刀心点样条控制点
      CPtArray m_TCtrlPtArr;//方位样条控制点
      CDoubleArray m_RArr;//权重
      CDoubleArray m_KArr;};//节点矢量
由语法词法分析模块1,读取含有五轴样条指令的数控加工程序并进行语法词法检查后,将与样条插补相关的所有数据信息按上述数据结构存储。
如果布尔变量TLC=TRUE则进行以下刀具长度补偿,如果TLC=FALSE则跳过。
刀具长度补偿的方向为由刀心点指向刀轴上一点的方向,参见图2。
图2中刀具长度为l,刀心点坐标及刀轴方位角分别为(x,y,z,a,c),首先计算刀轴方位单位矢量T。初始状态下,刀具竖直向下且平行于机床坐标系的Z轴,即T0={0,0,1}。刀具分别绕X轴和Z轴旋转a角和c角后刀轴单位矢量为T,由坐标变换原理有
T = cos c - sin c 0 sin c cos c 0 0 0 1 1 0 0 0 cos a - sin a 0 sin a cos a 0 0 1
= sin a sin c - sin a cos c cos a ,
进而可以计算出刀具长度补偿后新控制点坐标
x s = x + l sin a sin c y s = y - l sin a cos c . z s = z + l cos a
对所有刀心点和刀轴上一点重复进行上述计算过程,将计算得到的新的控制点坐标替换语法词法分析模块1中原来存储的控制点坐标数据。
按照下述方法计算减速点的参数值ud
用户指定的加速度为a,加速结束时的目标速度为F,参见图3。设
f ( u ) = &Integral; u l | | dC 0 ( u ) du | | du - F 2 2 a ,
则有f(ud)=0;设ud的初始值为u0=1,此时
Figure BSA00000426739300035
则有如下迭代计算式
ui+1=ui-f(ui)/f′(ui)   i≥0,
停止迭代计算的条件为(设迭代误差为η)
| u i + 1 - u i | u i < &eta; ui≠0。
在刀位生成模块4,结合加速度a、切削进给速度、插补周期Ts、长度补偿后的控制点坐标、样条阶次、节点矢量等,利用一阶泰勒方法将刀心点样条C0(u)展开后得出下一个插补周期刀位点参数值uk+1与当前插补周期刀位点参数值uk的递推表达式
u k + 1 = u k + T s V ( k ) | | d C 0 ( u k ) du | | ,
上式中
Figure BSA00000426739300041
进而得出插补后的刀心点坐标C0(uk+1)和刀轴上一点的坐标C1(uk+1),计算刀轴方位矢量并进行单位化
T k + 1 = C 1 ( u k + 1 ) - C 0 ( u k + 1 ) | | C 1 ( u k + 1 ) - C 0 ( u k + 1 ) | | .
提取出下一个插补周期刀心点坐标C0(uk+1)的数据信息和刀轴方位单位矢量Tk+1的数据信息分别定义为(ox,oy,oz)和(tx,ty,tz),同时提取五轴数控机床结构参数信息a和b。
用如下步骤计算两个转动坐标Ak+1、Ck+1
(1)当tz=1时,tx=ty=0,有
Xk+1=ox,Yk+1=oy,Zk+1=oz-a,Ak+1=0,Ck+1=0;
(2)当tz≠1时,sin Ak+1≠0,有
Ak+1=arccostz,Ak+1∈[-π/2,0)∪(0,π/2]
Ck+1=arctan(-tx/ty)。
采用克莱姆法则计算三个平动坐标Xk+1、Yk+1、Zk+1
X k + 1 = d 1 b 1 c 1 d 2 b 2 c 2 d 3 b 3 c 3 / a 1 b 1 c 1 a 2 b 2 c 2 a 3 b 3 c 3 ,
Y k + 1 = a 1 d 1 c 1 a 2 d 2 c 2 a 3 d 3 c 3 / a 1 b 1 c 1 a 2 b 2 c 2 a 3 b 3 c 3 ,
Z k + 1 = a 1 b 1 d 1 a 2 b 2 d 2 a 3 b 3 d 3 / a 1 b 1 c 1 a 2 b 2 c 2 a 3 b 3 c 3 ,
以上三式中
a 1 = cos C k + 1 b 1 = - cos A k + 1 sin C k + 1 c 1 = t x d 1 = o x - ( a + b ) t x a 2 = sin C k + 1 b 2 = cos A k + 1 cos C k + 1 c 2 = t y d 2 = o y - ( a + b ) t y a 3 = 0 b 3 = sin A k + 1 c 3 = t z d 3 = o z + b - ( a + b ) t z .
将得到的下一插补周期的三个平动坐标和两个转动坐标与当前插补周期相应的坐标作差后计算机床运动坐标增量(ΔXk+1,ΔYk+1,ΔZk+1,ΔAk+1,ΔCk+1),将该增量坐标采用数据通信的方式输入给五轴数控机床的数字式伺服驱动器以完成各运动坐标轴的进给伺服,参见图4。

Claims (2)

1.具有刀具长度补偿功能的五轴样条插补器,它由以下五个功能模块组成:语法词法分析模块(1)、刀具长度补偿模块(2)、加减速处理模块(3)、刀位生成模块(4)以及刀位后处理模块(5);
其技术特征在于:
语法词法分析模块(1),读取含有五轴样条指令的数控加工程序、语法检查、词法检查、存储布尔变量TLC、控制点坐标、长度补偿值、节点矢量、样条阶次以及切削进给速度F等数据信息;
刀具长度补偿模块(2),如果语法词法分析模块(1)中的布尔变量TLC=TRUE则按照以下过程进行刀具长度补偿:将语法词法分析模块(1)中寄存器H存储的刀具长度值或者长度补偿值、形成刀心点样条和方位样条的所有控制点坐标等数据信息经数学计算后得出补偿后新的控制点坐标、更新语法词法分析模块(1)中原来存储的控制点坐标数据;
加减速处理模块(3),采用梯形加减速控制规律、根据用户事先设定的加速度值a和插补周期Ts应用数值分析方法预测出减速点的参数值ud和减速时插补周期个数kd、将ud和kd写入语法词法分析模块(1)的数据区、进给速度V关于刀位点参数u对应的插补周期序号k的函数关系;
Figure FSA00000426739200011
刀位生成模块(4),结合加速度a、切削进给速度、插补周期Ts、长度补偿后的控制点坐标、样条阶次、节点矢量等,利用一阶泰勒方法将刀心点样条C0(u)展开后得出下一个插补周期刀位点参数值uk+1与当前插补周期刀位点参数值uk的递推表达式
u k + 1 = u k + T s V ( k ) | | dC 0 ( u k ) du | | ,
进而得出插补后的刀心点坐标C0(uk+1)和刀轴上一点的坐标C1(uk+1),计算刀轴方位矢量并进行单位化
T k + 1 = C 1 ( u k + 1 ) - C 0 ( u k + 1 ) | | C 1 ( u k + 1 ) - C 0 ( u k + 1 ) | | ;
刀位后处理模块(5),结合刀位生成模块(4)得到的刀心点坐标C0(uk+1)和刀轴方位单位矢量Tk+1,将刀轴方位单位矢量Tk+1转化为五轴数控机床的两个转动坐标、利用克莱姆法则解出五轴数控机床的三个平动坐标。
2.根据权利要求1所述的具有刀具长度补偿功能的五轴样条插补器,得到的下一插补周期的三个平动坐标和两个转动坐标与当前插补周期相应的坐标作差后计算机床运动坐标增量,将该增量坐标采用数据通信的方式输入给五轴数控机床的数字式伺服驱动器以完成各运动坐标轴的进给伺服。
CN2011100275782A 2011-01-21 2011-01-21 具有刀具长度补偿功能的五轴样条插补器 Expired - Fee Related CN102073301B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100275782A CN102073301B (zh) 2011-01-21 2011-01-21 具有刀具长度补偿功能的五轴样条插补器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100275782A CN102073301B (zh) 2011-01-21 2011-01-21 具有刀具长度补偿功能的五轴样条插补器

Publications (2)

Publication Number Publication Date
CN102073301A true CN102073301A (zh) 2011-05-25
CN102073301B CN102073301B (zh) 2012-07-25

Family

ID=44031874

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100275782A Expired - Fee Related CN102073301B (zh) 2011-01-21 2011-01-21 具有刀具长度补偿功能的五轴样条插补器

Country Status (1)

Country Link
CN (1) CN102073301B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102393678A (zh) * 2011-08-25 2012-03-28 沈阳高精数控技术有限公司 适用于五轴数控装置的轨迹平滑处理方法
CN103324141A (zh) * 2013-06-14 2013-09-25 浙江工业大学 一种高精度变插补周期的多轴联动运动控制方法
CN104020717A (zh) * 2014-05-28 2014-09-03 上海工程技术大学 基于参数同步的五轴等距双nurbs刀具路径插补方法
CN104850059A (zh) * 2014-02-13 2015-08-19 发那科株式会社 防止加工条件变更后的误加工的数值控制装置
CN109933009A (zh) * 2019-04-15 2019-06-25 天津工业大学 一种基于刀触点路径段的五轴线性插补方法
CN112346411A (zh) * 2020-11-23 2021-02-09 清华大学 使用双nurbs刀具轨迹叶轮数控加工速度控制方法及装置
CN113031525A (zh) * 2021-03-03 2021-06-25 福州大学 一种应用于数控加工的多项式加减速运动控制方法及设备
CN115599043A (zh) * 2022-12-12 2023-01-13 北京精雕科技集团有限公司(Cn) 用于提升五轴机床加工精度的方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5723961A (en) * 1995-07-17 1998-03-03 Mitsubishi Denki Kabushiki Kaisha Numerical control apparatus having spline interpolating function
CN101497140A (zh) * 2009-02-26 2009-08-05 上海交通大学 五轴数控侧铣加工进给率离线规划方法
CN101738984A (zh) * 2009-12-23 2010-06-16 哈尔滨工业大学 一种基于四元数的五坐标样条插补控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5723961A (en) * 1995-07-17 1998-03-03 Mitsubishi Denki Kabushiki Kaisha Numerical control apparatus having spline interpolating function
CN101497140A (zh) * 2009-02-26 2009-08-05 上海交通大学 五轴数控侧铣加工进给率离线规划方法
CN101738984A (zh) * 2009-12-23 2010-06-16 哈尔滨工业大学 一种基于四元数的五坐标样条插补控制方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102393678A (zh) * 2011-08-25 2012-03-28 沈阳高精数控技术有限公司 适用于五轴数控装置的轨迹平滑处理方法
CN102393678B (zh) * 2011-08-25 2014-05-07 沈阳高精数控技术有限公司 适用于五轴数控装置的轨迹平滑处理方法
CN103324141A (zh) * 2013-06-14 2013-09-25 浙江工业大学 一种高精度变插补周期的多轴联动运动控制方法
CN103324141B (zh) * 2013-06-14 2015-04-29 浙江工业大学 一种高精度变插补周期的多轴联动运动控制方法
US9715228B2 (en) 2014-02-13 2017-07-25 Fanuc Corporation Numerical controller capable of preventing wrong machining after machining condition change
CN104850059A (zh) * 2014-02-13 2015-08-19 发那科株式会社 防止加工条件变更后的误加工的数值控制装置
CN104020717B (zh) * 2014-05-28 2016-08-24 上海工程技术大学 基于参数同步的五轴等距双nurbs刀具路径插补方法
CN104020717A (zh) * 2014-05-28 2014-09-03 上海工程技术大学 基于参数同步的五轴等距双nurbs刀具路径插补方法
CN109933009A (zh) * 2019-04-15 2019-06-25 天津工业大学 一种基于刀触点路径段的五轴线性插补方法
CN112346411A (zh) * 2020-11-23 2021-02-09 清华大学 使用双nurbs刀具轨迹叶轮数控加工速度控制方法及装置
CN113031525A (zh) * 2021-03-03 2021-06-25 福州大学 一种应用于数控加工的多项式加减速运动控制方法及设备
CN113031525B (zh) * 2021-03-03 2022-06-21 福州大学 一种应用于数控加工的多项式加减速运动控制方法及设备
CN115599043A (zh) * 2022-12-12 2023-01-13 北京精雕科技集团有限公司(Cn) 用于提升五轴机床加工精度的方法及装置

Also Published As

Publication number Publication date
CN102073301B (zh) 2012-07-25

Similar Documents

Publication Publication Date Title
CN102073301B (zh) 具有刀具长度补偿功能的五轴样条插补器
CN102591257B (zh) 面向参数曲线刀具轨迹的数控系统轮廓误差控制方法
CN103064344B (zh) 一种基于nurbs曲线插补的速度平滑控制方法
CN102147600B (zh) 实时生成曲率连续路径的数控插补系统
CN104714477B (zh) 一种加工文件规划系统与加工文件的生成方法
CN103934528B (zh) 一种用于电火花加工的六轴联动插补方法
CN102023613A (zh) 一种五轴联动数控加工后置处理器及其处理方法
CN105425727A (zh) 五轴侧铣加工刀具路径光顺方法
CN102866671B (zh) 大圆弧直纹面数控加工刀具轨迹规划方法
CN103801981A (zh) 一种用于样条插补的四次多项式速度规划算法
CN104460515B (zh) 一种基于后置处理五轴刀具长度补偿方法
CN103176428B (zh) 基于球坐标的cnc系统插补算法及实现该算法的装置
CN104678888B (zh) 基于机床运动约束环形刀多轴加工复杂曲面刀轨生成方法
CN102069419B (zh) 一种主动误差抵消的自由曲面超精密车削方法
CN105234743A (zh) 一种五轴加工中心刀具变形误差补偿方法
CN101870073B (zh) 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法
CN106647623A (zh) 一种几何精度及衔接速度最优化的五轴联动平滑插补方法
CN102789199A (zh) 数控系统中多次误差加权叠加的补偿系统
CN106547251A (zh) 一种基于低速敏感区域插补信息反馈的五轴刀路生成方法
CN106202822A (zh) 面向叶片自适应加工的b样条曲面模型重构方法
CN113741342A (zh) 一种五轴联动轨迹误差溯源方法
CN202120088U (zh) 数控系统中多次误差加权叠加的补偿系统
Xu et al. An interpolation method based on adaptive smooth feedrate scheduling and parameter increment compensation for NURBS curve
Wei et al. Modeling and machining of integral impeller based on NURBS curve
CN102402198B (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
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120725

Termination date: 20130121

CF01 Termination of patent right due to non-payment of annual fee