CN116690567B - 一种基于拟人神经模型的下肢外骨骼驱动方法 - Google Patents
一种基于拟人神经模型的下肢外骨骼驱动方法 Download PDFInfo
- Publication number
- CN116690567B CN116690567B CN202310742393.2A CN202310742393A CN116690567B CN 116690567 B CN116690567 B CN 116690567B CN 202310742393 A CN202310742393 A CN 202310742393A CN 116690567 B CN116690567 B CN 116690567B
- Authority
- CN
- China
- Prior art keywords
- muscle
- model
- anthropomorphic
- gait
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 210000003141 lower extremity Anatomy 0.000 title claims abstract description 30
- 210000005036 nerve Anatomy 0.000 title claims abstract description 28
- 210000003205 muscle Anatomy 0.000 claims abstract description 84
- 230000005021 gait Effects 0.000 claims abstract description 65
- 238000005457 optimization Methods 0.000 claims abstract description 24
- 230000008602 contraction Effects 0.000 claims abstract description 17
- 230000004118 muscle contraction Effects 0.000 claims abstract description 17
- 230000015654 memory Effects 0.000 claims abstract description 10
- 230000000694 effects Effects 0.000 claims abstract description 6
- 238000003062 neural network model Methods 0.000 claims abstract description 6
- 210000003423 ankle Anatomy 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000033001 locomotion Effects 0.000 claims description 22
- 210000000544 articulatio talocruralis Anatomy 0.000 claims description 18
- 210000002027 skeletal muscle Anatomy 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 12
- 210000002414 leg Anatomy 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 9
- 210000002683 foot Anatomy 0.000 claims description 8
- 208000010428 Muscle Weakness Diseases 0.000 claims description 6
- 206010028372 Muscular weakness Diseases 0.000 claims description 6
- 230000007787 long-term memory Effects 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 6
- 210000002303 tibia Anatomy 0.000 claims description 6
- 239000000523 sample Substances 0.000 claims description 4
- 230000006403 short-term memory Effects 0.000 claims description 4
- 210000001519 tissue Anatomy 0.000 claims description 4
- 238000002604 ultrasonography Methods 0.000 claims description 4
- 206010062575 Muscle contracture Diseases 0.000 claims description 3
- 208000006111 contracture Diseases 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 210000002435 tendon Anatomy 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 2
- 208000006011 Stroke Diseases 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 4
- 210000003414 extremity Anatomy 0.000 description 3
- 208000020431 spinal cord injury Diseases 0.000 description 3
- 206010008190 Cerebrovascular accident Diseases 0.000 description 2
- 208000010886 Peripheral nerve injury Diseases 0.000 description 2
- 230000002490 cerebral effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000001537 neural effect Effects 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 210000001142 back Anatomy 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1602—Programme controls characterised by the control system, structure, architecture
- B25J9/161—Hardware, e.g. neural networks, fuzzy logic, interfaces, processor
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/0006—Exoskeletons, i.e. resembling a human figure
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1656—Programme controls characterised by programming, planning systems for manipulators
- B25J9/1661—Programme controls characterised by programming, planning systems for manipulators characterised by task planning, object-oriented languages
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Robotics (AREA)
- Mechanical Engineering (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Rehabilitation Tools (AREA)
Abstract
本发明提供了一种基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于,包括以下步骤:步骤1、根据用户个性化参数构造其对应拟人神经反馈模型;步骤2,通过传感器测量肌肉生物力学特性,解算肌肉收缩速度、肌肉力大小;步骤3,以肌肉收缩速度为输入,经由神经网络模型优化拟人神经反馈模型参数步速∝、步幅μ、步频ω;步骤4,计算步态识别结果,并通过长短期记忆权重优化方法,输出优化后步态识别结果;步骤5,根据步态事件得分识别对应事件,预测将到达的步态事件;步骤6,计算最大扭矩补偿值,电机输出助力;步骤7,每次助力完成后,将新解算的对应肌肉力对比正常人体步态事件肌肉力大小,输出助力效果。
Description
技术领域
本发明涉及康复器械技术领域,具体涉及个性化拟人神经模型,开发一种下肢外骨骼驱动方法。
背景技术
由脑卒中、脊髓损伤与周围神经损伤而导致的下肢肌无力患者往往需要康复器具进行其肢体运动的辅助。临床医学表明,脑卒中,脊髓损伤等疾病或生理功能的退化都可由步态特征信息所反映,准确地检测行走步态轨迹能更好地分析患者的步态能力与康复状况。
下肢外骨骼在助行方面不但需要协同患者下肢的运动规律,而且还需要准确的提供合适助力以修正患者步态。当前,下肢柔性外骨骼因其结构简单轻便、易于穿戴等优势已成为继刚性外骨骼后的潮流与方向,前期下肢外骨骼的控制策略系由正常人体运动作为标定进行患者步态动作的矫正,即按照正常人体肢体运动的角度与行程进行患者动作的矫正,该控制策略的人机协同性较差,忽视了患者个体的不同,外骨骼在助力步频、步幅、步速等方面完全按照预设指令实施,造成患者体验感下降。
目前的市场上,下肢柔性外骨应用广泛,以一种下肢助力外骨骼机器人的步态控制方法和控制系统(CN109172287B)专利为代表的下肢柔性外骨骼控制策略以对比人体模型实现了实时步态检测,能够识别出跖屈、背屈等动作进而实现助力,然而当实时采集步态进而驱动电机输出存在一定的延时,也造成了助力的滞后性,使助力效果下降。以(CN115294653A)为代表的大多数现有技术系基于惯性传感器(IMU)进行人体姿态的捕捉,IMU是间接测量的,需要复杂的身体参数进行后计算建模,以计算出运动,并不时进行额外的修正,以减少随时间增长的不对准误差;再以(CN110141239A)为代表的的现有技术进一步基于表面肌电的测量开展运动意图识别工作。此外,现有技术由步态开始初期开始进行反馈信号采集,而由于人体步态前几步的稳定性缺乏,而康复患者步态从开始达到稳定则需要的时间更长,故初期步态数据作为采集信号将对后续步态的识别,助力输出产生一定影响,目前此问题尚缺乏研究与探索。
综上所述,当前下肢柔性外骨骼的控制策略在人机协同性及个性化助力方面尚欠缺,因此,对于下肢柔性助力外骨骼而言,其核心问题是个性化人体运动意图识别、预测与助力驱动控制,建立基于用户个性化参数的驱动模型,进而实现外骨骼的个性化助力。
发明内容
本发明旨在至少解决现有技术或相关技术中存在的技术问题之一。
本发明的有益效果是本发明创造相较以往的技术,构建以人体个性化参数为指引的步态模型,精准识别用户步态,定义最佳优化驱动力,进行个性化步态精准识别预测,满足个性化助力需求。本发明将独特的个性化驱动模型与下肢柔性外骨骼结合,对具有不同体重、足码、步幅、步频、步速的下肢肌无力患者提供了个性化助力方案。
本发明的技术方案在于提供了一种基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于,包括以下步骤:
步骤1、即根据用户个性化参数体重、足码、步幅、步频、步速、助力位置,构造其对应拟人神经反馈模型,拟人神经反馈模型包含步态识别优化模型M1、事件模型M2、电机补偿力矩输出模型M3;
步骤2,通过传感器测量肌肉生物力学特性,解算肌肉收缩速度、肌肉力大小;
步骤21、通过贴附在胫骨前肌、腓肠肌的超声探头及采集卡收集超声反馈信号;
步骤22、基于超声反馈信号开展肌肉厚度解算;
步骤23、基于肌肉厚度T变化计算肌肉收缩速度;
步骤24、基于肌肉收缩速度计算肌肉力大小;
步骤25、基于肌肉力大小计算踝关节实际力矩;
步骤3,以肌肉收缩速度为输入,经由神经网络模型优化拟人神经反馈模型参数步速∝、步幅μ、步频ω;
步骤4,计算步态识别结果,并通过长短期记忆权重优化方法,输出优化后步态识别结果;即将优化的步速∝、步幅μ、步频ω代入至步态识别优化模型M1,经由踝关节角度计算模型M11与长短期记忆权重优化M12计算出步态识别参数
步骤5,根据步态事件得分识别对应事件,预测将到达的步态事件;
步骤6,计算最大扭矩补偿值,电机输出助力;即当电机收到输出信号时,结合输出模型M3计算最大扭矩补偿值,输出跖屈助力扭矩Tz,背屈助力扭矩Tb;
步骤7,每次助力完成后,将新解算的对应肌肉力对比正常人体步态事件肌肉力大小,输出助力效果。
进一步地,步态识别优化模型M1包含踝关节角度计算模型M11与长短期记忆权重优化M12,其中:
踝关节角度计算模型M11的数学模型为:(以左腿为例)
其中,θ为踝关节运动角度,参数σ、∝、μ、ω分别对应用户体重权重参数、收敛速度,幅值,以及周期,即患者个性化的步速、步幅、步频;
长短期记忆权重优化M12解算步骤为:
将踝关节运动角度划分为n个步态周期,以n个关节运动角度θ左为采集样本,按采集的时间先后顺序命名为为θ左1、θ左2、…、θ左n,通过长短期记忆权重优化步态识别参数并不断循环,其公式为:
进一步地,事件模型M2包含开始行走、跖屈、背屈、停止四个事件,其数学模型为:
其中,P表示踝关节动作事件,其根据不同的步态识别参数进行划分。
进一步地,电机补偿力矩输出模型输M3基于预测事件发生,结合人体开展最大助力扭矩不超过踝关节峰值扭矩的20%的助力输出,建立电机补偿力矩输出模型
其中,W为用户体重千克数,T实际MAX为基于测量模块计算的实际踝关节力矩峰值。
进一步地,步骤22中,
反馈信号的界面离超声发射器距离d计算公式为:
其中,t为脉冲信号从发射到接收的时间,c为超声在人体组织介质中的传播速度,为1400m/s;
故肌肉厚度T为BC界面之间的距离
T=dC-dB=1400(tC-tB)
其中,tC为界面C回声返回接收的时间,tB为界面B回声返回接收的时间。
进一步地,步骤23中,
基于连续采集的肌肉厚度T,以肌肉厚度的变化表征肌肉收缩速度,则肌肉收缩速度v为肌肉厚度的变化速率,对肌肉厚度求导得:
得到左右腿胫骨前肌、腓肠肌的收缩速度分别为:va左、va右、vs左、vs右。
进一步地,步骤24中,
根据Hill方程进行对应肌肉力的解算,则此时的肌肉力F大小为:
(v+b)(F+a)=b(F0+a)
其中,F0为肌肉挛缩时的张力,a为肌肉张力常量,b为肌肉收缩速度单位常量,由患者肌无力等级及个体详细情况定义;
得到左右腿胫骨前肌、腓肠肌的肌肉力分别为Fa左、Fa右、Fs右、Fs右。
进一步地,步骤25中,
背屈主要由胫骨前肌收缩完成,其踝关节力矩计算为
M实际=Fa左la
跖屈主要由腓肠肌收缩完成,其踝关节力矩计算为
M实际=Fs左ls
其中,la、ls分别为胫骨前肌肌腱及跟腱至踝关节旋转中心的垂直距离本发明的有益效果在于:
1、本发明建立了个性化拟人步态模型,能够有效识别对应步态;
2、本发明建立识别结果长短期记忆权重优化方法,能够精准识别步态;
3、本发明通过肌肉厚度变化表征了肌肉力大小变化及踝关节力矩变化;
4、本发明建立电机补偿力矩输出模型,设计峰值助力扭矩,加强装置安全性能。
附图说明
图1为发明方法流程示意图;
图2为拟人神经反馈模型示意图;
图3为肌肉收缩速度、肌肉力大小解算流程示意;
图4为超声传感器与下肢贴附示意;
图5为超声肌肉反馈信号示意;
图6踝关节实际力矩解算示意;
图7为神经网络模型优化人神经反馈模型参数步速∝、步幅μ、步频ω示意;
图8计算最大扭矩补偿值,电机输出助力示意图;
图9为下肢柔性外骨骼示意图;
图10为电机线轮示意图;
具体实施方式
为了可以更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本发明的实施方式及实施方式中的特征可以相互组合。
为了实现助力脑中风、脊髓损伤与周围神经损伤引起的肌无力等患者,本发明采用下肢柔性外骨骼进行患者在行走时跖屈、背屈的个性化助力。该技术方案将人体步态分为开始行走、跖屈、背屈、停止四个事件,基于建立用户个性化参数构造其对应拟人神经反馈模型及步态识别模型长短期记忆权重优化对患者个体步态的精准识别,再构建最大助力峰值计算模型计算电机助力输出结合传感采集数据进行骨肌动力学解算,并优化上述拟人神经反馈模型参数,开展预测人体运动意图,进而实时输出助力。
该实施例中,选定的患者体重65Kg,BMI指数为20.5,体重权重参数设置为0.25,屈助力提拉点至足踝中心点距离为l1=0.1m,背屈助力提拉点至足踝中心点距离为l2=0.15m,电机线轮的半径r=0.025m;以采集患者5个正常步态作为建立拟人神经反馈模型的基础开展助力。
基于该实施例,提供了一种基于个性化拟人神经模型的下肢外骨骼驱动方法,具体包括以下步骤:
步骤1,根据用户个性化参数构造其对应拟人神经反馈模型;
该步骤中,如图2所示,该实施例首先根据用户个性化参数,体重、足码、步幅、步频、步速、助力位置,构造其对应拟人神经反馈模型,其包含步态识别优化模型M1、事件模型M2、电机补偿力矩输出模型M3;
步态识别优化模型M1包含踝关节角度计算模型M11与长短期记忆权重优化M12,其中:
踝关节角度计算模型M11的数学模型为:(以左腿为例)
其中,θ为踝关节运动角度,参数σ、∝、μ、ω分别对应用户体重权重参数、收敛速度,幅值,以及周期,即患者个性化的步速、步幅、步频;
长短期记忆权重优化M12解算步骤为:
采集患者运动的5个步态,将踝关节运动角度划分为5个步态周期,以5个关节运动角度θ为采集样本,以左腿为例,按采集的时间先后顺序命名为θ左1、θ左2、…、θ左5,通过长短期记忆权重优化步态识别参数ε表示权重参数,赋值以最近采集的步态权重占比最高,逐渐降低,并不断循环,其公式为:
令ε5=0.5,ε4=0.25,ε3=0.125,ε2=0.07,ε1=0.55,
则
如图2,事件模型M2包含开始行走、跖屈、背屈、停止四个事件,其数学模型为:
其中,P表示踝关节动作事件,其根据不同的步态识别参数进行划分。
如图2、图8所示,电机补偿力矩输出模型输M3基于预测事件发生,结合人体开展最大助力扭矩不超过踝关节峰值扭矩的20%的助力输出,建立电机补偿力矩输出模型
其中,65为用户体重千克数,T实际MAX为基于测量模块计算的实际踝关节力矩峰值;
如图8、图9、图10所示,根据足长及提拉点位置计算电机在各事件中需输出的助力扭矩,跖屈助力时,提拉点至足踝中心点距离为l1=0.1m,电机线轮的半径r=0.025m,扭矩Tz计算为:
背屈助力时,提拉点至足踝中心点距离为l2=0.15m,电机线轮的半径r=0.025m扭矩Tb计算为:
步骤2,通过传感器测量肌肉生物力学特性,解算肌肉收缩速度、肌肉力大小;
如图3所示,该步骤通过以下方式实现:
步骤21、通过贴附在胫骨前肌、腓肠肌的超声探头及采集卡收集超声反馈信号;
如图4所示,将传感器贴附在患者胫骨前肌、腓肠肌的外侧,进行信号的发射与采集;
如图5所示,反馈回声的采集信号在骨肌组织的交界面发生明显的变化,根据反馈信号将不同深度的区域划分为皮肤界面A、脂肪-肌肉界面B、肌肉-骨骼界面C,界面BC之间的距离即为肌肉厚度T;
步骤22、基于超声反馈信号开展肌肉厚度解算;
如图5,一组反馈信号的界面离超声发射器距离d计算公式为:
其中,t为脉冲信号从发射到接收的时间,c为超声在人体组织介质中的传播速度,约为1400m/s;
故肌肉厚度T为BC界面之间的距离
T=dC-dB=1400(tC-tB)
其中,tC为界面C回声返回接收的时间,tB为界面B回声返回接收的时间;
步骤23、基于肌肉厚度T变化计算肌肉收缩速度;
基于连续采集的肌肉厚度T,以肌肉厚度的变化表征肌肉收缩速度,则肌肉收缩速度v为肌肉厚度的变化速率,对肌肉厚度求导得
上述方法得到左右腿胫骨前肌、腓肠肌的收缩速度分别为:va左、va右、vs左、vs右;
步骤24、基于肌肉收缩速度计算肌肉力大小;
根据Hill方程进行对应肌肉力的解算,则此时的肌肉力F大小为:
(v+b)(F+a)=b(F0+a)
其中,F0为肌肉挛缩时的张力,a为肌肉张力常量,b为肌肉收缩速度单位常量,由患者肌无力等级及个体详细情况定义;
上述方法得到左右腿胫骨前肌、腓肠肌的肌肉力分别为Fa左、Fa右、Fs右、Fs右;
步骤25、基于肌肉力大小计算踝关节实际力矩;
以左腿为例,背屈主要由胫骨前肌收缩完成,其踝关节力矩计算为
M实际=Fa左la
跖屈主要由腓肠肌收缩完成,其踝关节力矩计算为
M实际=Fs左ls
如图5,其中,la、ls分别为胫骨前肌肌腱及跟腱至踝关节旋转中心的垂直距离;
当患者足码为42码时,la约为0.055m,ls约为0.015m
步骤3,以肌肉收缩速度为输入,经由神经网络模型优化拟人神经反馈模型参数步速∝、步幅μ、步频ω
如图6,以步骤23所计算的左右腿胫骨前肌、腓肠肌的收缩速度为输入量,构建神经网络模型,设置合理的隐藏层及神经元数量,优化拟人神经反馈模型参数步速∝、步幅μ、步频ω;
为了使得系统计算量导致的延迟与优化效果达到较好的平衡,在此选择隐藏层为3层,各层神经元数量为4个,经测试得到系统计算延时为20ms;
步骤4,计算步态识别结果,并通过长短期记忆权重优化方法,输出优化后步态识别结果;
将优化的步速∝、步幅μ、步频ω代入至步态识别优化模型M1,经由踝关节角度计算模型M11与长短期记忆权重优化M12计算出步态识别参数
步骤5,根据步态事件得分至划分,预测将到达的步态事件;
事件模型M2将根据步态识别参数值的大小进行踝关节运动事件的判断,并将以预测将要发生事件的时机,提前系统延时时间20ms驱动电机进行助力输出;
步骤6,计算最大扭矩补偿值,电机输出助力;
如图5、图8所示,当电机收到输出信号时,结合输出模型M3计算最大扭矩补偿值,输出跖屈助力扭矩Tz,背屈助力扭矩Tb;
如图9所示,实施例还提供了下肢柔性外骨骼,肢柔性外骨骼装置由左电机、右电机、左线轮、右线轮、左跖屈鲍登线、左背屈鲍登线、右跖屈鲍登线、右背屈鲍登线、左鞋垫、右鞋垫组成;
左线轮连接在左电机上,顺时针旋转为提拉左跖屈鲍登线,逆时针旋转为提拉左背屈鲍登线;
左跖屈鲍登线与鞋垫鞋跟部连接;
左背屈鲍登线与鞋垫鞋弓部连接。
下肢柔性外骨骼装置由左电机、右电机、左线轮、右线轮、左跖屈鲍登线、左背屈鲍登线、右跖屈鲍登线、右背屈鲍登线、左鞋垫、右鞋垫组成;
左线轮连接在左电机上,顺时针旋转为提拉左跖屈鲍登线,逆时针旋转为提拉左背屈鲍登线;
左跖屈鲍登线与鞋垫鞋跟部连接;左背屈鲍登线与鞋垫鞋弓部连接。
如图9、图10所示,当外骨骼执行左腿跖屈助力时,左电机以扭矩Tz执行顺时针旋转,驱动鲍登线进行足根部的提拉助力;
如图9、图10所示,当外骨骼执行左腿背屈助力时,左电机以扭矩Tb执行逆时针旋转,驱动鲍登线进行足弓部的提拉助力;
步骤7,每次助力完成后,将新解算的对应肌肉肌肉力对比正常人体步态事件肌肉力大小,输出助力效果
如患者在行走初期,步态中对应胫骨前肌、腓肠肌肌肉力(F原)与正常人体步态下对应肌肉力(F正)差距较大,再一段时间的助力后,实测(F测)该差距减小,则助力百分比f为:
本发明的描述中,需要理解的是,术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或单元必须具有特定的方向、以特定的方位构造和操作,因此,不能理解为对本发明的限制。
在本说明书的描述中,术语“一个实施例/方式”、“一些实施例/方式”、“具体实施例/方式”等的描述意指结合该实施例/方式或示例描述的具体特征、结构、材料或特点包含于本发明的至少一个实施例/方式或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例/方式或实例。而且,描述的具体特征、结构、材料或特点可以在任何的一个或多个实施例/方式或示例中以合适的方式结合。
以上仅为本发明的优选实施例/方式而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于,包括以下步骤:
步骤1、即根据用户个性化参数体重、足码、步幅、步频、步速、助力位置,构造其对应拟人神经反馈模型,拟人神经反馈模型包含步态识别优化模型M1、事件模型M2、电机补偿力矩输出模型M3;
步骤2,通过传感器测量肌肉生物力学特性,解算肌肉收缩速度、肌肉力大小;
步骤21、通过贴附在胫骨前肌、腓肠肌的超声探头及采集卡收集超声反馈信号;
步骤22、基于超声反馈信号开展肌肉厚度解算;其中:根据反馈信号将不同深度的区域划分为皮肤界面A、脂肪-肌肉界面B、肌肉-骨骼界面C,界面BC之间的距离即为肌肉厚度T;
步骤23、基于肌肉厚度T变化计算肌肉收缩速度;
步骤24、基于肌肉收缩速度计算肌肉力大小;
步骤25、基于肌肉力大小计算踝关节实际力矩;
步骤3,以肌肉收缩速度为输入,经由神经网络模型优化拟人神经反馈模型参数步速、步幅/>、步频/>;
步骤4,计算步态识别结果,并通过长短期记忆权重优化方法,输出优化后步态识别结果;即将优化的步速、步幅/>、步频/>代入至步态识别优化模型M1,经由踝关节角度计算模型M11与长短期记忆权重优化M12计算出步态识别参数/>;其中:步态识别优化模型M1包含踝关节角度计算模型M11与长短期记忆权重优化M12:
踝关节角度计算模型M11的数学模型为:
其中,为踝关节运动角度,参数/>、/>、/>、/>分别对应用户体重权重参数、收敛速度,幅值,以及周期,即患者个性化的步速、步幅、步频;
长短期记忆权重优化M12解算步骤为:
采集患者运动的5个步态,将踝关节运动角度划分为5个步态周期,以5个关节运动角度为采集样本,按采集的时间先后顺序命名为为/>、/>、…/>,通过长短期记忆权重优化步态识别参数/>,并不断循环,其公式为:
令=0.5,/>=0.25,/>=0.125,/>=0.07,/>=0.55,
则
步骤5,根据步态事件得分识别对应事件,预测将到达的步态事件;
步骤6,计算最大扭矩补偿值,电机输出助力;即当电机收到输出信号时,结合输出模型M3计算最大扭矩补偿值,输出跖屈助力扭矩Tz,背屈助力扭矩Tb;
步骤7,每次助力完成后,将新解算的对应肌肉力对比正常人体步态事件肌肉力大小,输出助力效果。
2.根据权利要求1所述的基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于:事件模型M2包含开始行走、跖屈、背屈、停止四个事件,其数学模型为:
其中,P表示踝关节动作事件,其根据不同的步态识别参数进行划分。
3.根据权利要求1所述的基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于:
电机补偿力矩输出模型输M3基于预测事件发生,结合人体开展最大助力扭矩不超过踝关节峰值扭矩的20%的助力输出,建立电机补偿力矩输出模型:
其中,W为用户体重千克数,为基于测量模块计算的实际踝关节力矩峰值。
4.根据权利要求1所述的基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于:步骤22中,
反馈信号的界面离超声发射器距离d计算公式为:
其中,t为脉冲信号从发射到接收的时间,c为超声在人体组织介质中的传播速度,为1400m/s;
故肌肉厚度T为BC界面之间的距离
其中,为界面C回声返回接收的时间,/>为界面B回声返回接收的时间。
5.根据权利要求1所述的基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于:步骤23中,
基于连续采集的肌肉厚度T,以肌肉厚度的变化表征肌肉收缩速度,则肌肉收缩速度v为肌肉厚度的变化速率,对肌肉厚度求导得:
得到左右腿胫骨前肌、腓肠肌的收缩速度分别为:、/>、/>、/>。
6.根据权利要求1所述的基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于:步骤24中,
根据Hill方程进行对应肌肉力的解算,则此时的肌肉力F大小为:
其中,为肌肉挛缩时的张力,a为肌肉张力常量,b为肌肉收缩速度单位常量,由患者肌无力等级及个体详细情况定义;
得到左右腿胫骨前肌、腓肠肌的肌肉力分别为、/>、/>、/>。
7.根据权利要求1所述的基于个性化拟人神经模型的下肢外骨骼驱动方法,其特征在于:步骤25中,
背屈主要由胫骨前肌收缩完成,其踝关节力矩计算为
跖屈主要由腓肠肌收缩完成,其踝关节力矩计算为
其中,、/>分别为胫骨前肌肌腱及跟腱至踝关节旋转中心的垂直距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310742393.2A CN116690567B (zh) | 2023-06-21 | 2023-06-21 | 一种基于拟人神经模型的下肢外骨骼驱动方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310742393.2A CN116690567B (zh) | 2023-06-21 | 2023-06-21 | 一种基于拟人神经模型的下肢外骨骼驱动方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116690567A CN116690567A (zh) | 2023-09-05 |
CN116690567B true CN116690567B (zh) | 2024-03-12 |
Family
ID=87825478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310742393.2A Active CN116690567B (zh) | 2023-06-21 | 2023-06-21 | 一种基于拟人神经模型的下肢外骨骼驱动方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116690567B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103610568A (zh) * | 2013-12-16 | 2014-03-05 | 哈尔滨工业大学 | 一种拟人化的下肢助力外骨骼机器人 |
JP2014104549A (ja) * | 2012-11-28 | 2014-06-09 | Advanced Telecommunication Research Institute International | 外骨格型ロボットおよびリハビリテーション装置 |
CN106217352A (zh) * | 2016-08-17 | 2016-12-14 | 尖叫智能科技(上海)有限公司 | 基于压力传感器的外骨骼机器人步态控制方法 |
CN107397649A (zh) * | 2017-08-10 | 2017-11-28 | 燕山大学 | 一种基于径向基神经网络的上肢外骨骼康复机器人控制方法 |
CN107423729A (zh) * | 2017-09-20 | 2017-12-01 | 湖南师范大学 | 一种面向复杂视觉场景下的远距离类脑三维步态识别系统及实现方法 |
CN109940584A (zh) * | 2019-03-25 | 2019-06-28 | 杭州程天科技发展有限公司 | 一种外骨骼机器人及其检测人体运动意图的检测方法 |
CN110653817A (zh) * | 2019-08-20 | 2020-01-07 | 南京航空航天大学 | 一种基于神经网络的外骨骼机器人助力控制系统与方法 |
CN114043459A (zh) * | 2021-11-25 | 2022-02-15 | 湖南大学 | 柔性下肢外骨骼控制方法、外骨骼控制系统及使用方法 |
EP4011347A1 (en) * | 2020-12-14 | 2022-06-15 | Able Human Motion, S.L. | Powered-knee exoskeleton system |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9221177B2 (en) * | 2012-04-18 | 2015-12-29 | Massachusetts Institute Of Technology | Neuromuscular model-based sensing and control paradigm for a robotic leg |
-
2023
- 2023-06-21 CN CN202310742393.2A patent/CN116690567B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014104549A (ja) * | 2012-11-28 | 2014-06-09 | Advanced Telecommunication Research Institute International | 外骨格型ロボットおよびリハビリテーション装置 |
CN103610568A (zh) * | 2013-12-16 | 2014-03-05 | 哈尔滨工业大学 | 一种拟人化的下肢助力外骨骼机器人 |
CN106217352A (zh) * | 2016-08-17 | 2016-12-14 | 尖叫智能科技(上海)有限公司 | 基于压力传感器的外骨骼机器人步态控制方法 |
CN107397649A (zh) * | 2017-08-10 | 2017-11-28 | 燕山大学 | 一种基于径向基神经网络的上肢外骨骼康复机器人控制方法 |
CN107423729A (zh) * | 2017-09-20 | 2017-12-01 | 湖南师范大学 | 一种面向复杂视觉场景下的远距离类脑三维步态识别系统及实现方法 |
CN109940584A (zh) * | 2019-03-25 | 2019-06-28 | 杭州程天科技发展有限公司 | 一种外骨骼机器人及其检测人体运动意图的检测方法 |
CN110653817A (zh) * | 2019-08-20 | 2020-01-07 | 南京航空航天大学 | 一种基于神经网络的外骨骼机器人助力控制系统与方法 |
EP4011347A1 (en) * | 2020-12-14 | 2022-06-15 | Able Human Motion, S.L. | Powered-knee exoskeleton system |
CN114043459A (zh) * | 2021-11-25 | 2022-02-15 | 湖南大学 | 柔性下肢外骨骼控制方法、外骨骼控制系统及使用方法 |
Non-Patent Citations (2)
Title |
---|
下肢康复外骨骼机器人结构设计与运动学分析;曹慧林等;《机械设计与研究》;20200820;第36卷(第4期);第12-17页 * |
基于人体模型映射的步态识别及意图感知方法;贾晓辉等;《仪器仪表学报》;第41卷(第12期);第236-244页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116690567A (zh) | 2023-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101061984B (zh) | 利用肌电信号提供机械帮助的康复机器人系统 | |
Ishikawa et al. | Muscle-tendon interaction and elastic energy usage in human walking | |
Buchanan et al. | Estimation of muscle forces and joint moments using a forward-inverse dynamics model | |
Dugan et al. | Biomechanics and analysis of running gait | |
JP4541867B2 (ja) | 外力制御方法、外力制御システム及び外力制御プログラム | |
CN113101134B (zh) | 一种基于动力外骨骼的儿童下肢运动辅助康复系统 | |
CN105213153A (zh) | 基于脑肌信息变阻抗的下肢康复机器人控制方法 | |
JP4138814B2 (ja) | 運動管理システム | |
CN105288933A (zh) | 并联下肢康复机器人自适应训练控制方法及康复机器人 | |
WO2003002054A1 (fr) | Systeme de transmission de couple | |
US20140100494A1 (en) | Smart gait rehabilitation system for automated diagnosis and therapy of neurologic impairment | |
KR20080024695A (ko) | 보행훈련제어시스템 | |
CN112842824B (zh) | 一种用于下肢康复恢复的训练方法 | |
US20230330835A1 (en) | System and process for estimation of gait phase for use with a wearable robot | |
CN116690567B (zh) | 一种基于拟人神经模型的下肢外骨骼驱动方法 | |
JP2019084130A (ja) | 歩行動作評価装置、歩行動作評価方法及びプログラム | |
CN113995629A (zh) | 基于镜像力场的上肢双臂康复机器人导纳控制方法及系统 | |
Wang et al. | A real-time walking pattern recognition method for soft knee power assist wear | |
CN115544777B (zh) | 一种表征关节助力补偿值的方法及系统 | |
Kim et al. | Gait training algorithm of an end-effector typed hybrid walking rehabilitation robot | |
WO2023035457A1 (zh) | 助行器的控制方法、系统和驱动设备 | |
TWI687215B (zh) | 下肢外骨骼機器人及其輔助方法 | |
Safarov et al. | Robotic Device for Active-Passive Mechanotherapy of the Ankle Joint | |
KR102630930B1 (ko) | 웨어러블 로봇의 동력 보조 방법, 장치 및 컴퓨터프로그램 | |
JP2002165799A (ja) | 健康状態診断方法およびその方法を用いた健康状態診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |