机器人杂志 1687 - 9619 1687 - 9600 Hindawi 10.1155 / 2020/3481056 3481056 研究文章 仿生假肢的设计和实验建模研究气动人工肌肉 https://orcid.org/0000 - 0002 - 0247 - 9516 华龙 1 https://orcid.org/0000 - 0002 - 2932 - 1823 Zhijie 1 2 命运 l 1 先进制造与自动化学院 机械工程与自动化学院 东北大学 沈阳 中国 neu.edu.cn 2 信息科学与工程学院 沈阳理工大学沈阳 沈阳 中国 2020年 28 2 2020年 2020年 19 11 2019年 08年 01 2020年 28 2 2020年 2020年 版权©2020华龙谢et al。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

在智能假肢的研究和开发,一些性能测试的实验是必需的。为了提供一个理想的智能假肢的性能测试实验平台,提出了一种异构双足步行机器人模型。假肢是异构双足步行机器人的一个重要组成部分,和它的主要功能是模拟残疾人健康正常的步态,提供智能仿生腿步态跟随目标轨迹。气动人工肌肉(PAM)有良好应用的假肢。假肢的仿生设计主要包括结构的髋关节,膝关节,踝关节,采用四杆机构作为膝关节的机械结构,和PAM作为驱动源的膝关节。其次,PAM性能测试平台的构建建立输出的力之间的关系,收缩率,输入压力测量等压条件下,建立了数学模型,PAM。最后,虚拟样机技术是用来建立一个联合仿真平台,和PID控制算法用于验证仿真。结果表明,人工腿可以按照目标轨迹。

中国国家自然科学基金 61803272 中央大学基础研究基金 中华人民共和国教育部的 N170313025
1。介绍

为了为下肢截肢者提供便利,使用智能假肢代替残疾人肢体生物领域的康复,和智能假肢的开发和研究需要理想的人工腿( 1]。大多数二足机器人使用直流(DC)伺服电机作为驱动装置,也和一些使用气动和液压设备。尽管液压装置驱动的位置精度可以保证,其刚度很大,可能发生泄漏和污染。气动设备驱动主要用于关节位置控制,和密封筒的要求很高。作为一种新型的气动元件,PAM具有结构简单、灵活性好,输出力大,重量轻,和类似的机械性能生物肌肉( 2]。因此,采用气动人工肌肉作为膝关节的传动装置,结合仿人机器人和智能假肢的研究研究假肢在异构双腿行走机器人模式 3),使其模拟健康腿的残疾人的步态。

周和Hannaford详细研究了PAM从生物学的角度来看 4, 5]。根据能量守恒的原理,气动人工肌肉的静态数学模型推断,和它的静态特性进行了测试。范Damme等人提出了一个Preisach-based模型来描述磁滞褶PAM ( 6]。钟发现Bouc-Wen模型的参数进行了优化,并与压力的关系,讨论了PAM。王等人提出了一种串级控制策略基于鲁棒建模简单、高效的轨迹跟踪控制的仿生关节由一个PAM ( 7, 8]。

La Rosa等人形成了一个链被禁止的化学突触的修改版本Hindmarsh-Rose slow-regular振荡模型,讨论了它的转变和混沌振荡 9]。基于神经网络(NN)的建模是最常见的类型。它有许多优点,如大规模并行处理、分布式存储的知识,自学能力强。由于神经网络技术不需要任何先验知识的过程的一个模型,估计是广泛应用于工业过程变量在实时监测和控制 10]。江等人提出了一个人形下肢肌肉的并联机构的形式被PAM(分布不均 11]。刘等人提出了一种仿生肩关节机器人一个球铰链是由一群PAM ( 12]。

要比人类假肢具有良好的灵活性和步态,开发一种改进假肢。总之,本文设计了机械系统的人工腿,建立PAM的实验平台,进行数学建模。此外,仿真验证虚拟样机的假肢,并通过实验研究得出结论。

2。结构设计和模型建立 2.1。机械系统设计的人工腿

针对假肢的总体结构设计,它使用PAM作为膝关节的传动装置和直流伺服电机的传动装置髋关节和踝关节。简化人工腿结构模型包括髋关节、膝关节、踝关节,每一个都有一个自由度。四连杆膝关节机制有很多优势( 13]。PAM驱动可以实现弯曲和扩展。图 1是一个膝关节的三维模型。髋关节控制旋角通过谐波减速器大腿的酒吧。踝关节是由谐波减速器和常见的直齿圆柱齿轮实现弯曲和脚趾弯曲。

