A Monte Carlo analysis of muscle force estimation sensitivity to muscle-tendon properties using a Hill-based muscle model

• Published in: Journal of biomechanics Vol. 79; pp. 67 - 77
• Main Authors: Bujalski, P., Martins, J., Stirling, L.
• Format: Journal Article
• Language: English
• Published: United States Elsevier Ltd 05.10.2018; Elsevier Limited
• Subjects: Computer simulation; Confidence intervals; Contraction; Curvature; Electromyography; Force; Humans; Isometric Contraction - physiology; Mathematical models; Models, Biological; Monte Carlo Method; Monte Carlo simulation; Morphology; Muscle contraction; Muscle Contraction - physiology; Muscle, Skeletal - physiology; Muscles; Muscles - physiology; Parameter sensitivity; Parameter uncertainty; Physiology; Response surface methodology; Running - physiology; Sensitivity analysis; Shape factor; Tendons - physiology; Velocity
• ISSN: 0021-9290; 1873-2380; 1873-2380
• DOI: 10.1016/j.jbiomech.2018.07.045

Surface electromyography driven models are desirable for estimating subject-specific muscle forces. However, these models include parameters that come from an array of sources, thus creating uncertainty in the model-estimated force. In this study, we used Monte-Carlo simulations to evaluate the sensitivity of Hill-based model muscle forces to changes in 11 parameters in the muscle-tendon unit morphological properties and in the model force-length and force-velocity relationships. We decomposed the force variability and ranked the sensitivity of the model to the underlying parameters using the Variogram Analysis of Response Surfaces. For the analyzed running experiments and the adopted Hill model structure, our results show that the parameters are separable into four groups, where the parameters in each group have a synergistic contribution to the model global sensitivity. The first group consists of the maximum isometric force and the pennation angle. The second group contains the optimal fiber length, the tendon slack length, the tendon reference strain and the tendon shape factor. The third group contains the width and shape at the extremities of the active contractile element, along with the maximum contraction velocity and the curvature constant in the force-velocity curve. The fourth group consisted only of the force enhancement during eccentric contraction. The first two groups revealed the largest influence on the output force sensitivity. As many input parameters are difficult to measure and impact estimated forces, we propose that model estimates be presented with confidence intervals as well as inter-parameter relationships, to encourage users to explicitly consider the model uncertainty.