目前,我国人口老龄化的发展趋势不容乐观,预计老年人口数量在2035年将突破3亿,2053年将达到峰值4.87亿,占届时全球老年人口的1/4[1 ] 。随着人口老龄化的加剧,因机体老化、脑卒中等原因造成下肢运动功能障碍或残疾的患者日益增多。因此,研究下肢康复机器人对实现患者下肢运动功能恢复以及解决我国“患者多、康复医师少”的社会问题具有重要意义。
国内外学者在康复机器人研究领域取得了丰硕的成果。瑞士Hocoma医疗器械公司与苏黎世Balgrist医学院联合研制了Lokomat步行康复训练机器人[2 ] ,它可以辅助患者模拟健康人的行走状态。饶玲军等[3 ] 设计的外骨骼康复机器人包括外骨骼本体和1副可实现身体重心控制和平衡的拐杖,其能使患者实现站立与行走。Aggogeri等[4 ] 研发的六自由度并联步态训练康复机器人Rutgers Ankle通过调节2个踏板的相对位置与姿态来协助患者进行平地行走、转弯行走以及平衡能力康复训练。Cudby等[5 ] 研制的可穿戴式下肢康复机器人REX可以帮助患者外出行走。Chen等[6 ] 研制了一款名为CUHK-EXO的下肢外骨骼康复机器人,它可以协助患者站立、坐下和行走等。上述穿戴式下肢康复机器人主要服务于康复后期的患者,能够辅助他们恢复行走能力,而对处于康复初期的患者存在一定局限性。
针对上述问题,部分学者提出了末端牵引式下肢康复机器人,其在辅助处于康复初期的患者恢复运动功能方面具有较大潜力。例如:Bouri等[7 ] 研制的牵引式下肢康复机器人Lambda采用左右对称的并联机械结构,能够以水平或垂直的姿态布置,适用于辅助康复初期患者进行矢状面内的康复训练;Monaco等[8 ] 研制的NEUROBike由2个相互垂直安装的导轨和踝部的转动装置组成,能够辅助卧床的患者模拟健康人的行走动作。相较于穿戴式下肢康复机器人,末端牵引式下肢康复机器人还具有轻量化与小型化的特点,且机械结构简单,易于控制,在训练过程中能够保证患者的安全性与舒适性。然而,上述学者研制的末端牵引式下肢康复机器人的训练方式较为单一,并不十分完善。
为满足患者在不同康复时期的训练需求,针对现有下肢康复机器人训练方式单一的问题,从人体下肢运动机理出发,基于仿生学提出一种可实现卧姿、坐姿训练模式的末端牵引式下肢康复机器人,其兼具并联机构和串联机构的优势:结构简单、承载能力强、正解简单、运动灵活和易于控制等。此外,相对于穿戴式下肢康复机器人,所设计的机器人可克服因重力作用下机器人零部件与人体发生相对滑动而导致的人机运动不一致的缺点,进而降低对穿戴性和舒适性的要求。然后,通过对机器人的运动学正、逆解进行分析,并以其安全工作空间为目标函数,通过遗传算法对其机构参数进行优化。接着,根据康复护理相关理论对机器人的康复训练轨迹进行规划。最后,搭建实验平台并开展运动捕捉实验,以验证理论分析的正确性。
1 牵引式下肢康复机器人机构设计及自由度计算
1.1 人体下肢运动机理分析
下肢是人体结构的重要组成部分,其由骨骼、关节、肌肉以及韧带组成,主要起支撑人体重量和辅助人体行走的作用。对人体下肢运动机理进行分析,可为机器人机构设计提供理论基础。
根据人体解剖学原理,利用假想的3个相互垂直的平面将人体分为8个部分。鉴于拟设计的牵引式下肢康复机器人仅针对患者的髋、膝关节进行康复训练,故根据基本切面和基本轴线对人体下肢髋、膝关节的运动进行描述。
髋关节是典型的杵臼关节,从形态上看其近似为球窝关节,其运动形式为矢状面内的屈伸运动、水平面内的内外旋运动以及冠状面内的外展内收运动。膝关节属于滑车关节,其运动形式为矢状面内的屈伸运动。
根据上述分析,在构建人体下肢运动等效模型时,可将髋关节等效为虎克铰链,将膝关节等效为回转副,将足部视作下肢末端。
1.2 机器人机构设计
基于人体下肢运动等效模型,对末端牵引式下肢康复机器人的机构进行设计。由生物学理论可知,人体下肢运动是靠骨骼肌的收缩与舒张实现的,骨骼肌长度的改变可实现相邻2个骨骼的相对转动。故根据仿生学原理,在虎克铰链处布置2条UPS支链来模拟人体大腿的骨骼肌肉,在大腿与小腿之间加入1条RPR支链来模拟人体小腿的肌肉[9 ] ,且加入UPS支链和RPR支链可以提高机器人的承载能力,上述支链均以移动副作为驱动副,则以上部分组成(2-UPS&U)+(R&RPR)机构。然后,在(2-UPS&U)+(R&RPR)机构末端加入2条RR支链,通过伺服电机驱动齿轮箱内齿轮的传动实现左、右脚踏板的镜像开合,以辅助患者进行髋关节内外旋的康复训练;脚踏板跟随人体足部运动,不设置驱动。最终组成(2-UPS&U)+(R&RPR)+RR混联机构,如图1 (a)所示,对应的机器人三维模型如图1 (b)所示。
图1
图1
牵引式下肢康复机器人机构简图及三维模型
Fig.1
Mechanism sketch and three-dimensional model of traction lower limb rehabilitation robot
1.3 机器人机构自由度计算
利用修正的G-K公式计算(2-UPS&U)+(R&RPR)机构的自由度,计算式如下:
M = d n - g - 1 + ∑ i = 1 g f i + v - η (1)
式中:M 为机构的自由度;d 为机构的阶数;n 为机构中构件的数量;g 为机构中运动副的数量;fi 为第i 个运动副的自由度;v 为机构的冗余约束数;η 为机构的局部自由度。
由于所设计的机器人并联部分中的UPS支链的自由度为6,不会对OD 杆产生约束,故该支链为无约束支链;RPR支链始终与OD 杆、DF 杆处于同一平面,故该支链产生的约束为冗余约束。由此可知,(2-UPS&U)+(R&RPR)机构的阶数d =6,构件数n =9,运动副数g =11,所有运动副自由度的代数和为20,冗余约束数v =1,且该机构不存在局部自由度,故η = 0 。将上述参数代入式(1)即可求得该机构的自由度:
M = 6 × 9 - 11 - 1 + 20 + 1 - 0 = 3 (2)
此外,2条RR支链的运动具有一致性,故这2条支链整体的自由度为2。综上,所设计的牵引式下肢康复机器人机构的自由度为5。
2 牵引式下肢康复机器人运动学分析
为计算所设计牵引式下肢康复机器人的运动学正、逆解,构建如图2 所示的机构坐标系。其中:固定坐标系O-XYZ 的原点为虎克铰链的2条转动轴线的交点,X 轴与U支链中虎克铰链的第一转动轴线重合,沿OA 2 方向;Y 轴与U支链中虎克铰链的第二转动轴线重合,沿OA 1 方向;Z 轴由右手定则确定。动坐标系B-XB YB ZB 的原点位于OD 杆上某固定位置处,XB 轴沿BB 2 方向,YB 轴沿BB 1 方向,ZB 由右手定则确定。动坐标系D-XD YD ZD 的XD 轴与其原点D 处的回转副轴线平行,ZD 轴沿OD 方向。动坐标系F-XF YF ZF 的XF 轴沿FM 2 方向,ZF 轴沿DF 方向。动坐标系F 1 -XF 1 YF 1 ZF 1 的XF 1 轴垂直于原点F 1 处的回转副轴线,YF 1 轴与原点F 1 处的回转副轴线平行。动坐标系E 1 -XE 1 YE 1 ZE 1 的XE 1 轴垂直于原点E 1 处的回转副轴线,YE 1 轴与原点E 1 处的回转副轴线平行。动坐标系F 2 -XF 2 YF 2 ZF 2 与E 2 -XE 2 YE 2 ZE 2 同理,不再赘述。
图2
图2
牵引式下肢康复机器人机构运动学坐标系
Fig.2
Kinematics coordinate system of traction lower limb rehabilitation robot mechanism
2.1 运动学正解
运动学正解是指在已知各个关节输入参数的情况下求解的末端执行器位置参数。为表述方便,定义OD 杆的长度为l 1 ,DF 杆的长度为l 2 ,FM 1 杆、FM 2 杆的长度为l 3 ,M 1 F 1 杆、M 2 F 2 杆的长度为l 4 ,F 1 E 1 杆、F 2 E 2 杆的长度为l 5 。首先,分析左踏板连接处E 1 点相对于固定坐标系原点O 的位置。动坐标系D -XD YD Z D 相对于固定坐标系O -XYZ 的位姿可看作先绕X 轴旋转α ,然后绕Y 轴旋转β ,最后沿OD 杆方向平移l 1 得到。动坐标系F -XF YF ZF 相对于动坐标系D -XD YD ZD 的位姿可看作先绕XD 轴旋转γ ,再沿DF 方向平移l 2 得到。动坐标系F 1 -XF 1 YF 1 ZF 1 可由动坐标系F -XF YF ZF 平移得到。动坐标系E 1 -XE 1 YE 1 ZE 1 相对于动坐标系F 1 -XF 1 YF 1 ZF 1 的位姿可看作先绕YF 1 轴旋转- θ ,再沿ZF 1 轴平移l 5 得到。由此可得,动坐标系E 1 -XE 1 YE 1 ZE 1 相对于固定坐标系O -XYZ 的齐次变换矩阵E 1 O T 为:
E 1 O T = T D O T F D T F 1 F T E 1 F 1 = s θ s β c α + γ + c β c θ s α + γ s β s β c θ c α + γ E 1 x - s α + γ s θ c α + γ - s α + γ c θ E 1 y s θ c β c α + γ - s β c θ s α + γ c β c θ c β c α + γ + s β s θ E 1 z 0 0 0 1 (3)
E 1 x = l 5 c θ + l 2 + l 4 s β c α + γ - l 3 c β + l 1 c α s β - l 5 s θ E 1 y = - l 2 + l 4 + l 5 c θ s α + γ - l 1 s α E 1 z = l 5 c θ + l 2 + l 4 c β c α + γ + l 3 s β + l 1 c α c β + l 5 s β s θ
式中: T 为各坐标系之间的变换矩阵;c表示cos,s表示sin,下文矩阵中同。
同理可得,动坐标系E 2 -XE 2 YE 2 ZE 2 相对于固定坐标系O-XYZ 的齐次变换矩阵E 2 O T 为:
E 2 O T = T D O T F D T F 2 F T E 2 F 2 = c β c θ - s θ s β c α + γ s α + γ s β c β s θ + s β c θ c α + γ E 2 x s α + γ s θ c α + γ - s α + γ c θ E 2 y - s θ c β c α + γ - s β c θ s α + γ c β c θ c β c α + γ E 2 z 0 0 0 1 (4)
E 2 x = l 5 c θ + l 2 + l 4 s β c α + γ + l 3 c β + l 1 c α s β + l 5 c β s θ E 2 y = - l 2 + l 4 + l 5 c θ s α + γ - l 1 s α E 2 z = l 5 c θ + l 2 + l 4 c β c α + γ - l 3 s β + l 1 c α c β - l 5 s β s θ
当并联部分(2-UPS&U)的静、动坐标系重合时,动平台B 点相对于静平台OA 1 A 2 的转动可看作动坐标系B -XB YB ZB 绕固定坐标系O-XYZ 的X 轴、Y 轴转动,则其旋转矩阵B O R 为:
B O R = c β s α s β c α s β 0 c α - s α - s β s α c β c α c β (5)
A 1 、A 2 、O 点在固定坐标系O-XYZ 中的位置矢量 m 1 、 m 2 、 m 3 分别为:
m 1 = 0 a 1 0 T m 2 = a 2 0 0 T m 3 = 0 0 0 T (6)
B 1 、B 2 、B 点在动坐标系B-XB YB ZB 中的位置矢量 n 1 、 n 2 、 n 3 分别为:
n 1 = 0 b 1 0 T n 2 = b 2 0 0 T n 3 = 0 0 0 T (7)
在固定坐标系O-XYZ 下,建立并联部分的位置约束方程:
l O B + B O R n j = m j + A j B j ( j = 1 , 2 ) (8)
l O B = l O B c α s β - l O B s α l O B c α c β T
式中: l OB 为OB 杆的长度; A j 、 B j 分别为Aj 、Bj 点在固定坐标系O -XYZ 中的位置矢量。
代入具体数值并化简,可得电动推杆1,2的长度l 01 、l 02 :
l 01 = a 1 2 + b 1 2 + l O B 2 - 2 a 1 b 1 c o s α + 2 l O B a 1 s i n α l 02 = a 2 2 + b 2 2 + l O B 2 - 2 a 2 b 2 c o s β - 2 l O B a 2 c o s α s i n β (9)
c o s ∠ C 2 D C 3 = b 3 2 + m 2 2 + m 1 2 + a 3 2 - l 03 2 2 b 3 2 + m 2 2 m 1 2 + a 3 2 (10)
式中:l 03 为电动推杆3的长度;a 3 为CC 3 杆的长度;b 3 为C 1 C 2 杆的长度,m 1 为CD 杆的长度,m 2 为C 1 D 杆的长度。
γ = π - a r c c o s b 3 2 + m 2 2 + a 3 2 + m 1 2 - l 03 2 2 b 3 2 + m 2 2 a 3 2 + m 1 2 (11)
l 03 = a 3 2 + b 3 2 + m 1 2 + m 2 2 + K (12)
K = 2 m 1 2 + a 3 2 m 2 2 + b 3 2 c o s γ + a r c t a n b 3 m 2 + a r c t a n a 3 m 1
2.2 运动学逆解
通过变换将动坐标系E 1 -XE 1 YE 1 ZE 1 与E 2 -XE 2 YE 2 ZE 2 转换为F-XF YF ZF ,得到F 点在固定坐标系O-XYZ 中的位置(Fx , Fy , Fz ),从而计算得到α 、β 与γ :
α = a r c c o s l 1 2 - F y 2 - C 2 - l 2 2 2 l 1 F y 2 + C 2 - a r c t a n F y c o s β F z β = a r c t a n F x F z γ = a r c c o s F y + l 1 s i n α l 2 - a r c t a n c o t α (13)
C = F z c o s β
将求得的3个转角的值代入式(9)和式(12),可得到3个电动推杆的长度。
3 牵引式下肢康复机器人机构参数优化
为了使所设计的牵引式下肢康复机器人可以更好地满足患者的康复训练需求,以及使患者在训练过程中更加舒适,对其机构参数进行优化。由于康复初期人体关节的活动度较小[10 ] ,故以康复中后期坐姿训练状态下人体下肢末端与机器人末端工作空间重合度(即机器人的安全工作空间[11 ] )为目标函数,以机器人机构参数为设计变量,利用MATLAB遗传算法工具箱[12 ] 进行优化计算。
由于该机器人机构中的2条RR支链左右对称,故以人体左下肢与机器人左踏板为对象进行分析。
3.1 人体左下肢运动学分析
将机器人与坐姿状态下的人体下肢看作人机系统,以机器人固定坐标系O-XYZ 为人机系统的固定坐标系。根据人体下肢运动等效模型,建立如图3 所示的人机系统坐标系。其中:O 1 -XO 1 YO 1 ZO 1 表示髋关节坐标系,N-XN YN ZN 表示膝关节坐标系,M-XM YM ZM 表示足底坐标系;O 1 表示人体髋关节,N 表示人体膝关节,O 1 N 表示人体大腿,其长度为l 6 ;NM 表示人体小腿与足部,其长度为l 7 。
图3
图3
人机系统坐标系
Fig.3
Human-machine system coordinate system
令xO 1 表示人体髋部的横向宽度,yO 1 表示髋关节在固定坐标系中的Y 向位置,zO 1 表示坐姿状态下人体髋关节距离地面的高度。则动坐标系O 1 -XO 1 YO 1 ZO 1 的原点O 1 相对于固定坐标系O-XYZ 的原点O 的位置为(-xO 1 , -yO 1 , zO 1 )。
动坐标系N-XN YN ZN 相对于动坐标系O 1 -XO 1 YO 1 ZO 1 的位姿可看作先绕XO 1 轴旋转θ 1 ,再绕ZO 1 轴旋转θ 2 ,最后沿O 1 N 方向平移l 6 得到。动坐标系M-XM YM ZM 相对于动坐标系N-XN YN ZN 的位姿可看作先绕XN 轴旋转-θ 3 ,再沿NM 方向平移l 7 得到。故动坐标系M-XM YM ZM 相对于人机系统固定坐标系的齐次变换矩阵M O T 为:
M O T = c θ 2 - c θ 3 s θ 2 - s θ 2 s θ 3 M x c θ 1 s θ 2 s θ 1 s θ 3 + c θ 1 c θ 2 c θ 3 c θ 1 c θ 2 s θ 3 - c θ 3 s θ 1 M y s θ 1 s θ 2 c θ 2 c θ 3 s θ 1 - c θ 1 s θ 3 c θ 1 c θ 3 + c θ 2 s θ 1 s θ 3 M z 0 0 0 1 (14)
M x = - x 1 - l 6 s θ 2 - l 7 c θ 3 s θ 2 M y = - y 1 + l 6 c θ 1 c θ 2 + l 7 s θ 1 s θ 3 + l 7 c θ 1 c θ 2 c θ 3 M z = z 1 + l 6 c θ 2 s θ 1 - l 7 c θ 1 s θ 3 + l 7 c θ 2 c θ 3 s θ 1
参考《中国成年人人体尺寸》[13 ] ,以身高为1 754 mm、体重为75 kg的男性模型为参考,其对应的下肢参数如下:大腿长度为496 mm、小腿长度为396 mm,足部厚度为75 mm,胯部宽度为335 mm。故取xO 1 =167.5 mm,zO 1 =440 mm,l 6 =496 mm,l 7 =471 mm。根据人机工程学选取yO 1 =1 420 mm。规定坐姿状态下人体下肢各关节的康复训练角度,如表1 所示。
0 ° ≤ θ 1 ≤ 50 ° - 20 ° ≤ θ 2 ≤ 20 ° 0 ° ≤ θ 3 ≤ 90 ° (15)
将上述参数值代入式(14),可得人体左下肢末端在人机系统固定坐标系中的位置。
令Mx 、My 、Mz 与E 1 x 、E 1 y 、E 1 z 分别对应相等,求得人体左下肢的运动学逆解:
θ 1 = 2 a r c t a n 2 c l 6 + P b 2 + 2 b l 6 + c 2 + l 6 2 - l 7 2 θ 2 = a r c s i n a l 6 + l 7 c o s θ 3 θ 3 = 2 a r c t a n Q b 2 + c 2 - l 6 2 + 2 l 6 l 7 - l 7 2 (16)
a = x O 1 + E 1 x b = y O 1 + E 1 y c = E 1 z - z O 1 P = - ( b 2 + c 2 ) 2 + 2 ( b 2 + c 2 ) ( l 6 2 + l 7 2 ) - ( l 6 2 - l 7 2 ) 2 Q = ( b 2 + c 2 - l 6 2 + 2 l 6 l 7 - l 7 2 ) ( - b 2 - c 2 + l 6 2 + 2 l 6 l 7 + l 7 2 )
当机器人的输入已知时,可通过式(16)求得人体下肢髋、膝关节的角度。
3.2 设计变量选取
根据左踏板连接处的运动学正解,以机器人机构的主要杆长作为设计变量,即:
x = x 1 x 2 x 3 x 4 x 5 T = l 1 l 2 l 3 l 4 l 5 T (17)
3.3 目标函数确定
机器人末端与人体下肢末端工作空间的重合度越高,表明康复训练效果越好。故以机器人末端与人体下肢末端工作空间对应散点的距离之和最小为目标,建立目标函数S :
S = ∑ r = 1 k X r - x r 2 + Y r - y r 2 + Z r - z r 2 (18)
式中:Xr 、Yr 、Zr 和xr 、yr 、zr 分别为机器人末端和人体下肢末端在固定坐标系O -XYZ 中的坐标;k 为工作空间的散点数量。
3.4 约束条件
机构参数的变化会引起机器人工作空间的变化。首先,研究虎克铰链以及回转副的转角范围对机器人工作空间的影响,以“点集数”作为工作空间的评价指标[14 ] :点集数越大,工作空间越大。
按照固定步长依次改变4个转角的范围,其他参数值全部采用初选值,然后通过单一变量法来研究各转角范围对机器人工作空间的影响。最终得到机器人工作空间随各转角范围的变化曲线,如图4 所示。
图4
图4
各转角范围对机器人工作空间的影响
Fig.4
Influence of each angle range on robot workspace
由图4 可知,虎克铰链与回转副的转角范围越大,机器人的工作空间越大。另外,根据人机工程学考虑训练过程中患者的安全性以及电动推杆的技术要求,最终规定虎克铰链及回转副的转角范围为:
0 ° ≤ α ≤ 40 ° - 20 ° ≤ β ≤ 20 ° 47 ° ≤ γ ≤ 101 ° 0 ° ≤ θ ≤ 60 ° (19)
390 ≤ l 1 ≤ 420 m m 340 ≤ l 2 ≤ 400 m m 90 ≤ l 3 ≤ 110 m m 30 ≤ l 4 ≤ 45 m m 150 ≤ l 5 ≤ 180 m m (20)
采用蒙特卡洛法[15 ] 求解优化后的人机系统工作空间,结果如图5 所示。其中:红色部分代表人体左下肢末端的工作空间,蓝色部分代表机器人左踏板的工作空间,两者相交的部分为机器人的安全工作空间。由此可得,矢状面内机器人的安全工作空间区域如图6 所示。
图5
图5
人机系统工作空间
Fig.5
Human-machine system workspace
图6
图6
人机系统矢状面工作空间区域
Fig.6
Sagittal workspace area of human-machine system
由于下肢康复机器人以辅助患者进行矢状面内的康复训练为主[16 ] ,利用式(21)求解人体下肢矢状面内的有效工作空间比ξ [17 ] :
ξ = S t S w (21)
式中:S t 为机器人安全工作空间区域;S w 为人体下肢等效模型的工作空间区域。
根据图6 和式(21),求得ξ =0.71,说明优化结果良好。此外,由于本文所研究的机器人为末端牵引式下肢康复机器人,其安全工作空间无法与人体下肢等效模型的工作空间完全一致,而根据康复护理相关理论[18 ] ,在图6 所示的人体末端与机器人末端工作空间的重合区域内,机器人可以帮助人体下肢进行多种康复训练,即机器人对人体有较好的适应性,可以满足其康复需求。
4 牵引式下肢康复机器人轨迹规划
在康复训练的过程中,要求机器人尽量作速度、加速度连续平滑的无冲击运动。五次多项式插值能够解决三次多项式插值的速度变化不平滑且加速度存在突变的问题,故本文采用五次多项式插值法对所设计的牵引式下肢康复机器人进行轨迹规划。
4.1 卧姿康复训练模式轨迹规划
患者在康复初期身体虚弱,无法支撑其体重,处于这个时期的患者适合卧姿康复训练,通过保持正常的关节活动来恢复患者肢体的运动能力,以避免关节变形和肌肉萎缩。根据康复护理相关理论,采用CPM(continous passive motion,持续被动运动)方法[18 ] 进行康复训练。在康复训练过程中,小腿始终保持水平,髋关节与膝关节的角度变化一致。
为了确保运动的灵活性和舒适性,通常要求髋关节与膝关节的初始夹角为10 ° [18 ] ,规定人体髋、膝关节的运动角度为10 ° ~ 50 ° ,运动起始与终止处的速度与加速度均为0,故髋关节角度的五次插值多项式为:
φ t = 0.174 5 + 0.069 t 3 - 0.001 t 4 + 0.000 041 88 t 5 (22)
则可得在人机系统固定坐标系中,人体下肢末端Y 、Z 方向的位移为:
y t = - 147 0 + l 6 c o s φ t + l 7 z t = l 6 s i n φ t (23)
4.2 坐姿康复训练模式轨迹规划
4.2.1 圆周运动训练轨迹规划
当患者处于康复中后期时,适合坐姿康复训练模式。圆周运动训练即人体矢状面内的康复训练,与脚踏车运动类似,能够帮助患者的髋、膝关节进行屈伸运动[19 ] 。以工作空间中某一圆周轨迹为例,其圆心相对人机系统固定坐标系原点的位置(x 0 , y 0 , z 0 )=(-90, -665, 470) mm,半径R =115 mm,规定运动时间为10 s。令圆周运动的起始角度为360°,终止角度为0°,起始、终止处的速度与加速度均为0,由此得到圆心角位移的五次插值多项式:
ε ( t ) = 6.282 3 - 0.062 8 t 3 + 0.009 4 t 4 - 0.000 38 t 5 (24)
同理可得,在人机系统固定坐标系中,人体下肢末端Y、Z 方向的位移为:
y t = R c o s α t + y 0 z t = R s i n α t + z 0 (25)
4.2.2 螺旋运动训练轨迹规划
螺旋运动由直线运动与圆周运动组合形成,相较于圆周运动,其可以帮助患者实现空间内的康复训练[20 ] 。定义在人机系统固定坐标系中,圆周分运动的圆心坐标为(x 0 , y 0 , z 0 )=(-200, -630, 600) mm,半径R =70 mm;令圆周分运动的起始角度为720°,终止角度为0°。直线分运动的起始位置为(-200, -630, 600) mm,终止位置为(0, -630, 600) mm;2个分运动的起始与终止处的速度、加速度均为0。规定螺旋运动时间为15 s。由此可得,螺旋运动的2个分运动的五次插值多项式为:
δ t = 12.57 + 0.001 9 t 3 + 0.001 2 t 4 + 0.000 031 6 t 5 l t = - 200 + 0.592 6 t 3 - 0.068 1 t 4 + 0.001 8 t 5 (26)
同理可得,在人机系统固定坐标系中,人体下肢末端X 、Y、Z 方向的位移为:
x t = l t y t = R c o s δ t + y 0 z t = R s i n δ t + z 0 (27)
根据式(9)、式(12)、式(13)可得上述康复训练轨迹下电动推杆长度的变化情况。再根据式(16)可进一步计算出上述康复训练轨迹下人体髋、膝关节运动角度的理论变化值。
5 牵引式下肢康复机器人训练实验
为了验证牵引式下肢康复机器人设计与优化的合理性以及轨迹规划的正确性,根据轨迹规划结果对其电动推杆进行控制,并通过Xsens MVN系统采集康复训练过程中人体髋、膝关节的运动角度数据。由于目前机器人样机的底部没有设置升降装置,故现阶段开展卧姿康复训练实验时,被试者平躺在地面上,以模拟卧床期的患者。
机器人样机选用带上下限位开关的伺服电动推杆作为驱动,采用Copley公司的Accelnet伺服电机驱动器模块驱动电动推杆,选用基于STM32H750的主控板作为主驱动器,通过所搭建的伺服电机驱动系统来实现电动推杆位移的精密控制。考虑到安全性问题,在实验开始前将被试者的足部与机器人样机的脚踏板绑在一起,从而避免两者脱离。
图7 和图8 所示分别为被试者在卧姿和坐姿状态下进行康复训练的实验现场。
图7
图7
卧姿康复训练实验现场
Fig.7
Experimental site of rehabilitation training with lying posture
图8
图8
坐姿康复训练实验现场
Fig.8
Experimental site of rehabilitation training with sitting posture
3种康复训练轨迹下人体髋、膝关节运动角度的测量值与理论值对比分别如图9 至图11 所示。
图9
图9
CPM训练轨迹下人体髋、膝关节运动角度对比
Fig.9
Comparison of motion angles of human hip and knee joints under CPM training trajectory
图10
图10
圆周运动训练轨迹下人体髋、膝关节运动角度对比
Fig.10
Comparison of motion angles of human hip and knee joints under circular exercise training trajectory
图11
图11
螺旋运动训练轨迹下人体髋、膝关节运动角度对比
Fig.11
Comparison of motion angles of human hip and knee joints under spiral training trajectory
从图9 至图11 中可知,髋、膝关节运动角度的理论值与测量值基本一致;在3种训练轨迹下,髋关节运动角度的最大误差约为2.8°,膝关节运动角度的最大误差约为4.2°。产生误差的原因是在实验过程中被试者的腿部会不自主地抖动,且被试者与标准人体尺寸存在差异。另外,被试者的足底与机器人脚踏板之间可能会发生微小的相对移动,也会导致实验结果产生偏差。但总体来说,误差在许用范围之内,且人体髋、膝关节的活动度符合人体康复运动参数,验证了机器人的设计及其优化后的安全工作空间的合理性以及轨迹规划的正确性。
6 结 论
1)为满足下肢运动功能障碍患者在不同阶段的康复训练需求,针对现有下肢康复机器人训练方式单一的问题,根据人体下肢运动机理以及仿生学,设计了一种可实现卧姿、坐姿训练模式的牵引式下肢康复机器人。
2)基于对牵引式下肢康复机器人的运动学分析,将人体下肢末端与机器人末端工作空间的重合部分作为机器人的安全工作空间,并以此为目标函数,通过遗传算法对机器人的机构参数进行了优化;同时,根据优化结果绘制人机系统工作空间散点图,得到矢状面内人体下肢的有效工作空间比为0.71,表明优化结果良好。
3)根据康复护理相关理论,在机器人安全工作空间内分别开展卧姿、坐姿康复训练轨迹规划,并通过实验对理论分析结果进行了验证。结果表明,机器人设计与参数优化合理以及轨迹规划正确,可为人机交互系统的研究奠定基础。
4)本文仅针对被动训练进行了研究,后续将通过在人体与机器人的接触部位安装力/力矩传感器来获取人机之间的交互力,以进行人体运动意图识别,并在此基础上制定控制策略,以实现机器人助力训练以及主动训练。此外,还将通过理论分析与实验建立患者舒适度与康复效果的评价指标体系,以使机器人能够更好地辅助患者进行康复训练。
参考文献
View Option
[2]
TABEA A S , ANJA G , ROB L . Correction to: the FreeD module for the Lokomat facilitates a physiological movement pattern in healthy people: a proof of concept study
[J]. Journal of Neuroengineering & Rehabilitation , 2019 , 16 (1 ): 1 -13 .
[本文引用: 1]
[4]
AGGOGERI F , PELLEGRINI N , ADAMINI R . Functional design in rehabilitation: modular mechanisms for ankle complex
[J]. Applied Bionics and Biomechanics , 2016 , 2016 : 9707801 .
[本文引用: 1]
[5]
CUDBY K . Liberty autonomy independence
[J]. Engineering Insight , 2011 , 12 (1 ): 8 -14 .
[本文引用: 1]
[6]
CHEN B , ZHONG C H , ZHAO X , et al . A wearable exoskeleton suit for motion assistance to paralysed patients
[J]. Journal of Orthopaedic Translation , 2017 , 11 (3 ): 7 -18 .
[本文引用: 1]
[7]
BOURI M , LE GALL B , CLAVEL R . A new concept of parallel robot for rehabilitation and fitness: the Lambda
[C]//Proceedings of the 2009 International Conference on Robotics and Biomimetics . Guilin, Dec . 19 -23 , 2009 .
[本文引用: 1]
[8]
MONACO V , GALARDI G , COSCIA M , et al . Design and evaluation of NEUROBike: a neurorehabilitative platform for bedridden post-stroke patients
[J]. IEEE Transactions on Neural System and Rehablitation Engineering , 2012 , 20 (6 ): 845 -852 .
[本文引用: 1]
[9]
牛建业 .基于串并混联机构的四轮足步行机器人研究
[D].秦皇岛 :燕山大学 ,2018 :21 -23 .
[本文引用: 1]
NIU Jian-ye . Research on wheel-legged quadruped walking robot based on series-parallel mechanism
[D]. Qinhuangdao : Yanshan University , 2018 : 21 -23 .
[本文引用: 1]
[10]
[本文引用: 1]
WEI Yan-fen , CHEN Min-ru , DENG Min . Analysis of physical function status and influencing factors in patients with hip joint replacement at early rehabilitation
[J]. Journal of Nursing Administration , 2018 , 18 (8 ): 567 -571 .
DOI:10.3969/j.issn.1671-315x.2018.08.008
[本文引用: 1]
[11]
[本文引用: 1]
XIA Peng-peng , HAN Jian-hai , LI Liao-yuan , et al . A study on trajectory planning of an upper limb rehabilitation robot with distal traction in three-dimensional space
[J]. Chinese Journal of Rehabilitation Medicine , 2021 , 36 (11 ): 1414 -1419 .
DOI:10.3969/j.issn.1001-1242.2021.11.015
[本文引用: 1]
[12]
雷英杰 ,张善文 .MATLAB遗传算法工具箱及应用 [M].西安 :西安电子科技大学出版社 ,2014 :147 -155 .
[本文引用: 1]
LEI Ying-jie , ZHANG Shan-wen . MATLAB genetic algorithm toolbox and application [M]. Xi’an : Xidian University Press , 2014 : 147 -155 .
[本文引用: 1]
[13]
国家技术监督局 . 中国成年人人体尺寸国家标准 : [S].北京 :中国标准出版社 ,1988 : 6 -8 .
[本文引用: 1]
National Technical Supervision Bureau . Human dimensions of Chinese adults : [S]. Beijing : China Standard Press , 1988 : 6 -8 .
[本文引用: 1]
[14]
莫贤 ,陈文家 ,陈淑艳 .一种3-PRS并联机器人的姿态空间分析
[J].扬州大学学报(自然科学版) ,2009 ,12 (3 ):57 -61 .
[本文引用: 1]
MO Xian , CHEN Wen-jia , CHEN Shu-yan . Analysis on working space of a 3-PRS parallel manipulator
[J]. Journal of Yangzhou University (Natural Science Edition) , 2009 , 12 (3 ): 57 -61 .
[本文引用: 1]
[16]
[本文引用: 1]
GUO Bing-jing , HAN Jian-hai , LI Xiang-pan , et al . Personalized gait planning method for the lower-limb rehabilitation training robot with the physiotherapist interaction
[J]. Robot , 2018 , 40 (4 ): 479 -490 , 499 .
DOI:10.13973/j.cnki.robot.180139
[本文引用: 1]
[17]
陈佳丽 ,许勇 ,刘文彩 .线性驱动并联机构基于有效工作空间比的尺度优化
[J].机械设计与研究 ,2019 ,35 (1 ):42 -46 ,52 .
[本文引用: 1]
CHEN Jia-li , XU Yong , LIU Wen-cai . Geometric solution and scale optimization of parallel mechanism workspace
[J]. Mechine Design and Research , 2019 , 35 (1 ): 42 -46 , 52 .
[本文引用: 1]
[18]
GUO B , HAN J , LI X , et al . Research and design of a new horizontal lower limb rehabilitation training robot
[J]. International Journal of Advanced Robotic Systems , 2016 , 13 (1 ): 1 -10 .
[本文引用: 3]
[19]
黄高 ,张伟民 , MARCO Ceccarelli ,等 .一种新的康复与代步外骨骼机器人研究
[J].自动化学报 ,2016 ,42 (12 ):1933 -1942 .
[本文引用: 1]
HUANG Gao , ZHANG Wei-min , MARCO Ceccarelli , et al . Research of a new rehabilitation and assisting robot
[J]. Journal of Automation , 2016 , 42 (12 ): 1933 -1942 .
[本文引用: 1]
[20]
靳振楠 .空间四自由度下肢康复机器人设计与分析
[D].秦皇岛 :燕山大学 ,2019 :42 -43 .
[本文引用: 1]
JIN Zhen-nan . Design and analysis of spatial four-degree-of freedom lower limb rehabilitation robot
[D]. Qinhuangdao : Yanshan University , 2019 : 42 -43 .
[本文引用: 1]
中国老龄化的现状与积极应对
1
2017
... 目前,我国人口老龄化的发展趋势不容乐观,预计老年人口数量在2035年将突破3亿,2053年将达到峰值4.87亿,占届时全球老年人口的1/4[1 ] .随着人口老龄化的加剧,因机体老化、脑卒中等原因造成下肢运动功能障碍或残疾的患者日益增多.因此,研究下肢康复机器人对实现患者下肢运动功能恢复以及解决我国“患者多、康复医师少”的社会问题具有重要意义. ...
中国老龄化的现状与积极应对
1
2017
... 目前,我国人口老龄化的发展趋势不容乐观,预计老年人口数量在2035年将突破3亿,2053年将达到峰值4.87亿,占届时全球老年人口的1/4[1 ] .随着人口老龄化的加剧,因机体老化、脑卒中等原因造成下肢运动功能障碍或残疾的患者日益增多.因此,研究下肢康复机器人对实现患者下肢运动功能恢复以及解决我国“患者多、康复医师少”的社会问题具有重要意义. ...
Correction to: the FreeD module for the Lokomat facilitates a physiological movement pattern in healthy people: a proof of concept study
1
2019
... 国内外学者在康复机器人研究领域取得了丰硕的成果.瑞士Hocoma医疗器械公司与苏黎世Balgrist医学院联合研制了Lokomat步行康复训练机器人[2 ] ,它可以辅助患者模拟健康人的行走状态.饶玲军等[3 ] 设计的外骨骼康复机器人包括外骨骼本体和1副可实现身体重心控制和平衡的拐杖,其能使患者实现站立与行走.Aggogeri等[4 ] 研发的六自由度并联步态训练康复机器人Rutgers Ankle通过调节2个踏板的相对位置与姿态来协助患者进行平地行走、转弯行走以及平衡能力康复训练.Cudby等[5 ] 研制的可穿戴式下肢康复机器人REX可以帮助患者外出行走.Chen等[6 ] 研制了一款名为CUHK-EXO的下肢外骨骼康复机器人,它可以协助患者站立、坐下和行走等.上述穿戴式下肢康复机器人主要服务于康复后期的患者,能够辅助他们恢复行走能力,而对处于康复初期的患者存在一定局限性. ...
下肢外骨骼行走康复机器人研究与设计
1
2012
... 国内外学者在康复机器人研究领域取得了丰硕的成果.瑞士Hocoma医疗器械公司与苏黎世Balgrist医学院联合研制了Lokomat步行康复训练机器人[2 ] ,它可以辅助患者模拟健康人的行走状态.饶玲军等[3 ] 设计的外骨骼康复机器人包括外骨骼本体和1副可实现身体重心控制和平衡的拐杖,其能使患者实现站立与行走.Aggogeri等[4 ] 研发的六自由度并联步态训练康复机器人Rutgers Ankle通过调节2个踏板的相对位置与姿态来协助患者进行平地行走、转弯行走以及平衡能力康复训练.Cudby等[5 ] 研制的可穿戴式下肢康复机器人REX可以帮助患者外出行走.Chen等[6 ] 研制了一款名为CUHK-EXO的下肢外骨骼康复机器人,它可以协助患者站立、坐下和行走等.上述穿戴式下肢康复机器人主要服务于康复后期的患者,能够辅助他们恢复行走能力,而对处于康复初期的患者存在一定局限性. ...
下肢外骨骼行走康复机器人研究与设计
1
2012
... 国内外学者在康复机器人研究领域取得了丰硕的成果.瑞士Hocoma医疗器械公司与苏黎世Balgrist医学院联合研制了Lokomat步行康复训练机器人[2 ] ,它可以辅助患者模拟健康人的行走状态.饶玲军等[3 ] 设计的外骨骼康复机器人包括外骨骼本体和1副可实现身体重心控制和平衡的拐杖,其能使患者实现站立与行走.Aggogeri等[4 ] 研发的六自由度并联步态训练康复机器人Rutgers Ankle通过调节2个踏板的相对位置与姿态来协助患者进行平地行走、转弯行走以及平衡能力康复训练.Cudby等[5 ] 研制的可穿戴式下肢康复机器人REX可以帮助患者外出行走.Chen等[6 ] 研制了一款名为CUHK-EXO的下肢外骨骼康复机器人,它可以协助患者站立、坐下和行走等.上述穿戴式下肢康复机器人主要服务于康复后期的患者,能够辅助他们恢复行走能力,而对处于康复初期的患者存在一定局限性. ...
Functional design in rehabilitation: modular mechanisms for ankle complex
1
2016
... 国内外学者在康复机器人研究领域取得了丰硕的成果.瑞士Hocoma医疗器械公司与苏黎世Balgrist医学院联合研制了Lokomat步行康复训练机器人[2 ] ,它可以辅助患者模拟健康人的行走状态.饶玲军等[3 ] 设计的外骨骼康复机器人包括外骨骼本体和1副可实现身体重心控制和平衡的拐杖,其能使患者实现站立与行走.Aggogeri等[4 ] 研发的六自由度并联步态训练康复机器人Rutgers Ankle通过调节2个踏板的相对位置与姿态来协助患者进行平地行走、转弯行走以及平衡能力康复训练.Cudby等[5 ] 研制的可穿戴式下肢康复机器人REX可以帮助患者外出行走.Chen等[6 ] 研制了一款名为CUHK-EXO的下肢外骨骼康复机器人,它可以协助患者站立、坐下和行走等.上述穿戴式下肢康复机器人主要服务于康复后期的患者,能够辅助他们恢复行走能力,而对处于康复初期的患者存在一定局限性. ...
Liberty autonomy independence
1
2011
... 国内外学者在康复机器人研究领域取得了丰硕的成果.瑞士Hocoma医疗器械公司与苏黎世Balgrist医学院联合研制了Lokomat步行康复训练机器人[2 ] ,它可以辅助患者模拟健康人的行走状态.饶玲军等[3 ] 设计的外骨骼康复机器人包括外骨骼本体和1副可实现身体重心控制和平衡的拐杖,其能使患者实现站立与行走.Aggogeri等[4 ] 研发的六自由度并联步态训练康复机器人Rutgers Ankle通过调节2个踏板的相对位置与姿态来协助患者进行平地行走、转弯行走以及平衡能力康复训练.Cudby等[5 ] 研制的可穿戴式下肢康复机器人REX可以帮助患者外出行走.Chen等[6 ] 研制了一款名为CUHK-EXO的下肢外骨骼康复机器人,它可以协助患者站立、坐下和行走等.上述穿戴式下肢康复机器人主要服务于康复后期的患者,能够辅助他们恢复行走能力,而对处于康复初期的患者存在一定局限性. ...
A wearable exoskeleton suit for motion assistance to paralysed patients
1
2017
... 国内外学者在康复机器人研究领域取得了丰硕的成果.瑞士Hocoma医疗器械公司与苏黎世Balgrist医学院联合研制了Lokomat步行康复训练机器人[2 ] ,它可以辅助患者模拟健康人的行走状态.饶玲军等[3 ] 设计的外骨骼康复机器人包括外骨骼本体和1副可实现身体重心控制和平衡的拐杖,其能使患者实现站立与行走.Aggogeri等[4 ] 研发的六自由度并联步态训练康复机器人Rutgers Ankle通过调节2个踏板的相对位置与姿态来协助患者进行平地行走、转弯行走以及平衡能力康复训练.Cudby等[5 ] 研制的可穿戴式下肢康复机器人REX可以帮助患者外出行走.Chen等[6 ] 研制了一款名为CUHK-EXO的下肢外骨骼康复机器人,它可以协助患者站立、坐下和行走等.上述穿戴式下肢康复机器人主要服务于康复后期的患者,能够辅助他们恢复行走能力,而对处于康复初期的患者存在一定局限性. ...
A new concept of parallel robot for rehabilitation and fitness: the Lambda
1
2009
... 针对上述问题,部分学者提出了末端牵引式下肢康复机器人,其在辅助处于康复初期的患者恢复运动功能方面具有较大潜力.例如:Bouri等[7 ] 研制的牵引式下肢康复机器人Lambda采用左右对称的并联机械结构,能够以水平或垂直的姿态布置,适用于辅助康复初期患者进行矢状面内的康复训练;Monaco等[8 ] 研制的NEUROBike由2个相互垂直安装的导轨和踝部的转动装置组成,能够辅助卧床的患者模拟健康人的行走动作.相较于穿戴式下肢康复机器人,末端牵引式下肢康复机器人还具有轻量化与小型化的特点,且机械结构简单,易于控制,在训练过程中能够保证患者的安全性与舒适性.然而,上述学者研制的末端牵引式下肢康复机器人的训练方式较为单一,并不十分完善. ...
Design and evaluation of NEUROBike: a neurorehabilitative platform for bedridden post-stroke patients
1
2012
... 针对上述问题,部分学者提出了末端牵引式下肢康复机器人,其在辅助处于康复初期的患者恢复运动功能方面具有较大潜力.例如:Bouri等[7 ] 研制的牵引式下肢康复机器人Lambda采用左右对称的并联机械结构,能够以水平或垂直的姿态布置,适用于辅助康复初期患者进行矢状面内的康复训练;Monaco等[8 ] 研制的NEUROBike由2个相互垂直安装的导轨和踝部的转动装置组成,能够辅助卧床的患者模拟健康人的行走动作.相较于穿戴式下肢康复机器人,末端牵引式下肢康复机器人还具有轻量化与小型化的特点,且机械结构简单,易于控制,在训练过程中能够保证患者的安全性与舒适性.然而,上述学者研制的末端牵引式下肢康复机器人的训练方式较为单一,并不十分完善. ...
基于串并混联机构的四轮足步行机器人研究
1
2018
... 基于人体下肢运动等效模型,对末端牵引式下肢康复机器人的机构进行设计.由生物学理论可知,人体下肢运动是靠骨骼肌的收缩与舒张实现的,骨骼肌长度的改变可实现相邻2个骨骼的相对转动.故根据仿生学原理,在虎克铰链处布置2条UPS支链来模拟人体大腿的骨骼肌肉,在大腿与小腿之间加入1条RPR支链来模拟人体小腿的肌肉[9 ] ,且加入UPS支链和RPR支链可以提高机器人的承载能力,上述支链均以移动副作为驱动副,则以上部分组成(2-UP S&U)+(R&RP R)机构.然后,在(2-UP S&U)+(R&RP R)机构末端加入2条RR支链,通过伺服电机驱动齿轮箱内齿轮的传动实现左、右脚踏板的镜像开合,以辅助患者进行髋关节内外旋的康复训练;脚踏板跟随人体足部运动,不设置驱动.最终组成(2-UP S&U)+(R&RP R)+R R混联机构,如图1 (a)所示,对应的机器人三维模型如图1 (b)所示. ...
基于串并混联机构的四轮足步行机器人研究
1
2018
... 基于人体下肢运动等效模型,对末端牵引式下肢康复机器人的机构进行设计.由生物学理论可知,人体下肢运动是靠骨骼肌的收缩与舒张实现的,骨骼肌长度的改变可实现相邻2个骨骼的相对转动.故根据仿生学原理,在虎克铰链处布置2条UPS支链来模拟人体大腿的骨骼肌肉,在大腿与小腿之间加入1条RPR支链来模拟人体小腿的肌肉[9 ] ,且加入UPS支链和RPR支链可以提高机器人的承载能力,上述支链均以移动副作为驱动副,则以上部分组成(2-UP S&U)+(R&RP R)机构.然后,在(2-UP S&U)+(R&RP R)机构末端加入2条RR支链,通过伺服电机驱动齿轮箱内齿轮的传动实现左、右脚踏板的镜像开合,以辅助患者进行髋关节内外旋的康复训练;脚踏板跟随人体足部运动,不设置驱动.最终组成(2-UP S&U)+(R&RP R)+R R混联机构,如图1 (a)所示,对应的机器人三维模型如图1 (b)所示. ...
髋关节置换术后康复初期患者的躯体活动功能及其影响因素研究
1
2018
... 为了使所设计的牵引式下肢康复机器人可以更好地满足患者的康复训练需求,以及使患者在训练过程中更加舒适,对其机构参数进行优化.由于康复初期人体关节的活动度较小[10 ] ,故以康复中后期坐姿训练状态下人体下肢末端与机器人末端工作空间重合度(即机器人的安全工作空间[11 ] )为目标函数,以机器人机构参数为设计变量,利用MATLAB遗传算法工具箱[12 ] 进行优化计算. ...
髋关节置换术后康复初期患者的躯体活动功能及其影响因素研究
1
2018
... 为了使所设计的牵引式下肢康复机器人可以更好地满足患者的康复训练需求,以及使患者在训练过程中更加舒适,对其机构参数进行优化.由于康复初期人体关节的活动度较小[10 ] ,故以康复中后期坐姿训练状态下人体下肢末端与机器人末端工作空间重合度(即机器人的安全工作空间[11 ] )为目标函数,以机器人机构参数为设计变量,利用MATLAB遗传算法工具箱[12 ] 进行优化计算. ...
三维空间末端牵引上肢康复机器人轨迹规划研究
1
2021
... 为了使所设计的牵引式下肢康复机器人可以更好地满足患者的康复训练需求,以及使患者在训练过程中更加舒适,对其机构参数进行优化.由于康复初期人体关节的活动度较小[10 ] ,故以康复中后期坐姿训练状态下人体下肢末端与机器人末端工作空间重合度(即机器人的安全工作空间[11 ] )为目标函数,以机器人机构参数为设计变量,利用MATLAB遗传算法工具箱[12 ] 进行优化计算. ...
三维空间末端牵引上肢康复机器人轨迹规划研究
1
2021
... 为了使所设计的牵引式下肢康复机器人可以更好地满足患者的康复训练需求,以及使患者在训练过程中更加舒适,对其机构参数进行优化.由于康复初期人体关节的活动度较小[10 ] ,故以康复中后期坐姿训练状态下人体下肢末端与机器人末端工作空间重合度(即机器人的安全工作空间[11 ] )为目标函数,以机器人机构参数为设计变量,利用MATLAB遗传算法工具箱[12 ] 进行优化计算. ...
1
2014
... 为了使所设计的牵引式下肢康复机器人可以更好地满足患者的康复训练需求,以及使患者在训练过程中更加舒适,对其机构参数进行优化.由于康复初期人体关节的活动度较小[10 ] ,故以康复中后期坐姿训练状态下人体下肢末端与机器人末端工作空间重合度(即机器人的安全工作空间[11 ] )为目标函数,以机器人机构参数为设计变量,利用MATLAB遗传算法工具箱[12 ] 进行优化计算. ...
1
2014
... 为了使所设计的牵引式下肢康复机器人可以更好地满足患者的康复训练需求,以及使患者在训练过程中更加舒适,对其机构参数进行优化.由于康复初期人体关节的活动度较小[10 ] ,故以康复中后期坐姿训练状态下人体下肢末端与机器人末端工作空间重合度(即机器人的安全工作空间[11 ] )为目标函数,以机器人机构参数为设计变量,利用MATLAB遗传算法工具箱[12 ] 进行优化计算. ...
1
1988
... 参考《中国成年人人体尺寸》[13 ] ,以身高为1 754 mm、体重为75 kg的男性模型为参考,其对应的下肢参数如下:大腿长度为496 mm、小腿长度为396 mm,足部厚度为75 mm,胯部宽度为335 mm.故取xO 1 =167.5 mm,zO 1 =440 mm,l 6 =496 mm,l 7 =471 mm.根据人机工程学选取yO 1 =1 420 mm.规定坐姿状态下人体下肢各关节的康复训练角度,如表1 所示. ...
1
1988
... 参考《中国成年人人体尺寸》[13 ] ,以身高为1 754 mm、体重为75 kg的男性模型为参考,其对应的下肢参数如下:大腿长度为496 mm、小腿长度为396 mm,足部厚度为75 mm,胯部宽度为335 mm.故取xO 1 =167.5 mm,zO 1 =440 mm,l 6 =496 mm,l 7 =471 mm.根据人机工程学选取yO 1 =1 420 mm.规定坐姿状态下人体下肢各关节的康复训练角度,如表1 所示. ...
一种3-PRS并联机器人的姿态空间分析
1
2009
... 机构参数的变化会引起机器人工作空间的变化.首先,研究虎克铰链以及回转副的转角范围对机器人工作空间的影响,以“点集数”作为工作空间的评价指标[14 ] :点集数越大,工作空间越大. ...
一种3-PRS并联机器人的姿态空间分析
1
2009
... 机构参数的变化会引起机器人工作空间的变化.首先,研究虎克铰链以及回转副的转角范围对机器人工作空间的影响,以“点集数”作为工作空间的评价指标[14 ] :点集数越大,工作空间越大. ...
一种双机协同作业机器人的协作工作空间分析
1
2020
... 采用蒙特卡洛法[15 ] 求解优化后的人机系统工作空间,结果如图5 所示.其中:红色部分代表人体左下肢末端的工作空间,蓝色部分代表机器人左踏板的工作空间,两者相交的部分为机器人的安全工作空间.由此可得,矢状面内机器人的安全工作空间区域如图6 所示. ...
一种双机协同作业机器人的协作工作空间分析
1
2020
... 采用蒙特卡洛法[15 ] 求解优化后的人机系统工作空间,结果如图5 所示.其中:红色部分代表人体左下肢末端的工作空间,蓝色部分代表机器人左踏板的工作空间,两者相交的部分为机器人的安全工作空间.由此可得,矢状面内机器人的安全工作空间区域如图6 所示. ...
理疗师交互下的下肢康复训练机器人个性化步态规划方法
1
2018
... 由于下肢康复机器人以辅助患者进行矢状面内的康复训练为主[16 ] ,利用式(21) 求解人体下肢矢状面内的有效工作空间比ξ [17 ] : ...
理疗师交互下的下肢康复训练机器人个性化步态规划方法
1
2018
... 由于下肢康复机器人以辅助患者进行矢状面内的康复训练为主[16 ] ,利用式(21) 求解人体下肢矢状面内的有效工作空间比ξ [17 ] : ...
线性驱动并联机构基于有效工作空间比的尺度优化
1
2019
... 由于下肢康复机器人以辅助患者进行矢状面内的康复训练为主[16 ] ,利用式(21) 求解人体下肢矢状面内的有效工作空间比ξ [17 ] : ...
线性驱动并联机构基于有效工作空间比的尺度优化
1
2019
... 由于下肢康复机器人以辅助患者进行矢状面内的康复训练为主[16 ] ,利用式(21) 求解人体下肢矢状面内的有效工作空间比ξ [17 ] : ...
Research and design of a new horizontal lower limb rehabilitation training robot
3
2016
... 根据图6 和式(21) ,求得ξ =0.71,说明优化结果良好.此外,由于本文所研究的机器人为末端牵引式下肢康复机器人,其安全工作空间无法与人体下肢等效模型的工作空间完全一致,而根据康复护理相关理论[18 ] ,在图6 所示的人体末端与机器人末端工作空间的重合区域内,机器人可以帮助人体下肢进行多种康复训练,即机器人对人体有较好的适应性,可以满足其康复需求. ...
... 患者在康复初期身体虚弱,无法支撑其体重,处于这个时期的患者适合卧姿康复训练,通过保持正常的关节活动来恢复患者肢体的运动能力,以避免关节变形和肌肉萎缩.根据康复护理相关理论,采用CPM(continous passive motion,持续被动运动)方法[18 ] 进行康复训练.在康复训练过程中,小腿始终保持水平,髋关节与膝关节的角度变化一致. ...
... 为了确保运动的灵活性和舒适性,通常要求髋关节与膝关节的初始夹角为10 ° [18 ] ,规定人体髋、膝关节的运动角度为10 ° ~ 50 ° ,运动起始与终止处的速度与加速度均为0,故髋关节角度的五次插值多项式为: ...
一种新的康复与代步外骨骼机器人研究
1
2016
... 当患者处于康复中后期时,适合坐姿康复训练模式.圆周运动训练即人体矢状面内的康复训练,与脚踏车运动类似,能够帮助患者的髋、膝关节进行屈伸运动[19 ] .以工作空间中某一圆周轨迹为例,其圆心相对人机系统固定坐标系原点的位置(x 0 , y 0 , z 0 )=(-90, -665, 470) mm,半径R =115 mm,规定运动时间为10 s.令圆周运动的起始角度为360°,终止角度为0°,起始、终止处的速度与加速度均为0,由此得到圆心角位移的五次插值多项式: ...
一种新的康复与代步外骨骼机器人研究
1
2016
... 当患者处于康复中后期时,适合坐姿康复训练模式.圆周运动训练即人体矢状面内的康复训练,与脚踏车运动类似,能够帮助患者的髋、膝关节进行屈伸运动[19 ] .以工作空间中某一圆周轨迹为例,其圆心相对人机系统固定坐标系原点的位置(x 0 , y 0 , z 0 )=(-90, -665, 470) mm,半径R =115 mm,规定运动时间为10 s.令圆周运动的起始角度为360°,终止角度为0°,起始、终止处的速度与加速度均为0,由此得到圆心角位移的五次插值多项式: ...
空间四自由度下肢康复机器人设计与分析
1
2019
... 螺旋运动由直线运动与圆周运动组合形成,相较于圆周运动,其可以帮助患者实现空间内的康复训练[20 ] .定义在人机系统固定坐标系中,圆周分运动的圆心坐标为(x 0 , y 0 , z 0 )=(-200, -630, 600) mm,半径R =70 mm;令圆周分运动的起始角度为720°,终止角度为0°.直线分运动的起始位置为(-200, -630, 600) mm,终止位置为(0, -630, 600) mm;2个分运动的起始与终止处的速度、加速度均为0.规定螺旋运动时间为15 s.由此可得,螺旋运动的2个分运动的五次插值多项式为: ...
空间四自由度下肢康复机器人设计与分析
1
2019
... 螺旋运动由直线运动与圆周运动组合形成,相较于圆周运动,其可以帮助患者实现空间内的康复训练[20 ] .定义在人机系统固定坐标系中,圆周分运动的圆心坐标为(x 0 , y 0 , z 0 )=(-200, -630, 600) mm,半径R =70 mm;令圆周分运动的起始角度为720°,终止角度为0°.直线分运动的起始位置为(-200, -630, 600) mm,终止位置为(0, -630, 600) mm;2个分运动的起始与终止处的速度、加速度均为0.规定螺旋运动时间为15 s.由此可得,螺旋运动的2个分运动的五次插值多项式为: ...