四连杆膝关节的结构。

驾驶模式是一种链轮。如图 2(一个),最初的和对手的肌肉通过链轮链条连接。原来的肌肉产生动力,对手肌肉产生敌对的力量,和链式驱动链轮转动,所以联合可以完成拉伸和弯曲的动作。这种驱动方式的输出转矩较大,和膝关节可以满足负荷下运动的范围。

驾驶原则:(a)联合的方式由PAM和(b) PAM的安装位置。

考虑到仿人机器人的结构尺寸和所需的各关节的时刻,PAM的模型由费斯托是费斯托- mas - 20 - 200 - n - aa - mc。安装位置如图 2 (b)。PAM是倾斜的,和链轮的包角超过120°。虚拟样机的假肢的整体结构如图 3

整个结构的人工腿。

4假肢的控制系统图。运动控制卡与电脑通过ISA接口和一个直流伺服电机由PAMC控制电机驱动放大器,以髋关节和踝关节的运动控制。气动人工肌肉的运动控制卡控制通胀压力通过电动控制比例阀的压力信号,以控制膝关节的运动。在假肢的运动,编码器和电位计的反馈信息是通过D / A与工业控制计算机。

假肢的控制系统。

2.2。PAM和性能测试实验平台的建设

PAM的静态特性测试方案如图 5。主要由步进电机、压力比例阀、位移传感器、压力传感器、压力传感器、数据采集仪等。PAM的长度精确控制通过控制步进电机的转动和使用滚珠丝杆。压力比例阀控制气动人工肌肉的充气压力。位移传感器、压力传感器和压力传感器是用来测量收缩内部压力和输出气动人工肌肉的力量在实验中,分别。显示输出信号的数据采集仪器。图 6显示实际的实验。

PAM的实验程序的静态特征。1、空气压缩机;2、压力比例阀;3、压力传感器;4、张力传感器;5,气动人工肌肉;6,位移传感器;7,步进电机;8,滚珠丝杆。

的物理图的实验。

PAM的数学模型建立了通过实验建模。等压特性实验研究之间的关系的输出力和收缩速度气动人工肌肉在不同压力下通过加载条件下的步进电机的输入压力气动人工肌肉保持不变。

在均衡特性实验中,气动人工肌肉水平放置。一端是与张力传感器来测量输出的力的气动人工肌肉在实验。另一端连接步进电机的应用气动人工肌肉的外力。在实验中,首先,压力比例阀的输出压力调整,气动人工肌肉的内部压力达到设定值和仍然是固定的。此时,气动人工肌肉最大收缩状态。其次,旋转步进电机发送一个脉冲信号,线性位移和角位移转换为滚珠螺杆,和气动人工肌肉萎缩了通胀的作用下逐渐扩展外部负载。然后,停止旋转步进电机的原始长度。然后,改变步进电机的旋转方向,气动人工肌肉逐渐返回到最大收缩状态,完成一个工作循环。张力传感器和位移传感器分别测量输出的力和长度变化的气动人工肌肉在实验中,分别记录和数据后,系统是稳定的。最后,通胀压力气动人工肌肉的改变,重复上述实验过程,和气动人工肌肉的均衡特性实验在不同的压力下可以完成,从而获得PAM在不同的收缩力和收缩压力。

从图 7,可以看出内部摩擦的气动肌肉不同充气压力下几乎没有变化。如果充气气动人工肌肉的压力被认为是生物肌肉的激活水平,之间的特性曲线气动人工肌肉的收缩力和收缩类似于生物肌肉。

收缩力和收缩的关系平等条件下的压力。

2.3。气动人工肌肉的数学建模

气动人工肌肉的静态模型确定的等压特性试验。考虑弹性变形的橡皮管和网间的摩擦力的影响和橡胶层气动人工肌肉的数学模型,建立了气动人工肌肉的数学模型通过使用公式( 1)和( 2)参数识别( 14]。 (1) F = k 0 π r 0 2 p 一个 1 k 1 ε 2 b + k 2 F t , (2) f t = 2 3 一个 3 π r 0 l 0 t E 1 ε 2 1 3 3 b 一个 1 ε 2 , 在哪里 一个 = 3 / 棕褐色 2 θ 0 , b = 1 / 2 θ 0 , r0气动人工肌肉是初始半径, p是相对压力气动人工肌肉内腔, ε气动人工肌肉的收缩速度, l0是最初的气动人工肌肉的长度, θ0初始角的气动人工肌肉, t是橡皮管的壁厚对气动人工肌肉,然后呢 E气动人工肌肉橡皮管的弹性模量。假设PAM始终是一个理想的气缸,但在实际应用中,气动人工肌肉的结束将膨胀后半球形。PAM的收缩速率越大,收缩部分效果越小,并最终实际最大收缩速率要小于理论值。为了解决这个问题,身体萎缩率转化为有效的气动人工肌肉的收缩速度参数 k1。之间的静摩擦橡胶和纤维编织层和不考虑理论公式,但它是修改参数 k2。考虑到上述参数将不可避免地产生一些错误的理论公式修正,这将影响输出气动人工肌肉的力量,由参数补偿的因素 k0

