CN106503364B - 一种不确定性条件下的下肢外骨骼时变可靠性分析方法 - Google Patents

一种不确定性条件下的下肢外骨骼时变可靠性分析方法 Download PDF

Info

Publication number
CN106503364B
CN106503364B CN201610953616.XA CN201610953616A CN106503364B CN 106503364 B CN106503364 B CN 106503364B CN 201610953616 A CN201610953616 A CN 201610953616A CN 106503364 B CN106503364 B CN 106503364B
Authority
CN
China
Prior art keywords
joint
lower limb
hip
limb exoskeleton
ankle
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
Application number
CN201610953616.XA
Other languages
English (en)
Other versions
CN106503364A (zh
Inventor
汪忠来
俞水
李文生
孟德彪
张小玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610953616.XA priority Critical patent/CN106503364B/zh
Publication of CN106503364A publication Critical patent/CN106503364A/zh
Application granted granted Critical
Publication of CN106503364B publication Critical patent/CN106503364B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Evolutionary Computation (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Pathology (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Prostheses (AREA)
  • Manipulator (AREA)

Abstract

本发明公开一种不确定性条件下的下肢外骨骼时变可靠性分析方法，通过简化下肢外骨骼模型，建立独立的髋关节、膝关节和踝关节三个关节的简化模型；并充分考虑下肢外骨骼机械结构存在的不确定性因素，建立在不确定性条件下关节角的数学模型；利用运动学正问题求解方法建立末端轨迹关于不确定条件下关节角角度的数学模型；通过运动精度分析得到在不确定条件下关节角角度和末端轨迹的均值与方差；考虑下肢外骨骼失效时序，将下肢外骨骼从髋关节到末端轨迹看作自上而下的串联系统，分析每个单元的失效概率，实现下肢外骨骼系统时变可靠性计算；结果对全面提高下肢外骨骼的设计水平具有较高的理论支撑和工程实践意义。

Description

一种不确定性条件下的下肢外骨骼时变可靠性分析方法
技术领域
本发明属于可靠性工程技术领域，具体是一种面向下肢外骨骼的时变可靠性分析方法。
背景技术
可穿戴下肢外骨骼机器人是一种能够使人在负重条件下进行长距离行走的人机一体化系统。由于它的可穿戴性，这就要求下肢外骨骼必须通过模拟人体下肢的骨骼结构、肌肉运动和关节转动给人体提供额外的动力，可以给穿戴者起到支撑、助力和保护作用。同时由于它的助力作用，这使得下肢外骨骼在医疗健康、军事作战和助老助残等方面广泛应用，目前下肢外骨骼的研究已成为了国际前沿研究方向。美国、日本等国家的一些科研机构的研究成果已经在单兵军事作战装备、辅助医疗设备、助老助残等领域获得了实际应用。我国在“十一五”国家科技支撑计划“残疾人康复技术及设备研发”重点项目中提出了运动功能康复方面的重点课题“肢体康复医疗机器人研发”，它标志着下肢外骨骼技术在我国备受重视，并会在不久的将来推向市场。
下肢外骨骼在设计、制造及其使用过程中将受到众多不确定性因素的影响，例如：制造精度误差、运动副间隙误差、液压缸驱动源误差和载荷、客户使用和工作环境的随机变化。这些不确定性将对下肢外骨骼与人体配合的一致性、人体穿戴的舒适性、人体穿戴的安全性等性能产生严重的影响。因此，对于具有杆件精度误差、运动副间隙误差以及液压缸驱动源误差等不确定性的下肢外骨骼，结合下肢外骨骼的尺寸、D-H转换矩阵和不确定性分析，需要一种不确定性条件下的下肢外骨骼时变可靠性分析方法，以对其进行可靠性分析。
发明内容
本发明为解决上述技术问题，提出了一种不确定性条件下的下肢外骨骼时变可靠性分析方法，通过对下肢外骨骼进行简化，考虑下肢外骨骼存在的不确定性因素，从而获取时变可靠性分析模型，为下肢外骨骼一致性、安全性分析奠定理论基础和技术支撑。
本发明采用的技术方案是：一种不确定性条件下的下肢外骨骼时变可靠性分析方法，包括：
S1、建立下肢外骨骼包含液压缸的简化模型，量化下肢外骨骼各关节杆件、运动副和液压缸的不确定性，计算在不确定性条件下各关节角角度以及末端轨迹的均值和方差；
S2、考虑下肢外骨骼失效时序，将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，分析每个单元的失效概率，并确定下肢外骨骼失效概率；所述的四单元依次为：髋关节、膝关节、踝关节以及末端轨迹；
S3、根据步骤S2所得的下肢外骨骼失效概率确定下肢外骨骼的可靠度。
进一步地，所述步骤S1包括：
S11、分析下肢外骨骼，建立下肢外骨骼髋关节、膝关节和踝关节对应的关节角角度在各自不确定性条件下的数学模型；并确定各自关节角角度的均值和方差；
S12、通过D-H转换矩阵，建立末端轨迹关于髋关节、膝关节和踝关节对应的关节角角度的数学模型，并求解下肢外骨骼末端轨迹的上、下限。
更进一步地，所述步骤S11具体包括以下分步骤：
S111、建立下肢外骨骼髋关节、膝关节和踝关节各自包含液压缸的简化模型，得到在满足外骨骼正常步态下的液压缸理想位移；
S112、分别量化髋关节、膝关节和踝关节的不确定性；
S113、根据步骤S111得到的液压缸理想位移，建立下肢外骨骼髋关节、膝关节和踝关节各自对应的关节角角度在步骤S112得到的不确定性条件下的数学模型，并确定各自关节角角度的均值和方差。
进一步地，步骤S112所述的髋关节不确定性至少包括：髋关节杆件尺寸误差、髋关节液压缸误差以及髋关节运动副间隙误差；
所述的膝关节不确定性至少包括：膝关节杆件尺寸误差、膝关节液压缸误差以及膝关节运动副间隙误差；
所述的踝关节不确定性至少包括：踝关节杆件尺寸误差、踝关节液压缸误差以及踝关节运动副间隙误差。
进一步地，所述步骤S12具体包括以下分步骤：
S121、通过D-H转换矩阵，建立下肢外骨骼末端轨迹关于髋关节、膝关节和踝关节对应的关节角角度的数学模型；
S122、根据步骤S113得到的髋关节、膝关节和踝关节各自对应的关节角角度的均值和方差，通过仿真确定下肢外骨骼末端轨迹的均值和方差；
S123、以CGA数据中关节角角度值的上、下限为约束，分别以下肢外骨骼末端轨迹的最小和最大为目标，建立求解下肢外骨骼末端轨迹上、下限的优化模型；
S124、求解步骤S123建立的优化模型，得到下肢外骨骼末端轨迹的上、下限。
进一步地，所述步骤S2包括：
S21、考虑下肢外骨骼失效时序，将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，计算髋关节失效概率；
S22、对四单元时序相关串联系统的膝关节，计算髋关节可靠情况下膝关节的失效概率；
S23、对四单元时序相关串联系统的踝关节，计算在髋关节和膝关节都可靠情况下踝关节的失效概率；
S24、对四单元时序相关串联系统的末端轨迹，计算在髋关节、膝关节和踝关节都可靠情况下末端轨迹的失效概率。
更进一步地，所述步骤S21具体包括：
S211、将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，所述的四时序单元依次为：髋关节、膝关节、踝关节以及末端轨迹；
S212、根据步骤S11得到的髋关节所对应的关节角角度的均值和方差，以及CGA数据中的髋关节所对应关节角角度的上、下限，得到髋关节失效概率。
进一步地，所述步骤S22具体包括以下分步骤：
S221、对步骤S211确定的四单元时序相关串联系统中的膝关节单元；
S222、根据步骤S11得到的髋关节和膝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节和膝关节所对应的各自关节角角度的上、下限，得到在髋关节可靠情况下膝关节的失效概率。
进一步地，所述步骤S23具体包括以下分步骤：
S231、对步骤S211确定的四单元时序相关串联系统中的踝关节单元；
S232、根据步骤S11得到的髋关节、膝关节和踝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节、膝关节和踝关节所对应的各自关节角角度的上、下限，根据容斥原理，忽略三个单元同时失效的情况，得到在髋关节和膝关节都可靠情况下踝关节的失效概率。
进一步地，所述步骤S24具体包括以下分步骤：
S241、对步骤S211确定的四单元时序相关串联系统中的踝关节单元；
S242、根据步骤S11得到的髋关节、膝关节和踝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节、膝关节和踝关节所对应的各自关节角角度的上、下限，以及步骤S2得到的下肢外骨骼末端轨迹的上、下限，根据容斥原理，忽略三个或三个以上单元同时失效的情况，得到在髋关节、膝关节和踝关节可靠情况下末端轨迹的失效概率。
进一步地，所述步骤S3具体为：根据步骤S2得到的髋关节失效概率、膝关节失效概率、踝关节失效概率以及末端轨迹失效概率，得到在不确定性条件下的下肢外骨骼时变可靠度。
本发明的有益效果：本发明通过简化下肢外骨骼模型，建立了独立的髋关节、膝关节和踝关节三个关节的简化模型，并充分考虑了下肢外骨骼存在的不确定性因素，建立了在不确定性条件下关节角的数学模型，通过分析得到了只考虑关节角情况下的下肢外骨骼时变可靠性模型，同时结合末端轨迹满足误差范围，实现下肢外骨骼时变可靠性计算，同时对全面提高下肢外骨骼的设计水平具有较高的理论支撑和工程实践意义。
附图说明
图1是本发明下肢外骨骼时变可靠性分析方法的流程图。
图2是本发明针对下肢外骨骼的简化模型。
图3是本发明实施例中在不确定条件下的关节角均值。
图4是本发明实施例中在不确定条件下的关节角方差。
图5是本发明实施例中在不确定条件下的末端轨迹均值。
图6是本发明实施例中在不确定条件下的末端轨迹方差。
图7是本发明实施例中的末端轨迹上、下限。
图8是本发明实施例中的髋关节失效概率。
图9是本发明实施例中的在髋关节可靠情况下膝关节的失效概率。
图10是本发明实施例中的在髋关节和膝关节都可靠情况下踝关节的失效概率。
图11是本发明实施例中的在髋关节、膝关节和踝关节都可靠情况下末端轨迹的失效概率。
图12是本发明实施例中的下肢外骨骼时变可靠度。
具体实施方式
为了对本发明做进一步的详细说明，下面再结合一具体实例对本发明的方案做一个说明。本实例以伯克利下肢外骨骼为实施例，在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
如表1所示，为本发明的方法所针对的伯克利下肢外骨骼三个关节液压缸杆件的安装尺寸。
表1伯克利下肢外骨骼三个关节液压缸杆件的安装尺寸
Figure BDA0001142025880000051
如图1所示为本发明的方案流程图，本发明采用的技术方案是：一种不确定性条件下的下肢外骨骼时变可靠性分析方法，针对如表1所示尺寸的伯克利下肢外骨骼液压缸杆件的安装尺寸，包括以下步骤：
S1、建立下肢外骨骼包含液压缸的简化模型，量化下肢外骨骼各关节杆件、运动副和液压缸的不确定性，计算在不确定性条件下各关节角角度、末端轨迹的均值和方差；具体包括：
S11、如图2所示，为本发明下肢外骨骼简化模型；图中H，K和A分别是下肢外骨骼髋关节、膝关节和踝关节，图中Lh、Lk和La分别是三个关节的液压缸简化模型，图中βh、βk和βa分别是下肢外骨骼三个关节的关节角，分析下肢外骨骼，建立下肢外骨骼髋、膝和踝三个关节对应的三个关节角角度各自在不确定性条件下的数学模型；确定各自关节角角度的均值和方差；
所述步骤S11具体包括以下分步骤：
S111、建立下肢外骨骼髋关节、膝关节和踝关节各自包含液压缸的简化模型，得到在满足外骨骼正常步态下的液压缸理想位移；
S112、分别量化髋关节、膝关节和踝关节的不确定性；所述的髋关节不确定性至少包括：髋关节杆件尺寸误差、髋关节液压缸误差以及髋关节运动副间隙误差；所述的膝关节不确定性至少包括：膝关节杆件尺寸误差、膝关节液压缸误差以及膝关节运动副间隙误差；所述的踝关节不确定性至少包括：踝关节杆件尺寸误差、踝关节液压缸误差以及踝关节运动副间隙误差。具体各关节不确定性因素如表2所示。
表2下肢外骨骼三个关节不确定性因素参数
不确定性因素 均值(mm) 方差(mm<sup>2</sup>) 分布类型
杆件尺寸误差 0 1/6 正态分布
髋关节液压缸误差 0 40 正态分布
膝关节液压缸误差 0 16 正态分布
踝关节液压缸误差 0 8 正态分布
运动副间隙误差 0.2 0.1 正态分布
S113、根据步骤S111得到的液压缸理想位移，建立下肢外骨骼髋关节、膝关节和踝关节各自对应的关节角角度在步骤S112得到的不确定性条件下的数学模型，并计算各自关节角角度的均值和方差，考虑如表2所示的不确定性因素，各自关节角角度的均值和方差分别如图3和图4所示。
关节角β(t)的计算公式如公式(1)所示：
Figure BDA0001142025880000061
公式(1)中，li(i＝1,2,3,4)表示液压缸杆件的安装尺寸，本申请实施例为便于理解，将各个关节的杆件统一简化为4部分，当i＝1,2,3,4，表示分别取不同部分杆件的尺寸大小；髋关节、膝关节和踝关节各自对应的液压缸安装尺寸如表1所示，当l1，l2，l3，l4取表1中的不同值时可对应得到βh或βk或βa；arctan(*)表示反正切函数；arccos(*)表示反余弦函数；L0表示液压缸的初始长度；ΔL(t)表示t时刻关节角角度满足人体正常步态对应的液压缸理想位移。
三个关节角的均值和方差可采用统一的公式计算，各关节角角度的均值统一表示为：μβ(t)；各关节角角度的方差统一表示为：
Figure BDA0001142025880000062
如公式(2)所示：
Figure BDA0001142025880000071
公式(2)中，
Figure BDA0001142025880000072
表示液压缸杆件安装尺寸误差的均值；
Figure BDA0001142025880000073
表示在运动副间隙误差的方差；
Figure BDA0001142025880000074
表示在运动副间隙误差的均值；
Figure BDA0001142025880000075
表示液压缸杆件安装尺寸误差的方差；
Figure BDA0001142025880000076
表示液压缸驱动源误差的均值；
Figure BDA0001142025880000077
表示液压缸驱动源误差的方差；β(l1,l2,l3,l4,L0,ΔL(t))表示在t时刻通过公式(1)求得的关节角数学模型在各个杆件取原始长度以及液压缸位移为理想位移时的取值；μβ(t)与
Figure BDA0001142025880000078
的公式中对β(l1,l2,l3,l4,L0,ΔL(t))进行简化表达为β，即
Figure BDA0001142025880000079
实际是对β(l1,l2,l3,l4,L0,ΔL(t))求偏导。
S12、通过D-H转换矩阵，建立下肢外骨骼末端轨迹关于髋关节、膝关节和踝关节对应的关节角角度的数学模型，并求解下肢外骨骼末端轨迹的上、下限；具体包括以下分步骤：
S121、通过D-H转换矩阵，建立下肢外骨骼末端轨迹关于髋关节、膝关节和踝关节对应的关节角角度的数学模型；末端轨迹的s(β(t))的计算公式如公式(3)所示，是关于三个关节的关节角角度函数，也就是一个三元函数：
s(β(t))＝x2(β(t))+y2(β(t))+z2(β(t)) 公式(3)
公式(3)中，x(β(t))表示由D-H转换矩阵确定的下肢外骨骼末端t时刻在基坐标系x方向的投影；y(β(t))表示由D-H转换矩阵确定的下肢外骨骼末端t时刻在基坐标系y方向的投影；z(β(t))表示由D-H转换矩阵确定的下肢外骨骼末端t时刻在基坐标系z方向的投影。
S122、根据步骤S113得到的髋关节、膝关节和踝关节对应的关节角角度的均值和方差，通过仿真确定下肢外骨骼末端轨迹的均值和方差；所述仿真为Matlab仿真，得到末端轨迹的均值和方差分别如图5和图6所示。
S123、以CGA(人体临床步态分析：Clinical gait analysis,CGA)数据中关节角角度的上、下限为约束，分别以下肢外骨骼末端轨迹的最小和最大为目标，建立求解下肢外骨骼末端轨迹上、下限的优化模型；末端轨迹上、下限的优化模型为：
Figure BDA0001142025880000081
Figure BDA0001142025880000082
公式(4)和(5)中，βmin(t)表示关节角角度β(t)在t时刻的下限；βmax(t)表示关节角角度β(t)在t时刻的上限。
S124、求解步骤S123建立的优化模型，得到下肢外骨骼末端轨迹的上、下限；结合人体步态的CGA数据，通过优化模型，得到如图7所示的末端轨迹上、下限。
S2、考虑下肢外骨骼失效时序，将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，分析每个单元的失效概率，并确定下肢外骨骼失效概率；所述的四单元依次为：髋关节、膝关节、踝关节以及末端轨迹；包括：
S21、计算髋关节失效概率；所述步骤S21具体包括以下分步骤：
S211、将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，所述的四时序单元依次为：髋关节、膝关节、踝关节以及末端轨迹；
S212、根据步骤S11得到的髋关节所对应的关节角角度的均值和方差，以及CGA数据中的髋关节对应的关节角角度的上、下限，得到髋关节失效概率，如图8所示；髋关节失效概率P1(t)的计算公式如公式(6)所示：
Figure BDA0001142025880000083
公式(6)中，Rh(t)表示t时刻髋关节的可靠度；
Figure BDA0001142025880000084
表示t时刻髋关节所对应关节角角度的上限；
Figure BDA0001142025880000085
表示t时刻髋关节所对应关节角角度的下限；fh(x,t)表示t时刻
Figure BDA0001142025880000086
的一维高斯分布概率密度函数，μh(t)表示髋关节的均值，
Figure BDA0001142025880000087
表示髋关节的方差。
S22、计算髋关节可靠情况下膝关节的失效概率；
S221、对步骤S211确定的四单元时序相关串联系统中的膝关节单元；
S222、根据步骤S11得到的髋关节和膝关节所对应的关节角的均值、方差，以及CGA数据中的髋关节和膝关节各自对应的关节角角度的上、下限，得到在髋关节可靠情况下膝关节的失效概率，如图9所示；在髋关节可靠情况下膝关节失效概率P2(t)的计算公式如公式(7)所示：
Figure BDA0001142025880000091
公式(7)中，R(βk|βh)表示t时刻在髋关节可靠情况下膝关节的可靠度；Rhk(t)表示在t时刻髋关节和膝关节同时可靠的概率，计算公式为
Figure BDA0001142025880000092
Figure BDA0001142025880000093
表示t时刻膝关节对应的关节角角度的上限；
Figure BDA0001142025880000094
表示t时刻髋关节对应的关节角角度的下限；fhk(x,y,t)表示t时刻
Figure BDA0001142025880000095
的二维高斯分布概率密度函数，μk(t)表示膝关节的均值，
Figure BDA0001142025880000096
表示膝关节的方差；rhk(t)表示髋关节和膝关节联合分布fhk(x,y,t)的相关系数；rhk(t)是通过在不确定条件下髋关节和膝关节各自方差和均值产生的两列随机数求得的。
S23、计算在髋关节和膝关节都可靠情况下踝关节的失效概率；
S231、对步骤S211确定的四单元时序相关串联系统中的踝关节单元；
S232、根据步骤S11得到的髋关节、膝关节和踝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节、膝关节和踝关节各自对应的关节角角度的上、下限，根据容斥原理，忽略三个单元同时失效的情况，得到在髋关节和膝关节都可靠情况下踝关节的失效概率，如图10所示；在髋关节和膝关节都可靠情况下踝关节的失效概率P3(t)的计算公式如公式(8)所示：
Figure BDA0001142025880000097
公式(8)中，R(βa|βhβk)表示t时刻在髋关节和膝关节都可靠的情况下踝关节的可靠度；Rha(t)表示在t时刻髋关节和踝关节同时可靠的概率；Rka(t)表示在t时刻膝关节和踝关节同时可靠的概率；Rk(t)表示t时刻髋关节的可靠度；Ra(t)表示t时刻踝关节的可靠度。
S24、计算在髋关节、膝关节和踝关节都可靠情况下末端轨迹的失效概率；
S241、对步骤S211确定的四单元时序相关串联系统中的踝关节单元；
S242、根据步骤S11得到的髋关节、膝关节和踝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节、膝关节和踝关节各自对应的关节角角度的上、下限，以及步骤S2得到的下肢外骨骼末端轨迹的上、下限，根据容斥原理，忽略三个或三个以上单元同时失效的情况，得到在髋关节、膝关节和踝关节可靠情况下末端轨迹的失效概率，如图11所示；在髋关节、膝关节和踝关节都可靠情况下末端轨迹的失效概率P4(t)的计算公式如公式(9)所示：
Figure BDA0001142025880000101
公式(9)中，R(s|βhβkβa)表示t时刻在髋关节、髋关节以及踝关节都可靠的情况下末端轨迹的可靠度；Rhs(t)表示在t时刻髋关节和末端轨迹同时可靠的概率；Rks(t)表示在t时刻膝关节和末端轨迹同时可靠的概率；Ras(t)表示在t时刻踝关节和末端轨迹同时可靠的概率；Rs(t)表示t时刻末端轨迹的可靠度。
S3、根据步骤S21得到的髋关节失效概率P1(t)、步骤S22得到的膝关节失效概率P2(t)、步骤S23得到的踝关节失效概率P3(t)以及步骤S24得到的末端轨迹失效概率P4(t)，得到在不确定性条件下的下肢外骨骼时变可靠度，如图12所示。
下肢外骨骼的可靠度计算公式如公式(10)所示：
R(t)＝1-[P1(t)+P2(t)+P3(t)+P4(t)] 公式(10)
公式(10)中，R(t)表示t时刻下肢外骨骼的可靠度。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。

Claims (8)

1.一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，包括：
S1、建立下肢外骨骼包含液压缸的简化模型，量化下肢外骨骼各关节杆件、运动副和液压缸的不确定性，计算在不确定性条件下各关节角角度以及末端轨迹的均值和方差；所述步骤S1包括：
S11、分析下肢外骨骼，建立下肢外骨骼髋关节、膝关节和踝关节对应的关节角角度在各自不确定性条件下的数学模型；并确定各自关节角角度的均值和方差；
S12、通过D-H转换矩阵，建立末端轨迹关于髋关节、膝关节和踝关节对应的关节角角度的数学模型，并求解下肢外骨骼末端轨迹的上、下限；所述步骤S12具体包括以下分步骤：
S121、通过D-H转换矩阵，建立下肢外骨骼末端轨迹关于髋关节、膝关节和踝关节对应的关节角角度的数学模型；
S122、根据步骤S113得到的关节角的均值和方差，通过仿真确定下肢外骨骼末端轨迹的均值和方差；
S123、以CGA数据中关节角角度的上、下限为约束，分别以下肢外骨骼末端轨迹的最小和最大为目标，建立求解下肢外骨骼末端轨迹上、下限的优化模型；
S124、求解步骤S123建立的优化模型，得到下肢外骨骼末端轨迹的上、下限；
S2、考虑下肢外骨骼失效时序，将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，分析每个单元的失效概率，并确定下肢外骨骼失效概率；所述的四单元依次为：髋关节、膝关节、踝关节以及末端轨迹；
S3、根据步骤S2所得的下肢外骨骼失效概率确定下肢外骨骼的可靠度。
2.根据权利要求1所述的一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，所述步骤S11具体包括以下分步骤：
S111、建立下肢外骨骼髋关节、膝关节和踝关节各自包含液压缸的简化模型，得到在满足外骨骼正常步态下的液压缸理想位移；
S112、分别量化髋关节、膝关节和踝关节的不确定性；
S113、根据步骤S111得到的液压缸理想位移，建立下肢外骨骼髋关节、膝关节和踝关节各自对应的关节角角度在步骤S112得到的不确定性条件下的数学模型，并确定各自关节角角度的均值和方差。
3.根据权利要求2所述的一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，步骤S112所述的髋关节不确定性至少包括：髋关节杆件尺寸误差、髋关节液压缸误差以及髋关节运动副间隙误差；
所述的膝关节不确定性至少包括：膝关节杆件尺寸误差、膝关节液压缸误差以及膝关节运动副间隙误差；
所述的踝关节不确定性至少包括：踝关节杆件尺寸误差、踝关节液压缸误差以及踝关节运动副间隙误差。
4.根据权利要求1所述的一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，所述步骤S2包括：
S21、考虑下肢外骨骼失效时序，将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，计算髋关节失效概率；
S22、对四单元时序相关串联系统的膝关节，计算髋关节可靠情况下膝关节的失效概率；
S23、对四单元时序相关串联系统的踝关节，计算在髋关节和膝关节都可靠情况下踝关节的失效概率；
S24、对四单元时序相关串联系统的末端轨迹，计算在髋关节、膝关节和踝关节都可靠情况下末端轨迹的失效概率。
5.根据权利要求4所述的一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，所述步骤S21具体包括：
S211、将下肢外骨骼从髋关节到末端轨迹看作自上而下的四单元时序相关串联系统，所述的四单元依次为：髋关节、膝关节、踝关节以及末端轨迹；
S212、根据步骤S11得到的髋关节所对应的关节角角度的均值和方差，以及CGA数据中的髋关节所对应关节角角度的上、下限，得到髋关节失效概率。
6.根据权利要求4所述的一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，所述步骤S22具体包括以下分步骤：
S221、对步骤S211确定的四单元时序相关串联系统中的膝关节单元；
S222、根据步骤S11得到的髋关节和膝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节和膝关节所对应的各自关节角角度的上、下限，得到在髋关节可靠情况下膝关节的失效概率。
7.根据权利要求4所述的一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，所述步骤S23具体包括以下分步骤：
S231、对步骤S211确定的四单元时序相关串联系统中的踝关节单元；
S232、根据步骤S11得到的髋关节、膝关节和踝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节、膝关节和踝关节所对应的各自关节角角度的上、下限，根据容斥原理，忽略三个单元同时失效的情况，得到在髋关节和膝关节都可靠情况下踝关节的失效概率。
8.根据权利要求4所述的一种不确定性条件下的下肢外骨骼时变可靠性分析方法，其特征在于，所述步骤S24具体包括以下分步骤：
S241、对步骤S211确定的四单元时序相关串联系统中的踝关节单元；
S242、根据步骤S11得到的髋关节、膝关节和踝关节所对应的关节角角度的均值、方差，以及CGA数据中的髋关节、膝关节和踝关节所对应的各自关节角角度值的上、下限，以及步骤S2得到的下肢外骨骼末端轨迹的上、下限，根据容斥原理，忽略三个或三个以上单元同时失效的情况，得到在髋关节、膝关节和踝关节可靠情况下末端轨迹的失效概率。
CN201610953616.XA 2016-11-03 2016-11-03 一种不确定性条件下的下肢外骨骼时变可靠性分析方法 Active CN106503364B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610953616.XA CN106503364B (zh) 2016-11-03 2016-11-03 一种不确定性条件下的下肢外骨骼时变可靠性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610953616.XA CN106503364B (zh) 2016-11-03 2016-11-03 一种不确定性条件下的下肢外骨骼时变可靠性分析方法

Publications (2)

Publication Number Publication Date
CN106503364A CN106503364A (zh) 2017-03-15
CN106503364B true CN106503364B (zh) 2020-07-31

Family

ID=58322297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610953616.XA Active CN106503364B (zh) 2016-11-03 2016-11-03 一种不确定性条件下的下肢外骨骼时变可靠性分析方法

Country Status (1)

Country Link
CN (1) CN106503364B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108634934B (zh) * 2018-05-07 2021-01-29 北京长木谷医疗科技有限公司 对脊柱矢状位图像进行处理的方法和设备
CN110502866A (zh) * 2019-08-30 2019-11-26 电子科技大学 一种考虑输入误差的连杆机构运动可靠性分析方法
CN113012812B (zh) * 2021-02-03 2023-02-21 上海橙捷健康科技有限公司 一种膝、踝关节以及足底压力数据整合方法及系统
CN115401699B (zh) * 2022-10-31 2023-03-24 广东隆崎机器人有限公司 工业机器人精度可靠性分析方法、装置、设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150001658A (ko) * 2013-06-26 2015-01-06 다솔 시스템즈 시뮬리아 코포레이션 유한 요소 분석, 프로세스 통합, 및 설계 최적화를 이용한 근골격계 모델링
CN205198395U (zh) * 2015-12-22 2016-05-04 中国科学院沈阳自动化研究所 一种可穿戴式下肢外骨骼助力行走机器人机构
CN105662789A (zh) * 2016-03-30 2016-06-15 北京精密机电控制设备研究所 一种基于健全肢体运动监测、控制残肢运动的外骨骼系统
CN205459236U (zh) * 2016-01-08 2016-08-17 中国人民解放军理工大学 可折叠便携式下肢外骨骼
CN105902368A (zh) * 2016-03-01 2016-08-31 芜湖安普机器人产业技术研究院有限公司 连杆关节一体化液压驱动外骨骼控制系统及其控制方法
CN105943316A (zh) * 2016-05-23 2016-09-21 成都润惠科技有限公司 一种具有回力结构的人体下肢外骨骼

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150001658A (ko) * 2013-06-26 2015-01-06 다솔 시스템즈 시뮬리아 코포레이션 유한 요소 분석, 프로세스 통합, 및 설계 최적화를 이용한 근골격계 모델링
CN205198395U (zh) * 2015-12-22 2016-05-04 中国科学院沈阳自动化研究所 一种可穿戴式下肢外骨骼助力行走机器人机构
CN205459236U (zh) * 2016-01-08 2016-08-17 中国人民解放军理工大学 可折叠便携式下肢外骨骼
CN105902368A (zh) * 2016-03-01 2016-08-31 芜湖安普机器人产业技术研究院有限公司 连杆关节一体化液压驱动外骨骼控制系统及其控制方法
CN105662789A (zh) * 2016-03-30 2016-06-15 北京精密机电控制设备研究所 一种基于健全肢体运动监测、控制残肢运动的外骨骼系统
CN105943316A (zh) * 2016-05-23 2016-09-21 成都润惠科技有限公司 一种具有回力结构的人体下肢外骨骼

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
仿生下肢外骨骼机器人结构设计及可靠性分析;陈伟;《中国优秀硕士学位论文全文数据库医药卫生科技辑》;20160215(第02期);第E080-82页 *

Also Published As

Publication number Publication date
CN106503364A (zh) 2017-03-15

Similar Documents

Publication Publication Date Title
CN106503364B (zh) 一种不确定性条件下的下肢外骨骼时变可靠性分析方法
Kong et al. Quantifying the physical intensity of construction workers, a mechanical energy approach
Blana et al. Real-time simulation of hand motion for prosthesis control
Li et al. Zeroing dynamics method for motion control of industrial upper-limb exoskeleton system with minimal potential energy modulation
Kim et al. Prediction and analysis of human motion dynamics performing various tasks
Sun et al. Sensor reduction, estimation, and control of an upper-limb exoskeleton
Ayusawa et al. Simultaneous global inverse kinematics and geometric parameter identification of human skeletal model from motion capture data
Sun et al. A novel RBF neural network-based iterative learning control for lower limb rehabilitation robot with strong robustness
Vatankhah et al. Intermittent control model for ascending stair biped robot using a stable limit cycle model
Asghari et al. A planar neuro-musculoskeletal arm model in post-stroke patients
Wang et al. A forward-inverse dynamics modeling framework for human musculoskeletal multibody system
CN112650079B (zh) 一种七自由度机器人的运动学反解方法
Pontonnier et al. Inverse dynamics method using optimization techniques for the estimation of muscles forces involved in the elbow motion
KR20120070660A (ko) 다물체 모델링과 통계 방법론을 이용한 인체 하지부 구속력의 불확실성 추정 방법
CN109887570B (zh) 一种基于rgb-d相机与imu传感器的机器人辅助康复训练方法
Pereira et al. Implementation of an efficient muscle fatigue model in the framework of multibody systems dynamics for analysis of human movements
Xie et al. A new criterion for redundancy resolution of human arm in reaching tasks
Gragg et al. Optimization-based posture reconstruction for digital human models
Bueno et al. Neuromusculoskeletal model self-calibration for on-line sequential Bayesian moment estimation
Wei et al. Development of an upper limb rehabilitation robot system for bilateral training
CN111772638A (zh) 步态控制方法、装置、电子设备及存储介质
Jacquier-Bret et al. Effect of spinal cord injury at C6–C7 on global upper-limb coordination during grasping: manipulability approach
Xie et al. [Retracted] Sensor‐Based Exercise Rehabilitation Robot Training Method
Xi et al. Finite Element Analysis of Series Elastic Actuator for Exoskeleton Robot Joint
Lura et al. Robotic model for simulating upper body movement

Legal Events

Date Code Title Description
C06 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