一种绳索康复机器人运动学参数标定和解算方法及系统与流程

1.本发明属于康复机器人技术领域，涉及一种基于编码器和imu的绳索康复机器人运动学参数标定和解算方法及系统。


背景技术：

2.由于社会老龄化的加剧，我国脑卒中的发病率不断上升，且大部分脑卒中患者在术后会遗留功能障碍。偏瘫是脑卒中后常见残疾，其主要表现为患侧肌力减退、肢体痉挛以及活动受限，对患者日常生活造成严重影响。
3.研究表明，通过康复训练能够促进大脑损伤区域或冗余神经重新学习，恢复人体运动功能，近年来，机器人行业发展迅猛，使用机器人进行康复训练也愈发普及，相关研究也表明机器人康复训练具有较好的效果且能显著降低治疗成本。
4.上肢康复训练机器人经过多年发展，大体可以分为末端式和外骨骼式两种，末端牵引式机器人，通过与人手臂末端相固定，进而牵引手臂进行运动康复，这类机器人结构通常较为简单且不会引起训练过程中的关节错位问题，但是其工作范围有限且运动灵活性相对较差；外骨骼机器人通过与人上肢多个关节绑定来驱动手臂运动，相对末端机器人具有较好的灵活性和较大的工作空间，但是其结构通常比较复杂，并且，由于需要让机器人关节旋转中心与人体关节旋转中心相对准，增大了设计难度，在肩关节等非中心固定关节处，更是较难满足要求。
5.上述问题在传统的刚性机器人上多有发生。近年来，柔性外骨骼机器人得到了广泛的发展。与刚性机器人相比，柔性机器人使用绳索或人工肌肉模拟人体肌肉收缩，因而对人体的适应性更强，且结构更加简单轻便。同时由于柔性外骨骼机器人常将人体关节作为机器人结构的一部分，因此有效避免了传统外骨骼机器人关节轴对准和关节长度调整的问题。
6.而对于此类柔性外骨骼机器人来说，其柔性的结构通常会导致不确定的运动学参数，如穿戴时患者的肩关节中心位置及设备穿戴位置等运动学参数会相对于预设值发生偏移。同时考虑到人肩关节的运动不能简单理解为一个球状关节的运动，在实际的肩关节运动中，上臂的抬举动作会使得肩关节的旋转中心同步抬升，给机器人的运动学建模和运动控制带来影响。目前国内外的柔性外骨骼上肢康复机器人基本上将肩关节中心视为定点进行参数解算，会给上肢的运动学模型引入误差，因此研究柔性外骨骼机器人的运动学参数的标定和实时解算问题具有实际的意义以及较高的价值。


技术实现要素：