使用最小二乘法确定气动人工肌肉的参数。的参数的估计价值气动人工肌肉在收缩和扩展了表 1 2

参数的估计值与PAM萎缩。

相对压力 p (MPa) k0 k1 k2
0.5 1.01 1.547 120.2
0.4 0.8576 2.183 303.7
0.3 0.7599 3.1 343.7
0.2 0.7616 4.461 298.9
0.1 0.9846 7.331 179.2

参数的估计值与PAM伸展。

相对压力 p (MPa) k0 k1 k2
0.5 0.8466 1.988 300.2
0.4 0.7655 2.433 326.8
0.3 0.6977 3.241 316.3
0.2 0.6804 4.676 266.2
0.1 0.761 11.06 203.2

确定参数带入修改模型。数学模型的气动人工肌肉在收缩和扩展。理论模型的比较,如图 8,当 p = 0.5 MPa

数学模型和实验数据之间的比较与PAM的收缩和伸展 p = 0.5 MPa

修改后的气动人工肌肉的数学模型可以准确反映时的静态特征 p = 0.1 , 0.2 , 0.3 , 0.4 , 0.5 MPa 。在实际应用程序中,气动人工肌肉的许多结构参数难以精确测量。因此,建立一个简单的数学模型来代替复杂的修改模型将有助于对气动人工肌肉的控制。

等压实验的特点,考虑到的输出力和位移的非线性气动人工肌肉在恒压下,拟合函数的数学模型,提出了气动人工肌肉如下: (3) F = 一个 b ε + c 经验值 μ ε , 在哪里 一个, b, c拟合系数, µ是非线性衰减系数描述收缩, μ= 45。

拟合函数的输入压力 p 一个( p), B( p), C( p假)气动肌肉在收缩和伸展模式,如图 9。在此基础上,两种模式的参数识别结果分别进行评估,如表所示 3

的拟合曲线 一个, b, c在压力 p 与PAM萎缩(a)和(b)。

的拟合相关系数 一个, b, c在压力 p 与PAM收缩和伸展。

系数 ρ(减少) ρ(拉伸)
一个( p) 0.9997 0.9995
b( p) 0.9647 0.9937
c( p) 0.9758 0.9548

最后,获得的气动人工肌肉的数学模型如下:

收缩状态: (4) F p , ε = 2442年 p 197.8 5464年 p + 1735年 ε + 688.3 p 9.3 经验值 45 ε

扩展状态: (5) F p , ε = 2554年 p 344.5 8430年 p 43.2 ε + 899.2 p + 43.24 经验值 45 ε

气动人工肌肉的静态数学模型可以改写如下: (6) F p , ε = f 0 ε + f 1 ε p , 在哪里 (7) f 0 ε = 197.8 1735年 ε 9.3 经验值 45 ε , 施工状态, 344.5 + 43.2 ε + 43.24 经验值 45 ε , 扩展的状态, f 1 ε = 2442年 5464年 ε + 688.3 经验值 45 ε , 施工状态 , 2554年 8403年 ε + 899.2 经验值 45 ε , 扩展状态

如图 10,数学模型可以准确反映PAM的静态特征。它是相对简单和控制仿真奠定了良好的基础。

静态数学模型与实验数据的比较,当PAM收缩和伸展。

3所示。实验仿真和结果分析 3.1。动态分析的驱动关节

两个气动人工肌肉组合驱动关节的连接链和链轮的气动人工肌肉。原理图所示 11

联合由PAM的原则。

如图 11, F1 F2分别输出屈肌和伸肌的力量。 R是力的手臂关节由气动人工肌肉驱动的, T驱动力矩, J的惯性负载, C是关节的阻尼, θ是关节的旋转角度。建立驱动关节的动力学方程如下: (8) J θ ¨ + C θ ˙ = T , T = F 1 F 2 R , T = φ 0 ε 1 , ε 2 + φ 1 ε 1 , ε 2 Δ p , Δ p = J θ ¨ + C θ ˙ φ 0 ε 1 , ε 2 φ 1 ε 1 , ε 2

