CN109375579B - 一种基于运动学的五轴数控加工刀具姿态规划方法 - Google Patents
一种基于运动学的五轴数控加工刀具姿态规划方法 Download PDFInfo
- Publication number
- CN109375579B CN109375579B CN201811525179.7A CN201811525179A CN109375579B CN 109375579 B CN109375579 B CN 109375579B CN 201811525179 A CN201811525179 A CN 201811525179A CN 109375579 B CN109375579 B CN 109375579B
- Authority
- CN
- China
- Prior art keywords
- cutter
- axis
- vector
- tool
- space
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical 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/19—Numerical 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 positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/45—Nc applications
- G05B2219/45136—Turning, lathe
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
本发明公开了一种基于运动学的五轴数控加工刀具姿态规划方法,将刀轴矢量定义为刀具轴线方向在工件坐标系下的矢量描述,采用同一刀路上刀位点PL,i处刀轴矢量vi与其相邻的刀位点PL,i‑1和PL,i+1处刀轴矢量vi‑1和vi+1的变化作为矢量描述的度量,具体采用相邻刀轴矢量之间夹角的平方之和然后开方作为度量指标,基于通用结构机床运动学模型和刀轴矢量光顺度量指标,即可建立考虑运动学空间映射的刀轴矢量整体光顺模型,完成刀具姿态规划,本发明大大提高多轴数控机床加工的精度和加工效率,减少刀具和工件之间的干扰,防止事故发生。
Description
技术领域
本发明属于机床加工领域,涉及一种刀具姿态规划方法,具体涉及一种基于运动学的五轴数控加工刀具姿态规划方法。
背景技术
五轴加工引入的两个旋转轴使得避免干涉、获得更好的表面质量和更高的加工效率成为可能。基于可行空间的规划可避免刀具与工件等的各种干涉碰撞,基于刀具偏离量子空间的规划可提高加工精度等。然而,旋转轴可实现更多的姿态也对刀具路径规划提出更高要求,单个刀位处的刀轴矢量规划可以追求尽可能高的加工效率或精度,而实际生产中,零件表面上离散出大量的网格刀位点,面临的新问题是在刀轴矢量规划时需要考虑相邻刀位点处的刀轴矢量之间的变化关系。刀轴矢量光顺性就是这样一个对工件表面质量、加工效率和加工过程稳定性有重要影响的规划要求。它反映的是整个加工表面上所有相邻刀位点处刀轴矢量的变化情况,因此对刀轴矢量光顺性进行度量建模即可完成对数控加工刀具姿态规划,然而现有技术中还缺乏相关成熟的技术。
发明内容
本发明的目的是提供一种基于运动学的五轴数控加工刀具姿态规划方法,解决现有技术中多轴机床工具姿态规划精度不高,控制困难的问题。
为了解决上述技术问题,本发明采用的技术方案是:
一种基于运动学的五轴数控加工刀具姿态规划方法,其特征在于:将刀轴矢量定义为刀具轴线方向在工件坐标系下的矢量描述,采用同一刀路上刀位点PL,i处刀轴矢量vi与其相邻的刀位点PL,i-1和PL,i+1处刀轴矢量vi-1和vi+1的变化作为矢量描述的度量,具体采用相邻刀轴矢量之间夹角的平方之和然后开方作为度量指标上述度量指标Δv1也即为刀轴矢量光顺度量指标;
广义刀具圆周包络曲面STs与工件实时曲面SW之间的干涉为全局干涉,二者之间保留一定的安全距离δTs,建立的刀轴矢量子空间表达式为:
VTs(Pi,j)={v|dmin(STs(Pi,j,v),SW)>δTs}
dmin表示刀具包络曲面与工件实时曲面的最小距离;
基于通用结构机床运动学模型和刀轴矢量光顺度量指标,即可建立考虑运动学空间映射的刀轴矢量整体光顺模型如下:
minΔv1(α,β)
s.t.
v(Pi,j,α,β)∈VTs(Pi,j)
通过上述模型即完成五轴数控加工刀具姿态规划,上式中,NL为同一刀路上刀位点总数,i为刀位点序列,α为刀位点表达参数,β为刀轴矢量表达参数,Pi,j为机床加工过程中,第j个刀路上的第i个刀位点,s.t为数学公式符号,是使得满足,v是刀轴矢量,VTs是刀轴矢量集合。
作为优选,考虑相邻刀路之间刀轴矢量变化情况,将上述度量指标Δv1更新为
Np为刀路总条数,刀轴矢量整体光顺模型更新为
minΔv2(α,β)
s.t.。
v(Pi,j,α,β)∈VTs(Pi,j)
作为优选,Δv1和Δv2均仅反映了工件坐标系下在刀具空间内的变化情况,要进行空间映射的运动学建模,将机床轴空间的变化情况纳入度量指标;
五轴机床轴空间的传动链包括的运动学元素主要有三个平动轴(x,y,z)和两个转动轴(a,c),可用PM表达为
PM=(x,y,z,a,c)T
PM映射到刀具空间内即为刀位点PL和刀轴矢量v,T表示转置矩阵;
以下建立两个空间的运动学映射关系,首先根据多体运动学的基本理论给出各种运动学变换的矩阵公式,设在五轴机床轴空间坐标系下绕x和z两个坐标轴分别旋转角度为a和c,及沿x、y和z三个坐标轴分别平移距离,这五种运动学变换元素的复合矩阵描述分别为:
将上述五种运动学变换元素合成,即同时绕两个坐标轴做两种旋转和沿各个坐标轴做三种平移,复合变换矩阵描述为
基于五轴机床轴空间与刀具空间的运动学映射关系,建立通用结构五轴机床的运动学映射模型,即可提出刀轴矢量变化映射到机床轴空间下的度量指标,将度量指标Δv2更新为:
刀轴矢量整体光顺模型更新为
minΔv3(α,β)
s.t.。
v(Pi,j,α,β)∈VTs(Pi,j)
作为优选,考虑相邻刀路之间的刀轴矢量变化在轴空间的映射,将度量指标Δv3更新为:
式中,Δxi,j、Δyi,j、Δzi,j、Δai,j和Δci,j依次为第i行第j列的刀位点对应的平动轴和转动轴解变化量;
刀轴矢量整体光顺模型更新为
minΔv4(α,β)
s.t.。
v(Pi,j,α,β)∈VTs(Pi,j)
5.如权利要求4所述的五轴数控加工刀具姿态规划方法,其特征在于:综合考虑Δv3和Δv1可以将度量指标Δv4更新为
Δv5=μW(Δv1)+μM(Δv3)
上式中,μW和μM分别为刀具空间和机床轴空间下矢量描述度量指标的权值,μW和μM均在[0,1]范围内取值,且μW+μM=1;
相应的刀轴矢量整体光顺模型更新为
minΔv5(α,β)
s.t.。
v(Pi,j,α,β)∈VTs(Pi,j)
本发明有益效果是:
本发明提供了一种刀轴矢量光顺性进行度量建模方法,可以进行多轴数控加工刀具姿态规划,通过本发明方法可以大大提高多轴数控机床加工的精度,提高加工效率,减少刀具和工件之间的干扰,防止事故发生。
附图说明
图1为本发明实施例中五轴机床轴空间和刀具空间映射关系示意图。
具体实施方式
本文将刀轴矢量定义为刀具轴线方向在工件坐标系下的矢量描述。将刀轴矢量定义为刀具轴线方向在工件坐标系下的矢量描述,采用同一刀路上刀位点PL,i处刀轴矢量vi与其相邻的刀位点PL,i-1和PL,i+1处刀轴矢量vi-1和vi+1的变化作为矢量描述的度量,具体采用相邻刀轴矢量之间夹角的平方之和然后开方作为度量指标,上述度量指标Δv1也即为刀轴矢量光顺度量指标;
公式(1)中,i表示刀位文件中同一刀路上的刀位点序列,NL为同一刀路上刀位点的总数。
现有技术表明相邻刀路之间较大的刀轴矢量变化会引起降低工件加工表面质量。基于此结论,考虑相邻刀路之间刀轴矢量的变化情况,将公式(1)中的指标更新为
公式(2)中,i表示刀位文件中第j条刀路上的刀位点序列,NL为同一刀路上刀位点的总数,NP为刀路总条数。
进一步地,Δv1和Δv2均仅反映了工件坐标系下在刀具空间内的变化情况,未能反映机床轴空间内因刀轴矢量变化引起的各轴运动量变化情况。为此,需要进行空间映射的运动学建模,将机床轴空间的变化情况纳入度量指标。
如图1所示,五轴机床轴空间的传动链包括的运动学元素主要有三个平动轴(x,y,z)和两个转动轴(a,c),可用PM表达为
PM=(x,y,z,a,c)T 公式(3)
它映射到刀具空间内即为刀位点PL和刀轴矢量v。
以下建立两个空间的运动学映射关系,首先根据多体运动学的基本理论给出各种运动学变换的矩阵公式,设在五轴机床轴空间坐标系下绕x和z两个坐标轴分别旋转角度为a和c,及沿x、y和z三个坐标轴分别平移距离,这五种运动学变换元素的复合矩阵描述分别为:
将这五种运动学变换元素合成,即同时绕两个坐标轴做两种旋转和沿各个坐标轴做三种平移,复合变换矩阵描述为
基于五轴机床轴空间与刀具空间的运动学映射关系,建立通用结构五轴机床的运动学映射模型,即可提出刀轴矢量变化映射到机床轴空间下的度量指标如下。
式中,Δxi、Δyi、Δzi、Δai和Δci分别是x、y、z、a和c五个轴的轴解变化量,μ1、μ2、μ3、μ4和μ5分别是Δxi、Δyi、Δzi、Δai和Δci的度量权重,在[0,1]内取值,且μ1+μ2+μ3+μ4+μ5=1,数值取值需要根据机床机构参数和工件特点进行设置。
考虑相邻刀路之间的刀轴矢量变化在轴空间的映射,公式(6)中度量指标v3更新为
公式(8)中,Δxi,j、Δyi,j、Δzi,j、Δai,j和Δci,j依次为第i行第j列的刀位点对应的各平动轴和转动轴的轴解变化量。
由公式(1)和公式(6),可得同时考虑刀具空间和机床轴空间内的刀轴矢量变化的度量指标为
Δv5=μW(Δv1)+μM(Δv3) 公式(10)
式中,μW和μM分别为刀具空间和机床轴空间下度量指标的权值,均在[0,1]范围内取值,且μW+μM=1。
公式(1)、公式(2)、公式(6)、公式(8)和公式(10)分别给出了不同考虑下的刀轴矢量光顺度量指标。对于大型曲面离散点数较多且对光顺的要求不是很高时,可选择公式(1)仅考虑同一刀路上相邻刀位点处刀轴矢量的几何光顺。对于结构复杂的机床,由刀具空间到机床轴空间的映射计算量较大,此时,可根据需要仅考虑刀具空间的几何光顺,而如果对几何整体光顺性要求高,可同时考虑相邻刀路上刀轴矢量的光顺,选择公式(2)作为光顺度量指标。对于机床运动轴的光顺的必要性较大时,比如复杂曲面工件高速加工,需要控制各轴运动量的变化,则选择公式(6)、公式(8)和公式(10)对应的包含运动学光顺的度量指标。具体选择还需根据机床结构和尺寸特点、工件设计曲面特点、加工工艺、切削参数、公差要求和任务时间等进行综合考量。
整体光顺建模与算例
基于通用结构机床运动学模型和整体光顺度量指标,即可建立考虑运动学空间映射的刀轴矢量整体光顺模型如下
minΔv5(α,β)
s.t. 公式(11)
v(Pi,j,α,β)∈VTs(Pi,j)
其中,优化目标为同时考虑刀具空间和机床轴空间内的刀轴矢量变化及相邻刀路之间的的刀轴矢量变化的度量指标,由公式(10)确定。可根据需要进行调整,比如只考虑单独空间,或只考虑所在刀路,可将其简化为公式公式(1)、公式(2)、公式(6)和公式(8)中的指标,其中约束为考虑干涉避免的子空间约束。亦可根据需要添加上一章中提出的考虑刀具偏离量与工件表面成形误差的子空间约束。
考虑干涉避免的刀轴矢量子空间
广义刀具圆周包络曲面STs与工件实时曲面SW之间的干涉,这种干涉简称为全局干涉。此处,广义刀具圆周曲面是指刀具、刀柄与主轴头部分组成的包络体曲面,工件实时曲面为加工过程中的判断瞬时工件的加工后曲面。避免全局干涉,需要控制广义刀具圆周曲面与工件实时曲面之间的距离,使二者之间保留一定的安全距离δTs,建立的刀轴矢量子空间表达式为
VTs(Pi,j)={v|dmin(STs(Pi,j,v),SW)>δTs} 公式(12)
dmin表示刀具包络曲面与工件实时曲面的最小距离。
Claims (3)
1.一种基于运动学的五轴数控加工刀具姿态规划方法,其特征在于:将刀轴矢量定义为刀具轴线方向在工件坐标系下的矢量描述,采用同一刀路上刀位点PL,i处刀轴矢量vi与其相邻的刀位点PL,i-1和PL,i+1处刀轴矢量vi-1和vi+1的变化作为矢量描述的度量,具体采用相邻刀轴矢量之间夹角的平方之和然后开方作为度量指标上述度量指标Δv1也即为刀轴矢量光顺度量指标;
广义刀具圆周包络曲面STs与工件实时曲面SW之间的干涉为全局干涉,二者之间保留一定的安全距离δTs,建立的刀轴矢量子空间表达式为:
VTs(Pi,j)={v|dmin(STs(Pi,j,v),SW)>δTs}
dmin表示刀具包络曲面与工件实时曲面的最小距离;
基于通用结构机床运动学模型和刀轴矢量光顺度量指标,即可建立考虑运动学空间映射的刀轴矢量整体光顺模型如下:
minΔv1(α,β)
s.t.
v(Pi,j,α,β)∈VTs(Pi,j)
通过上述模型即完成五轴数控加工刀具姿态规划,上式中,NL为同一刀路上刀位点总数,i为刀位点序列,α为刀位点表达参数,β为刀轴矢量表达参数,Pi,j为机床加工过程中,第j个刀路上的第i个刀位点,s.t为数学公式符号,是使得满足,v是刀轴矢量,VTs是刀轴矢量集合;
考虑相邻刀路之间刀轴矢量变化情况,将上述度量指标Δv1更新为
Np为刀路总条数,刀轴矢量整体光顺模型更新为
minΔv2(α,β)
s.t.
v(Pi,j,α,β)∈VTs(Pi,j)
Δv1和Δv2均仅反映了工件坐标系下在刀具空间内的变化情况,要进行空间映射的运动学建模,将机床轴空间的变化情况纳入度量指标;
五轴机床轴空间的传动链包括的运动学元素主要有三个平动轴(x,y,z)和两个转动轴(a,c),可用PM表达为
PM=(x,y,z,a,c)T
PM映射到刀具空间内即为刀位点PL和刀轴矢量v,T表示转置矩阵;
以下建立两个空间的运动学映射关系,首先根据多体运动学的基本理论给出各种运动学变换的矩阵公式,设在五轴机床轴空间坐标系下绕x和z两个坐标轴分别旋转角度为a和c,及沿x、y和z三个坐标轴分别平移距离,这五种运动学变换元素的复合矩阵描述分别为:
将上述五种运动学变换元素合成,即同时绕两个坐标轴做两种旋转和沿各个坐标轴做三种平移,复合变换矩阵描述为
基于五轴机床轴空间与刀具空间的运动学映射关系,建立通用结构五轴机床的运动学映射模型,即可提出刀轴矢量变化映射到机床轴空间下的度量指标,将度量指标Δv2更新为:
上式中,Δxi、Δyi、Δzi、Δai和Δci分别是x、y、z、a和c五个轴的轴解变化量,μ1、μ2、μ3、μ4和μ5分别是Δxi、Δyi、Δzi、Δai和Δci的度量权重,在[0,1]内取值,且μ1+μ2+μ3+μ4+μ5=1;
刀轴矢量整体光顺模型更新为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811525179.7A CN109375579B (zh) | 2018-12-13 | 2018-12-13 | 一种基于运动学的五轴数控加工刀具姿态规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811525179.7A CN109375579B (zh) | 2018-12-13 | 2018-12-13 | 一种基于运动学的五轴数控加工刀具姿态规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109375579A CN109375579A (zh) | 2019-02-22 |
CN109375579B true CN109375579B (zh) | 2021-02-02 |
Family
ID=65373597
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811525179.7A Active CN109375579B (zh) | 2018-12-13 | 2018-12-13 | 一种基于运动学的五轴数控加工刀具姿态规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109375579B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110209123B (zh) * | 2019-06-10 | 2020-05-29 | 哈尔滨工业大学 | 一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法 |
CN111665784B (zh) * | 2020-05-15 | 2021-08-03 | 成都飞机工业(集团)有限责任公司 | 一种基于西门子系统的空间定位误差补偿方法 |
CN111857047B (zh) * | 2020-06-15 | 2021-11-12 | 广东长盈精密技术有限公司 | 四轴联动加工方法、装置、计算机设备和存储介质 |
CN112904796B (zh) * | 2021-01-15 | 2022-05-06 | 西北工业大学 | 一种五轴铣削中刀具底刃刮蹭现象预测方法 |
CN113189937B (zh) * | 2021-04-28 | 2023-03-24 | 成都永峰科技有限公司 | 航空类零件自动生产线刀具集成管理方法及其系统和应用 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103592891B (zh) * | 2013-09-28 | 2015-12-23 | 大连理工大学 | 运动学约束的复杂曲面五轴数控加工刀矢光顺方法 |
CN105425727B (zh) * | 2015-12-08 | 2018-11-16 | 上海交通大学 | 五轴侧铣加工刀具路径光顺方法 |
CN106933190B (zh) * | 2017-04-26 | 2019-07-02 | 山东理工大学 | 基于bc型五轴数控机床环形刀加工刀轴矢量光顺方法 |
CN108415374B (zh) * | 2018-02-01 | 2019-08-13 | 大连理工大学 | 基于机床旋转进给轴运动学特性的刀轴矢量光顺方法 |
-
2018
- 2018-12-13 CN CN201811525179.7A patent/CN109375579B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109375579A (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109375579B (zh) | 一种基于运动学的五轴数控加工刀具姿态规划方法 | |
Wang et al. | Trajectory planning and optimization for robotic machining based on measured point cloud | |
CN111367236B (zh) | 一种面向加工过程的移动机器人系统标定方法及系统 | |
CN111702762B (zh) | 一种工业机器人作业姿态优化方法 | |
CN109765848B (zh) | 五轴机床运动轴几何误差贡献值计算及影响评估方法 | |
Liao et al. | Optimization of robot posture and workpiece setup in robotic milling with stiffness threshold | |
CN106873522B (zh) | 一种非轴对称扫掠曲面的数控车削刀具轨迹规划方法 | |
CN111958611B (zh) | 考虑轮廓误差最小的铣削加工机器人姿态优化方法及设备 | |
Li et al. | Interference-free inspection path generation for impeller blades using an on-machine probe | |
CN105467927A (zh) | 一种基于配准对齐的免划线锻铸件加工定位方法 | |
CN109176487A (zh) | 一种联动关节段标定方法、系统、设备、存储介质 | |
Li et al. | Optimal measurement configurations for kinematic calibration of six-DOF serial robot | |
CN110989490A (zh) | 一种基于轮廓误差的工件最优安装位置的获取方法 | |
CN111452033A (zh) | 一种工业机器人双nurbs曲线铣削轨迹规划方法 | |
CN111679629A (zh) | 一种多主轴头加工的空行程无干涉轨迹规划方法 | |
CN115113582A (zh) | 一种面向零件轮廓误差的五轴位置环增益离线矫正方法 | |
CN113778018B (zh) | 一种基于R-test的五轴机床刀轴矢量误差测量方法 | |
Guo et al. | A robotic boring system for intersection holes in aircraft assembly | |
CN113917888B (zh) | 一种基于固定角向标定及补偿的加工精度提升方法 | |
CN113985813B (zh) | 一种基于在机检测的机床原点误差补偿方法 | |
CN113325802B (zh) | 一种五轴机床几何误差补偿方法 | |
Yu et al. | Post-processing algorithm of a five-axis machine tool with dual rotary tables based on the TCS method | |
CN106959665B (zh) | 基于大数据的机床产品多权值融合的几何精度维护方法 | |
CN116834009A (zh) | 一种分布式虚实孪生的机器人位姿误差在线补偿系统 | |
CN108549323B (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: 20220526 Address after: 430080 No. 1212, iron and steel building, Wuhan University of science and technology, No. 743, Heping Avenue, Qingshan District, Wuhan City, Hubei Province Patentee after: Wuhan Wuke servo Electromechanical Equipment Co.,Ltd. Address before: 430081 No. 947 Heping Avenue, Qingshan District, Hubei, Wuhan Patentee before: WUHAN University OF SCIENCE AND TECHNOLOGY |