本发明涉及一种多轴机器人7r机械臂逆解建模与解算方法,属于机器人技术领域。
背景技术
自主机器人研究的一个重要方面是需要解决变拓扑结构机器人的运动学建模问题。在mas中,具有动态的图结构(dynamicgraphstructure),可以动态地建立基于运动轴的有向span树,为研究可变拓扑结构(variabletopologystructure)的机器人建模与控制奠定了基础。为此,需要提出基于轴不变量的通用机械臂逆解原理,既要建立包含坐标系、极性、结构参数、关节变量的完全参数化的正运动学模型,又要实时地计算位姿方程;一方面,可以提高机器人的自主性,另一方面,可以提高机器人位姿控制的绝对精度。
由于通用6r机械臂不存在共点约束,现有技术中其逆解计算已经十分困难,在工程上不得不屈从于解耦约束,该约束既增加了机械臂加工及装配难度,又降低了机械臂绝对定位精度。通用7r机械臂逆解由于计算复杂度极高,在现有技术条件下根本无法实现。
技术实现要素:
本发明所要解决的技术问题是提供一种基于轴不变量的通用7r机械臂逆解建模与解算方法,解决了现有技术中无法计算7r机械臂逆解的问题。
为解决上述技术问题,本发明采用以下技术方案:
一种基于轴不变量的通用7r机械臂逆解建模与解算方法,其特征是,
设定有7个转动轴,拾取点位于第7轴轴线上,且第5轴与第6轴不共轴的机械臂为通用7r机械臂;通用7r机械臂通过前6轴控制第7轴与期望位置及姿态对齐,使第7轴无限转动或控制第7轴满足径向对齐;
将给定第7轴期望位置
基于轴不变量的通用机械臂运动规划,建立通用7r机械臂增量位姿方程,第6轴取距拾取点一定距离的点为名义拾取点,先计算通用6r机械臂的逆解,再应用数值迭代法,完成通用7r机械臂的运动规划与逆解计算。
若给定7r轴链il7=(i,1:7],il1=03,期望位置矢量及ju-gibbs四元数分别记为
其中:
系统结构参数及期望ju-gibbs姿态四元数构成的矩阵表示为
式中,\为续行符;
消去τ5及τ6后的位姿方程(200),是4个“4元2阶”多项式方程。
基于轴不变量的通用机械臂运动规划,具体包括以下步骤:
【1】建立通用7r轴链的增量位姿方程;
【2】对基于偏速度迭代的通用7r机械臂进行运动规划。
定义ju-gibbs增量四元数:
ju-gibbs规范四元数为
其中:轴不变量
定义居-吉布斯增量四元数
其中:
步骤【1】中,若给定6r轴链il7=(i,1:7],il1=03,位置矢量及ju-gibbs增量四元数分别记为
其中:
由行四元数构成的系统结构参数矩阵表示为
其中,
分析通用7r机械臂的增量位姿逆解,式(239)是关于{εl|l∈[1:4}的线性方程。将式(239)重新表示为
a·[ε1ε2ε3ε4]t=b;(247)
若a-1存在,解式(245)得
[ε1ε2ε3ε4]t=a-1·b;(248)
至此,得到全部逆解。
通过增量位置矢量
步骤【2】中,基于偏速度迭代优化的步骤如下:
记运动链为
记
由式(252)应用梯度下降法得
其中:步长step>0,step→0;有
选择step步长,由初态
(1)确定目标函数
goal表示
(2)选择步长;
(3)进行迭代计算,得到稳态解,即为通用机械臂的位姿逆解。
应用构造法确定步长
仅当
其中:ε(i,6]=[ε1ε2…ε6];由式(257)及式(253)确定步长
其中:
迭代过程
取式(258)之步长step,由式(259)完成迭代计算:
其中:当
对于式(259)的迭代过程,则必有
δgoal≤0,(260)
即式(259)的迭代过程一定收敛;
当
本发明所达到的有益效果:
本发明建立了通用7r机械臂增量位姿方程,并进行实时逆解计算。特征在于:
具有简洁、优雅的运动链符号系统,具有伪代码的功能,具有迭代式结构,保证系统实现的可靠性及机械化演算。
具有基于轴不变量的迭代式,保证计算的实时性;实现坐标系、极性及系统结构参量的完全参数化,基于轴不变量的可逆解运动学具有统一的表达及简洁的结构化层次模型,保证位姿分析逆解的通用性。
直接应用激光跟踪仪精密测量获得的基于固定轴不变量的结构参数,保证位姿逆解的准确性;从而,使系统的绝对定位与定姿精度接近重复精度。
附图说明
图1自然坐标系与轴链;
图2固定轴不变量;
图3为定轴转动示意图;
图4为轴不变量的导出不变量。
具体实施方式
下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
定义1自然坐标轴:称与运动轴或测量轴共轴的,具有固定原点的单位参考轴为自然坐标轴,亦称为自然参考轴。
定义2自然坐标系:如图1所示,若多轴系统d处于零位,所有笛卡尔体坐标系方向一致,且体坐标系原点位于运动轴的轴线上,则该坐标系统为自然坐标系统,简称自然坐标系。
自然坐标系优点在于:(1)坐标系统易确定;(2)零位时的关节变量为零;(3)零位时的系统姿态一致;(4)不易引入测量累积误差。
由定义2可知,在系统处于零位时,所有杆件的自然坐标系与底座或世界系的方向一致。系统处于零位即
由上式知,
轴不变量与坐标轴具有本质区别:
(1)坐标轴是具有零位及单位刻度的参考方向,可以描述沿该方向平动的位置,但不能完整描述绕该方向的转动角度,因为坐标轴自身不具有径向参考方向,即不存在表征转动的零位。在实际应用时,需要补充该轴的径向参考。例如:在笛卡尔系f[l]中,绕lx转动,需以ly或lz为参考零位。坐标轴自身是1d的,3个正交的1d参考轴构成3d的笛卡尔标架。
(2)轴不变量是3d的空间单位参考轴,其自身就是一个标架。其自身具有径向参考轴,即参考零位。空间坐标轴及其自身的径向参考轴可以确定笛卡尔标架。空间坐标轴可以反映运动轴及测量轴的三个基本参考属性。
已有文献将无链指标的轴矢量记为
【1】给定旋转变换阵
【2】是3d参考轴,不仅具有轴向参考方向,而且具有径向参考零位,将在3.3.1节予以阐述。
【3】在自然坐标系下:
对轴不变量而言,其绝对导数就是其相对导数。因轴不变量是具有不变性的自然参考轴,故其绝对导数恒为零矢量。因此,轴不变量具有对时间微分的不变性。有:
【4】在自然坐标系统中,通过轴矢量
【5】应用轴矢量
因基矢量el是与f[l]固结的任一矢量,基矢量
在系统处于零位时,以自然坐标系为参考,测量得到坐标轴矢量
因此,应用自然坐标系统,当系统处于零位时,只需确定一个公共的参考系,而不必为系统中每一杆件确定各自的体坐标系,因为它们由轴不变量及自然坐标唯一确定。当进行系统分析时,除底座系外,与杆件固结的其它自然坐标系只发生在概念上,而与实际的测量无关。自然坐标系统对于多轴系统(mas)理论分析及工程作用在于:
(1)系统的结构参数测量需要以统一的参考系测量;否则,不仅工程测量过程烦琐,而且引入不同的体系会引入更大的测量误差。
(2)应用自然坐标系统,除根杆件外,其它杆件的自然坐标系统由结构参量及关节变量自然确定,有助于mas系统的运动学与动力学分析。
(3)在工程上,可以应用激光跟踪仪等光学测量设备,实现对固定轴不变量的精确测量。
(4)由于运动副r及p、螺旋副h、接触副o是圆柱副c的特例,可以应用圆柱副简化mas运动学及动力学分析。
定义3不变量:称不依赖于一组坐标系进行度量的量为不变量。
定义4转动坐标矢量:绕坐标轴矢量
定义5平动坐标矢量:沿坐标轴矢量
定义6自然坐标:以自然坐标轴矢量为参考方向,相对系统零位的角位置或线位置,记为ql,称为自然坐标;称与自然坐标一一映射的量为关节变量;其中:
定义7机械零位:对于运动副
故关节
定义8自然运动矢量:将由自然坐标轴矢量
自然运动矢量实现了轴平动与转动的统一表达。将由自然坐标轴矢量及关节确定的矢量,例如
定义9关节空间:以关节自然坐标ql表示的空间称为关节空间。
定义10位形空间:称表达位置及姿态(简称位姿)的笛卡尔空间为位形空间,是双矢量空间或6d空间。
定义11自然关节空间:以自然坐标系为参考,通过关节变量
如图2所示,给定链节
轴矢量
定义12自然坐标轴空间:以固定轴不变量作为自然参考轴,以对应的自然坐标表示的空间称为自然坐标轴空间,简称自然轴空间。它是具有1个自由度的3d空间。
如图2所示,
给定多轴系统d={t,a,b,k,f,nt},在系统零位时,只要建立底座系或惯性系,以及各轴上的参考点ol,其它杆件坐标系也自然确定。本质上,只需要确定底座系或惯性系。
给定一个由运动副连接的具有闭链的结构简图,可以选定回路中任一个运动副,将组成该运动副的定子与动子分割开来;从而,获得一个无回路的树型结构,称之为span树。t表示带方向的span树,以描述树链运动的拓扑关系。
i为结构参数;a为轴序列,f为杆件参考系序列,b为杆件体序列,k为运动副类型序列,nt为约束轴的序列即非树。
描述运动链的基本拓扑符号及操作是构成运动链拓扑符号系统的基础,定义如下:
【1】运动链由偏序集合(]标识。
【2】a[l]为取轴序列a的成员;因轴名l具有唯一的编号对应于a[l]的序号,故a[l]计算复杂度为o(1)。
【3】
【4】
【5】llk为取由轴l至轴k的运动链,输出表示为
【6】ll为取轴l的子。该操作表示在
【7】ll表示获得由轴l及其子树构成的闭子树,ll为不含l的子树;递归执行ll,计算复杂度为
【8】支路、子树及非树弧的增加与删除操作也是必要的组成部分;从而,通过动态span树及动态图描述可变拓扑结构。在支路llk中,若
定义以下表达式或表达形式:
轴与杆件具有一一对应性;轴间的属性量
约定:
本申请中约定:在运动链符号演算系统中,具有偏序的属性变量或常量,在名称上包含表示偏序的指标;要么包含左上角及右下角指标,要么包含右上角及右下角指标;它们的方向总是由左上角指标至右下角指标,或由右上角指标至右下角指标,本申请中为叙述简便,有时省略方向的描述,即使省略,本领域技术人员通过符号表达式也可以知道,本申请中采用的各参数,对于某种属性符,它们的方向总是由偏序指标的左上角指标至右下角指标,或由右上角指标至右下角指标。例如:
本申请的符号规范与约定是根据运动链的偏序性、链节是运动链的基本单位这两个原则确定的,反映了运动链的本质特征。链指标表示的是连接关系,右上指标表征参考系。采用这种符号表达简洁、准确,便于交流与书面表达。同时,它们是结构化的符号系统,包含了组成各属性量的要素及关系,便于计算机处理,为计算机自动建模奠定基础。指标的含义需要通过属性符的背景即上下文进行理解;比如:若属性符是平动类型的,则左上角指标表示坐标系的原点及方向;若属性符是转动类型的,则左上角指标表示坐标系的方向。
(1)ls-杆件l中的点s;而s表示空间中的一点s。
(2)
(3)
(4)
(5)
(6)
(7)左下角指标为0时,表示机械零位;如:
(8)0-三维零矩阵;1-三维单位矩阵;
(9)约定:“\”表示续行符;
幂符
叉乘符运算的优先级高于投影符
(10)单位矢量在大地坐标系的投影矢量
(11)
(12)iql,相对绝对空间的旋转变换阵;
(13)以自然坐标轴矢量为参考方向,相对系统零位的角位置或线位置,记为ql,称为自然坐标;关节变量
(14)对于一给定有序的集合r=[1,4,3,2]t,记r[x]表示取集合r的第x行元素。常记[x]、[y]、[z]及[ω]表示取第1、2、3及4列元素。
(15)ilj表示由i到j的运动链;llk为取由轴l至轴k的运动链;
给定运动链
(16)rodrigues四元数表达形式:
欧拉四元数表达形式:
不变量的四元数(也称为轴四元数)表达形式
分块矩阵的高维行列式计算:
记<1:n>表示自然数[1:n]的全排列,共有n!个实例。给定属于数域的大小为n×n的矩阵m,其j行i列元素记为
其中:i[i1,…in]表示排列<i1,…in>的逆序个数。式(1)计算复杂度为:n!次n个数积及n!次加法,具有指数计算复杂度,只能适用于维度较小的行列式。对于维度较大的行列式,通常应用laplace公式进行递规运算,记
更简单的算法通常应用高斯消去法或lu分解法,先通过初等变换将矩阵变为三角阵或三角阵的乘积,后计算行列式。上述针对数域的行列式计算方法不适用于高维度的多项式矩阵,需要引入分块矩阵的行列式计算方法。计算矢量多项式(vectorpolynomial)的行列式是一个特定的分块矩阵行列式的计算问题,它在矢量层次上表达了矢量与行列式的内在联系。而分块矩阵行列式计算则从矩阵层次上表达分块矩阵与行列式的内在规律。
若给定矢量多项式
则有
上式的推导步骤为:因
故式(4)成立。
式(3)及式(4)可以推广至n维空间。式(3)有助于从矢量层次上分析行列式的内在规律;比如,当任两矢量平行或三矢量共面时,对应的行列式为零。式(4)表明:矢量多项式的行列式易导致“组合爆炸”。
实施例1
给定2个2维行矢量多项式
另一方面,
上面的结果验证了式(4)的正确性。
给出分块矩阵的行列式计算定理:
若记大小为(n+m)·(n+m)的方阵为m,大小为n·n的矩阵
对行列式进行行阶梯化计算原理:
对于s×s矩阵,其每一项是关于τ1的n阶多项式。计算该矩阵的行列式时,可通过初等行变换将原行列式变为上三角行列式,再将非零的对角线元素相乘,得到行列式的多项式表达式。因该式为0,求得到τ1的所有解。
阶梯化的具体方法为,先对行列式第一列的最高阶次由高到低进行排序,再进行最多(s-1)×n次初等行变换消元,得到第一列只有第一个元素不为0的行列式。再对该行列式第1行及1列的余子式进行初等行变换消元,依次迭代求解。
实施例2
通过矩阵的初等行变换,得到
步骤为:rk代表第k行。得
则得
基于“n进位字”的n阶多项式系统:
若n个“n元1阶”多项式幂积
n个“n元n阶”多项式系统的dixon多项式:
引入辅助变量[y2,y3,…,yn],且有
在多元多重多项式(8)中,用辅助变量ym的前m个依次替换原变量(originalvariables)xn中的m个变量,记“|”为替换操作符,得到增广的(extended)多项式
式中右上角标α、α表示幂;
由式(6)及式(12)得
其中:
定义可分离组合变量
由式(14)及式(15)可知:替换式
给定n个“n元n阶”多项式系统
由式(17)得
式(15)中分离变量与文献不同:原变量xn-1被辅助变量yn-1替换的次序不同,dixon多项式也不同。考虑式(13)及式(18)得该多项式的dixon行列式
在笛卡尔空间下,由位置矢量或转动矢量构成的行列式表示矢量张成空间的容积(volume);在不同笛卡尔空间下具有容积的不变性。其中:
n个“n元n阶”多项式的dixon行列式的阶次及替换变量项数分别为:
n个“n元n阶”dixon矩阵:
给定n个“n元n阶”多项式系统fn(yn-1|xn-1),n≥2;存在与消去变量x2,…,xn无关的dixon矩阵sθs(x1),其dixon多项式
α[l]∈[0,n·(n-l+1)-1]],l∈[2:n](23)
其中:
若
则有
n个“n元n阶”多项式的dixon消元与求解
考虑式(22),若
det(sθs(x1))=0。(28)
称式(28)中“n个n元”为dixon消元的必要条件,从而获得可行解。若sθs存在零行或零列向量,则无法建立x1的多项式方程;此时,通过除标量积之外的初等变换,将sθs变为行阶梯(rowechelon)矩阵ech(sθs);在计算该矩阵的杻轴(pivot)的积之后得方阵
任一个n个“n元n阶”多项式系统
确定双线性型
其中:
称其为结式或消去式。式(32)是单变量x1的多项式方程;消去了n-1个未知量;从而,可以获得单变量x1的可行解。若x1同时满足
则x1为正确解。将已解的x1代入式(34),因式(32)成立且
即有
若有必要条件
成立,解式(35),得被消去变量
若同时满足
则由式(35)能解得
给定n个“n元n阶”多项式
①确定系统结构。方程数及独立变量数记为n;独立变量记为xn;多项式复合变量记为
②由式(8)得xα与
③根据式(19)及sarrus规则,计算dixon(fn(yn-1|xn-1));根据
④dixon矩阵成员如式(32)所示,由式(32)计算dixon矩阵sθs的(n+1)·s2个系数。
⑤当满足式(37)及式(38)直接解判别准则时,由式(34)及式(35)得全部数值解。
实施例3
对多项式系统(39)进行dixon消元。
步骤为:该式是4个“4元1阶”多项式系统,满足dixon消元条件。由式(19)及式(22),得
其中:
由式(34)及式(40)得5个解:
其中:
解得:τ3=1,τ4=-2。将
基于轴不变量的定轴转动
如图3所示,给定轴矢量
轴矢量
因单位矢量
若
式(43)是关于
对于轴链
由式(44)及式(43)得
称(45)为改进的cayley变换。即有
规范的位置方程为
“居-吉布斯”四元数的确定:
对于任意杆件
其中:
其中:
显然,
其中:
由式(52)得
习惯上,单关节及运动链的期望姿态以规范的ju-gibbs四元数(简称规范ju-gibbs四元数,即“标部”为1的四元数)表示;但是它们积运算通常是不规范的,即其标部不为1。由式(53)可知:只有给定轴l及
由式(53)得
由四维复数性质得
记
故
由式(48)至式(50)及式(55)得
由式(50)、式(54)及式(57)得
类dcm及性质:
对于轴链
由式(59)得
式中,
其中:
由式(61)可知:iqn及
显然,类dcm可以通过ju-gibbs四元数表达。因此,式(59)姿态方程及式(47)位置方程是关于ju-gibbs四元数的表达式。
分块方阵的逆:
若给定可逆方阵k、b及c,其中b及c分别为l×l、c×c的方阵;a、d分别为l×c、c×l的矩阵,且
则有
基于ju-gibbs四元数的指向对齐原理
考虑轴链ill,其中
其中:
且有
上式的具体建立步骤为:
定轴转动的cayley正变换为
由式(69)得
由式(70)得
及
因
式(73)表明
及
由式(74)及(75)得式(66)。若
由式(76)得
因
及
由式(78)及式(79)分别得式(67)及式(68)。证毕。
基于ju-gibbs四元数的指向对齐原理表明:至少存在一个期望的ju-gibbs四元数
实施例4
考虑轴链il6,由基于ju-gibbs四元数的指向对齐原理得
基于ju-gibbs四元数的2r方向逆解
以ju-gibbs四元数指向对齐为基础,阐述2r方向逆解。
若给定6r转动链il6=(i,1:6]],记第5轴关节ju-gibbs方向四元数期望为
其中:
ju-gibbs方向四元数
上式的具体建立步骤为:
首先,考虑基于欧拉四元数的姿态对齐。由
得
由式(86)得
其中:
由式(64)及式(88)得
由式(87)得
其中:
由式(89)及式(90)得
由式(92)得
若
其次,考虑ju-gibbs四元数的指向对齐。因
以规范ju-gibbs四元数表征关节变量,由式(53)得
由式(96)得
因3n4及4n5独立,由式(83)可知,3e5必存在。显然,
若
将式(99)代入式(98)得
由式(94)及式(100)可知两种原理等价。由式(100)第2、3行得
由式(101)可知式(80)成立。因式存(100)在4个等式,2个独立变量,由式(101)及式(100)中第4行得约束方程
若
由式(103)得
显然,
若
由式(106)可知:要么
ju-gibbs方向四元数
由式(80),(81)和(99)得
且
把
实施例5
继实施例4,考虑轴链il6,且有3n4=1[x],4n5=1[y],由式(83)得3e5=1。由式(80)及式(81)得
由式(81)得:
基于ju-gibbs四元数的2r方向逆解原理表明ju-gibbs四元数与欧拉四元数同构;同时,式(63)表明以ju-gibbs四元数表示的类dcm与dcm同构。因此,应用ju-gibbs四元数可以完整表达位姿关系。
若给定运动链iln,k∈iln,期望规范ju-gibbs四元数
且具有模不变性
与欧拉四元数及对偶四元数相比,ju-gibbs四元数表征的位姿对齐不存在冗余方程;通过指向对齐,可以求解第4轴及第5轴的关节变量,为6r及7r机械臂逆解奠定了基础。
基于轴不变量的dixon行列式计算原理:
下面基于轴不变量,提出径向不变量及运动链的dixon行列式基本性质,为基于轴不变量的机器人逆运动学分析奠定基础。
【1】轴不变量
首先,轴不变量与坐标轴具有本质区别:坐标轴是具有零位及单位刻度的参考方向,可以描述沿轴向平动的线位置,但不能完整描述绕轴向的角位置,因为坐标轴自身不具有径向参考方向,即不存在表征转动的零位。在实际应用时,需要补充坐标轴的径向参考。坐标轴自身是1d的,3个正交的坐标轴构成3d的笛卡尔标架;轴不变量是3d空间单位参考轴(简称3d参考轴),具有径向参考零位。“3d参考轴”及其径向参考零位可以确定对应的笛卡尔系。以自然坐标系为基础的轴不变量可以准确地反映运动轴及测量轴的“共轴性”、“极性”与“零位”三个基本属性。
其次,轴不变量与欧拉轴具有本质的区别:方向余弦矩阵(dcm)是实矩阵,轴矢量是dcm的特征值1对应的特征矢量,是不变量;固定轴不变量是“3d参考轴”,不仅具有原点及轴向,也有径向参考零位;在自然坐标系下,轴不变量不依赖于相邻固结的自然坐标系,即在相邻固结的自然坐标系下具有不变的自然坐标;轴不变量具有幂零特性等优良的数学操作功能;在自然坐标系统中,通过轴不变量及关节坐标,可以唯一确定dcm及参考极性;没有必要为每一个杆件建立各自的体系,可以极大地简化建模的工作量。
同时,以唯一需要定义的笛卡尔直角坐标系为参考,测量轴不变量,可以提高结构参数的测量精度。基于轴不变量的优良操作及属性,可以建立包含拓扑结构、坐标系、极性、结构参量及动力学参量的迭代式的运动学及动力学方程。
由式(59)及式(47)可知:多轴系统的姿态及位置方程本质上是多元二阶多项式方程,其逆解本质上归结于多元二阶多项式的消元问题,包含dixon矩阵及dixon行列式计算的两个子问题。用式(47)的表达3r机械臂位置方程,是3个“3元2阶”多项式,应用dixon消元方法计算逆解,有两个替换变量,在计算8×8的dixon行列式时,最大可能的阶次为16。由式(4)可知:行列式计算是一个排列过程,面临着“组合爆炸”的难题。
所有的不在确定的多项式时间内可解的问题称为np问题。非确定性算法将问题分解为“猜测”与“验证”两个阶段:算法的“猜测”阶段具有非确定性,算法的“验证”阶段具有确定性,通过验证来确定猜测的解是否正确。假如可以在多项式时间内计算出来,就称为多项式非确定性问题。多元多项式的消元通常被认为是np问题。通常应用
【2】径向不变量
结构参数
任一个矢量可以分解为零位矢量及轴向矢量,故有
其中:
考虑链节
显然,
若两个零位不变量
因此,称式(110)所示的轴向不变量及零位不变量是结构参数矢量对自然轴的分解。
由式(113)及式(114)可知:同一个轴的三个径向矢量的行列式为零;同一个轴的任意两个轴向矢量的行列式为零。可以用轴不变量及其导出的不变量来简化dixon行列式计算。
由轴不变量导出的零位矢量、径向矢量及轴向矢量具有以下关系:
称式(115)为零位矢量的反转公式;称式(116)为零位矢量与径向矢量的互换公式;称式(117)为径向矢量不变性公式。由式(109)、式(115)至式(117)得
由式(118)得
因
【3】运动链dixon行列式性质
定义
由式(52)得
其中:
考虑式(122),若m为4·4的矩阵,则有
且有
由式(62)及式(110)得
由式(126)证得
由式(129)得
式(127)可以将
式(129)由三个导出的独立结构参量
由式(127)及式(131)得
由式(127)及式(132)得
基于类dcm的2r方向逆解
给定6r轴链il6=(i,1:6],轴矢量3n4及4n5,期望第5轴的dcm为
上述方程的建立步骤为:
方向矢量5l6与期望方向
由式(60)得
即
式(137)重新表示为式(134)。
基于轴不变量的cayley变换
当给定角度
由式(138)得
定义
故有
一、基于轴不变量的3r机械臂位置逆解方法
给定3r转动链
【1】根据机械臂n元3d矢量位姿方程,获得n个“n元2阶”多项式方程。
由式(47)得3r运动学方程
由式(142)得
由式(143)得
若记
则由式(61)及式得式(145)
由式(144)及式(145)得
下面,阐述3r机械臂运动学方程的dixon行列式的结构模型及特点。由式(147)得3r运动学多项式方程
多项式系统f3(y2|t2),根据双线性型行列式通式
则有
其中:
由式(18)、式(147)及式(148)得
由式(22)及式(153)可知式(151)成立。由式(127)及式(145)得
由式(145)、式(154)及式(155)得
其中:应用式(133)计算
显然,式(156)中的y2阶次β2∈[0:3]及y3阶次β3∈[0:1]。考虑式(153)后三项:
由式(145)、式(153)至式(156)可知:
【2】应用“基于轴不变量的dixon行列式计算”方法及“分块矩阵的高维行列式计算”方法简化行列式计算。
根据运动链dixon行列式性质,由式(127)及式(145)得
由式(157)得
由式(158)得
由式(159)得
由式(153)得
将式(160)至式(162)代入式(163)得
【3】应用n个“n元n阶”多项式的dixon消元与求解原理完成位姿逆解计算,其中:根据dixon矩阵的行列式为0,得到一元高阶多项式方程,应用基于友阵的一元高阶多项式方程求解一元高阶多项式方程的解。
一元n阶多项式p(x)=a0+a1x+…an-1xn-1+xn具有n个解。若能找到一个矩阵a,满足|a-λl·1n|·vl=0,其中:l∈[1:n],λl为该矩阵的特征值,vl为对应的特征矢量。若矩阵a的特征方程为
若多项式p(x)的友阵为
则由矩阵a的特征向量构成的矩阵为范德蒙德(vandermonde)矩阵为
且有
p(λl)=|a-λl·1n|=0(167)。
由式(28)、式(151)及式(152)得
因s=8,应用式(1)计算
二、基于轴不变量的通用6r机械臂位姿逆解方法
设定有6个转动轴,拾取点位于第6轴轴线上,且第4轴与第5轴不共轴的机械臂为通用6r机械臂。通过前5轴控制第6轴与期望的位置及指向对齐,第6轴可以无限转动或控制第6轴满足径向对齐。
基于轴不变量的通用6r机械臂位姿逆解:
给定第6轴期望位置矢量
若给定6r轴链il6=(i,1:6],il1=03,第6轴期望位置矢量为
其中
式中,\为续行符;
系统结构参数及期望ju-gibbs姿态四元数构成的矩阵表示为
其中,
并进行如下定义:
下面给出上式的建立过程:
期望
由式(175)得
由式(52)得式(170),由式(176)得
由式(53)得
其中:
由式(177)及式(178)得
由式(64)得
将式(181)代入式(177)得
其中:
由式(183)得式(172)。由式(172)及式(182)得
式(184)是关于
一方面,由式(184)、式(185)及式(186)得
另一方面,由式(176)、(184)及式(187)得
故得
由式(185)得(174)及
由式(184)及式(174)得
式(184)至式(190)用于后续方程简化,带有上下标的c是结构常数矩阵。考虑式(182)两边2范数得
考虑当il1=03时的位置矢量对齐关系得
由式(47)及式(193)得
进而,得
即有
显然,有
由式(110)、式(176)及式(195)得式(194)之左式
结构参数
由式(110)、式(185)、式(190)及式(188)得
由式(190)及式(197)得
由式(184)、式(190)及式(198)得
将式(196)及式(199)代入(194),且消去两边的
消去τ4及τ5后的位置方程(169)是3个“3元2阶”多项式方程,可以等同于3r机械臂问题,采用上述基于轴不变量的3r机械臂位置逆解方法进行求解,为实时计算通用6r轴机械臂的逆解奠定了基础。一方面,将有利于提高6r机械臂的绝对定位精度;另一方面,在结构上可以使传统解耦机械臂的第4轴及第5轴向根方向移动,不仅可以优化机械臂的结构,而且有利于提高6r机械臂避让障碍的灵活性。
三、基于轴不变量的通用7r机械臂位姿逆解
称有7个转动轴,拾取点位于第7轴轴线上,且第5轴与第6轴不共轴的机械臂为通用7r机械臂。通用7r机械臂通过前6轴控制第7轴与期望位置及姿态对齐,第7轴无限转动或控制第7轴满足径向对齐。因此,7r通用机械臂比6r机械臂具有更强的空间操作的灵活性。
给定期望位置
若给定7r轴链il7=(i,1:7],il1=03,期望位置矢量及ju-gibbs四元数分别记为
其中:
系统结构参数及期望ju-gibbs姿态四元数构成的矩阵表示为
上式建立步骤为:
期望
由式(52)得式(201),由式(206)得
由式(53)得
其中:
由式(207)及式(208)得
由式(64)得
将式(211)代入式(210)得
其中:
由式(213)得式(203)。由式(203)及式(212)得
式(214)是关于
由式(214)及式(205)得
式(214)至式(217)用于后续方程简化,c是结构常数矩阵。考虑式(216)两边2范数得
由式(62)得
一方面,由式(216)、式(215)及式(219)得
另一方面,由式(206)、(214)及式(220)得
故得
由式(222)得式(200)中的姿态方程。
考虑当il1=03时的位置矢量对齐关系得
由式(47)及式(223)得
进而,得
即有
显然,有
由式(206)及式(225)得式(224)之左式
由式(110)、式(215)、式(216)及式(221)得
由式(216)及式(227)得
由式(214)、式(217)及式(228)得
将式(226)及式(229)代入(224),且消去两边的
消去τ5及τ6后的位姿方程(200),是4个“4元2阶”多项式方程。
1.ju-gibbs增量四元数及性质
通用7r机械臂位姿逆解的计算复杂度较高,需要解决实时计算的技术问题。在工程上,计算精度是一个相对的概念,只需要保证数值计算精度高于系统结构参数精度的4至6倍就能满足工程要求。式(59)所示的姿态方程及式(47)所示的位置方程本质上是ju-gibbs四元数的表达式。只要保证式(59)及式(47)具有足够的计算精度就能满足工程精度要求。下面,首先提出“居-吉布斯”增量四元数(delta-quaternion),再建立通用7r机械臂增量位姿方程,并进行实时逆解计算。
【1】“居-吉布斯”增量四元数定义
定义“居-吉布斯”增量四元数
其中:
显然,“居-吉布斯”增量四元数是四维复数,且有
【2】“居-吉布斯”增量四元数性质
由式(48)及式(230)得
由式(51)得
由式(234)得
由式(50)得
由式(236)及式(54)得
由式(62)得
2.基于轴不变量的通用7r机械臂运动规划
通用7r机械臂逆解由于计算复杂度极高,在现有技术条件下无法实现。但是,通常位于第7轴上的拾取点与第6轴的距离很小。因此,在第6轴取距拾取点较近的点为名义拾取点,先计算通用6r机械臂的逆解;以之为基础,再应用数值迭代法,完成通用7r机械臂的运动规划与逆解计算。下面,探讨通用7r机械臂的增量(delta)位姿方程建立及逆解的问题。
【1】通用7r轴链的增量位姿方程
首先陈述由gibbs增量四元数表征的通用7r轴链的增量位姿方程,最后求其逆解。
若给定6r轴链il7=(i,1:7],il1=03,位置矢量及ju-gibbs增量四元数分别记为
其中:
由行四元数构成的系统结构参数矩阵表示为
下面,分析通用7r机械臂的增量位姿逆解。显然,式(239)是关于{εl|l∈[1:4}的线性方程。将式(239)重新表示为
a·[ε1ε2ε3ε4]t=b(247)
若a-1存在,解式(245)得
[ε1ε2ε3ε4]t=a-1·b。(248)
由式(235)及式(248)得
通用7r机械臂位姿逆解的作用在于:可以通过增量位置矢量
【2】基于偏速度迭代的通用7r机械臂运动规划
当k∈iln时,由式(3.304)至式(3.306)可知
记运动链为
记
由式(252)应用梯度(gradientdescentmethod/gdm)下降法得
其中:步长step>0,step→0。显然有
选择step步长,由初态
(1)确定目标函数
显然,goal表示
(2)选择步长
一方面,应用构造法得确定步长
则由式(253)及式(256)可知:仅当
另一方面,由式(231)及式(252)得
其中:ε(i,6]=[ε1ε2…ε6]t。由式(257)及式(253)确定步长
(3)迭代过程
一方面,若取式(258)之步长step,由式(259)完成迭代计算:
其中:当
δgoal≤0,(260)
即式(259)的迭代过程一定收敛。
迭代计算式的推导过程为:由式(252)及式(255)得
由式(253)、式(256)及式(261)得
因式(256)与式(258)在理论上等价,故可用式(258)替代式(256)。但是式(256)与式(258)的计算过程不同:因计算机字长有限,当
当
基于偏速度的通用7r机械臂运动规划特点在于:通过迭代,逐步逼近期望的位姿,可以得到一条由初始位姿至期望位姿的路径。因为该方法是趋向目标的优化过程,所以实时性较差。若在迭代过程中,控制关节增量,则可以满足关节速度的约束;因此,在完成运动规划的同时,也获得对应期望位姿的一组逆解。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。