Compliant Vs. rigid tendon models: A simulation study on precision, computational efficiency and numerical stability

Providing an efficient mathematical model of the skeletal muscles which takes both computational efficiency and accuracy into account is a crucial factor in the simulations of multiple-muscle problems. Previous studies stated that ignoring the elastic characteristics of the tendon can reduce the time cost of simulations at the expense of introducing some minor errors if the ratio of tendon slack length to muscle optimum length is less than or equal to unity. The purpose of this paper was to test the precision, efficiency and numerical stability of this criterion for the muscles of the human body in their usual length excursions. In this regard two muscles of the upper extremity (Brachioradialis (BRD) long head of biceps (BICL)) and one from the lower extremity (soleus (SOL)) of the human body have been chosen. Two variations of a general Hill-based musculotendon model have been considered in this study. In the first one, using a nonlinear spring the elastic properties of the tendon has been incorporated into the model and in the second one, ignoring this properties, a constant length for the tendon has been assumed. The mean absolute error between the force profiles of the two models for BRD, BICL and SOL were 4.2, 12 and 13.1 respectively. Also rigid-tendon model was 7.3 to 9.5 times faster than compliant-tendon model using the implicit integrator. For BRD the outcomes of the two models, have similar trends and the discrepancies between the force profiles are negligible. However, the results obtained from the compliant-tendon model illustrate some numerical stability problems. In the second muscle, i.e. BICL, likewise BRD the trends of the force profiles are the same; however, the disparity among the outcomes of the two models have escalated. Likewise BRD, the rigid-tendon model required less computational time. Inspecting the results obtained for SOL, one can easily spot the significant differences between the outcomes of the two models. Considering the tendon slack length to the optimum muscle length ratio for the three mentioned musculotendon units, one can draw this conclusion that, in case this value is less than unity using the rigid-tendon model is recommended. If this value is not much greater than unity, like BICL, exploiting the rigid-tendon model will increase the computational efficiency in expense of contaminating the outcomes with some amounts of error. However, if this ratio is far from unity, like SOL, ignoring the length alterations in the tendon is not recommended.

[1]  Ayman Habib,et al.  OpenSim: Open-Source Software to Create and Analyze Dynamic Simulations of Movement , 2007, IEEE Transactions on Biomedical Engineering.

[2]  M. Pandy,et al.  Dynamic optimization of human walking. , 2001, Journal of biomechanical engineering.

[3]  John H Challis,et al.  The influence of an elastic tendon on the force producing capabilities of a muscle during dynamic movements , 2007, Computer methods in biomechanics and biomedical engineering.

[4]  Germund Dahlquist,et al.  Numerical methods in scientific computing , 2008 .

[5]  F. Zajac Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. , 1989, Critical reviews in biomedical engineering.

[6]  David G Lloyd,et al.  Neuromusculoskeletal modeling: estimation of muscle forces and joint moments and movements from measurements of neural command. , 2004, Journal of applied biomechanics.

[7]  Matthew Millard,et al.  Flexing computational muscle: modeling and simulation of musculotendon dynamics. , 2013, Journal of biomechanical engineering.

[8]  Joseph Hamill,et al.  Evaluation of the minimum energy hypothesis and other potential optimality criteria for human running , 2012, Proceedings of the Royal Society B: Biological Sciences.

[9]  I. Matsubara,et al.  X-ray diffraction studies on skinned single fibres of frog skeletal muscle. , 1972, Journal of molecular biology.

[10]  S L Delp,et al.  A graphics-based software system to develop and analyze models of musculoskeletal structures. , 1995, Computers in biology and medicine.

[11]  Jeremy D Wong,et al.  The central nervous system does not minimize energy cost in arm movements. , 2010, Journal of neurophysiology.