人体下肢外骨骼康复机器人的仿真分析
2022-10-18
来源:好走旅游网
第34卷第6期 陕西科技大学学报 Journal of Shaanxi University of Science&Technology Vo1.34 No.6 Dec.2O16 2016年12月 文章编号:1000.5811(2016)06—0154—06 人体下肢外骨骼康复机器人的仿真分析 夏 田,桓 茜,陈 宇,徐建林 (陕西科技大学机电工程学院,陕西西安710021) 摘 要:采用电机驱动方式,设计出一种用于辅助老年人和患者生活的服务型机器人——人体 下肢外骨骼康复机器人,通过建立人体行走姿态时的【)IH数学模型,推导出了髋关节、膝关节 和踝关节行走姿态中的坐标方程.在ADAMS环境下进行人机耦合的仿真,结果表明:下肢外 骨骼各关节在空间中具有连续的运动轨迹,满足各关节运动的可达性,符合人体下肢的运动特 性.搭建样机试验平台,测取了三个关节角度随周期的变化曲线,验证了下肢外骨骼装置的运 动轨迹,并且试验运动曲线和仿真曲线基本一致. 关键词:下肢外骨骼;运动学;D—H模型;耦合仿真;ADAMS 中图分类号:TP242 文献标志码:A Simulation analysis of human lower limb exoskeleton rehabiljtati0n robot XIA Tian,HUAN Xi,CHEN Yu,XU Jian—lin (College of Mechanical and Electrical Engineering,Shaanxi University of Science 8L Technology,Xi an 710021,China) Abstract:The kind of device of human lower limb exoskeleton rehabilitation robot iS a kind of service robot,which is designed to be used for assisting the life of old people and sufferer, based on the motor drive.By establishing the D—H mathematical model of human walking posture,the coordinate equations of hip j oint,knee j oint and ankle j oint are derived.Under the ADAMS environment,the simulation results show that the joints of the lower limb exo— skeleton have a continuous traj ectory in space,which can meet the motility of the j oint movement,and meet the motion characteristics of the human lower 1imb.A prototype test platform is set up to test the angle curve which changes with the cycle.The motion traj ectory of the 1ower limb exoskeleton iS verified.and the experimental curves and the simulation curves are basically consistent. Key words:lower limb exoskeleton;kinematics;D—H model;coupling simulation;ADAMS 收稿日期:2016-10—22 基金项目:陕西省科技厅科技计划项目(2013K07—08) 作者简介:夏田(1962一),女,陕西咸阳人,教授,研究方向:数控装备、机器臂及AVG小车 第6期 夏 田等:人体下肢外骨骼康复机器人的仿真分析 ・ 155 ・ 尼龙;各关节连接构件采用高强度不锈钢材料.各 0 引言 个构件的不同连接和配合方式,可以满足各关节的 自由度,各构件的连接关系如表2所示. 随着全球人类老龄化的到来,老年人的生活运 动障碍已经成为社会运作的一个巨大压力,老年人 正常行走及上下楼梯等运动都需要人帮扶,给家人 造成不小的生活负担和压力,另外,在医疗领域,下 肢瘫痪或运动功能障碍患者进行康复训练也需要 机械装置提供支撑、保护和助力 ]. 在国外,斯坦福大学在2000年推出了手臂康 复训练机器人样机,并且己打入了市场,同时日本 也研制了旋转式康复机器人,柏林大学研制出了康 复机器人腿部样机;在国内,清华大学在国内率先 研制了卧式下肢康复训练机器人样机,哈尔滨工程 大学的下肢康复训练机器人是通过步态和姿态机 器人控制系统模拟正常人行走的步态、踩关节的运 动姿态,浙江大学流体与传动控制国家重点实验室 采用气动系统设计了一种新型的可穿戴式下肢步 行外骨骼机构 ]. 人体下肢外骨骼康复机器人是一种可穿戴的 仿人型动力系统装置,其将人体和外骨骼康复装置 高度耦合在一起,辅助人体运动,这种人机耦合交 互系统,需要各个构件与身体的运动姿态保持一 致,即外骨骼自由度的形式和构成需要和人体基本 吻合 .通过对下肢外骨骼机械结构的设计,并建 立了下肢外骨骼的运动学方程,用ADAMS进行 仿真分析,搭建试验平台进行运动分析,为样机制 作提供一定的理论支撑. 1 下肢外骨骼机械结构 根据GB/T 10O0—1988《中国成年人人体尺 寸》提供的立姿人体尺寸基础数据 ],按照人体身 高175 mm的比例,将下肢外骨骼的基础尺寸设计 为:髋部尺寸z 一334 mm,大腿杆尺寸z :428 mm,小腿杆尺寸 3=430 mm,脚部尺寸z4—68 mm.人体下肢各关节的自由度如表I所示. 表I下肢各关节自由度 图1 下肢外骨骼机械结构 表2各构件连接关系 名称及序号 背架i 支撑架2 连接方式 与2固定连接 与i固连,与3铰接 髋部电机3 与2固连.与5减速器联接 分别与2、5铰接 与3减速器联接.与6螺纹联接 分别与5、7螺纹联接 分别与8、9铰接,与6螺纹联接 髋部弹簧阻尼4 髋部转动轴5 大腿高度调节6 大腿连杆7 膝部电机8 小腿连杆9 小腿高度调节Io 踝部转动轴Il 与7固连,与9减速器联接 与7铰接,与Io螺纹联接 分别与9、l1螺纹联接 与Io螺纹联接,与I2铰接 分别与I1、14铰接 与I4减速器联接 与i3减速器联接,与I2铰接 如图1所示,采用电机驱动方式,设计了人体 下肢外骨骼康复机器人,各连接部件采用铝合金材 料;背部、腰部和腿部的束紧带为柔性构件,材料为 弹簧阻尼I2 脚部电机13 脚部构件I4 ・ 156 ・ 陕西科技大学学报 第34卷 此外,各关节处的电机驱动系统,可为外骨骼 提供动力以辅助穿戴者运动,大小腿高度调节装置 通过螺纹进行调节,以适应不同的穿戴者使用,弹 簧阻尼装置具有缓冲和限位的功能,起到储能和助 力的作用. 2下肢外骨骼的运动学建模 2.1 I)-H模型 在机器人运动学研究中,I)_H模型是最常用 的运动学模型,从结构设计角度可知,下肢外骨骼 是由多个刚体串联构成,将各组成构件简化为连 杆,将多自由度关节分解成几个单自由度关节,对 下肢外骨骼进行运动学建模分析[7]. 图2所示为人体直线行走姿态时单侧腿的I)_ H数学模型.其中,0 、O 、O。分别为髋关节、膝关 节和踝关节的坐标原点,O。一x。Y。Z。是建立在髋 部和躯干之间的绝对坐标系,各连杆的动坐标系为 O1一X1YlZ1、O2一x2Y2Z2和O3一X3Y3Z3,0f为轴X 1嘲 o oa S 绕轴Z 旋转变换到轴x 的旋转角,即0 、0 、0。, 各个连杆长度为l 、z 、z。. 咖o o 图2 D-H模型简图 根据下肢外骨骼的运动情况,I)_H模型中变 化矩阵的运动学参数如表3所示. 表3 H模型参数 根据表3中的I)_H运动学参数,可求解出躯 干和大腿杆、大腿杆和小腿杆、小腿杆和脚构件三 个齐次坐标变换矩阵,即髋关节、膝关节和踝关节 三个旋转关节的坐标变换,见公式(1)、公式(2)和 公式(3). 0 Z1 cos01 0 Z1 sin01 Mo1 1 O 0 1 0 Z2COS02 0 Z2 sin02 M12= .C S o o嘲 o o ..C S (2) 1 O O 1 0 Z3 cos03 c藿}0 o 咖0 o 0 0 Z3 sin03 M2。= (3) 1 O 0 1 2.2运动学正问题 外骨骼机器人是已知关节的角度变量与时间 的函数关系,然后求解各活动杆件的位姿关系,采 用运动学的正问题求解. 按照齐次变换矩阵的右乘法则,分别求得小腿 连杆和脚部构件的齐次坐标变换矩阵,见公式(4) 和公式(5). Mo2一Mo1×M12 (4) Mo3一Mol×M12×M23 (5) 齐次变换矩阵建立两个相邻连杆的位姿变换 关系,由此推导出髋关节、膝关节和踝关节直线行 走姿态中的坐标方程分别为: =Mo1 (6) 一Mo2 (7) 第6期 夏 田等:人体下肢外骨骼康复机器人的仿真分析 ・ 157 ・ 3 3 :Mo3 (8) 0 1 3人机耦合仿真分析 人体与外骨骼康复装置的耦合性,影响到装置 可穿戴性和运动关联性,如果人机耦合设计合理得 当,可使身体各个部位感觉舒适且没有约束感,而 且装置结构紧凑,运动协调,与外界干涉可能性小, 通过Solidworks建立人体和外骨骼康复装置的耦 合虚拟模型,导入ADAMS中对其进行仿真分析. 3.1人机耦合系统模型 3.1.1 人体模型 按照GB/T 1000—1988提供的人体站立时姿 态的身体各部位相对于身高的比例关系,以身高为 1 75 mm的中国成年人为例,其身体各部位的尺寸 比例如表4所示__6].通过实际测量成年人的围度、 宽度和高度尺寸,在Solidworks中建立人体的近 似模型 ],如图3所示. 表4身体各部位比例关系 图3耦合模型 3.2 ADAMS环境下耦合仿真分析 将人机耦合模型导人到ADAMS中,如图4 所示,包括人体模型,下肢外骨骼模型和地面模型, 进行相关的设置建立一个近似实际的仿真系统,其 主要分析过程如下: (1)设置材料属性 由于人体靠骨架成型和支撑,可以将骨头作为 人体模型的成分: ,在此将人体模型的材料属性设 置为杨氏模量值为1.3×10 MPa,泊松比为 0.031,密度为1.057×10。Kg/m。. 外骨骼模型的刚性构件选用铝合金材料,杨氏 3.1.2下肢外骨骼仿真模型 为了便于动态仿真分析,可以对系统进行简化 处理为一个负重且满足各部位相对运动关系的仿 真模型,如图3所示.其中背负装置1,支撑架2,大 模量为7.17×10 MPa,泊松比为0.33,密度为 2.71X10。Kg/m。.柔性构件选用尼龙材料,杨氏 模量为2.62×10 Mpa,泊松比为0.34,密度为 1.12×10。Kg/m . 腿杆件3,小腿杆件4和脚部5为刚性构件,主要 起支撑和传递运动的作用;背部束紧带6、腰部束 紧带7和腿部束紧带8为柔性构件,主要起连接作 用.为了降低仿真分析的难度,只在髋关节、膝关节 和踝关节处分别设置一个屈伸运动的自由度. 实际中,地面为混凝土,由于脚部材料为铝合 金,为了不影响仿真的结果.将地板的材料设置为 橡胶.其杨氏模量为6.11 MPa,泊松比为0.49,密 度为6.7×10 Kg/m . (2)添加约束 为了保证人机耦合的正确性,对模型添加合适 的约束是非常关键的 . ・ 158 ・ 陕西科技大学学报 下肢的运动特性. 1.2 1.0 第34卷 ①在人体模型中,将各关节设置为旋转副,将 人体模型和各束紧带之间设置为固定副,将人体躯 干部分设置为移动副; g 0.8 \ ②在外骨骼模型中,将各关节设置为旋转副, 将刚性构件和柔性构件之间设置为固定副. (3)添加接触力 专0.6 。o 暑0 4 0.2 0 0 人行走过程中主要产生重力、支撑力和摩擦 力,将整体模型添加竖直向下重力加速度;在行走 200 0 150 100 50 / 、. ,, 、、 、 过程中.脚部构件是唯一与地面接触的部分,在脚 部构件和地面之间设置为接触力. (4)添加关节驱动 g 0 -, 、\/ .一一一一 、、、、 \\ 、,50 / / 。 、 , 一100 200 踝关节 ’、、\ 一 / 鼍.150 .膝 ①将髋关节、膝关节和踝关节的关节运动数据 保存为.txt格式,导入ADAMS中,生成6个关节 的Spline曲线函数¨ll; ②在6个关节的旋转副处设置旋转驱动,定义 为AKISPI (time.0,SPI 1NEj,0)等. 《 30 2O 10 0 .10 .20 30 40 .(5)仿真及结果分析 ..正常人行走的一个周期是1.5 S,设置仿真时 间为1.5 S。为了便于观察将仿真步长设置为 2 0O0,人机耦合系统的仿真动画如图4所示,可知 50 (c)角度曲线 图5耦合模型的仿真结果 4样机试验平台 人机耦合模型行走的运动轨迹和实际人体行走基 本一致,证明了仿真的正确性和合理性. 搭建样机试验平台是下肢外骨骼康复机器人 诸多试验中的一种,其目的是步态预编程控制(被 动控制)下外骨骼装置的运动轨迹验证_】 .如图6 所示,试验的硬件组成部分包括下肢外骨骼装置、 步进电机、驱动器、驱动电源、信号发生器和控制器 等. 图4 仿真动画结果 以下肢外骨骼中的右侧下肢为例,在后处理模 块中进行数据处理,选择各关节的位移、角速度和 角度随着运动周期的关系曲线,滤波处理掉跳跃式 信号,如图5所示. 可知,下肢外骨骼中髋关节、膝关节和踝关节 图6样机行走试验 用图像分析法对行走试验进行分析,拍摄样机 试验画面,然后绘制样机各关节运行的轨迹,用关 键点连接法求解出试验结果如图7所示,测取了三 个关节角度随周期的变化曲线. 由图7可知,试验运动曲线的平滑性较差,波 动性较大,局部变形突出,由于试验过程中电机的 在空间中具有连续的运动轨迹。可以满足运动的可 达性.同时,在初始阶段,脚部构件位移最 .大腿 构件位移最小,并在运动过程中.大腿和脚部的位 移交替变化,小腿的位移处于中间状态,符合人体 振动所造成.对比仿真曲线图5(c),试验运动曲线 第6期 夏 田等:人体下肢外骨骼康复机器人的仿真分析 ・ 159 ・ 30 20 \ —0 . 10 =cP 0 、 / 襁.10 、 、' .} 20 \ .30 0 0.3 0.6 0.9 1.2 1.5 时间/S (a)髋关节角度变化图 厂\ g-- \ f ’; ±仁 } 、 壁 八 { { . . 时IWs (b)膝关节角度变化图 10 。g--0 、 / \ 接 拉一10 V 、 ;7 f :H< 髅 .20 f { } 时I司/S (c)踝关节角度变化图 图7样机试验结果 的轨迹和仿真曲线走向一致,这也在一定程度上验 证了虚拟耦合仿真的正确合理,但由于方向设置问 题,试验曲线与仿真曲线的方向相反. 此外,试验中未在外骨骼关节及末端安装角度 和位移传感器,而且步进电机没有反馈系统,以后 的试验需要搭建反馈系统,并且考虑电机的动态指 标——共振点,避免电机作为激振源引起外骨骼的 共振. 5 结论 (1)采用电机驱动方式,设计了人体下肢外骨 骼康复机器人装置,并且设计有缓冲和限位功能, 大小腿的长度均可调节. (2)建立了人体行走姿态时的D—H数学模型, 推导出髋关节、膝关节和踝关节行走姿态中的坐标 方程. (3)在ADAMS环境下进行人机耦合的仿真, 结果表明下肢外骨骼各关节在空间中具有连续的 运动轨迹,满足运动的可达性,符合人体下肢的运 动特性,搭建样机试验平台验证了下肢外骨骼装置 的运动轨迹. 参考文献 [1]Meng W,Liu Q,Zhou Z D,et a1.Recent development of mechanisms and control strategies for robot—assisted lower limb rehabilitation[J].Mechatronics,2015,31(4):132— 145. [2]Viteckova S,Kutilek P,Jirina M.Wearable lower limb ro— botics:A review[J'].Biocybernetics and Biomedical Engi— neering,2013,33(2):96—105. [3]Huang C E,Chen J S.On the implementation and control of a pneumatic power active lower—limb orthosis[J]. Mechatronics,2013,23(5):505—517. [4]Galle S,Malcolm P,Derave W,et a1.Adaptation to walk— ing with an exoskeleton that assists ankle extension[J]. Gait&Posture,2013,38(3):495—499. E5]Li Z Q,Xie H Q.Proceeding of human exoskeleton tech— nology and discussions on future research[J].Chinese Journal of Mechanical Engineering,2014,27(3):437—442. [6]孟庆鑫.机器人技术基础[M].哈尔滨:哈尔滨工业大学出 版社,2006. [7]李团结.机器人技术[M].北京:电子工业出版社,2009. [8]王 阳,宋道志,王晓光.下肢外骨骼机器人人机腿部约束 分析[J].现代制造技术与装备,2016(1):51—52. [9]Robert B.Exoskeletons and robotic prosthetics:A review of recent developments[J ̄.Industrial Robot,2009,36(5): 421—427. [1O]高东强,王 伟.摆线钢球行星传动机构设计与仿真分析 [J].陕西科技大学学报(自然科学版),2014,32(6):139— 144. [¨]杨宗林,曾亿山,王善杰.下肢外骨骼机器人结构设计和 动力学仿真[J].机械设计与制造,2016(1):91-97. [12]Arnold E,Ward S,Lieber R.A model of the lower limb for analysis of human movement[J].Annals of Biomed— ical Engineering,2010,38(2):269—279. 【责任编辑:陈佳】