CN111024292B - 一种连续步态下单步作用力状态测试提取方法 - Google Patents
一种连续步态下单步作用力状态测试提取方法 Download PDFInfo
- Publication number
- CN111024292B CN111024292B CN201911306438.1A CN201911306438A CN111024292B CN 111024292 B CN111024292 B CN 111024292B CN 201911306438 A CN201911306438 A CN 201911306438A CN 111024292 B CN111024292 B CN 111024292B
- Authority
- CN
- China
- Prior art keywords
- force
- gait
- threshold
- search
- value
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/16—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring several components of force
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2455—Query execution
- G06F16/24564—Applying rules; Deductive queries
- G06F16/24565—Triggers; Constraints
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/283—Multi-dimensional databases or data warehouses, e.g. MOLAP or ROLAP
Abstract
本发明公开了一种连续步态下单步作用力状态测试提取方法,包括步骤:同步采集并按时间顺序存储各维力及其力矩值;搜索查找所存储的地面反应力Fz的最大值Fp和所对应的存储位置Ip;判定Fp是否大于已设定的触发门限Fmin,如大于进入下一步,否则结束;设置搜索门限FSmin,从Ip开始,分别向前、向后遍历搜索,直到其Fz<FSmin,其对应点位置分别确定为Ia和Ib;提取所有维度测试力和测试力矩在Ia和Ib间数据作为此次单步态值,再把所有在Ia和Ib间的测试数据删去;重复前述搜索查找地面反应力Fz的最大值,直至跳出循环。本发明利用所测竖向地面反应力Fz的非负特征,对单个步态力状态进行有效提取,该技术在解决人体步态测试力的科研和工程问题上具有非常实用的价值。
Description
技术领域
本发明涉及作用力状态测试提取方法,特别是涉及一种连续步态下单步作用力状态测试提取方法。
背景技术
随着社会的发展,建筑师在土木工程设计过程中不断追求结构纤细美观的外表,其带来的问题是,结构体系的尺度在不断增大,结构的谐振频率和阻尼耗能能力在不断减小,当频率减小到一定程度,且接近于人的激励频率时,容易导致共振发生,由此带来很多结构使用的不适定以及结构被破坏的问题。为探索结构在行人作用下的振动机理,工程师和研究学者需要测试人的双足对结构所产生的反作用力,以提取单个步态下人对结构的荷载激励特征。在体育运动学领域,对运动员单步作用力状态测试研究,是改进训练方法重要依据,从而对提高运动员训练效率和提升竞技水平产生积极影响。
要获得单步作用力状态准确的测量结果,就必须要求测试仪器具备高灵敏度和大的测量动态范围。高灵敏度用于满足精度要求,大的测量动态范围用于满足体型差异很大的各类人员步态测试需求。然而,受周围毗连空间结构以及周围电磁环境等的影响,必然导致高灵敏度和大的测量动态范围的测试仪器的测量结果中含有许多噪声。比如,当没有荷载作用时候,测量仪器依然显示有微量的荷载波动,这种波动使得在提取单步态信号起始时刻点产生困难,而载荷因人而异,差异也较大,有关作用力状态提取的准确性就更不好保证,使得依据作用力状态开展的研究变得困难。
发明内容
发明目的:针对现有技术不足,本发明提出了一种连续步态下单步作用力状态测试提取方法,利用所测单步作用力状态中的竖向地面反应力Fz的非负特征,并利用空间单个步态所产生的空间三维力和对应的力矩具有时间同步特征,实现对单个步态值进行科学快速的有效提取。通过关联反应力Fz,自适应设置搜索门限值,解决了测量载荷因人的差异造成的测量误差。
技术方案:本发明所述一种连续步态下单步作用力状态测试提取方法,包括以下步骤:
(1)利用测量装置同步采集并按时间顺序存储各维力和各维力矩,其中必包含竖向地面反应力Fz;
(2)从存储的数据中搜索查找竖向地面反应力Fz的最大值Fp和所对应的位置Ip;
(3)判定Fp是否满足已设定的触发门限Fmin,如满足进入下一步,否则结束;
(4)设置搜索门限FSmin,从Ip开始,向前进行遍历搜索,直到Fz≤FSmin,把该点位置确定为Ia;
(5)从Ip开始,向后进行遍历搜索,直到Fz≤FSmin,把该点位置确定为Ib;
(6)提取所有维度力和力矩在Ia和Ib间的采集数据作为此次单步态值并存储和显示;
(7)把所有在Ia和Ib间的采集数据删除,重复步骤(2)直至跳出循环,结束。
进一步地,所述步骤(1)还包括:同步采集并存储垂直于Fz的且彼此垂直的作用力Fx、Fy以及对应的空间三维力矩Mx、My和Mz,构成单个步态所产生的空间六维力;其中,所述Fx为步态力在行走方向上的力;Mx、My和Mz分别表示x、y、z方向上的力矩,x方向为行走方向。这就让更多型号的力和力矩测量装置可以被采用,也可以根据科研需要选择测量的维度。
进一步地,在步骤(3)中,所述触发门限Fmin的设定是通过实验获得的,选择应用空间中体重最轻的人产生作用力获取其采集数据波形,并综合考虑波形噪声的影响来确定触发门限的值。这就保证了分析采集数据过程的科学性和数据分析的有效性。
进一步地,所述搜索门限FSmin理想条件下取值为0。
进一步地,所述搜索门限FSmin根据竖向地面反应力Fz的最大值Fp自适应选取。
进一步地,步骤(4)与步骤(5)中的搜索门限FSmin的取值不一样,步骤(4)中所述搜索门限FSmin与Fp的百分比取值较高,步骤(5)中所述搜索门限FSmin与Fp的百分比取值较低。
进一步地,步骤(4)中所述搜索门限FSmin与Fp的百分比取值范围为15%-30%。
进一步地,步骤(5)中所述搜索门限FSmin与Fp的百分比取值范围为0-10%。
事实上,地面反应力Fz的起始波形陡峭,而终止的波形较平缓。所以,上述门限值设置的改进进一步提高了有关作用力及力矩信息提取的准确性。
有益效果:本发明利用所测竖向地面反应力Fz的非负特征,并利用空间单个步态所产生的空间多维力及力矩所具有的时间上的同步特征,能够对单个步态力特征进行科学快速的有效提取,可实现不同人连续步态下单步作用力状态自动测试提取,减少了测试提取时间与步骤,测试提取的精度高。本发明中的搜索门限FSmin是可自适应控制的,可以根据不同体重人在连续步态下单步作用力状态产生合适搜索门限,进一步提高测试提取精确度。该技术在解决人体步态力特征测试的科研问题和工程问题上具有非常实用的价值。
附图说明
图1是本发明的方法流程图;
图2是本发明的测试结构模型图;
图3是本发明的部分地面反应力及力矩原始采集波形图;
图4是本发明的某测试数据的单步态提取结果图。
具体实施方式
下面结合附图和具体实施方式,对本发明的技术方案作进一步的介绍。
图1是本发明的方法流程图,包括以下步骤:
a、应用相关测量装置同步采集并存储各维力和力矩值,竖向地面反应力Fz必包含其中;
b、从存储的数据中搜索查找地面反应力Fz的最大值Fp和所对应的位置Ip;
c、判定Fp是否满足已设定的触发门限Fmin,如满足进入下一步,否则结束;
在一个实施例中,所述触发门限Fmin的设定是通过实验获得的,选择应用空间中体重最轻的人产生作用力获取其采集数据波形,并综合考虑波形噪声的影响来确定触发门限的值。
d、设置搜索门限FSmin,从Ip开始,向前进行遍历搜索,直到Fz≤FSmin,把该点位置确定为Ia;
e、从Ip开始,向后进行遍历搜索,直到Fz≤FSmin,把该点位置确定为Ib;
理想条件下,搜索门限FSmin取值为0。而考虑到其它影响因素,显示表现为波形图存在噪声,搜索门限FSmin应当根据竖向地面反应力Fz的最大值Fp自适应选取。步骤d与步骤e中的搜索门限FSmin的取值不一样,步骤d中所述搜索门限FSmin与Fp的百分比取值较高,步骤e中所述搜索门限FSmin与Fp的百分比取值较低。
f、提取所有维度测试力和测试力矩在Ia和Ib间采集数据作为此次单步态力特征值并存储和显示;
g、把所有在Ia和Ib间采集的数据删除,重复步骤b直至跳出循环,结束。本发明利用所测竖向地面反应力Fz的非负特征,并利用空间单个步态所产生的空间多维力具有同步特征,对单个步态力特征进行科学快速的有效提取,并能够自动完成。
图2是本发明的测试结构模型图,其中包括:展示了行人在测力板上行走时的空间测试反应力状态:分别为空间三维力Fx,Fy和Fz,以及空间三维力矩Mx,My和Mz。
图3是本发明的部分地面反应力及力矩原始采集波形图。其中,横坐标表示时间,单位s为秒。纵坐标Fx(N)、Fy(N)、Fz(N)分别表示x、y、z方向上的力,单位N为牛顿;纵坐标Mx(Nm)、My(Nm)、Mz(Nm)分别表示x、y、z方向上的力矩,单位Nm为牛顿·米。
图4是本发明的某测试数据的单步态提取结果图。表述了一位重69公斤、高1.82米的男人,在测力板上行走时的空间测试反应力状态提取结果。本实施例考虑的是理想情况,即FSmin=0。其中,横坐标表示的时间,单位s为秒。纵坐标Fx(N)、Fy(N)、Fz(N)分别表示x、y、z方向上的力,单位N为牛顿;纵坐标Mx(Nm)、My(Nm)、Mz(Nm)分别表示x、y、z方向上的力矩,单位Nm为牛顿·米。由图可知本方法能够很好的提取连续行走下的单步态激励荷载。由于地面反应力Fz的起始波形陡峭,而终止的波形较平缓,因此建议步骤d中搜索门限FSmin与Fp的百分比取值范围为15%-30%,步骤e中所述搜索门限FSmin与Fp的百分比取值范围为0-10%。上述门限值设置的改进进一步提高了有关作用力及力矩信息提取的准确性。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内还可以对其作出种种变化。例如,在上述实施方式中,测量装置在测试反应力状态时各种关于采集和存储参数设置的改变,以及对所测试反应力状态统计分析等。
Claims (7)
1.一种连续步态下单步作用力状态测试提取方法,其特征在于,包括以下步骤:
(1)利用测量装置同步采集并按时间顺序存储各维力和各维力矩,其中必包含竖向地面反应力Fz;
(2)从存储的数据中搜索查找竖向地面反应力Fz的最大值Fp和所对应的位置Ip;
(3)判定Fp是否满足已设定的触发门限Fmin,如满足进入下一步,否则结束;
(4)设置搜索门限FSmin,从Ip开始,向前进行遍历搜索,直到Fz≤FSmin,把该点位置确定为Ia;
(5)从Ip开始,向后进行遍历搜索,直到Fz≤FSmin,把该点位置确定为Ib;
(6)提取所有维度力和力矩在Ia和Ib间的采集数据作为此次单步态值并存储和显示;
(7)把所有在Ia和Ib间的采集数据删去,重复步骤(2)直至跳出循环,结束。
2.根据权利要求1所述的连续步态下单步作用力状态测试提取方法,其特征在于,所述步骤(1)包括:同步采集并存储垂直于Fz的且彼此垂直的作用力Fx、Fy以及对应的空间三维力矩Mx、My和Mz,构成单个步态所产生的空间六维力;其中,所述Fx为步态力在行走方向上的力;Mx、My和Mz分别表示x、y、z方向上的力矩,x方向为行走方向。
3.根据权利要求1所述的连续步态下单步作用力状态测试提取方法,其特征在于,在步骤(3)中,所述触发门限Fmin的设定是通过实验获得的,选择应用空间中体重最轻的人产生作用力获取其采集数据波形,并综合考虑波形噪声的影响来确定触发门限的值。
4.根据权利要求1所述的连续步态下单步作用力状态测试提取方法,其特征在于,所述搜索门限FSmin理想条件下取值为0。
5.根据权利要求1所述的连续步态下单步作用力状态测试提取方法,其特征在于,所述搜索门限FSmin根据竖向地面反应力Fz的最大值Fp自适应选取。
6.根据权利要求5所述的连续步态下单步作用力状态测试提取方法,其特征在于,步骤(4)中所述搜索门限FSmin与Fp的百分比取值范围为15%-30%。
7.根据权利要求5所述的连续步态下单步作用力状态测试提取方法,其特征在于,步骤(5)中所述搜索门限FSmin与Fp的百分比取值范围为0-10%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911306438.1A CN111024292B (zh) | 2019-12-18 | 2019-12-18 | 一种连续步态下单步作用力状态测试提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911306438.1A CN111024292B (zh) | 2019-12-18 | 2019-12-18 | 一种连续步态下单步作用力状态测试提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111024292A CN111024292A (zh) | 2020-04-17 |
CN111024292B true CN111024292B (zh) | 2021-07-09 |
Family
ID=70210239
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911306438.1A Active CN111024292B (zh) | 2019-12-18 | 2019-12-18 | 一种连续步态下单步作用力状态测试提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111024292B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113273980A (zh) * | 2021-06-07 | 2021-08-20 | 深圳麦格米特电气股份有限公司 | 一种获取生命体征数据的方法、光纤传感器及智能床垫 |
CN115553757B (zh) * | 2022-12-06 | 2023-04-07 | 深圳市第二人民医院(深圳市转化医学研究院) | 步态过程中的有效地面反力智能筛选方法及设备 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101347333A (zh) * | 2007-07-17 | 2009-01-21 | 中国科学院理化技术研究所 | 记录人体步态信息的微型测量装置 |
CN102175356A (zh) * | 2011-02-22 | 2011-09-07 | 中国科学院合肥物质科学研究院 | 一种步态触觉信息测量装置及测量方法 |
CN104015829A (zh) * | 2014-06-10 | 2014-09-03 | 西北工业大学 | 一种爬壁机器人的负压吸附式双足 |
CN104921851A (zh) * | 2015-05-25 | 2015-09-23 | 河北工业大学 | 主动型膝上假肢膝关节的预测控制方法 |
CN105608309A (zh) * | 2015-12-11 | 2016-05-25 | 杭州南江机器人股份有限公司 | 一种双足机器人步行规划和控制方法 |
CN106020215A (zh) * | 2016-05-09 | 2016-10-12 | 北京航空航天大学 | 一种基于单步预测矩阵博弈的近距空战自主决策方法 |
CN106595912A (zh) * | 2016-11-01 | 2017-04-26 | 中国海洋大学 | 一种人足动态力学实时检测评估系统及方法 |
CN107095677A (zh) * | 2016-02-19 | 2017-08-29 | 财团法人资讯工业策进会 | 步态分析系统及方法 |
JP2018017585A (ja) * | 2016-07-27 | 2018-02-01 | 富士通株式会社 | 誘導システム及び誘導制御装置 |
CN107703747A (zh) * | 2017-10-09 | 2018-02-16 | 东南大学 | 一种面向搅拌摩擦焊应用的重载机器人动力学参数自标定方法 |
US9916011B1 (en) * | 2015-08-22 | 2018-03-13 | Bertec Corporation | Force measurement system that includes a force measurement assembly, a visual display device, and one or more data processing devices |
CN108058758A (zh) * | 2017-12-07 | 2018-05-22 | 东莞深圳清华大学研究院创新中心 | 基于六自由度腿的机器人双足行走仿真评估系统及评估方法 |
CN109514602A (zh) * | 2018-11-08 | 2019-03-26 | 东南大学 | 一种基于负载自适应辨识的工业机器人力矩补偿控制方法 |
CN109846672A (zh) * | 2019-03-07 | 2019-06-07 | 合肥工业大学 | 可变刚度的踝足康复矫形器及其运动控制方法 |
-
2019
- 2019-12-18 CN CN201911306438.1A patent/CN111024292B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101347333A (zh) * | 2007-07-17 | 2009-01-21 | 中国科学院理化技术研究所 | 记录人体步态信息的微型测量装置 |
CN102175356A (zh) * | 2011-02-22 | 2011-09-07 | 中国科学院合肥物质科学研究院 | 一种步态触觉信息测量装置及测量方法 |
CN104015829A (zh) * | 2014-06-10 | 2014-09-03 | 西北工业大学 | 一种爬壁机器人的负压吸附式双足 |
CN104921851A (zh) * | 2015-05-25 | 2015-09-23 | 河北工业大学 | 主动型膝上假肢膝关节的预测控制方法 |
US9916011B1 (en) * | 2015-08-22 | 2018-03-13 | Bertec Corporation | Force measurement system that includes a force measurement assembly, a visual display device, and one or more data processing devices |
CN105608309A (zh) * | 2015-12-11 | 2016-05-25 | 杭州南江机器人股份有限公司 | 一种双足机器人步行规划和控制方法 |
CN107095677A (zh) * | 2016-02-19 | 2017-08-29 | 财团法人资讯工业策进会 | 步态分析系统及方法 |
CN106020215A (zh) * | 2016-05-09 | 2016-10-12 | 北京航空航天大学 | 一种基于单步预测矩阵博弈的近距空战自主决策方法 |
JP2018017585A (ja) * | 2016-07-27 | 2018-02-01 | 富士通株式会社 | 誘導システム及び誘導制御装置 |
CN106595912A (zh) * | 2016-11-01 | 2017-04-26 | 中国海洋大学 | 一种人足动态力学实时检测评估系统及方法 |
CN107703747A (zh) * | 2017-10-09 | 2018-02-16 | 东南大学 | 一种面向搅拌摩擦焊应用的重载机器人动力学参数自标定方法 |
CN108058758A (zh) * | 2017-12-07 | 2018-05-22 | 东莞深圳清华大学研究院创新中心 | 基于六自由度腿的机器人双足行走仿真评估系统及评估方法 |
CN109514602A (zh) * | 2018-11-08 | 2019-03-26 | 东南大学 | 一种基于负载自适应辨识的工业机器人力矩补偿控制方法 |
CN109846672A (zh) * | 2019-03-07 | 2019-06-07 | 合肥工业大学 | 可变刚度的踝足康复矫形器及其运动控制方法 |
Non-Patent Citations (4)
Title |
---|
A Three-Dimensional Pedestrian–Structure Interaction Model for;Yan An Gao;《International Journal of Structural Stability and Dynamics》;20181231;正文第1-28页 * |
Estimation of Ground Reaction Forces and Moments;Angelos Karatsidis;《Sensors》;20171231;正文第1-22页 * |
基于步态触觉特征的生物特征识别;徐苏;《中国博士学位论文全文数据库 信息科技辑》;20100930(第9期);第13-81页 * |
板结构在行人激励作用下的动力行为研究;高延安;《第十三届全国振动理论及应用学术会议论文集》;20191109;正文第1-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111024292A (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111024292B (zh) | 一种连续步态下单步作用力状态测试提取方法 | |
Rainieri et al. | Development and validation of an automated operational modal analysis algorithm for vibration-based monitoring and tensile load estimation | |
Dossett et al. | Spatial curvature and cosmological tests of general relativity | |
CN104983489B (zh) | 下肢假肢行走时的路况识别方法 | |
CN104834004B (zh) | 基于峰值前后波形斜率的矿山微震与爆破信号识别方法 | |
CN105910601A (zh) | 一种基于隐马尔科夫模型的室内地磁定位方法 | |
US11789173B1 (en) | Real-time microseismic magnitude calculation method and device based on deep learning | |
CN108875178A (zh) | 用于减小结构模态识别不确定性的传感器布置方法 | |
Li et al. | An automated operational modal analysis algorithm and its application to concrete dams | |
CN110232412A (zh) | 一种基于多模态深度学习的人体步态预测方法 | |
CN110322509A (zh) | 基于层级类激活图的目标定位方法、系统及计算机设备 | |
CN106842161B (zh) | 基于置信度评价的微动周期提取方法 | |
CN102539537B (zh) | 基于附加虚拟质量的框架结构损伤方法 | |
Gong et al. | Empirical formula of fundamental period for steel structure based on seismic response record | |
JP2016161463A (ja) | 液状化強度曲線生成方法、装置及びプログラム | |
Hölter et al. | Optimal sensor location for parameter identification in soft clay | |
CN102551731B (zh) | 一种基于数据曲线比较的跌倒活动检测方法 | |
Ping et al. | Pedestrian stride-length estimation based on bidirectional LSTM network | |
Kalinowska et al. | Behaviour of cohesive soil subjected to low-frequency cyclic loading in strain-controlled tests | |
TWI467212B (zh) | 建築物樓層之地震即時分析系統及其方法與儲存媒體 | |
CN105233459B (zh) | 跳远训练装置及其方法 | |
CN106033490B (zh) | 一种基于卫星遥感技术的岩石风化程度提取方法 | |
Guoqing et al. | A Study of Fall Risk Assessment in the Elderly | |
Balachandra et al. | Validating SSI Response of Buildings on Liquefiable Deposits using Centrifuge Testing and Numerical Modelling | |
Özçelik et al. | Modal identification results of RC frames at different damage levels |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20200417 Assignee: HUAI'AN XINYUANMING NEW MATERIAL Co.,Ltd. Assignor: HUAIYIN INSTITUTE OF TECHNOLOGY Contract record no.: X2021980012804 Denomination of invention: A test and extraction method of single step force state under continuous gait Granted publication date: 20210709 License type: Common License Record date: 20211122 |
|
EE01 | Entry into force of recordation of patent licensing contract |