Summary

动态运动过程中骨应变的学科特异性肌肉模型研究

Published: April 11, 2018
doi:

Summary

在着陆期间, 下半身的骨骼经历了巨大的机械负荷, 并且变形了。测量骨的变形以更好地了解与撞击有关的骨应激损伤的机制是非常必要的。一种新的方法, 结合学科特定的肌肉骨骼模型和有限元分析是用来测量胫骨应变在动态运动。

Abstract

骨应激损伤在运动和军事训练中普遍存在。在训练过程中, 重复的大地面撞击力可能是原因。重要的是要确定高地冲击力对下肢骨变形的影响, 以更好地了解骨应激损伤的机制。传统的应变测量方法用于研究体内胫骨变形。这种方法与限制, 包括程序的侵袭性, 少数人类的参与, 和有限的应变数据从小骨表面积。本研究旨在介绍一种新的方法来研究胫骨骨应变在高冲击载荷条件下。一个主题特定的肌肉骨骼模型被创建来代表一个健康的男性 (19 岁, 80 公斤, 1800 毫米)。基于 ct 扫描的对象右胫骨, 建立了一种灵活的有限元胫骨模型。进行了实验室运动捕获, 以获得不同高度 (26、39、52厘米) 降落的运动学和地面反应力。采用多体动态计算机模拟方法, 结合弹性胫骨的模态分析, 对下落着陆时胫骨应变进行量化。计算出的胫骨应变数据与先前的体内研究吻合良好。很明显, 这种非侵入性的方法可以应用于研究大型队列的高冲击活动胫骨骨应变, 这将有助于更好地理解胫骨应力性骨折的损伤机制。

Introduction

骨应激损伤, 如应力性骨折, 是严重的过度使用性损伤, 需要长时间的恢复, 并招致重大医疗费用1,2。应力性骨折在运动和军事人口中都很常见。在所有运动相关的损伤中, 应力性骨折占总3的10%。特别是, 田径运动员在 20%4时面临较高的伤害率。士兵也经历了高压力骨折的发病率。例如, 据报道, 美国陆军的6% 受伤率为1 , 在以色列军队的5中报告了31% 的伤害率。在所有报告的应力性骨折中, 胫骨应力骨折是最常见的一个6,7,8

运动和体育训练, 有较高的胫骨应力骨折风险通常与高地面影响 (例如, 跳跃, 着陆和切割)。在运动过程中, 当脚接触地面时, 地面撞击力会被应用到身体上。这种撞击力被肌肉骨骼系统和鞋类所消散。骨骼系统作为一系列杠杆, 允许肌肉施加力量来吸收地面撞击9。当腿部肌肉不能充分减少地面撞击时, 下肢骨骼必须吸收残余力。骨结构在这个过程中会经历变形。残余冲击力的重复吸收可能导致骨骼中的 microdamages, 这将积聚并成为应力性骨折。到目前为止, 与外部地面撞击力有关的骨骼反应的信息是有限的。研究胫骨骨在动态运动过程中如何响应高冲击力引入的机械载荷是很重要的。在高撞击活动中检查胫骨骨的变形, 可以更好地了解胫骨应力性骨折的发生机制。

用于测量骨骼变形的常规技术在体内依赖于仪器应变表10,11,12,13,14,15。在骨表面植入应变片需要手术治疗。由于侵入性的性质,在体内研究受志愿者的小样本的限制。此外, 应变计只能监测骨骼表面的一个小区域。本文介绍了一种利用计算机模拟分析骨应变的非侵入性方法1617。这种方法可以将肌肉骨骼模型和计算模拟结合起来, 研究人体运动过程中的骨骼应变。

肌肉骨骼模型由骨骼和骨骼肌肉代表。骨架由骨段组成, 它们是刚性或非变形体。骨骼肌肉被建模为控制器使用渐进积分-导数 (PID) 算法。三期 PID 控制使用估计中的错误来提高输出精度18。从本质上讲, 代表肌肉的 PID 控制器试图通过发展必要的力量来复制身体运动, 从而在一段时间内产生肌肉的长度变化。PID 控制器使用长度/时间曲线中的误差来修改复制运动的力。这个模拟过程创造了一个可行的解决方案, 协调所有的肌肉一起工作, 以移动骨骼和产生身体运动。

肌肉骨骼模型骨架中的一个或多个段可以被建模为柔性体, 以允许测量变形。例如, 胫骨骨可以分解成有限数量的元素, 由数以千计的元素和节点组成。通过有限元分析, 可以研究机械载荷对挠性胫骨的影响。有限元分析计算了各个元素随时间推移的加载响应。随着骨元素和节点数量的增加, 有限元分析的计算时间将显著增加。

