CN109313428B - 运动轨迹平滑转接的方法和装置及其相关设备 - Google Patents

运动轨迹平滑转接的方法和装置及其相关设备 Download PDF

Info

Publication number
CN109313428B
CN109313428B CN201780034014.7A CN201780034014A CN109313428B CN 109313428 B CN109313428 B CN 109313428B CN 201780034014 A CN201780034014 A CN 201780034014A CN 109313428 B CN109313428 B CN 109313428B
Authority
CN
China
Prior art keywords
track
arc
curve
point
circular arc
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
CN201780034014.7A
Other languages
English (en)
Other versions
CN109313428A (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.)
Shenzhen A&E Intelligent Technology Institute Co Ltd
Original Assignee
Shenzhen A&E Intelligent Technology Institute Co Ltd
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 Shenzhen A&E Intelligent Technology Institute Co Ltd filed Critical Shenzhen A&E Intelligent Technology Institute Co Ltd
Publication of CN109313428A publication Critical patent/CN109313428A/zh
Application granted granted Critical
Publication of CN109313428B publication Critical patent/CN109313428B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/02Programme-controlled manipulators characterised by movement of the arms, e.g. cartesian coordinate type
    • 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/41Numerical 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 interpolation, e.g. the computation of intermediate points between programmed end points to define the path to be followed and the rate of travel along that path
    • G05B19/4103Digital interpolation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Numerical Control (AREA)
  • Manipulator (AREA)

Abstract

本发明公开了一种运动轨迹平滑转接的方法和装置及其相关设备,该方法包括:获取待平滑转接的第一和第二轨迹以及预设轮廓误差;根据第一和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一和第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线。通过上述方式,本发明能够准确的拟合出实现曲率和速度均连续的插补曲线且能够保证设备的轮廓误差。

Description