3.2。PAM通货膨胀的动态模型

假设没有热量交换之间的气动人工肌肉和外部在膨胀过程中,这个过程可以被视为绝热膨胀过程。根据能量守恒方程的常数气体容积绝热充电源变量,气动人工肌肉充电过程的能量方程可以得到 15]: (9) d p d t = κ RT 年代 1 V κ p V d V d t , 在哪里 1 气体的质量流量,单位是千克/秒; p 是人造肌肉的内部压力,单位是爸爸; R是一种气体常数, RN = 287 (米/公斤 K); κ 是等熵指数和 κ = 1.4 ;和 T 年代 是气体源温度。

气动人工肌肉在曝气的质量流量方程如下: (10) 1 = 一个 e p 年代 1 b RT 年代 ω , 在哪里 一个 e 是管道系统的有效面积, p 年代 气源压力, T 年代 是气体源温度, b 压力比是至关重要的。

最后,动态模型的关节驱动人工肌肉得到如下: (11) ω = 1 σ b 1 b 2 , σ = p p 年代 > b , 1 , σ = p p 年代 b

3.3。虚拟样机联合仿真平台的建设

Pro / E、ADAMS和MATLAB / Simulink是用来建立一个虚拟样机联合仿真平台人工腿,取代物理模型来解决和控制仿真。建立三维模型导入ADAMS进行动态分析。主要的模型导入到ADAMS部分,关节,司机,地面,等等。每个部分由一对运动对连接。假肢的虚拟样机模型导入ADAMS如图 12,每个关节的约束和驾驶条件。

ADAMS虚拟样机的假肢。

假肢的虚拟样机模型,在ADAMS导入MATLAB / Simulink仿真,仿真系统框图和相应的控制是建立基于PID控制算法。联合控制假肢进行虚拟样机的仿真,以及理想步态之间的误差的原因和实际步态进行了分析。其中,步态数据来自人类的膝盖运动轨迹披露APAS软件。采样间隔0.005秒,仿真时间是0.2秒,仿真时间间隔是0.005秒。

膝关节,PID控制流图所示 13。理想的步态数据之间的误差和编码器测量值作为反馈信号,和压力比例阀的输出压力是由PID控制器调节输出气动人工肌肉的力量。在实际应用中,一个两个气动人工肌肉的膨胀和收缩,另一个是收缩和伸展,所以他们的输出力是不同的,这将产生驱动力矩的传动轴气动人工肌肉,并使用这一时刻作为输入控制模块。

膝关节由PAM的过程。

控制仿真、人工腿的虚拟样机由气动人工肌肉驱动的生成到模块my_test亚当斯/控制,及其结构如图 14。Control_torque膝盖关节的驱动力矩,和angle_1 angle_2是膝关节的旋转角度。以它为控制对象,系统控制仿真框图是由连接它与其他控制在MATLAB和Simulink模块。

控制植物由亚当斯。

基于PID控制算法的框图气动人工muscle-driven假肢的控制系统虚拟样机是建立在MATLAB / Simulink仿真,如图 15。在仿真过程中,理想步态轨迹之间的误差和实际步态作为反馈信号,控制模块的输入是纠正通过调整PID控制器,它是充气气动人工肌肉的压力,和膝关节的驱动力矩。作为输入ADAMS虚拟样机模型的包装,膝关节的旋转角度,这是实际的输出控制仿真系统。

仿真软件模型的PID控制算法。

使用“试错”的方法,实际输出曲线与理想曲线可以安装很多次的试错。

跟踪曲线和膝关节的跟踪误差,如图 16。下面的稳态误差相对较小,和下面的目标更好的实现。

跟踪曲线和膝关节的跟踪误差。

错误的原因如下:

结构设计:结构设计的过程中人工腿和膝关节,尽管机制进行了优化,它不同于人类下肢膝关节的结构

控制算法:选择PID控制算法的主要参数有很大的影响控制系统的性能,和“试错”的方法必然会对结果产生影响

的非线性和滞后的影响气动人工肌肉的控制系统

4所示。结论