为降低计算成本, 并对柔性体的变形进行精确评价, 采用模态有限元分析方法, 在汽车和航天工业中应用了19,20。该过程不是分析单个 FE 元素对时间域内机械载荷的响应, 而是根据频率域中不同的振动频率来评估物体的机械响应。此方法使计算时间大大减少, 同时提供了对变形20的精确测量。虽然模态有限元分析已广泛应用于汽车和航天领域的机械疲劳研究, 但这种方法在人类运动科学中的应用却十分有限。Al Nazer et, 使用一种模态有限元分析, 检查人体步态中的胫骨变形, 并报告令人鼓舞的结果16,17。然而, 只有利用有限的运动学数据进行计算机仿真, 才对其方法产生很大的影响;没有真正的地面撞击力用于辅助模拟。这种方法可能是合理的研究低冲击慢动作, 如步行, 但它不是一个可行的解决方案, 研究高地面撞击运动。因此, 为了检查在动态高影响活动期间的下肢骨反应, 必须开发一种创新的方法来解决与以前报告的方法相关的限制。具体来说, 必须开发一种利用精确的实验运动数据和实际地面冲击力的方法。因此, 本研究的目的是建立一个特定主题的肌肉骨骼模型, 进行多体动力学模拟, 用模态有限元分析, 检查胫骨应变在高冲击活动。选择了由不同高度的降落所表示的动态高冲击运动来测试该方法。

Protocol

这项试验是根据《赫尔辛基宣言》进行的。在收集数据之前, 该主题审查并签署了大学机构审查委员会批准的同意表格, 然后才参加研究。 1. CT 成像协议 把参与者带到一个有 CT 扫描仪的设施。在 ct 扫描之前, 配置 ct 机具有以下参数: ct 切片厚度为0.625 毫米, ct 视场为15厘米 x 15 厘米, 并采用机器算法自动设置峰值电压 (kVp) 和毫安秒 (mAs) 参数。 让参与者躺在一个…

Representative Results

一个健康的高加索男性 (19 岁, 身高1800毫米, 质量80公斤) 自愿为研究。在收集数据之前, 该主题审查并签署了大学机构审查委员会批准的同意表格, 然后才参加研究。这项试验是根据《赫尔辛基宣言》进行的。实验是根据下面的协议进行的。 为了验证正向动态仿真的正确性, 将模拟中的低体关节角与由生物力学分析程序处理的运…

Discussion

本研究的目的是发展一种非侵入性的方法, 以确定胫骨的变形在高影响活动。对胫骨负荷的影响进行量化, 将有助于更好地理解胫骨应力性骨折。在本研究中, 开发了一个特定于主题的肌肉骨骼模型, 并运行计算机模拟, 以复制在实验室环境中进行的降落运动。研究了降落地高度对胫骨应变的影响。在本研究中, 我们观察到, 随着降落高度的增加, 峰值最大的主要菌株也是如此。此外, 在三个着陆条件?…

Divulgazioni

The authors have nothing to disclose.

Acknowledgements

陆军部 #W81XWH-08-1-0587, #W81XWH-15-1-0006;球州立大学2010立志补助金。

Materials

CT Scanner GE Medical System N/A Light Speed VCT. For performing tibia CT scan.
Motion Capture System Vicon Inc N/A Vicon FX40 high speed cameras. For performing 3D motion capture.
Force plates AMTI Inc N/A Collecting 3D ground reaction forces
Vicon Nexus Vicon Inc N/A Motion capture software program. For processing visual marker trajectory data.
Visual 3D C-Motion Inc N/A Biomechanics analysis software. For computing 3D kinematics and kinetics of human movements.
MATLAB Mathworks Inc N/A Computer programming software. For performing raw data filtering, data conversion, and data processing.
ADAMS 2012 MSC Software Inc N/A Multibody dynamic computer simulation program.
LifeMOD Lifemodeler Inc N/A A software Plug-in in ADAMS. For building human body musculo-skeletal models.
MIMICS 13 Materialise Inc N/A Image processing program. A 3D modeling tool to process imaging data. For creating 3D tibia model from CT scans.
MARC 2012 MSC Software Inc N/A Finite element analysis software. For performing volumn meshing, generating tibia FE model, and running modal FE analysis.
SPSS 19 IBM Inc N/A Statistical analysis software.