运动轨迹平滑转接的方法和装置及其相关设备
技术领域
本发明涉及运动控制技术领域,特别是涉及一种运动轨迹平滑转接的方法和装置及其相关设备。
背景技术
在很多情况下,相邻两段运动轨迹间通常会存在一个夹角,运动到该处时,进给速度的方向发生突变,若保持速度大小不变,则会产生一个较大的加速度,这个加速度可能超出设备的加速能力,导致对应的设备震动。为了让此加速度小于***最大加速度,需要减小两段运动轨迹连接点处的速度,大量的运动轨迹累积下来,使得运动效率受到严重影响。
为此,有人提出在存在夹角的两段轨迹之间引入一段插补曲线,用“圆角”近似“尖角”,从而避免在尖角处速度需要降低很多的做法。不同的插补曲线都有一个最基本的统一要求,即插补曲线与上述两段轨迹的连接点位置保证C1连续(一阶导数连续),从而速度连续。进一步,我们希望在连接点处的曲率也连续,即C2连续(二阶导数连续),这样可以保证加速度不会出现跳变引起设备震动,在保证加工效率的同时保证了加工的平稳性,可以进一步提升***性能。
现有的插补曲线通常为近似的圆弧,圆弧计算方式虽然简单但是无法保证在圆弧与两段轨迹的连接点处的C2连续,只能保证C1连续。
在利用其它曲线对上述两段轨迹进行插补时,由于计算的条件不充分,无法准确的确定曲线的方程,充其量也是只能近似的算出插补曲线。
另外,设备通常会设定一个轮廓误差,目前在计算插补曲线时通常无法保证轮廓误差,导致实际按照计算得到的插补曲线加工时,仍然会由于不平滑出现运动不平稳运动效率低的问题。
因此,需要提供运动轨迹平滑转接的方法和装置及其相关设备,以解决上述技术问题。
发明内容
本发明主要解决的技术问题是提供一种运动轨迹平滑转接的方法和装置及其相关设备,能够准确的拟合出实现曲率和速度均连续的插补曲线且能够保证设备的轮廓误差。
为解决上述技术问题,本发明采用的一个技术方案是:提供一种运动轨迹平滑转接的方法,该方法包括:获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等;利用插补曲线平滑转接第一轨迹和第二轨迹。
为解决上述技术问题,本发明采用的一个技术方案是:一种运动轨迹平滑转接的装置,装置包括处理器和与处理器连接的存储器,存储器用于存储第一轨迹、第二轨迹、预设轮廓误差,处理器用于从存储器获取第一轨迹、第二轨迹以及预设轮廓误差;根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等。
为解决上述技术问题,本发明采用的又一个技术方案是:提供一种存储装置,该存储装置存储有程序数据,程序数据能够被执行以实现上述的方法。
为解决上述技术问题,本发明采用的又一个技术方案是:提供一种机器人***,该机器人***包括机器人以及与机器人连接的控制器,控制器用于获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等;以及根据第一轨迹、第二轨迹以及插补曲线控制机器人或者机器人的部件的运动。
为解决上述技术问题,本发明采用的又一个技术方案是:提供一种数控机床,该数控机床包括机床本体和数控设备,数控设备用于获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等;根据第一轨迹、第二轨迹以及插补曲线控制机床本体或者机床本体的部件的运动。
本发明的有益效果是:区别于现有技术的情况,本发明通过将轮廓误差作为拟合插补曲线的条件,根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置,再根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,一方面可以准确的拟合出插补曲线,另一方面拟合出的插补曲线可以保证设备运动时的轮廓误差。
附图说明
图1是本发明第一实施例运动轨迹平滑转接的方法流程图;
图2是本发明第二实施例运动轨迹平滑转接的方法流程图;
图3是线段、圆弧以及插补曲线的几何示意图;
图4是本发明第三实施例运动轨迹平滑转接的方法流程图;
图5是第一圆弧、第二圆弧以及插补曲线的几何示意图;
图6是本发明实施例的运动轨迹平滑转接装置的模块示意图;
图7是本发明实施例的存储装置的示意图;
图8是本发明实施例的机器人***的模块示意图;
图9是本发明实施例的数控机床的模块示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的说明。
请参阅图1,图1是本发明第一实施例运动轨迹平滑转接的方法流程图。在本实施例中,运动轨迹平滑转接的方法包括以下步骤:
步骤S11:获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差。
在本实施例中,所述第一轨迹为圆弧或线段,第二轨迹为圆弧,具体请参见下文的描述。
步骤S12:根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置。
在本实施例中,所述根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置可以为:在第一轨迹和第二轨迹二者的连接点处的切线夹角的角平分线上选择距离第一轨迹和第二轨迹的连接点为预设轮廓误差的点作为特征插补点,具体请参见下文的描述。
步骤S13:根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等,具体请参见下文的描述。
在本实施例中,所述根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线的步骤包括:拟合出连接第一轨迹和特征插补点的第一子曲线以及连接特征插补点和第二轨迹之间的第二子曲线,第一子曲线和第二子曲线分别满足插补曲线通用方程且具有不同的曲线参数,第一子曲线和第二子曲线在特征插补点处的曲率与切向角均相等,具体请参见下文的描述。
在本实施例中,采用分为两段子曲线的方式进行拟合。在其他实施例中,也可以采用分为三段或三段以上子曲线的方式进行拟合,也可以采用一段曲线的方式进行拟合。
在本实施例中,所述第一子曲线为第一回旋(clothoid)曲线,所述第二子曲线为第二回旋曲线。在其他实施例中,所述第一子曲线和第二子曲线还可以采用其他的特性曲线。
请参阅图2,图2是本发明第二实施例运动轨迹平滑转接的方法流程图。
在本实施例中,优选地,根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线的步骤包括:拟合出连接第一轨迹和特征插补点的第一子曲线以及连接特征插补点和第二轨迹之间的第二子曲线,第一子曲线和第二子曲线分别满足插补曲线通用方程且具有不同的曲线参数,第一子曲线和第二子曲线在特征插补点处的曲率与切向角均相等。第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等。第一轨迹为线段,第二轨迹为圆弧。第一子曲线为第一回旋(clothoid)曲线,第二子曲线为第二回旋曲线。
在本实施例中,运动轨迹平滑转接的方法包括以下步骤:
步骤S21:获取待平滑转接的线段和圆弧以及预设轮廓误差。
在本实施例中,所述获取待平滑转接的线段和圆弧包括:获取线段的位置和圆弧的位置以及二者之间的位置关系。
步骤S22:根据线段和圆弧以及预设轮廓误差确定特征插补点的位置。
在本实施例中,所述特征插补点的确定条件为:圆弧和线段连接点处的切线与线段的夹角的角平分线上选择距离线段和圆弧的连接点为预设轮廓误差的点作为特征插补点。也即图3或者图4中的Q点,具体参见下文的描述。
步骤S23:将线段设置在直角坐标系的一坐标轴上,且线段的未与圆弧连接的一端设置成与直角坐标系的原点重合;
也即,以线段的未与圆弧连接的一端为原点,以线段为一个坐标轴建立一个坐标系,具体参见下文的描述。
步骤S24:通过联立以下两组方程求解出第一回旋曲线和第二回旋曲线的曲线参数:
Figure BDA0001887568300000041
Figure BDA0001887568300000042
其中,x0为第一回旋曲线与线段的连接点在坐标轴上的坐标,c0为第一回旋曲线的曲率变化率,s0为第一回旋曲线从其与线段的连接点到特征插补点之间的弧长,c1为第二回旋曲线的曲率变化率,s1为第二回旋曲线从特征插补点到其与圆弧的连接点之间的弧长,α为圆弧从第二回旋曲线与圆弧的连接点到圆弧与线段的连接点之间的圆弧片段的圆心角,Tccw为旋转矩阵,Tccw、x0、c0、s0、c1、s1以及α为需要通过两组方程进行求解的未知量;
其中,L为线段的长度,E为预设轮廓误差,β为线段与圆弧的夹角,R为圆弧的半径,xc为圆弧的圆心在坐标轴上的坐标,yc为圆弧的圆心在直角坐标系另一坐标轴上的坐标,
Figure BDA0001887568300000051
为圆弧的圆心到圆弧与线段的连接点的方向向量,xc、yc
Figure BDA0001887568300000052
L、β和R为能够通过线段和圆弧及二者的相对位置关系获得的已知量。
下面说明推导过程。
首先介绍插补曲线的通用方程,插补曲线为回旋(clothoid)曲线,其具有曲率随路径弧长线性变化的特点。
插补曲线的通用方程表达形式有三种:
k(s)=k0+cs(3)
Figure BDA0001887568300000053
Figure BDA0001887568300000054
其中,s为曲线路径弧长,x(s)为曲线路径弧长为s处对应的横坐标,x0为曲线起点处的横坐标,y0为曲线起点处的横坐标,y(s)为曲线路径弧长为s处对应的纵坐标,θ(s)为曲线路径弧长为s处对应的切向角,θ0为曲线起点的切向角,k0为曲线起点的曲率,c为曲线的曲率变化率,k(s)为曲线路径弧长为s处对应的曲率。
其次,进一步说明特征插补点的位置的确定以及方程组(1)和方程组(2)的推导。记第一轨迹与插补曲线在二者连接点的曲率与切向角相等为第一平滑转接条件,第二轨迹与插补曲线在二者连接点的曲率与切向角相等为第二平滑转接条件,第一子曲线和第二子曲线在特征插补点处的曲率与切向角均相等为第三平滑转接条件。
请参阅图3,图3是线段、圆弧以及插补曲线的几何示意图。
下面以线段NP和圆弧PM的平滑转接为例进行说明。I点为线段NP与第一回旋曲线IQ的连接点,O点为圆弧PM与第二回旋曲线QO的连接点,P点为线段NP与圆弧PM的连接点,C为圆弧PM的圆心。关于运动方向,以N点→I点→Q点→O点→M点为例。
设P点和O点之间的圆弧片段OP对应的圆心角为α。
在步骤S21中,获取待平滑转接的线段和圆弧以及预设轮廓误差包括:
获取预设轮廓误差E、线段NP的长度L、线段NP与圆弧PM的连接点P、线段NP与圆弧PM的夹角β、圆心C的位置以及圆弧的圆心到圆弧与线段的连接点的方向向量。
在步骤S22中,根据前文描述的特征插补点的确定条件为特征插补点Q在线段NP与圆弧PM的切线的夹角的角平分线上且离P点距离QP为E以及根据预设轮廓误差E、夹角β、线段NP的长度L确定特征插补点Q的位置。
在步骤S23和步骤S24中:
1、先建立直角坐标系,优选地,将线段NP设置在x轴上,将N点设置为与坐标原点重合,结合步骤S22可以得出特征插补点Q在建立的坐标系中的坐标为:
Figure BDA0001887568300000061
在其他实施例中,x轴与y轴的位置可以调换,例如将线段NP设置在y轴上,将N点设置为与坐标原点重合。
2、然后根据第一平滑转接条件和第一轨迹确定第一轨迹与第一回旋曲线的连接点的部分参数和第一回旋曲线在第一轨迹与第一回旋曲线的连接点的部分参数:
由线段NP在x轴上,可以确定I点的纵坐标为0;
由线段NP的曲率和切向角处处为0,且由第一平滑转接条件可得第一回旋曲线在点I处曲率为和切向角均为0。
3、根据第一回旋曲线满足插补曲线通用方程将第一轨迹与第一回旋曲线连接点的部分参数、第一回旋曲线在第一轨迹与第一回旋曲线连接点的部分参数以及特征插补点的位置带入插补曲线通用方程得到方程组(1):
将I点纵坐标为0,第一回旋曲线在I点处曲率为和切向角均为0以及插补特征点的坐标带入插补曲线通用方程得到方程组(1)
Figure BDA0001887568300000062
4、根据第二平滑转接条件和第二轨迹确定第二轨迹与第二回旋曲线的连接点的部分参数和第二回旋曲线在第二轨迹与第二回旋曲线的连接点的部分参数:
根据建立的直角坐标系和获取的圆心的位置、圆弧的圆心到圆弧与线段的连接点的方向向量得到圆心C在该直角坐标系的坐标为(xc,yc),C点到P点的方向向量
Figure BDA0001887568300000063
将点O的坐标表示为
Figure BDA0001887568300000064
5、根据第二回旋曲线满足插补曲线通用方程将特征插补点的位置、第二轨迹与第二回旋曲线的连接点的部分参数和第二回旋曲线在第二轨迹与第二回旋曲线的连接点的部分参数以及第三平滑转接条件带入插补曲线通用方程得到方程组(2):
根据第三平滑转接条件将第二回旋曲线在特征插补点Q的曲率表示为
k1=c0·s0
根据第三平滑转接条件将第二回旋曲线在特征插补点Q的切向角表示为
Figure BDA0001887568300000065
将O点的坐标、特征插补点Q的坐标、第二回旋曲线在特征插补点Q的曲率和切向角带入插补曲线通用方程得到方程组(2)
Figure BDA0001887568300000071
根据方程组(1)和(2)求解Tccw、x0、c0、s0、c1、s1以及α;然后将Tccw、x0、c0、s0、c1、s1以及α带回到插补曲线通用方程得到第一回旋曲线和第二回旋曲线的方程。
请参阅图4,图4是本发明第三实施例运动轨迹平滑转接的方法流程图。在本实施例中,优选地,根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线的步骤包括:拟合出连接第一轨迹和特征插补点的第一子曲线以及连接特征插补点和第二轨迹之间的第二子曲线,第一子曲线和第二子曲线分别满足插补曲线通用方程且具有不同的曲线参数,第一子曲线和第二子曲线在特征插补点处的曲率与切向角均相等。第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等。第一子曲线为第一回旋(clothoid)曲线,第二子曲线为第二回旋曲线。第一轨迹为第一圆弧,第二轨迹为第二圆弧。
在本实施例中,运动轨迹平滑转接的方法包括以下步骤:
步骤S31:获取待平滑转接的第一圆弧和第二圆弧以及预设轮廓误差。
在步骤S31中,例如获取待平滑转接的第一圆弧和第二圆弧包括:获取第一圆弧和第二圆弧以及二者的位置关系。
步骤S32:根据第一圆弧和第二圆弧以及预设轮廓误差确定特征插补点的位置。
在步骤S32中,例如,特征插补点的确定条件为:第一圆弧与第二圆弧在二者连接点处的切线的夹角的角平分线上选择距离第一圆弧和第二圆弧的连接点为预设轮廓误差的点作为特征插补点。
步骤S33:将第一圆弧在未与第二圆弧连接的一端处的切线设置在直角坐标系的一坐标轴上,将第一圆弧未与第二圆弧连接的一端设为与直角坐标系的原点重合。
步骤S34:通过联立以下两组方程求解出第一回旋曲线和第二回旋曲线的曲线参数:
Figure BDA0001887568300000081
Figure BDA0001887568300000082
其中,c0为第一回旋曲线的曲率变化率,s0为第一回旋曲线从其与第一圆弧的连接点到特征插补点之间的弧长,c1为第二回旋曲线的曲率变化率,s1为第二回旋曲线从特征插补点到其与第二圆弧的连接点之间的弧长,α1为第一圆弧从第一圆弧未与第二圆弧连接的一端到第一圆弧与第一回旋曲线的连接点之间的圆弧片段的圆心角,α2为第二圆弧从第二回旋曲线与第二圆弧的连接点到第二圆弧与第一圆弧的连接点之间的圆弧片段的圆心角,Tcw为旋转矩阵;c0、s0、c1、s1、α1、α2、Tcw为需要通过两组方程进行求解的未知量;
其中,R1为第一圆弧的半径,R2为第二圆弧的半径,E为预设轮廓误差,x1为第一圆弧与第二圆弧的连接点在坐标轴上的坐标,y1为第一圆弧与第二圆弧的连接点在直角坐标系另一坐标轴上的坐标,x2为第二圆弧的圆心在坐标轴上的坐标,y2为第二圆弧的圆心在直角坐标系另一坐标轴上的坐标,
Figure BDA0001887568300000083
为第一圆弧在第一圆弧与第二圆弧的连接点的切向量,
Figure BDA0001887568300000084
为第二圆弧在第一圆弧与第二圆弧的连接点的切向量,
Figure BDA0001887568300000085
为第二圆弧的圆心到第一圆弧与第二圆弧的连接点的方向向量,R1、R2、x1、y1、x2、y2
Figure BDA0001887568300000086
为能够通过第一圆弧和第二圆弧及二者的相对位置关系获得的已知量。
下面说明推导过程。
首先,插补曲线的通用方程与第二实施例相同,三种表达方式分别为方程(3)、(4)、(5),此处不再赘述。
其次,进一步说明特征插补点的位置的确定以及方程组(6)和方程组(7)的推导。
记第一轨迹与插补曲线在二者连接点的曲率与切向角相等为第一平滑转接条件,第二轨迹与插补曲线在二者连接点的曲率与切向角相等为第二平滑转接条件,第一子曲线和第二子曲线在特征插补点处的曲率与切向角均相等为第三平滑转接条件。
请参阅图5,图5是第一圆弧、第二圆弧以及插补曲线的几何示意图。
下面以第一圆弧NP和第二圆弧PM的平滑转接为例进行说明。I点为第一圆弧NP与第一回旋曲线IQ的连接点,O点为第二圆弧PM与第二回旋曲线QO的连接点,P点为第一圆弧NP与第二圆弧PM的连接点。关于运动方向,以N点→I点→Q点→O点→M点为例。
设第一圆弧NP上的圆弧片段NI对应的圆心角为α1,设第二圆弧PM上的圆弧片段PO对应的圆心角为α2
在步骤S31中,获取待平滑转接的第一圆弧和第二圆弧以及预设轮廓误差包括:
获取预设轮廓误差E、第一圆弧NP与第二圆弧PM的连接点P、第一圆弧NP在连接点P的切向量
Figure BDA0001887568300000091
第二圆弧PM在连接点P的切向量
Figure BDA0001887568300000092
第一圆弧NP与第二圆弧PM的夹角β、第二圆弧PM的圆心O2、圆心O2到连接点P的方向向量
Figure BDA0001887568300000093
第一圆弧NP的半径R1,第二圆弧PM的半径R2
在步骤S32中,根据前文描述的特征插补点的确定条件为特征插补点Q在第一圆弧NP与第二圆弧PM的切线的夹角的角平分线上且离P点距离QP为E以及根据预设轮廓误差E、夹角β、
Figure BDA0001887568300000094
确定特征插补点Q的位置。
在步骤S33和步骤S34中:
1、先建立直角坐标系,优选地,将第一圆弧NP的切线设置在x轴上,将N点设置为与坐标原点重合,结合步骤S32可以得出在获取的P点在建立的直角坐标系中坐标为(x1,y1)的情况下,特征插补点Q在建立的坐标系中的坐标为:
Figure BDA0001887568300000095
在其他实施例中,直角坐标系x轴与y轴的设置可以调换,例如,将第一圆弧NP的切线设置在y轴上,将N点设置为与坐标原点重合。
2、然后根据第一平滑转接条件和第一轨迹确定第一轨迹与第一回旋曲线的连接点的部分参数和第一回旋曲线在第一轨迹与第一回旋曲线的连接点的部分参数:
由第一圆弧NP的半径R1,可知第一圆弧NP上所有点曲率均为
Figure BDA0001887568300000096
由第一平滑转接条件确定第一回旋曲线在点I的曲率等于第一圆弧NP在I点处的曲率
Figure BDA0001887568300000097
由第一平滑转接条件确定第一回旋曲线在点I的切向角等于第一圆弧NP在I点处的切向角将第一回旋曲线在I点的切向角表示为
θ0=α1
由几何关系将I点的坐标表示为
(R1sinα1,R1(1-cosα1))
3、根据第一回旋曲线满足插补曲线通用方程将第一轨迹与第一回旋曲线连接点的部分参数、第一回旋曲线在第一轨迹与第一回旋曲线连接点的部分参数以及特征插补点的位置带入插补曲线通用方程得到方程组(6):
将第一回旋曲线在点I的曲率和切向角、点I的坐标以及特征插补点的坐标带入插补曲线通用方程得到方程组(6)
Figure BDA0001887568300000101
4、根据第二平滑转接条件和第二轨迹确定第二轨迹与第二回旋曲线的连接点的部分参数和第二回旋曲线在第二轨迹与第二回旋曲线的连接点的部分参数:
根据第二平滑转接条件将第二回旋曲线在O点的切向角表示为
Figure BDA0001887568300000102
其中,atan2表示求反正切函数,只是将角度限制在[-π,π]区间;
将第二回旋曲线在O点的曲率表示为
Figure BDA0001887568300000103
结合步骤S32可以得出在获取的圆心O2在建立的直角坐标系中坐标为(x2,y2)的情况下将点O的坐标表示为
Figure BDA0001887568300000104
其中,T为旋转矩阵
Figure BDA0001887568300000105
5、根据第二回旋曲线满足插补曲线通用方程将特征插补点的位置、第二轨迹与第二回旋曲线的连接点的部分参数和第二回旋曲线在第二轨迹与第二回旋曲线的连接点的部分参数以及第三平滑转接条件带入插补曲线通用方程得到方程组(7):
根据第三平滑转接条件将第二回旋曲线在特征插补点Q的曲率表示为
k1=k0+c0·s0
Figure BDA0001887568300000106
根据第三平滑转接条件将第二回旋曲线在特征插补点Q的切向角表示为
Figure BDA0001887568300000107
Figure BDA0001887568300000111
将点O的坐标、特征插补点Q的坐标、第二回旋曲线在特征插补点Q的曲率和切向角带入插补曲线通用方程得到方程组(7)
Figure BDA0001887568300000112
根据方程组(6)和(7)求解c0、s0、c1、s1、α1、α2、Tcw,然后将c0、s0、c1、s1、α1、α2、Tcw带回插补曲线通用方程得到第一回旋曲线和第二回旋曲线的方程。
请参阅图6,图6是本发明实施例的运动轨迹平滑转接装置的模块示意图。在本实施例中,运动轨迹平滑转接的装置包括处理器61和与处理器61连接的存储器62,存储器62用于存储第一轨迹、第二轨迹、预设轮廓误差,处理器61用于从存储器获取第一轨迹、第二轨迹以及预设轮廓误差;根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等。
其中,处理器61执行的上述各个步骤请参见前文任意一个实施例所述的运动轨迹平滑转接的方法中的说明,此处不再赘述。
请参阅图7,图7是本发明实施例的存储装置的示意图。在本实施例中,存储装置70存储有程序数据,程序数据能够被执行以实现上述任意一个实施例所述的运动轨迹平滑转接的方法。
该存储装置70可以为U盘、光盘、硬盘、移动硬盘、服务器等,当然,存储装置也可以为上述实施例中的存储器62。
请参阅图8,图8是本发明实施例的机器人***的模块示意图。在本实施例中,机器人***包括机器人81以及与机器人81连接的控制器82,控制器82用于获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等。
在一种实施例中,机器人81可包括工具,控制器82可根据第一轨迹、第二轨迹以及插补曲线控制机器人及该工具运动实现加工、运输等操作。在另一实施例中,控制器82可根据第一轨迹、第二轨迹以及插补曲线控制机器人81运动。
其中,控制器82执行的上述各个步骤请参见前文任意一个实施例所述的运动轨迹平滑转接的方法中的说明,此处不再赘述。
请参阅图9,图9是本发明实施例的数控机床的模块示意图。在本实施例中,数控机床包括机床本体91以及与机床本体91连接的数控设备92,数控设备92用于获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置;根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,以使得第一轨迹与插补曲线在二者连接点的曲率与切向角相等,第二轨迹与插补曲线在二者连接点的曲率与切向角相等。
在一种实施例中,机床本体91可包括加工工具,数控设备92可根据第一轨迹、第二轨迹以及插补曲线控制机床本土及该加工工具运动实现加工等操作。
其中,数控设备92执行的上述各个步骤请参见前文任意一个实施例所述的运动轨迹平滑转接的方法中的说明,此处不再赘述。
本发明通过将轮廓误差作为拟合插补曲线的条件,根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置,再根据第一轨迹、第二轨迹、特征插补点以及插补曲线通用方程拟合出连接第一轨迹、特征插补点和第二轨迹的插补曲线,一方面可以准确的拟合出插补曲线,另一方面拟合出的插补曲线可以保证设备加工时的轮廓误差。
以上所述仅为本发明的实施方式,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (9)

