CN112883502A - *速度曲线的设计方法及基于*速度曲线的五轴轨迹加工方法 - Google Patents

*速度曲线的设计方法及基于*速度曲线的五轴轨迹加工方法 Download PDF

Info

Publication number
CN112883502A
CN112883502A CN201911203626.1A CN201911203626A CN112883502A CN 112883502 A CN112883502 A CN 112883502A CN 201911203626 A CN201911203626 A CN 201911203626A CN 112883502 A CN112883502 A CN 112883502A
Authority
CN
China
Prior art keywords
speed
velocity
track
critical
curve
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
CN201911203626.1A
Other languages
English (en)
Other versions
CN112883502B (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.)
Xian Jiaotong University
Original Assignee
Xian 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201911203626.1A priority Critical patent/CN112883502B/zh
Publication of CN112883502A publication Critical patent/CN112883502A/zh
Application granted granted Critical
Publication of CN112883502B publication Critical patent/CN112883502B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)

Abstract

本发明公开了
Figure DDA0002296462940000011
速度曲线的设计方法及基于
Figure DDA0002296462940000012
速度曲线的五轴轨迹加工方法,本发明的方法可以将S2速度曲线的时间延长至任意值的速度规划方法;基于此,根据五轴加工轨迹的刀尖和刀轴方向子轨迹的运动学约束生成各子轨迹的最短运动时间,并通过
Figure DDA0002296462940000014
速度曲线调整较快子轨迹的最大速度或末端速度,实现两子轨迹的运动时间同步,最终实现五轴机床在给定的时间内沿指令轨迹运动到指定位置和位姿。本发明的方法综合考虑了刀尖和刀轴方向子轨迹的运动学约束,通过提出的
Figure DDA0002296462940000013
速度曲线能够实现在相同时间内沿指令轨迹运动到指定位置和刀具位姿,避免了任何子轨迹的驱动约束饱和,获得了连续平稳的运动过程。

Description

