CN114770519B - 一种关节力估计方法 - Google Patents
一种关节力估计方法 Download PDFInfo
- Publication number
- CN114770519B CN114770519B CN202210569906.XA CN202210569906A CN114770519B CN 114770519 B CN114770519 B CN 114770519B CN 202210569906 A CN202210569906 A CN 202210569906A CN 114770519 B CN114770519 B CN 114770519B
- Authority
- CN
- China
- Prior art keywords
- joint
- robot
- auxiliary
- standing
- track
- 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
Images
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
-
- 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/1612—Programme controls characterised by the hand, wrist, grip control
-
- 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/1664—Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
Landscapes
- Engineering & Computer Science (AREA)
- Robotics (AREA)
- Mechanical Engineering (AREA)
- Automation & Control Theory (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Manipulator (AREA)
Abstract
本发明属于关节力估计技术领域,尤其涉及一种关节力估计方法。其保证了老人起立过程的平稳性同时使辅助机器人提供适当的接触力,使辅助机器人的功能性和可靠性得到有效的提升。包括:步骤1、采集人体在起立动作时的姿态信息;并通过运动过程中重心的移动轨迹判断稳定性;步骤2、在确定起立过程处于稳定状态后,将人体视为左右对称的简化模型,并根据机器人正向运动学计算该过程中肩关节的移动轨迹,将该移动轨迹作为辅助机器人的参考轨迹;步骤3、在辅助机器人与用户接触时,基于参考轨迹估计辅助机器人肩关节力矩,控制机器人按照该轨迹完成辅助任务。
Description
技术领域
本发明属于关节力估计技术领域,尤其涉及一种关节力估计方法。
背景技术
人口老龄化是现代社会的一个重要人口特征,随着年龄的增长老人难以独自完成由坐到站的起立运动,越来越多的老人需要辅助机器人辅助起立。由于下肢肌肉无力,人们在使用进行机器人辅助起立时容易失去平衡,同时机器人与用户间的交互力过大或过小可能会造成二次伤害。
发明内容
本发明就是针对现有技术存在的缺陷,提供一种关节力估计方法。其基于服务机器人,针对现有辅助起立方法的不足,提出一种辅助机器人的起立辅助关节力估计方法,目的是保证老人起立过程的平稳性同时使辅助机器人提供适当的接触力,使辅助机器人的功能性和可靠性得到有效的提升。
为实现上述目的,本发明采用如下技术方案,其特征在于,包括:
步骤1、采集人体在起立动作时的姿态信息;并通过运动过程中重心的移动轨迹判断稳定性;
步骤2、在确定起立过程处于稳定状态后,将人体视为左右对称的简化模型,并根据机器人正向运动学计算该过程中肩关节的移动轨迹,将该移动轨迹作为辅助机器人的参考轨迹;
步骤3、在辅助机器人与用户接触时,基于参考轨迹估计辅助机器人肩关节力矩,控制机器人按照该轨迹完成辅助任务。
进一步地,所述姿态信息包括踝关节角度记为θa、膝关节角度记为θk、髋关节角度记为θx。
进一步地,步骤1中,所述通过运动过程中重心的移动轨迹判断稳定性的稳定性判断方法为:根据运动中人体重心轨迹是否保持在脚掌的范围内;
人体重心的运动轨迹计算方法:
式中,i为身体各部位的编号,xi是部位i重心在x坐标系的位置,mi是身体部位i的质量;若x<xd即重心始终在脚掌范围内,则起立过程中人体保持稳定,该起立轨迹可以作为辅助机器人的参考轨迹。
进一步地,获得身体数据,包括身体各部位的长度、质量计算公式如下:
其中,H为身高,L1为小腿长度,L2为大腿长度,L3为躯干长度,C1为小腿重心到踝关节的距离,C2为大腿重心到膝关节的距离,C3为躯干重心到髋关节的距离,W1为小腿的质量,W2为大腿的质量,W3为躯干的质量。
进一步地,起立运动中肩关节的移动轨迹求解方法:首先,通过对稳定性判断的分析选择能帮助用户平稳起立的辅助动作,然后根据姿态传感器记录的该动作姿态信息,包括踝关节角度记为θa,膝关节角度记为θk,髋关节角度记为θx,以及人体简化模型求解辅助起立过程中,以脚踝为原点的坐标系下肩关节的移动轨迹;计算过程如下:
其中,x10为膝关节在基坐标系下的横坐标,y10为基坐标系下的纵坐标,x20为髋关节在基坐标系下的横坐标,y20为髋关节在基坐标系下的纵坐标,x30为肩关节在基坐标系下的横坐标,y30为基坐标系下的纵坐标。
更进一步地,步骤3包括:根据线性二次调节器对机器人的相关状态变量及其微分引入状态反馈,完成对服务机器人的肩关节夹角及接触点位置和接触力大小的控制。具体过程如下:
机器人与用户交互时的受力分析:
状态空间表达式:
U=[Fx Fz]T (12)
U=-KX (13)
其中,M为机器人重量,g为重力加速度,F为机器人对用户施加的力,R为用户对机器人施加的力,q为机器人肩关节角度,I为机器人转动惯量,L为人机接触点距肩关节的长度。
更进一步地,不使用力传感器而通过仿真实验训练关于接触力的神经网络模型;首先,基于现有的服务机器人,在SolidWorks中绘制好机器人的实体模型并导入SimscapeMultibody工具箱中获得服务机器人仿真模型,再结合空间接触力模块实现对接触力的仿真,建立服务机器人助立仿真模型,利用该模型完成服务机器人辅助人体起立的仿真实验;仿真过程中通过模拟机器人与人的交互情况,记录辅助机器人受到的接触力法向分量和切向分量以及肩关节角度作为训练数据去训练神经网络;
实验内容如下:首先令服务机器人机械臂以相同的动作完成预定的轨迹,然后通过在机械臂上施加不同的接触力来模拟机器人扶人起立的过程,并记录此时机器人肩关节角度和机器人施加的关节力矩;
最后以记录的服务机器人辅助起立时接触力与肩关节力矩作为数据训练神经网络模型并得到其近似函数:
τ=h(q,Fx,Fz) (14)
式中:τ为辅助机器人肩关节的力矩,q为肩关节角度,Fx,Fz分别为辅助机器人与用户切向力和法向力;
神经网络模型的结构:
y=b2+wL tansig(b1+wIx) (15)
通过仿真获取1000组数据,将70%的数据用于培训,15%用于验证,15%用于测试,将机器人肩关节角度,接触力的切向法向分量作为神经网络输入量,它们的一维时间序列数据形成3×1的矩阵,共计1000个这样的矩阵。每个矩阵经过隐藏层输出机器人肩关节的力矩。当达到最大迭代次数、μ超过最大值或梯度低于最小值时,训练停止;
神经网络模型根据输入的肩关节角度和接触力信息,估算出服务机器人肩关节的输出力矩,并传输给机器人控制末端执行器完成辅助运动。最后通过随机选取人机间的接触力,对训练出的神经网络预测的机器人肩关节力矩与实际的关节力矩进行对比,以验证神经网络模型的有效性。最终的结果显示该神经网络模型能够根据输入的数据输出我们需求的关节力矩数据,且精确度满足要求。
与现有技术相比本发明有益效果。
本发明辅助机器人可以准确按照设计的平稳轨迹帮助用户起立,提高了人机接触过程中的安全性。在辅助机器人帮助下起立时,辅助机器人能够根据当前的位置信息和关节角度提供适当的辅助力。
附图说明
下面结合附图和具体实施方式对本发明做进一步说明。本发明保护范围不仅局限于以下内容的表述。
图1为人体起立过程示意图。
图2为本发明传感器安装位置示意图。
图3为本发明机器人辅助起立示意图。
图4为本发明人体简化模型示意图。
图5为本发明神经网络结构示意图。
图6为本发明参考轨迹求解示意图。
图7为本发明机器人关节力矩求解示意图。
图8为本发明角度信息示意图。
图9为本发明某实验者重心移动轨迹示意图。
图10为本发明人体起立肩关节移动轨迹示意图。
图11为本发明辅助机器人末端执行器轨迹跟踪结果示意图。
具体实施方式
如图1-11所示,本发明包括:
步骤1、采集人体在起立动作时的姿态信息;并通过运动过程中重心的移动轨迹判断稳定性。
步骤2、在确定起立过程处于稳定状态后,将人体视为左右对称的简化模型,并根据机器人正向运动学计算该过程中肩关节的移动轨迹,将该移动轨迹作为辅助机器人的参考轨迹。
步骤3、在辅助机器人与用户接触时,基于参考轨迹估计辅助机器人肩关节力矩,控制机器人按照该轨迹完成辅助任务。
优选地,所述姿态信息包括踝关节角度记为θa、膝关节角度记为θk、髋关节角度记为θx。
优选地,步骤1中,所述通过运动过程中重心的移动轨迹判断稳定性的稳定性判断方法为:根据运动中人体重心轨迹是否保持在脚掌的范围内;
人体重心的运动轨迹计算方法:
式中,i为身体各部位的编号,xi是部位i重心在基坐标系x轴的位置,mi是身体部位i的质量;若x<xd即重心始终在脚掌范围内,则起立过程中人体保持稳定,该起立轨迹可以作为辅助机器人的参考轨迹。
优选地,获得身体数据,包括身体各部位的长度、质量计算公式如下:
其中,H为身高,L1为小腿长度,L2为大腿长度,L3为躯干长度,C1为小腿重心到踝关节的距离,C2为大腿重心到膝关节的距离,C3为躯干重心到髋关节的距离,W1为小腿的质量,W2为大腿的质量,W3为躯干的质量。
优选地,起立运动中肩关节的移动轨迹求解方法:通过姿态传感器记录起立的姿态信息,包括踝关节角度记为θa,膝关节角度记为θk;髋关节角度记为θx,根据关节角度信息以及人体的简化模型计算在以脚踝为原点的坐标系中肩关节的移动轨迹。
其中,x10为膝关节在基坐标系下的横坐标,y10为基坐标系下的纵坐标,x20为髋关节在基坐标系下的横坐标,y20为髋关节在基坐标系下的纵坐标,x30为膝关节在基坐标系下的横坐标,y30为基坐标系下的纵坐标,θa、θk、θx分别踝关节、膝关节、肩关节的角度。
更优选地,步骤3包括:根据线性二次调节器对机器人的相关状态变量及其微分引入状态反馈,得到状态反馈控制矩阵K,基于机器人和用户接触位置控制接触力,具体过程如下:
机器人与用户交互时的受力分析:
其中,M为机器人重量,g为重力加速度,F为机器人对用户施加的力,R为用户对机器人施加的力,q为机器人肩关节角度,I为机器人转动惯量,L为人机接触点距肩关节的长度。
状态空间表达式:
U=[Fx Fz]T (12)
U=-KX (13)
更优选地,所述起立运动辅助机器人肩关节力矩估算方法包括:
基于现有的服务机器人,在SolidWorks中绘制好机器人的实体模型,接下来利用该模型进行服务机器人辅助人体起立的仿真实验。实验内容如下:首先令服务机器人机械臂以相同的动作完成预定的轨迹,然后通过在机械臂上施加不同的接触力来模拟机器人扶人起立的过程,并记录此时机器人肩关节角度和施加的力矩。
最后以记录的服务机器人辅助起立时接触力与肩关节力矩作为数据训练神经网络模型并得到其近似函数:
τ=h(q,Fx,Fz) (14)
式中:τ为辅助机器人肩关节的力矩,q为肩关节角度,Fx,Fz分别为辅助机器人与用户切向力和法向力;
神经网络模型的结构:
y=b2+wL tansig(b1+wIx) (15)
通过仿真获取1000组数据,将70%的数据用于培训,15%用于验证,15%用于测试,当达到最大迭代次数、μ超过最大值或梯度低于最小值时,训练停止;
神经网络模型根据输入的肩关节角度和接触力信息,估算出辅助机器人肩关节力矩,并传输给机器人控制末端执行器完成辅助运动。
下面结合附图对本发明做进一步详细说明:
本发明使用了服务机器人作为实验平台、姿态传感器用于采集数据。姿态传感器用弹性绷带将姿态传感器固定在小腿、大腿和躯干上,分别采集踝关节、膝关节和肩关节的角度。根据多次实验可知,人在起立过程中脚的位置保持不变,因此,以脚后跟作为坐标原点建立坐标系。
设计起立实验,为了能够得到用户身体保持稳定不失衡的轨迹,我们招募了10名无骨骼、肌肉、神经病史的健康受试者。受试者的平均年龄为24岁,平均身高为170cm,平均体重为60kg,所有受试者实验环境相同,实验在室内进行,每组实验耗时5s。每名受试者都被要求在同样搀扶人员帮助下完成10组起立实验,且每天只进行一组实验。
接下来根据运动中重心的移动轨迹判断用户起立过程的稳定性。首先根据实验者的身高H和体重W,计算出人体各部位的近似长度和质量。然后基于姿态传感器采集的角度信息结合身体各部位数据计算重心的移动轨迹。
身体各部位的长度、质量计算方法如下:
其中,H为身高,L1为小腿长度,L2为大腿长度,L3为躯干长度,C1为小腿重心到踝关节的距离,C2为大腿重心到膝关节的距离,C3为躯干重心到髋关节的距离,W1为小腿的质量,W2为大腿的质量,W3为躯干的质量。
通过计算出的重心轨迹在起立过程中是否保持在脚掌范围内判断起立动作的稳定性。若稳定,则根据关节角度信息以及人体简化模型计算肩关节的移动轨迹。
基于参考轨迹,在辅助机器人与用户接触时,估计辅助机器人肩关节力矩,进而控制辅助机器人按照期望轨迹完成辅助起立任务。首先利用线性二次调节器对机器人的相关状态变量及其微分引入状态反馈,得到状态反馈控制矩阵K,基于机器人和用户接触位置控制接触力,具体过程如下:
机器人与用户交互时的受力分析:
状态空间表达式:
U=[Fx Fz]T
U=-KX
然后根据接触力通过神经网络模型估算机器人肩关节的力矩。为了训练神经网络,使用Simscape Multibody建立机器人辅助用户起立时的简化模型,并根据简化模型的仿真结果训练神经网络模型,获得h的近似函数:
τ=h(q,Fx,Fz)
式中:τ为辅助机器人肩关节的力矩,q为肩关节角度,Fx,Fz分别为辅助机器人与用户切向力和法向力。
神经网络模型的结构:
y=b2+wL tansig(b1+wIx)
通过仿真获取1000组数据,将70%的数据用于培训,15%用于验证,15%用于测试,当达到最大迭代次数、μ超过最大值或梯度低于最小值时,训练停止。
神经网络模型根据输入的肩关节角度和接触力信息,估算出辅助机器人肩关节力矩,并传输给机器人控制末端执行器完成辅助运动。
为了验证方法的有效性,我们根据起立实验获得的参考轨迹进行了仿真。受试者A和B的其中一组实验数据如图7所示,其中受试者B重心的移动轨迹如图8所示,在1.5s时人体离开椅面,此时重心在0-17cm的范围内波动,没有超出脚掌的范围,则此起立动作是稳定的。因此基于该过程计算出肩关节的移动轨迹如图9所示。
考虑辅助机器人和用户的相对位置,将肩关节的轨迹换算成机器人末端执行器的参考轨迹,并通过仿真,对参考轨迹进行跟踪,结果如图10所示。由此可以看出该方法能够基于给定的轨迹为机器人提供适当的关节力矩,完成辅助起立任务。
可以理解的是,以上关于本发明的具体描述,仅用于说明本发明而并非受限于本发明实施例所描述的技术方案,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换,以达到相同的技术效果;只要满足使用需要,都在本发明的保护范围之内。
Claims (4)
1.一种关节力估计方法,其特征在于:包括:
步骤1、采集人体在起立动作时的姿态信息;并通过运动过程中人体重心的移动轨迹判断稳定性;
步骤2、在确定起立过程处于稳定状态后,将人体视为左右对称的简化模型,并根据辅助机器人正向运动学计算该过程中肩关节的移动轨迹,将该移动轨迹作为辅助机器人的辅助起立轨迹;
步骤3、基于辅助起立轨迹,将辅助机器人关节角度数据和接触力切向力分量、法向力分量作为输入,辅助机器人肩关节输出力矩作为输出训练神经网络;使辅助机器人沿着辅助起立轨迹运行并提供适当的关节力矩完成辅助起立任务;
基于人体统计学获得身体数据,包括身体各部位的长度、质量及重心位置,计算公式如下:
其中,H为身高,L1为小腿长度,L2为大腿长度,L3为躯干长度,C1为小腿重心到踝关节的距离,C2为大腿重心到膝关节的距离,C3为躯干重心到髋关节的距离,W1为小腿的质量,W2为大腿的质量,W3为躯干的质量;
步骤3包括:对辅助机器人辅助起立过程受力分析建立状态空间方程,根据线性二次调节器对辅助机器人的相关状态变量及其微分引入状态反馈,求解状态反馈控制矩阵K完成对辅助机器人与用户接触力的控制;具体过程如下:
辅助机器人与用户交互过程进行受力分析:
状态空间表达式:
U=[Fx Fz]T
U=-KX
其中,M为辅助机器人重量,g为重力加速度,F为辅助机器人对用户施加的力,R为用户对辅助机器人施加的力,q为辅助机器人肩关节角度,I为辅助机器人转动惯量,l为人机接触点距肩关节的长度。
2.根据权利要求1所述的一种关节力估计方法,其特征在于:所述姿态信息包括踝关节角度记为θa、膝关节角度记为θk、髋关节角度记为θx。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210569906.XA CN114770519B (zh) | 2022-05-24 | 2022-05-24 | 一种关节力估计方法 |
ZA2023/00417A ZA202300417B (en) | 2022-05-24 | 2023-01-10 | Method for estimating joint forces |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210569906.XA CN114770519B (zh) | 2022-05-24 | 2022-05-24 | 一种关节力估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114770519A CN114770519A (zh) | 2022-07-22 |
CN114770519B true CN114770519B (zh) | 2023-03-31 |
Family
ID=82408659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210569906.XA Active CN114770519B (zh) | 2022-05-24 | 2022-05-24 | 一种关节力估计方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN114770519B (zh) |
ZA (1) | ZA202300417B (zh) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101422907B (zh) * | 2008-12-16 | 2011-01-19 | 吉林大学 | 欠驱动双足步行机器人髋关节机构 |
CN102922508B (zh) * | 2012-09-21 | 2015-01-07 | 西安交通大学 | 一种电动汽车电池换装外骨骼机器人系统 |
CN109124916B (zh) * | 2018-06-30 | 2024-01-30 | 源珈力医疗器材国际贸易(上海)有限公司 | 一种辅助站立椅及其运动轨迹研究方法 |
CN111515929A (zh) * | 2020-04-15 | 2020-08-11 | 深圳航天科技创新研究院 | 人体运动状态估计方法、装置、终端及计算机可读存储介质 |
WO2021213639A1 (en) * | 2020-04-22 | 2021-10-28 | Abb Schweiz Ag | Method of controlling industrial robot, control system and robot system |
CN112494284A (zh) * | 2020-12-18 | 2021-03-16 | 深圳市迈步机器人科技有限公司 | 一种髋关节外骨骼的控制方法及装置 |
CN113081582B (zh) * | 2021-03-18 | 2022-06-28 | 上海交通大学 | 一种机器人辅助站立轨迹生成方法 |
-
2022
- 2022-05-24 CN CN202210569906.XA patent/CN114770519B/zh active Active
-
2023
- 2023-01-10 ZA ZA2023/00417A patent/ZA202300417B/en unknown
Also Published As
Publication number | Publication date |
---|---|
ZA202300417B (en) | 2023-04-26 |
CN114770519A (zh) | 2022-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ferrati et al. | Virtual modelling of a real exoskeleton constrained to a human musculoskeletal model | |
Toussaint et al. | Controlling the ground reaction force during lifting | |
Vashista et al. | Force adaptation in human walking with symmetrically applied downward forces on the pelvis | |
Mankala et al. | Novel swing-assist un-motorized exoskeletons for gait training | |
CN114770519B (zh) | 一种关节力估计方法 | |
Dinovitzer et al. | Accurate real-time joint torque estimation for dynamic prediction of human locomotion | |
Demir et al. | Human impedance parameter estimation using artificial neural network for modelling physiotherapist motion | |
Dainis | Cinematographic analysis of the handspring vault | |
Konosu et al. | Upper limb joint torques during performances of kicking pullovers | |
Bejczy | Towards development of robotic aid for rehabilitation of locomotion-impaired subjects | |
Rejman et al. | An evaluation of kinesthetic differentiation ability in monofin swimmers | |
Lee et al. | Estimation of unmeasured golf swing of arm based on the swing dynamics | |
Felton et al. | Are planar simulation models affected by the assumption of coincident joint centers at the hip and shoulder? | |
González et al. | Center of mass estimation for rehabilitation in a multi-contact environment: A simulation study | |
Kamnik et al. | Human voluntary activity integration in the control of a standing-up rehabilitation robot: A simulation study | |
Afzal et al. | Simulation of a patient driven strategy for FES supported sit-to-stand movement | |
Han et al. | Estimation of the center of bodymass during forward stepping using body acceleration | |
Hirata et al. | Control of wearable walking helper on slope based on integration of acceleration and GRF information | |
Hnat et al. | Estimating center of mass kinematics during perturbed human standing using accelerometers | |
Miyawaki et al. | Evaluation of the gait of elderly people using an assisting cart: gait on flat surface | |
Zhang et al. | Gait Prediction and Assist Control of Lower Limb Exoskeleton Based on Inertia Measurement Unit | |
Maddalena et al. | An optimized design of a parallel robot for gait training | |
Chambers et al. | A Model-Based Analysis of The Effect of Repeated Unilateral Low Stiffness Perturbations on Human Gait: Toward Robot-Assisted Rehabilitation | |
Chen et al. | Towards Human-like Walking with Biomechanical and Neuromuscular Control Features: Personalized Attachment Point Optimization Method of Cable-Driven Exoskeleton | |
Kondo et al. | Algorithm of pattern generation for mimicking disabled person’s gait |
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 |