1.一种运动轨迹平滑转接的方法,其特征在于,所述方法包括:
获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;
根据所述第一轨迹和所述第二轨迹以及所述预设轮廓误差确定特征插补点的位置;
根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,以使得所述第一轨迹与所述插补曲线在二者连接点的曲率与切向角相等,所述第二轨迹与所述插补曲线在二者连接点的曲率与切向角相等;
其中,所述根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线的步骤包括:拟合出连接所述第一轨迹和所述特征插补点的第一子曲线以及连接所述特征插补点和所述第二轨迹之间的第二子曲线,所述第一子曲线和所述第二子曲线分别满足所述插补曲线通用方程且具有不同的曲线参数,所述第一子曲线和所述第二子曲线在所述特征插补点处的曲率与切向角均相等;所述第一子曲线为第一回旋曲线,所述第二子曲线为第二回旋曲线;
其中,所述第一轨迹为线段,所述第二轨迹为圆弧,其中所述根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线的步骤包括:
将所述线段设置在直角坐标系的一坐标轴上,且所述线段的未与所述圆弧连接的一端设置成与所述直角坐标系的原点重合;
通过联立以下两组方程求解出所述第一回旋曲线和所述第二回旋曲线的曲线参数:
Figure FDA0003504343590000011
Figure FDA0003504343590000021
其中,x0为所述第一回旋曲线与所述线段的连接点在所述坐标轴上的坐标,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述线段的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述圆弧的连接点之间的弧长,α为所述圆弧从所述第二回旋曲线与所述圆弧的连接点到所述圆弧与所述线段的连接点之间的圆弧片段的圆心角,Tccw为旋转矩阵,Tccw、x0、c0、s0、c1、s1以及α为需要通过所述两组方程进行求解的未知量;
其中,L为所述线段的长度,E为所述预设轮廓误差,β为所述线段与所述圆弧的夹角,R为所述圆弧的半径,xc为所述圆弧的圆心在所述坐标轴上的坐标,yc为所述圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000022
为所述圆弧的圆心到所述圆弧与所述线段的连接点的方向向量,xc、yc
Figure FDA0003504343590000023
L、β和R为能够通过所述线段和所述圆弧及二者的相对位置关系获得的已知量;或者,
所述第一轨迹为第一圆弧,所述第二轨迹为第二圆弧,其中所述根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线的步骤包括:
将所述第一圆弧在未与所述第二圆弧连接的一端处的切线设置在直角坐标系的一坐标轴上,将所述第一圆弧未与所述第二圆弧连接的一端设为与所述直角坐标系的原点重合;
通过联立以下两组方程求解出所述第一回旋曲线和所述第二回旋曲线的曲线参数:
Figure FDA0003504343590000024
Figure FDA0003504343590000031
其中,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述第一圆弧的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述第二圆弧的连接点之间的弧长,α1为所述第一圆弧从所述第一圆弧未与所述第二圆弧连接的一端到所述第一圆弧与所述第一回旋曲线的连接点之间的圆弧片段的圆心角,α2为所述第二圆弧从所述第二回旋曲线与所述第二圆弧的连接点到所述第二圆弧与所述第一圆弧的连接点之间的圆弧片段的圆心角,Tcw为旋转矩阵;c0、s0、c1、s1、α1、α2、Tcw为需要通过所述两组方程进行求解的未知量;
其中,R1为所述第一圆弧的半径,R2为所述第二圆弧的半径,E为所述预设轮廓误差,x1为所述第一圆弧与所述第二圆弧的连接点在所述坐标轴上的坐标,y1为所述第一圆弧与所述第二圆弧的连接点在所述直角坐标系另一坐标轴上的坐标,x2为所述第二圆弧的圆心在所述坐标轴上的坐标,y2为所述第二圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000032
为所述第一圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000033
为所述第二圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000034
为所述第二圆弧的圆心到所述第一圆弧与所述第二圆弧的连接点的方向向量,R1、R2、x1、y1、x2、y2
Figure FDA0003504343590000035
为能够通过所述第一圆弧和所述第二圆弧及二者的相对位置关系获得的已知量。
2.根据权利要求1所述的方法,其特征在于,所述根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置包括:
在所述第一轨迹和第二轨迹二者的连接点处的切线夹角的角平分线上选择距离所述第一轨迹和第二轨迹的连接点为所述预设轮廓误差的点作为所述特征插补点。
3.一种运动轨迹平滑转接的装置,其特征在于,所述装置包括处理器和与所述处理器连接的存储器,所述存储器用于存储第一轨迹、第二轨迹、预设轮廓误差,所述处理器用于从所述存储器获取所述第一轨迹、所述第二轨迹以及所述预设轮廓误差;根据所述第一轨迹和所述第二轨迹以及所述预设轮廓误差确定特征插补点的位置;根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,以使得所述第一轨迹与所述插补曲线在二者连接点的曲率与切向角相等,所述第二轨迹与所述插补曲线在二者连接点的曲率与切向角相等;
其中,所述处理器根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,包括:拟合出连接所述第一轨迹和所述特征插补点的第一子曲线以及连接所述特征插补点和所述第二轨迹之间的第二子曲线,所述第一子曲线和所述第二子曲线分别满足所述插补曲线通用方程且具有不同的曲线参数,所述第一子曲线和所述第二子曲线在所述特征插补点处的曲率与切向角均相等;所述第一子曲线为第一回旋曲线,所述第二子曲线为第二回旋曲线;
其中,所述第一轨迹为线段,所述第二轨迹为圆弧时,所述处理器根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,包括:
将所述线段设置在直角坐标系的一坐标轴上,且所述线段的未与所述圆弧连接的一端设置成与所述直角坐标系的原点重合;以及通过联立以下两组方程求解出第一回旋曲线和第二回旋曲线的曲线参数:
Figure FDA0003504343590000041
其中,x0为所述第一回旋曲线与所述线段的连接点在所述坐标轴上的坐标,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述线段的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述圆弧的连接点之间的弧长,α为所述圆弧从所述第二回旋曲线与所述圆弧的连接点到所述圆弧与所述线段的连接点之间的圆弧片段的圆心角,Tccw为旋转矩阵,Tccw、x0、c0、s0、c1、s1以及α为需要通过所述两组方程进行求解的未知量;
其中,L为所述线段的长度,E为所述预设轮廓误差,β为所述线段与所述圆弧的夹角,R为所述圆弧的半径,xc为所述圆弧的圆心在所述坐标轴上的坐标,yc为所述圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000051
为所述圆弧的圆心到所述圆弧与所述线段的连接点的方向向量,xc、yc
Figure FDA0003504343590000052
L、β和R为能够通过所述线段和所述圆弧及二者的相对位置关系获得的已知量;或者,
所述第一轨迹为第一圆弧,所述第二轨迹为第二圆弧时,所述处理器根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,包括:
将所述第一圆弧在未与所述第二圆弧连接的一端处的切线设置在直角坐标系的一坐标轴上,将所述第一圆弧未与所述第二圆弧连接的一端设为与所述直角坐标系的原点重合;通过联立以下两组方程求解出所述第一回旋曲线和所述第二回旋曲线的曲线参数:
Figure FDA0003504343590000053
Figure FDA0003504343590000054
其中,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述第一圆弧的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述第二圆弧的连接点之间的弧长,α1为所述第一圆弧从所述第一圆弧未与所述第二圆弧连接的一端到所述第一圆弧与所述第一回旋曲线的连接点之间的圆弧片段的圆心角,α2为所述第二圆弧从所述第二回旋曲线与所述第二圆弧的连接点到所述第二圆弧与所述第一圆弧的连接点之间的圆弧片段的圆心角,Tcw为旋转矩阵;c0、s0、c1、s1、α1、α2、Tcw为需要通过所述两组方程进行求解的未知量;
其中,R1为所述第一圆弧的半径,R2为所述第二圆弧的半径,E为所述预设轮廓误差,x1为所述第一圆弧与所述第二圆弧的连接点在所述坐标轴上的坐标,y1为所述第一圆弧与所述第二圆弧的连接点在所述直角坐标系另一坐标轴上的坐标,x2为所述第二圆弧的圆心在所述坐标轴上的坐标,y2为所述第二圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000055
为所述第一圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000061
为所述第二圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000062
为所述第二圆弧的圆心到所述第一圆弧与所述第二圆弧的连接点的方向向量,R1、R2、x1、y1、x2、y2
Figure FDA0003504343590000063
为能够通过所述第一圆弧和所述第二圆弧及二者的相对位置关系获得的已知量。
4.根据权利要求3所述的装置,其特征在于,所述处理器根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置包括:
在所述第一轨迹和第二轨迹二者的连接点处的切线夹角的角平分线上选择距离所述第一轨迹和第二轨迹的连接点为所述预设轮廓误差的点作为所述特征插补点。
5.一种存储装置,其特征在于,所述存储装置存储有程序数据,所述程序数据能够被执行以实现如权利要求1-2任意一项所述的方法。
6.一种机器人***,其特征在于,所述机器人***包括机器人以及与所述机器人连接的控制器,所述控制器用于获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;根据所述第一轨迹和所述第二轨迹以及所述预设轮廓误差确定特征插补点的位置;根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,以使得所述第一轨迹与所述插补曲线在二者连接点的曲率与切向角相等,所述第二轨迹与所述插补曲线在二者连接点的曲率与切向角相等;
其中,所述控制器根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线包括:拟合出连接所述第一轨迹和所述特征插补点的第一子曲线以及连接所述特征插补点和所述第二轨迹之间的第二子曲线,所述第一子曲线和所述第二子曲线分别满足所述插补曲线通用方程且具有不同的曲线参数,所述第一子曲线和所述第二子曲线在所述特征插补点处的曲率与切向角均相等;所述第一子曲线为第一回旋曲线,所述第二子曲线为第二回旋曲线;其中,所述第一轨迹为线段,所述第二轨迹为圆弧时,所述控制器根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,包括:
将所述线段设置在直角坐标系的一坐标轴上,且所述线段的未与所述圆弧连接的一端设置成与所述直角坐标系的原点重合;以及通过联立以下两组方程求解出第一回旋曲线和第二回旋曲线的曲线参数:
Figure FDA0003504343590000064
Figure FDA0003504343590000071
其中,x0为所述第一回旋曲线与所述线段的连接点在所述坐标轴上的坐标,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述线段的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述圆弧的连接点之间的弧长,α为所述圆弧从所述第二回旋曲线与所述圆弧的连接点到所述圆弧与所述线段的连接点之间的圆弧片段的圆心角,Tccw为旋转矩阵,Tccw、x0、c0、s0、c1、s1以及α为需要通过所述两组方程进行求解的未知量;
其中,L为所述线段的长度,E为所述预设轮廓误差,β为所述线段与所述圆弧的夹角,R为所述圆弧的半径,xc为所述圆弧的圆心在所述坐标轴上的坐标,yc为所述圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000072
为所述圆弧的圆心到所述圆弧与所述线段的连接点的方向向量,xc、yc
Figure FDA0003504343590000073
L、β和R为能够通过所述线段和所述圆弧及二者的相对位置关系获得的已知量;或者,
所述第一轨迹为第一圆弧,所述第二轨迹为第二圆弧时,所述控制器根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,包括:
将所述第一圆弧在未与所述第二圆弧连接的一端处的切线设置在直角坐标系的一坐标轴上,将所述第一圆弧未与所述第二圆弧连接的一端设为与所述直角坐标系的原点重合;通过联立以下两组方程求解出所述第一回旋曲线和所述第二回旋曲线的曲线参数:
Figure FDA0003504343590000074
Figure FDA0003504343590000081
其中,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述第一圆弧的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述第二圆弧的连接点之间的弧长,α1为所述第一圆弧从所述第一圆弧未与所述第二圆弧连接的一端到所述第一圆弧与所述第一回旋曲线的连接点之间的圆弧片段的圆心角,α2为所述第二圆弧从所述第二回旋曲线与所述第二圆弧的连接点到所述第二圆弧与所述第一圆弧的连接点之间的圆弧片段的圆心角,Tcw为旋转矩阵;c0、s0、c1、s1、α1、α2、Tcw为需要通过所述两组方程进行求解的未知量;
其中,R1为所述第一圆弧的半径,R2为所述第二圆弧的半径,E为所述预设轮廓误差,x1为所述第一圆弧与所述第二圆弧的连接点在所述坐标轴上的坐标,y1为所述第一圆弧与所述第二圆弧的连接点在所述直角坐标系另一坐标轴上的坐标,x2为所述第二圆弧的圆心在所述坐标轴上的坐标,y2为所述第二圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000082
为所述第一圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000083
为所述第二圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000084
为所述第二圆弧的圆心到所述第一圆弧与所述第二圆弧的连接点的方向向量,R1、R2、x1、y1、x2、y2
Figure FDA0003504343590000085
为能够通过所述第一圆弧和所述第二圆弧及二者的相对位置关系获得的已知量。
7.根据权利要求6所述的机器人***,其特征在于,所述控制器根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置包括:
在所述第一轨迹和第二轨迹二者的连接点处的切线夹角的角平分线上选择距离所述第一轨迹和第二轨迹的连接点为所述预设轮廓误差的点作为所述特征插补点。
8.一种数控机床,其特征在于,所述数控机床包括机床本体和数控设备,所述数控设备用于获取待平滑转接的第一轨迹和第二轨迹以及预设轮廓误差;根据所述第一轨迹和所述第二轨迹以及所述预设轮廓误差确定特征插补点的位置;根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,以使得所述第一轨迹与所述插补曲线在二者连接点的曲率与切向角相等,所述第二轨迹与所述插补曲线在二者连接点的曲率与切向角相等;
其中,所述数控设备根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线包括:拟合出连接所述第一轨迹和所述特征插补点的第一子曲线以及连接所述特征插补点和所述第二轨迹之间的第二子曲线,所述第一子曲线和所述第二子曲线分别满足所述插补曲线通用方程且具有不同的曲线参数,所述第一子曲线和所述第二子曲线在所述特征插补点处的曲率与切向角均相等;所述第一子曲线为第一回旋曲线,所述第二子曲线为第二回旋曲线;
其中,所述第一轨迹为线段,所述第二轨迹为圆弧时,所述数控设备根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,包括:
将所述线段设置在直角坐标系的一坐标轴上,且所述线段的未与所述圆弧连接的一端设置成与所述直角坐标系的原点重合;以及通过联立以下两组方程求解出第一回旋曲线和第二回旋曲线的曲线参数:
Figure FDA0003504343590000091
Figure FDA0003504343590000092
其中,x0为所述第一回旋曲线与所述线段的连接点在所述坐标轴上的坐标,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述线段的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述圆弧的连接点之间的弧长,α为所述圆弧从所述第二回旋曲线与所述圆弧的连接点到所述圆弧与所述线段的连接点之间的圆弧片段的圆心角,Tccw为旋转矩阵,Tccw、x0、c0、s0、c1、s1以及α为需要通过所述两组方程进行求解的未知量;
其中,L为所述线段的长度,E为所述预设轮廓误差,β为所述线段与所述圆弧的夹角,R为所述圆弧的半径,xc为所述圆弧的圆心在所述坐标轴上的坐标,yc为所述圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000093
为所述圆弧的圆心到所述圆弧与所述线段的连接点的方向向量,xc、yc
Figure FDA0003504343590000094
L、β和R为能够通过所述线段和所述圆弧及二者的相对位置关系获得的已知量;或者,
所述第一轨迹为第一圆弧,所述第二轨迹为第二圆弧时,所述数控设备根据所述第一轨迹、第二轨迹、所述特征插补点以及插补曲线通用方程拟合出连接所述第一轨迹、所述特征插补点和所述第二轨迹的插补曲线,包括:
将所述第一圆弧在未与所述第二圆弧连接的一端处的切线设置在直角坐标系的一坐标轴上,将所述第一圆弧未与所述第二圆弧连接的一端设为与所述直角坐标系的原点重合;通过联立以下两组方程求解出所述第一回旋曲线和所述第二回旋曲线的曲线参数:
Figure FDA0003504343590000101
Figure FDA0003504343590000102
其中,c0为所述第一回旋曲线的曲率变化率,s0为所述第一回旋曲线从其与所述第一圆弧的连接点到所述特征插补点之间的弧长,c1为所述第二回旋曲线的曲率变化率,s1为所述第二回旋曲线从所述特征插补点到其与所述第二圆弧的连接点之间的弧长,α1为所述第一圆弧从所述第一圆弧未与所述第二圆弧连接的一端到所述第一圆弧与所述第一回旋曲线的连接点之间的圆弧片段的圆心角,α2为所述第二圆弧从所述第二回旋曲线与所述第二圆弧的连接点到所述第二圆弧与所述第一圆弧的连接点之间的圆弧片段的圆心角,Tcw为旋转矩阵;c0、s0、c1、s1、α1、α2、Tcw为需要通过所述两组方程进行求解的未知量;
其中,R1为所述第一圆弧的半径,R2为所述第二圆弧的半径,E为所述预设轮廓误差,x1为所述第一圆弧与所述第二圆弧的连接点在所述坐标轴上的坐标,y1为所述第一圆弧与所述第二圆弧的连接点在所述直角坐标系另一坐标轴上的坐标,x2为所述第二圆弧的圆心在所述坐标轴上的坐标,y2为所述第二圆弧的圆心在所述直角坐标系另一坐标轴上的坐标,
Figure FDA0003504343590000103
为所述第一圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000104
为所述第二圆弧在所述第一圆弧与所述第二圆弧的连接点的切向量,
Figure FDA0003504343590000105
为所述第二圆弧的圆心到所述第一圆弧与所述第二圆弧的连接点的方向向量,R1、R2、x1、y1、x2、y2
Figure FDA0003504343590000111
为能够通过所述第一圆弧和所述第二圆弧及二者的相对位置关系获得的已知量。
9.根据权利要求8所述的数控机床,其特征在于,所述数控设备根据第一轨迹和第二轨迹以及预设轮廓误差确定特征插补点的位置包括:
在所述第一轨迹和第二轨迹二者的连接点处的切线夹角的角平分线上选择距离所述第一轨迹和第二轨迹的连接点为所述预设轮廓误差的点作为所述特征插补点。
CN201780034014.7A 2017-05-12 2017-05-12 运动轨迹平滑转接的方法和装置及其相关设备 Active CN109313428B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2017/084199 WO2018205276A1 (zh) 2017-05-12 2017-05-12 运动轨迹平滑转接的方法和装置及其相关设备