*速度曲线的设计方法及基于*速度曲线的五轴轨迹加工 方法
技术领域
本发明属于运动规划技术领域,具体涉及
Figure BDA0002296462920000013
速度曲线的设计方法及基于
Figure BDA0002296462920000014
速度曲线的五轴轨迹加工方法。
背景技术
五轴机床轨迹由刀尖子轨迹和刀轴方向子轨迹组成,其中刀尖子轨迹位于笛卡尔平面坐标系而刀轴方向子轨迹位于球面坐标系,两轨迹所处空间是独立的,实现刀尖和刀轴子轨迹的运动同步是五轴插补中的关键技术之一。文献1“Beudaert X,Lavernhe S,Tournier C.5-axis local corner rounding of linear tool path discontinuities[J].International Journal of Machine Tools and Manufacture,2013,73:9-16.”和文献2“Tulsyan S,Altintas Y.Local toolpath smoothing for five-axis machine tools[J].International Journal of Machine Tools and Manufacture,2015,96:15-26.”采用了线性参数同步的方法,通过将刀尖子轨迹作为主运动,刀轴方向子轨迹作为从运动,并通过参数同时实现两轨迹的同步运动。速度规划过程中,先依据刀尖运动学约束对刀尖子轨迹进行速度规划生成刀尖速度曲线,并通过参数同步生成刀轴方向速度曲线,但是,该过程忽略了刀轴方向的运动学约束。在某些情况下,线性参数同步方法会造成刀轴子轨迹的速度规划违背其运动学约束,造成机床物理轴轴的饱和,进而破坏五轴加工质量。
发明内容
针对现有技术中的技术问题,本发明提供了
Figure BDA0002296462920000015
速度曲线的设计方法及基于
Figure BDA0002296462920000016
速度曲线的五轴轨迹加工方法,综合考虑了刀尖和刀轴方向子轨迹的运动学约束,通过提出的
Figure BDA0002296462920000017
速度曲线能够实现在相同时间内沿指令轨迹运动到指定位置和刀具位姿。
为解决上述技术问题,本发明通过以下技术方案予以实现:
一种
Figure BDA0002296462920000026
速度曲线的设计方法,包括以下步骤:
步骤1:确定
Figure BDA0002296462920000027
速度曲线的匀速速度上限值;
步骤2:计算临界时间,用于在步骤3中确定
Figure BDA0002296462920000028
速度曲线的类型;
步骤3:根据步骤1确定的匀速速度上限值和步骤2计算的临界时间确定
Figure BDA0002296462920000029
速度曲线类型。
进一步地,所述步骤1的具体方法为:
给定参数,所述参数包括目标位移L、初始速度vs、终止速度ve、匀速速度vm、加速度极限值Amax和跃度极限值Jmax;根据给定的所述参数假设运动过程包括加速阶段、匀速阶段和减速阶段,其中,加速阶段的持续时间ta为:
Figure BDA0002296462920000021
减速阶段的持续时间td为:
Figure BDA0002296462920000022
Figure BDA0002296462920000023
时,假设成立,运动过程存在匀速阶段,且匀速速度为vm
Figure BDA0002296462920000024
时,假设不成立,运动过程包括加速阶段和减速阶段,且速度曲线的最大速度无法取匀速速度vm,实际最大速度vma借助Newtown-Raphson方法在max(vs,ve)和vm之间寻找,直至匀速阶段的持续时间tc介于0与插补周期Ts之间,迭代过程结束,匀速阶段的持续时间tc为0,匀速速度vm为实际最大速度vma
Figure BDA0002296462920000025
本步骤获取的匀速速度vm
Figure BDA00022964629200000210
速度曲线的匀速速度上限值vmm
进一步地,所述步骤2的具体方法为:
定义临界时间T1,T1为在给定的初始速度vs、终止速度ve和目标位移L下,采用加速度极限值Amax和跃度极限值Jmax进行加速运动和匀速运动的最短持续时间,表示为:
Figure BDA0002296462920000031
其中,T1a为加速时间,表示为
Figure BDA0002296462920000032
定义临界时间T2,T2为仅进行加速运动时,从初始速度vs到终止速度ve的时间,表示为:
Figure BDA0002296462920000033
定义临界时间T3,T3为仅以初始速度vs进行匀速运动走完目标位移L的时间,表示为:
Figure BDA0002296462920000034
进一步地,所述步骤3的具体方法为:
1)若目标时间Te小于临界时间T1,则
Figure BDA0002296462920000039
速度曲线包含加速阶段、匀速阶段和减速阶段,定义为类型A1,具体如下:
加速阶段的持续时间ta为:
Figure BDA0002296462920000035
减速阶段的持续时间td为:
Figure BDA0002296462920000036
匀速阶段的持续时间tc为Te-ta-td,匀速速度vm为:
Figure BDA0002296462920000037
Figure BDA0002296462920000038
速度曲线的指令位移为:
Figure BDA0002296462920000041
2)若T1<Te<T2
Figure BDA0002296462920000042
速度曲线包含加速阶段和匀速阶段,定义为类型A2,具体如下:加速阶段的持续时间ta为:
Figure BDA0002296462920000043
匀速阶段的持续时间tc为Te-ta,终止速度ve为:
Figure BDA0002296462920000044
Figure BDA0002296462920000045
速度曲线的指令位移为:
Figure BDA0002296462920000046
3)若T2<Te<T3
Figure BDA0002296462920000047
速度曲线包含加速阶段,定义为类型A3,具体如下:
加速阶段的持续时间ta为Te,匀速阶段的终止速度ve为:
Figure BDA0002296462920000048
Figure BDA0002296462920000049
速度曲线的指令位移为:
Figure BDA00022964629200000410
4)若T3<Te
Figure BDA00022964629200000411
速度曲线先考虑先减速阶段后匀速阶段,定义为类型A4,具体如下:减速阶段的持续时间td为:
Figure BDA00022964629200000412
匀速阶段的持续时间tc为Te-td,匀速阶段的终止速度ve为:
Figure BDA0002296462920000051
当ve≥0时,
Figure BDA00022964629200000510
速度曲线的指令位移为:
Figure BDA0002296462920000052
若ve<0时,
Figure BDA0002296462920000053
速度曲线需要进一步修改成减速阶段I、匀速阶段和减速阶段II,定义为类型A5,此时,采用Newtown-Raphson方法寻找保证ve≥0的匀速速度vm,并确定各阶段的持续时间,匀速速度vm的搜索范围为0和vs,终止条件为|Tev-Te|≤Ts
Figure BDA0002296462920000054
其中,td1和td2为减速阶段I和减速阶段II的持续时间,Tev为整个变化过程的时间,为实现变化时间等于目标时间Te,匀速速度vm为:
Figure BDA0002296462920000055
Figure BDA0002296462920000056
速度曲线的指令位移为:
Figure BDA0002296462920000057
基于
Figure BDA0002296462920000058
速度曲线的五轴轨迹加工方法,
Figure BDA0002296462920000059
速度曲线应用所述的设计方法;
预处理阶段:对五轴加工轨迹进行预处理,待插补的五轴加工轨迹包含刀尖子轨迹和刀轴方向子轨迹;首先,计算各子轨迹的曲率,根据曲率极值将各子轨迹分成离散的N个曲线段,曲率极值点为曲线段的衔接点,定义为临界点;然后,计算各曲线段的弧长;将临界点坐标、曲率极值和弧长存储在计算机中;
双向极值扫描阶段:对预处理后的五轴加工轨迹进行双向极值扫描,通过双向极值扫描更新临界点速度,双向极值扫描阶段包含向后看阶段和向前看阶段,扫描过程包含N个曲线段;向后看从轨迹终止点开始向起点反向扫描,向前看从轨迹起始点开始向终点正向扫描;
向后看的具体步骤如下:
1)令i=N,向后看获取的第i个临界点处的刀尖速度极值
Figure BDA0002296462920000061
刀轴方向速度极值
Figure BDA0002296462920000062
2)若i=1,执行3);
否则,令当前模块的刀尖初始速度
Figure BDA0002296462920000063
刀轴方向初始速度
Figure BDA0002296462920000064
执行2a);
2a)结合预处理阶段计算的刀尖子轨迹第i段曲线的轨迹长度
Figure BDA0002296462920000065
和切向跃度
Figure BDA0002296462920000066
将刀尖临界点速度记为vsj,由下式计算临界点速度vsj
Figure BDA0002296462920000067
结合刀尖子轨迹的切向加速度
Figure BDA00022964629200000619
将刀尖临界点速度记为vsa,vsa满足如下不等式:
Figure BDA0002296462920000068
结合预处理阶段计算的刀尖子轨迹第i段曲线的临界点曲率极值
Figure BDA0002296462920000069
预设的弓差误差δ、预设的法向加速度
Figure BDA00022964629200000610
和预设的法向跃度
Figure BDA00022964629200000611
将刀尖临界点速度记为ven,其表达为:
Figure BDA00022964629200000612
其中,Ts为插补周期;
刀尖子轨迹的临界点速度
Figure BDA00022964629200000613
由切向约束和法向约束共同决定,表示为:
Figure BDA00022964629200000614
其中,vm为刀尖匀速速度;
2b)根据第i段刀尖子轨迹的初始速度
Figure BDA00022964629200000615
终止速度
Figure BDA00022964629200000616
轨迹长度
Figure BDA00022964629200000617
切向跃度
Figure BDA00022964629200000618
和切向加速度
Figure BDA0002296462920000076
结合
Figure BDA0002296462920000077
速度曲线的生成方法确定最短持续时间
Figure BDA0002296462920000078
2c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的轨迹长度
Figure BDA0002296462920000079
和切向跃度
Figure BDA00022964629200000710
将刀轴方向临界点速度记为ωsj,由下式计算临界点速度ωsj
Figure BDA0002296462920000071
结合刀轴方向子轨迹的切向加速度
Figure BDA00022964629200000711
将刀轴方向临界点速度记为ωsa,ωsa满足如下不等式:
Figure BDA0002296462920000072
结合预处理阶段计算的刀轴方向子轨迹第i段曲线的临界点曲率极值
Figure BDA00022964629200000712
法向加速度
Figure BDA00022964629200000713
和法向跃度
Figure BDA00022964629200000714
将刀轴方向临界点速度记为ωsn,其表示为:
Figure BDA0002296462920000073
刀轴方向临界点速度
Figure BDA00022964629200000715
由切向约束和法向约束共同决定,表示为:
Figure BDA0002296462920000074
其中,ωm为刀轴方向匀速速度;
2d)根据第i段刀轴方向子轨迹的初始速度
Figure BDA00022964629200000716
终止速度
Figure BDA00022964629200000717
轨迹长度
Figure BDA00022964629200000718
切向跃度
Figure BDA00022964629200000719
和切向加速度
Figure BDA00022964629200000720
结合
Figure BDA00022964629200000721
速度曲线的生成方法确定最短持续时间
Figure BDA00022964629200000722
2e)若
Figure BDA00022964629200000723
将刀轴方向子轨迹以
Figure BDA00022964629200000724
速度曲线修改并延长其持续时间至
Figure BDA00022964629200000725
先计算刀轴方向子轨迹的
Figure BDA00022964629200000726
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200000727
与临界时间比较并确定类型,记录修改后的临界速度为
Figure BDA00022964629200000728
Figure BDA00022964629200000729
将刀尖子轨迹以
Figure BDA00022964629200000730
速度曲线修改并延长其持续时间至
Figure BDA00022964629200000731
先计算刀尖子轨迹的
Figure BDA00022964629200000732
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200000734
与临界时间比较并确定类型,记录修改后的临界速度为
Figure BDA00022964629200000733
2f)令
Figure BDA0002296462920000075
i=i-1,执行2);
3)将
Figure BDA0002296462920000087
Figure BDA0002296462920000088
存储在数控单元缓存中;
向前看的具体步骤如下:
11)令i=1,向前看获取的第i个临界点处的刀尖速度极值vi=0,刀轴方向速度极值ωi=0;
22)若i=N,执行33);
否则,令当前模块的刀尖初始速度
Figure BDA0002296462920000089
刀轴方向初始速度
Figure BDA00022964629200000810
执行22a);
22a)结合预处理阶段计算的刀尖子轨迹第i段曲线的轨迹长度
Figure BDA00022964629200000811
和切向跃度
Figure BDA00022964629200000812
将刀尖临界点速度记为vej,由下式计算临界点速度vej
Figure BDA0002296462920000081
结合刀尖子轨迹的切向加速度
Figure BDA00022964629200000813
将刀尖临界点速度记为vea,vea满足如下不等式
Figure BDA0002296462920000082
刀尖子轨迹的临界点速度
Figure BDA00022964629200000814
由切向约束和向后看的临界速度共同决定,表示为:
Figure BDA0002296462920000083
22b)根据第i段刀尖子轨迹的初始速度
Figure BDA00022964629200000815
终止速度
Figure BDA00022964629200000816
轨迹长度
Figure BDA00022964629200000817
切向跃度
Figure BDA00022964629200000818
和切向加速度
Figure BDA00022964629200000819
结合
Figure BDA00022964629200000820
速度曲线的生成方法确定最短持续时间
Figure BDA00022964629200000821
22c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的轨迹长度
Figure BDA00022964629200000822
和切向跃度
Figure BDA00022964629200000823
将刀轴方向临界点速度记为ωej,由下式计算临界点速度ωsj
Figure BDA0002296462920000084
结合刀轴方向子轨迹的切向加速度
Figure BDA00022964629200000824
将刀轴方向临界点速度记为ωea,ωsa满足如下不等式:
Figure BDA0002296462920000085
刀轴方向子轨迹的临界点速度
Figure BDA00022964629200000825
由切向约束和向后看的临界速度共同决定,表示为
Figure BDA0002296462920000086
22d)根据第i段刀轴方向子轨迹的初始速度
Figure BDA0002296462920000091
终止速度
Figure BDA0002296462920000092
轨迹长度
Figure BDA0002296462920000093
切向跃度
Figure BDA0002296462920000094
和切向加速度
Figure BDA0002296462920000095
结合
Figure BDA0002296462920000096
速度曲线的生成方法确定最短持续时间
Figure BDA0002296462920000097
22e)若
Figure BDA0002296462920000098
将刀轴方向子轨迹以
Figure BDA0002296462920000099
速度曲线修改并延长其持续时间至
Figure BDA00022964629200000910
先计算刀轴方向子轨迹的
Figure BDA00022964629200000911
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200000912
与临界时间比较并确定类型,记录修改后的临界速度
Figure BDA00022964629200000913
和最大速度v′m
Figure BDA00022964629200000914
将刀尖子轨迹以
Figure BDA00022964629200000915
速度曲线修改并延长其持续时间至
Figure BDA00022964629200000916
先计算刀尖子轨迹的
Figure BDA00022964629200000917
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200000918
与临界时间比较并确定类型,记录修改后的临界速度
Figure BDA00022964629200000919
和最大速度ω′m
22f)令
Figure BDA00022964629200000920
i=i+1,执行22);
33)存储vi,ωi
Figure BDA00022964629200000921
和Ti(i=1,...,N)至数控单元缓存。
与现有技术相比,本发明至少具有以下有益效果:本发明的方法可以在满足速度、加速度和跃度极限值的前提下实现维持目标时间的速度规划;基于此,根据五轴加工轨迹的刀尖和刀轴方向子轨迹的运动学约束生成各子轨迹的最短运动时间,并通过
Figure BDA00022964629200000922
速度曲线调整较快子轨迹的最大速度或末端速度,实现两子轨迹的运动时间同步,最终实现五轴机床在给定的时间内沿指令轨迹运动到指定位置和位姿。本发明的方法将刀尖和刀轴方向子轨迹看作地位相等的运动,通过本发明提出的
Figure BDA00022964629200000923
速度曲线能够实现在相同时间内沿指令轨迹运动到指定位置和刀具位姿,避免了任何子轨迹的驱动约束饱和,且实现了刀尖和刀轴方向子轨迹的运动同步,获得了连续平稳的运动过程。本发明具有计算效率高和柔性运动能力优势,更加适用于五轴轨迹速度规划算法的在线应用。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式中的技术方案,下面将对具体实施方式描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1说明采用S2速度曲线生成的速度曲线示意图;
图2说明四种类型的
Figure BDA0002296462920000102
速度曲线示意图;
图3说明本实施例对应的工件坐标系下的刀尖子轨迹和机床坐标系下的刀轴方向子轨迹数据图;
图4说明时间同步算法前后的各位置点对应的速度曲线图;;
图5说明参数线性同步和时间同步算法对应的刀轴方向速度曲线图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
作为本发明的某一具体实施方式,一种
Figure BDA0002296462920000103
速度曲线的设计方法,包括以下步骤:
步骤1:确定
Figure BDA0002296462920000104
速度曲线的匀速速度上限值,具体方法为:
给定参数,所述参数包括目标位移L、初始速度vs、终止速度ve、匀速速度vm、加速度极限值Amax和跃度极限值Jmax;根据给定的所述参数假设运动过程包括加速阶段、匀速阶段和减速阶段,其中,加速阶段的持续时间ta为:
Figure BDA0002296462920000101
减速阶段的持续时间td为:
Figure BDA0002296462920000111
Figure BDA0002296462920000112
时,假设成立,运动过程存在匀速阶段,且匀速速度为vm
Figure BDA0002296462920000113
时,假设不成立,运动过程包括加速阶段和减速阶段,且速度曲线的最大速度无法取匀速速度vm,实际最大速度vma借助Newtown-Raphson方法在max(vs,ve)和vm之间寻找,直至匀速阶段的持续时间tc介于Ts与0之间,迭代过程结束,匀速阶段的持续时间tc为0,匀速速度vm为实际最大速度vma
Figure BDA0002296462920000114
此时,满足加速度极限值和跃度极限值获得的最短持续时间为ta+tc+td
本步骤获取的匀速速度vm
Figure BDA0002296462920000119
速度曲线的匀速速度上限值vmm
步骤2:计算临界时间,用于在步骤3中确定
Figure BDA00022964629200001110
速度曲线的类型,具体方法为:
定义临界时间T1,T1为在给定的初始速度vs、终止速度ve和目标位移L下,采用加速度极限值Amax和跃度极限值Jmax进行加速运动和匀速运动的最短持续时间,表示为:
Figure BDA0002296462920000115
其中,T1a为加速时间,表示为
Figure BDA0002296462920000116
定义临界时间T2,T2为仅进行加速运动时,从初始速度vs到终止速度ve的时间,表示为:
Figure BDA0002296462920000117
定义临界时间T3,T3为仅以初始速度vs进行匀速运动走完目标位移L的时间,表示为:
Figure BDA0002296462920000118
步骤3:根据步骤1确定的匀速速度上限值和步骤2计算的临界时间确定
Figure BDA00022964629200001210
速度曲线类型,具体方法为:
1)若目标时间Te小于临界时间T1,则
Figure BDA00022964629200001211
速度曲线包含加速阶段、匀速阶段和减速阶段,定义为类型A1,具体如下:
加速阶段的持续时间ta为:
Figure BDA0002296462920000121
减速阶段的持续时间td为:
Figure BDA0002296462920000122
匀速阶段的持续时间tc为Te-ta-td,匀速速度vm为:
Figure BDA0002296462920000123
Figure BDA0002296462920000124
速度曲线的指令位移为:
Figure BDA0002296462920000125
2)若T1<Te<T2
Figure BDA0002296462920000126
速度曲线包含加速阶段和匀速阶段,定义为类型A2,具体如下:
加速阶段的持续时间ta为:
Figure BDA0002296462920000127
匀速阶段的持续时间tc为Te-ta,终止速度ve为:
Figure BDA0002296462920000128
Figure BDA0002296462920000129
速度曲线的指令位移为:
Figure BDA0002296462920000131
3)若T2<Te<T3
Figure BDA0002296462920000132
速度曲线包含加速阶段,定义为类型A3,具体如下:
加速阶段的持续时间ta为Te,匀速阶段的终止速度ve为:
Figure BDA0002296462920000133
Figure BDA0002296462920000134
速度曲线的指令位移为:
Figure BDA0002296462920000135
4)若T3<Te
Figure BDA0002296462920000136
速度曲线先考虑先减速阶段后匀速阶段,定义为类型A4,具体如下:
减速阶段的持续时间td为:
Figure BDA0002296462920000137
匀速阶段的持续时间tc为Te-td,匀速阶段的终止速度ve为:
Figure BDA0002296462920000138
当ve≥0时,
Figure BDA0002296462920000139
速度曲线的指令位移为:
Figure BDA00022964629200001310
若ve<0时,
Figure BDA00022964629200001311
速度曲线需要进一步修改成减速阶段I、匀速阶段和减速阶段II,定义为类型A5,此时,采用Newtown-Raphson方法寻找保证ve≥0的匀速速度vm,并确定各阶段的持续时间,匀速速度vm的搜索范围为0和vs,终止条件为|Tev-Te|≤Ts
Figure BDA0002296462920000141
其中,td1和td2为减速阶段I和减速阶段II的持续时间,Tev为整个变化过程的时间,为实现变化时间等于目标时间Te,匀速速度vm为:
Figure BDA0002296462920000142
Figure BDA0002296462920000143
速度曲线的指令位移为:
Figure BDA0002296462920000144
基于
Figure BDA0002296462920000145
速度曲线的五轴轨迹加工方法,
Figure BDA0002296462920000146
速度曲线应用本发明所述的
Figure BDA0002296462920000147
速度曲线的设计方法;
预处理阶段:对五轴加工轨迹进行预处理,待插补的五轴加工轨迹包含刀尖子轨迹和刀轴方向子轨迹;首先,计算各子轨迹的曲率,根据曲率极值将各子轨迹分成离散的N个曲线段,曲率极值点为曲线段的衔接点,定义为临界点;然后,计算各曲线段的弧长;将临界点坐标、曲率极值和弧长存储在计算机中;
双向极值扫描阶段:对预处理后的五轴加工轨迹进行双向极值扫描,通过双向极值扫描更新临界点速度,双向极值扫描阶段包含向后看阶段和向前看阶段,扫描过程包含N个曲线段;向后看从轨迹终止点开始向起点反向扫描,向前看从轨迹起始点开始向终点正向扫描;
向后看的具体步骤如下:
1)令i=N,向后看获取的第i个临界点处的刀尖速度极值
Figure BDA0002296462920000156
刀轴方向速度极值
Figure BDA0002296462920000157
2)若i=1,执行3);
否则,令当前模块的刀尖初始速度
Figure BDA0002296462920000158
刀轴方向初始速度
Figure BDA0002296462920000159
执行2a);
2a)结合预处理阶段计算的刀尖子轨迹第i段曲线的轨迹长度
Figure BDA00022964629200001510
和切向跃度
Figure BDA00022964629200001511
将刀尖临界点速度记为vsj,由下式计算临界点速度vsj
Figure BDA0002296462920000151
结合刀尖子轨迹的切向加速度
Figure BDA00022964629200001513
将刀尖临界点速度记为vsa,vsa满足如下不等式:
Figure BDA0002296462920000152
结合预处理阶段计算的刀尖子轨迹第i段曲线的临界点曲率极值
Figure BDA00022964629200001512
预设的弓差误差d、预设的法向加速度
Figure BDA00022964629200001514
和预设的法向跃度
Figure BDA00022964629200001515
将刀尖临界点速度记为ven,其表达为:
Figure BDA0002296462920000153
其中,Ts为插补周期;
刀尖子轨迹的临界点速度
Figure BDA00022964629200001516
由切向约束和法向约束共同决定,表示为:
Figure BDA0002296462920000154
其中,vm为刀尖匀速速度;
2b)根据第i段刀尖子轨迹的初始速度
Figure BDA00022964629200001517
终止速度
Figure BDA00022964629200001518
轨迹长度
Figure BDA00022964629200001519
切向跃度
Figure BDA00022964629200001520
和切向加速度
Figure BDA00022964629200001521
结合
Figure BDA00022964629200001522
速度曲线的生成方法确定最短持续时间
Figure BDA00022964629200001523
2c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的轨迹长度
Figure BDA00022964629200001524
和切向跃度
Figure BDA00022964629200001525
将刀轴方向临界点速度记为ωsj,由下式计算临界点速度ωsj
Figure BDA0002296462920000155
结合刀轴方向子轨迹的切向加速度
Figure BDA00022964629200001526
将刀轴方向临界点速度记为ωsa,ωsa满足如下不等式:
Figure BDA0002296462920000161
结合预处理阶段计算的刀轴方向子轨迹第i段曲线的临界点曲率极值
Figure BDA0002296462920000164
法向加速度
Figure BDA0002296462920000165
和法向跃度
Figure BDA0002296462920000166
将刀轴方向临界点速度记为ωsn,其表示为:
Figure BDA0002296462920000162
刀轴方向临界点速度
Figure BDA0002296462920000167
由切向约束和法向约束共同决定,表示为:
Figure BDA0002296462920000163
其中,ωm为刀轴方向匀速速度;
2d)根据第i段刀轴方向子轨迹的初始速度
Figure BDA0002296462920000168
终止速度
Figure BDA0002296462920000169
轨迹长度
Figure BDA00022964629200001610
切向跃度
Figure BDA00022964629200001611
和切向加速度
Figure BDA00022964629200001612
结合
Figure BDA00022964629200001613
速度曲线的生成方法确定最短持续时间
Figure BDA00022964629200001614
2e)若
Figure BDA00022964629200001615
将刀轴方向子轨迹以
Figure BDA00022964629200001616
速度曲线修改并延长其持续时间至
Figure BDA00022964629200001617
先计算刀轴方向子轨迹的
Figure BDA00022964629200001618
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200001620
与临界时间比较并确定类型,记录修改后的临界速度为
Figure BDA00022964629200001619
Figure BDA00022964629200001621
将刀尖子轨迹以
Figure BDA00022964629200001622
速度曲线修改并延长其持续时间至
Figure BDA00022964629200001623
先计算刀尖子轨迹的
Figure BDA00022964629200001624
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200001625
与临界时间比较并确定类型,记录修改后的临界速度为
Figure BDA00022964629200001626
2f)令
Figure BDA00022964629200001627
i=i-1,执行2);
3)将
Figure BDA00022964629200001628
Figure BDA00022964629200001629
存储在数控单元缓存中;
向前看的具体步骤如下:
11)令i=1,向前看获取的第i个临界点处的刀尖速度极值vi=0,刀轴方向速度极值ωi=0;
22)若i=N,执行33);
否则,令当前模块的刀尖初始速度
Figure BDA00022964629200001630
刀轴方向初始速度
Figure BDA00022964629200001631
执行22a);
22a)结合预处理阶段计算的刀尖子轨迹第i段曲线的轨迹长度
Figure BDA00022964629200001632
和切向跃度
Figure BDA00022964629200001633
将刀尖临界点速度记为vej,由下式计算临界点速度vej
Figure BDA0002296462920000171
结合刀尖子轨迹的切向加速度
Figure BDA0002296462920000177
将刀尖临界点速度记为vea,vea满足如下不等式
Figure BDA0002296462920000172
刀尖子轨迹的临界点速度
Figure BDA0002296462920000178
由切向约束和向后看的临界速度共同决定,表示为:
Figure BDA0002296462920000173
22b)根据第i段刀尖子轨迹的初始速度
Figure BDA0002296462920000179
终止速度
Figure BDA00022964629200001710
轨迹长度
Figure BDA00022964629200001711
切向跃度
Figure BDA00022964629200001712
和切向加速度
Figure BDA00022964629200001713
结合
Figure BDA00022964629200001714
速度曲线的生成方法确定最短持续时间
Figure BDA00022964629200001715
22c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的轨迹长度
Figure BDA00022964629200001716
和切向跃度
Figure BDA00022964629200001717
将刀轴方向临界点速度记为ωej,由下式计算临界点速度ωsj
Figure BDA0002296462920000174
结合刀轴方向子轨迹的切向加速度
Figure BDA00022964629200001718
将刀轴方向临界点速度记为ωea,ωsa满足如下不等式:
Figure BDA0002296462920000175
刀轴方向子轨迹的临界点速度
Figure BDA00022964629200001719
由切向约束和向后看的临界速度共同决定,表示为
Figure BDA0002296462920000176
22d)根据第i段刀轴方向子轨迹的初始速度
Figure BDA00022964629200001720
终止速度
Figure BDA00022964629200001721
轨迹长度
Figure BDA00022964629200001722
切向跃度
Figure BDA00022964629200001723
和切向加速度
Figure BDA00022964629200001724
结合
Figure BDA00022964629200001725
速度曲线的生成方法确定最短持续时间
Figure BDA00022964629200001726
22e)若
Figure BDA00022964629200001727
将刀轴方向子轨迹以
Figure BDA00022964629200001728
速度曲线修改并延长其持续时间至
Figure BDA00022964629200001729
先计算刀轴方向子轨迹的
Figure BDA00022964629200001730
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200001731
与临界时间比较并确定类型,记录修改后的临界速度
Figure BDA00022964629200001732
和最大速度
Figure BDA00022964629200001733
Figure BDA00022964629200001734
将刀尖子轨迹以
Figure BDA00022964629200001735
速度曲线修改并延长其持续时间至
Figure BDA00022964629200001736
先计算刀尖子轨迹的
Figure BDA0002296462920000183
速度曲线类型的临界时间,然后根据目标时间
Figure BDA0002296462920000184
与临界时间比较并确定类型,记录修改后的临界速度
Figure BDA0002296462920000185
和最大速度
Figure BDA0002296462920000186
22f)令
Figure BDA0002296462920000187
i=i+1,执行22);
33)存储vi
Figure BDA0002296462920000188
和Ti(i=1,...,N)至数控单元缓存。
根据附图1所示的五轴加工轨迹说明本发明的具体实施方法,其中刀尖子轨迹匀速速度20mm/s,刀尖切向和法向加速度极限值为100mm/s2,刀尖切向和法向跃度极限值为1000mm/s3,弓差误差为2μm,刀轴方向子轨迹匀速速度为0.1rad/s,刀轴切向和法向加速度极限值为1rad/s2,刀轴切向和法向跃度极限值为10rad/s3,插补周期为1ms。
第一步,根据刀尖和刀轴方向子轨迹计算曲率并获得曲率极值,定义该点为临界点并将轨迹分成4段,计算各子曲线段的弧长;
第二步,双向极值扫描阶段:
1)令i=4,
Figure BDA0002296462920000189
2)若i=1,执行3)。否则,令
Figure BDA00022964629200001810
Figure BDA00022964629200001811
执行2a);
2a)结合预处理阶段计算的刀尖子轨迹第i段曲线的曲线长度
Figure BDA00022964629200001812
和切向跃度
Figure BDA00022964629200001813
由下式计算临界点速度vsj
Figure BDA0002296462920000181
其中:
Figure BDA00022964629200001814
为第i段曲线起点处的速度;
结合刀尖子轨迹的切向加速度
Figure BDA00022964629200001815
临界点速度vsa满足如下不等式
Figure BDA0002296462920000182
结合预处理阶段计算的刀尖子轨迹第i段曲线的临界点曲率极值
Figure BDA00022964629200001816
预设的弓差误差、预设的法向加速度、预设的法向跃度,刀尖临界点速度ven
Figure BDA0002296462920000191
其中:Ts为插补周期,
Figure BDA0002296462920000197
为刀尖方向子轨迹的法向加速度,
Figure BDA0002296462920000198
为刀尖子轨迹的法向跃度,d为刀尖子轨迹的弓差;
刀尖子轨迹的临界点速度由切向约束和法向约束共同决定,表示为
Figure BDA0002296462920000192
其中:vm为刀尖匀速速度;
2b)根据反向阶段的第i段刀尖子轨迹的初始速度
Figure BDA0002296462920000199
终止速度
Figure BDA00022964629200001910
轨迹长度
Figure BDA00022964629200001911
切向跃度
Figure BDA00022964629200001912
和切向加速度
Figure BDA00022964629200001913
结合
Figure BDA00022964629200001914
速度曲线的生成方法获得持续时间
Figure BDA00022964629200001915
2c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的曲线长度
Figure BDA00022964629200001916
和切向跃度
Figure BDA00022964629200001917
由下式计算临界点速度ωsj
Figure BDA0002296462920000193
其中:
Figure BDA00022964629200001918
为第i段曲线起点处的速度;
结合刀轴方向子轨迹的切向加速度
Figure BDA00022964629200001919
临界点速度ωsa满足如下不等式
Figure BDA0002296462920000194
结合预处理阶段计算的刀尖子轨迹第i段曲线的临界点曲率极值
Figure BDA00022964629200001920
法向加速度和法向跃度,临界点速度ωsn
Figure BDA0002296462920000195
其中,
Figure BDA00022964629200001922
为刀轴方向子轨迹的法向加速度,
Figure BDA00022964629200001921
为刀轴方向子轨迹的法向跃度;
刀轴方向子轨迹的临界点速度由切向约束和法向约束共同决定,表示为
Figure BDA0002296462920000196
其中:ωm为刀轴方向匀速速度;
2d)根据反向阶段的第i段刀轴方向子轨迹的初始速度
Figure BDA00022964629200002023
终止速度
Figure BDA00022964629200002024
轨迹长度
Figure BDA00022964629200002025
切向跃度
Figure BDA00022964629200002020
和切向加速度
Figure BDA00022964629200002019
结合
Figure BDA00022964629200002026
速度曲线的生成方法获得持续时间
Figure BDA00022964629200002027
2e)若
Figure BDA00022964629200002021
将刀轴方向子轨迹以
Figure BDA00022964629200002022
速度曲线修改并延长其持续时间至
Figure BDA00022964629200002028
先计算刀轴方向子轨迹的
Figure BDA00022964629200002029
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200002018
与临界时间比较并确定类型;
Figure BDA00022964629200002016
将刀尖子轨迹以
Figure BDA00022964629200002031
速度曲线修改并延长其持续时间至
Figure BDA00022964629200002030
先计算刀尖子轨迹的
Figure BDA00022964629200002017
速度曲线类型的临界时间,然后根据目标时间To i与临界时间比较并确定类型;
记录修改后的临界速度
Figure BDA00022964629200002014
Figure BDA00022964629200002015
2f)令
Figure BDA00022964629200002013
i=i-1,执行2);
3)将
Figure BDA00022964629200002011
Figure BDA00022964629200002012
存储在计算机中。
向前看从轨迹起始点开始向终点正向扫描,向前看具体执行流程如下:
11)令i=1,vi=0,ωi=0;
22)若i=4,转向33);
否则,令
Figure BDA0002296462920000209
Figure BDA00022964629200002010
转向22a);
22a)结合预处理阶段计算的刀尖子轨迹第i段曲线的曲线长度
Figure BDA00022964629200002032
和切向跃度
Figure BDA00022964629200002033
由下式计算临界点速度vej
Figure BDA0002296462920000201
其中:
Figure BDA00022964629200002034
为第i段曲线起点处的速度;
结合刀尖子轨迹的切向加速度
Figure BDA0002296462920000208
临界点速度vea满足如下不等式
Figure BDA0002296462920000202
刀尖子轨迹的临界点速度由切向约束和向后看的临界速度共同决定,表示为
Figure BDA0002296462920000203
22b)根据第i段刀尖子轨迹的起始速度
Figure BDA0002296462920000204
终止速度
Figure BDA0002296462920000205
轨迹长度
Figure BDA0002296462920000206
切向跃度
Figure BDA0002296462920000207
和切向加速度
Figure BDA0002296462920000214
结合
Figure BDA0002296462920000215
速度曲线的生成方法获得持续时间
Figure BDA0002296462920000216
22c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的曲线长度
Figure BDA0002296462920000217
和切向跃度
Figure BDA0002296462920000218
由下式计算临界点速度ωej
Figure BDA0002296462920000211
其中:
Figure BDA0002296462920000219
为第i段曲线起点处的速度;
结合刀轴方向子轨迹的切向加速度
Figure BDA00022964629200002110
临界点速度ωea满足如下不等式
Figure BDA0002296462920000212
刀轴方向子轨迹的临界点速度由切向约束和向后看的临界速度共同决定,表示为
Figure BDA0002296462920000213
22d)根据第i段刀轴方向子轨迹的初始速度
Figure BDA00022964629200002111
终止速度
Figure BDA00022964629200002112
轨迹长度
Figure BDA00022964629200002113
切向跃度
Figure BDA00022964629200002114
和切向加速度
Figure BDA00022964629200002115
结合
Figure BDA00022964629200002116
速度曲线的生成方法获得持续时间
Figure BDA00022964629200002117
22e)若
Figure BDA00022964629200002118
将刀轴方向子轨迹以
Figure BDA00022964629200002119
速度曲线修改并延长其持续时间至
Figure BDA00022964629200002120
先计算刀轴方向子轨迹的
Figure BDA00022964629200002121
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200002122
与临界时间比较并确定类型;
Figure BDA00022964629200002123
将刀尖子轨迹以
Figure BDA00022964629200002124
速度曲线修改并延长其持续时间至
Figure BDA00022964629200002125
先计算刀尖子轨迹的
Figure BDA00022964629200002134
速度曲线类型的临界时间,然后根据目标时间
Figure BDA00022964629200002126
与临界时间比较并确定类型;
记录修改后的临界速度
Figure BDA00022964629200002127
Figure BDA00022964629200002128
最大速度
Figure BDA00022964629200002129
Figure BDA00022964629200002130
22f)令
Figure BDA00022964629200002131
i=i+1,转向22);
33)存储vi,ωi
Figure BDA00022964629200002132
和Ti(i=1,...,4)至计算机缓存。
刀尖和刀轴子轨迹在引入时间同步算法前后的持续时间如下表1所示,附图4说明刀尖和刀轴子轨迹在引入时间同步算法前后的速度曲线,其中曲线1为引入时间同步算法前的速度曲线,曲线2为引入时间同步算法后的速度曲线。在第1和第3个曲线段,刀尖子轨迹的运动时间大于刀轴方向子轨迹时间,引入
Figure BDA00022964629200002133
速度曲线的类型A2和类型A4来延长刀轴方向子轨迹的运动时间,曲线如附图4(a)所示。在第2个和第4个曲线段,刀轴方向子轨迹的运动时间大于刀尖方向子轨时间引入
Figure BDA0002296462920000222
速度曲线的类型A3和类型A4延长刀尖字轨迹的运动时间,曲线如附图4(b)所示。从表1可知,引入时间同步算法后,两子轨迹在各阶段的时间完全相同,说明时间同步算法能够有效的实现两子轨迹的同步运动。
表1时间同步算法前后刀尖和刀轴方向子轨迹的运动时间
Figure BDA0002296462920000221
附图5为刀轴方向速度曲线图,曲线1为参数线性同步算法对应的刀轴方向速度曲线图,曲线2为本发明提出的时间同步算法对应的刀轴方向速度曲线图。图中可见,参数线性同步算法忽略了刀轴方向的运动学约束,导致刀轴方向速度超过了预设的速度极限值,而本发明提出的时间同步算法有效地遵守了刀轴方向的运动学约束,刀轴方向的速度曲线没有超过预设的极限值。

