CN102402198A - 一种多轴数控机床通用后置处理方法 - Google Patents

一种多轴数控机床通用后置处理方法 Download PDF

Info

Publication number
CN102402198A
CN102402198A CN2011103243887A CN201110324388A CN102402198A CN 102402198 A CN102402198 A CN 102402198A CN 2011103243887 A CN2011103243887 A CN 2011103243887A CN 201110324388 A CN201110324388 A CN 201110324388A CN 102402198 A CN102402198 A CN 102402198A
Authority
CN
China
Prior art keywords
machine tool
motion
cutter
cutter location
location
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
CN2011103243887A
Other languages
English (en)
Other versions
CN102402198B (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN 201110324388 priority Critical patent/CN102402198B/zh
Publication of CN102402198A publication Critical patent/CN102402198A/zh
Application granted granted Critical
Publication of CN102402198B publication Critical patent/CN102402198B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种多轴数控机床通用后置处理方法,包括如下步骤:(1)根据机床结构建立机床运动变换链(2)根据步骤(1)中所建立的机床运动变换链,建立从刀具到工件的变换矩阵QWT(3)输入初始刀位点CL1的机床各平动轴运动坐标Δ1(4)计算出所述初始刀位点CL1后续的任一刀位点CLi+1的平动轴运动坐标Δi+1(5)判断i的值是否小于n,如果小于n则返回步骤(4.1),将i增加1后继续计算,否则,计算结束,输出各刀位点对应的机床平动轴运动坐标和旋转轴运动坐标。本发明克服了一般后置处理方法中机床各轴运动坐标计算公式需要手工推导的缺点,可以满足各类多轴数控机床的后置处理需求,具有求解速度快、求解精度高的优点。

Description

一种多轴数控机床通用后置处理方法
技术领域:
本发明涉及多轴数控加工领域,具体是一种多轴数控机床通用后置处理方法。
背景技术:
随着新能源、航空航天、造船等行业的发展,对主要功能部件如风电叶片、飞机进气道、螺旋桨等的加工质量要求越来越高。这类零件外形及工作型面形状复杂,加工精度高,只有采用多轴数控机床才能完成加工。
在五轴数控加工中,后置处理的主要任务是将CAM软件生成的加工刀位轨迹源文件转换成为特定机床可接受的数控NC代码,目前多轴数控机床后置求解,主要包含手工推导和数值计算这两种方法。手工推导方法不仅过程复杂、容易出错,而且得到的计算公式只能用于特定结构的机床,不具备通用性;数值方法具有通用性,但是计算时间较长,在某些情况下求解精度不高,造成加工误差。随着加工零件日益复杂,机床结构不断更新,对后置处理的求解效率和求解精度提出了更高的要求,因此具有快速性和精确性的通用后置求解算法具有重要意义。
发明内容:
针对现有的通用后置处理方法的不足,本发明提出一种多轴机床通用后置处理方法,根据机床结构建立机床运动变换链,再采用数值方法和全微分方法计算刀位源文件中其他刀位点对应的机床各轴运动量,进而获得数控机代码。本发明为工艺人员提供一种适合于各种机床结构的通用后置方法,减少了人工推导计算公式的环节,有利于提高后置处理的速度和效率。
实现本发明的目的所采用的具体技术方案如下:
一种多轴数控机床通用后置处理方法,包含如下步骤:
(一)根据机床结构建立机床运动变换链,机床运动链包括刀具到机床运动链QMT和工件到机床运动链QWM,刀具刀位点方向矢量和位置经运动变换得到工件上加工点的方向矢量和位置信息;
(二)根据步骤(一)中所建立的机床运动变换链,建立从刀具到工件的变换矩阵QWT
(三)输入初始刀位点CL1的机床各平动轴运动坐标Δ1
(四)计算出初始刀位点CL1后续的各刀位点CLi+1的平动轴运动坐标Δi+1,具体为:
1.输入刀位点CLi+1的位置矩阵该矩阵表明该刀位点的位置坐标。第一次进入步骤(四)中的子步骤(1)时i的取值为1,其中i=1,2…,n;
2.设定全微分算法初始迭代条件,即取CLi+1,1为迭代计算初始刀位点,通过迭代计算,最终使其与目标刀位点CLi+1重合,自然,刀位点CLi+1,1对应的机床运动坐标与目标刀位点CLi+1相同,即Δi+1,1=Δi(下标中逗号前数值表示第i+1个刀位点,逗号后的数值表示进行第一次迭代)
3.令k=1,利用全微分算法对每个刀位点CLi+1进行迭代计算,建立第k+1次迭代计算的中间刀位点CLi+1,k+1与第k次迭代计算的中间刀位点CLi+1,k之间机床各平动轴和转动轴的运动坐标增量dΔk的关系,求解出dΔk。第一次进入(四)中3时k的取值为1;
4.根据步骤(四)中子步骤(3)得到的机床各平动轴和转动轴运动坐标增量dΔk,求出第k+1次迭代计算的中间刀位点CLi+1,k+1的机床运动坐标Δi+1,k+1
5.根据机床运动变换链,通过所述第k+1次迭代计算中间刀位点CLi+1,k+1对应的机床运动坐标Δi+1,k+1,得到刀位点CLi+1,k+1对应的位置矩阵
Figure BDA0000101307700000031
6.如果第k+1次迭代计算中间刀位点CLi+1,k+1对应的机床运动坐标Δi+1,k+1对应的位置矩阵
Figure BDA0000101307700000032
与目标刀位点CLi+1对应的位置矩阵小于设定的误差δ之内,或者k的值大于设定的迭代次数t,则此时目标刀位点CLi+1所对应的机床运动坐标Δi+1=Δi+1,k+1,否则,将k的值增加1,返回步骤(四)中子步骤(3),,进行再次计算;
(五)判断i的值是否小于n,如果小于n则返回步骤(四)中子
步骤1,将i增加1后继续计算,否则计算结束,输出各刀位点对应的机床平动轴运动坐标和旋转轴运动坐标。
本发明通过机床结构建立机床运动变换链,根据微分算法求解刀位源文件中各刀位点对应的机床各轴运动量,得到机床机代码。
本发明所建立的多轴数控机床通用后置处理方法,其有益成果在于:本发明消除了手工推导方法中机床各轴运动量计算公式的推导环节,简化了机床后置处理的流程,提高后置处理的速度;本发明消除了数值处理方法计算时间长,求解精度不稳定的缺点,提高了后置处理的精度;适用于任意结构多轴机床后置处理,具有求解速度快,求解精度高的特点。
附图说明:
图1为多轴机床通用后置处理流程图
具体实施方式:
下面结合附图和具体实施例对本发明作进一步说明。
多轴机床通用后置处理方法,包括如下步骤:
(一)对多轴数控机床各运动部件之间相对位置关系进行分析测量,建立机床运动变换链,机床运动链包括刀具到机床运动链QMT和工件到机床运动链TWM,刀具刀位点方向矢量和位置经运动变换得到工件上加工点的方向矢量和位置信息,关系如下面方程组所示:
Q MW U W 0 T = Q MT 0 0 1 0 T Q MW P W 1 T = Q MT 0 0 0 1 T
其中QMT是刀具相对于机床床身的变换矩阵,QMW是工件相对于机床床身的变换矩阵。[0010]T是刀具坐标系中刀位点位置矩阵,[0001]T是刀具坐标系中刀位点方向矢量,UW是工件坐标系中刀位点位置矩阵,该矢量通过UX,UY,UZ三个数值分别描述该位置矩阵在X,Y,Z方向上的分量,PW是工件坐标系中刀位点方向矢量,该矢量通过PX,PY,PZ三个数值表示该方向矢量在X,Y,Z方向上的分量。
具体实现步骤如下:
1.建立刀具到机床运动链的变换矩阵QMT,该变换矩阵计算公式如下:
QMT=QM,n×…×Qi,j×…Q1,T
其中QM,n是刀具到机床运动链中第n个运动部件到机床床身的运动变换矩阵,Qi,j是刀具到机床运动链中第j个运动部件到第i个运动部件的运动变换矩阵,Q1,T是刀具到机床运动链中刀具到第1个运动链的变换矩阵,上述QM,n、Qi,j和Q1,T三类变换矩阵反映机床不同运动部件相对位置关系,该类矩阵可以通过机床运动模型导出。
2.建立从工件到机床运动链的变换矩阵QMW,该变换矩阵计算公式如下:
QMW=QM,n×…×Qi,j×…Q1,W
其中QM,n是工件到机床运动链中第n个运动部件到机床床身的运动变换矩阵,Qi,j是工件到机床运动链中第j个运动部件到第i个运动部件的运动变换矩阵,Q1,W是工件到机床运动链中刀具到第1个运动链的变换矩阵。
上述QM,n、Qi,j和Q1,W三类变换矩阵反映机床不同运动部件相对位置关系,该类矩阵可以通过机床运动模型导出。
3.将1和2步骤中得到的从刀具到机床运动链的变换矩阵QMT和从工件到机床运动链的变换矩阵QMW带入如下方程组:
Q MW U W 0 T = Q MT 0 0 1 0 T Q MW P W 1 T = Q MT 0 0 0 1 T
即可建立机床运动变换链。
(二)根据(一)中所建立的机床运动变换链,可以建立从刀具到工件的变换矩阵QWT
QWT=QMW -1×QMT
QMW -1为(一)中所建立的工件相对于机床床身的变换矩阵的逆矩阵,而QMT为(一)中所建立的刀具相对于机床床身的变换矩阵。
(三)输入初始刀位点CL1的机床各平动轴运动坐标Δ1
(四)计算出刀位点CLi+1的平动轴运动坐标Δi+1
1.输入刀位点CLi+1的位置矩阵
Figure BDA0000101307700000061
该矩阵表明该刀位点的位置坐标。第一次进入(四)时i的取值为1,其中i=1,2…,n;
2.设定全微分算法初始迭代条件,取CLi+1,1为迭代计算初始刀位点,通过迭代计算,最终使其与目标刀位点CLi+1重合,自然,刀位点CLi+1,1对应的机床运动坐标与目标刀位点CLi+1相同,即Δi+1,i=Δi
3.根据全微分方法,可以建立第k+1次迭代计算中间刀位点CLi+1,k+1与第k次迭代计算中间刀位点CLi+1,k之间机床各平动轴和转动轴运动坐标增量dΔk的关系,求解出dΔk。第一次进入(四)中3时k的取值为1;
4.根据(四)中3得到的机床各平动轴和转动轴运动坐标增量dΔk,即可求出第k+1次迭代计算中间刀位点CLi+1,k+1的机床运动坐标Δi+1,k+1
5.根据机床运动变换链,通过第k+1次迭代计算中间刀位点CLi+1,k+1对应的机床运动坐标Δi+1,k+1,得到刀位点CLi+1,k+1对应的位置矩阵
Figure BDA0000101307700000071
6.如果第k+1次迭代计算中间刀位点CLi+1,k+1对应的机床运动坐标对Δi+1,k+1应的位置矩阵
Figure BDA0000101307700000072
与目标刀位点CLi+1对应的位置矩阵
Figure BDA0000101307700000073
小于设定的误差δ(如1e-3)之内,或者k的值大于设定的迭代次数t(如10000),则此时目标刀位点CLi+1所对应的机床运动坐标Δi+1=Δi+1,k+1,否则,返回(六)中,将k的值增加1,再次计算。
上述1-6步骤如图2所示,其中虚线框内的内容为(四)步骤中的输入量,而实线框内的内容则对应1-6步骤。
(五)判断i的值是否小于n,如果小于n则返回(四)中,将i增加1后继续计算,否则计算结束,输出各刀位点对应的机床平动轴运动坐标和旋转轴运动坐标。
整体流程图如图1表示。

Claims (1)

1.一种多轴数控机床通用后置处理方法,包括如下步骤:
(1)根据机床结构建立机床运动变换链,其中,所述机床运动链包括刀具到机床运动链QMT,工件到机床运动链QWM,以及刀具刀位点方向矢量和位置经运动变换得到的工件上加工点的方向矢量和位置;
(2)根据步骤(1)中所建立的机床运动变换链,建立从刀具到工件的变换矩阵QWT
(3)输入初始刀位点CL1的机床各平动轴运动坐标Δ1
(4)计算出所述初始刀位点CL1后续的任一刀位点CLi+1的平动轴运动坐标Δi+1,其中i=1,2…,n,n为正整数,具体为:
(4.1)输入所述刀位点CLi+1的位置矩阵
Figure FDA0000101307690000011
该矩阵表明该刀位点的位置坐标,其中第一次进入步骤(4.1)时i的取值为1;
(4.2)利用全微分算法对刀位点CLi+1进行迭代计算,建立第k+1次迭代计算的中间刀位点CLi+1,k+1与前一次迭代计算的中间刀位点CLi+1,k之间机床各平动轴和转动轴的运动坐标增量dΔk的关系,求解出dΔk,其中第一次进入步骤(4.2)时k的取值为1;
(4.3)根据所述的运动坐标增量dΔk,求出第k+1次迭代计算的中间刀位点CLi+1,k+1的机床运动坐标Δi+1,k+1
(4.4)根据机床运动变换链,通过所述第k+1次迭代计算中间刀位点CLi+1,k+1对应的机床运动坐标Δi+1,k+1,得到刀位点CLi+1,k+1对应的位置矩阵
Figure FDA0000101307690000012
(4.5)如果第k+1次迭代计算中间刀位点CLi+1,k+1对应的机床运动坐标Δi+1,k+1对应的位置矩阵与目标刀位点CLi+1对应的位置矩阵小于设定的误差δ之内,或者k的值大于设定的迭代次数t,则此时目标刀位点CLi+1所对应的机床运动坐标Δi+1=Δi+1,k+1,否则,将k的值增加1,跳转到步骤(4.2)重新计算。
(5)判断i的值是否小于n,如果小于n则返回步骤(4.1),将i增加1后循环计算,否则,计算结束,输出各刀位点对应的机床平动轴运动坐标和旋转轴运动坐标。
CN 201110324388 2011-10-24 2011-10-24 一种多轴数控机床通用后置处理方法 Active CN102402198B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110324388 CN102402198B (zh) 2011-10-24 2011-10-24 一种多轴数控机床通用后置处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110324388 CN102402198B (zh) 2011-10-24 2011-10-24 一种多轴数控机床通用后置处理方法

Publications (2)

Publication Number Publication Date
CN102402198A true CN102402198A (zh) 2012-04-04
CN102402198B CN102402198B (zh) 2013-03-27

Family

ID=45884510

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110324388 Active CN102402198B (zh) 2011-10-24 2011-10-24 一种多轴数控机床通用后置处理方法

Country Status (1)

Country Link
CN (1) CN102402198B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106292542A (zh) * 2016-09-26 2017-01-04 北京航空航天大学 一种带c轴车加工的通用后置处理方法和装置
CN106802628A (zh) * 2017-02-22 2017-06-06 西安应用光学研究所 五轴孔加工通用后置处理方法
CN104345687B (zh) * 2013-08-09 2017-10-20 中国科学院沈阳计算技术研究所有限公司 一种刀具扫描体建模方法
CN107703744A (zh) * 2017-09-19 2018-02-16 浙江大学 考虑非线性误差和广义模态的数控机床运动链设计方法
CN108549323A (zh) * 2018-04-18 2018-09-18 西南交通大学 一种五轴数控工具磨床通用后置处理方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001105278A (ja) * 1999-10-05 2001-04-17 Toyoda Mach Works Ltd 工作機械における先使用・後使用工具の刃先位置整合方法
CN1624613A (zh) * 2004-12-11 2005-06-08 东方汽轮机厂 一种四轴联动数控加工后置处理方法及其控制系统
WO2006031821A2 (en) * 2004-09-14 2006-03-23 Ugs Corp. System, method, and computer program product for machine tool programming
CN101859126A (zh) * 2010-04-30 2010-10-13 西北工业大学 五坐标数控加工刀具路径转换方法
CN102023613A (zh) * 2010-12-29 2011-04-20 广西工学院 一种五轴联动数控加工后置处理器及其处理方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001105278A (ja) * 1999-10-05 2001-04-17 Toyoda Mach Works Ltd 工作機械における先使用・後使用工具の刃先位置整合方法
WO2006031821A2 (en) * 2004-09-14 2006-03-23 Ugs Corp. System, method, and computer program product for machine tool programming
CN1624613A (zh) * 2004-12-11 2005-06-08 东方汽轮机厂 一种四轴联动数控加工后置处理方法及其控制系统
CN101859126A (zh) * 2010-04-30 2010-10-13 西北工业大学 五坐标数控加工刀具路径转换方法
CN102023613A (zh) * 2010-12-29 2011-04-20 广西工学院 一种五轴联动数控加工后置处理器及其处理方法

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104345687B (zh) * 2013-08-09 2017-10-20 中国科学院沈阳计算技术研究所有限公司 一种刀具扫描体建模方法
CN106292542A (zh) * 2016-09-26 2017-01-04 北京航空航天大学 一种带c轴车加工的通用后置处理方法和装置
CN106292542B (zh) * 2016-09-26 2019-02-22 北京航空航天大学 一种带c轴车加工的通用后置处理方法和装置
CN106802628A (zh) * 2017-02-22 2017-06-06 西安应用光学研究所 五轴孔加工通用后置处理方法
CN106802628B (zh) * 2017-02-22 2019-03-29 西安应用光学研究所 五轴孔加工通用后置处理方法
CN107703744A (zh) * 2017-09-19 2018-02-16 浙江大学 考虑非线性误差和广义模态的数控机床运动链设计方法
CN107703744B (zh) * 2017-09-19 2020-07-14 浙江大学 考虑非线性误差和广义模态的数控机床运动链设计方法
CN108549323A (zh) * 2018-04-18 2018-09-18 西南交通大学 一种五轴数控工具磨床通用后置处理方法

Also Published As

Publication number Publication date
CN102402198B (zh) 2013-03-27

Similar Documents

Publication Publication Date Title
CN102591257B (zh) 面向参数曲线刀具轨迹的数控系统轮廓误差控制方法
CN101980091B (zh) 双转台五轴联动数控加工旋转刀具中心点补偿方法
CN102402198B (zh) 一种多轴数控机床通用后置处理方法
CN112558547B (zh) 一种五轴数控机床平动轴几何误差补偿数据快速优化方法
CN102023613A (zh) 一种五轴联动数控加工后置处理器及其处理方法
CN101653921A (zh) 凸轮轴数控磨削轮廓误差补偿方法
Lei et al. NURBS-based fast geometric error compensation for CNC machine tools
CN101639681B (zh) 一种电子装备运动机构性能参数优化方法
CN101870073B (zh) 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法
CN104460515A (zh) 一种基于后置处理五轴刀具长度补偿方法
CN102073301B (zh) 具有刀具长度补偿功能的五轴样条插补器
CN110597183B (zh) 一种磨齿关键误差高效补偿方法
CN106125674A (zh) 一种高精度实时轮廓误差估计方法
CN103760827A (zh) 跃度约束的数控加工进给率的离线规划方法
Liu et al. A method of energy-optimal trajectory planning for palletizing robot
Wang et al. Design and implementation of five-axis transformation function in CNC system
Li et al. Off-line optimization on NC machining based on virtual machining
Li et al. Simulation and comparison research of Lagrange and Kane dynamics modeling for the 4-DOF modular industrial robot
CN104155909B (zh) 基于pc+fpga的带速度观测器的有限时间稳定性控制系统
CN102426436A (zh) 考虑机床结构误差的多轴数控加工后置处理方法
CN202985116U (zh) 五轴联动强力数控机床
CN108415371A (zh) 一种通用五轴数控机床几何误差补偿方法
Wang et al. Ball End Mill—Tool Radius Compensation of Complex NURBS Surfaces for 3-Axis CNC Milling Machines
CN115822879B (zh) 采用三维拟合调试叶片组装工装的方法、设备及存储介质
Lang et al. Simulation of Kinematic Accuracy of a 6-DOF Serial Robotic Mechanism Based on Screw Theory

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