Publications (2)

Publication Number Publication Date
CN109313428A CN109313428A (zh) 2019-02-05
CN109313428B true CN109313428B (zh) 2022-04-29

Family

ID=64104275

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201780034014.7A Active CN109313428B (zh) 2017-05-12 2017-05-12 运动轨迹平滑转接的方法和装置及其相关设备

Country Status (2)

Country Link
CN (1) CN109313428B (zh)
WO (1) WO2018205276A1 (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111216119B (zh) * 2018-11-23 2022-04-15 深圳市优必选科技有限公司 机器人关节运动控制方法、装置及终端设备
TWI716127B (zh) * 2019-10-01 2021-01-11 東元電機股份有限公司 機器人及其路徑插值規劃命令產生系統
CN112621739B (zh) * 2019-10-08 2022-03-15 东元电机股份有限公司 机器人及其路径插值规划命令产生***
CN113450262A (zh) * 2020-03-25 2021-09-28 武汉金山办公软件有限公司 一种实现图像处理的方法、装置、计算机存储介质及终端
CN111958611B (zh) * 2020-07-31 2021-10-08 华中科技大学 考虑轮廓误差最小的铣削加工机器人姿态优化方法及设备
CN113103240B (zh) * 2021-04-29 2022-08-23 哈尔滨工业大学 用于实现c2连续的机器人轨迹规划方法、装置及***
CN113504764B (zh) * 2021-06-30 2024-07-05 浙江大学 基于位置矢量加权积分的连续线段数控加工路径平滑方法
CN113791578B (zh) * 2021-08-23 2023-05-02 五邑大学 基于数控加工***的轨迹滤波方法、装置和电子设备
CN113741338B (zh) * 2021-09-10 2023-08-22 深圳市英威腾电气股份有限公司 一种运动设备的运行路径控制方法、***及装置
CN114434241B (zh) * 2021-12-30 2023-06-02 西门子(中国)有限公司 控制石磨机横梁运动的方法、装置、电子设备及存储介质
CN115128996A (zh) * 2022-07-04 2022-09-30 深圳市汇川技术股份有限公司 数控加工的误差补偿方法、装置、设备及介质
CN116501071B (zh) * 2023-06-30 2023-09-05 苏州铼钠克信息技术有限公司 近似回折路径的规划方法和装置、电子设备和存储介质
CN117944057B (zh) * 2024-03-26 2024-06-21 北京云力境安科技有限公司 一种机械臂轨迹规划方法、装置、设备及介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1398680A2 (en) * 1999-11-19 2004-03-17 Siemens Energy & Automation, Inc. Apparatus and method for smooth cornering in a motion control system
CN103699056A (zh) * 2013-12-02 2014-04-02 嘉兴学院 高速高精度数控加工的小线段实时平滑过渡插补方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0764622A (ja) * 1993-08-26 1995-03-10 Matsushita Electric Ind Co Ltd ロボットの軌道補間装置
JPH10320026A (ja) * 1997-05-20 1998-12-04 Mitsubishi Electric Corp 数値制御装置及び方法
JP5653972B2 (ja) * 2012-07-19 2015-01-14 ファナック株式会社 コーナ複数曲線挿入部を有する数値制御装置
CN106527351B (zh) * 2015-09-10 2019-01-22 沈阳高精数控智能技术股份有限公司 一种过渡曲线长度动态确定的小线段实时插补方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1398680A2 (en) * 1999-11-19 2004-03-17 Siemens Energy & Automation, Inc. Apparatus and method for smooth cornering in a motion control system
CN103699056A (zh) * 2013-12-02 2014-04-02 嘉兴学院 高速高精度数控加工的小线段实时平滑过渡插补方法

Also Published As

Publication number Publication date
WO2018205276A1 (zh) 2018-11-15
CN109313428A (zh) 2019-02-05

Similar Documents

Publication Publication Date Title
CN109313428B (zh) 运动轨迹平滑转接的方法和装置及其相关设备
CN107615194B (zh) 加工轨迹平滑转接的方法及加工装置
US8831777B2 (en) Robot program changing device
CN107980109B (zh) 机器人运动轨迹规划方法及相关装置
CN108594815B (zh) 一种分阶段的轮式机器人移动路径规划方法
CN109605369B (zh) 机械臂奇异点控制方法及***
WO2017113195A1 (zh) 一种加工路径规划方法、加工路径规划装置及数控机床
CN104635762A (zh) 一种面向srs仿人手臂的自运动角计算方法
CN109794943B (zh) 一种拐角过渡路径及确定方法
CN115502966B (zh) 用于机器人的可变导纳控制方法
CN108062071B (zh) 参数曲线轨迹伺服轮廓误差的实时测定方法
US10807238B2 (en) Robot system and method for controlling robot
CN103454979B (zh) 一种封装成PLCOpen指令的变速曲线圆弧快速插补方法
CN108803480B (zh) 轨迹间g2连续光顺转接方法、设备及计算机可读存储介质
JP4577247B2 (ja) 移動体及びその制御方法
CN111670093B (zh) 机器人的运动控制方法、控制***和存储装置
WO2023207164A1 (zh) 一种机器人作业控制方法及装置
WO2023124339A1 (zh) 路径规划方法、运动控制方法及计算机程序产品
CN105676787A (zh) 一种椭圆弧插补算法
CN109648563B (zh) 串联机器人运动控制方法及计算机存储介质
JPS5916286B2 (ja) 工業用ロボツトの動作制御方法
CN111405966B (zh) 用于控制机器人组的方法和控制装置
CN116494250A (zh) 基于速度补偿的机械臂控制方法、控制器、介质及***
US20210260761A1 (en) Method And Control System For Controlling An Industrial Actuator
CN110625617B (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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 518063 23 Floor (Room 2303-2306) of Desai Science and Technology Building, Yuehai Street High-tech Zone, Nanshan District, Shenzhen City, Guangdong Province

Applicant after: Shenzhen AANDE Intelligent Technology Research Institute Co., Ltd.

Address before: 518104 Shajing Industrial Co., Ltd. No. 3 Industrial Zone, Hexiang Road, Shajing Street, Baoan District, Shenzhen City, Guangdong Province

Applicant before: Shenzhen AANDE Intelligent Technology Research Institute Co., Ltd.

GR01 Patent grant
GR01 Patent grant