一种基于激励轨迹的机器人末端负载快速辨识方法

1.本发明属于机器人控制技术领域，具体是涉及一种基于激励轨迹的机器人末端负载快速辨识方法。


背景技术：

2.为了实现各种柔顺/力控时，在机器人机械臂末端加装6维力传感器是比较典型的做法，其可以感知其后方工具重力以及真实施加的有效负载力，然后作为反馈可以实现各种力控，比如力控笛卡尔空间拖动示教、末端力控等。在机器人接触作业中，依靠安装在末端多维力传感器感知外部的接触力。但当机械臂的有效载荷和速度较大时，其产生的重力/惯性力对力传感器的数据输出将产生不可忽视的影响，从而严重影响力控的精度和效果。
3.针对这些问题，目前有不少的文献提出了重力补偿算法，但目前的大部分算法都只采集了空间有限点进行重力补偿的计算，补偿的精度相对较差且运行的效率也非常低，已经无法满足现有的工业现场的高精度高速率的力控要求，将很大程度地影响力控作业的效率。


技术实现要素：

4.发明目的：本发明目的在于针对现有技术的不足，提供一种基于激励轨迹的机器人末端负载快速辨识方法，设计基于激励轨迹的最优空间辨识位姿，各关节的激励轨迹采用有限项的傅里叶级数表示，采用最小二乘法计算得到负载的辨识参数，解决机器人末端负载辨识效率低和精度差的问题。
5.技术方案：本发明所述基于激励轨迹的机器人末端负载快速辨识方法，包括如下步骤：
6.s10、建立机器人和机器人末端负载的空间坐标系，并根据建立的坐标系关系，分别推导负载辨识计算公式和惯性力影响公式，从而获得重力和惯性力的综合表达式；
7.s20、基于优化目标函数及各限制条件，设计机器人末端各关节激励轨迹，以此改变机器人末端姿态；
8.s30、在采集机器人末端姿态基于激励轨迹改变时，采集至少12组机器人在不同姿态下安装在机器人末端的六维力传感器的数据，最后基于重力和惯性力的综合表达式通过最小二乘法得到机器人末端负载的各项参数。
9.本发明进一步优选地技术方案为，步骤s10获得重力和惯性力的综合表达式的具体方法为：
10.s1001、首先建立参考坐标系，设定世界坐标系为w，机器人基坐标系为b，机器人末端坐标系为e，六维力传感器自身的坐标系为ft；其中六维力传感器坐标系的原点为o，坐标系的方向和传感器厂家给定的坐标系方向一致；
11.然后确定坐标系之间的变换矩阵，通过所建立的坐标系，得到六维力传感器坐标系相对于机器人末端坐标系的齐次变换矩阵为对于任意给定的末端位子末端负载
的重力在ft坐标系中的各个轴上均有分量；
12.s1002、假设六维力传感器的三个力的初始值为和三个力矩的初始值为和负载重力为g以及负载重心在力传感器坐标系中的位置为x,y和z，则负载重力在传感器坐标系下的作用力分量为t
gx
、t
gy
和t
gz
，作用力矩分量为t
gx
、t
gy
和t
gz
，六维力传感器实测的作用力反馈值为f
x
、fy和fg，作用力矩为t
x
、ty和tz；
13.s1003、由力与力矩的关系以及右手定则可得公式(1)：
[0014][0015]
其中：
[0016][0017]
通过以上公式，并对多组数据采样，最终通过最小二乘法可以得到负载质心坐标：
[0018]
p＝(f
t
f)-1
·ft
τ#(3)；
[0019]
其中：p＝[x y z k
1 k
2 k3]
t
[0020][0021]
在机器人基坐标系的z轴和重力的方向不一致，存在安装倾角时，假设机器人基座与世界坐标系的x轴和y轴方向的偏角即安装倾角分别为u和v，且重力在世界坐标系中的方向向量为gw＝[0 0
ꢀ‑
1]
t
，设负载重力为g，得到：
[0022][0023]
s1004、机器人在快速运动中，在负载的作用下，机器人末端产生惯性f
inert
和惯性力矩t
inert
，则传感器的实测值f
meas
与t
meas
由期望跟踪力与力矩f
des
与t
des
以及惯性力与力矩组成，即：
[0024][0025]
s1005、通过公式(4)和公式(5)得到重力和惯性力的综合表达式：
[0026][0027]
优选地，步骤s20中设计机器人末端各关节激励轨迹的具体方法为：
[0028]
s2001、找到与关节参数相关的方程组，将其作为激励方程，即考虑惯性力下的力矩方程组：
[0029][0030]
由于负载的加速度为负载绕末端坐标系旋转产生的加速度与坐标系本身加速度的矢量和，则有：
[0031][0032]
其中α是机器人末端坐标系的角加速度，是末端坐标系的加速度；得到：
[0033]
ta＝a
·
p#(9)；
[0034]
其中：
[0035]
a是关于θ4、θ5、θ6的函数，记为a＝f(q4,q5,q6)；
[0036]
s2002、各个关节的激励轨迹用有限项的傅里叶级数表示，则机器人第i关节的角位移qi、角速度以及角加速度表示为：
[0037][0038][0039][0040]
式中的n表示正弦和余弦的项数，wf＝2πff是基波角频率，a
i,l
和b
i,l
是待确定的参数；
[0041]
为使机器人末端尽可能的充满整个空间，设优化的目标函数为minf(a)；由于重力补偿只涉及机械臂末端的末三轴，因此共有6n个参数待优化；为满足机器人关节限制条件，给出如下关节约束表达式：
[0042][0043]
启停约束：
[0044][0045]
s2003、考虑辨识取点姿态“不具有代表性”及传感器数据浮动等因素，ma、a和p均受到扰动，则扰动后的方程组为：
[0046]
ta+δta＝(a+δa)
·
(p+δp)#(15)；
[0047]
由此得到：
[0048]
δp＝a-(δt
a-δa
·
p-δa
·
δp)#(16)；
[0049]
由||δp||≤||a-||(||δta||+||δa||
·
||p||+||δa||
·
||δp||)，解出：
[0050][0051]
由于ta＝a
·
p，则||ta||≤||a||
·
||p||，因此式(16)写成：
[0052][0053]
则辨识参数b的相对误差为：
[0054][0055]
由于是相对扰动项，辨识参数矩阵p的误差主要由||a-||
·
||a||引起，在此引入条件数为激励轨迹的优化目标函数：
[0056]
cond(a)＝||a-||
·
||a||#(20)；
[0057]
因此得到：
[0058][0059]
则优化目标函数为：
[0060]
minf(a)＝||a-||
·
||a||#(22)；
[0061]
在满足机器人各关节位置、速度、加速度条件下，当优化目标函数minf(a)值最小时即表示外界扰动对辨识参数矩阵p的扰动最小，对应的参数a
i,l
、b
i,l
即激励轨迹对应参数；先假定工件质心在机器人末端坐标系中的位置矢量为一已知量p，并将其代入式(9)，从而激励出具有代表性的机器人位姿。
[0062]
优选地，步骤s30获得机器人末端负载的各项参数的具体方法为：
[0063]
通过重力和惯性力的综合表达式，可以得到如下公式：
[0064]
f＝rl#(23)；
[0065]
其中：
[0066][0067]
通过最小二乘法可以得到：
[0068]
l＝(r
t
r)-1
·rt
f#(24)；
[0069]
h＝(r
t
r)-1
·rt
t#(25)；
[0070]
方程(25)中：
[0071][0072][0073][0074]
[0075][0076]
采集至少12组机器人在不同姿态下安装在机器人末端的六维力传感器的数据，带入公式(24)和公式(25)，求出负载重力g，基座倾角u和v，传感器零点值入公式(24)和公式(25)，求出负载重力g，基座倾角u和v，传感器零点值和负载质心x、y和z共12个参数的数值。
[0077]
有益效果：本发明通过设计一条尽可能代表空间预期姿态的激励轨迹，并采用有限项的傅里叶级数表示各个关节的激励轨迹，然后通过最小二乘法拟合n组机器人末端不同姿态下的数据并结合参数化模型对机器人末端负载进行快速辨识。本发明可快速辨识机器人末端的负载参数，弥补了辨识效率低和精度差的问题，使之能适应当前高速与高精度的力控工业场景中。
附图说明
[0078]
图1为本发明实施例中坐标系和负载重力示意图；
[0079]
图2为本发明实施例中末三轴的关节角度运动轨迹；
[0080]
图3为本发明实施例中末端负载位置曲线；
[0081]
图4为本发明实施例中末端负载速度曲线；
[0082]
图5为本发明实施例中末端负载加速度曲线。
具体实施方式
[0083]
下面通过附图对本发明技术方案进行详细说明，但是本发明的保护范围不局限于所述实施例。
[0084]
实施例：一种基于激励轨迹的机器人末端负载快速辨识方法，包括如下步骤：
[0085]
s10、选定一款六维力传感器安装在机器人末端，建立机器人和机器人末端负载的空间坐标系，并根据建立的坐标系关系，分别推导负载辨识计算公式和惯性力影响公式，从而获得重力和惯性力的综合表达式；
[0086]
s1001、首先建立参考坐标系，设定世界坐标系为w，机器人基坐标系为b，机器人末端坐标系为e，六维力传感器自身的坐标系为ft；其中六维力传感器坐标系的原点为o，坐标
系的方向和传感器厂家给定的坐标系方向一致；
[0087]
然后确定坐标系之间的变换矩阵，通过所建立的坐标系，得到六维力传感器坐标系相对于机器人末端坐标系的齐次变换矩阵为对于任意给定的末端位子末端负载的重力在ft坐标系中的各个轴上均有分量；
[0088]
s1002、假设六维力传感器的三个力的初始值为和三个力矩的初始值为和负载重力为g以及负载重心在力传感器坐标系中的位置为x,y和z，则负载重力在传感器坐标系下的作用力分量为f
gx
、f
gy
和f
gz
，作用力矩分量为t
gx
、t
gy
和t
gz
，六维力传感器实测的作用力反馈值为f
x
、fy和fg，作用力矩为t
x
、ty和tz；
[0089]
s1003、由力与力矩的关系以及右手定则可得公式(1)：
[0090][0091]
其中：
[0092][0093]
通过以上公式，并对多组数据采样，最终通过最小二乘法可以得到负载质心坐标：
[0094]
p＝(f
t
f)-1
·ft
τ#(3)；
[0095]
其中：p＝[x y z k
1 k
2 k3]
t
[0096][0097]
在机器人基坐标系的z轴和重力的方向不一致，存在安装倾角时，假设机器人基座与世界坐标系的x轴和y轴方向的偏角即安装倾角分别为u和v，且重力在世界坐标系中的方向向量为gw＝[0 0
ꢀ‑
1]
t
，设负载重力为g，得到：
[0098][0099]
s1004、机器人在快速运动中，在负载的作用下，机器人末端产生惯性f
inert
和惯性力矩t
inert
，则传感器的实测值f
meas
与t
meas
由期望跟踪力与力矩f
des
与t
des
以及惯性力与力矩组成，即：
[0100][0101]
s1005、通过公式(4)和公式(5)得到重力和惯性力的综合表达式：
[0102][0103]
s20、重力和惯性力的各个未知参数通过改变末端的姿态可辨识得到(即与末三轴有关)，为使得选取的姿态能尽可能代表辨识空间，需要设计一条最优的关节参考轨迹。将各关节的激励轨迹用有限项的傅里叶级数表示，设计一条最优关节参考轨迹，以此改变末端姿态辨识得到重力的各个未知参数。将傅里叶级数中的正余弦项数取为n＝5，基频f_f选定为0.1hz，采样频率取100hz，总共1000次采样来计算矩阵a。
[0104]
s2001、找到与关节参数相关的方程组，将其作为激励方程，即考虑惯性力下的力矩方程组：
[0105][0106]
由于负载的加速度为负载绕末端坐标系旋转产生的加速度与坐标系本身加速度的矢量和，则有：
[0107][0108]
其中α是机器人末端坐标系的角加速度，是末端坐标系的加速度；得到：
[0109]
ta＝a
·
p#(9)；
[0110]
其中：
[0111]
a是关于θ4、θ5、θ6的函数，记为a＝f(q4,q5,q6)；
[0112]
s2002、各个关节的激励轨迹用有限项的傅里叶级数表示，则机器人第i关节的角位移qi、角速度以及角加速度表示为：
[0113]
[0114][0115][0116]
式中的n表示正弦和余弦的项数，wf＝2πff是基波角频率，a
i,l
和b
i,l
是待确定的参数；
[0117]
为使机器人末端尽可能的充满整个空间，设优化的目标函数为minf(a)；由于重力补偿只涉及机械臂末端的末三轴，因此共有6n个参数待优化；为满足机器人关节限制条件，给出如下关节约束表达式：
[0118][0119]
启停约束：
[0120][0121]
s2003、考虑辨识取点姿态“不具有代表性”及传感器数据浮动等因素，ma、a和p均受到扰动，则扰动后的方程组为：
[0122]
ta+δta＝(a+δa)
·
(p+δp)#(15)；
[0123]
由此得到：
[0124]
δp＝a-(δt
a-δa
·
p-δa
·
δp)#(16)；
[0125]
由||δp||≤||a-|(||δta||+||δa||
·
||p||+||δa||
·
||δp||)，解出：
[0126][0127]
由于ta＝a
·
p，则||ta||≤||a||
·
||p||，因此式(16)写成：
[0128]
[0129]
则辨识参数b的相对误差为：
[0130][0131]
由于是相对扰动项，辨识参数矩阵p的误差主要由||a-||
·
||a||引起，在此引入条件数为激励轨迹的优化目标函数：
[0132]
cond(a)＝||a-||
·
||a||#(20)；
[0133]
因此得到：
[0134][0135]
则优化目标函数为：
[0136]
minf(a)＝||a-||
·
||a||#(22)；
[0137]
在满足机器人各关节位置、速度、加速度条件下，当优化目标函数minf(a)值最小时即表示外界扰动对辨识参数矩阵p的扰动最小，对应的参数a
i,l
、b
i,l
即激励轨迹对应参数；先假定工件质心在机器人末端坐标系中的位置矢量为一已知量p，并将其代入式(9)，从而激励出具有代表性的机器人位姿。
[0138]
s30、在采集机器人末端姿态基于激励轨迹改变时，采集12组机器人在不同姿态下安装在机器人末端的六维力传感器的数据，最后基于重力和惯性力的综合表达式通过最小二乘法得到机器人末端负载的各项参数。
[0139]
通过重力和惯性力的综合表达式，可以得到如下公式：
[0140]
f＝rl#(23)；
[0141]
其中：
[0142][0143]
通过最小二乘法可以得到：
[0144]
l＝(r
t
r)-1
·rt
f#(24)；
[0145]
h＝(r
t
r)-1
·rt
t#(25)；
[0146]
方程(25)中：
[0147][0148]
[0149][0150][0151][0152]
采集至少12组机器人在不同姿态下安装在机器人末端的六维力传感器的数据，带入公式(24)和公式(25)，求出负载重力g，基座倾角u和v，传感器零点值入公式(24)和公式(25)，求出负载重力g，基座倾角u和v，传感器零点值和负载质心x、y和z共12个参数的数值。
[0153]
如上所述，尽管参照特定的优选实施例已经表示和表述了本发明，但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下，可对其在形式上和细节上作出各种变化。