-
摘要: 针对康复机器人运动过程中的人机交互性问题，提出一种下肢康复机器人自适应人机交互控制策略.提取伸屈运动中下肢表面肌电信号（Surface electromyography，sEMG）和足底压力特征，分别用于表征下肢运动意图和人机交互力（Interaction force，IF）信息，建立基于sEMG-IF的人机交互信息融合模型，实现下肢康复机器人运动轨迹的在线规划；考虑主动康复运动过程中的人机交互作用，建立具有时变动态特性的人机系统动力学模型，设计间接模糊自适应控制器对期望轨迹进行跟踪控制，实现下肢康复机器人自适应人机交互控制.通过对5名被试者进行下肢康复机器人运动控制实验研究，验证所提方法的可行性和有效性.Abstract: Aiming at the problem of human-machine interaction in rehabilitation robot's movement, we propose an adaptive control strategy for lower limb rehabilitation robots. During flexion and extension, the surface electromyography (sEMG) of lower limbs and plantar pressure features are extracted respectively to represent lower limbs' motion intention and interaction force (IF). An sEMG-IF based human-machine interaction and information fusion model is established to program the motion trails of the rehabilitation robot online. Considering the human-machine interaction in active rehabilitation, a man-machine system dynamic model with time-varying dynamic characteristics is established. An indirect fuzzy adaptive controller is designed to trace and control the desired trajectory, and achieve adaptive interactive control of the lower limb rehabilitation robot. Validity and feasibility of the proposed strategy are verified by analysis of the data from 5 subjects under limb movement with the rehabilitation robot.
-
由脑卒中、外伤等导致的运动功能障碍康复问题, 已经在世界范围内引起广泛关注.康复机器人由于能够协助或替代医师进行患者康复训练, 近年来得到迅速发展并在临床康复中得到了广泛应用[1-2].但现有的康复机器人还存在控制策略不够智能, 自适应调整能力和人机交互性水平不高等问题[3].因此, 研究有效的下肢康复机器人自适应人机交互控制方法, 体现患者的运动意图和人机交互的同时自适应控制康复机器人运动, 成为康复机器人研究的关键问题.
康复机器人的人机交互可以通过人体生理信号和力/位置信息两种途径体现[4].生理电信号中的表面肌电信号(Surface electromyography, sEMG)能够反映肌肉状态和人体运动意图[5].有研究将sEMG引入康复机器人, 如Pittaccio提出了一种基于sEMG的踝关节康复机器人[6], 利用胫骨前肌sEMG触发康复机器人按照预定义轨迹运动, 但未考虑人机交互问题; Song等开发基于sEMG连续控制的康复机器人[7], 基于sEMG幅值按照比例控制关节扭矩, 具有了简单程度的人机交互.基于sEMG的康复机器人控制一定程度上体现了人体的运动意图, 但sEMG易受干扰, 采集精度无法保证[4].利用力/位置传感器可以获取人机之间的作用力和位置信息, 并实现康复机器人的阻抗控制和人机交互[8].如Saglia等利用阻抗控制调整脚踝康复机器人的柔顺性实现患者的辅助训练[9]; Duschauwicke等则利用阻抗控制调节患者的步行速度和牵引力[10]; Koopman等通过阻抗控制实现"patient-cooperative"和"human-centred"控制以提高患者的积极参与度[11].上述康复机器人阻抗控制方法具有一定人机交互能力, 但其阻抗控制参数固定不变, 应用于具有个体差异性的康复训练中略有不足.同时, 人机之间的交互作用使得康复机器人具有时变不确定性[3], 导致人机系统动力学模型无法精确建立, 造成传统的基于模型的控制方法无法实现有效控制.
为此, 本文考虑人体主动运动的因素, 基于拉格朗日法建立下肢康复机器人动力学模型, 提取人体下肢伸屈运动中相关拮抗肌的sEMG信号识别人体运动意图; 引入力反馈构建sEMG与人机交互力(sEMG-IF (Interaction force))信息融合模型, 在线规划下肢康复机器人运动轨迹, 并设计间接自适应模糊控制器, 实现下肢康复机器人运动轨迹跟踪控制.
1. 基于sEMG和IF的下肢康复机器人控制
1.1 人机系统动力学建模
以膝关节下肢康复机器人为研究对象, 针对膝关节康复, 实现被试者下肢水平伸屈运动训练.康复机器人电机驱动安装在髋关节处.
由于下肢康复机器人运动过程中存在人机作用力, 需将二者看作一个整体进行动力学分析[4].其人机系统可简化为曲柄滑块机构, 如图 1所示.
图 1中, $\theta $为$X$轴与大腿连杆的转角(髋关节角度), 为$X$轴与小腿连杆的转角, $M_{i}$为第$i$杆的质量, 即人体下肢第$i$段质量$m_{i}$与康复机器人第$i$杆质量$m'_{i}$之和, $L_{i}$为第$i$杆的长度, $d_{i}$为第$i$杆质心与关节的距离, $I_{i}$为人体下肢第$i$段的转动惯量.由几何法得到与$\theta $的关系为
$ \begin{equation} \label{eq1}\varphi =\arcsin \bigg(\frac{L_1 }{L_2 }\sin \theta \bigg) \end{equation} $
(1) 利用拉格朗日法进行人机动力学建模:
$ \begin{equation} \label{eq2}\tau =D(\theta )\ddot {\theta }+H(\theta, \dot {\theta })\dot {\theta }+G(\theta ) \end{equation} $
(2) 式中, $\tau $为系统驱动力, $D(\theta )$为系统的转动惯量, 为科氏力及离心相关项, $G(\theta )$为重力项.
$ \begin{align} \label{eq3} D(\theta )=\, &I_1 \theta +d_1^2 M_1 \theta +[I_2 +M_2 (L_2-d_2 )^2]\times\nonumber\\ &\frac{L_1^2 \cos ^2\theta }{c^2}+M_2 L_1^2 \sin ^2\theta +\nonumber\\ & M_2 \bigg(\frac{2d_2 }{L_2 }-1\bigg)L_1^4 \frac{\cos ^2\theta \sin ^2\theta }{c^2}+\nonumber\\ & M_2 \frac{d_2 }{L_2 }L_1^3 \frac{\cos \theta \sin ^2\theta }{c} +M_3 L_1^2 \sin ^2\theta +\nonumber\\ &2M_3 L_1^3 \frac{\sin ^2\theta \cos \theta }{c}+M_3 L_1^4 \frac{\sin ^2\theta \cos ^2\theta }{c^2} \end{align} $
(3) 式中, $c=\sqrt {L_2^2 -L_1^2 \sin ^2\theta } $.
$ \begin{align} \label{eq4} H(\theta, \dot {\theta })=\, &\bigg[\frac{1}{2}I_2 +\frac{1}{2}M_2 d_2 ^2L_2^2 \bigg]\times\nonumber\\ &\frac{2L_1^2 \cos \theta \sin \theta c^2+2L_1^4 \cos ^3\theta \sin \theta }{c^4}\dot {\theta } +\nonumber\\ &M_2 L_1^2 \sin \theta \cos \theta \cdot \dot {\theta }+\frac{1}{2}M_2 \frac{d_2 }{L_2 }L_1^3 \times\nonumber\\ &\bigg(\frac{3\cos ^2\theta \sin \theta -\sin \theta }{c}+\nonumber\\ & \frac{L_1^2 \cos ^2\theta \sin ^3\theta }{c^3}\bigg)\dot {\theta }+ M_2 \bigg(\frac{2d_2 }{L_2 }-1\bigg)\times\nonumber\\ & L_1^4 \bigg(\frac{\sin \theta \cos \theta -2\sin ^3\theta \cos \theta }{c^2}+\nonumber\\ &\frac{L_1^2 \cos ^2\theta \sin ^3\theta }{c^4}\bigg)\dot {\theta }+ M_3 L_1^2 \sin \theta \cos \theta +\nonumber\\ &M_3 L_1^3 \bigg(\frac{3\cos ^2\theta \sin \theta -\sin \theta }{c}+\nonumber\\ &\frac{L_1^2 \cos ^2\theta \sin ^3\theta }{c^3}\bigg)\dot {\theta }+\nonumber\\ & M_3 L_1^4 \bigg(\frac{\sin \theta \cos \theta -2\sin ^3\theta \cos \theta }{c^2}+\nonumber\\ &\frac{L_1^2 \cos ^2\theta \sin ^3\theta }{c^4}\bigg)\dot {\theta } \end{align} $
(4) $ \begin{equation} \label{eq5} G(\theta )=\left[\frac{d_1 }{L_1 }M_1 +M_2 \left(1-\frac{d_2 }{L_2 }\right)\right]gL_1 \cos \theta \end{equation} $
(5) 式中, $g$为重力常数.
如上文所述, 本文将人体下肢和康复机构看作一个整体进行动力学分析.当患者进行被动康复训练时, 人机系统第$i$个杆件的质量$M_i$为常数, 即人体下肢第$i$段质量与康复机器人第$i$杆质量之和.而当患者进行主动训练时, 康复机器人杆件会受到人体运动及人机交互力的作用, 主动力会在杆件上形成等效质量$\hat{M_i}(t)$, 且施加的力的大小不同, 等效质量也不同, 如式(6)所示:
$ \begin{equation} \label{eq6} M_i'(t)=M_i+\hat{M_i}(t) \end{equation} $
(6) 式中, $M_i'(t)$为患者主动训练时, 第$i$杆件上的质量.因此, 在人体主动运动时, 式(3)~(5)中$M_1$、$M_2$变为随时间变化的$M_1'(t)$和$M_2'(t)$.同时, 人机系统第$i$段杆件质心的位置也会因人体的主动运动而发生变化, 从而导致$d_i$的值不再为固定的常数.因此, 人机系统动力学系数具有时变特性, 其数值会随着人体下肢的主动运动而变化, 分别用、${H}'(\theta )$、${G}'(\theta )$表示.此时式(2)可表示为
$ \begin{equation} \label{eq7} \tau ={D}'(\theta )\ddot {\theta }+{H}'(\theta, \dot {\theta })\dot {\theta }+{G}'(\theta ) \end{equation} $
(7) 即人机交互时的动力学方程.
1.2 基于间接自适应模糊控制器的轨迹跟踪
由式(7)可见, 人机系统动力学模型具有时变不确定性, 从而使得基于模型的控制方法无法有效跟踪期望轨迹.因此, 本文设计模糊自适应控制方法根据人机系统运动状态自适应调整控制参数, 实现轨迹跟踪.
对式(7)进行变换可得:
$ \begin{equation} \label{eq8} \ddot {\theta }=f(\theta, \dot {\theta })+g(\theta , \dot {\theta })u \end{equation} $
(8) 式中, $\theta $为髋关节角度, 和均为非线性函数, $u$为输入驱动力$\tau $.
为表征变量$f({\pmb{\theta}} )$和的时变不确定性, 利用能够反映人机系统状态的模糊变量$\tilde{f}({\pmb{\theta}})$和$\tilde {g}({\pmb{\theta}})$进行代换, $\tilde{f}({\pmb{\theta}}\vert{\pmb{\alpha _f}})$模糊系统的构造步骤如下:
1) 对变量$\theta _i~ (i=1, 2, \cdots, n)$, 定义个模糊集合$A_i^{l_i }~ (l_i =1, 2, \cdots, p_i )$;
2) 采用$\prod_{i=1}^n {p_i } $条模糊规则来构造模糊系统:
If $\theta _1 $ is $A_1^{l_1 }$ and $\cdots$ and $\theta _n$ is $A_n^{l_n }$, then $\hat{f}$ is $E^{l_1 \cdots \ l_n }$.
其中, ${\pmb{\alpha _f}} $为集合, 为下肢康复机器人期望轨迹$\theta _m $与跟踪轨迹之间的误差$e=\theta _m -\theta $, $\theta _2 $为误差变化量$ec$, 模糊集合$A_i^{l_i } $定义为{NB, NS, O, PS, PB}, 即$p_i =5$, $\prod_{i=1}^n {p_i } =25$.
运用乘积推理机、单值模糊器和中心平均解模糊器可得模糊变量为
$ \begin{equation} \label{eq9} \tilde {f}({\pmb{\theta}} \vert {\pmb{\alpha _f}} )=\frac{\sum\limits_{l_1 =1}^{p_1 } \cdots \sum\limits_{l_n =1}^{p_n } {\bar {y}_f^{l_1 \cdots l_n } } (\prod\limits_{i=1}^n {\mu A_i^{l_i } (\theta _i )} )}{\sum\limits_{l_1 =1}^{p_1 } \cdots \sum\limits_{l_n =1}^{p_n } {(\prod\limits_{i=1}^n {\mu A_i^{l_i } (\theta _i )} )} } \end{equation} $
(9) 式中, $\mu A_i^j (\theta _i )$为隶属度函数, 为自由参数.将其放入集合中, 同时引入$\prod_{i=1}^n {p_i } $维向量$\xi ({\pmb{\theta}} )$, 则式(9)变为
$ \begin{equation} \label{eq10} \tilde {f}({\pmb{\theta}} \vert {\pmb{\alpha _g}} )={\pmb\alpha} _{\pmb f}^{\rm T} \xi ({\pmb{\theta}} ) \end{equation} $
(10) 其第$l_1 \cdots l_n $个元素为
$ \begin{equation} \label{eq11} \xi _{l_1 \cdots l_n } ({\pmb{\theta}} )=\frac{\prod\limits_{i=1}^n {\mu A_i^{l_i } (\theta _i )} }{\sum\limits_{l_1 =1}^{p_1 } \cdots \sum\limits_{l_n =1}^{p_n } {(\prod\limits_{i=1}^n {\mu A_i^{l_i } (\theta _i )} )} } \end{equation} $
(11) 同理可得, .进一步, 设计自适应律, 用于调整集合参数和${\pmb\alpha}_{\pmb g}^{\rm T} $:
$ \begin{equation} \label{eq12} {\pmb{\dot {\alpha }_f}} =-\gamma _1 {\pmb{e}}^{\rm T}{U}{\pmb{b}}\xi ({\pmb{\theta}} ) \end{equation} $
(12) $ \begin{equation} \label{eq13} {\pmb{\dot {\alpha }_g}} =-\gamma _2 {\pmb{e}}^{\rm T} {U}{\pmb{b}}\eta ({\pmb{\theta}})u \end{equation} $
(13) 式中, $\gamma _1 $、$\gamma _2 $为正常数.本文中取、$\gamma _2 =0.5$、${\pmb{b}}=[0, 1]^{\rm T}$, $U$为满足Lyapunov方程的正定矩阵:
$ \begin{equation} \label{eq14} {\Lambda}^{\rm T} {U}+{U}{\Lambda} =-{Q} \end{equation} $
(14) 式中, , ${Q}$为任意$n\times n$矩阵.本文中取$k_1 =2$, $k_2 =1$, .
以下肢康复机器人运动轨迹误差${\pmb{e}}$为控制目标, 令
$ \begin{equation} \label{eq15} {\pmb{e}}=(e, \dot {e}, \cdots e^{(n-1)})^{\rm T} \end{equation} $
(15) 取控制律为
$ \begin{equation} \label{eq16} u^\ast =\frac{1}{g({\pmb{\theta}})}(-f({\pmb{\theta}} )+\ddot {\theta}_m +{ {\pmb K}}^{\rm T}{\pmb{e}}) \end{equation} $
(16) 式中, ${\pmb K}=(k_2, k_1 )^{\rm T}$为使得多项式所有根都位于复平面左半平面内的向量.将式(16)带入式(8), 得到系统的闭环控制方程为
$ \begin{equation} \label{eq17} \ddot {e}+k_1 \dot {e}+k_2 e=0 \end{equation} $
(17) 经模糊化后得到模糊变量$\tilde {f}({\pmb{\theta}} )$和, 则控制律变为
$ \begin{equation} \label{eq18} u=\frac{1}{\tilde {g}({\pmb{\theta}} \vert {\pmb{\alpha _g}} )}(-\tilde {f}({\pmb{\theta}} \vert {\pmb{\alpha _f}} )+\ddot {\theta}_m +{\pmb K}^{\rm T}{\pmb{e}}) \end{equation} $
(18) 将式(18)代入式(8), 得到间接自适应控制器模型:
$ \begin{equation} \label{eq19} \ddot {e}=-{\pmb K}^{\rm T}{\pmb{e}}+[\tilde {f}({\pmb{\theta}} \vert {\pmb{\alpha _f}} )-f({\pmb{\theta}} )]+[\tilde {g}({\pmb{\theta}} \vert {\pmb{\alpha _g}} )-g({\pmb{\theta}} )]u \end{equation} $
(19) 通过选取${\pmb K}$值, 使$t\to \infty $时, $e(t)\to 0$, 即系统的输出$\theta $渐近收敛于期望轨迹$\theta _m $.间接自适应控制器的稳定性利用Lyapunov定理[12]得到了证明.
1.3 基于sEMG-IF信息融合的运动轨迹在线规划
为在下肢康复机器人运动过程中体现人体运动意图和人机交互性, 提高患者的主动参与意识和积极性, 本文构建基于表面肌电和人机交互力(sEMG-IF)的信息融合模型实现运动轨迹在线规划.
1.3.1 人体运动意图识别
在下肢伸屈运动过程中, 拮抗肌(股二头肌和股外侧肌) sEMG特征能够体现人体运动意图.但由于sEMG易受噪声干扰, 需要进行信号预处理[13], 具体过程如图 2.
首先, 利用5~200 Hz带通滤波器对原始信号进行滤波; 其次, 采用自适应高通滤波去除sEMG中的基线漂移; 再次, 运用自适应50 Hz陷波滤除sEMG中的工频及谐波干扰; 最后, 对sEMG信号进行全波整流, 得到预处理后的sEMG信号.
进一步, 运用文献[14]中方法获取人体下肢运动意图.将预处理后的股二头肌和股外侧肌sEMG进行归一化处理.
$ \begin{equation} \label{eq20} \bar {x}(t)=\frac{\left| {x(t)-x_r } \right|} {x_m -x_r } \end{equation} $
(20) 式中, $x_r $为静息状态下的sEMG值, 为伸屈运动中最大自主收缩下的sEMG值, $x(t)$为预处理后的sEMG值, $\bar {x}(t)$为归一化后的sEMG值.
由于人体下肢伸屈运动主要受拮抗肌的协同收缩控制, 计算归一化后sEMG的包络, 并取其包络差值识别下肢伸屈运动特征.
$ \begin{equation} \label{eq21} \tilde {x}(t)=\tilde {x}_e (t)-\tilde {x}_f (t) \end{equation} $
(21) 式中, $\tilde {x}_e (t)$为股二头肌包络, 为股外侧肌包络. $\tilde {x}(t)>0$表征下肢伸展运动, 表征下肢屈曲运动.
1.3.2 sEMG-IF信息融合模型
在基于肌电信号定性辨识出人体运动意图基础上, 为进一步提高下肢康复机器人的人机交互能力, 利用足底压力传感器获取人机交互力信息, 建立sEMG-IF信息融合模型实现运动意图量化推理, 进而在线规划下肢康复机器人运动轨迹.
综合考虑被试者运动过程中的参与程度、运动状态的差异以及干扰等因素, 利用人机交互力前后采样时刻的差值刻画被试者的运动意图.当被试者处于下肢伸展运动状态($P=1$)时, 若前后采样时刻人机交互力差值大于阈值$\varepsilon $, 则被试者具有加速运动意图; 若其差值小于阈值$-{\varepsilon }$, 则被试者具有减速运动意图; 否则, 被试者维持当前运动速度.由此, 可以构造出如下模型:
$ \begin{equation} \label{eq22} \begin{cases} P({F}'_R -F_R )>\varepsilon, &\Delta x={1 } \\ \left| {P({F}'_R -F_R )} \right|\le \varepsilon, &\Delta x={0 } \\ P({F}'_R -F_R )<-\varepsilon, &\Delta x={-1} \\ \end{cases} \end{equation} $
(22) 式中, ${F}'_R $为当前时刻的人机交互力, 为前一采样时刻的人机交互力, $P$为第1.3.1节中人体运动意图识别结果($P=1$表示伸展运动, $P=-1$表示屈曲运动), $\varepsilon$为设定的下肢伸屈运动速度改变阈值, 为加减速运动标志量($\Delta x=1$表示加速运动, 表示减速运动).
进一步, 通过下肢康复机器人末端运动速度的变化, 定量描述人体运动意图的改变, 建立如下模型:
$ \begin{equation} \label{eq23} \left| {{F}'_R -F_R } \right|\Delta x=K\frac{{\rm d}v}{{\rm d}t} \end{equation} $
(23) 式中, $K$为根据被试者个体差异性进行调整的下肢康复机器人助力系数, $v$为下肢康复机器人末端沿$X$轴的运动速度.由人机交互力的变化量(和标志量$\Delta x$的乘积体现人体下肢伸屈运动意图的变化.当$\left| {{F}'_R -F_R } \right|$一定时, $K$值越大, 下肢康复机器人运动速度改变越小, 即下肢康复机器人提供的驱动力越小.
为将下肢康复机器人运动速度$v$转换为驱动髋关节角度, 基于下肢康复机器人运动学分析, 得到如下关系式:
$\begin{equation} \label{eq24} x_B =L_1 \cos \theta +L_2 \cos \varphi \end{equation} $
(24) 式中, $x_{B}$为下肢康复机器人机构末端位置, 、$\varphi $为时间变量.对上式求导可得下肢康复机器人末端运动速度:
$ \begin{equation} \label{eq25} v=-L_1 \dot {\theta }\sin \theta -L_2 \dot {\varphi }\sin \varphi \end{equation} $
(25) 将式(1)、式(25)代入式(23)即可得到髋关节角度的期望轨迹$\theta _m $.
2. 实验结果分析
2.1 研究对象与实验流程
采集5名健康被试者(4名男生、1名女生, 年龄: 25±2), 要求被试者实验前无剧烈运动、精神状态良好, 且熟悉实验流程.利用16通道无线表面肌电采集仪, 采样率设为1 kHz, 带通频率设为5 Hz~200 Hz, 输入方式采用三点式差动, AgCl主动电极贴于被试者下肢股二头肌和股外侧肌肌腹处, 参考电极贴于脚踝处; 下肢康复机器人驱动器选用ZM-6405E, 通过FSR402薄膜压力传感器采集足底压力信号.实验过程如下:被试者平躺于下肢康复机器人平台上, 踝关节固定于踏板, 采集被试者下肢屈曲/伸展运动的sEMG信号, 每20次伸展/屈曲动作为一组, 每名被试者采集3组, 每组动作后休息10分钟, 并利用第1节中方法进行下肢康复机器人自适应人机交互控制, 如图 3所示.
2.2 下肢康复机器人运动轨迹在线规划及自适应控制
以其中一名被试者为例, 采集下肢伸屈运动过程中股二头肌和股外侧sEMG, 利用第1.3.1节中方法对原始sEMG进行预处理, 如图 4、图 5所示.由图可见, 被试者进行了3次完整的下肢伸屈运动, 且预处理后的sEMG基线漂移得到了明显抑制.
如图 6所示, 对归一化后的sEMG进一步提取包络值. 图 7为利用第1.3.1节中方法计算股二头肌和股外侧肌的包络差值.由图 7可见, 当包络差值大于零时, 被试者进行伸展动作; 当包络差值小于零时, 被试者进行屈曲动作.
图 7 股二头肌与股外侧肌包络差值Fig. 7 The envelope difference of biceps and lateral muscles of femoral head将上述分析过程应用于下肢康复机器人自适应人机交互控制中的人体运动意图在线识别.由表 1可见, 由于存在个体差异性, 使得5名被试者伸展/屈曲运动完成时间不同, 造成被试者运动意图识别率不同, 其识别正确率均在96%左右, 表明利用股二头肌和股外侧肌sEMG包络差值, 能够用于被试者运动意图在线识别.
表 1 人体运动意图识别结果Table 1 Results of human motion intent recognition被试者 识别总数(个) 识别正确数(个) 识别率(%) A 348 332 95.40 B 315 301 95.55 C 296 285 96.28 D 308 291 94.48 E 286 277 96.85 图 8为伸展/屈曲运动中人机交互力(IF)和运动意图($P$)变化曲线, 图 9为利用本文提出的sEMG-IF信息融合模型得到的下肢康复机器人运动速度变化曲线.由图 8、图 9可见, 下肢康复机器人运动速度与人机作用力变化趋势基本一致:在第一个伸展阶段($P=1$), 下肢康复机器人运动速度随着IF差值的变化而相应的加速和减速; 在屈曲阶段($P=-1$), 下肢康复机器人运动速度随着IF差值的变化出现反向加速和减速.
图 10为基于运动速度并利用第1.3.2节中方法得到的下肢康复机器人期望轨迹, 以及运用间接自适应模糊控制康复机器人跟踪轨迹变化曲线; 图 11为两者误差变化曲线.由图 10可见, 期望轨迹(髋关节角度)出现先减小后增大再减小的变化趋势, 与被试者伸展/屈曲运动相对应, 同样符合第1.3.2节中式(25)所描述的规律.由图 11可见, 轨迹跟踪误差曲线在刚开始时达到最大值后迅速减小, 与之对应的是间接自适应模糊控制器的建立过程; 在7 s和12 s附近误差突然变大, 这是由于为验证所设计控制器的有效性, 人为引入了干扰项所导致的; 整个过程误差保持在$\pm $0.04范围内, 能够满足下肢康复机器人在线自适应控制的要求.
3. 结论
本文提出了一种基于sEMG和力反馈信息融合的下肢康复机器人自适应人机交互控制策略.在建立人机系统动力学模型基础上, 利用提取下肢运动过程中拮抗肌包络差值的方法识别出人体运动意图, 结合人机交互力进行sEMG-IF信息融合, 在线规划下肢康复机器人运动轨迹, 进一步利用间接自适应模糊控制器跟踪期望轨迹, 实现了下肢康复机器人自适应人机交互控制.通过5名被试者伸展/屈曲运动实验, 采集运动相关的股二头肌和股外侧肌sEMG以及足底压力信息, 并运用所提方法进行下肢康复机器人自适应人机交互控制.结果表明, 本文方法能够在线规划出体现人体运动意图和人机交互力信息的运动轨迹, 并以$\pm $0.04范围内的误差自适应跟踪下肢康复机器人运动轨迹.
-
图 7 股二头肌与股外侧肌包络差值
Fig. 7 The envelope difference of biceps and lateral muscles of femoral head
表 1 人体运动意图识别结果
Table 1 Results of human motion intent recognition
被试者 识别总数(个) 识别正确数(个) 识别率(%) A 348 332 95.40 B 315 301 95.55 C 296 285 96.28 D 308 291 94.48 E 286 277 96.85 
下载: 导出CSV
-
[1] Lo A C, Guarino P D, Richards L G, Haselkorn J K, Wittenberg G F, Federman D G, Ringer R J, Wagner T H, Krebs H I, Volpe B T, Bever C T, Bravata D M, Duncan P W, Corn B H, Maffucci A D, Nadeau S E, Conroy S S, Powell J M, Huang G D, Peduzzi P. Robot-assisted therapy for long-term upper-limb impairment after stroke. New England Journal of Medicine, 2010, 362(19):1772-1783 doi: 10.1056/NEJMoa0911341 [2] Klamroth-Marganska V, Blanco J, Campen K, Curt A, Dietz V, Ettlin T, Felder M, Fellinghauer B, Guidali M, Kollmar A, Luft A, Nef T, Schuster-Amft C, Stahel W, Riener R. Three-dimensional, task-specific robot therapy of the arm after stroke:a multicentre, parallel-group randomised trial. The Lancet Neurology, 2014, 13(2):159-166 doi: 10.1016/S1474-4422(13)70305-3 [3] Meng W, Liu Q, Zhou Z D, Ai Q S, Sheng B, Xie S Q. Recent development of mechanisms and control strategies for robot-assisted lower limb rehabilitation. Mechatronics, 2015, 31:132-145 doi: 10.1016/j.mechatronics.2015.04.005 [4] 彭亮, 侯增广, 王卫群.康复机器人的同步主动交互控制与实现.自动化学报, 2015, 41(11):1837-1846 http://www.aas.net.cn/CN/abstract/abstract18759.shtmlPeng Liang, Hou Zeng-Guang, Wang Wei-Qun. Synchronous active interaction control and its implementation for a rehabilitation robot. Acta Automatica Sinica, 2015, 41(11):1837-1846 http://www.aas.net.cn/CN/abstract/abstract18759.shtml [5] 佟丽娜, 侯增广, 彭亮, 王卫群, 陈翼雄, 谭民.基于多路sEMG时序分析的人体运动模式识别方法.自动化学报, 2014, 40(5):810-821 http://www.aas.net.cn/CN/abstract/abstract18349.shtmlTong Li-Na, Hou Zeng-Guang, Peng Liang, Wang Wei-Qun, Chen Yi-Xiong, Tan Min. Multi-channel sEMG time series analysis based human motion recognition method. Acta Automatica Sinica, 2014, 40(5):810-821 http://www.aas.net.cn/CN/abstract/abstract18349.shtml [6] Pittaccio S, Viscuso S. An EMG-controlled SMA device for the rehabilitation of the ankle joint in post-acute stroke. Journal of Materials Engineering and Performance, 2011, 20(4-5):666-670 doi: 10.1007/s11665-010-9826-7 [7] Song R, Tong K Y, Hu X L, Zhou W. Myoelectrically controlled wrist robot for stroke rehabilitation. Journal of Neuroengineering and Rehabilitation, 2013, 10(1):Article No.52 doi: 10.1186/1743-0003-10-52 [8] 胡进, 侯增广, 陈翼雄, 张峰, 王卫群.下肢康复机器人及其交互控制方法.自动化学报, 2014, 40(11):2377-2390 http://www.aas.net.cn/CN/abstract/abstract18514.shtmlHu Jin, Hou Zeng-Guang, Chen Yi-Xiong, Zhang Feng, Wang Wei-Qun. Lower limb rehabilitation robots and interactive control methods. Acta Automatica Sinica, 2014, 40(11):2377-2390 http://www.aas.net.cn/CN/abstract/abstract18514.shtml [9] Saglia J A, Tsagarakis N G, Dai J S, Caldwell D G. Control strategies for patient-assisted training using the ankle rehabilitation robot (ARBOT). IEEE/ASME Transactions on Mechatronics, 2013, 18(6):1799-1808 doi: 10.1109/TMECH.2012.2214228 [10] Duschau-Wicke A, Von Zitzewitz J, Caprez A, Lunenburger L, Riener R. Path control:a method for patient-cooperative robot-aided gait rehabilitation. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2010, 18(1):38-48 doi: 10.1109/TNSRE.2009.2033061 [11] Koopman B, van Asseldonk E H F, van der Kooij H. Selective control of gait subtasks in robotic gait training:foot clearance support in stroke survivors with a powered exoskeleton. Journal of Neuroengineering and Rehabilitation, 2013, 10(1):Article No.3 doi: 10.1186/1743-0003-10-3 [12] 刘金琨.机器人控制系统的设计与MATLAB仿真.北京:清华大学出版社, 2008. 169-172Liu Jin-Kun. Design and Simulation of Robot Control System Based on MATLAB. Beijing:Tsinghua University Press, 2008. 169-172 [13] 谢平, 宋妍, 郭子晖, 陈晓玲, 吴晓光, 苏玉萍, 杜义浩.中风康复运动中肌肉异常耦合分析.生物医学工程学杂志, 2016, 33(2):244-254 doi: 10.7507/1001-5515.20160043Xie Ping, Song Yan, Guo Zi-Hui, Chen Xiao-Ling, Wu Xiao-Guang, Su Yu-Ping, Du Yi-Hao. Analysis of abnormal muscular coupling during rehabilitation after stroke. Journal of Biomedical Engineering, 2016, 33(2):244-254 doi: 10.7507/1001-5515.20160043 [14] Hu J, Hou Z G, Peng L, Peng L, Gu N. sEMG-Based single-joint active training with iLeg-a horizontal exoskeleton for lower limb rehabilitation. In: Proceedings of the 21st International Conference. Kuching, Malaysia: Springer, 2014. 535-542 期刊类型引用(20)
1. 邢艳丽. 下肢康复训练机器人迭代学习控制探究. 上海轻工业. 2024(01): 186-188 . 
百度学术2. 董泽萍，仇大伟，刘静. 表面肌电人体下肢动作识别预测进展. 计算机工程与应用. 2023(14): 39-50 . 百度学术3. 高一聪，王彦坤，郑浩，魏喆，谭建荣. 下肢康复外骨骼的步态交互设计研究进展. 机械工程学报. 2023(17): 175-188 . 百度学术4. 邢艳丽，亢瑾. 下肢康复训练机器人膝踝协调运动控制方法. 机械设计与制造. 2022(03): 287-291 . 百度学术5. 刘克平，滕召纬，孙中波，李婉婷. 基于自适应模糊神经网络的下肢关节角度估计. 计算机仿真. 2022(09): 456-461 . 百度学术6. 黎颢宇，冯永飞，刘国伟，彭旭，牛建业，尹东颖. 坐卧式下肢康复机器人的机构设计与仿真. 机械制造. 2022(07): 35-40 . 百度学术7. 程帆，董宇欣. 基于力阻抗模型的上肢康复机器人交互控制系统设计. 计算机测量与控制. 2021(01): 111-114+125 . 百度学术8. 刘景森，吉宏远，李煜. 基于改进蝙蝠算法和三次样条插值的机器人路径规划. 自动化学报. 2021(07): 1710-1719 . 本站查看9. 梁旭，王卫群，苏婷婷，侯增广，何广平，任士鑫，石伟国. 下肢康复机器人的主动柔顺自适应交互控制. 机器人. 2021(05): 547-556 . 百度学术10. 程洪，黄瑞，邱静，王艺霖，邹朝彬，施柯丞. 康复机器人及其临床应用综述. 机器人. 2021(05): 606-619 . 百度学术11. 刘佳，温宇强. 基于双模糊PID算法的凹版印刷机多轴控制方法. 数字印刷. 2021(06): 57-64 . 百度学术12. 雷东威，陈彬，雷静桃. 穿戴式外骨骼机器人肌电信号采集与处理. 计量与测试技术. 2021(12): 11-15 . 百度学术13. 黄智晖. 基于ADAMS的康复机器人建模与运动仿真的研究. 电子测量技术. 2020(09): 127-132 . 百度学术14. 王永琦，江潇潇. 基于混合灰狼算法的机器人路径规划. 计算机工程与科学. 2020(07): 1294-1301 . 百度学术15. 滕召纬，孙中波，刘克平. 基于自适应模糊神经网络的下肢关节运动意图估计. 长春工业大学学报. 2020(06): 558-562 . 百度学术16. 张弼，姚杰，赵新刚，谈晓伟. 一种基于肌电信号的自适应人机交互控制方法. 控制理论与应用. 2020(12): 2560-2570 . 百度学术17. 吴佳佳，刘箴，许辉煌，刘婷婷，王瑾，柴艳杰. 面向康复训练的严肃游戏人机交互设计. 系统仿真学报. 2019(05): 909-918 . 百度学术18. 刘中华. 下肢康复机器人的交互控制研究. 人人健康. 2019(17): 282-283 . 百度学术19. 李剑锋，张子康，张雷雨，董明杰，左世平，张凯. 并联踝康复机器人的系统搭建与运动控制策略. 中南大学学报(自然科学版). 2019(11): 2753-2762 . 百度学术20. 刘中华. 基于主动交互控制的家用康复治疗床自适应阻抗控制实现. 人人健康. 2019(23): 468 . 百度学术其他类型引用(32)
-
下载:
计量
- 文章访问数: 2482
- HTML全文浏览量: 1300
- PDF下载量: 1058
- 被引次数: 52