Claims (5)

1.一种
Figure FDA0002296462910000011
速度曲线的设计方法,其特征在于,包括以下步骤:
步骤1:确定
Figure FDA0002296462910000012
速度曲线的匀速速度上限值;
步骤2:计算临界时间,用于在步骤3中确定
Figure FDA0002296462910000013
速度曲线的类型;
步骤3:根据步骤1确定的匀速速度上限值和步骤2计算的临界时间确定
Figure FDA0002296462910000014
速度曲线类型。
2.根据权利要求1所述的一种
Figure FDA0002296462910000015
速度曲线的设计方法,其特征在于,所述步骤1的具体方法为:
给定参数,所述参数包括目标位移L、初始速度vs、终止速度ve、匀速速度vm、加速度极限值Amax和跃度极限值Jmax;根据给定的所述参数假设运动过程包括加速阶段、匀速阶段和减速阶段,其中,加速阶段的持续时间ta为:
Figure FDA0002296462910000016
减速阶段的持续时间td为:
Figure FDA0002296462910000017
Figure FDA0002296462910000018
时,假设成立,运动过程存在匀速阶段,且匀速速度为vm
Figure FDA0002296462910000019
时,假设不成立,运动过程包括加速阶段和减速阶段,且速度曲线的最大速度无法取匀速速度vm,实际最大速度vma借助Newtown-Raphson方法在max(vs,ve)和vm之间寻找,直至匀速阶段的持续时间tc介于0与插补周期Ts之间,迭代过程结束,匀速阶段的持续时间tc为0,匀速速度vm为实际最大速度vma
Figure FDA00022964629100000110
本步骤获取的匀速速度vm
Figure FDA00022964629100000111
速度曲线的匀速速度上限值vmm
3.根据权利要求2所述的一种
Figure FDA0002296462910000021
速度曲线的设计方法,其特征在于,所述步骤2的具体方法为:
定义临界时间T1,T1为在给定的初始速度vs、终止速度ve和目标位移L下,采用加速度极限值Amax和跃度极限值Jmax进行加速运动和匀速运动的最短持续时间,表示为:
Figure FDA0002296462910000022
其中,T1a为加速时间,表示为
Figure FDA0002296462910000023
定义临界时间T2,T2为仅进行加速运动时,从初始速度vs到终止速度ve的时间,表示为:
Figure FDA0002296462910000024
定义临界时间T3,T3为仅以初始速度vs进行匀速运动走完目标位移L的时间,表示为:
Figure FDA0002296462910000025
4.根据权利要求3所述的一种
Figure FDA0002296462910000026
速度曲线的设计方法,其特征在于,所述步骤3的具体方法为:
1)若目标时间Te小于临界时间T1,则
Figure FDA0002296462910000027
速度曲线包含加速阶段、匀速阶段和减速阶段,定义为类型A1,具体如下:
加速阶段的持续时间ta为:
Figure FDA0002296462910000028
减速阶段的持续时间td为:
Figure FDA0002296462910000029
匀速阶段的持续时间tc为Te-ta-td,匀速速度vm为:
Figure FDA00022964629100000210
Figure FDA0002296462910000031
速度曲线的指令位移为:
Figure FDA0002296462910000032
2)若T1<Te<T2
Figure FDA0002296462910000033
速度曲线包含加速阶段和匀速阶段,定义为类型A2,具体如下:加速阶段的持续时间ta为:
Figure FDA0002296462910000034
匀速阶段的持续时间tc为Te-ta,终止速度ve为:
Figure FDA0002296462910000035
Figure FDA0002296462910000036
速度曲线的指令位移为:
Figure FDA0002296462910000037
3)若T2<Te<T3
Figure FDA0002296462910000038
速度曲线包含加速阶段,定义为类型A3,具体如下:
加速阶段的持续时间ta为Te,匀速阶段的终止速度ve为:
Figure FDA0002296462910000039
Figure FDA00022964629100000310
速度曲线的指令位移为:
Figure FDA00022964629100000311
4)若T3<Te
Figure FDA00022964629100000312
速度曲线先考虑先减速阶段后匀速阶段,定义为类型A4,具体如下:减速阶段的持续时间td为:
Figure FDA00022964629100000313
匀速阶段的持续时间tc为Te-td,匀速阶段的终止速度ve为:
Figure FDA0002296462910000041
当ve≥0时,
Figure FDA0002296462910000042
速度曲线的指令位移为:
Figure FDA0002296462910000043
若ve<0时,
Figure FDA0002296462910000044
速度曲线需要进一步修改成减速阶段I、匀速阶段和减速阶段II,定义为类型A5,此时,采用Newtown-Raphson方法寻找保证ve≥0的匀速速度vm,并确定各阶段的持续时间,匀速速度vm的搜索范围为0和vs,终止条件为|Tev-Te|≤Ts
Figure FDA0002296462910000045
其中,td1和td2为减速阶段I和减速阶段II的持续时间,Tev为整个变化过程的时间,为实现变化时间等于目标时间Te,匀速速度vm为:
Figure FDA0002296462910000046
Figure FDA0002296462910000047
速度曲线的指令位移为:
Figure FDA0002296462910000048
5.基于
Figure FDA0002296462910000049
速度曲线的五轴轨迹加工方法,其特征在于:
Figure FDA00022964629100000410
速度曲线应用权利要求1至4任一项所述的设计方法;
预处理阶段:对五轴加工轨迹进行预处理,待插补的五轴加工轨迹包含刀尖子轨迹和刀轴方向子轨迹;首先,计算各子轨迹的曲率,根据曲率极值将各子轨迹分成离散的N个曲线段,曲率极值点为曲线段的衔接点,定义为临界点;然后,计算各曲线段的弧长;将临界点坐标、曲率极值和弧长存储在计算机中;
双向极值扫描阶段:对预处理后的五轴加工轨迹进行双向极值扫描,通过双向极值扫描更新临界点速度,双向极值扫描阶段包含向后看阶段和向前看阶段,扫描过程包含N个曲线段;向后看从轨迹终止点开始向起点反向扫描,向前看从轨迹起始点开始向终点正向扫描;
向后看的具体步骤如下:
1)令i=N,向后看获取的第i个临界点处的刀尖速度极值
Figure FDA0002296462910000051
刀轴方向速度极值
Figure FDA0002296462910000052
2)若i=1,执行3);
否则,令当前模块的刀尖初始速度
Figure FDA0002296462910000053
刀轴方向初始速度
Figure FDA0002296462910000054
执行2a);
2a)结合预处理阶段计算的刀尖子轨迹第i段曲线的轨迹长度
Figure FDA0002296462910000055
和切向跃度
Figure FDA0002296462910000056
将刀尖临界点速度记为vsj,由下式计算临界点速度vsj
Figure FDA0002296462910000057
结合刀尖子轨迹的切向加速度
Figure FDA0002296462910000058
将刀尖临界点速度记为vsa,vsa满足如下不等式:
Figure FDA0002296462910000059
结合预处理阶段计算的刀尖子轨迹第i段曲线的临界点曲率极值
Figure FDA00022964629100000510
预设的弓差误差δ、预设的法向加速度
Figure FDA00022964629100000511
和预设的法向跃度
Figure FDA00022964629100000512
将刀尖临界点速度记为ven,其表达为:
Figure FDA00022964629100000513
其中,Ts为插补周期;
刀尖子轨迹的临界点速度
Figure FDA00022964629100000514
由切向约束和法向约束共同决定,表示为:
Figure FDA00022964629100000515
其中,vm为刀尖匀速速度;
2b)根据第i段刀尖子轨迹的初始速度
Figure FDA0002296462910000061
终止速度
Figure FDA0002296462910000062
轨迹长度
Figure FDA0002296462910000063
切向跃度
Figure FDA0002296462910000064
和切向加速度
Figure FDA0002296462910000065
结合
Figure FDA0002296462910000066
速度曲线的生成方法确定最短持续时间
Figure FDA0002296462910000067
2c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的轨迹长度
Figure FDA0002296462910000068
和切向跃度
Figure FDA0002296462910000069
将刀轴方向临界点速度记为ωsj,由下式计算临界点速度ωsj
Figure FDA00022964629100000610
结合刀轴方向子轨迹的切向加速度
Figure FDA00022964629100000611
将刀轴方向临界点速度记为ωsa,ωsa满足如下不等式:
Figure FDA00022964629100000612
结合预处理阶段计算的刀轴方向子轨迹第i段曲线的临界点曲率极值
Figure FDA00022964629100000613
法向加速度
Figure FDA00022964629100000614
和法向跃度
Figure FDA00022964629100000615
将刀轴方向临界点速度记为ωsn,其表示为:
Figure FDA00022964629100000616
刀轴方向临界点速度
Figure FDA00022964629100000617
由切向约束和法向约束共同决定,表示为:
Figure FDA00022964629100000618
其中,ωm为刀轴方向匀速速度;
2d)根据第i段刀轴方向子轨迹的初始速度
Figure FDA00022964629100000619
终止速度
Figure FDA00022964629100000620
轨迹长度
Figure FDA00022964629100000621
切向跃度
Figure FDA00022964629100000622
和切向加速度
Figure FDA00022964629100000623
结合
Figure FDA00022964629100000624
速度曲线的生成方法确定最短持续时间
Figure FDA00022964629100000625
2e)若
Figure FDA00022964629100000626
将刀轴方向子轨迹以
Figure FDA00022964629100000627
速度曲线修改并延长其持续时间至
Figure FDA00022964629100000628
先计算刀轴方向子轨迹的
Figure FDA00022964629100000629
速度曲线类型的临界时间,然后根据目标时间
Figure FDA00022964629100000630
与临界时间比较并确定类型,记录修改后的临界速度为
Figure FDA00022964629100000631
Figure FDA00022964629100000632
将刀尖子轨迹以
Figure FDA00022964629100000633
速度曲线修改并延长其持续时间至
Figure FDA00022964629100000634
先计算刀尖子轨迹的
Figure FDA00022964629100000635
速度曲线类型的临界时间,然后根据目标时间
Figure FDA00022964629100000636
与临界时间比较并确定类型,记录修改后的临界速度为
Figure FDA0002296462910000071
2f)令
Figure FDA0002296462910000072
i=i-1,执行2);
3)将
Figure FDA0002296462910000073
Figure FDA0002296462910000074
存储在数控单元缓存中;
向前看的具体步骤如下:
11)令i=1,向前看获取的第i个临界点处的刀尖速度极值vi=0,刀轴方向速度极值ωi=0;
22)若i=N,执行33);
否则,令当前模块的刀尖初始速度
Figure FDA0002296462910000075
刀轴方向初始速度
Figure FDA0002296462910000076
执行22a);
22a)结合预处理阶段计算的刀尖子轨迹第i段曲线的轨迹长度
Figure FDA0002296462910000077
和切向跃度
Figure FDA0002296462910000078
将刀尖临界点速度记为vej,由下式计算临界点速度vej
Figure FDA0002296462910000079
结合刀尖子轨迹的切向加速度
Figure FDA00022964629100000710
将刀尖临界点速度记为vea,vea满足如下不等式
Figure FDA00022964629100000711
刀尖子轨迹的临界点速度
Figure FDA00022964629100000712
由切向约束和向后看的临界速度共同决定,表示为:
Figure FDA00022964629100000713
22b)根据第i段刀尖子轨迹的初始速度
Figure FDA00022964629100000714
终止速度
Figure FDA00022964629100000715
轨迹长度
Figure FDA00022964629100000716
切向跃度
Figure FDA00022964629100000717
和切向加速度
Figure FDA00022964629100000718
结合
Figure FDA00022964629100000719
速度曲线的生成方法确定最短持续时间
Figure FDA00022964629100000720
22c)结合预处理阶段计算的刀轴方向子轨迹第i段曲线的轨迹长度
Figure FDA00022964629100000721
和切向跃度
Figure FDA00022964629100000722
将刀轴方向临界点速度记为ωej,由下式计算临界点速度ωsj
Figure FDA00022964629100000723
结合刀轴方向子轨迹的切向加速度
Figure FDA00022964629100000724
将刀轴方向临界点速度记为ωea,ωsa满足如下不等式:
Figure FDA00022964629100000725
刀轴方向子轨迹的临界点速度
Figure FDA0002296462910000081
由切向约束和向后看的临界速度共同决定,表示为
Figure FDA0002296462910000082
22d)根据第i段刀轴方向子轨迹的初始速度
Figure FDA0002296462910000083
终止速度
Figure FDA0002296462910000084
轨迹长度
Figure FDA0002296462910000085
切向跃度
Figure FDA0002296462910000086
和切向加速度
Figure FDA0002296462910000087
结合
Figure FDA0002296462910000088
速度曲线的生成方法确定最短持续时间
Figure FDA0002296462910000089
22e)若
Figure FDA00022964629100000810
将刀轴方向子轨迹以
Figure FDA00022964629100000811
速度曲线修改并延长其持续时间至
Figure FDA00022964629100000812
先计算刀轴方向子轨迹的
Figure FDA00022964629100000813
速度曲线类型的临界时间,然后根据目标时间
Figure FDA00022964629100000814
与临界时间比较并确定类型,记录修改后的临界速度
Figure FDA00022964629100000815
和最大速度v′m
Figure FDA00022964629100000816
将刀尖子轨迹以
Figure FDA00022964629100000817
速度曲线修改并延长其持续时间至
Figure FDA00022964629100000818
先计算刀尖子轨迹的
Figure FDA00022964629100000819
速度曲线类型的临界时间,然后根据目标时间
Figure FDA00022964629100000820
与临界时间比较并确定类型,记录修改后的临界速度
Figure FDA00022964629100000821
和最大速度ω′m
22f)令
Figure FDA00022964629100000822
i=i+1,执行22);
33)存储vi,ωi
Figure FDA00022964629100000823
和Ti(i=1,...,N)至数控单元缓存。
CN201911203626.1A 2019-11-29 2019-11-29 St2速度曲线的设计方法及基于sst2速度曲线的五轴轨迹加工方法 Active CN112883502B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911203626.1A CN112883502B (zh) 2019-11-29 2019-11-29 St2速度曲线的设计方法及基于sst2速度曲线的五轴轨迹加工方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911203626.1A CN112883502B (zh) 2019-11-29 2019-11-29 St2速度曲线的设计方法及基于sst2速度曲线的五轴轨迹加工方法

