第10卷第4期 智能系统学报 Vol.10 No.4 2015年8月 CAAI Transactions on Intelligent Systems Aug.2015 D0:10.3969/j.issn.1673-4785.201503039 网络出版地址：http://www.cnki.net/kcms/detail/23.1538.tp.20150702.1104.001.html 人体下肢生物力学建模研究进展 邵明旭，王斐，般腾龙，刘健 (东北大学信息科学与工程学院，辽宁沈阳110819) 摘要：人体下肢生物力学建模与仿真是穿戴式外骨骼机器人系统开发的一个重要内容。对其研究所获得相关的 理论与技术方法对生物动力学、康复医学、假肢及运动康复器械设计等领域的发展具有促进作用。本文以人体下肢 生物力系统为研究对象，概括和总结了国内外下肢建模与仿真技术的研究现状，就目前普遍采用的基于Lagrange方 程和角动量定理的多刚体模型法、仿真软件建模法、H三元素法、黑箱训练等方法进行了详细的分析，并对研究趋 势进行了展望。本文所综述的动力学建模与仿真验证方法对实现穿戴式外骨骼机器人和谐自然人机交互设计具有 重要的指导意义。 关键词：下肢生物力学建模；Lagrange方程；角动量定理；Hill模型；仿真软件建模；黑箱训练法 中图分类号：TP391文献标志码：A文章编号：1673-4785(2015)04-0518-10 中文引用格式：邵明旭，王斐，殷腾龙，等.人体下肢生物力学建模研究进展[J].智能系统学报，2015,10(4)：518-527. 英文引用格式：SHAO Mingxu,WANG Fei,YIN Tenglong,etal.Progress on the research of human lower limb biomechanical modeling[J].CAAI Transactions on Intelligent Systems,2015,10(4):518-527. Research progress on the human lower limb biomechanical modeling SHAO Mingxu,WANG Fei,YIN Tenglong,LIU Jian (College of Information Science and Engineering,Northeastern University,Shenyang 110819,China) Abstract:The research on the biomechanical modeling and simulation of human lower limbs is an important content in the development of wearable exoskeleton robots.Theoretical and technical methods derived from this research can promote the process of biomechanics,rehabilitation medicine and prosthetic/orthotic devices.This work reviews the state-of-the-art techniques for modeling and simulating biomechanics of human lower limbs and makes analysis of popular methods,such as multi-body modeling,simulation software modeling,Hill three elements modeling and black box training modeling based on Lagrange equation and theorem of angular momentum.The future prospects in this research field are also provided in this paper.The biomechanical modeling and simulating methods discussed is of great significance to the design of naturally harmonious human-robot interaction of wearable exoskeleton robots. Keywords:biomechanical modeling of human lower limb;Lagrange equation;theorem of angular momentum;Hill model;simulation software modeling;black box training modeling 随着MEMS、材料、控制和计算机技术的飞速发围内引起了越来越多学者的关注。但目前的动力外 展，以及我国快速进入老龄化社会后对康复辅助设 骨骼和假肢研究存在的一个共性问题是对人体感觉 备的旺盛需求，穿戴式外骨骼机器人技术在世界范 运动控制、生物力学等知识了解不够，这就导致了 虽然设计出了先进的机器人硬件和软件，但由于没 收稿日期：2015-03-24.网络出版日期：2015-07-02 有考虑人如何作用，因此设备在与人交互和协同方 基金项目：中央高校基础科研业务费资助项目(120124002)：辽宁省自 面效果欠佳。由于人体结构和功能的特殊性和复杂 然科学基金资助项目(2013020040). 通信作者：王斐.E-mail:wangfei@(ise.ncu.cd.cn 性，建立一个能完全反映人体状态的生物力学模型
第 １０ 卷第 ４ 期 智 能 系 统 学 报 Ｖｏｌ．１０ №．４ ２０１５ 年 ８ 月 ＣＡＡＩ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｓｙｓｔｅｍｓ Ａｕｇ． ２０１５ ＤＯＩ：１０．３９６９ ／ ｊ．ｉｓｓｎ．１６７３⁃４７８５．２０１５０３０３９ 网络出版地址：ｈｔｔｐ： ／ ／ ｗｗｗ．ｃｎｋｉ．ｎｅｔ ／ ｋｃｍｓ／ ｄｅｔａｉｌ ／ ２３．１５３８．ｔｐ．２０１５０７０２．１１０４．００１．ｈｔｍｌ 人体下肢生物力学建模研究进展 邵明旭，王斐，殷腾龙，刘健 （东北大学 信息科学与工程学院，辽宁 沈阳 １１０８１９） 摘 要：人体下肢生物力学建模与仿真是穿戴式外骨骼机器人系统开发的一个重要内容。 对其研究所获得相关的 理论与技术方法对生物动力学、康复医学、假肢及运动康复器械设计等领域的发展具有促进作用。 本文以人体下肢 生物力系统为研究对象，概括和总结了国内外下肢建模与仿真技术的研究现状，就目前普遍采用的基于 Ｌａｇｒａｎｇｅ 方 程和角动量定理的多刚体模型法、仿真软件建模法、Ｈｉｌｌ 三元素法、黑箱训练等方法进行了详细的分析，并对研究趋 势进行了展望。 本文所综述的动力学建模与仿真验证方法对实现穿戴式外骨骼机器人和谐自然人机交互设计具有 重要的指导意义。 关键词：下肢生物力学建模；Ｌａｇｒａｎｇｅ 方程；角动量定理；Ｈｉｌｌ 模型；仿真软件建模；黑箱训练法 中图分类号： ＴＰ３９１ 文献标志码：Ａ 文章编号：１６７３⁃４７８５（２０１５）０４⁃０５１８⁃１０ 中文引用格式：邵明旭，王斐，殷腾龙，等． 人体下肢生物力学建模研究进展［Ｊ］． 智能系统学报， ２０１５， １０（４）： ５１８⁃５２７． 英文引用格式：ＳＨＡＯ Ｍｉｎｇｘｕ， ＷＡＮＧ Ｆｅｉ， ＹＩＮ Ｔｅｎｇｌｏｎｇ， ｅｔ ａｌ． Ｐｒｏｇｒｅｓｓ ｏｎ ｔｈｅ ｒｅｓｅａｒｃｈ ｏｆ ｈｕｍａｎ ｌｏｗｅｒ ｌｉｍｂ ｂｉｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｉｎｇ［Ｊ］． ＣＡＡＩ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｓｙｓｔｅｍｓ， ２０１５， １０（４）： ５１８⁃５２７． Ｒｅｓｅａｒｃｈ ｐｒｏｇｒｅｓｓ ｏｎ ｔｈｅ ｈｕｍａｎ ｌｏｗｅｒ ｌｉｍｂ ｂｉｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｉｎｇ ＳＨＡＯ Ｍｉｎｇｘｕ， ＷＡＮＧ Ｆｅｉ， ＹＩＮ Ｔｅｎｇｌｏｎｇ， ＬＩＵ Ｊｉａｎ （Ｃｏｌｌｅｇｅ ｏｆ Ｉｎｆｏｒｍａｔｉｏｎ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ， Ｎｏｒｔｈｅａｓｔｅｒｎ Ｕｎｉｖｅｒｓｉｔｙ， Ｓｈｅｎｙａｎｇ １１０８１９， Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｅ ｒｅｓｅａｒｃｈ ｏｎ ｔｈｅ ｂｉｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｉｎｇ ａｎｄ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｈｕｍａｎ ｌｏｗｅｒ ｌｉｍｂｓ ｉｓ ａｎ ｉｍｐｏｒｔａｎｔ ｃｏｎｔｅｎｔ ｉｎ ｔｈｅ ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ ｗｅａｒａｂｌｅ ｅｘｏｓｋｅｌｅｔｏｎ ｒｏｂｏｔｓ． Ｔｈｅｏｒｅｔｉｃａｌ ａｎｄ ｔｅｃｈｎｉｃａｌ ｍｅｔｈｏｄｓ ｄｅｒｉｖｅｄ ｆｒｏｍ ｔｈｉｓ ｒｅｓｅａｒｃｈ ｃａｎ ｐｒｏｍｏｔｅ ｔｈｅ ｐｒｏｃｅｓｓ ｏｆ ｂｉｏｍｅｃｈａｎｉｃｓ， ｒｅｈａｂｉｌｉｔａｔｉｏｎ ｍｅｄｉｃｉｎｅ ａｎｄ ｐｒｏｓｔｈｅｔｉｃ ／ ｏｒｔｈｏｔｉｃ ｄｅｖｉｃｅｓ． Ｔｈｉｓ ｗｏｒｋ ｒｅｖｉｅｗｓ ｔｈｅ ｓｔａｔｅ⁃ｏｆ⁃ｔｈｅ⁃ａｒｔ ｔｅｃｈｎｉｑｕｅｓ ｆｏｒ ｍｏｄｅｌｉｎｇ ａｎｄ ｓｉｍｕｌａｔｉｎｇ ｂｉｏｍｅｃｈａｎｉｃｓ ｏｆ ｈｕｍａｎ ｌｏｗｅｒ ｌｉｍｂｓ ａｎｄ ｍａｋｅｓ ａｎａｌｙｓｉｓ ｏｆ ｐｏｐｕｌａｒ ｍｅｔｈｏｄｓ， ｓｕｃｈ ａｓ ｍｕｌｔｉ⁃ｂｏｄｙ ｍｏｄｅｌｉｎｇ ， ｓｉｍｕｌａｔｉｏｎ ｓｏｆｔｗａｒｅ ｍｏｄｅｌｉｎｇ， Ｈｉｌｌ ｔｈｒｅｅ ｅｌｅｍｅｎｔｓ ｍｏｄｅｌｉｎｇ ａｎｄ ｂｌａｃｋ ｂｏｘ ｔｒａｉｎｉｎｇ ｍｏｄｅｌｉｎｇ ｂａｓｅｄ ｏｎ Ｌａｇｒａｎｇｅ ｅｑｕａｔｉｏｎ ａｎｄ ｔｈｅｏｒｅｍ ｏｆ ａｎｇｕｌａｒ ｍｏｍｅｎｔｕｍ． Ｔｈｅ ｆｕｔｕｒｅ ｐｒｏｓｐｅｃｔｓ ｉｎ ｔｈｉｓ ｒｅｓｅａｒｃｈ ｆｉｅｌｄ ａｒｅ ａｌｓｏ ｐｒｏｖｉｄｅｄ ｉｎ ｔｈｉｓ ｐａｐｅｒ． Ｔｈｅ ｂｉｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｉｎｇ ａｎｄ ｓｉｍｕｌａｔｉｎｇ ｍｅｔｈｏｄｓ ｄｉｓｃｕｓｓｅｄ ｉｓ ｏｆ ｇｒｅａｔ ｓｉｇｎｉｆｉｃａｎｃｅ ｔｏ ｔｈｅ ｄｅｓｉｇｎ ｏｆ ｎａｔｕｒａｌｌｙ ｈａｒｍｏｎｉｏｕｓ ｈｕｍａｎ⁃ｒｏｂｏｔ ｉｎｔｅｒａｃｔｉｏｎ ｏｆ ｗｅａｒａｂｌｅ ｅｘｏｓｋｅｌｅｔｏｎ ｒｏｂｏｔｓ． Ｋｅｙｗｏｒｄｓ：ｂｉｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｉｎｇ ｏｆ ｈｕｍａｎ ｌｏｗｅｒ ｌｉｍｂ； Ｌａｇｒａｎｇｅ ｅｑｕａｔｉｏｎ； ｔｈｅｏｒｅｍ ｏｆ ａｎｇｕｌａｒ ｍｏｍｅｎｔｕｍ； Ｈｉｌｌ ｍｏｄｅｌ； ｓｉｍｕｌａｔｉｏｎ ｓｏｆｔｗａｒｅ ｍｏｄｅｌｉｎｇ； ｂｌａｃｋ ｂｏｘ ｔｒａｉｎｉｎｇ ｍｏｄｅｌｉｎｇ 收稿日期：２０１５⁃０３⁃２４． 网络出版日期：２０１５⁃０７⁃０２． 基金项目：中央高校基础科研业务费资助项目（ １２０１２４００２）；辽宁省自 然科学基金资助项目（２０１３０２００４０）． 通信作者：王斐． Ｅ⁃ｍａｉｌ： ｗａｎｇｆｅｉ＠ ｉｓｅ．ｎｅｕ．ｅｄｕ．ｃｎ． 随着 ＭＥＭＳ、材料、控制和计算机技术的飞速发 展，以及我国快速进入老龄化社会后对康复辅助设 备的旺盛需求，穿戴式外骨骼机器人技术在世界范 围内引起了越来越多学者的关注。 但目前的动力外 骨骼和假肢研究存在的一个共性问题是对人体感觉 运动控 制、生物力学等知识了解不够，这就导致了 虽然设计出了先进的机器人硬件和软件，但由于没 有考虑人如何作用，因此设备在与人交互和协同方 面效果欠佳。 由于人体结构和功能的特殊性和复杂 性，建立一个能完全反映人体状态的生物力学模型
第4期 邵明旭，等：人体下肢生物力学建模研究进展 ·519 非常困难。为此，根据研究的目的，对人体或人体的 医学伦理性问题，肌肉的生物动力学参数测量具有 一部分作适当的简化，在一些必要的假设条件下，从 一定困难，现有的大部分下肢建模方式均未考虑人 不同角度出发，建立人体下肢的生物动力学模型有 体下肢的生物模型，只考虑下肢的典型结构，或是建 着巨大的意义山。 立简化后的下肢多刚体模型9】，或是用3D建模技 本文旨在对国内外人体下肢建模和仿真的最新 术模拟简易的人体下肢[1013)，或是采用黑箱训练等 研究状况进行综述，并对今后的研究趋势进行阐述 方法跨过人体结构[)，进而通过动力学分析方法描 分析，为进一步深入研究下肢生物动力学建模提供 述模型，并对模型进行求解。下面分别介绍几种典 一定的参考依据。 型的下肢建模方法。 1人体下肢结构概述 3.1 Lagrange多刚体模型 国内外的许多学者针对人体不同的运动或动作 从解剖学的角度看，人的躯体是由骨、骨连接和 情形，考虑不同情况下模型的具体要求，建立出了相 骨骼肌组成，全身骨骨间借关节连成骨骼，构成主体 应的多刚体模型，其中最简单的方法是用单质点模 人体支架。在运动过程中，骨、骨连接主要起杠杆作 拟人体运动，不过模型过于简单，实际价值较低。 用，而具体的肢体关节是运动的枢纽，骨骼肌是运动 Hanavan于1964年提出一个由15个多刚体环节组 的动力来源和主要执行部件。 成的，通过球铰系统连接的人体力学模型，这种模型 在人体下肢运动中起主要参与的骨骼如下：髋 将人体主要的生物结构简化为几何形状的均质刚 骨、股骨、腓骨、胫骨等。骨与骨之间主要由结缔组织 体，建立并计算出基于人体惯性参数的回归方程，是 连接，共同形成人体骨架，骨头之间连接的地方称为 一套比较完整的具有个性化参数特征的多刚体人体 关节。而在人体下肢运动中起重要作用的关节主要 模型，这个模型也在一段时间内被连续引用)。波 有3个：髋关节、膝关节和踝关节。髋关节为多轴性 兰学者Blajera等为了研究杂技演员在蹦床上腾空 关节，可做屈、伸、内收、外展、内旋、外旋和旋转运动， 翻跟斗的动作过程，建立了一个自由度为10的多刚 由于髋关节起到连接上下肢和承重的作用，其骨密度 体人体模型io。Anderson等在研究人垂直起跳直 较大，故而稳固性较好。膝关节是人体负重比例最 到达到最大高度过程中肌肉运动协调模式时，把人 大、结构最复杂的关节，其主要负责大腿和小腿的屈、 体简化为10段、具有54块肌肉的模型，并在此模型 伸运动，和某些特定位置下的小幅度旋转运动。踝关 的基础上进一步建立人在行走过程中行走距离和步 节主要负责脚掌的背屈、伸、内收和外展运动。 态与能量消耗间的关系[)。Pejhan等在基于步态 与此同时，对人体下肢运动充当动力元件的肌 动力学仿真分析的人体下肢假肢优化系统的设计研 肉主要有股直肌、股外侧肌、股内侧肌、腓肠肌、胫骨 究中，将人体下肢简化为6刚体的二维力学模型， 前肌等。其中股直肌、股内侧肌、股外侧肌主要负责 大腿伸和屈、小腿伸和屈、膝关节的伸和髋关节的 并使用经典的Lagrange动力学方程对其建模[1)。 屈，并保持人体直立姿势：腓肠肌和胫骨前肌主要负 多刚体建模方法的核心是不考虑肌肉组成结 责小腿屈伸和人体的直立、行走[24]。 构，把大腿、小腿、上肢等结构视作一个个刚体，从系 统能量或动力学角度，对人体下肢进行建模，其中主 2人体下肢力学建模方法分类 要有Newton-Euler方法、Lagrange方法、角动量方 从关节和肌肉角度2个出发，可以得到不同的 法、Kane方法等[i9-2],其中Lagrange建模法应用最 下肢模型。从关节角度出发可以得到膝关节、髋关 为广泛，下面以这种方法为例介绍典型多刚体建模 节和踝关节力矩：从肌肉角度出发可以得到肌肉力， 方法[24」 而且肌肉力进一步可合成关节力矩。考虑到外骨骼 Lagrange建模法的核心是建立人体下肢多刚体 发展趋势，操作者的舒适性和安全性变得越来越重 模型的Lagrange动力学方程。其中，多轴加速度传 要。从个性差异化骨骼肌模型出发的肌肉力模型因 感器用来检测大腿和小腿的倾角并获取下肢运动过 与人体下肢结构更为相似，符合人体下肢建模的仿 程中的姿态以及膝关节角度和角速度；足底压力传 生设计的出发点，具有先天优势，已成为近几年外骨 感器被用来检测下肢的脚与地面之间接触状态和作 骼的发展热点[s6]」 用力等信息。结合人体结构和运动学的分析方法， 3基于关节力矩的下肢建模方法 建立人体下肢运动数学模型。从系统能量角度出发 构建动力学模型，并对得到的下肢模型进行动力学 由于人体下肢有着复杂的生物结构，并且限于 分析，进而得到关节力矩
非常困难。 为此，根据研究的目的，对人体或人体的 一部分作适当的简化，在一些必要的假设条件下，从 不同角度出发，建立人体下肢的生物动力学模型有 着巨大的意义［１］ 。 本文旨在对国内外人体下肢建模和仿真的最新 研究状况进行综述，并对今后的研究趋势进行阐述 分析，为进一步深入研究下肢生物动力学建模提供 一定的参考依据。 １ 人体下肢结构概述 从解剖学的角度看，人的躯体是由骨、骨连接和 骨骼肌组成，全身骨骨间借关节连成骨骼，构成主体 人体支架。 在运动过程中，骨、骨连接主要起杠杆作 用，而具体的肢体关节是运动的枢纽，骨骼肌是运动 的动力来源和主要执行部件。 在人体下肢运动中起主要参与的骨骼如下：髋 骨、股骨、腓骨、胫骨等。 骨与骨之间主要由结缔组织 连接，共同形成人体骨架，骨头之间连接的地方称为 关节。 而在人体下肢运动中起重要作用的关节主要 有 ３ 个：髋关节、膝关节和踝关节。 髋关节为多轴性 关节，可做屈、伸、内收、外展、内旋、外旋和旋转运动， 由于髋关节起到连接上下肢和承重的作用，其骨密度 较大，故而稳固性较好。 膝关节是人体负重比例最 大、结构最复杂的关节，其主要负责大腿和小腿的屈、 伸运动，和某些特定位置下的小幅度旋转运动。 踝关 节主要负责脚掌的背屈、伸、内收和外展运动。 与此同时，对人体下肢运动充当动力元件的肌 肉主要有股直肌、股外侧肌、股内侧肌、腓肠肌、胫骨 前肌等。 其中股直肌、股内侧肌、股外侧肌主要负责 大腿伸和屈、小腿伸和屈、膝关节的伸和髋关节的 屈，并保持人体直立姿势；腓肠肌和胫骨前肌主要负 责小腿屈伸和人体的直立、行走［２⁃４］ 。 ２ 人体下肢力学建模方法分类 从关节和肌肉角度 ２ 个出发，可以得到不同的 下肢模型。 从关节角度出发可以得到膝关节、髋关 节和踝关节力矩；从肌肉角度出发可以得到肌肉力， 而且肌肉力进一步可合成关节力矩。 考虑到外骨骼 发展趋势，操作者的舒适性和安全性变得越来越重 要。 从个性差异化骨骼肌模型出发的肌肉力模型因 与人体下肢结构更为相似，符合人体下肢建模的仿 生设计的出发点，具有先天优势，已成为近几年外骨 骼的发展热点［５⁃６］ 。 ３ 基于关节力矩的下肢建模方法 由于人体下肢有着复杂的生物结构，并且限于 医学伦理性问题，肌肉的生物动力学参数测量具有 一定困难，现有的大部分下肢建模方式均未考虑人 体下肢的生物模型，只考虑下肢的典型结构，或是建 立简化后的下肢多刚体模型［７⁃９］ ，或是用 ３Ｄ 建模技 术模拟简易的人体下肢［１０⁃１３］ ，或是采用黑箱训练等 方法跨过人体结构［１４］ ，进而通过动力学分析方法描 述模型，并对模型进行求解。 下面分别介绍几种典 型的下肢建模方法。 ３．１ Ｌａｇｒａｎｇｅ 多刚体模型 国内外的许多学者针对人体不同的运动或动作 情形，考虑不同情况下模型的具体要求，建立出了相 应的多刚体模型，其中最简单的方法是用单质点模 拟人体运动，不过模型过于简单，实际价值较低。 Ｈａｎａｖａｎ 于 １９６４ 年提出一个由 １５ 个多刚体环节组 成的，通过球铰系统连接的人体力学模型，这种模型 将人体主要的生物结构简化为几何形状的均质刚 体，建立并计算出基于人体惯性参数的回归方程，是 一套比较完整的具有个性化参数特征的多刚体人体 模型，这个模型也在一段时间内被连续引用［１５］ 。 波 兰学者 Ｂｌａｊｅｒａ 等为了研究杂技演员在蹦床上腾空 翻跟斗的动作过程，建立了一个自由度为 １０ 的多刚 体人体模型［１６］ 。 Ａｎｄｅｒｓｏｎ 等在研究人垂直起跳直 到达到最大高度过程中肌肉运动协调模式时，把人 体简化为 １０ 段、具有 ５４ 块肌肉的模型，并在此模型 的基础上进一步建立人在行走过程中行走距离和步 态与能量消耗间的关系［１７］ 。 Ｐｅｊｈａｎ 等在基于步态 动力学仿真分析的人体下肢假肢优化系统的设计研 究中，将人体下肢简化为 ６ 刚体的二维力学模型， 并使用经典的 Ｌａｇｒａｎｇｅ 动力学方程对其建模［１８］ 。 多刚体建模方法的核心是不考虑肌肉组成结 构，把大腿、小腿、上肢等结构视作一个个刚体，从系 统能量或动力学角度，对人体下肢进行建模，其中主 要有 Ｎｅｗｔｏｎ⁃Ｅｕｌｅｒ 方法、 Ｌａｇｒａｎｇｅ 方法、角动量方 法、Ｋａｎｅ 方法等［１９⁃２３］ ，其中 Ｌａｇｒａｎｇｅ 建模法应用最 为广泛，下面以这种方法为例介绍典型多刚体建模 方法［２４］ 。 Ｌａｇｒａｎｇｅ 建模法的核心是建立人体下肢多刚体 模型的 Ｌａｇｒａｎｇｅ 动力学方程。 其中，多轴加速度传 感器用来检测大腿和小腿的倾角并获取下肢运动过 程中的姿态以及膝关节角度和角速度；足底压力传 感器被用来检测下肢的脚与地面之间接触状态和作 用力等信息。 结合人体结构和运动学的分析方法， 建立人体下肢运动数学模型。 从系统能量角度出发 构建动力学模型，并对得到的下肢模型进行动力学 分析，进而得到关节力矩。 第 ４ 期 邵明旭，等：人体下肢生物力学建模研究进展 ·５１９·
.520. 智能系统学报 第10卷 考虑到人体下肢在行走和跳跃等主要活动中踝 骤发展都较为成熟。在人体下肢生物动力学模型研 关节转动轻微，因此忽略踝关节的转动，将人体下肢 究的早期，研究者们多用此方法来描述人体的运动 结构简化为平面的二刚体模型如图1所示。其中， 规律。 H表示髋关节，K为膝关节，A为踝关节。T和S分 1)选取关节广义变量及广义力 别代表大腿和小腿的质心。 首先选取笛卡尔坐标系，将大腿和小腿的关节 变量分别定义为运动过程中大腿和小腿的转角 H 和02，M。和M4分别代表髋关节和膝关节的力矩。 m1和m2分别表示大腿和小腿的质量，大、小腿长分 别为l和L2,质心分别用T和S代替，与关节中心的 距离分别为P,和p2。如上所述，大腿质心T的位置 坐标可用一组坐标描述： mg X=pisin 0 Y1=-P1c0s6 (3) 小腿质心S的位置坐标描述如下： X2=l1sin01+P2sin(01+02) Y2=-l1cos01-p2cos(91+02) (4) 2)系统动能K Ek=Ek +Eg a=ni诚Ba= 1 1 2m,14e+2m(0+0,)2+ m2l2p2(0+0,02)cos02 (5) 图1下肢二刚体二自由度模型 3)系统势能P Fig.1 The two rigid body and two de- E。=Eni+Ep2 grees of freedom model of lower E=migp (1 cos 0)E2 limb m282[1-c0s(01+02)]+ 在分析力学中，拉格朗日方程是解决受约束的 m2gl (1 -cos 0,) 质点系动力学问题的一个有力工具。拉格朗日函数 4)拉格朗日函数定义为 L被定义为系统动能K和势能P之差，即 L=Ek-E。 L=K-P (1) 5)系统动力学方程 其中系统动能K和势能P可以用任何方便的坐标 根据拉格朗日方程定义式计算各关节上的力 系来表示。用拉格朗日方程表示系统动力学方程式 矩，得到系统动力学方程。髋关节H上的力矩M 如下： 计算如下： d aLaL F:= (2) dt aq:aqi =(m,p+m,l9,+m,lp,(28,+,)cos8,+ a0, 式中：9：为表示动能和位能的坐标，9：为相应质点 m2P22(01+02) 的速度，而F:为作用在第i个坐标上的力或力矩。T aL 是作用在第i个坐标上的力或力矩，具体计算公式 a01 =-(m1P1+m2l1)gsin81-m2g2sin(81+02) 由9：为直线坐标或角坐标决定。这些力、力矩和坐 (7) 标分别称为广义力、广义力矩和广义坐标，n为连杆 得到 总数目。 d aL aL M6 (8) 用Lagrange方程解决实际力学问题是经典力 d业a8,a9, 学分析中一种非常成功的方法，具体内容和解题步 进一步简化为
考虑到人体下肢在行走和跳跃等主要活动中踝 关节转动轻微，因此忽略踝关节的转动，将人体下肢 结构简化为平面的二刚体模型如图 １ 所示。 其中， Ｈ 表示髋关节， Ｋ 为膝关节， Ａ 为踝关节。 Ｔ 和 Ｓ 分 别代表大腿和小腿的质心。 图 １ 下肢二刚体二自由度模型 Ｆｉｇ．１ Ｔｈｅ ｔｗｏ ｒｉｇｉｄ ｂｏｄｙ ａｎｄ ｔｗｏ ｄｅ⁃ ｇｒｅｅｓ ｏｆ ｆｒｅｅｄｏｍ ｍｏｄｅｌ ｏｆ ｌｏｗｅｒ ｌｉｍｂ 在分析力学中，拉格朗日方程是解决受约束的 质点系动力学问题的一个有力工具。 拉格朗日函数 Ｌ 被定义为系统动能 Ｋ 和势能 Ｐ 之差，即 Ｌ ＝ Ｋ － Ｐ （１） 其中系统动能 Ｋ 和势能 Ｐ 可以用任何方便的坐标 系来表示。 用拉格朗日方程表示系统动力学方程式 如下： Ｆｉ ＝ ｄ ｄｔ ∂Ｌ ∂ｑ · ｉ － ∂Ｌ ∂ｑｉ （２） 式中： ｑｉ 为表示动能和位能的坐标， ｑｉ ˙ 为相应质点 的速度，而 Ｆｉ 为作用在第 ｉ 个坐标上的力或力矩。 Ｔｉ 是作用在第 ｉ 个坐标上的力或力矩，具体计算公式 由 ｑｉ 为直线坐标或角坐标决定。 这些力、力矩和坐 标分别称为广义力、广义力矩和广义坐标， ｎ 为连杆 总数目。 用 Ｌａｇｒａｎｇｅ 方程解决实际力学问题是经典力 学分析中一种非常成功的方法，具体内容和解题步 骤发展都较为成熟。 在人体下肢生物动力学模型研 究的早期，研究者们多用此方法来描述人体的运动 规律。 １）选取关节广义变量及广义力 首先选取笛卡尔坐标系，将大腿和小腿的关节 变量分别定义为运动过程中大腿和小腿的转角 θ１ 和 θ２ ， Ｍｈ 和 Ｍｋ 分别代表髋关节和膝关节的力矩。 ｍ１ 和 ｍ２ 分别表示大腿和小腿的质量，大、小腿长分 别为 ｌ １ 和 ｌ ２ ，质心分别用 Ｔ 和 Ｓ 代替，与关节中心的 距离分别为 ｐ１ 和 ｐ２ 。 如上所述，大腿质心 Ｔ 的位置 坐标可用一组坐标描述： Ｘ１ ＝ ｐ１ ｓｉｎ θ１ Ｙ１ ＝ － ｐ１ ｃｏｓ θ１ （３） 小腿质心 Ｓ 的位置坐标描述如下： Ｘ２ ＝ ｌ １ ｓｉｎ θ１ ＋ ｐ２ ｓｉｎ θ１ ＋ θ２ ( ) Ｙ２ ＝ － ｌ １ ｃｏｓ θ１ － ｐ２ ｃｏｓ θ１ ＋ θ２ ( ) （４） ２）系统动能 Ｋ Ｅｋ ＝ Ｅｋ１ ＋ Ｅｋ２ Ｅｋ１ ＝ １ ２ ｍ１ ｐ ２ １ θ ·２ １Ｅｋ２ ＝ １ ２ ｍ２ ｌ ２ １ θ ·２ １ ＋ １ ２ ｍ２ ｐ ２ ２ θ · １ ＋ θ · ２ ( ) ２ ＋ ｍ２ ｌ ２ ｐ２ θ ·２ １ ＋ θ · １ θ · ２ ( ) ｃｏｓθ２ （５） ３）系统势能 Ｐ Ｅｐ ＝ Ｅｐ１ ＋ Ｅｐ２ Ｅｐ１ ＝ ｍ１ ｇｐ１ １ － ｃｏｓ θ１ ( ) Ｅｐ２ ＝ ｍ２ ｇｐ２ １ － ｃｏｓ θ１ ＋ θ２ [ ( ) ] ＋ ｍ２ ｇｌ １ １ － ｃｏｓ θ１ ( ) （６） ４）拉格朗日函数定义为 Ｌ ＝ Ｅｋ － Ｅｐ ５）系统动力学方程 根据拉格朗日方程定义式计算各关节上的力 矩，得到系统动力学方程。 髋关节 Ｈ 上的力矩 Ｍｈ 计算如下 ： ∂Ｌ ∂θ · １ ＝ ｍ１ ｐ ２ １ ＋ ｍ２ ｌ ２ １ ( ) θ · １ ＋ ｍ２ ｌ １ ｐ２ ２θ · １ ＋ θ · ２ ( ) ｃｏｓ θ２ ＋ ｍ２ ｐ２ ２ θ · １ ＋ θ · ２ ( ) ∂Ｌ ∂θ１ ＝ － ｍ１ ｐ１ ＋ ｍ２ ｌ １ ( ) ｇｓｉｎθ１ － ｍ２ ｇｐ２ ｓｉｎ θ１ ＋ θ２ ( ) （７） 得到 Ｍｈ ＝ ｄ ｄｔ ∂Ｌ ∂θ · １ － ∂Ｌ ∂θ１ （８） 进一步简化为 ·５２０· 智 能 系 统 学 报 第 １０ 卷
第4期 邵明旭，等：人体下肢生物力学建模研究进展 ·521· Mg=Du01+D1202+D20102+D12+D1(9) 系，但定义中不涉及力偶或力矩，并且运动量也不涉 膝关节K上的力矩M计算如下： 及角速度和角加速度。而角动量定理则描述了质点 =m,p(0,+a,)+m4p:8,es8, 系或刚体的运动量和力矩之间的关系。结合这2个 a02 定理就可完整地描述外力系与受力体运动参数之间 的定量关系。 aL a02 =-m2lP2(0+0,02)sin02- m2g2sin(81+02) (10) 进而得到 d aL aL M=a0,a0: (11) 进一步简化为 M=D101+Dn02+D20,82+ 速度运量运量定量 D26,2+D2 (12) (a)质点运动 最终得到人体下肢模型如下，带入大腿小腿质 量、关节角度等参数即可求得所需要的关节力矩： +C F(++dr) (13) d ) 式中： 面积速度→角动量→角动量定理 A-2 D.=Ds D B= (b)旋转状态 矩阵中详细参数如下： 图2质点运动及旋转状态示意图 Du mipi map2 m2l+2mzl pzcos 0 Fig.2 Sketch of the motion and rotating state of the particle D12=m2P3+m2l1P2c0s02 如图2所示，质点的运动状态以矢量形式定义： D21=m2P2+m1l1P2c0s62 (r、),相对某参考点的位置矢量为r,速度为v。 D22=m2P2 则动量定理为 D11=0 d(2mv)F·d (14) D122 =-m2liP2sin 02 Dau malipasin 02 惯性系$中的一个质点在相对某参考点O的运 D22=0 动过程中，其相对参考点0的径矢会相应的旋转。 Du=-2m2l Pzsin 02 假设在dt时间质点位移为vdt,转过的角度为d9, D121=0 r便会扫过面积dS: D22 =-m2l pasin 02 m2l pzsin 02 dsrxl (15) D =(mip m)gsin 0 m2p2gsin(0 +02) 面积速度为 D2=m2gP2sin(81+02) 3.2角动量定理建模法 di=rxv (16) 经典动力学研究中的3个定理分别为动量定 质点在惯性系S中相对参考点O的角动量L: 理、角动量定理和动能定理。Lagrange方程从能量 L=rXmw=r×p (17) 角度出发，应用动能定理和能量守恒建立和求解下 肢动力方程，而从角动量定理和动量定理角度出发 角动量L随时间的变化如下： 也可以建立人体多质点模型和多刚体模型的动力学 dL_d(r×p)_dr (18) dt dt dt xp+rx史 dt 方程252，求出关节力矩。 式中： 动量定理描述了物体的运动状态和力之间的关
Ｍｈ ＝ Ｄ１１ θ ¨ １ ＋ Ｄ１２ θ ¨ ２ ＋ Ｄ１１２ θ · １ θ · ２ ＋ Ｄ１２２ θ ·２ ２ ＋ Ｄ１ （９） 膝关节 Ｋ 上的力矩 Ｍｋ 计算如下： ∂Ｌ ∂θ · ２ ＝ ｍ２ ｐ ２ ２ θ · １ ＋ θ · ２ ( ) ＋ ｍ２ ｌ １ ｐ２ θ · １ ｃｏｓθ２ ∂Ｌ ∂θ２ ＝ － ｍ２ ｌ １ ｐ２ θ ·２ １ ＋ θ · １ θ · ２ ( ) ｓｉｎθ２ － ｍ２ ｇｐ２ ｓｉｎ θ１ ＋ θ２ ( ) （１０） 进而得到 Ｍｋ ＝ ｄ ｄｔ ∂Ｌ ∂θ · ２ － ∂Ｌ ∂θ２ （１１） 进一步简化为 Ｍｋ ＝ Ｄ２１ θ ¨ １ ＋ Ｄ２２ θ ¨ ２ ＋ Ｄ２１２ θ · １ θ · ２ ＋ Ｄ２１１ θ · １ ２ ＋ Ｄ２ （１２） 最终得到人体下肢模型如下，带入大腿小腿质 量、关节角度等参数即可求得所需要的关节力矩： Ｍｈ Ｍｋ é ë ê ê ù û ú ú ＝ Ａ θ ¨ １ θ ¨ ２ é ë ê ê ù û ú ú ＋ Ｂ θ ·２ １ θ ·２ ２ é ë ê ê ê ù û ú ú ú ＋ Ｃ θ · １ θ · ２ θ · ２ θ · １ é ë ê ê ê ù û ú ú ú ＋ Ｄ１ Ｄ２ é ë ê ê ù û ú ú （１３） 式中： Ａ ＝ Ｄ１１ Ｄ１２ Ｄ２１ Ｄ２２ é ë ê ê ù û ú ú ，Ｂ ＝ Ｄ１１１ Ｄ１２２ Ｄ２１１ Ｄ２２２ é ë ê ê ù û ú ú ，Ｃ ＝ Ｄ１２２ Ｄ１２１ Ｄ２１２ Ｄ２２１ é ë ê ê ù û ú ú 。 矩阵中详细参数如下： Ｄ１１ ＝ ｍ１ ｐ ２ １ ＋ ｍ２ ｐ ２ ２ ＋ ｍ２ ｌ ２ １ ＋ ２ｍ２ ｌ １ ｐ２ ｃｏｓ θ２ Ｄ１２ ＝ ｍ２ ｐ ２ ２ ＋ ｍ２ ｌ １ ｐ２ ｃｏｓ θ２ Ｄ２１ ＝ ｍ２ ｐ ２ ２ ＋ ｍ１ ｌ １ ｐ２ ｃｏｓ θ２ Ｄ２２ ＝ ｍ２ ｐ ２ ２ Ｄ１１１ ＝ ０ Ｄ１２２ ＝ － ｍ２ ｌ １ ｐ２ ｓｉｎ θ２ Ｄ２１１ ＝ ｍ２ ｌ １ ｐ２ ｓｉｎ θ２ Ｄ２２２ ＝ ０ Ｄ１１２ ＝ － ２ｍ２ ｌ １ ｐ２ ｓｉｎ θ２ Ｄ１２１ ＝ ０ Ｄ２１２ ＝ － ｍ２ ｌ １ ｐ２ ｓｉｎ θ２ ＋ ｍ２ ｌ １ ｐ２ ｓｉｎ θ２ Ｄ１ ＝ ｍ１ ｐ１ ＋ ｍ１ ｌ １ ( ) ｇｓｉｎ θ１ ＋ ｍ２ ｐ２ ｇｓｉｎ θ１ ＋ θ２ ( ) Ｄ２ ＝ ｍ２ ｇｐ２ ｓｉｎ θ１ ＋ θ２ ( ) ３．２ 角动量定理建模法 经典动力学研究中的 ３ 个定理分别为动量定 理、角动量定理和动能定理。 Ｌａｇｒａｎｇｅ 方程从能量 角度出发，应用动能定理和能量守恒建立和求解下 肢动力方程，而从角动量定理和动量定理角度出发 也可以建立人体多质点模型和多刚体模型的动力学 方程［２５⁃２８］ ，求出关节力矩。 动量定理描述了物体的运动状态和力之间的关 系，但定义中不涉及力偶或力矩，并且运动量也不涉 及角速度和角加速度。 而角动量定理则描述了质点 系或刚体的运动量和力矩之间的关系。 结合这 ２ 个 定理就可完整地描述外力系与受力体运动参数之间 的定量关系。 （ａ） 质点运动 （ｂ） 旋转状态 图 ２ 质点运动及旋转状态示意图 Ｆｉｇ．２ Ｓｋｅｔｃｈ ｏｆ ｔｈｅ ｍｏｔｉｏｎ ａｎｄ ｒｏｔａｔｉｎｇ ｓｔａｔｅ ｏｆ ｔｈｅ ｐａｒｔｉｃｌｅ 如图 ２ 所示，质点的运动状态以矢量形式定义： (ｒ、ｖ) ，相对某参考点的位置矢量为 ｒ ，速度为 ｖ 。 则动量定理为 ｄ（ １ ２ ｍｖ·ｖ） «Ｆ·ｄｒ （１４） 惯性系 Ｓ 中的一个质点在相对某参考点 Ｏ 的运 动过程中，其相对参考点 Ｏ 的径矢 ｒ 会相应的旋转。 假设在 ｄｔ 时间质点位移为 ｖｄｔ ， 转过的角度为 ｄθ ， ｒ 便会扫过面积 ｄＳ ： ｄＳ ＝ １ ２ ｒ × ｖ ｄｔ （１５） 面积速度为 Ｋ ＝ ｄＳ ｄｔ ＝ １ ２ ｒ × ｖ （１６） 质点在惯性系 Ｓ 中相对参考点 Ｏ 的角动量 Ｌ ： Ｌ ＝ ｒ × ｍｖ ＝ ｒ × ｐ （１７） 角动量 Ｌ 随时间的变化如下： ｄＬ ｄｔ ＝ ｄ（ｒ × ｐ） ｄｔ ＝ ｄｒ ｄｔ × ｐ ＋ ｒ × ｄｐ ｄｔ （１８） 式中： 第 ４ 期 邵明旭，等：人体下肢生物力学建模研究进展 ·５２１·
·522. 智能系统学报 第10卷 dr 骨骼康复机器人的动力学模型，以角度、角速度和角 Xp=v×p=0 dt 加速度作为输入信号，输出信号为髋、膝关节力矩， 二F (19) 进行逆动力学仿真分析，获得了外骨骼逆动力学动 dt 态数学模型[35]。 于是 3.4SVR黑箱训练法 d业=rF (20) 下肢多刚体模型实质上采用传统的动力学方程 d 求解，具有模型精准度不高，计算关节力矩受参数影 质点所受到相对参考点O的力矩为 响大，不利于实时控制等缺点。针对这些缺点，Wi M= =rxF (21) Meng等提出了SVR黑箱训练法对下肢进行建模， dt 彻底放弃考虑人体下肢结构，并且证明该方法有利 根据下肢的平面二刚体模型，由角动量定理可 于实时控制。该方法从4个下肢肌肉提取的EMG 以得出： 信号均方根值(RMS)作为输入向量，以力传感器采 髋关节力矩为 集到的下肢关节力矩作为训练数据，将其输入支持 M=J 0 mip 01 migpisin 0 (22) 向量回归(SVR)模型来估算下肢关节力矩。使用自 膝关节力矩为 回归算法来构建EMG信号和下肢关节力矩之间的 Me=J202 m2p02 -m2gp2sin 62 (23) 关系8】。 式中：J和J2分别代表实验对象的大、小腿固有转 在本方法中，$VR算法被用来预测下肢力量与 动惯量。 肌电信号之间的关系。SVR黑箱算法的一个关键 3.3仿真软件建模法 问题是用函数来最大限度的接近模型，即找到一个 仿真软件建模属于多刚体建模法中的一种，将 确定性的函数g(x)代替模型本身非确定性函数 人体下肢主要结构在ADAMS或SimMechanics等仿 f(x),可用式(24)进行说明： 真分析或工业设计软件中用一些环节和约束函数代 替，再由视频采集软件或运动姿态采集设备提供肢 R(f.g)=L(f.g)dx (24) 体角度、大小腿长度等输入参数，通过软件强大的仿 其中，L(f,g)代表惩罚函数。函数f(x)表示样本 真功能，直接输出想要得到的关节力矩93。 序列，g(x)可以通过回归分析来获得。最后，SVR 在国内，朱昌义等在研究人体在做单杠摆运动 模型的求解可被转换为一组拉格朗日对偶问题，最 过程特点时，建立了一个5环节组成的多刚体力学 优化目标函数和约束条件如下： 模型，并且结合经典力学分析中的Kane方程，推导 min(O++03) 出人体行走步态规划过程中下肢摆动的动力学方 程)。刘明辉等在对外骨骼机器的研究中把人体 2a,-a)g-g)(x) Q1= 21 简化成一个九连杆的多刚体模型，并在机械系统动 力学自动分析软件(automatic dynamic analysis of Q2= ae-) i=1 mechanical system,ADAMS)中建立了简化后的人体 多刚体模型，对人体行走模型及步态分析进行了仿 0-Za(e+y (25) i=1 真分析，该研究为在虚拟环境中模拟现实情况下人 约束如下： 机外骨骼刚柔耦合系统提供了实验依据别。沈凌 等在对人体下肢外骨骼假肢的研究过程中，建立了 (a:-a)=0 i=1 (26) 基于四连杆系统的膝关节模型，并利用Lagrange方 a:,a∈[0，C] 法推导出下肢动力学方程，最后在Pro/E(Prol Engi- 然后，可以得到模型的回归方程： neer)软件中建立了下肢假肢结构的仿真模型，通过 系统集成的动力学仿真功能，得出的膝关节角度变 )=∑a-a)(0+b(27) 化和输出力矩之间的关系曲线，为下肢动力学建模 式中：SV代表支持向量空间，a:、a;是拉格朗日乘 提供了一定的参考数据。 子系数向量，b为偏移量。利用核函数来描述线性 洪晓明等人通过Solidworks软件建立精确的三 SVR向量到非线性SVR向量的映射。它们之间的 维实体模型，联合Matlab/SimMechanics建立下肢外 映射关系如下：
ｄｒ ｄｔ × ｐ ＝ ｖ × ｐ ＝ ０ ｄｐ ｄｔ ＝ Ｆ （１９） 于是 ｄＬ ｄｔ ＝ ｒ × Ｆ （２０） 质点所受到相对参考点 Ｏ 的力矩为 Ｍ ＝ ｄＬ ｄｔ ＝ ｒ × Ｆ （２１） 根据下肢的平面二刚体模型，由角动量定理可 以得出： 髋关节力矩为 Ｍｈ ＝ Ｊｘ１ θ ¨ １ ＋ ｍ ２ １ ｐ１ θ ¨ １ － ｍ１ ｇｐ１ ｓｉｎ θ１ （２２） 膝关节力矩为 Ｍｋ ＝ Ｊｘ２ θ ¨ ２ ＋ ｍ２ ｐ ２ ２ θ ¨ ２ － ｍ２ ｇｐ２ ｓｉｎ θ２ （２３） 式中： Ｊｘ１ 和 Ｊｘ２ 分别代表实验对象的大、小腿固有转 动惯量。 ３．３ 仿真软件建模法 仿真软件建模属于多刚体建模法中的一种，将 人体下肢主要结构在 ＡＤＡＭＳ 或 ＳｉｍＭｅｃｈａｎｉｃｓ 等仿 真分析或工业设计软件中用一些环节和约束函数代 替，再由视频采集软件或运动姿态采集设备提供肢 体角度、大小腿长度等输入参数，通过软件强大的仿 真功能，直接输出想要得到的关节力矩［２９⁃３１］ 。 在国内，朱昌义等在研究人体在做单杠摆运动 过程特点时，建立了一个 ５ 环节组成的多刚体力学 模型，并且结合经典力学分析中的 Ｋａｎｅ 方程，推导 出人体行走步态规划过程中下肢摆动的动力学方 程［３２］ 。 刘明辉等在对外骨骼机器的研究中把人体 简化成一个九连杆的多刚体模型，并在机械系统动 力学自动分析软件 （ ａｕｔｏｍａｔｉｃ ｄｙｎａｍｉｃ ａｎａｌｙｓｉｓ ｏｆ ｍｅｃｈａｎｉｃａｌ ｓｙｓｔｅｍ， ＡＤＡＭＳ）中建立了简化后的人体 多刚体模型，对人体行走模型及步态分析进行了仿 真分析，该研究为在虚拟环境中模拟现实情况下人 机外骨骼刚柔耦合系统提供了实验依据［３３］ 。 沈凌 等在对人体下肢外骨骼假肢的研究过程中，建立了 基于四连杆系统的膝关节模型，并利用 Ｌａｇｒａｎｇｅ 方 法推导出下肢动力学方程，最后在 Ｐｒｏ ／ Ｅ（Ｐｒｏｌ Ｅｎｇｉ⁃ ｎｅｅｒ）软件中建立了下肢假肢结构的仿真模型，通过 系统集成的动力学仿真功能，得出的膝关节角度变 化和输出力矩之间的关系曲线，为下肢动力学建模 提供了一定的参考数据［３４］ 。 洪晓明等人通过 Ｓｏｌｉｄｗｏｒｋｓ 软件建立精确的三 维实体模型，联合 Ｍａｔｌａｂ ／ ＳｉｍＭｅｃｈａｎｉｃｓ 建立下肢外 骨骼康复机器人的动力学模型，以角度、角速度和角 加速度作为输入信号，输出信号为髋、膝关节力矩， 进行逆动力学仿真分析，获得了外骨骼逆动力学动 态数学模型［３５⁃３７］ 。 ３．４ ＳＶＲ 黑箱训练法 下肢多刚体模型实质上采用传统的动力学方程 求解，具有模型精准度不高，计算关节力矩受参数影 响大，不利于实时控制等缺点。 针对这些缺点，Ｗｅｉ Ｍｅｎｇ 等提出了 ＳＶＲ 黑箱训练法对下肢进行建模， 彻底放弃考虑人体下肢结构，并且证明该方法有利 于实时控制。 该方法从 ４ 个下肢肌肉提取的 ＥＭＧ 信号均方根值（ＲＭＳ）作为输入向量，以力传感器采 集到的下肢关节力矩作为训练数据，将其输入支持 向量回归（ＳＶＲ）模型来估算下肢关节力矩。 使用自 回归算法来构建 ＥＭＧ 信号和下肢关节力矩之间的 关系［３８］ 。 在本方法中，ＳＶＲ 算法被用来预测下肢力量与 肌电信号之间的关系。 ＳＶＲ 黑箱算法的一个关键 问题是用函数来最大限度的接近模型，即找到一个 确定性的函数 ｇ(ｘ) 代替模型本身非确定性函数 ｆ(ｘ) ，可用式（２４）进行说明： Ｒ(ｆ，ｇ) ＝ ∫Ｌ(ｆ，ｇ) ｄｘ （２４） 其中， Ｌ(ｆ，ｇ) 代表惩罚函数。 函数 ｆ(ｘ) 表示样本 序列， ｇ(ｘ) 可以通过回归分析来获得。 最后，ＳＶＲ 模型的求解可被转换为一组拉格朗日对偶问题，最 优化目标函数和约束条件如下： ｍｉｎ Ｑ１ ＋ Ｑ２ ＋ Ｑ３ { } Ｑ１ ＝ １ ２ ∑ ｍ ｉ ＝ １ ａｉ － ａ ∗ ｉ ( ) ａｊ － ａ ∗ ｊ ( ) 〈ｘｉ，ｘｊ〉 Ｑ２ ＝ ∑ ｍ ｉ ＝ １ ａｉ ε － ｙｉ ( ) Ｑ３ ＝ ∑ ｍ ｉ ＝ １ ａ ∗ ｉ ε ＋ ｙｉ ( ) （２５） 约束如下： ∑ ｍ ｉ ＝ １ ａｉ － ａ ∗ ｉ ( ) ＝ ０ ａｉ，ａ ∗ ｉ ∈ [０，Ｃ] ì î í ï ï ïï （２６） 然后，可以得到模型的回归方程： ｆ(ｘ) ＝ ∑ＳＶ ａｉ － ａ ∗ ｉ ( ) 〈ｘｉ·ｘ〉 ＋ ｂ （２７） 式中：ＳＶ 代表支持向量空间， ａｉ 、 ａ ∗ ｉ 是拉格朗日乘 子系数向量， ｂ 为偏移量。 利用核函数来描述线性 ＳＶＲ 向量到非线性 ＳＶＲ 向量的映射。 它们之间的 映射关系如下： ·５２２· 智 能 系 统 学 报 第 １０ 卷
第4期 邵明旭，等：人体下肢生物力学建模研究进展 .523. f)=∑(a,-a)Ke,)+b(28) 所示。 i=1 式中：z:和z。是训练样本，K和径向基函数(RBF) 作为核函数，用来逼近真实函数。 CE 最常用的径向基函数是高斯核函数[o]: Kc(x-x)=exp{-llx-x2/(2*)2) PE (29) 式中：x。为核函数中心。σ作为函数的宽度参数控 SE W 制了函数的径向作用范围。核函数中心可以通过 1 cos a K-Means聚类算法来得到聚类的中心向量，而σ其 实就是聚类中心的半径了。这个半径等于聚类中心 图3Hl三元素法模型 向量和属于该类的样本之间的距离的平均值。 Fig.3 Hill Three-element model 考虑EMG信号与肢体的力之间的关系，从不同 H山三元素模型考虑肌肉的本质结构，将肌肉 情况出发，根据经验公式或者通过对实验数据回归 简化为3种元素，即收缩元PE、串联弹性元SE和并 分析得到拉格朗日乘子、网络节点数、宽度参数等训 联弹性元CE来反映肌肉的功能。用一个弹簧阻尼 练参数，即可由SVR模型直接推导出关节力矩与 系统代表肌肉模型，如图4所示。 EMG信号的关系。 4基于肌肉力的下肢建模方法 从肌肉的本质模型角度出发，骨骼肌力预测有 3种代表性的肌肉模型一Hil模型、Huxley模型和 流变学模型。其中，应用在下肢建模方面最广泛且 最成功的是Hil三元素模型。Nikos Karavas等人提 出了一种基于Hl三元素法的下肢模型，下面详细 阐述Hill建模法的基本原理[o42] 肌肉的基本机能是把化学能通过做功形式转化 (a)超微结构 为机械能，表征肌肉活动的基本生物力学指标主要 有2个：1)肌肉端部测量出来的肌肉对骨骼的拉 肌动蛋白单体原肌球蛋白质肌原蛋白 力，即肌张力：2)运动过程中肌肉纤维长度变化的 0 3 Thl TnT 细肌丝 TnC 速度，即收缩速度。这2个指标是与肌肉实际活动 8888888文 8888的 密切相关且可测的宏观量。H在1938年以青蛙 0 的缝匠肌为试样，做了快速电刺激肌肉释放实验，利 用能量原理及经验假设，给出了肌肉收缩力F与收 肌球蛋白杆乃 0 粗肌丝 缩速度V间的关系式： 肌球蛋白头（横桥） (F+a)(e+b)=(F。+a)b (30) (b)肌丝结构 式中：F。是肌肉最大等长收缩力，a、b是力学常 图4基于超微结构和肌丝结构的骨骼肌模型 数。这就是著名的经典Hl方程。F。曲线及常数 Fig.4 Skeletal muscle model based on ultra- a、b可以通过纯力学的方法来确定，而不需要额外 structure and muscle wire structure 再做实验测定。H山方程是由快速电刺激释放痉挛 该模型中，各种元素的意义和肌肉的超微结构 态的青蛙缝匠肌的实验推导而来的，它仅仅揭示了 一一对应。收缩元代表着滑动的肌动蛋白和肌球蛋 肌肉收到电刺激后快速收缩时产生的肌肉力与肌肉 白分子，串联弹性元代表肌节和结缔组织的固有弹 长度和收缩速度之间的关系，并不能描述肌肉缓慢 性，并联弹性元代表连接肌肉骨骼的结缔组织。 释放时的肌张力与速度关系，也不能描述不受刺激 肌肉力是由于神经电信号对肌肉收缩单元的刺 处于正常状态下肌肉的力学性质。基于这些缺陷， 激产生的，因此Hl肌力模型应考虑2个关键问题， H山之后提出了改进后的肌肉三元素模型。如图3 分别为“EG-激活程度”和“激活程度-肌力”之间
ｆ ｚ０ ( ) ＝ ∑ ｌ ｉ ＝ １ ａｉ － ａ ∗ ｉ ( ) Ｋ ｚｉ，ｚ０ ( ) ＋ ｂ （２８） 式中： ｚｉ 和 ｚｏ 是训练样本， Ｋ 和径向基函数（ＲＢＦ） 作为核函数，用来逼近真实函数。 最常用的径向基函数是高斯核函数［４０］ ： ＫＣ ｘ － ｘｃ ( ) ＝ ｅｘｐ － ‖ｘ － ｘｃ‖２ ／ (２∗σ) ２ { } （２９） 式中： ｘｃ 为核函数中心。 σ 作为函数的宽度参数控 制了函数的径向作用范围。 核函数中心可以通过 Ｋ⁃Ｍｅａｎｓ 聚类算法来得到聚类的中心向量，而 σ 其 实就是聚类中心的半径了。 这个半径等于聚类中心 向量和属于该类的样本之间的距离的平均值。 考虑 ＥＭＧ 信号与肢体的力之间的关系，从不同 情况出发，根据经验公式或者通过对实验数据回归 分析得到拉格朗日乘子、网络节点数、宽度参数等训 练参数，即可由 ＳＶＲ 模型直接推导出关节力矩与 ＥＭＧ 信号的关系。 ４ 基于肌肉力的下肢建模方法 从肌肉的本质模型角度出发，骨骼肌力预测有 ３ 种代表性的肌肉模型—Ｈｉｌｌ 模型、Ｈｕｘｌｅｙ 模型和 流变学模型。 其中，应用在下肢建模方面最广泛且 最成功的是 Ｈｉｌｌ 三元素模型。 Ｎｉｋｏｓ Ｋａｒａｖａｓ 等人提 出了一种基于 Ｈｉｌｌ 三元素法的下肢模型，下面详细 阐述 Ｈｉｌｌ 建模法的基本原理［４０⁃４２］ 。 肌肉的基本机能是把化学能通过做功形式转化 为机械能，表征肌肉活动的基本生物力学指标主要 有 ２ 个：１） 肌肉端部测量出来的肌肉对骨骼的拉 力，即肌张力；２） 运动过程中肌肉纤维长度变化的 速度，即收缩速度。 这 ２ 个指标是与肌肉实际活动 密切相关且可测的宏观量。 Ｈｉｌｌ 在 １９３８ 年以青蛙 的缝匠肌为试样，做了快速电刺激肌肉释放实验，利 用能量原理及经验假设，给出了肌肉收缩力 Ｆ 与收 缩速度 Ｖ 间的关系式 ： (Ｆ ＋ ａ) (ｖ ＋ ｂ) ＝ Ｆ０ ( ＋ ａ) ｂ （３０） 式中： Ｆ０ 是肌肉最大等长收缩力， ａ 、 ｂ 是力学常 数。 这就是著名的经典 Ｈｉｌｌ 方程。 Ｆ０ 曲线及常数 ａ 、 ｂ 可以通过纯力学的方法来确定，而不需要额外 再做实验测定。 Ｈｉｌｌ 方程是由快速电刺激释放痉挛 态的青蛙缝匠肌的实验推导而来的，它仅仅揭示了 肌肉收到电刺激后快速收缩时产生的肌肉力与肌肉 长度和收缩速度之间的关系，并不能描述肌肉缓慢 释放时的肌张力与速度关系，也不能描述不受刺激 处于正常状态下肌肉的力学性质。 基于这些缺陷， Ｈｉｌｌ 之后提出了改进后的肌肉三元素模型。 如图 ３ 所示。 图 ３ Ｈｉｌｌ 三元素法模型 Ｆｉｇ．３ Ｈｉｌｌ Ｔｈｒｅｅ⁃ｅｌｅｍｅｎｔ ｍｏｄｅｌ Ｈｉｌｌ 三元素模型考虑肌肉的本质结构，将肌肉 简化为 ３ 种元素，即收缩元 ＰＥ、串联弹性元 ＳＥ 和并 联弹性元 ＣＥ 来反映肌肉的功能。 用一个弹簧阻尼 系统代表肌肉模型，如图 ４ 所示。 （ａ） 超微结构 （ｂ） 肌丝结构 图 ４ 基于超微结构和肌丝结构的骨骼肌模型 Ｆｉｇ． ４ Ｓｋｅｌｅｔａｌ ｍｕｓｃｌｅ ｍｏｄｅｌ ｂａｓｅｄ ｏｎ ｕｌｔｒａ⁃ ｓｔｒｕｃｔｕｒｅ ａｎｄ ｍｕｓｃｌｅ ｗｉｒｅ ｓｔｒｕｃｔｕｒｅ 该模型中，各种元素的意义和肌肉的超微结构 一一对应。 收缩元代表着滑动的肌动蛋白和肌球蛋 白分子，串联弹性元代表肌节和结缔组织的固有弹 性，并联弹性元代表连接肌肉骨骼的结缔组织。 肌肉力是由于神经电信号对肌肉收缩单元的刺 激产生的，因此 Ｈｉｌｌ 肌力模型应考虑 ２ 个关键问题， 分别为“ＥＭＧ⁃激活程度” 和“激活程度⁃肌力” 之间 第 ４ 期 邵明旭，等：人体下肢生物力学建模研究进展 ·５２３·
.524. 智能系统学报 第10卷 的关系描述。前者表示了神经电刺激与实验测得的 肌肉力的预测计算，具体流程如下。 表面肌电信号之间的潜在关系，后者则表示了肌肉 EMG 模糊 的激活程度同所产生的肌张力之间的关系。在H山 原始 滤波器 信号 推理机 Hill 三元素模型图中，为肌腱和肌纤维的总长度，L, 肌肉 模型 肌肉力 为肌腱长度，1为肌纤维长度，α为肌肉纤维束与 肌束和肌腱长度 肌腱的夹角。 图5Hl肌力预测模型流程图 根据上诉的三元素结构，肌肉力的生成用弹簧 Fig.5 Flow chart of Hill muscle strength prediction model 阻尼系统动力学方程来描述，其中肌腱单元所产生 模糊逻辑推理模型主要有Mamdani和Sugeno 的合力如下： 两种方法。首先，将采样得到的EMG原始信号进行 F()=(E+F)cos(9,()) (31) 滤波和平滑处理：其次对EMG原始信号进行滤波和 其中，F由式(32)求出： 平滑处理，将其用模糊集进行表示。然后建立模糊 =b ,°+a,i9 推理规则，根据处理后的EMG信号幅值，推理出作 (32) b-i理 为肌肉激活程度的隶属度函数。然后结合肌束和肌 式中：α：为神经单元作用在肌肉上的激活程度，b是 腱的长度，通过H肌力模型，从而计算出肌肉力。 在Matlab中的模糊工具箱中设计了模糊推理规则 肌肉固有的常数，和L分别表示肌肉的长度和 和隶属度函数，通过Matlab编程，实现对Hil肌力 收缩速度，F等张力则由式(33)给出： 预测模型的计算。 E,°=1- (33) 计算出肌肉力之后，就可以从力矩的定义出发， 将肌肉力与各自的力臂相乘，之后对关节进行简单 式中：F“是肌肉的最大等张力，能通过实验测量得 的受力分析，得出所需要的关节力矩4。 到。心则是表征肌肉宽度的形参，对人体骨骼肌应 肌肉力臂，：(t)可通过对骨骼肌长度和形状求 取0.56L0r。而肌纤维夹角9，()则是随着肌纤 导得到，公式如下： 维长度的变化而不断改变的，能通过式计算： L (t)=Co+C0mee(t)+C20(t) (t)=arcsin (34) al (t) L"(t) r:(t)=- (35) d0tnce 式中：L(t)是t时刻的肌纤维长度，9则是在优 代入后得到力臂：(t),其中，C。、C,和C,均 化肌纤维长度处的肌腱-肌纤维夹角。 为常数。 肌电信号是目前测量肌肉激活程度的唯一手 段。但由于SEMG信号在实验测量中容易受到表面 (0）=[C,+2C,90]180 (36) 湿度、皮下脂肪等因素的干扰，原始的EMG信号往 关节力矩M(t)计算如下： 往很难直接应用于肌力预测，必须经过一系列的数 M)= 含含0 (37) 字信号处理才能使用。依据H肌力模型，肌肉激 活程度表现为0~1的某个数值。0表示肌肉没有任 式中：M表示所有与关节相连的主动肌力矩，M表 何激活，1则表示肌肉完全激活。这是2种极端情 示所有拮抗肌力矩，如下： 况，实际中肌肉激活程度是介于0~1的常数，其值 M:(t)=F,(t)r:(t) 大小可通过归一化实验处理得到或者通过模糊推理 M,(t)=F()r(t) (38) 机得到。本文以模糊推理机为例介绍基于H模型 下肢建模仿真技术手段 的下肢建模方法[3】 模糊理论是随着模糊数学的出现而蓬勃发展 下肢建模的仿真验证是建模过程中的重要一 的，其核心是确定0~1的隶属度函数值，非常适合 环，也因为真实值测量困难、参数复杂多样等因素一 于估计SEMG的肌肉激活常数。本方法根据采样得 直是建模过程中的难点。目前国内外对下肢建模的 到的EMG信号，通过对其进行模糊化处理，经过一 仿真验证主要采取以下几种方法。 系列模糊规则推理，采用模糊隶属度函数来表示肌 孙棕檀等人将基于Lagrange方程或Kane方程 肉激活程度，从而作为H山肌肉模型的输入，再根据 等传统动力学分析的结果，与Simulink等仿真软件 个人差异化肌肉参数，输入到H山三元素模型，进行 输出的结果进行比较分析，证明了传统力学分析应
的关系描述。 前者表示了神经电刺激与实验测得的 表面肌电信号之间的潜在关系，后者则表示了肌肉 的激活程度同所产生的肌张力之间的关系。 在 Ｈｉｌｌ 三元素模型图中， ｌｍｔ 为肌腱和肌纤维的总长度， ｌ ｔ 为肌腱长度， ｌｍ 为肌纤维长度， α 为肌肉纤维束与 肌腱的夹角。 根据上诉的三元素结构，肌肉力的生成用弹簧 阻尼系统动力学方程来描述，其中肌腱单元所产生 的合力如下： Ｆ ｍｔ ｊ (ｔ) ＝ Ｆ ＣＥ ｊ ＋ Ｆ ＰＥ ｊ ( ) ｃｏｓ φｊ ( (ｔ) ) （３１） 其中， Ｆ ＣＥ ｊ 由式（３２）求出： Ｆ ＣＥ ｊ ＝ ｂＦｊ ０ ＋ ａｊＬ ·ＣＥ ｊ ｂ － Ｌ ·ＣＥ ｊ （３２） 式中： ａｊ 为神经单元作用在肌肉上的激活程度， ｂ 是 肌肉固有的常数， Ｌ ＣＥ ｊ 和 Ｌ ＣＥ ｊ ˙ 分别表示肌肉的长度和 收缩速度， Ｆ ０ ｊ 等张力则由式（３３）给出： Ｆｊ ０ ＝ Ｆ ｍａｘ ｊ １ － Ｌ ＣＥ ｊ － Ｌ ＣＥＯＰＴ ｊ ｗ æ è ç ö ø ÷ ２ é ë ê ê ù û ú ú （３３） 式中： Ｆ ｍａｘ ｊ 是肌肉的最大等张力，能通过实验测量得 到。 ｗ 则是表征肌肉宽度的形参，对人体骨骼肌应 取 ０．５６ Ｌ ＣＥＯＰＴ ｊ 。 而肌纤维夹角 φｊ (ｔ) 则是随着肌纤 维长度的变化而不断改变的，能通过式计算： φｊ (ｔ) ＝ ａｒｃｓｉｎ Ｌ ０ ｍｊ ｓｉｎ φ ０ ｊ Ｌ ｍ ｊ (ｔ) æ è ç ö ø ÷ （３４） 式中： Ｌ ｍ ｊ (ｔ) 是 ｔ 时刻的肌纤维长度， φ ０ ｊ 则是在优 化肌纤维长度 Ｌ ０ ｍｊ 处的肌腱－肌纤维夹角。 肌电信号是目前测量肌肉激活程度的唯一手 段。 但由于 ＳＥＭＧ 信号在实验测量中容易受到表面 湿度、皮下脂肪等因素的干扰，原始的 ＥＭＧ 信号往 往很难直接应用于肌力预测，必须经过一系列的数 字信号处理才能使用。 依据 Ｈｉｌｌ 肌力模型，肌肉激 活程度表现为 ０～１ 的某个数值。 ０ 表示肌肉没有任 何激活，１ 则表示肌肉完全激活。 这是 ２ 种极端情 况，实际中肌肉激活程度是介于 ０ ～ １ 的常数，其值 大小可通过归一化实验处理得到或者通过模糊推理 机得到。 本文以模糊推理机为例介绍基于 Ｈｉｌｌ 模型 的下肢建模方法［ ４３ ］ 。 模糊理论是随着模糊数学的出现而蓬勃发展 的，其核心是确定 ０ ～ １ 的隶属度函数值，非常适合 于估计 ＳＥＭＧ 的肌肉激活常数。 本方法根据采样得 到的 ＥＭＧ 信号，通过对其进行模糊化处理，经过一 系列模糊规则推理，采用模糊隶属度函数来表示肌 肉激活程度，从而作为 Ｈｉｌｌ 肌肉模型的输入，再根据 个人差异化肌肉参数，输入到 Ｈｉｌｌ 三元素模型，进行 肌肉力的预测计算，具体流程如下。 图 ５ Ｈｉｌｌ 肌力预测模型流程图 Ｆｉｇ．５ Ｆｌｏｗ ｃｈａｒｔ ｏｆ Ｈｉｌｌ ｍｕｓｃｌｅ ｓｔｒｅｎｇｔｈ ｐｒｅｄｉｃｔｉｏｎ ｍｏｄｅｌ 模糊逻辑推理模型主要有 Ｍａｍｄａｎｉ 和 Ｓｕｇｅｎｏ 两种方法。 首先，将采样得到的 ＥＭＧ 原始信号进行 滤波和平滑处理；其次对 ＥＭＧ 原始信号进行滤波和 平滑处理，将其用模糊集进行表示。 然后建立模糊 推理规则，根据处理后的 ＥＭＧ 信号幅值，推理出作 为肌肉激活程度的隶属度函数。 然后结合肌束和肌 腱的长度，通过 Ｈｉｌｌ 肌力模型，从而计算出肌肉力。 在 Ｍａｔｌａｂ 中的模糊工具箱中设计了模糊推理规则 和隶属度函数，通过 Ｍａｔｌａｂ 编程，实现对 Ｈｉｌｌ 肌力 预测模型的计算。 计算出肌肉力之后，就可以从力矩的定义出发， 将肌肉力与各自的力臂相乘，之后对关节进行简单 的受力分析，得出所需要的关节力矩［４４⁃４６］ 。 肌肉力臂 ｒｉ (ｔ) 可通过对骨骼肌长度和形状求 导得到，公式如下： ｌｍｔ (ｔ) ＝ Ｃ０ ＋ Ｃ１ θｋｎｅｅ（ｔ） ＋ Ｃ２ θ ２ ｋｎｅｅ (ｔ) ｒｉ (ｔ) ＝ ∂ｌｍｔ (ｔ) ∂θｋｎｅｅ （３５） 代入后得到力臂 ｒｉ (ｔ) ，其中， Ｃ０ 、 Ｃ１ 和 Ｃ２ 均 为常数。 ｒｉ (ｔ) ＝ Ｃ１ ＋ ２Ｃ２ θｋｎｅｅ [ (ｔ) ] １８０ π （３６） 关节力矩 Ｍ(ｔ) 计算如下： Ｍ(ｔ) ＝ ∑ ｎ ｉ ＝ １ Ｍｉ (ｔ) － ∑ ｋ ｊ ＝ １ Ｍｊ (ｔ) （３７） 式中： Ｍｉ 表示所有与关节相连的主动肌力矩， Ｍｊ 表 示所有拮抗肌力矩，如下： Ｍｉ (ｔ) ＝ Ｆｉ (ｔ) ｒｉ (ｔ) Ｍｊ (ｔ) ＝ Ｆｊ (ｔ) ｒｊ (ｔ) （３８） ５ 下肢建模仿真技术手段 下肢建模的仿真验证是建模过程中的重要一 环，也因为真实值测量困难、参数复杂多样等因素一 直是建模过程中的难点。 目前国内外对下肢建模的 仿真验证主要采取以下几种方法。 孙棕檀等人将基于 Ｌａｇｒａｎｇｅ 方程或 Ｋａｎｅ 方程 等传统动力学分析的结果，与 Ｓｉｍｕｌｉｎｋ 等仿真软件 输出的结果进行比较分析，证明了传统力学分析应 ·５２４· 智 能 系 统 学 报 第 １０ 卷
第4期 邵明旭，等：人体下肢生物力学建模研究进展 ·525· 用在下肢建模具有一定的可行性4]」 LYU Houshan.Artificial Joint Surgery[M].Beijing:Sci- Wei Meng等将SVR黑箱训练计算出的关节力 ence Press,2001:1-799. 矩与实验中力矩检测装置反馈回的真实值进行比较 [3]MOEINZADEH M H,ENGIN A E,AKKAS N.Two-dimen- 分析，结果证明模型准确度较高，EMG信号与关节力 sional dynamic modeling of human knee joint[Journal of 矩匹配情况较好，模型具有较高的应用意义[8]。陈 Biomechanics,1983,16(4):253-264. [4]ROYER T,KOENIG M.Joint loading and bone mineral 贵亮等人使用实验数据经过BP神经网络训练，对模 density in persons with unilateral,trans-tibial amputation 型进行辨别，确定模型中的未知参数，并使用额外的 [J].Clinical Biomechanics,2005,20(10):1119-1125. 数据去验证模型的准确度，有一定参考价值9]。 [5]YANG C J,ZHANG J F,CHEN Y,et al.A review of exo- Nikos Karavas等从表面肌电信号出发，验证得 skeleton-type systems and their key technologies[J].Pro- 到的肌肉力或转矩与真实人体行走过程中记录下的 ceedings of the Institution of Mechanical Engineers,Part C: EMG和扭矩波形进行趋势比较和均方根误差分析， Journal of Mechanical Engineering Science,2008,222(8): 证明模型的精准度较高，对仿真验证算法有着较高 1599-1612. 的指导意义[s0]。 [6]KAZEROONI H.Human augmentation and exoskeleton sys- tems in BerkeleyJ.International Journal of Humanoid Ro- 6结论 botics,2007.4(3):575-605. [7BRAND R A.The biomechanics and motor control of human 人体是一个非常复杂的多体模型。国内外研究 gait:normal,elderly,and pathological[M].Waterloo:U- 者针对人体不同的运动或动作，对人体做出了相应 niversity of Waterloo Press,1991:14-15,1-384. 的多刚体简化模型，并使用不同的建模方法建立了 [8]TUMER S T,ENGIN A E.Three-body segment dynamic 人体运动的动力学方程，取得了一些有意义的成果。 model of the human knee J.Journal of Biomechanical En- 但将人体简化为多刚体模型，忽略了肌肉和神经系 gineering,1993,115(4A):350-356. [9]MOEINZADEH M H.ENGIN A E.AKKAS N.Two-dimen- 统的作用，模拟的人体运动系统缺乏灵活性和自适 sional dynamic modeling of human knee joint[J.Journal of 应性。因此，所建的模型与人体实际运动情况具有 Biomechanics,1982,15(4):346. 较大差异。 [10]AYOUB MM.A 2-D simulation model for lifting activities SVR和仿真软件建模法也各有利弊。SVM模 [J].Computers Industrial Engineering,1998,35(3- 型法的优点是模型简单，但核心参数的选取还没有 4):619-622. 统一的理论公式，准确度有待提高。仿真软件建模 [11]PANDY M G,ZAJAC F E,SIM E,et al.An optimal con- 优点是下肢参数考虑较为完全，但缺点是数据需要 trol model for maximum-height human jumping[J].Journal 引出，并且计算量较大不利于实时控制。 of Biomechanics,1990.23(12):1185-1198. 相比较而言基于H山三元素模型建立的人体下 [12]SPAGELE T,KISTNER A,GOLLHOFER A.Modelling, simulation and optimisation of a human vertical jump[J]. 肢模型更为贴近人体的实际下肢结构，对人体运动 Journal of Biomechanics,1999,32(5):521-530. 的建模和仿真更符合人体实际运动状态。由于人体 [13]Anolg Devices,Inc.ADXL203 Data Sheet DB/OL]. 的运动过程非常复杂，而且每一时刻都会有不同的 [2006-04-27].http://www.analog.com/UploadedFiles/ 状态，应该根据不同的需求建立不同的模型。为了 Data_Sheets/ADXL103_203.pdf. 得到更加符合实际的模型，我们需要在现有建模技 [14]FAN Yuanjie,YIN Yuehong.Active and progressive exo- 术的基础上，更加深入的考虑人体的生理结构，将外 skeleton rehabilitation using multisource information fusion 骨骼与操作者的行走习惯有机结合在一起，兼顾人 from EMG and force-position EPP[J].IEEE Transactions 的意愿和外界作用力的影响，采用合理的控制策略， on Biomedical Engineering,2013,60(12):3314-3321. [15]HANAVAN E P.A mathematical model of the human 提高系统的鲁棒性，达到最优的控制目的。 body.AD608463[R].[S.1],1964. 参考文献： [16]BLAJER W,CZAPLICKI A.Modeling and inverse simula- tion of somersaults on the trampoline[J].Journal of Biome- [1]DOLLAR A M,HERR H.Lower extremity exoskeletons and chanics,2001,34(12):1619-1629. active orthoses:challenges and state-of-the-art[J].IEEE [17]ANDERSON F C,PANDY M G.Dynamic optimization of Transactions on Robotics,2008,24(1):144-158. human walking[J].Journal of Biomechanical Engineering, [2]吕厚山.人工关节外科学[M].北京：科学出版社， 2001,123(5):381-390. 2001:1-799 [18]PEJHAN S,FARAHMAND F.PARNIANPOUR M.De-
用在下肢建模具有一定的可行性［ ４７ ］ 。 Ｗｅｉ Ｍｅｎｇ 等将 ＳＶＲ 黑箱训练计算出的关节力 矩与实验中力矩检测装置反馈回的真实值进行比较 分析，结果证明模型准确度较高，ＥＭＧ 信号与关节力 矩匹配情况较好，模型具有较高的应用意义［ ４８ ］ 。 陈 贵亮等人使用实验数据经过 ＢＰ 神经网络训练，对模 型进行辨别，确定模型中的未知参数，并使用额外的 数据去验证模型的准确度，有一定参考价值［ ４９ ］ 。 Ｎｉｋｏｓ Ｋａｒａｖａｓ 等从表面肌电信号出发，验证得 到的肌肉力或转矩与真实人体行走过程中记录下的 ＥＭＧ 和扭矩波形进行趋势比较和均方根误差分析， 证明模型的精准度较高，对仿真验证算法有着较高 的指导意义［５０］ 。 ６ 结论 人体是一个非常复杂的多体模型。 国内外研究 者针对人体不同的运动或动作，对人体做出了相应 的多刚体简化模型，并使用不同的建模方法建立了 人体运动的动力学方程，取得了一些有意义的成果。 但将人体简化为多刚体模型，忽略了肌肉和神经系 统的作用，模拟的人体运动系统缺乏灵活性和自适 应性。 因此，所建的模型与人体实际运动情况具有 较大差异。 ＳＶＲ 和仿真软件建模法也各有利弊。 ＳＶＭ 模 型法的优点是模型简单，但核心参数的选取还没有 统一的理论公式，准确度有待提高。 仿真软件建模 优点是下肢参数考虑较为完全，但缺点是数据需要 引出，并且计算量较大不利于实时控制。 相比较而言基于 Ｈｉｌｌ 三元素模型建立的人体下 肢模型更为贴近人体的实际下肢结构，对人体运动 的建模和仿真更符合人体实际运动状态。 由于人体 的运动过程非常复杂，而且每一时刻都会有不同的 状态，应该根据不同的需求建立不同的模型。 为了 得到更加符合实际的模型，我们需要在现有建模技 术的基础上，更加深入的考虑人体的生理结构，将外 骨骼与操作者的行走习惯有机结合在一起，兼顾人 的意愿和外界作用力的影响，采用合理的控制策略， 提高系统的鲁棒性，达到最优的控制目的。 参考文献： ［１］ＤＯＬＬＡＲ Ａ Ｍ， ＨＥＲＲ Ｈ． Ｌｏｗｅｒ ｅｘｔｒｅｍｉｔｙ ｅｘｏｓｋｅｌｅｔｏｎｓ ａｎｄ ａｃｔｉｖｅ ｏｒｔｈｏｓｅｓ： ｃｈａｌｌｅｎｇｅｓ ａｎｄ ｓｔａｔｅ⁃ｏｆ⁃ｔｈｅ⁃ａｒｔ ［ Ｊ］． ＩＥＥＥ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｒｏｂｏｔｉｃｓ， ２００８， ２４（１）： １４４⁃１５８． ［２］吕厚山． 人工关节外科学 ［ Ｍ］． 北京： 科学出版社， ２００１： １⁃７９９． ＬＹＵ Ｈｏｕｓｈａｎ． Ａｒｔｉｆｉｃｉａｌ Ｊｏｉｎｔ Ｓｕｒｇｅｒｙ ［ Ｍ］． Ｂｅｉｊｉｎｇ： Ｓｃｉ⁃ ｅｎｃｅ Ｐｒｅｓｓ， ２００１： １⁃７９９． ［３］ＭＯＥＩＮＺＡＤＥＨ Ｍ Ｈ， ＥＮＧＩＮ Ａ Ｅ， ＡＫＫＡＳ Ｎ． Ｔｗｏ⁃ｄｉｍｅｎ⁃ ｓｉｏｎａｌ ｄｙｎａｍｉｃ ｍｏｄｅｌｉｎｇ ｏｆ ｈｕｍａｎ ｋｎｅｅ ｊｏｉｎｔ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃｓ， １９８３， １６（４）： ２５３⁃２６４． ［４］ ＲＯＹＥＲ Ｔ， ＫＯＥＮＩＧ Ｍ． Ｊｏｉｎｔ ｌｏａｄｉｎｇ ａｎｄ ｂｏｎｅ ｍｉｎｅｒａｌ ｄｅｎｓｉｔｙ ｉｎ ｐｅｒｓｏｎｓ ｗｉｔｈ ｕｎｉｌａｔｅｒａｌ， ｔｒａｎｓ⁃ｔｉｂｉａｌ ａｍｐｕｔａｔｉｏｎ ［Ｊ］． Ｃｌｉｎｉｃａｌ Ｂｉｏｍｅｃｈａｎｉｃｓ， ２００５， ２０（１０）： １１１９⁃１１２５． ［５］ＹＡＮＧ Ｃ Ｊ， ＺＨＡＮＧ Ｊ Ｆ， ＣＨＥＮ Ｙ， ｅｔ ａｌ． Ａ ｒｅｖｉｅｗ ｏｆ ｅｘｏ⁃ ｓｋｅｌｅｔｏｎ⁃ｔｙｐｅ ｓｙｓｔｅｍｓ ａｎｄ ｔｈｅｉｒ ｋｅｙ ｔｅｃｈｎｏｌｏｇｉｅｓ ［ Ｊ］． Ｐｒｏ⁃ ｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ Ｉｎｓｔｉｔｕｔｉｏｎ ｏｆ Ｍｅｃｈａｎｉｃａｌ Ｅｎｇｉｎｅｅｒｓ， Ｐａｒｔ Ｃ： Ｊｏｕｒｎａｌ ｏｆ Ｍｅｃｈａｎｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｓｃｉｅｎｃｅ， ２００８， ２２２（８）： １５９９⁃１６１２． ［６］ＫＡＺＥＲＯＯＮＩ Ｈ． Ｈｕｍａｎ ａｕｇｍｅｎｔａｔｉｏｎ ａｎｄ ｅｘｏｓｋｅｌｅｔｏｎ ｓｙｓ⁃ ｔｅｍｓ ｉｎ Ｂｅｒｋｅｌｅｙ［Ｊ］． Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｈｕｍａｎｏｉｄ Ｒｏ⁃ ｂｏｔｉｃｓ， ２００７， ４（３）： ５７５⁃６０５． ［７］ＢＲＡＮＤ Ｒ Ａ． Ｔｈｅ ｂｉｏｍｅｃｈａｎｉｃｓ ａｎｄ ｍｏｔｏｒ ｃｏｎｔｒｏｌ ｏｆ ｈｕｍａｎ ｇａｉｔ： ｎｏｒｍａｌ， ｅｌｄｅｒｌｙ， ａｎｄ ｐａｔｈｏｌｏｇｉｃａｌ［Ｍ］． Ｗａｔｅｒｌｏｏ： Ｕ⁃ ｎｉｖｅｒｓｉｔｙ ｏｆ Ｗａｔｅｒｌｏｏ Ｐｒｅｓｓ， １９９１： １４⁃１５， １⁃３８４． ［８］ ＴＵＭＥＲ Ｓ Ｔ， ＥＮＧＩＮ Ａ Ｅ． Ｔｈｒｅｅ⁃ｂｏｄｙ ｓｅｇｍｅｎｔ ｄｙｎａｍｉｃ ｍｏｄｅｌ ｏｆ ｔｈｅ ｈｕｍａｎ ｋｎｅｅ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃａｌ Ｅｎ⁃ ｇｉｎｅｅｒｉｎｇ， １９９３， １１５（４Ａ）： ３５０⁃３５６． ［９］ＭＯＥＩＮＺＡＤＥＨ Ｍ Ｈ， ＥＮＧＩＮ Ａ Ｅ， ＡＫＫＡＳ Ｎ． Ｔｗｏ⁃ｄｉｍｅｎ⁃ ｓｉｏｎａｌ ｄｙｎａｍｉｃ ｍｏｄｅｌｉｎｇ ｏｆ ｈｕｍａｎ ｋｎｅｅ ｊｏｉｎｔ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃｓ， １９８２， １５（４）： ３４６． ［１０］ＡＹＯＵＢ Ｍ Ｍ． Ａ ２⁃Ｄ ｓｉｍｕｌａｔｉｏｎ ｍｏｄｅｌ ｆｏｒ ｌｉｆｔｉｎｇ ａｃｔｉｖｉｔｉｅｓ ［Ｊ］． Ｃｏｍｐｕｔｅｒｓ ＆ Ｉｎｄｕｓｔｒｉａｌ Ｅｎｇｉｎｅｅｒｉｎｇ， １９９８， ３５ （ ３⁃ ４）： ６１９⁃６２２． ［１１］ＰＡＮＤＹ Ｍ Ｇ， ＺＡＪＡＣ Ｆ Ｅ， ＳＩＭ Ｅ， ｅｔ ａｌ． Ａｎ ｏｐｔｉｍａｌ ｃｏｎ⁃ ｔｒｏｌ ｍｏｄｅｌ ｆｏｒ ｍａｘｉｍｕｍ⁃ｈｅｉｇｈｔ ｈｕｍａｎ ｊｕｍｐｉｎｇ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃｓ， １９９０， ２３（１２）： １１８５⁃１１９８． ［１２］ＳＰÄＧＥＬＥ Ｔ， ＫＩＳＴＮＥＲ Ａ， ＧＯＬＬＨＯＦＥＲ Ａ． Ｍｏｄｅｌｌｉｎｇ， ｓｉｍｕｌａｔｉｏｎ ａｎｄ ｏｐｔｉｍｉｓａｔｉｏｎ ｏｆ ａ ｈｕｍａｎ ｖｅｒｔｉｃａｌ ｊｕｍｐ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃｓ， １９９９， ３２（５）： ５２１⁃５３０． ［１３ ］ Ａｎｏｌｇ Ｄｅｖｉｃｅｓ， Ｉｎｃ． ＡＤＸＬ２０３ Ｄａｔａ Ｓｈｅｅｔ ［ ＤＢ／ ＯＬ］． ［２００６⁃０４⁃２７］． ｈｔｔｐ： ／ ／ ｗｗｗ． ａｎａｌｏｇ． ｃｏｍ／ ＵｐｌｏａｄｅｄＦｉｌｅｓ／ Ｄａｔａ＿Ｓｈｅｅｔｓ／ ＡＤＸＬ１０３＿２０３．ｐｄｆ． ［１４］ ＦＡＮ Ｙｕａｎｊｉｅ， ＹＩＮ Ｙｕｅｈｏｎｇ． Ａｃｔｉｖｅ ａｎｄ ｐｒｏｇｒｅｓｓｉｖｅ ｅｘｏ⁃ ｓｋｅｌｅｔｏｎ ｒｅｈａｂｉｌｉｔａｔｉｏｎ ｕｓｉｎｇ ｍｕｌｔｉｓｏｕｒｃｅ ｉｎｆｏｒｍａｔｉｏｎ ｆｕｓｉｏｎ ｆｒｏｍ ＥＭＧ ａｎｄ ｆｏｒｃｅ⁃ｐｏｓｉｔｉｏｎ ＥＰＰ［ Ｊ］． ＩＥＥＥ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｂｉｏｍｅｄｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ， ２０１３， ６０（１２）： ３３１４⁃３３２１． ［１５］ ＨＡＮＡＶＡＮ Ｅ Ｐ． Ａ ｍａｔｈｅｍａｔｉｃａｌ ｍｏｄｅｌ ｏｆ ｔｈｅ ｈｕｍａｎ ｂｏｄｙ．ＡＤ６０８４６３［Ｒ］． ［Ｓ．ｌ．］， １９６４． ［１６］ＢＬＡＪＥＲ Ｗ， ＣＺＡＰＬＩＣＫＩ Ａ． Ｍｏｄｅｌｉｎｇ ａｎｄ ｉｎｖｅｒｓｅ ｓｉｍｕｌａ⁃ ｔｉｏｎ ｏｆ ｓｏｍｅｒｓａｕｌｔｓ ｏｎ ｔｈｅ ｔｒａｍｐｏｌｉｎｅ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅ⁃ ｃｈａｎｉｃｓ， ２００１， ３４（１２）： １６１９⁃１６２９． ［１７］ＡＮＤＥＲＳＯＮ Ｆ Ｃ， ＰＡＮＤＹ Ｍ Ｇ． Ｄｙｎａｍｉｃ ｏｐｔｉｍｉｚａｔｉｏｎ ｏｆ ｈｕｍａｎ ｗａｌｋｉｎｇ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ， ２００１， １２３（５）： ３８１⁃３９０． ［１８］ ＰＥＪＨＡＮ Ｓ， ＦＡＲＡＨＭＡＮＤ Ｆ， ＰＡＲＮＩＡＮＰＯＵＲ Ｍ． Ｄｅ⁃ 第 ４ 期 邵明旭，等：人体下肢生物力学建模研究进展 ·５２５·
.526. 智能系统学报 第10卷 sign optimization of an above-knee prosthesis based on the control of intelligent lower limb prosthesis using electro- kinematics of gait[C]//30th Annual International Confer- rheological fluid[J].Journal of Tsinghua University:Sci- ence of the IEEE Engineering in Medicine and Biology So- ence and Technology,1998,38(2):40-43. ciety.Vancouver,Canada,2008:4274-4277. [30]MARTINEZ-VILLALPANDO E C,MOONEY L,ELLIOTT [19]张国伟，宋伟刚.并联机器人动力学问题的Kane方法 G,et al.Antagonistic active knee prosthesis.A metabolic [J].系统仿真学报，2004,16(7)：1386-1391. cost of walking comparison with a variable-damping pros- ZHANG Guowei,SONG Weigang.A Kane formulation for thetic knee[C]//Proceedings of the 2011 Annual Interna- the inverse dynamic of Stewart platform manipulator[J]. tional Conference of the IEEE Engineering in Medicine and Journal of System Simulation,2004,16(7):1386-1391. Biology Society.Boston,USA,2011:8519-8522 [20]刘延斌，韩秀英，薛玉君，等.3-RRRT并联机器人动力 [31]王斌锐，金英连，徐心和仿生膝关节虚拟样机与协同仿 学仿真[J].系统仿真学报，2006,18(7)：1962-1965. 真方法研究[J].系统仿真学报，2006,18(6)：1554- LIU Yanbin,HAN Xiuying,XUE Yujun,et al.Dynamics 1557. simulation of a 3-RRRT parallel robot[J].Journal of Sys- WANG Binrui,JIN Yinglian,XU Xinhe.Virtual prototype tem Simulation,2006,18(7):1962-1965. and collaborative simulation of bionic knee joint[J].Jour- [21]TIBERT G.Deployable tenegrity structures for space appli- nal of System Simulation,2006,18(6):1554-1557. cation[R].Stockholm:Royal Institute of Technology, [32]朱昌义.单杠上人体摆动的凯恩动力学模型[J].成都 2002. 体育学院学报，2000,26(6)：71-74. [22]JENKYN T R,NICOL A C.A multi-segment kinematic ZHU Changyi.Kaine dynamics model of human swing on model of the foot with a novel definition of forefoot motion the horizontal bar Journal of Chengdu Sport Universi- for use in clinical gait analysis during walking[J].Journal ty,2000.26(6):71-74. of Biomechanics,2007,40(14):3271-3278. [33]刘明辉，顾文锦，陈占伏.基于骨骼服的虚拟人体建模 [23]SMITH S L.Application of high-speed videography in 与仿真[J].海军航空工程学院学报，2009,24(2)： sports analysis [C]//Proceedings of SPIE 1757,Ultra- 157-161. high-and High-Speed Photography,Videography,and LIU Minghui,GU Wenjin,CHEN Zhanfu.Virtual human Photonics.San Diego,USA,1992:118. body modeling and simulation based on skeletal services [24]NARUSE K,KAWAI S,KUKICHI T.Three-dimensional [J].Journal of Naval Aeronautical and Astronautical Uni- lifting-up motion analysis for wearable power assist device versity,2009,24(2):157-161. of lower back support[C]//IEEE/RSJ Intemational Con- [34]沈凌，孟青云，喻洪流.基于虚拟样机技术的下肢假肢 ference on Intelligent Robots and Systems.Edmonton,Can- 结构设计与仿真[J].工程设计学报，2011,18(1)：34- ada,2005:2959-2964. 38. [25]KIGUCHI K,HAYASHI Y.Estimation of joint torque for a SHEN Ling,MENG Qingyun,YU Hongliu.Design and myoelectric arm by genetic programming based on EMG simulation of leg prosthesis structure based on virtual proto- signals[C]//IEEE World Automation Congress (WAC). type technology[J].Journal of Engineering Design,2011, Puerto Vallarta,Mexico,2012:1-4. 18(1):34-38. [26SARTORI M,REGGIANI M,PAGELLO E,et al.Model- [35]MOURAGNON E,LHUILLIER M,DHOME M,et al.Mo- ing the human knee for assistive technologies[J].IEEE nocular vision based SLAM for mobile robots [C]//IEEE Transactions on Biomedical Engineering,2012,59(9): International Conference on Pattern Recognition.Hong 2642-2649. Kong,China,2006,3:1027-1031. [27]PATRINOS P,ALEXANDRIDIS A,NINOS K,et al.Var- [36]PRATT J,CHEW C M,TORRES A,et al.Virtual model iable selection in nonlinear modeling based on RBF net- control:An intuitive approach for bipedal locomotion[J]. works and evolutionary computation[J].Intemational Jour- The International Journal of Robotics Research,2001,20 nal of Neural Systems,2010,20(5):365-379. (2):129-143. [28 ]SONG R,TONG K Y.Using recurrent artificial neural net- [37]YUAN Shaoqiang,LIU Zhong,LI Xingshan.Modeling and work model to estimate voluntary elbow torque in dynamic simulation of robot based on Matlab/SimMechanics[C]// situations[J].Medical and Biological Engineering and Proceedings of the 27th Chinese Control Conference.Kun- Computing,2005,43(4):473-480. ming,China,2008:161-165. [29]金德闻，王人成，白彩勤.电流变液智能下肢假肢摆动 [38]XU Guozheng,SONG Aiguo,LI Huijun.Control system 相控制原理与方法[J]清华大学学报：自然科学版， design for an upper-limb rehabilitation robot[J].Advanced 1998,38(2):40-43. Robotics.,2011,25(1-2):229-251. JIN Dewen,WANG Rencheng,BAI Caiqin.Swing phase [39]RARNJANI A N,CORSS P A.A Kalman filter model of
ｓｉｇｎ ｏｐｔｉｍｉｚａｔｉｏｎ ｏｆ ａｎ ａｂｏｖｅ⁃ｋｎｅｅ ｐｒｏｓｔｈｅｓｉｓ ｂａｓｅｄ ｏｎ ｔｈｅ ｋｉｎｅｍａｔｉｃｓ ｏｆ ｇａｉｔ［Ｃ］ ／ ／ ３０ｔｈ Ａｎｎｕａｌ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒ⁃ ｅｎｃｅ ｏｆ ｔｈｅ ＩＥＥＥ Ｅｎｇｉｎｅｅｒｉｎｇ ｉｎ Ｍｅｄｉｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ Ｓｏ⁃ ｃｉｅｔｙ． Ｖａｎｃｏｕｖｅｒ， Ｃａｎａｄａ， ２００８： ４２７４⁃４２７７． ［１９］张国伟，宋伟刚．并联机器人动力学问题的 Ｋａｎｅ 方法 ［Ｊ］．系统仿真学报， ２００４， １６（７）： １３８６⁃１３９１． ＺＨＡＮＧ Ｇｕｏｗｅｉ， ＳＯＮＧ Ｗｅｉｇａｎｇ． Ａ Ｋａｎｅ ｆｏｒｍｕｌａｔｉｏｎ ｆｏｒ ｔｈｅ ｉｎｖｅｒｓｅ ｄｙｎａｍｉｃ ｏｆ Ｓｔｅｗａｒｔ ｐｌａｔｆｏｒｍ ｍａｎｉｐｕｌａｔｏｒ ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｓｙｓｔｅｍ Ｓｉｍｕｌａｔｉｏｎ， ２００４， １６（７）： １３８６⁃１３９１． ［２０］刘延斌， 韩秀英，薛玉君，等．３⁃ＲＲＲＴ 并联机器人动力 学仿真［Ｊ］．系统仿真学报， ２００６， １８（７）： １９６２⁃１９６５． ＬＩＵ Ｙａｎｂｉｎ， ＨＡＮ Ｘｉｕｙｉｎｇ， ＸＵＥ Ｙｕｊｕｎ， ｅｔ ａｌ． Ｄｙｎａｍｉｃｓ ｓｉｍｕｌａｔｉｏｎ ｏｆ ａ ３⁃ＲＲＲＴ ｐａｒａｌｌｅｌ ｒｏｂｏｔ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｓｙｓ⁃ ｔｅｍ Ｓｉｍｕｌａｔｉｏｎ， ２００６， １８（７）： １９６２⁃１９６５． ［２１］ＴＩＢＥＲＴ Ｇ． Ｄｅｐｌｏｙａｂｌｅ ｔｅｎｅｇｒｉｔｙ ｓｔｒｕｃｔｕｒｅｓ ｆｏｒ ｓｐａｃｅ ａｐｐｌｉ⁃ ｃａｔｉｏｎ ［ Ｒ］． Ｓｔｏｃｋｈｏｌｍ： Ｒｏｙａｌ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｔｅｃｈｎｏｌｏｇｙ， ２００２． ［２２］ ＪＥＮＫＹＮ Ｔ Ｒ， ＮＩＣＯＬ Ａ Ｃ． Ａ ｍｕｌｔｉ⁃ｓｅｇｍｅｎｔ ｋｉｎｅｍａｔｉｃ ｍｏｄｅｌ ｏｆ ｔｈｅ ｆｏｏｔ ｗｉｔｈ ａ ｎｏｖｅｌ ｄｅｆｉｎｉｔｉｏｎ ｏｆ ｆｏｒｅｆｏｏｔ ｍｏｔｉｏｎ ｆｏｒ ｕｓｅ ｉｎ ｃｌｉｎｉｃａｌ ｇａｉｔ ａｎａｌｙｓｉｓ ｄｕｒｉｎｇ ｗａｌｋｉｎｇ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃｓ， ２００７， ４０（１４）： ３２７１⁃３２７８． ［２３ ］ ＳＭＩＴＨ Ｓ Ｌ． Ａｐｐｌｉｃａｔｉｏｎ ｏｆ ｈｉｇｈ⁃ｓｐｅｅｄ ｖｉｄｅｏｇｒａｐｈｙ ｉｎ ｓｐｏｒｔｓ ａｎａｌｙｓｉｓ ［ Ｃ］ ／ ／ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＳＰＩＥ １７５７， Ｕｌｔｒａ⁃ ｈｉｇｈ⁃ ａｎｄ Ｈｉｇｈ⁃Ｓｐｅｅｄ Ｐｈｏｔｏｇｒａｐｈｙ， Ｖｉｄｅｏｇｒａｐｈｙ， ａｎｄ Ｐｈｏｔｏｎｉｃｓ． Ｓａｎ Ｄｉｅｇｏ， ＵＳＡ， １９９２： １１８． ［２４］ＮＡＲＵＳＥ Ｋ， ＫＡＷＡＩ Ｓ， ＫＵＫＩＣＨＩ Ｔ． Ｔｈｒｅｅ⁃ｄｉｍｅｎｓｉｏｎａｌ ｌｉｆｔｉｎｇ⁃ｕｐ ｍｏｔｉｏｎ ａｎａｌｙｓｉｓ ｆｏｒ ｗｅａｒａｂｌｅ ｐｏｗｅｒ ａｓｓｉｓｔ ｄｅｖｉｃｅ ｏｆ ｌｏｗｅｒ ｂａｃｋ ｓｕｐｐｏｒｔ［ Ｃ］ ／ ／ ＩＥＥＥ／ ＲＳＪ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎ⁃ ｆｅｒｅｎｃｅ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｒｏｂｏｔｓ ａｎｄ Ｓｙｓｔｅｍｓ． Ｅｄｍｏｎｔｏｎ， Ｃａｎ⁃ ａｄａ， ２００５： ２９５９⁃２９６４． ［２５］ＫＩＧＵＣＨＩ Ｋ， ＨＡＹＡＳＨＩ Ｙ． Ｅｓｔｉｍａｔｉｏｎ ｏｆ ｊｏｉｎｔ ｔｏｒｑｕｅ ｆｏｒ ａ ｍｙｏｅｌｅｃｔｒｉｃ ａｒｍ ｂｙ ｇｅｎｅｔｉｃ ｐｒｏｇｒａｍｍｉｎｇ ｂａｓｅｄ ｏｎ ＥＭＧ ｓｉｇｎａｌｓ［Ｃ］ ／ ／ ＩＥＥＥ Ｗｏｒｌｄ Ａｕｔｏｍａｔｉｏｎ Ｃｏｎｇｒｅｓｓ （ ＷＡＣ）． Ｐｕｅｒｔｏ Ｖａｌｌａｒｔａ， Ｍｅｘｉｃｏ， ２０１２： １⁃４． ［２６］ＳＡＲＴＯＲＩ Ｍ， ＲＥＧＧＩＡＮＩ Ｍ， ＰＡＧＥＬＬＯ Ｅ， ｅｔ ａｌ． Ｍｏｄｅｌ⁃ ｉｎｇ ｔｈｅ ｈｕｍａｎ ｋｎｅｅ ｆｏｒ ａｓｓｉｓｔｉｖｅ ｔｅｃｈｎｏｌｏｇｉｅｓ ［ Ｊ］． ＩＥＥＥ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｂｉｏｍｅｄｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ， ２０１２， ５９ （ ９）： ２６４２⁃２６４９． ［２７］ＰＡＴＲＩＮＯＳ Ｐ， ＡＬＥＸＡＮＤＲＩＤＩＳ Ａ， ＮＩＮＯＳ Ｋ， ｅｔ ａｌ． Ｖａｒ⁃ ｉａｂｌｅ ｓｅｌｅｃｔｉｏｎ ｉｎ ｎｏｎｌｉｎｅａｒ ｍｏｄｅｌｉｎｇ ｂａｓｅｄ ｏｎ ＲＢＦ ｎｅｔ⁃ ｗｏｒｋｓ ａｎｄ ｅｖｏｌｕｔｉｏｎａｒｙ ｃｏｍｐｕｔａｔｉｏｎ［Ｊ］． Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒ⁃ ｎａｌ ｏｆ Ｎｅｕｒａｌ Ｓｙｓｔｅｍｓ， ２０１０， ２０（５）： ３６５⁃３７９． ［２８］ＳＯＮＧ Ｒ， ＴＯＮＧ Ｋ Ｙ． Ｕｓｉｎｇ ｒｅｃｕｒｒｅｎｔ ａｒｔｉｆｉｃｉａｌ ｎｅｕｒａｌ ｎｅｔ⁃ ｗｏｒｋ ｍｏｄｅｌ ｔｏ ｅｓｔｉｍａｔｅ ｖｏｌｕｎｔａｒｙ ｅｌｂｏｗ ｔｏｒｑｕｅ ｉｎ ｄｙｎａｍｉｃ ｓｉｔｕａｔｉｏｎｓ ［ Ｊ ］． Ｍｅｄｉｃａｌ ａｎｄ Ｂｉｏｌｏｇｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ Ｃｏｍｐｕｔｉｎｇ， ２００５， ４３（４）： ４７３⁃４８０． ［２９］金德闻，王人成，白彩勤． 电流变液智能下肢假肢摆动 相控制原理与方法［ Ｊ］． 清华大学学报：自然科学版， １９９８， ３８（２）： ４０⁃４３． ＪＩＮ Ｄｅｗｅｎ， ＷＡＮＧ Ｒｅｎｃｈｅｎｇ， ＢＡＩ Ｃａｉｑｉｎ． Ｓｗｉｎｇ ｐｈａｓｅ ｃｏｎｔｒｏｌ ｏｆ ｉｎｔｅｌｌｉｇｅｎｔ ｌｏｗｅｒ ｌｉｍｂ ｐｒｏｓｔｈｅｓｉｓ ｕｓｉｎｇ ｅｌｅｃｔｒｏ⁃ ｒｈｅｏｌｏｇｉｃａｌ ｆｌｕｉｄ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｔｓｉｎｇｈｕａ Ｕｎｉｖｅｒｓｉｔｙ： Ｓｃｉ⁃ ｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ， １９９８， ３８（２）： ４０⁃４３． ［３０］ＭＡＲＴＩＮＥＺ⁃ＶＩＬＬＡＬＰＡＮＤＯ Ｅ Ｃ， ＭＯＯＮＥＹ Ｌ， ＥＬＬＩＯＴＴ Ｇ， ｅｔ ａｌ． Ａｎｔａｇｏｎｉｓｔｉｃ ａｃｔｉｖｅ ｋｎｅｅ ｐｒｏｓｔｈｅｓｉｓ． Ａ ｍｅｔａｂｏｌｉｃ ｃｏｓｔ ｏｆ ｗａｌｋｉｎｇ ｃｏｍｐａｒｉｓｏｎ ｗｉｔｈ ａ ｖａｒｉａｂｌｅ⁃ｄａｍｐｉｎｇ ｐｒｏｓ⁃ ｔｈｅｔｉｃ ｋｎｅｅ［Ｃ］ ／ ／ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ２０１１ Ａｎｎｕａｌ Ｉｎｔｅｒｎａ⁃ ｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｆ ｔｈｅ ＩＥＥＥ Ｅｎｇｉｎｅｅｒｉｎｇ ｉｎ Ｍｅｄｉｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ Ｓｏｃｉｅｔｙ． Ｂｏｓｔｏｎ， ＵＳＡ， ２０１１： ８５１９⁃８５２２． ［３１］王斌锐，金英连，徐心和．仿生膝关节虚拟样机与协同仿 真方法研究［ Ｊ］．系统仿真学报， ２００６， １８（ ６）： １５５４⁃ １５５７． ＷＡＮＧ Ｂｉｎｒｕｉ， ＪＩＮ Ｙｉｎｇｌｉａｎ， ＸＵ Ｘｉｎｈｅ． Ｖｉｒｔｕａｌ ｐｒｏｔｏｔｙｐｅ ａｎｄ ｃｏｌｌａｂｏｒａｔｉｖｅ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｂｉｏｎｉｃ ｋｎｅｅ ｊｏｉｎｔ［ Ｊ］． Ｊｏｕｒ⁃ ｎａｌ ｏｆ Ｓｙｓｔｅｍ Ｓｉｍｕｌａｔｉｏｎ， ２００６， １８（６）： １５５４⁃１５５７． ［３２］朱昌义． 单杠上人体摆动的凯恩动力学模型［ Ｊ］． 成都 体育学院学报， ２０００， ２６（６）： ７１⁃７４． ＺＨＵ Ｃｈａｎｇｙｉ． Ｋａｉｎｅ ｄｙｎａｍｉｃｓ ｍｏｄｅｌ ｏｆ ｈｕｍａｎ ｓｗｉｎｇ ｏｎ ｔｈｅ ｈｏｒｉｚｏｎｔａｌ ｂａｒ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｃｈｅｎｇｄｕ Ｓｐｏｒｔ Ｕｎｉｖｅｒｓｉ⁃ ｔｙ， ２０００， ２６（６）： ７１⁃７４． ［３３］刘明辉， 顾文锦， 陈占伏． 基于骨骼服的虚拟人体建模 与仿真［Ｊ］． 海军航空工程学院学报， ２００９， ２４ （ ２）： １５７⁃１６１． ＬＩＵ Ｍｉｎｇｈｕｉ， ＧＵ Ｗｅｎｊｉｎ， ＣＨＥＮ Ｚｈａｎｆｕ． Ｖｉｒｔｕａｌ ｈｕｍａｎ ｂｏｄｙ ｍｏｄｅｌｉｎｇ ａｎｄ ｓｉｍｕｌａｔｉｏｎ ｂａｓｅｄ ｏｎ ｓｋｅｌｅｔａｌ ｓｅｒｖｉｃｅｓ ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｎａｖａｌ Ａｅｒｏｎａｕｔｉｃａｌ ａｎｄ Ａｓｔｒｏｎａｕｔｉｃａｌ Ｕｎｉ⁃ ｖｅｒｓｉｔｙ， ２００９， ２４（２）： １５７⁃１６１． ［３４］沈凌， 孟青云， 喻洪流． 基于虚拟样机技术的下肢假肢 结构设计与仿真［Ｊ］． 工程设计学报， ２０１１， １８（１）： ３４⁃ ３８． ＳＨＥＮ Ｌｉｎｇ， ＭＥＮＧ Ｑｉｎｇｙｕｎ， ＹＵ Ｈｏｎｇｌｉｕ． Ｄｅｓｉｇｎ ａｎｄ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｌｅｇ ｐｒｏｓｔｈｅｓｉｓ ｓｔｒｕｃｔｕｒｅ ｂａｓｅｄ ｏｎ ｖｉｒｔｕａｌ ｐｒｏｔｏ⁃ ｔｙｐｅ ｔｅｃｈｎｏｌｏｇｙ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｅｎｇｉｎｅｅｒｉｎｇ Ｄｅｓｉｇｎ， ２０１１， １８（１）： ３４⁃３８． ［３５］ＭＯＵＲＡＧＮＯＮ Ｅ， ＬＨＵＩＬＬＩＥＲ Ｍ， ＤＨＯＭＥ Ｍ， ｅｔ ａｌ． Ｍｏ⁃ ｎｏｃｕｌａｒ ｖｉｓｉｏｎ ｂａｓｅｄ ＳＬＡＭ ｆｏｒ ｍｏｂｉｌｅ ｒｏｂｏｔｓ［Ｃ］ ／ ／ ＩＥＥＥ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｐａｔｔｅｒｎ Ｒｅｃｏｇｎｉｔｉｏｎ． Ｈｏｎｇ Ｋｏｎｇ， Ｃｈｉｎａ， ２００６， ３： １０２７⁃１０３１． ［３６］ＰＲＡＴＴ Ｊ， ＣＨＥＷ Ｃ Ｍ， ＴＯＲＲＥＳ Ａ， ｅｔ ａｌ． Ｖｉｒｔｕａｌ ｍｏｄｅｌ ｃｏｎｔｒｏｌ： Ａｎ ｉｎｔｕｉｔｉｖｅ ａｐｐｒｏａｃｈ ｆｏｒ ｂｉｐｅｄａｌ ｌｏｃｏｍｏｔｉｏｎ［ Ｊ］． Ｔｈｅ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｒｏｂｏｔｉｃｓ Ｒｅｓｅａｒｃｈ， ２００１， ２０ （２）： １２９⁃１４３． ［３７］ＹＵＡＮ Ｓｈａｏｑｉａｎｇ， ＬＩＵ Ｚｈｏｎｇ， ＬＩ Ｘｉｎｇｓｈａｎ． Ｍｏｄｅｌｉｎｇ ａｎｄ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｒｏｂｏｔ ｂａｓｅｄ ｏｎ Ｍａｔｌａｂ ／ ＳｉｍＭｅｃｈａｎｉｃｓ［Ｃ］ ／ ／ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ２７ｔｈ Ｃｈｉｎｅｓｅ Ｃｏｎｔｒｏｌ Ｃｏｎｆｅｒｅｎｃｅ． Ｋｕｎ⁃ ｍｉｎｇ， Ｃｈｉｎａ， ２００８： １６１⁃１６５． ［３８］ ＸＵ Ｇｕｏｚｈｅｎｇ， ＳＯＮＧ Ａｉｇｕｏ， ＬＩ Ｈｕｉｊｕｎ． Ｃｏｎｔｒｏｌ ｓｙｓｔｅｍ ｄｅｓｉｇｎ ｆｏｒ ａｎ ｕｐｐｅｒ⁃ｌｉｍｂ ｒｅｈａｂｉｌｉｔａｔｉｏｎ ｒｏｂｏｔ［Ｊ］． Ａｄｖａｎｃｅｄ Ｒｏｂｏｔｉｃｓ， ２０１１， ２５（１⁃２）： ２２９⁃２５１． ［３９］ＲＡＲＮＪＡＮＩ Ａ Ｎ， ＣＯＲＳＳ Ｐ Ａ． Ａ Ｋａｌｍａｎ ｆｉｌｔｅｒ ｍｏｄｅｌ ｏｆ ·５２６· 智 能 系 统 学 报 第 １０ 卷
第4期 邵明旭，等：人体下肢生物力学建模研究进展 ·527. an integrated land vehicle navigation system[J].Journal of 48]MENG Wei,DING Bo,ZHOU Zude,et al.An EMG- Navigation,1995,48(2):293-302. based force prediction and control approach for robot-assis- [40]HARDYK A TT.Force and power-velocity relationships in ted lower limb rehabilitation[J].IEEE International Con- a multi-joint movement[D].Pennsylvania State:The ference on Systems,Man,and Cybernetics.San Diego. Pennsylvania State University,2000:102-123. USA,2014:2198-2203. [41]FENN W O,MARSH B S.Muscular force at different [49]陈贵亮，李长鹏，赵月，等.下肢外骨酪康复机器人的 speeds of shortening[J].The Journal of Physiology,1935, 动力学建模及神经网络辨识仿真[J].机械设计与制 85(3):277-297. 造，2013，(11)：197-200. [42]POLISSAR M J.Physical chemistry of contractile process CHEN Guiliang,LI Changpeng,ZHAO Yue,et al.Dy- in muscle.I.A physicochemical model of contractile mech- namic modeling and neural network identification simula- anism[J].The American Journal of Physiology,1952, tion for lower limbs exoskeletons rehabilitation robot[]. 168(3):766-811. Mechanical Design and Manufacture,2013,(11):197- [43]SHADMEHR R,ARBIB M.A mathematical analysis of the 200. force-stiffness characteristics of muscles in control of a sin- [50]KARAVAS N,AJOUDANI A,TSAGARAKIS N,et al. gle joint system[J].Biological Cybernetics,1992,66 Tele-impedance based assistive control for a compliant (6):463-477. knee exoskeleton[J].Robotics and Autonomous Systems, [44]VISSER JJ,HOOGKAMER J E,BOBBERT M F,et al. 2014,doi:10.1016/j.robot..2014.09.027. Length and moment arm of human leg muscles as a function 作者简介： of knee and hip-joint angles[J].European Journal of Ap- 邵明旭，男，1993年生，硕士研究 plied Physiology and Occupational Physiology,1990,61 生，主要研究方向为模式识别与智能机 (5-6)：453-460. 器人。 [45]AN K N,TAKAHASHI K,HARRIGAN T P,et al.Deter- mination of muscle orientations and moment arms [J]. Journal of Biomechanical Engineering,1984,106(3): 280-282. [46]OSU R,FRANKLIN D W,KATO H,et al.Short-and 王斐，男，1974年生，副教授，博士， long-term changes in joint co-contraction associated with 主要研究方向为模式识别与智能机器 motor learning as revealed from surface EMG[J].Journal 人。完成包括国家自然科学基金，辽宁 of Neurophysiology,2002,88(2):991-1004. 省自然科学基金、教育部博士点基金和 [47]孙棕檀.刚柔耦合系统分析动力学建模研究[D].哈尔 中国人民解放军总装备部预研项目在 滨：哈尔滨工程大学，2013：44-56 内的国家和省部级科研项目15项，获 SUN Zongtan.Research on dynamic modeling analysis of 得辽宁省自然科学学术成果一等奖1项。发表学术论文50 rigid flexible coupling system[D].Harbin:Harbin Engi- 余篇。 neering University,2013:44-56. [责任编辑：郑可为]
ａｎ ｉｎｔｅｇｒａｔｅｄ ｌａｎｄ ｖｅｈｉｃｌｅ ｎａｖｉｇａｔｉｏｎ ｓｙｓｔｅｍ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｎａｖｉｇａｔｉｏｎ， １９９５， ４８（２）： ２９３⁃３０２． ［４０］ＨＡＲＤＹＫ Ａ Ｔ Ｔ． Ｆｏｒｃｅ ａｎｄ ｐｏｗｅｒ⁃ｖｅｌｏｃｉｔｙ ｒｅｌａｔｉｏｎｓｈｉｐｓ ｉｎ ａ ｍｕｌｔｉ⁃ｊｏｉｎｔ ｍｏｖｅｍｅｎｔ ［ Ｄ ］． Ｐｅｎｎｓｙｌｖａｎｉａ Ｓｔａｔｅ： Ｔｈｅ Ｐｅｎｎｓｙｌｖａｎｉａ Ｓｔａｔｅ Ｕｎｉｖｅｒｓｉｔｙ， ２０００： １０２⁃１２３． ［４１］ ＦＥＮＮ Ｗ Ｏ， ＭＡＲＳＨ Ｂ Ｓ． Ｍｕｓｃｕｌａｒ ｆｏｒｃｅ ａｔ ｄｉｆｆｅｒｅｎｔ ｓｐｅｅｄｓ ｏｆ ｓｈｏｒｔｅｎｉｎｇ［Ｊ］． Ｔｈｅ Ｊｏｕｒｎａｌ ｏｆ Ｐｈｙｓｉｏｌｏｇｙ， １９３５， ８５（３）： ２７７⁃２９７． ［４２］ ＰＯＬＩＳＳＡＲ Ｍ Ｊ． Ｐｈｙｓｉｃａｌ ｃｈｅｍｉｓｔｒｙ ｏｆ ｃｏｎｔｒａｃｔｉｌｅ ｐｒｏｃｅｓｓ ｉｎ ｍｕｓｃｌｅ． Ｉ． Ａ ｐｈｙｓｉｃｏｃｈｅｍｉｃａｌ ｍｏｄｅｌ ｏｆ ｃｏｎｔｒａｃｔｉｌｅ ｍｅｃｈ⁃ ａｎｉｓｍ［ Ｊ］． Ｔｈｅ Ａｍｅｒｉｃａｎ Ｊｏｕｒｎａｌ ｏｆ Ｐｈｙｓｉｏｌｏｇｙ， １９５２， １６８（３）： ７６６⁃８１１． ［４３］ＳＨＡＤＭＥＨＲ Ｒ， ＡＲＢＩＢ Ｍ． Ａ ｍａｔｈｅｍａｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ ｆｏｒｃｅ⁃ｓｔｉｆｆｎｅｓｓ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｍｕｓｃｌｅｓ ｉｎ ｃｏｎｔｒｏｌ ｏｆ ａ ｓｉｎ⁃ ｇｌｅ ｊｏｉｎｔ ｓｙｓｔｅｍ ［ Ｊ ］． Ｂｉｏｌｏｇｉｃａｌ Ｃｙｂｅｒｎｅｔｉｃｓ， １９９２， ６６ （６）： ４６３⁃４７７． ［４４］ＶＩＳＳＥＲ Ｊ Ｊ， ＨＯＯＧＫＡＭＥＲ Ｊ Ｅ， ＢＯＢＢＥＲＴ Ｍ Ｆ， ｅｔ ａｌ． Ｌｅｎｇｔｈ ａｎｄ ｍｏｍｅｎｔ ａｒｍ ｏｆ ｈｕｍａｎ ｌｅｇ ｍｕｓｃｌｅｓ ａｓ ａ ｆｕｎｃｔｉｏｎ ｏｆ ｋｎｅｅ ａｎｄ ｈｉｐ⁃ｊｏｉｎｔ ａｎｇｌｅｓ［ Ｊ］． Ｅｕｒｏｐｅａｎ Ｊｏｕｒｎａｌ ｏｆ Ａｐ⁃ ｐｌｉｅｄ Ｐｈｙｓｉｏｌｏｇｙ ａｎｄ Ｏｃｃｕｐａｔｉｏｎａｌ Ｐｈｙｓｉｏｌｏｇｙ， １９９０， ６１ （５⁃６）： ４５３⁃４６０． ［４５］ＡＮ Ｋ Ｎ， ＴＡＫＡＨＡＳＨＩ Ｋ， ＨＡＲＲＩＧＡＮ Ｔ Ｐ， ｅｔ ａｌ． Ｄｅｔｅｒ⁃ ｍｉｎａｔｉｏｎ ｏｆ ｍｕｓｃｌｅ ｏｒｉｅｎｔａｔｉｏｎｓ ａｎｄ ｍｏｍｅｎｔ ａｒｍｓ ［ Ｊ ］． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｃｈａｎｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ， １９８４， １０６ （ ３ ）： ２８０⁃２８２． ［４６］ ＯＳＵ Ｒ， ＦＲＡＮＫＬＩＮ Ｄ Ｗ， ＫＡＴＯ Ｈ， ｅｔ ａｌ． Ｓｈｏｒｔ⁃ａｎｄ ｌｏｎｇ⁃ｔｅｒｍ ｃｈａｎｇｅｓ ｉｎ ｊｏｉｎｔ ｃｏ⁃ｃｏｎｔｒａｃｔｉｏｎ ａｓｓｏｃｉａｔｅｄ ｗｉｔｈ ｍｏｔｏｒ ｌｅａｒｎｉｎｇ ａｓ ｒｅｖｅａｌｅｄ ｆｒｏｍ ｓｕｒｆａｃｅ ＥＭＧ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｎｅｕｒｏｐｈｙｓｉｏｌｏｇｙ， ２００２， ８８（２）： ９９１⁃１００４． ［４７］孙棕檀． 刚柔耦合系统分析动力学建模研究［Ｄ］． 哈尔 滨： 哈尔滨工程大学， ２０１３： ４４⁃５６． ＳＵＮ Ｚｏｎｇｔａｎ． Ｒｅｓｅａｒｃｈ ｏｎ ｄｙｎａｍｉｃ ｍｏｄｅｌｉｎｇ ａｎａｌｙｓｉｓ ｏｆ ｒｉｇｉｄ ｆｌｅｘｉｂｌｅ ｃｏｕｐｌｉｎｇ ｓｙｓｔｅｍ［Ｄ］． Ｈａｒｂｉｎ： Ｈａｒｂｉｎ Ｅｎｇｉ⁃ ｎｅｅｒｉｎｇ Ｕｎｉｖｅｒｓｉｔｙ， ２０１３： ４４⁃５６． ［４８］ ＭＥＮＧ Ｗｅｉ， ＤＩＮＧ Ｂｏ， ＺＨＯＵ Ｚｕｄｅ， ｅｔ ａｌ． Ａｎ ＥＭＧ⁃ ｂａｓｅｄ ｆｏｒｃｅ ｐｒｅｄｉｃｔｉｏｎ ａｎｄ ｃｏｎｔｒｏｌ ａｐｐｒｏａｃｈ ｆｏｒ ｒｏｂｏｔ⁃ａｓｓｉｓ⁃ ｔｅｄ ｌｏｗｅｒ ｌｉｍｂ ｒｅｈａｂｉｌｉｔａｔｉｏｎ［ Ｊ］． ＩＥＥＥ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎ⁃ ｆｅｒｅｎｃｅ ｏｎ Ｓｙｓｔｅｍｓ， Ｍａｎ， ａｎｄ Ｃｙｂｅｒｎｅｔｉｃｓ． Ｓａｎ Ｄｉｅｇｏ， ＵＳＡ， ２０１４： ２１９８⁃２２０３． ［４９］陈贵亮， 李长鹏， 赵月， 等． 下肢外骨骼康复机器人的 动力学建模及神经网络辨识仿真［ Ｊ］． 机械设计与制 造， ２０１３， （１１）： １９７⁃２００． ＣＨＥＮ Ｇｕｉｌｉａｎｇ， ＬＩ Ｃｈａｎｇｐｅｎｇ， ＺＨＡＯ Ｙｕｅ， ｅｔ ａｌ． Ｄｙ⁃ ｎａｍｉｃ ｍｏｄｅｌｉｎｇ ａｎｄ ｎｅｕｒａｌ ｎｅｔｗｏｒｋ ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｓｉｍｕｌａ⁃ ｔｉｏｎ ｆｏｒ ｌｏｗｅｒ ｌｉｍｂｓ ｅｘｏｓｋｅｌｅｔｏｎｓ ｒｅｈａｂｉｌｉｔａｔｉｏｎ ｒｏｂｏｔ ［ Ｊ］． Ｍｅｃｈａｎｉｃａｌ Ｄｅｓｉｇｎ ａｎｄ Ｍａｎｕｆａｃｔｕｒｅ， ２０１３， （ １１）： １９７⁃ ２００． ［５０］ ＫＡＲＡＶＡＳ Ｎ， ＡＪＯＵＤＡＮＩ Ａ， ＴＳＡＧＡＲＡＫＩＳ Ｎ， ｅｔ ａｌ． Ｔｅｌｅ⁃ｉｍｐｅｄａｎｃｅ ｂａｓｅｄ ａｓｓｉｓｔｉｖｅ ｃｏｎｔｒｏｌ ｆｏｒ ａ ｃｏｍｐｌｉａｎｔ ｋｎｅｅ ｅｘｏｓｋｅｌｅｔｏｎ［Ｊ］． Ｒｏｂｏｔｉｃｓ ａｎｄ Ａｕｔｏｎｏｍｏｕｓ Ｓｙｓｔｅｍｓ， ２０１４， ｄｏｉ： １０．１０１６ ／ ｊ．ｒｏｂｏｔ．２０１４．０９．０２７． 作者简介： 邵明旭，男，１９９３ 年生，硕士研究 生，主要研究方向为模式识别与智能机 器人。 王斐，男，１９７４ 年生，副教授，博士， 主要研究方向为模式识别与智能机器 人。 完成包括国家自然科学基金、辽宁 省自然科学基金、教育部博士点基金和 中国人民解放军总装备部预研项目在 内的国家和省部级科研项目 １５ 项，获 得辽宁省自然科学学术成果一等奖 １ 项。 发表学术论文 ５０ 余篇。 ［责任编辑：郑可为］ 第 ４ 期 邵明旭，等：人体下肢生物力学建模研究进展 ·５２７·