摘要四杆机构作为膝关节的结构设计,和假肢的结构综合设计。它类似于人体的下肢关节,它可以更好地满足仿人假肢的要求。驱动,气动人工肌肉作为膝关节的驱动源,从而避免大的刚度等问题造成的使用伺服电机或液压设备。一个测试的实验平台建立了气动人工肌肉的性能,表明,气动人工肌肉的特性曲线类似于生物肌肉。建立更精确的数学模型,与实验数据相比。结果表明,它能准确反映静态气动人工肌肉的特点,为控制仿真奠定了基础。气动人工肌肉驱动的动态分析联合进行,和驾驶联合建立的数学模型。联合仿真平台构建和人工腿的虚拟样机模拟基于PID控制算法。仿真结果表明,人工腿气动人工肌肉驱动的可以按照目标轨迹。

4.1。限制

人体下肢有六个自由度,即简化为三个。更好的设计的假肢结构是必需的。PID控制算法具有适应性差的非线性和时变系统,需要更好的控制算法仿真。

数据可用性

使用的数据来支持本研究的发现可以从相应的作者。

的利益冲突

作者宣称没有利益冲突有关的出版。

确认

这项工作得到了国家自然科学基金(批准号61803272)和基础研究基金为中央大学中国教育部(批准号N170313025)。作者衷心感谢支持。

H。 K。 F。 仿生假肢设计和控制仿真基于模糊PID控制算法 学报》第五届IEEE国际会议网络技术在自动化、控制和智能系统 2015年6月 沈阳,中国 IEEE 2097年 2102年 10.1109 / cyber.2015.7288272 2 - s2.0 - 84962290745 Q。 X。 年代。 项目和气动人工肌肉在仿生关节的前景 中国液压和气动 2007年 4 50 51 F。 H。 X。 的双足机器人步态规划和控制实现异构双腿 东北大学学报:自然科学 2007年 28 1081年 1084年 悬疑类 g·K。 Czerntecki j . M。 赫纳夫德 B。 McKibben人造肌肉:气动执行机构与生物力学的情报 学报1999年IEEE / ASME国际会议上先进的智能机电一体化 1999年9月 亚特兰大,乔治亚州,美国 IEEE 221年 226年 10.1109 / aim.1999.803170 悬疑类 g·K。 赫纳夫德 B。 McKibben人工肌肉致动器的疲劳特征 学报1998年IEEE / RSJ智能机器人和系统国际会议 1998年10月 加拿大维多利亚 IEEE 1776年 1781年 10.1109 / iros.1998.724854 范Damme M。 Beyl P。 Vanderborght B。 建模滞后褶气动人工肌肉 诉讼IEEE国际会议上的机器人,自动化和机电一体化 2008年9月 成都,中国 IEEE 471年 476年 10.1109 / ramech.2008.4681431 2 - s2.0 - 58049103182 J。 风扇 J。 Y。 J。 W。 一个非线性PID控制来提高机械手的控制性能由气动肌肉驱动执行机构 机械工程的发展 2015年 6 ID 172782年 10.1155 / 2014/172782 2 - s2.0 - 84901761955 Y。 Q。 x H。 轨迹跟踪控制的仿生关节气动人工肌肉驱动的基于健壮建模 机器人 2016年 38 249年 256年 La Rosa M。 拉比诺维奇 m . I。 韦尔塔 R。 Abarbanel h . d . I。 命运 l 通过混沌振荡转移缓慢的正规化的单向链Hindmarsh-Rose模型 物理信 2000年 266年 1 88年 93年 10.1016 / s0375 - 9601 (00) 00015 - 3 2 - s2.0 - 0034645878 命运 l Giannone P。 Graziani 年代。 Xibilia m·G。 基于多层神经网络的虚拟仪器在炼油厂提高产品质量监控 IEEE仪表和测量 2007年 56 1 95年 101年 10.1109 / tim.2006.887331 2 - s2.0 - 33847381990 F。 G。 Q。 分析和控制一个平行的下肢基于气动人工肌肉 机械工程的发展 2017年 9 1 14 10.1177 / 1687814016685002 2 - s2.0 - 85013848492 Y。 T。 风扇 W。 静力学分析仿生肩关节机器人 北京理工的事务 2015年 35 607年 611年 X。 年代。 Q。 四假肢膝关节的优化和仿真 机械设计 2011年 28 42 45 P。 赫纳夫德 B。 Mckibben气动人工肌肉的测量和建模 IEEE机器人和自动 1996年 12 1 90年 102年 10.1109/70.481753 2 - s2.0 - 0030085935 J。 气力输送系统的动力学 1991年 广州,中国 华南理工大学出版社