CN108549319B - 一种双转台五轴数控机床通用后处理方法 - Google Patents

一种双转台五轴数控机床通用后处理方法 Download PDF

Info

Publication number
CN108549319B
CN108549319B CN201810282519.1A CN201810282519A CN108549319B CN 108549319 B CN108549319 B CN 108549319B CN 201810282519 A CN201810282519 A CN 201810282519A CN 108549319 B CN108549319 B CN 108549319B
Authority
CN
China
Prior art keywords
axis
machine tool
angle
tool
motion
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.)
Active
Application number
CN201810282519.1A
Other languages
English (en)
Other versions
CN108549319A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201810282519.1A priority Critical patent/CN108549319B/zh
Publication of CN108549319A publication Critical patent/CN108549319A/zh
Application granted granted Critical
Publication of CN108549319B publication Critical patent/CN108549319B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/404Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开了一种双转台五轴数控机床通用后处理方法,属于数控机床加工领域。包括:引入机床坐标系下旋转轴和工件的位置向量,建立双转台五轴数控机床运动链;结合各个运动轴相对于全局坐标系的位置向量,根据指数积理论建立刀具相对于工件的齐次运动矩阵;建立机床正向运动学方程;根据刀具姿态表达式,求解旋转轴角度的两组解;考虑C轴真实运动范围,与上一个刀具姿态的加工代码比较获得C轴真实角度,更新旋转轴角度的两组解;建立最短旋转轴角度路径原则,选择合适的旋转轴角度;根据刀具位置表达式,计算三个平动轴运动量;读取工件刀位文件,计算得到工件加工代码。该发明生成的加工代码反映了运动轴相对于其零位置的运动量,可直接用于加工。

Description