Publications (2)

Publication Number Publication Date
CN112883502A true CN112883502A (zh) 2021-06-01
CN112883502B CN112883502B (zh) 2022-12-09

Family

ID=76038868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911203626.1A Active CN112883502B (zh) 2019-11-29 2019-11-29 St2速度曲线的设计方法及基于sst2速度曲线的五轴轨迹加工方法

Country Status (1)

Country Link
CN (1) CN112883502B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114035582A (zh) * 2021-11-16 2022-02-11 深圳数马电子技术有限公司 S形轨迹规划方法、装置、计算机设备和存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0704781A2 (de) * 1994-09-29 1996-04-03 Siemens Aktiengesellschaft Geschwindigkeitsabhängige Beschleunigungskennlinie
US20030155881A1 (en) * 2002-01-10 2003-08-21 Siemens Aktiengesellschaft Method for adapting jerk profiles through jerk limitation for minimizing excitations of oscillations during movement of a machine element
CN107608313A (zh) * 2017-09-11 2018-01-19 大连理工大学 一种五轴双样条曲线插补速度规划方法
CN107765639A (zh) * 2017-09-30 2018-03-06 山东大学 一种s曲线加减速的圆整误差补偿方法
CN110134065A (zh) * 2019-05-30 2019-08-16 西安交通大学 一种基于正弦平方加速度前瞻的机床加工轨迹运动规划方法
CN110209173A (zh) * 2019-06-28 2019-09-06 中科新松有限公司 基于s型速度规划的时间最优在线轨迹规划方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0704781A2 (de) * 1994-09-29 1996-04-03 Siemens Aktiengesellschaft Geschwindigkeitsabhängige Beschleunigungskennlinie
US20030155881A1 (en) * 2002-01-10 2003-08-21 Siemens Aktiengesellschaft Method for adapting jerk profiles through jerk limitation for minimizing excitations of oscillations during movement of a machine element
CN107608313A (zh) * 2017-09-11 2018-01-19 大连理工大学 一种五轴双样条曲线插补速度规划方法
CN107765639A (zh) * 2017-09-30 2018-03-06 山东大学 一种s曲线加减速的圆整误差补偿方法
CN110134065A (zh) * 2019-05-30 2019-08-16 西安交通大学 一种基于正弦平方加速度前瞻的机床加工轨迹运动规划方法
CN110209173A (zh) * 2019-06-28 2019-09-06 中科新松有限公司 基于s型速度规划的时间最优在线轨迹规划方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
XAVIER BEUDAERT等: "5-axis local corner rounding of linear tool path discontinuities", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS AND MANUFACTURE》 *
XIAOYONG等: "A novel triple-stage friction compensation for a feed system based on electromechanical characteristics", 《PRECISION ENGINEERING》 *
孙竞潇等: "临界转速测量与实时显示", 《微电子技术》 *
张强等: "切向加速度与弦误差约束的CNC系统时间最优轨迹规划", 《系统工程理论与实践》 *
赵欢等: "五轴线性刀路的转接光顺及轨迹生成算法", 《机械工程学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114035582A (zh) * 2021-11-16 2022-02-11 深圳数马电子技术有限公司 S形轨迹规划方法、装置、计算机设备和存储介质