Riferimenti

  1. Brukner, P., Bennell, K., Matheson, G. . Stress fracture. , (1999).
  2. Zadpoor, A., Nikooyan, A. The relationship between lower-extremity stress fractures and the ground reaction force: A systematic review. Clin Biomech. 26, 23-28 (2011).
  3. Matheson, G. O., Clement, D. B., McKenzie, D. C., Taunton, J. E., Lioyd-Smith, D. R., Maclntyre, J. G. Stress fractures in athletes. A study of 320 cases. Am J Sports Med. 15, 46-58 (1987).
  4. Bennell, K., Grimston, S., Burr, D., Milgrom, C. Risk factors for developing stress fractures. Musculoskeletal fatigue and stress fractures. , 15-33 (2001).
  5. Milgrom, C., Giladi, M., Stein, M., Kashtan, H., Margulies, J. Y., Chisin, R., Stenberg, R., Aharonson, Z. Stress fractures in military recruits. A prospective study showing an unusually high incidence. J Bone Joint Surg Br. 67, 732-735 (1985).
  6. Almeida, S. A., Williams, K. M., Shaffer, R. A., Brodine, S. K. Epidemiological patterns of musculoskeletal injuries and physical training. Med Sci Sports Exerc. 31, 1176-1182 (1999).
  7. Jones, B. H., Knapik, J. J. Physical training and exercise-related injuries, surveillance, research and injury prevention in military populations. Sports Med. 27, 111-125 (1999).
  8. Jones, B. H., Thacker, S., Gilchrist, J., Kimsey, C. D., Sosin, D. M. Prevention of lower extremity stress fractures in athletes and soldiers: a systematic review. Epidemiol Rev. 24, 228-247 (2002).
  9. Voloshin, A., Wosk, J. An in vivo study of low back pain and shock absorption in the human locomotor system. J Biomech. 15, 21-27 (1982).
  10. Burr, D. B., Milgrom, C., Fyhrie, D., Forwood, M., Nyska, M., Finestone, A., Hoshaw, S., Saiag, E., Simkin, A. In vivo measurement of human tibial strains during vigorous activity. Bone. 18, 405-410 (1996).
  11. Ekenman, I., Halvorsen, K., Westblad, P., Fellander-Tsai, L., Rolf, C. The reliability and validity of an instrumented staple system for in vivo measurement of local bone deformation. An in vitro study. Scand J Med Sci Sports. 8, 172-176 (1998).
  12. Lanyon, L. E., Hampson, W. G., Goodship, A. E., Shah, J. S. Bone deformation recorded in vivo from strain gauges attached to the human tibial shaft. Acta Orthop Scand. 46, 256-268 (1975).
  13. Ekenman, I., Halvorsen, K., Westblad, P., Tsai, L. F., Rolf, C. Local bone deformation at two predominant sites for stress fractures of the tibia: an in vivo study. Foot Ankle Int. 19, 479-484 (1998).
  14. Milgrom, C., Finestone, A., Levi, Y., Simkin, A., Ekenman, I., Mendelson, S., Millgram, M., Nyska, M., Benjuya, N., Burr, D. Do high impact exercises produce higher tibial strains than running?. Br J Sports Med. 34, 195-199 (2000).
  15. Milgrom, C., Finestone, A., Simkin, A., Ekenman, I., Mendelson, S., Millgram, M., Nyska, M., Larsson, E., Burr, D. In-vivo strain measurements to evaluate the strengthening potential of exercises on the tibial bone. J Bone Joint Surg Br. 82, 591-594 (2000).
  16. Al Nazer, R., Rantalainen, T., Heinonen, A., Sievanen, H., Mikkola, A. Flexible multibody simulation approach in the analysis of tibial strain during walking. J Biomech. 41, 1036-1043 (2008).
  17. Al Nazer, R., Klodowski, A., Rantalainen, T., Heinonen, A., Sievanen, H., Mikkola, A. A full body musculoskeletal model based on flexible multibody simulation approach utilised in bone strain analysis during human locomotion. Comput Method Biomec. 14, 573-579 (2011).
  18. Johnson, M. A., Moradi, M. H., Crowe, J. . PID control: new identification and design methods. , 543 (2005).
  19. Craig, R. R., Bampton, M. C. C. Coupling of substructures for dynamics analysis. American Institute of Aeronautics and Astronautics Journal. 6, 1313-1319 (1968).
  20. Wasfy, T. M., Noor, A. K. Computational strategies for flexible multibody systems. Appl Mech Rev. 56, 553-613 (2003).
  21. Kadaba, M. P., Ramakrishnan, H. k., Wootten, M. E. Measurement of lower extremity kinematics during level walking. J Orthop Res. 8, 383-392 (1990).
  22. Schwartz, M. H., Rozumalski, A. A new method for estimating joint parameters from motion data. J Biomech. 38, 107-116 (2005).
  23. Devita, P., Skelly, W. A. Effect of landing stiffness on joint kenetics and energetic in the lower extremity. Med Sci Sports Exerc. 24, 108-115 (1992).
  24. Dong, X. N., Guo, X. E. The dependence of transversely isotropic elasticity of human femoral cortical bone on porosity. J Biomech. 37, 1281-1287 (2004).
  25. Schileo, E., Taddei, F., Malandrino, A., Cristofolini, L., Viceconti, M. Subject-specific finite element models can accurately predict strain levels in long bones. J Biomech. 40, 2982-2989 (2007).
  26. Pattin, C. A., Caler, W. E., Carter, D. R. Cyclic mechanical property degradation during fatigue loading of cortical bone. J Biomech. 29, 69-79 (1996).
  27. Lifemodeler, I. . Lifemod Manual. , (2010).
check_url/it/56759?article_type=t

Play Video

Citazione di questo articolo
Wang, H., Dueball, S. Subject-specific Musculoskeletal Model for Studying Bone Strain During Dynamic Motion. J. Vis. Exp. (134), e56759, doi:10.3791/56759 (2018).

View Video