一种双转台五轴数控机床通用后处理方法
技术领域
本发明涉及数控机床加工制造技术领域,尤其涉及五轴数控机床后处理技术。
背景技术
随着科技发展,五轴数控机床已经广泛使用于制造行业。五轴数控机床后处理是机床可加工的重要功能之一。机床CAM系统是产生工件加工路径,数控系统是根据加工代码驱动机床运动轴运动进行加工,而五轴数控机床后处理是CAM系统和数控系统的,其作用是根据CAM系统生成的加工路径产生相应的加工代码。
目前五轴数控机床后处理大多是通过齐次变换矩阵得到机床正向运动学方程,然后根据逆向运动学计算运动轴运动量。由于在采用齐次变换矩阵建立机床运动学过程中忽略了旋转轴在机床中位置,使得生产的平动轴运动量没有反映其真实位置,无法直接进行加工。另外,大多后处理探索了旋转轴角度计算机制,并没有解决旋转轴角度双解求解和选择的问题。由于现有五轴后处理生成的加工代码没有反映各个运动轴相对于其零位置的运动量,所生成的加工代码无法用于机床误差补偿。这种后处理需要借助机床RTCP功能来实现工件加工,这使得五轴机床不得不配置高端的五轴数控系统。另外,RTCP的代码生成功能并不对用户开放,用户无法获得机床加工工件时运动轴的运动量。
依据旋转轴在五轴数控机床中位置可将其分为三大类,双转台五轴数控机床是其中之一,在工业生产制造中占用重要地位。由于平动轴位置的不同,双转台五轴数控机床也有很多种。现有后处理研究从机理上探索了五轴机床的代码生成机制,并没有针对某类机床研究通用的后处理。随着用户对系统开放性需求的提高,研究可直接用于工件加工的通用的后处理具有重要意义。
发明内容
本发明的目的是提供一种双转台五轴数控机床通用后处理方法,它能有效地获得所有运动轴相对于其零位置的运动量,能够直接用于工件加工。
本发明的目的是通过以下技术方案来实现的:
一种双转台五轴数控机床通用后处理方法,包括如下步骤:
步骤1、引入机床坐标系下旋转轴和工件的位置向量,建立双转台五轴数控机床运动链,具体步骤包括:
步骤1.1、机床坐标系原点建立在双转台五轴数控机床三个平动轴零位置时刀尖点处;
步骤1.2、获得C轴旋转轴线在机床坐标系中的位置向量VC=[Cx,Cy,0]T,其中Cx表示C轴旋转轴线x方向坐标,Cy表示C轴旋转轴线y方向坐标;
步骤1.3、获得A轴旋转轴线在机床坐标系中的位置向量VA=[0,Ay,Az]T,其中Ay表示A轴旋转轴线y方向坐标,Az表示A轴旋转轴线z方向坐标;
步骤1.4、获得工件在机床坐标系中的位置向量Vw=[Wx,Wy,Wz]T,其中Wx表示工件x方向坐标,Wy表示工件y方向坐标,Wz表示工件z方向坐标;
步骤1.5、根据双转台五轴数控机床拓扑结构,得到运动轴拓扑关系,建立双转台五轴数控机床运动链。
步骤2、结合各个运动轴相对于全局坐标系的位置向量,根据指数积理论,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵;
步骤3、建立双转台五轴数控机床正向运动学方程,包括刀具位置表达式和刀具姿态表达式;
步骤4、根据刀具姿态表达式和刀具姿态,结合旋转轴运动范围,求解旋转轴角度的两组解;
步骤5、考虑C轴真实运动范围,将步骤4得到的旋转轴角度的两组解与上一个刀具姿态对应的加工代码比较,获得C轴真实角度,并更新旋转轴角度的两组解;
步骤6、建立最短旋转轴角度路径原则,将步骤5更新的旋转轴角度的两组解与上一个刀具姿态对应的加工代码比较,选择合适的旋转轴角度;
步骤7、根据刀具位置表达式,结合步骤6中选择的旋转轴角度,计算刀具位置对应的三个平动轴运动量:
Figure GDA0002596668250000021
其中[Px,Py,Pz]T表示刀具位置,[Ox,Oy,Oz]T表示刀具姿态,x,y,z分别表示X轴、Y轴和Z轴的相对于其零位置的运动量,α和γ分别表示A轴和C轴的旋转角度;
步骤8、获得双转台五轴数控机床的机床坐标系下旋转轴和工件的位置向量数值,读取工件刀位文件数据中刀具位置和刀具姿态,根据步骤4到步骤7计算每个刀具位置和刀具姿态的旋转轴角度和平动轴运动量,得到工件加工代码。
作为优选,所述步骤2中结合各个运动轴相对于全局坐标系的位置向量,根据指数积理论,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵,具体步骤包括:
步骤2.1、将双转台五轴数控全局坐标系固定在工件原点处;
步骤2.2、结合机床坐标系下旋转轴和工件的位置向量,获得各个运动轴相对于全局坐标系的位置向量,其中C轴旋转轴线相对于全局坐标系的位置向量为LC=-VW+VC,A轴旋转轴线相对于全局坐标系的位置向量为LA=-VW+VA,X轴相对于全局坐标系的位置向量为LX=[x-Wx,0,0]T,Y轴相对于全局坐标系的位置向量为LY=[0,y-Wy,0]T,Z轴相对于全局坐标系的位置向量为LZ=[0,0,z-Wz]T
步骤2.3、根据指数积理论,结合C轴旋转轴线相对于全局坐标系的位置向量,建立C轴运动旋量和运动指数矩阵;
步骤2.4、根据指数积理论,结合A轴旋转轴线相对于全局坐标系的位置向量,建立A轴运动旋量和运动指数矩阵;
步骤2.5、根据指数积理论,结合三个平动轴相对于全局坐标系的位置向量,建立X轴、Y轴和Z轴的运动旋量和运动指数矩阵;
步骤2.6、根据双转台五轴数控机床运动链,结合各个运动轴的运动指数矩阵,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵,计算表达式为:
Tt w=TC·TA·Tx·Ty·Tz
其中Tt w表示刀具相对于工件的齐次运动矩阵,TC、TA、Tx、Ty和Tz分别表示C轴、A轴、X轴、Y轴和Z轴的运动指数矩阵。
作为优选,所述步骤3中建立双转台五轴数控机床正向运动学方程,得到的刀具位置表达式为:
Figure GDA0002596668250000031
作为优选,所述步骤4中根据刀具姿态表达式和刀具姿态,结合旋转轴运动范围,求解旋转轴角度的两组解的具体方法为:
步骤4.1、根据刀具姿态表达式,暂时忽略旋转轴运动范围,得到旋转轴角度表达式为:
Figure GDA0002596668250000032
其中a和c分别表示满足刀具姿态表达式的A轴和C轴角度的一个值。
步骤4.2、根据刀具姿态表达式,考虑A轴运动范围为[-π/2,π/2],建立A轴角度表达式为:
Figure GDA0002596668250000033
步骤4.3、考虑C轴运动范围为[0,2π],结合反正弦函数的周期性,建立当Oy≠0时C轴角度两组解的表达式为:
Figure GDA0002596668250000034
其中γ1表示C轴角度两组解中第一组解,γ2表示C轴角度两组解中第二组解。
步骤4.4、结合旋转轴角度表达式、A轴角度表达式和当Oy≠0时C轴角度两组解的表达式,求解得到旋转轴角度的两组解为:
Figure GDA0002596668250000041
其中α1表示A轴角度两组解中第一组解,α2表示A轴角度两组解中第二组解.
作为优选,所述步骤5中C轴真实运动范围为-∞~∞;
作为优选,所述步骤5中将步骤4得到的旋转轴角度的两组解与上一个刀具姿态对应的加工代码比较,获得C轴真实角度为:
Figure GDA0002596668250000042
其中γr1表示C轴角度两组解中第一组解的真实角度,γr2表示C轴角度两组解中第二组解的真实角度,γpre表示上一个刀具姿态对应的加工代码中C轴角度,k1和k2为整数。
作为优选,所述步骤6中最短旋转轴角度路径原则为将得到的旋转轴角度两组解与上一个刀具姿态对应的加工代码中旋转轴角度比较,选择旋转轴角度路径变化最小的组合。旋转轴角度路径表达式为:min|α-αpre|+|γ-γpre|。其中αpre表示上一个刀具姿态对应的加工代码中A轴角度。
作为优选,所述步骤8中获得双转台五轴数控机床的机床坐标系下旋转轴的位置向量数值的方法为查询双转台五轴数控机床结构参数;
作为优选,所述步骤8中获得双转台五轴数控机床的机床坐标系下工件的位置向量的方法为安装好工件毛坯后通过机床对刀功能获得;
作为优选,所述步骤5和所述步骤6中第一个刀具位置和刀具姿态的所述的上一个刀具姿态对应的加工代码为零;
作为优选,所述步骤6中第一个刀具位置和刀具姿态的选择的合适的旋转轴角度为旋转轴角度的两组解中第一组解;
本发明是一种双转台五轴数控机床通用后处理方法,具体的有益效果是:本发明引入了旋转轴和工件在机床坐标系下的位置向量,使得得到的加工代码反映了所有运动轴相对于其零位置的真实运动量,处理了旋转轴角度两组解的求解和旋转问题,建立了真实可用的五轴机床后处理,可用于双转台五轴数控机床的工件加工。
附图说明
图1为本发明流程图;
图2为本发明双转台五轴数控机床运动链示意图;
图3为C轴角度与上一个刀具姿态对应的加工代码中C轴角度关系图;
图4为某工件刀具路径示意图;
图5为本发明双转台五轴数控机床通用后处理方法生成的工件加工代码;
图6为根据加工代码进行加工的工件效果图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步说明。
附图1所示为本发明双转台五轴数控机床通用后处理方法流程图。以某双转台五轴数控机床为例阐述双转台五轴数控机床通用后处理方法。
步骤1、引入机床坐标系下旋转轴和工件的位置向量,建立双转台五轴数控机床运动链,具体步骤包括:
步骤1.1、机床坐标系原点建立在双转台五轴数控机床三个平动轴零位置时刀尖点处;
步骤1.2、获得C轴旋转轴线在机床坐标系中的位置向量VC=[Cx,Cy,0]T,其中Cx表示C轴旋转轴线x方向坐标,Cy表示C轴旋转轴线y方向坐标;
步骤1.3、获得A轴旋转轴线在机床坐标系中的位置向量VA=[0,Ay,Az]T,其中Ay表示A轴旋转轴线y方向坐标,Az表示A轴旋转轴线z方向坐标;
步骤1.4、获得工件在机床坐标系中的位置向量Vw=[Wx,Wy,Wz]T,其中Wx表示工件x方向坐标,Wy表示工件y方向坐标,Wz表示工件z方向坐标;
步骤1.5、根据双转台五轴数控机床拓扑结构,得到运动轴拓扑关系,建立双转台五轴数控机床运动链。附图2所示的双转台五轴数控机床的运动轴拓扑关系为工作台→C轴→A轴→Y轴→床身→X轴→Z轴→刀具,从而得到该机床运动链,附图2所示为包含旋转轴和工件位置向量的双转台五轴数控机床运动链。
步骤2、结合各个运动轴相对于全局坐标系的位置向量,根据指数积理论,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵,具体步骤包括:
步骤2.1、将双转台五轴数控全局坐标系固定在工件原点处;
步骤2.2、结合机床坐标系下旋转轴和工件的位置向量,获得各个运动轴相对于全局坐标系的位置向量,其中C轴旋转轴线相对于全局坐标系的位置向量为LC=-VW+VC,A轴旋转轴线相对于全局坐标系的位置向量为LA=-VW+VA,X轴相对于全局坐标系的位置向量为LX=[x-Wx,0,0]T,Y轴相对于全局坐标系的位置向量为LY=[0,y-Wy,0]T,Z轴相对于全局坐标系的位置向量为LZ=[0,0,z-Wz]T
步骤2.3、根据指数积理论,结合C轴旋转轴线相对于全局坐标系的位置向量,建立C轴运动旋量和运动指数矩阵为:
Figure GDA0002596668250000051
Figure GDA0002596668250000052
其中γ表示C轴旋转角度,
Figure GDA0002596668250000053
表示C轴运动旋量,TC表示C轴运动指数矩阵。
步骤2.4、根据指数积理论,结合A轴旋转轴线相对于全局坐标系的位置向量,建立A轴运动旋量和运动指数矩阵为:
Figure GDA0002596668250000061
Figure GDA0002596668250000062
其中α表示C轴旋转角度,
Figure GDA0002596668250000063
表示A轴运动旋量,TA表示A轴运动指数矩阵。
步骤2.5、根据指数积理论,结合三个平动轴相对于全局坐标系的位置向量,建立X轴、Y轴和Z轴的运动旋量和运动指数矩阵为:
Figure GDA0002596668250000064
Figure GDA0002596668250000065
其中
Figure GDA0002596668250000066
Figure GDA0002596668250000067
分别表示X轴、Y轴和Z轴的运动旋量,Tx、Ty和Tz分别表示X轴、Y轴和Z轴的运动指数矩阵。
步骤2.6、根据双转台五轴数控机床运动链,结合各个运动轴的运动指数矩阵,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵。平动轴运动轴指数矩阵为平动齐次变换矩阵,平动齐次变换矩阵之间相乘与矩阵次序无关,所以双转台五轴数控机床刀具相对于工件的齐次运动矩阵表达式可统一表示为:
Tt w=TC·TA·Tx·Ty·Tz
其中Tt w表示刀具相对于工件的齐次运动矩阵;
步骤3、建立双转台五轴数控机床正向运动学方程,包括刀具位置表达式和刀具姿态表达式。刀具位置和刀具姿态可根据机床刀具相对于工件的齐次运动矩阵得到:
[Px,Py,Pz,1]T=Tt w·[0,0,0,1]T
[Ox,Oy,Oz,0]T=Tt w·[0,0,1,0]T
展开后获得刀具位置表达式和刀具姿态表达式为:
Figure GDA0002596668250000068
Figure GDA0002596668250000069
步骤4、根据刀具姿态表达式和刀具姿态,结合旋转轴运动范围,求解旋转轴角度的两组解,具体方法为:
步骤4.1、根据刀具姿态表达式,暂时忽略旋转轴运动范围,得到旋转轴角度表达式为:
Figure GDA0002596668250000071
其中a和c分别表示满足刀具姿态表达式的A轴和C轴角度的一个值。
步骤4.2、根据刀具姿态表达式,考虑A轴运动范围为[-π/2,π/2],建立A轴角度表达式为:
Figure GDA0002596668250000072
步骤4.3、考虑C轴运动范围为[0,2π],结合反正弦函数的周期性,建立当Oy≠0时C轴角度两组解的表达式为:
Figure GDA0002596668250000073
其中γ1表示C轴角度两组解中第一组解,γ2表示C轴角度两组解中第二组解。
步骤4.4、结合旋转轴角度表达式、A轴角度表达式和当Oy≠0时C轴角度两组解的表达式,求解得到旋转轴角度的两组解为:
Figure GDA0002596668250000074
其中α1表示A轴角度两组解中第一组解,α2表示A轴角度两组解中第二组解.
步骤5、考虑C轴真实运动范围,将步骤4得到的旋转轴角度的两组解与上一个刀具姿态对应的加工代码比较,获得C轴真实角度,并更新旋转轴角度的两组解。机床C轴可以沿着一个方向一直旋转,理论上C轴的真实运动范围为(-∞,∞),将步骤4得到旋转轴角度的两组解与与上一个刀具姿态对应的加工代码比较,附图3所示为C轴角度与上一个刀具姿态对应的加工代码中C轴角度关系图,从而得到C轴真实角度为:
Figure GDA0002596668250000081
其中γr1表示C轴角度两组解中第一组解的真实角度,γr2表示C轴角度两组解中第二组解的真实角度,γpre表示上一个刀具姿态对应的加工代码中C轴角度。从而对步骤4中的旋转轴角度的两组解进行更新,k1和k2为整数。
步骤6、建立最短旋转轴角度路径原则,表示为将得到的旋转轴角度与上一个刀具姿态对应的加工代码中旋转轴角度比较,选择旋转轴角度路径变化最小的组合,旋转轴角度路径表达式为:min|α-αpre|+|γ-γpre|。其中αpre表示上一个刀具姿态对应的加工代码中A轴角度。然后根据此原则,从步骤5更新的旋转轴角度的两组解中选择合适的旋转轴角度,
步骤7、根据刀具位置表达式,结合步骤6中选择的旋转轴角度,计算刀具位置对应的三个平动轴运动量:
Figure GDA0002596668250000082
其中[Px,Py,Pz]T表示刀具位置,[Ox,Oy,Oz]T表示刀具姿态,x,y,z分别表示X轴、Y轴和Z轴的相对于其零位置的运动量,α和γ分别表示A轴和C轴的旋转角度;
步骤8、根据双转台五轴数控机床结构参数获得机床坐标系下旋转轴位置向量数值,安装好工件毛坯后通过机床对刀功能得到工件位置向量对数。附图2所示的双转台五轴数控机床的C轴位置向量为VC=[244.5247,105.9248,0]Tmm,A轴位置向量为VA=[0,105.9232,-237.154]T mm,工件安装后得到的工件位置向量为Vw=[244.7576,55.3002,-153.3136]T mm。然后读取工件刀位文件数据中刀具位置和刀具姿态,根据步骤4到步骤7计算每个刀具位置和刀具姿态的旋转轴角度和平动轴运动量,得到工件加工代码。附图4所示某工件的刀具路径示意图,读取该工件的刀位文件,获得相应的刀具位置和刀具姿态数据,根据步骤4到步骤7计算每个刀具位置和刀具姿态的旋转轴角度和平动轴运动量,得到工件加工代码。对于第一个刀具位置和刀具姿态,在步骤5和步骤6中将其对应的上一个刀具姿态对应的加工代码设置为零,也可以直接在骤6中第选择旋转轴角度的两组解中第一组解为合适的旋转轴角度。附图5所示为生成的工件加工代码。为了验证本发明双转台五轴数控机床通用后处理方法的有效性,根据生成的加工代码进行Vericut模拟加工来进行验证。在vericut中根据A轴和C轴位置向量设置机床状态,并关闭数控系统RTCP功能,读入本发明双转台五轴数控机床通用后处理方法生成的加工代码,对工件进行加工。附图6所示为工的工件效果图,加工工件的刀具路径分布与附图4中设计的刀具路径相同,说明本发明双转台五轴数控机床通用后处理方法可以生产有效的加工代码,并可以直接用于工件加工。