Also Published As

Publication number Publication date
CN112883502B (zh) 2022-12-09

Similar Documents

Publication Publication Date Title
Tajima et al. Kinematic corner smoothing for high speed machine tools
Tajima et al. Global tool-path smoothing for CNC machine tools with uninterrupted acceleration
CN109571473B (zh) 一种误差可控的小线段轨迹光顺方法
CN108829045B (zh) 连续微直线段的衔接速度的优化方法及系统
CN109623166B (zh) 一种激光切割的拐角处理方法及系统
JP5840288B2 (ja) モータ制御装置
WO2019047458A1 (zh) 一种五轴双样条曲线插补速度规划方法
CN107765648B (zh) 一种cnc加工的进给速度规划方法及装置
CN108227630B (zh) 一种采用时间参数多项式插补的自由曲面数控加工方法
CN112486101B (zh) Nurbs曲线自适应前瞻插补方法
CN113156893B (zh) 一种基于s形加减速的五轴机床速度规划方法
CN109048091B (zh) 激光切割速度规划方法、装置、存储介质及计算机设备
CN111897290A (zh) 一种轴加加速度光滑的拐角过渡平滑方法
JP2010511919A (ja) 許容差ベースの経路設計と制御の方法
CN108303952B (zh) 一种适用于数控装置刀具的小线段的速度曲线规划方法
US20110062118A1 (en) Method and machine for machining parts using spark-erosion milling
CN105785909B (zh) 一种高速高精的路径动态前瞻规划方法
CN102707671A (zh) 应用于工具机的加工路径最佳化方法
CN109901518B (zh) 一种恒力约束条件下的数控机床加减速速度规划方法
CN111633668A (zh) 一种用于机器人加工三维自由曲面的运动控制方法
CN116300698A (zh) 一种基于动态规划的螺纹切削加工方法
CN112883502B (zh) St2速度曲线的设计方法及基于sst2速度曲线的五轴轨迹加工方法
TWI336822B (en) Method of tolerance-based trajectory planning and control
CN114019910A (zh) 一种小线段刀具轨迹实时全局光顺方法
CN110531700B (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