Summary

Subject-specific Musculoskeletal Model for Studying Bone Strain During Dynamic Motion

Published: April 11, 2018
doi:

Summary

During landing, lower-body bones experience large mechanical loads and are deformed. It is essential to measure bone deformation to better understand the mechanisms of bone stress injuries associated with impacts. A novel approach integrating subject-specific musculoskeletal modeling and finite element analysis is used to measure tibial strain during dynamic movements.

Abstract

Bone stress injuries are common in sports and military trainings. Repetitive large ground impact forces during training could be the cause. It is essential to determine the effect of high ground impact forces on lower-body bone deformation to better understand the mechanisms of bone stress injuries. Conventional strain gauge measurement has been used to study in vivo tibia deformation. This method is associated with limitations including invasiveness of the procedure, involvement of few human subjects, and limited strain data from small bone surface areas. The current study intends to introduce a novel approach to study tibia bone strain under high impact loading conditions. A subject-specific musculoskeletal model was created to represent a healthy male (19 years, 80 kg, 1,800 mm). A flexible finite element tibia model was created based on a computed tomography (CT) scan of the subject's right tibia. Laboratory motion capture was performed to obtain kinematics and ground reaction forces of drop-landings from different heights (26, 39, 52 cm). Multibody dynamic computer simulations combined with a modal analysis of the flexible tibia were performed to quantify tibia strain during drop-landings. Calculated tibia strain data were in good agreement with previous in vivo studies. It is evident that this non-invasive approach can be applied to study tibia bone strain during high impact activities for a large cohort, which will lead to a better understanding of injury mechanism of tibia stress fractures.

Introduction

Bone stress injuries, such as stress fractures, are severe overuse injuries requiring long periods of recovery and incurring significant medical cost1,2. Stress fractures are common both in athletic and military populations. Among all sports related injuries, stress fractures account for 10% of the total3. In particular, track athletes face a higher injury rate at 20%4. Soldiers also experience a high rate of stress fractures. For instance, a 6% injury rate was reported for the US Army1 and a 31% injury rate was reported in the Israeli Army5. Among all reported stress fractures, tibia stress fracture is the most common one6,7,8.

Sports and physical trainings with a higher risk of tibia stress fracture are normally associated with high ground impacts (e.g., jumping, landing, and cutting). During locomotion, a ground impact force is applied to the body when the foot contacts the ground. This impact force is dissipated by the musculoskeletal system and footwear. The skeletal system serves as a series of levers allowing muscles to apply forces to absorb the ground impact9. When the leg muscles cannot adequately reduce the ground impact, lower-body bones must absorb the residual force. Bone structure will experience deformation during this process. Repetitive absorption of residual impact force may result in microdamages in the bone, which will accumulate and become stress fractures. To date, information related to bone reaction to external ground impact forces is limited. It is important to study how the tibia bone responds to the mechanical load introduced by high impact forces during dynamic motions. Examining tibia bone deformation during high impact activities could lead to a better understanding of the mechanism of tibia stress fracture.

Conventional techniques used to measure bone deformation in vivo rely on instrumented strain gauges10,11,12,13,14,15. Surgical procedures are needed to implant strain gauges on bone surface. Due to the invasive nature, in vivo studies are limited by a small sample of volunteers. In addition, the strain gauge can only monitor a small region of the bone surface. Recently, a non-invasive method utilizing computer simulation to analyze bone strain was introduced16,17. This methodology allows for the ability to combine musculoskeletal modeling and computational simulations to study bone strain during human movement.

A musculoskeletal model is represented by a skeleton and skeletal muscles. The skeleton consists of bone segments, which are rigid or non-deformable bodies. Skeletal muscles are modeled as controllers using the progressive-integral-derivative (PID) algorithm. The three-term PID control uses errors in estimation to improve the output accuracy18. In essence, PID controllers representing muscles try to duplicate body movements by developing necessary forces to produce length changes of the muscles over time. The PID controller uses the error in the length/time curve to modify the force for reproducing the movement. This simulation process creates a feasible solution to coordinate all muscles to work together to move the skeleton and produce the body movement.

One or more segments in the skeleton of the musculoskeletal model can be modeled as flexible bodies to allow measurement of deformation. For instance, the tibia bone can be broken down into a finite number of elements, which consists of thousands of elements and nodes. The effect of mechanical loading on the flexible tibia can be examined through finite element (FE) analysis. The FE analysis calculates the loading response of individual elements over time. As the number of bone elements and nodes increase, the computation time of the FE analysis will significantly increase.

To reduce computational cost with accurate evaluation of flexible bodies' deformation, modal FE analysis has been developed and used within the automotive and aerospace industry19,20. Instead of analyzing individual FE elements' responses to mechanical load in the time domain, this procedure assesses an object's mechanical responses based on different vibrational frequencies in the frequency domain. This method results in a significant reduction in computation time while providing accurate measurement of deformation20. Although modal FE analysis has been widely used to study mechanical fatigue in automotive and aerospace areas, the application of this method has been very limited in human movement science. Al Nazer et al., used a modal FE analysis to examine tibial deformation during human gait and reported encouraging results16,17. However, their method was greatly affected by only using limited kinematic data from an experiment to drive the computer simulations; There were no real ground impact forces used to assist the simulations. This approach may be reasonable for studying low impact slow motions such as walking, but it is not a feasible solution to study high ground impact movements. Thus, in order to examine lower-body bone reactions during dynamic high impact activities, it is essential to develop an innovative approach to address the limitations associated with the previously reported method. Specifically, a method utilizing accurate experimental kinematic data and real ground impact forces must be developed. Therefore, the goal of this study was to develop a subject-specific musculoskeletal model to perform multibody dynamic simulations with modal FE analysis to examine tibial strain during high impact activities. A dynamic high impact movement represented by drop-landings from different heights was selected to test the method.

Protocol

The experiment was conducted under the Helsinki Declaration. Prior to data collection, the subject reviewed and signed the consent form approved by the University Institutional Review Board before participating in the study. 1. CT Imaging Protocol Take the participant to a facility where a CT scanner is housed. Prior to the CT scan, configure the CT machine with the following parameters: CT slice thickness of 0.625 mm, CT field of view of 15 cm x 15 cm, and auto setting for paramet…

Representative Results

A healthy Caucasian male (19 years, height 1,800 mm, mass 80 kg) volunteered for the study. Prior to data collection, the subject reviewed and signed the consent form approved by the University Institutional Review Board before participating in the study. The experiment was conducted under the Helsinki Declaration. The experiment was performed based on the following protocol. In order to verify the accuracy of the forward dynami…

Discussion

The purpose of this study was to develop a non-invasive method to determine tibia deformation during high impact activities. Quantifying tibia strain due to impact loading will lead to a better understanding of tibia stress fracture. In this study, a subject-specific musculoskeletal model was developed, and computer simulations were run to duplicate the drop-landing movements performed in a laboratory setting. The effect of drop-landing height on tibial strain was examined. In this study, we observed that as the drop-lan…

Disclosures

The authors have nothing to disclose.

Acknowledgements

Department of the Army #W81XWH-08-1-0587, #W81XWH-15-1-0006; Ball State University 2010 ASPiRE grant.

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.

References

  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).

Play Video

Cite This Article
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