7.为解决上述问题，本发明提供了一种基于编码器和imu的绳索康复机器人运动学参数标定和解算方法及系统。本发明在一种绳索牵引上肢康复机器人的基础上，利用机器人自带的编码器和imu(惯性测量单元)两种传感器实现，首先在使用绳索康复机器人进行康复训练前对患者的运动学参数进行标定，再在使用绳索康复机器人进行康复训练过程中
根据患者的手臂姿态实时解算出当前肩关节中心位置。
8.为了达到上述目的，本发明提供如下技术方案：
9.一种绳索康复机器人运动学参数标定和解算方法，包括如下步骤：
10.步骤一，运动学参数标定
11.(1)首先取患者手臂处于某一姿态下的状态为初始状态，此时的绳长即为初始绳长，亦是待求参数，并记此时手臂的旋转矩阵为即无相对旋转；
12.(2)患者进行手臂各种动作，记录患者运动过程中的编码器数据和imu数据，即可得到绳长变化量δl
i
和上肢旋转矩阵参数，患者手臂在任一位姿下的δl
i
和为一组数据，采集n组数据；
13.(3)对算法迭代量赋初值，即对待求解参数初始化，初始值通过手动测量近似确定，分别赋值为其中，为肩关节中心相对于系统坐标原点的偏移值初始值；为手臂处于某一姿态下的四根绳索初始绳长，其中i＝1,2,3,4,为上臂环中心到肩关节中心的距离；
14.(4)进行迭代，算得迭代参数的修正值，对迭代参数进行修正；
15.(5)将第m次迭代后的迭代参数d
m
与第m
‑
1次的迭代参数d
m
‑1进行比较，当||d
m
‑
d
m
‑1||
∞
＜ε时，停止迭代，待求解参数的解即为当前迭代参数，为＜ε时，停止迭代，待求解参数的解即为当前迭代参数，为否则返回步骤(4)，继续迭代；其中d
m
是第m次迭代后的迭代参数，ε是设定的误差限，m为迭代次数；
16.步骤二，运动学参数解算
17.对n组和进行多元线性回归
18.记标定过程中的数据和的关系系数
[0019][0020][0021]
则bj＝z
[0022]
根据最小二乘法可解得j
[0023]
j＝(b
t
b)
‑1b
t
y
[0024]
在实际训练过程中，通过实时的手臂旋转矩阵可算得患者肩关节中心偏移量
[0025]
进一步的，所述步骤(4)中，具体迭代过程如下：
[0026]
考虑任意其中一组数据，忽略下标j，则在该组数据中，第i根绳索的绳长相对于初始绳长的绳长差用如下方程表示
[0027]
[0028]
而l
i
＝||
o
s
io
u
i
||
ꢀꢀꢀꢀ
(2)
[0029]
将式(1)代入式(2)，可得
[0030][0031]
假定肩关节中心坐标估计值上臂环中心到肩关节中心的距离估计值值为第i根绳索的初始绳长估计值为
[0032]
在待求解参数附近进行泰勒展开，并取至一次项；则
[0033][0034]
将式(4)代入式(3)，展开后并写成
[0035][0036]
假设共采集了n组数据，则根据公式(5)可列出4n组方程，对应的待求解参数为3n组肩关节中心偏移[
o
g
jx
,
o
g
jy
,
o
g
jz
]
t
，4组初始绳长和1组上臂环中心到肩关节中心的距离共3n 5个变量，其中，j＝1,2,...,n，i＝1,2,3,4；
[0037]
因此至少需要采集5组数据方可解出方程，即n≥5：
[0038]
记
[0039]
其中
[0040][0041]
记δl
ij
＝δl
ij
‑
δl
0ij
[0042]
记
[0043]
记
[0044]
则ax＝y
[0045]
根据最小二乘法有
[0046]
x＝(a
t
a)
‑1a
t
y
ꢀꢀꢀꢀꢀ
(6)
[0047]
最后对待求解参数估计值进行迭代
[0048][0049]
上述各式中，
o
g
j
＝[
o
g
jx
,
o
g
jy
,
o
g
jz
]
t
为在坐标系{o}下第j组数据中肩关节中心坐标；表示第j组数据中通过imu测得的手臂的旋转矩阵；为在坐标系{g}下第i根绳索在上臂环上的连接点坐标初始值；
o
s
i
为在坐标系{o}下第i根绳索在肩部支架的连接点坐标；
o
u
ij
为在坐标系{o}下第j组数据中第i根绳索在上臂环上的连接点坐标；l
ij
为第j组数据中第i根绳索的长度，即||
o
s
i o
u
ij
||；δl
ij
为第j组数据中第i根绳索长度相对于初始绳长的绳长差。
[0050]
进一步的，所述步骤二进行解算时，参数视为常量。
[0051]
进一步的，所述步骤一中某一姿态指患者手臂自然下垂时的姿态。
[0052]
进一步的，所述步骤一中各种动作包括手臂外展/内收、前屈/后伸、水平屈/伸动作。
[0053]
本发明还提供了绳索康复机器人运动学参数标定和解算系统，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述计算机程序被加载至处理器时用于实现绳索康复机器人运动学参数标定和解算方法。
[0054]
进一步的，所述计算机程序包括运动学参数标定模块和运动学参数解算模块，运动学参数标定模块用于在使用绳索康复机器人进行康复训练前对患者的运动学参数进行标定；运动学参数解算模块用于在使用绳索康复机器人进行康复训练过程中根据患者的手臂姿态实时解算出当前肩关节中心位置。
[0055]
与现有技术相比，本发明具有如下优点和有益效果：
[0056]
本发明仅使用编码器和imu两类传感器即可有效标定绳索机器人的运动学参数，并在康复训练过程中实时解算相关参数，提高了运动学模型的准确性和控制精度，成本低，安装与操作简单。
附图说明
[0057]
图1为绳索上肢康复机器人示意图，图中，4根绳索(虚线表示)穿过肩部支架后连接到上臂环，上臂环与患者上肢结合形成一个完全约束的机械系统，可控制上肢肩关节3自由度运动。机器人绳索的长度通过编码器测量，上肢的3自由度旋转通过imu测量；
[0058]
图2为绳索上肢康复机器人标定及使用状态示意图，其中(a)为手臂的前屈/后伸动作，(b)为手臂的外展/内收动作，(c)为手臂的水平屈/伸动作。
[0059]
图3为本发明方法流程示意图。
具体实施方式
[0060]
以下将结合具体实施例对本发明提供的技术方案进行详细说明，应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
[0061]
图1为绳索上肢康复机器人示意图，图中，4根绳索(虚线表示)穿过肩部支架后连接到上臂环，上臂环是一个穿戴在手臂上的套筒，上臂环与患者上肢结合形成一个完全约束的机械系统，可控制上肢肩关节3自由度运动。机器人绳索的长度通过编码器测量，上肢的3自由度旋转通过imu测量。图2为绳索上肢康复机器人标定及使用状态示意图。
[0062]
将绳索上肢康复机器人穿戴在患者手臂上，应用本发明提供的基于编码器和imu的绳索康复机器人运动学参数标定和解算方法，如图3所示，本发明方法具体包括如下步骤：
[0063]
(一)运动学参数标定
[0064]
首先统一对本发明需要的参数进行说明如下：
[0065]
下标i:第i根绳索对应参数，i＝1,2,3,4；
[0066]
下标j:第j组数据对应参数，标定过程共采集n组数据，j对应其中某一组数据，j＝1,2,
…
,n；
[0067]
顶标∧：用迭代参数算得的估计量；
[0068]
表示第j组数据中，坐标系{g}相对于坐标系{o}的旋转矩阵，在本发明中主要指代通过imu测得的手臂的旋转矩阵；
[0069]
o
g
j
：
o
g
j
＝[
o
g
jx
,
o
g
jy
,
o
g
jz
]
t
，在坐标系{o}下，第j组数据中，肩关节中心坐标；
[0070]
在坐标系{g}下，第i根绳索在上臂环上的连接点坐标初始值；
[0071]
o
s
i
:在坐标系{o}下，第i根绳索在肩部支架的连接点坐标，相对于坐标系{o}静止；
[0072]
o
u
ij
:在坐标系{o}下，第j组数据中，第i根绳索在上臂环上的连接点坐标，会随着上肢运动相对于坐标系{o}运动，即
[0073]
l
ij
：第j组数据中，第i根绳索的长度，即||
o
s
i o
u
ij
||；
[0074]
第i根绳索的初始长度，i＝1,2,3,4；
[0075]
δl
ij
：第j组数据中，第i根绳索长度相对于初始绳长的绳长差；
[0076]
待求解的参数是g
jx
，g
jy
，g
jz
，其各自含义如下：
[0077]
g
jx
，g
jy
，g
jz
：肩关节中心相对于系统坐标原点的偏移
o
g
j
＝[
o
g
jx
,
o
g
jy
,
o
g
jz
]
t
，(j＝1,2,...,n)。考虑到正常人手臂在抬升过程中，肩关节的旋转中心会同步抬升，因此n组数据对应n组肩关节中心偏移；
[0078]
取手臂处于某一姿态下的四根绳索长度记为初始绳长通常取手臂正常下垂时的绳长值作为初始绳长。
[0079]
上臂环中心到肩关节中心的距离为
[0080]
具体标定过程如下：
[0081]
1.首先取患者手臂自然下垂时的状态为初始状态，此时的绳长即为初始绳长，亦是待求参数，并记此时手臂的旋转矩阵为即无相对旋转。
[0082]
2.如图2所示，患者在医师的帮助下进行手臂的外展/内收，前屈/后伸、水平屈/伸等动作，记录运动过程中的编码器数据和imu数据，即可得到绳长变化量δl
i
和上肢旋转矩阵等参数，患者手臂在任一位姿下的δl
i
和为一组数据，采集n组数据(n≥5)。
[0083]
3.对算法迭代量赋初值，即对待求解参数初始化，初始值通过手动测量近似确定，分别赋值为
[0084]
4.进行迭代，算得迭代参数的修正值，对迭代参数进行修正。
[0085]
具体迭代原理如下：
[0086]
考虑任意其中一组数据，忽略下标j,则在该组数据中，第i根绳索的绳长相对于初始绳长的绳长差可用如下方程表示
[0087][0088]
而l
i
＝||
o
s
i o
u
i
||
ꢀꢀꢀꢀ
(2)
[0089]
将式(1)代入式(2)，可得
[0090][0091]
假定肩关节中心坐标估计值上臂环中心到肩关节中心的距离估计值值为第i根绳索的初始绳长估计值为
[0092]
在待求解参数附近进行泰勒展开，并取至一次项。则
[0093][0094]
将式(4)代入式(3)，展开后并写成
[0095][0096]
假设共采集了n组数据，则根据公式(5)可列出4n组方程，对应的待求解参数为3n组肩关节中心偏移[
o
g
jx
,
o
g
jy
,
o
g
jz
]
t
(j＝1,2,...,n)，4组初始绳长和1组上臂环中心到肩关节中心的距离共3n 5个变量。
[0097]
因此至少需要采集5组数据方可解出方程。
[0098]
记
[0099]
其中
[0100][0101]
记δl
ij
＝δl
ij
‑
δl
0ij
[0102]
记
[0103]
记
[0104]
则ax＝y
[0105]
根据最小二乘法有
[0106]
x＝(a
t
a)
‑1a
t
y
ꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0107]
最后对待求解参数估计值进行迭代
[0108][0109]
5.将第m次迭代后的迭代参数d
m
与第m
‑
1次的迭代参数d
m
‑1进行比较，当||d
m
‑
d
m
‑1||
∞
＜ε时(ε是设定的误差限)，停止迭代，待求解参数的解为＜ε时(ε是设定的误差限)，停止迭代，待求解参数的解为(m为迭代次数)，否则返回步骤4，继续迭代。
[0110]
(二)运动学参数解算
[0111]
根据上述运动学参数标定方法，能够在机器人训练前对患者进行运动学参数标定，获得待求解运动学参数在迭代解算后的近似解(m为迭代次数)。
[0112]
其中参数仅与患者手臂长度，上臂环穿戴位置等因素有关，不随患者手臂在训练过程的运动而变化，因此将其视为常量。而参数表示患者肩关节旋转中心在坐标系{o}下的坐标，会随着患者手臂在训练过程的运动而发生改变，属于需要实时解算的参数。解算过程如下：
[0113]
在上述运动学参数标定过程中，对于n组数据，共有n组手臂旋转矩阵和n组肩关节中心偏移量标定值运动学解算的目标是找到标定过程中的数据和的关系系数j，并在训练过程中，根据关系系数j和实时测得的手臂旋转矩阵实时解算出患者肩关节中心偏移量
o
g。
[0114]
为便于计算，对n组和进行多元线性回归
[0115]
记
[0116][0117][0118]
则bj＝z
[0119]
根据最小二乘法可解得j
[0120]
j＝(b
t
b)
‑1b
t
y
[0121]
在实际训练过程中，通过实时的手臂旋转矩阵可算得患者肩关节中心偏移量
[0122]
基于相同的发明构思，本发明还提供了一种绳索康复机器人运动学参数标定和解算系统，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，该计算机程序被加载至处理器时实现前述的一种绳索康复机器人运动学参数标定和解算方法方法。计算机程序包括运动学参数标定模块和运动学参数解算模块，运动学参数标定模块用于在使用绳索康复机器人进行康复训练前对患者的运动学参数进行标定，具体实现上述步骤一内容；运动学参数解算模块用于在使用绳索康复机器人进行康复训练过程中根据患者的手臂姿态实时解算出当前肩关节中心位置，具体实现上述步骤二内容。
[0123]
本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段，还包括由以上技术特征任意组合所组成的技术方案。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。