CN110209123B - 一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法 - Google Patents
一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法 Download PDFInfo
- Publication number
- CN110209123B CN110209123B CN201910494747.XA CN201910494747A CN110209123B CN 110209123 B CN110209123 B CN 110209123B CN 201910494747 A CN201910494747 A CN 201910494747A CN 110209123 B CN110209123 B CN 110209123B
- Authority
- CN
- China
- Prior art keywords
- machine tool
- cutter
- axis
- vector
- coordinate system
- 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/416—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 control of velocity, acceleration or deceleration
-
- 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/35—Nc in input of data, input till input file format
- G05B2219/35215—Generate optimal nc program variant as function of cost, time, surface, energy
Abstract
本发明涉及切削加工领域,更具体的说是一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,可以解决机床加工过程中旋转轴的速度和加速度超过其电机最大速度和最大加速度的问题,充分发挥出五轴数控机床的加工特性,减少了机床加工过程产生的冲击,保证机床加工过程的稳定性,提高了机床的加工效率;并且解决了刀轴矢量全局优化模型计算效率低和难以求出满足机床驱动约束条件可行解的问题,提高了刀轴矢量优化的效率。
Description
技术领域
本发明涉及切削加工领域,更具体的说是一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法。
背景技术
五轴数控机床是加工航空发动机叶轮、汽车模具、光学零件等具有复杂曲面关键零件的重要设备。相较三轴数控机床,五轴数控机床增加了两个旋转轴,使刀具相对于工件的位置理论上任意可控,但因为增加了两个旋转轴,五轴数控机床的刀具轨迹规划变得更加复杂;特别是在加工复杂曲面刀轴矢量发生突变的区域,规划的刀具轨迹很可能会超过机床驱动电机的最大速度和最大加速度,产生冲击,导致加工不稳定,影响机床的加工效率和加工质量。以往的刀轴矢量优化方法是在机床驱动约束下对刀轴矢量进行全局优化,该方法计算效率较低,甚至无法求出满足机床驱动约束的可行解。
发明内容
本发明的目的是提供一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,该方法计算效率较高,而且满足机床驱动约束条件,保证机床加工过程的稳定性,对提高机床的加工效率具有重要意义。
本发明的目的通过以下技术方案来实现:
一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,包括以下步骤:
步骤一:建立加工曲面的驱动模型,加工曲面的刀具运动轨迹,根据机床结构将刀轴矢量在工件坐标系下转换到机床坐标系下,并通过差分计算出机床旋转轴的速度和加速度;
步骤二:对机床旋转轴的运动学特性进行分析,计算出旋转轴的速度超限区域,并根据其加速度小于设定阈值扩展成局部优化区域;
步骤三:以机床旋转轴最小的加速度为目标函数建立刀轴矢量局部优化模型,计算得到优化后的刀轴矢量。
作为本技术方案的进一步优化,本发明一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,所述步骤一中所述的将刀位点从工件坐标系下转换到机床坐标系下的公式为:
式中:Lcwx、Lcwy和Lcwz为工件在机床上的安装位置,x、y、z为刀尖点,i、j、k为刀轴矢量,X、Y、Z、A、C为机床各轴的位移。
作为本技术方案的进一步优化,本发明一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,所述步骤一中所述的刀轴矢量在工件坐标系下为:[i,j,k,0]=WTLCS(Trans(z,R)Rot(x,t)Rot(y,l))[0,0,1,0]T
式中:WTLCS是局部坐标系到工件坐标系的变换矩阵;x、y、z组成Tcpxyz坐标系,是局部坐标系从刀触点平移到刀具球心的坐标系,z与刀具与曲面接触点的法向矢量平行,x与刀具进给方向平行;l为前倾角,是刀轴矢量T在Tcpxyz坐标系下绕y轴旋转的角度;t为侧倾角,是刀轴矢量T在Tcpxyz坐标系下绕x轴旋转的角度,旋转角度的方向由右手螺旋定则确定。
作为本技术方案的进一步优化,本发明一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,所述步骤二中所述的确定局部优化区域的具体方法为:
步骤一:对机床旋转轴的运动学特性进行分析,计算出旋转轴的速度大于其驱动电机最大速度的区域vLimitm=[p,q],其中p为第p个刀位点,q为第q个刀位点;
步骤二:速度超限区域向两边扩展,扩展到其旋转轴的加速度小于设置的阈值athreshold,计算得到局部优化区域optSectionm。
作为本技术方案的进一步优化,本发明一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,所述步骤刀轴矢量局部优化模型为:
s.t.(li,ti)∈C_Space,i=p+2,p+3,...,q-3,q-2
(li,ti)=(linitial,tinitial),i=p,p+1,q-1,q
本发明一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法的有益效果为:
本发明一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,可以解决机床加工过程中旋转轴的速度和加速度超过其电机最大速度和最大加速度的问题,充分发挥出五轴数控机床的加工特性,减少了机床加工过程产生的冲击,保证机床加工过程的稳定性,提高了机床的加工效率;并且解决了刀轴矢量全局优化模型计算效率低和难以求出满足机床驱动约束条件可行解的问题,提高了刀轴矢量优化的效率。
附图说明
下面结合附图和具体实施方法对本发明做进一步详细的说明。
图1是本发明的刀轴矢量局部优化的流程框图;
图2是刀轴矢量在局部坐标系下的示意图;
图3是机床旋转轴的速度曲线的示意图;
图4是机床旋转轴的加速度曲线的示意图;
图5是本发明的刀轴矢量局部优化后机床旋转轴的速度曲线的示意图;
图6是本发明的刀轴矢量局部优化后机床旋转轴的加速度曲线的示意图;
图7是本发明的刀轴矢量局部优化后的刀具轨迹的示意图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
具体实施方式一:
下面结合图1-7说明本实施方式,一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,包括以下步骤:
步骤一:建立加工曲面的驱动模型,加工曲面的刀具运动轨迹,根据机床结构将刀轴矢量在工件坐标系下转换到机床坐标系下,并通过差分计算出机床旋转轴的速度和加速度;
步骤二:对机床旋转轴的运动学特性进行分析,计算出旋转轴的速度超限区域,并根据其加速度小于设定阈值扩展成局部优化区域;
步骤三:以机床旋转轴最小的加速度为目标函数建立刀轴矢量局部优化模型,计算得到优化后的刀轴矢量。
具体实施方式二:
下面结合图1-7说明本实施方式,本实施方式对实施方式一作进一步说明,所述步骤一中所述的将刀位点从工件坐标系下转换到机床坐标系下的公式为:
式中:Lcwx、Lcwy和Lcwz为工件在机床上的安装位置,x、y、z为刀尖点,i、j、k为刀轴矢量,X、Y、Z、A、C为机床各轴的位移。
具体实施方式三:
下面结合图1-7说明本实施方式,本实施方式对实施方式二作进一步说明,所述步骤一中所述的刀轴矢量在工件坐标系下为:
[i,j,k,0]=WTLCS(Trans(z,R)Rot(x,t)Rot(y,l))[0,0,1,0]T
式中:WTLCS是局部坐标系到工件坐标系的变换矩阵;x、y、z组成Tcpxyz坐标系,是局部坐标系从刀触点平移到刀具球心的坐标系,z与刀具与曲面接触点的法向矢量平行,x与刀具进给方向平行;l为前倾角,是刀轴矢量T在Tcpxyz坐标系下绕y轴旋转的角度;t为侧倾角,是刀轴矢量T在Tcpxyz坐标系下绕x轴旋转的角度,旋转角度的方向由右手螺旋定则确定。
具体实施方式四:
下面结合图1-7说明本实施方式,本实施方式对实施方式一作进一步说明,所述步骤二中所述的确定局部优化区域的具体方法为:
步骤一:对机床旋转轴的运动学特性进行分析,计算出旋转轴的速度大于其驱动电机最大速度的区域vLimitm=[p,q],其中p为第p个刀位点,q为第q个刀位点;
步骤二:速度超限区域向两边扩展,扩展到其旋转轴的加速度小于设置的阈值athreshold,计算得到局部优化区域optSectionm;若优化区域有相交情况,则对这些区域进行求并,如:optSectionm和optSectionm+1相交,则局部优化区域为optSectionm∪optSectionm+1。
具体实施方式五:
下面结合图1-7说明本实施方式,本实施方式对实施方式一至四任一项作进一步说明,所述步骤刀轴矢量局部优化模型为:
s.t.(li,ti)∈C_Space,i=p+2,p+3,...,q-3,q-2
(li,ti)=(linitial,tinitial),i=p,p+1,q-1,q
式中:l为刀具的前倾角,t为刀具的侧倾角,C_Space为刀轴矢量的可行域,和分别为机床J1旋转轴电机的最大速度和最大加速度,和分别为设定的机床J2旋转轴速度和加速度阈值,其中取值范围为取值范围为和分别为机床J2旋转轴电机的最大速度和最大加速度,linitial和tinitial分别为刀具前倾角l和侧倾角t的初值,其对局部优化区域[p,q]边界点p,p+1,q-1,q的刀具角度进行约束,保证刀轴矢量优化后,优化区域和非优化区域连接处机床旋转轴的速度和加速度小于其驱动电机的最大速度和大加速度。
本发明的一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,其工作原理为:
使用时可以采用AC双转台五轴数控机床,以直纹叶片模型为例,使用UG软件和Matlab软件说明本发明的实施过程,刀轴矢量局部优化的流程图见附图1;使用UG软件建立直纹叶片模型,以叶片的曲面为驱动曲面,生成该叶片的加工轨迹,并保存为刀位文件,根据机床结构将刀轴矢量在工件坐标系下转换到机床坐标系下,并通过差分计算出机床旋转轴的速度和加速度;
式中:Lcwx、Lcwy和Lcwz为工件在机床上的安装位置,x、y、z为刀尖点,i、j、k为刀轴矢量,X、Y、Z、A、C为机床各轴的位移。
刀轴矢量在刀具坐标系下的示意图见附图2,其在工件坐标系下为:
[i,j,k,0]=WTLCS(Trans(z,R)Rot(x,t)Rot(y,l))[0,0,1,0]T (2)
式中:WTLCS是局部坐标系到工件坐标系的变换矩阵;x、y、z组成Tcpxyz坐标系,是局部坐标系从刀触点平移到刀具球心的坐标系,z和刀具与曲面接触点的法向矢量平行,x和刀具进给方向平行;l为前倾角,是刀轴矢量T在Tcpxyz坐标系下绕y轴旋转的角度;t为侧倾角,是刀轴矢量T在Tcpxyz坐标系下绕x轴旋转的角度,旋转角度的方向由右手螺旋定则确定;根据公式(1)-(2),使用Matlab计算得到机床旋转轴的位移A、C;Matlab是MathWorks公司推出的一款商业数学软件;Matlab具有数值分析、数值和符号计算、工程与科学绘图、控制系统的设计与仿真、数字图像处理、数字信号处理、财务与金融工程等功能;对机床A、C轴位移进行差分得到A、C轴的速度和加速度:
式中:M为A轴或C轴的位移,fi是刀具的进给速率;根据步骤二对机床旋转轴进行运动学特性分析,其速度曲线和加速度曲线分别见附图3和附图4;其中,第一个超限区域对应的刀位点序列为[26,45],机床旋转轴驱动电机的最大速度为2.61799rad/s,最大加速度为3.49066rad/s2,速度超限区域向两边扩展,扩展到对应旋转轴的加速度小于0.01rad/s2,计算得到局部优化区optSection1=[10,58];建立刀轴矢量局部优化模型接着进行最优化模型的求解,求解后机床旋转轴的速度曲线和加速度曲线分别见附图5和附图6;其中,C轴的最大速度从10.7397rad/s下降到2.5029rad/s,C轴的最大加速度从92.323rad/s2下降到3.4907rad/s2,A轴的最大速度从0.0086rad/s下降到0.0012rad/s,A轴的最大加速度从1.1568rad/s2下降到0.0018rad/s2,所以优化后该区域满足机床驱动约束条件,使机床加工过程更加平稳,从而提高了复杂曲面的加工效率,优化后的刀具轨迹见附图7。
当然,上述说明并非对本发明的限制,本发明也不仅限于上述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也属于本发明的保护范围。
Claims (4)
1.一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,其特征在于:该方法包括以下步骤:
步骤一:建立加工曲面的驱动模型,加工曲面的刀具运动轨迹,根据机床结构将刀轴矢量在工件坐标系下转换到机床坐标系下,并通过差分计算出机床旋转轴的速度和加速度;
步骤二:对机床旋转轴的运动学特性进行分析,计算出旋转轴的速度超限区域,并根据其加速度小于设定阈值扩展成局部优化区域;
步骤三:以机床旋转轴最小的加速度为目标函数建立刀轴矢量局部优化模型,计算得到优化后的刀轴矢量;
所述步骤二中确定局部优化区域的具体方法为:
步骤一:对机床旋转轴的运动学特性进行分析,计算出旋转轴的速度大于其驱动电机最大速度的区域vLimitm=[p,q],其中p为第p个刀位点,q为第q个刀位点;
步骤二:速度超限区域向两边扩展,扩展到其旋转轴的加速度小于设置的阈值athreshold,计算得到局部优化区域optSectionm;
所述步骤三中刀轴矢量局部优化模型为:
s.t.(li,ti)∈C_Space,i=p+2,p+3,...,q-3,q-2
(li,ti)=(linitial,tinitial),i=p,p+1,q-1,q
3.根据权利要求2所述的一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法,其特征在于:所述步骤一中所述的刀轴矢量在工件坐标系下为:
[i,j,k,0]=WTLCS(Trans(z,R)Rot(x,t)Rot(y,l))[0,0,1,0]T
式中:WTLCS是局部坐标系到工件坐标系的变换矩阵;x、y、z组成Tcpxyz坐标系,是局部坐标系从刀触点平移到刀具球心的坐标系,z与刀具与曲面接触点的法向矢量平行,x与刀具进给方向平行;l为前倾角,是刀轴矢量T在Tcpxyz坐标系下绕y轴旋转的角度;t为侧倾角,是刀轴矢量T在Tcpxyz坐标系下绕x轴旋转的角度,旋转角度的方向由右手螺旋定则确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910494747.XA CN110209123B (zh) | 2019-06-10 | 2019-06-10 | 一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910494747.XA CN110209123B (zh) | 2019-06-10 | 2019-06-10 | 一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110209123A CN110209123A (zh) | 2019-09-06 |
CN110209123B true CN110209123B (zh) | 2020-05-29 |
Family
ID=67791492
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910494747.XA Active CN110209123B (zh) | 2019-06-10 | 2019-06-10 | 一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110209123B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112947296B (zh) * | 2021-03-10 | 2022-04-26 | 深圳市泰达智能装备有限公司 | 一种五轴速度规划方法、装置、计算机设备和存储介质 |
CN113359619B (zh) * | 2021-06-09 | 2022-09-16 | 意特利(上海)科技有限公司 | 一种立卧转换五轴机床的转台惯量自适应控制方法 |
CN114102262B (zh) * | 2021-12-27 | 2023-03-24 | 爱派尔(常州)数控科技有限公司 | 一种机床主轴安全检测方法、系统及机床 |
CN114217570B (zh) * | 2022-02-22 | 2022-06-14 | 成都飞机工业(集团)有限责任公司 | 一种五轴转角特征高效加工刀轨生成方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558808A (zh) * | 2013-09-28 | 2014-02-05 | 大连理工大学 | 复杂曲面五轴数控加工刀矢的运动学控制方法 |
CN104102171A (zh) * | 2014-07-04 | 2014-10-15 | 山东理工大学 | 一种球头刀多轴加工刀轴矢量优化方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10163503A1 (de) * | 2001-12-21 | 2003-07-10 | Siemens Ag | Polynom- und Spline-Interpolation von Werkzeugorientierungen |
KR100517880B1 (ko) * | 2003-11-20 | 2005-09-28 | 학교법인 포항공과대학교 | Cnc 가공에서 생산시간 단축과 가공 정밀도 향상을위한 이송속도 스케줄링 방법 |
JP5425342B1 (ja) * | 2012-05-22 | 2014-02-26 | 三菱電機株式会社 | 数値制御装置 |
CN103592891B (zh) * | 2013-09-28 | 2015-12-23 | 大连理工大学 | 运动学约束的复杂曲面五轴数控加工刀矢光顺方法 |
CN104635619A (zh) * | 2013-11-12 | 2015-05-20 | 沈阳高精数控技术有限公司 | 基于刀具矢量插补的五轴数控加工方法 |
CN105527927B (zh) * | 2016-01-29 | 2018-06-12 | 大连理工大学 | 一种基于机床各旋转轴角加速度分治优化的五轴加工刀轴矢量插值方法 |
CN106933189B (zh) * | 2017-04-26 | 2019-07-02 | 山东理工大学 | 基于ac型五轴数控机床环形刀加工刀轴矢量光顺方法 |
CN107608313B (zh) * | 2017-09-11 | 2019-09-27 | 大连理工大学 | 一种五轴双样条曲线插补速度规划方法 |
CN109375579B (zh) * | 2018-12-13 | 2021-02-02 | 武汉科技大学 | 一种基于运动学的五轴数控加工刀具姿态规划方法 |
-
2019
- 2019-06-10 CN CN201910494747.XA patent/CN110209123B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558808A (zh) * | 2013-09-28 | 2014-02-05 | 大连理工大学 | 复杂曲面五轴数控加工刀矢的运动学控制方法 |
CN104102171A (zh) * | 2014-07-04 | 2014-10-15 | 山东理工大学 | 一种球头刀多轴加工刀轴矢量优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110209123A (zh) | 2019-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110209123B (zh) | 一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法 | |
Ming et al. | Four-axis trochoidal toolpath planning for rough milling of aero-engine blisks | |
CN102147600B (zh) | 实时生成曲率连续路径的数控插补系统 | |
CN105676786B (zh) | 一种五轴数控加工中考虑各旋转轴角速度平滑特性的刀轴矢量插值方法 | |
CN104102171A (zh) | 一种球头刀多轴加工刀轴矢量优化方法 | |
CN104460516A (zh) | 一种基于后置处理五轴刀具半径补偿方法 | |
Meng et al. | Optimal barrel cutter selection for the CNC machining of blisk | |
CN109597357B (zh) | 一种面向叶片旋铣工艺的数控编程方法及装置 | |
CN112947298A (zh) | 一种机器人曲面加工轨迹优化生成方法、系统及终端 | |
CN102608952A (zh) | 对采用球头刀具的五轴联动机床平滑加工路径的方法 | |
JPWO2018179401A1 (ja) | 工具経路生成方法および装置 | |
CN110618659A (zh) | 一种面向五轴线性插补的平动轴旋转轴协调运动规划控制方法 | |
CN110362039B (zh) | 一种五轴加工工件摆放姿态优化方法 | |
Wang et al. | Reducing roughness of freeform surface through tool orientation optimization in multi-axis polishing of blisk | |
WO2010134532A1 (ja) | 数値制御装置 | |
CN103543691B (zh) | 适用于自由曲面高速高精加工的旋转轴位置优化方法 | |
Lu et al. | Smooth flank milling tool path generation for blade surfaces considering geometric constraints | |
CN113910021B (zh) | 异形外圆磨削方法及系统 | |
Lu et al. | Collision-free and smooth path generation method of five-axis ball-end milling considering effective tool diameter change | |
CN112147893B (zh) | 一种基于直纹面空间的五轴铣削刀轴矢量优化方法 | |
CN112464399B (zh) | 一种基于运动旋量的五轴加工刀具包络面计算方法 | |
Watanabe et al. | Development of tool collision avoidance method adapted to uncut workpiece shape | |
Qian et al. | Optimization of tool axis vector for mirror milling of thin-walled parts based on kinematic constraints | |
Zou et al. | Tool Orientation Optimization Based on Spatial Tractrix Method for Five-Axis CNC Machining with Ball End Cutters | |
CN113946136A (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 |