Claims (9)

1.一种双转台五轴数控机床通用后处理方法,包括如下步骤:
步骤1、引入机床坐标系下旋转轴和工件的位置向量,建立双转台五轴数控机床运动链,具体步骤包括:
步骤1.1、机床坐标系原点建立在双转台五轴数控机床三个平动轴零位置时刀尖点处;
步骤1.2、获得C轴旋转轴线在机床坐标系中的位置向量VC=[Cx,Cy,0]T;其中,Cx表示C轴旋转轴线x方向坐标,Cy表示C轴旋转轴线y方向坐标,上标T为矩阵转置;
步骤1.3、获得A轴旋转轴线在机床坐标系中的位置向量VA=[0,Ay,Az]T,其中Ay表示A轴旋转轴线y方向坐标,Az表示A轴旋转轴线z方向坐标;
步骤1.4、获得工件在机床坐标系中的位置向量Vw=[Wx,Wy,Wz]T,其中Wx表示工件x方向坐标,Wy表示工件y方向坐标,Wz表示工件z方向坐标;
步骤1.5、根据双转台五轴数控机床拓扑结构,得到运动轴拓扑关系,建立双转台五轴数控机床运动链;
步骤2、结合各个运动轴相对于全局坐标系的位置向量,根据指数积理论,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵;
步骤3、建立双转台五轴数控机床正向运动学方程,包括刀具位置表达式和刀具姿态表达式;
步骤4、根据刀具姿态表达式和刀具姿态,结合旋转轴运动范围,求解旋转轴角度的两组解,具体方法为:
步骤4.1、根据刀具姿态表达式,暂时忽略旋转轴运动范围,得到旋转轴角度表达式为:
Figure FDA0002596668240000011
其中[Ox,Oy,Oz]T表示刀具姿态,a和c分别表示满足刀具姿态表达式的A轴和C轴角度的一个值;
步骤4.2、根据刀具姿态表达式,考虑A轴运动范围为[-π/2,π/2],建立A轴角度表达式为:
Figure FDA0002596668240000012
步骤4.3、考虑C轴运动范围为[0,2π],结合反正弦函数的周期性,建立当Oy≠0时C轴角度两组解的表达式为:
Figure FDA0002596668240000013
其中γ1表示C轴角度两组解中第一组解,γ2表示C轴角度两组解中第二组解;
步骤4.4、结合旋转轴角度表达式、A轴角度表达式和当Oy≠0时C轴角度两组解的表达式,求解得到旋转轴角度的两组解为:
Figure FDA0002596668240000021
其中α1表示A轴角度两组解中第一组解,α2表示A轴角度两组解中第二组解;
步骤5、考虑C轴真实运动范围,将步骤4得到的旋转轴角度的两组解与上一个刀具姿态对应的加工代码比较,获得C轴真实角度,并更新旋转轴角度的两组解;
步骤6、建立最短旋转轴角度路径原则,将步骤5更新的旋转轴角度的两组解与上一个刀具姿态对应的加工代码比较并选择所需的旋转轴角度;
步骤7、根据刀具位置表达式,结合步骤6中选择的旋转轴角度,计算刀具位置对应的三个平动轴运动量:
Figure FDA0002596668240000022
其中[Px,Py,Pz]T表示刀具位置,x,y,z分别表示X轴、Y轴和Z轴的相对于其零位置的运动量,α和γ分别表示A轴和C轴的旋转角度;
步骤8、获得双转台五轴数控机床的机床坐标系下旋转轴和工件的位置向量数值,读取工件刀位文件数据中刀具位置和刀具姿态,根据步骤4到步骤7计算每个刀具位置和刀具姿态的旋转轴角度和平动轴运动量,得到工件加工代码。
2.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤2中结合各个运动轴相对于全局坐标系的位置向量,根据指数积理论,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵,具体步骤包括:
步骤2.1、将双转台五轴数控全局坐标系固定在工件原点处;
步骤2.2、结合机床坐标系下旋转轴和工件的位置向量,获得各个运动轴相对于全局坐标系的位置向量,其中C轴旋转轴线相对于全局坐标系的位置向量为LC=-VW+VC,A轴旋转轴线相对于全局坐标系的位置向量为LA=-VW+VA,X轴相对于全局坐标系的位置向量为LX=[x-Wx,0,0]T,Y轴相对于全局坐标系的位置向量为LY=[0,y-Wy,0]T,Z轴相对于全局坐标系的位置向量为LZ=[0,0,z-Wz]T,其中x,y,z分别表示X轴、Y轴和Z轴的相对于其零位置的运动量;
步骤2.3、根据指数积理论,结合C轴旋转轴线相对于全局坐标系的位置向量,建立C轴运动旋量和运动指数矩阵;
步骤2.4、根据指数积理论,结合A轴旋转轴线相对于全局坐标系的位置向量,建立A轴运动旋量和运动指数矩阵;
步骤2.5、根据指数积理论,结合三个平动轴相对于全局坐标系的位置向量,建立X轴、Y轴和Z轴的运动旋量和运动指数矩阵;
步骤2.6、根据双转台五轴数控机床运动链,结合各个运动轴的运动指数矩阵,建立双转台五轴数控机床刀具相对于工件的齐次运动矩阵,计算表达式为:
Tt w=TC·TA·Tx·Ty·Tz
其中Tt w表示刀具相对于工件的齐次运动矩阵,TC、TA、Tx、Ty和Tz分别表示C轴、A轴、X轴、Y轴和Z轴的运动指数矩阵。
3.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤3中建立双转台五轴数控机床正向运动学方程,得到的刀具位置和刀具姿态表达式为:
Figure FDA0002596668240000031
4.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤5中C轴真实运动范围为-∞~∞。
5.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤5中将步骤4得到的旋转轴角度的两组解与上一个刀具姿态对应的加工代码比较,获得C轴真实角度为:
Figure FDA0002596668240000032
其中γr1表示C轴角度两组解中第一组解的真实角度,γr2表示C轴角度两组解中第二组解的真实角度,γpre表示上一个刀具姿态对应的加工代码中C轴角度,k1和k2为整数。
6.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤6中最短旋转轴角度路径原则为将得到的旋转轴角度两组解与上一个刀具姿态对应的加工代码中旋转轴角度比较,选择旋转轴角度路径变化最小的组合;旋转轴角度路径表达式为:min|α-αpre|+|γ-γpre|;其中αpre表示上一个刀具姿态对应的加工代码中A轴角度。
7.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤8中获得双转台五轴数控机床的机床坐标系下旋转轴的位置向量数值的方法为查询双转台五轴数控机床结构参数。
8.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤8中获得双转台五轴数控机床的机床坐标系下工件的位置向量的方法为安装好工件毛坯后通过机床对刀功能获得。
9.根据权利要求1所述的一种双转台五轴数控机床通用后处理方法,其特征在于:所述步骤5和所述步骤6中第一个刀具位置和刀具姿态的所述的上一个刀具姿态对应的加工代码为零。
CN201810282519.1A 2018-04-02 2018-04-02 一种双转台五轴数控机床通用后处理方法 Active CN108549319B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810282519.1A CN108549319B (zh) 2018-04-02 2018-04-02 一种双转台五轴数控机床通用后处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810282519.1A CN108549319B (zh) 2018-04-02 2018-04-02 一种双转台五轴数控机床通用后处理方法

