CN107736890B - 不同行走任务下膝关节内侧负载的估值方法 - Google Patents
不同行走任务下膝关节内侧负载的估值方法 Download PDFInfo
- Publication number
- CN107736890B CN107736890B CN201711054022.6A CN201711054022A CN107736890B CN 107736890 B CN107736890 B CN 107736890B CN 201711054022 A CN201711054022 A CN 201711054022A CN 107736890 B CN107736890 B CN 107736890B
- Authority
- CN
- China
- Prior art keywords
- knee joint
- kam
- period
- kmf
- kfm
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/112—Gait analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/22—Ergometry; Measuring muscular strength or the force of a muscular blow
- A61B5/224—Measuring muscular strength
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/22—Ergometry; Measuring muscular strength or the force of a muscular blow
- A61B5/224—Measuring muscular strength
- A61B5/227—Measuring muscular strength of constricting muscles, i.e. sphincters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4528—Joints
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4538—Evaluating a particular part of the muscoloskeletal system or a particular medical condition
- A61B5/4585—Evaluating the knee
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
Abstract
本发明涉及医疗器械、生物力学测量，为确定不同不行任务下膝关节内侧负载与若干简单参数(如膝关节内收力矩、屈曲力矩和内收角度)之间的显著线性关系，用于膝关节内侧负载的简单估算。本发明，不同行走任务下膝关节内侧负载的估值方法，步骤如下：1)采集直行、踏步转弯和交叉转弯任务下下肢的运动学数据和地反作用力GRF；2)截取右腿支撑期作为分析的时间段；3)通过反向动力学可以由这段数据计算出右腿支撑期的膝关节内收力矩、屈曲力矩和内收角度；4)计算这段时间内的膝关节内侧负载；5)求出右腿支撑期膝关节内收力矩的最大值KAM以及出现时刻t；确定膝关节内侧负载估值公式。本发明主要应用于生物力学测量。
Description
技术领域
本发明涉及医疗器械、生物力学测量，具体讲，涉及一种针对不同行走任务的膝关节内侧负载估值方法。
背景技术
膝关节内侧负载是膝关节内侧骨关节炎发生和发展的重要影响因素。确定膝关节内侧负载最准确的方法是利用膝关节内部的传感器直接在体测量，但这种方法仅适用于进行膝关节置换的患者，不便于应用推广。为了弥补在体测量方法应用局限的缺点，研究人员通常使用肌骨模型估算肌肉的收缩力，进而无创计算膝关节内部的负载。
尽管肌骨模型可以较为准确地估算膝关节负载，但模型复杂，计算量大，耗时长。因此，一些研究学者尝试利用替代参数简单表征膝关节内测负载的变化。1991年，美国圣卢克斯医疗中心的Schipplein等人基于静定肌肉模型，首次提出膝关节内收力矩是膝关节内侧负载的主要决定参数。为了得到更准确的结果，美国佛罗里达大学的Zhao等人利用相关分析研究了膝关节内收力矩与在体测量的内侧负载之间的关系，确定了二者高度的相关性。虽然利用膝关节内收力矩代替膝关节内侧负载已经被大多数学者普遍接受，但仍有少数学者提出了不同的意见。美国佛罗里达大学的Walker等人研究了不同矫正步态下的膝关节内侧负载和内收力矩、屈曲力矩之间的关系，发现利用内收力矩峰值和屈曲力矩绝对值的峰值能够最为准确地估计内侧负载。日本学者Ogaya和德国学者Trepczynski等人也肯定了屈曲力矩在预测膝关节内侧负载中的作用。另外，加拿大学者Adouni和Shirazi-Adl认为膝关节负载的内外侧分布主要由膝关节内收角度决定。综上所述，膝关节内收力矩、屈曲力矩和内收角度都可以表征膝关节内侧负载的变化。
发明内容
为克服现有技术的不足，本发明旨在提出一种膝关节内侧负载估值方法。借助于反向动力学计算、肌骨模型和多元线性回归分析，确定不同不行任务下膝关节内侧负载与若干简单参数(如膝关节内收力矩、屈曲力矩和内收角度)之间的显著线性关系，该量化的线性关系可用于膝关节内侧负载的简单估算。该项发明简化了膝关节内侧负载的测量方法，为膝骨关节炎的临床治疗和康复评价提供了一种简便易行的方法。本发明采用的技术方案是，不同行走任务下膝关节内侧负载的估值方法，步骤如下：
1)利用运动捕捉系统和三维测力板，采集直行、踏步转弯和交叉转弯任务下下肢的运动学数据和地反作用力GRF，运动学数据是指下肢关键位置处标记点的三维位置坐标；
2)截取右腿支撑期作为分析的时间段，即右脚踏在测力板上的时间，也就是GRF竖直分量大于1N的一段时间，这段时间内的运动学数据和GRF为要分析的原始数据；
3)通过反向动力学可以由这段数据计算出右腿支撑期的膝关节内收力矩、屈曲力矩和内收角度；
4)利用肌骨模型计算这段时间内的膝关节内侧负载；
5)求出右腿支撑期膝关节内收力矩的最大值KAM以及出现时刻t，确定t时刻的膝关节屈曲力矩KFM、内收角度KAA和内侧负载KMF；
6)将KMF作为结果变量，KAM，KAM和KFM，KAM、KFM和KAA依次作为预测变量，进行线性回归分析，确定膝关节内侧负载估值公式，如表1：
表1.不同行走任务下膝关节内侧负载估值公式
本发明的特点及有益效果是：
本技术通过反向动力学计算、肌骨模型和线性回归分析，确定了不同行走任务下膝关节内侧负载的估值模型。估值模型的建立可以有效简化膝关节内侧负载的确定，在临床步态分析中，能够快速准确估计不同步态下膝关节内侧负载的变化，实际操作性强，对临床膝关节内侧骨关节炎的康复治疗有重要的指导作用。
附图说明：
图1.膝关节内侧负载估值方法技术流程。
图2.不同直行任务示意图。
图3.标记点位置。
具体实施方式
不同行走任务(直行、踏步转弯和交叉转弯)下的膝关节内侧负载估值方法大体的技术流程如图1所示。
其具体步骤为：
1)利用运动捕捉系统和三维测力板，采集直行、踏步转弯和交叉转弯任务下下肢的运动学数据和地反作用力GRF，运动学数据是指下肢关键位置处标记点的三维位置坐标(如图3)；
2)截取右腿支撑期作为分析的时间段，即右脚踏在测力板上的时间，也就是GRF竖直分量大于1N的一段时间。这段时间内的运动学数据和GRF为要分析的原始数据。
3)通过反向动力学可以由这段数据计算出右腿支撑期的膝关节内收力矩、屈曲力矩和内收角度；
4)利用肌骨模型计算这段时间内的膝关节内侧负载；
5)求出右腿支撑期膝关节内收力矩的最大值KAM以及出现时刻，确定该时刻的膝关节屈曲力矩KFM、内收角度KAA和内侧负载KMF。
6)将KMF作为结果变量，KAM，KAM和KFM，KAM、KFM和KAA依次作为预测变量，进行线性回归分析，确定膝关节内侧负载估值公式(如表1)。由于交叉转弯中，KAA与KMF没有显著的相关性，因此在该任务下，不存在以KAM，KFM和KAA作为预测变量的估值模型。
表1.不同行走任务下膝关节内侧负载估值公式
本发明确定了不同行走任务下膝关节内侧负载的估值模型。本技术的最佳实施方案拟采用专利转让、技术合作或产品开发。
Claims (1)
1.一种不同行走任务下膝关节内侧负载的估值方法，其特征是，步骤如下：
1)利用运动捕捉系统和三维测力板，采集直行、踏步转弯和交叉转弯任务下下肢的运动学数据和地反作用力GRF，运动学数据是指下肢关键位置处标记点的三维位置坐标；
2)截取右腿支撑期作为分析的时间段，即右脚踏在测力板上的时间，也就是GRF竖直分量大于1N的一段时间，这段时间内的运动学数据和GRF为要分析的原始数据；
3)通过反向动力学由所述这段时间内的运动学数据和GRF为要分析的原始数据计算出右腿支撑期的膝关节内收力矩、屈曲力矩和内收角度；
4)利用肌骨模型计算这段时间内的膝关节内侧负载；
5)求出右腿支撑期膝关节内收力矩的最大值KAM以及出现时刻t，确定t时刻的膝关节屈曲力矩KFM、内收角度KAA和内侧负载KMF；
6)将KMF作为结果变量，KAM，KAM和KFM，KAM、KFM和KAA依次作为预测变量，进行线性回归分析，确定膝关节内侧负载估值公式：
对于预测变量KAM，直行任务时，估值公式KMF＝0.704+36.622×KAM，决定系数为0.736，踏步转弯任务时，估值公式KMF＝0.560+41.054×KAM，决定系数为0.744，交叉转弯时，估值公式KMF＝1.568+19.527×KAM，决定系数为0.171；
对于预测变量KAM、KFM，直行任务时，估值公式KMF＝0.473+36.778×KAM+9.000×KFM，决定系数为0.795，踏步转弯任务时，估值公式KMF＝0.448+38.206×KAM+8.467×KFM，决定系数为0.793，交叉转弯时，估值公式KMF＝1.053+39.481×KAM-19.049×KFM，决定系数为0.519；
对于预测变量KAM、KFM、KAA，直行任务时，估值公式KMF＝0.387+38.355×KAM+9.345×KFM+0.036×KAA，决定系数为0.809，踏步转弯任务时，估值公式KMF＝0.371+38.178×KAM+10.098×KFM+0.037×KAA，决定系数为0.806。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711054022.6A CN107736890B (zh) | 2017-10-31 | 2017-10-31 | 不同行走任务下膝关节内侧负载的估值方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711054022.6A CN107736890B (zh) | 2017-10-31 | 2017-10-31 | 不同行走任务下膝关节内侧负载的估值方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107736890A CN107736890A (zh) | 2018-02-27 |
| CN107736890B true CN107736890B (zh) | 2020-08-28 |
Family
ID=61233516
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711054022.6A Active CN107736890B (zh) | 2017-10-31 | 2017-10-31 | 不同行走任务下膝关节内侧负载的估值方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107736890B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109793500B (zh) * | 2019-01-24 | 2021-10-01 | 河南省人民医院 | 膝关节负荷力学分析装置 |
| US11645565B2 (en) * | 2019-11-12 | 2023-05-09 | Optum Services (Ireland) Limited | Predictive data analysis with cross-temporal probabilistic updates |
| CN111728616B (zh) * | 2020-07-01 | 2021-07-20 | 北京航空航天大学 | 一种人体膝关节载荷谱确定方法及系统 |
| CN113100789B (zh) * | 2021-04-16 | 2022-10-21 | 西北工业大学 | 一种膝关节内外侧受力实时分析系统 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009105665A1 (en) * | 2008-02-20 | 2009-08-27 | Mako Surgical Corp. | Implant planning using corrected captured joint motion information |
| JP2012090735A (ja) * | 2010-10-26 | 2012-05-17 | Toyota Motor Corp | 膝関節抵抗計測装置 |
| US9521965B2 (en) * | 2011-08-09 | 2016-12-20 | The Regents Of The University Of California | Mechanical pivot shift measurement to evaluate joints |
| CN102521508B (zh) * | 2011-12-08 | 2014-12-24 | 天津大学 | 一种功能性电刺激下的自适应神经模糊肌肉建模方法 |
| KR102146363B1 (ko) * | 2013-10-31 | 2020-08-20 | 삼성전자주식회사 | 착용형 로봇 및 그 제어 방법 |
| CN103761392B (zh) * | 2014-01-23 | 2017-02-15 | 南京工程学院 | 类人机器人协同运动的肌肉力模型优化方法 |
| CN104921851B (zh) * | 2015-05-25 | 2016-09-07 | 河北工业大学 | 主动型膝上假肢膝关节的预测控制方法 |
| CN105852866B (zh) * | 2016-04-01 | 2018-08-21 | 浙江大学 | 用于膝关节内收力矩测量的可穿戴传感系统及测量方法 |
| CN106326643B (zh) * | 2016-08-16 | 2020-01-14 | 完美鲲鹏(北京)动漫科技有限公司 | 一种应用于类人形骨架的反向动力学解算方法及装置 |
-
2017
- 2017-10-31 CN CN201711054022.6A patent/CN107736890B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107736890A (zh) | 2018-02-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107736890B (zh) | 不同行走任务下膝关节内侧负载的估值方法 | |
| Heintz et al. | Static optimization of muscle forces during gait in comparison to EMG-to-force processing approach | |
| Stief et al. | Reliability and accuracy in three-dimensional gait analysis: a comparison of two lower body protocols | |
| CN107920782B (zh) | 步行分析方法及步行分析系统 | |
| Stagni et al. | Effects of hip joint centre mislocation on gait analysis results | |
| Leardini et al. | A new anatomically based protocol for gait analysis in children | |
| Sartori et al. | Modeling the human knee for assistive technologies | |
| Buchanan et al. | Estimation of muscle forces and joint moments using a forward-inverse dynamics model | |
| Siston et al. | Evaluation of a new algorithm to determine the hip joint center | |
| Wang et al. | Real-time estimation of knee adduction moment for gait retraining in patients with knee osteoarthritis | |
| Liu et al. | Lower extremity joint torque predicted by using artificial neural network during vertical jump | |
| de Rosario et al. | Kinematic description of soft tissue artifacts: quantifying rigid versus deformation components and their relation with bone motion | |
| Larsen et al. | Estimation of spinal loading during manual materials handling using inertial motion capture | |
| Lugrís et al. | Solution methods for the double-support indeterminacy in human gait | |
| Ackermann et al. | Dynamic analysis of human gait disorder and metabolical cost estimation | |
| Duerinck et al. | Forefoot deformation during stance: does the forefoot collapse during loading? | |
| Son et al. | An EMG-based muscle force monitoring system | |
| Wang et al. | A wearable real-time kinetic measurement sensor setup for human locomotion | |
| Dinovitzer et al. | Accurate real-time joint torque estimation for dynamic prediction of human locomotion | |
| Wolf et al. | Day-to-day consistency of lower extremity kinematics during walking and running | |
| Yuen et al. | A comparison of gastrocnemius muscle–tendon unit length during gait using anatomic, cadaveric and MRI models | |
| Ryu et al. | Soft tissue artifact compensation using displacement dependency between anatomical landmarks and skin markers–a preliminary study | |
| Vogt et al. | Overground and walkway ambulation with unilateral hip osteoarthritis: comparison of step length asymmetries and reproducibility of treadmill mounted force plate readings | |
| Sartori et al. | Electromyography‐driven modeling for simulating subject‐specific movement at the neuromusculoskeletal level | |
| JP6079585B2 (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 |