Publications (2)

Publication Number Publication Date
CN108549319A CN108549319A (zh) 2018-09-18
CN108549319B true CN108549319B (zh) 2020-09-15

Family

ID=63513785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810282519.1A Active CN108549319B (zh) 2018-04-02 2018-04-02 一种双转台五轴数控机床通用后处理方法

Country Status (1)

Country Link
CN (1) CN108549319B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109839920B (zh) * 2019-03-18 2020-06-23 西南交通大学 一种五轴机床运动轴灵敏度分析方法
CN111624949B (zh) * 2020-06-03 2024-05-17 浙江日发航空数字装备有限责任公司 一种基于3轴程序编制5轴机床数控程序的方法
CN113759829B (zh) * 2021-08-26 2023-09-29 超同步股份有限公司 五轴机床带刀尖跟踪功能的后置处理方法及装置
CN114425666A (zh) * 2022-02-24 2022-05-03 北京工业大学 一种基于异形结构的六轴激光切割机床的后处理方法
CN117311268B (zh) * 2023-11-30 2024-02-23 北京天圣华信息技术有限责任公司 一种制作后处理文件方法及相关装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104965517A (zh) * 2015-07-07 2015-10-07 张耀伦 一种机器人笛卡尔空间轨迹的规划方法
CN107038275A (zh) * 2016-12-19 2017-08-11 中国科学院沈阳自动化研究所 一种机械臂误差分析方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980091B (zh) * 2010-08-23 2012-08-15 西安交通大学苏州研究院 双转台五轴联动数控加工旋转刀具中心点补偿方法
KR101199037B1 (ko) * 2010-12-02 2012-11-07 대영코어텍(주) 5축 cνc 가공기계를 이용한 드럼 캠의 정삭가공을 위한 후처리방법
CN103592888B (zh) * 2013-01-28 2016-09-07 广西工学院 国产双摆头五轴机床的后置处理方法
CN103862326B (zh) * 2014-03-11 2016-04-13 南京高传四开数控装备制造有限公司 一种双旋转工作台五轴机床旋转刀具刀尖控制方法
CN103878641B (zh) * 2014-03-14 2016-07-06 浙江大学 一种五轴数控机床通用的旋转轴几何误差辨识方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104965517A (zh) * 2015-07-07 2015-10-07 张耀伦 一种机器人笛卡尔空间轨迹的规划方法
CN107038275A (zh) * 2016-12-19 2017-08-11 中国科学院沈阳自动化研究所 一种机械臂误差分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Analysis of Multi-Product Manufacturing Systems with Homogeneous Exponential Machines;Cong Zhao and Jingshan Li;《2013 IEEE International Conference on Automation Science and Engineering (CASE)》;20131231;第386-391页 *

Also Published As

Publication number Publication date
CN108549319A (zh) 2018-09-18

Similar Documents

Publication Publication Date Title
CN108549319B (zh) 一种双转台五轴数控机床通用后处理方法
Fu et al. Product-of-exponential formulas for precision enhancement of five-axis machine tools via geometric error modeling and compensation
CN112558547B (zh) 一种五轴数控机床平动轴几何误差补偿数据快速优化方法
CN103048947B (zh) 五坐标联动机床数控程序后置处理方法
CN108829038B (zh) 一种刀尖运动轨迹控制算法
CN101980091A (zh) 双转台五轴联动数控加工旋转刀具中心点补偿方法
CN109664296B (zh) 一种面向机器人砂带打磨的优化轨迹搜索方法
CN105302070A (zh) 一种非正交摆头转台类五轴机床后置处理方法
CN104210302A (zh) 三维立体雕刻机、方法及系统
CN104460515A (zh) 一种基于后置处理五轴刀具长度补偿方法
Ding et al. Novel method for position-independent geometric error compensation of five-axis orthogonal machine tool based on error motion
CN108608425B (zh) 六轴工业机器人铣削加工离线编程方法及系统
CN111367237A (zh) 一种非正交五轴立卧转换数控机床后处理方法
CN104321707A (zh) 数控装置
CN114131611A (zh) 机器人重力位姿分解的关节误差离线补偿方法、系统及终端
Wang et al. Design and implementation of five-axis transformation function in CNC system
Xi et al. Kinematics for a six-axis EDM machine by screw theory and its application in feedrate planning in EDM for shrouded blisks
CN106020114B (zh) 一种数控机床综合误差可视化方法
CN109839920B (zh) 一种五轴机床运动轴灵敏度分析方法
CN111413924A (zh) 装万向角度头非正交五轴数控机床后处理方法
Kanda et al. Tool Path Generation for Five-Axis Controlled Machining with Consideration of Structural Interference.
CN108549323B (zh) 一种五轴数控工具磨床通用后置处理方法
Yun et al. A geometric postprocessing method for 5-axis machine tools using locations of joint points
JP2003005811A (ja) 6軸制御ncプログラム生成方法及び生成装置、並びに6軸制御ncプログラム生成プログラム及びこのプログラムを記憶したコンピュータ読み取り可能な記録媒体
CN109773581B (zh) 一种机器人应用于再现机加工的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201203

Address after: 610031 north section of two ring road, Sichuan, Chengdu

Patentee after: SOUTHWEST JIAOTONG University

Patentee after: BIJIE POWER SUPPLY BUREAU OF GUIZHOU POWER GRID Co.,Ltd.

Address before: 610031 Sichuan City, Chengdu Province, No. two North Ring Road, No. 111

Patentee before: SOUTHWEST JIAOTONG University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220117

Address after: 610031 north section of two ring road, Sichuan, Chengdu

Patentee after: SOUTHWEST JIAOTONG University

Address before: 610031 north section of two ring road, Sichuan, Chengdu

Patentee before: SOUTHWEST JIAOTONG University

Patentee before: Bijie Power Supply Bureau of Guizhou Power Grid Co., Ltd