skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Numerical instability of Hill-type muscle models
Hill-type muscle models are highly preferred as phenomenological models for musculoskeletal simulation studies despite their introduction almost a century ago. The use of simple Hill-type models in simulations, instead of more recent cross-bridge models, is well justified since computationally ‘light-weight’—although less accurate—Hill-type models have great value for large-scale simulations. However, this article aims to invite discussion on numerical instability issues of Hill-type muscle models in simulation studies, which can lead to computational failures and, therefore, cannot be simply dismissed as an inevitable but acceptable consequence of simplification. We will first revisit the basic premises and assumptions on the force–length and force–velocity relationships that Hill-type models are based upon, and their often overlooked but major theoretical limitations. We will then use several simple conceptual simulation studies to discuss how these numerical instability issues can manifest as practical computational problems. Lastly, we will review how such numerical instability issues are dealt with, mostly in an ad hoc fashion, in two main areas of application: musculoskeletal biomechanics and computer animation.  more » « less
Award ID(s):
1846368
PAR ID:
10419461
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of The Royal Society Interface
Volume:
20
Issue:
199
ISSN:
1742-5662
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; (, bioRxiv)
    Abstract Musculoskeletal simulations can offer valuable insight into how the properties of our musculoskeletal system influence the biomechanics of our daily movements. One such property is muscle’s history-dependent initial resistance to stretch, also known as short-range stiffness, which is key to stabilizing movements in response to external perturbations. Short-range stiffness is poorly captured by existing musculoskeletal simulations since they employ phenomenological Hill-type muscle models that lack the mechanisms underlying short-range stiffness. While it has been previously shown that biophysical cross-bridge models can reproduce muscle short-range-stiffness, it is unclear which specific biophysical properties are necessary to capture history-dependent muscle force responses in behaviorally relevant conditions. Here, we tested the ability of various biophysical cross-bridge models to reproduce empirical short-range stiffness and its history-dependent changes across a broad range of behaviorally relevant length changes and activation levels, using an existing dataset on permeabilized rat soleus muscle fibers (N = 11). We found that a biophysical cross-bridge model with cooperative myofilament activation reproduced the effects of muscle activation (R2= 0.86), stretch amplitude (R2= 0.71) and isometric recovery time (R2= 0.79) on history-dependent changes in short-range stiffness after shortening. Similar results were obtained when the cross-bridge distribution of the biophysical model was approximated by a Gaussian (R2= 0.73 - 0.88), but at a 20 times lower computational cost. These effects could not be reproduced by either a biophysical cross-bridge model without cooperative myofilament activation or a Hill-type model (R2< 0.5). The reduced computational demand of the Gaussian-approximated models facilitates implementing biophysical cross-bridge models with cooperative myofilament activation in musculoskeletal simulations to improve the prediction of short-range stiffness during movements. 
    more » « less
  2. ; ; (, American Journal of Biological Anthropology)
    Abstract ObjectivesMusculoskeletal modeling is a powerful approach for studying the biomechanics and energetics of locomotion.Australopithecus (A.) afarensisis among the best represented fossil hominins and provides critical information about the evolution of musculoskeletal design and locomotion in the hominin lineage. Here, we develop and evaluate a three‐dimensional (3‐D) musculoskeletal model of the pelvis and lower limb ofA. afarensisfor predicting muscle‐tendon moment arms and moment‐generating capacities across lower limb joint positions encompassing a range of locomotor behaviors. Materials and MethodsA 3‐D musculoskeletal model of an adultA. afarensispelvis and lower limb was developed based primarily on the A.L. 288‐1 partial skeleton. The model includes geometric representations of bones, joints and 35 muscle‐tendon units represented using 43 Hill‐type muscle models. Two muscle parameter datasets were created from human and chimpanzee sources. 3‐D muscle‐tendon moment arms and isometric joint moments were predicted over a wide range of joint positions. ResultsPredicted muscle‐tendon moment arms generally agreed with skeletal metrics, and corresponded with human and chimpanzee models. Human and chimpanzee‐based muscle parameterizations were similar, with some differences in maximum isometric force‐producing capabilities. The model is amenable to size scaling from A.L. 288‐1 to the larger KSD‐VP‐1/1, which subsumes a wide range of size variation inA. afarensis. DiscussionThis model represents an important tool for studying the integrated function of the neuromusculoskeletal systems inA. afarensis. It is similar to current human and chimpanzee models in musculoskeletal detail, and will permit direct, comparative 3‐D simulation studies. 
    more » « less
  3. ; ; (, Journal of Experimental Biology)
    ABSTRACT Hill-type muscle models are widely used, even though they do not accurately represent the relationship between activation and force in dynamic contractions. We explored the use of neural networks as an alternative approach to capture features of dynamic muscle function, without a priori assumptions about force–length–velocity relationships. We trained neural networks using an existing dataset of two guinea fowl muscles to estimate muscle force from activation, fascicle length and velocity. Training data were recorded using sonomicrometry, electromyography and a tendon buckle. First, we compared the neural networks with Hill-type muscle models, using the same data for network training and model optimization. Second, we trained neural networks on larger datasets, in a more realistic machine learning scenario. We found that neural networks generally yielded higher coefficients of determination and lower errors than Hill-type muscle models. Neural networks performed better when estimating forces on the muscle used for training, but on another bird, than on a different muscle of the same bird, likely due to inaccuracies in activation and force scaling. We extracted force–length and force–velocity relationships from the trained neural networks and found that both effects were underestimated and the relationships were not well replicated outside the training data distribution. We discuss suggested experimental designs and the challenge of collecting suitable training data. Given a suitable training dataset, neural networks could provide a useful alternative to Hill-type muscle models, particularly for modeling muscle dynamics in faster movements; however, scaling of the training data should be comparable between muscles and animals. 
    more » « less
  4. ; ; ; ; (, BioMed)
    Background/Objectives: This study explores an optimization-based strategy for muscle force estimation by employing simplified cost functions integrated with physiologically relevant muscle models. Methods: Considering elbow flexion as a case study, we employ an inverse-dynamics approach to estimate muscle forces for the biceps brachii, brachialis, and brachioradialis, utilizing different combinations of cost functions and muscle constitutive models. Muscle force generation is modeled by accounting for active and passive contractile behavior to varying degrees using Hill-type models. In total, three separate cost functions (minimization of total muscle force, mechanical work, and muscle stress) are evaluated with each muscle force model to represent potential neuromuscular control strategies without relying on electromyography (EMG) data, thereby characterizing the interplay between muscle models and cost functions. Results: Among the evaluated models, the Hill-type muscle model that incorporates both active and passive properties, combined with the stress minimization cost function, provided the most accurate predictions of muscle activation and force production for all three arm flexor muscles. Our results, validated against existing biomechanical data, demonstrate that even simplified cost functions, when paired with detailed muscle models, can achieve high accuracy in predicting muscle forces. Conclusions: This approach offers a versatile, EMG-free alternative for estimating muscle recruitment and force production, providing a more accessible and adaptable tool for muscle force analysis. It has profound implications for enhancing rehabilitation protocols and athletic training, not only broadening the applicability of muscle force estimation in clinical and sports settings but also paving the way for future innovations in biomechanical research. 
    more » « less
  5. ; ; ; ; (, ACM Transactions on Graphics)
    We propose a simple and practical approach for incorporating the effects of muscle inertia, which has been ignored by previous musculoskeletal simulators in both graphics and biomechanics. We approximate the inertia of the muscle by assuming that muscle mass is distributed along the centerline of the muscle. We express the motion of the musculotendons in terms of the motion of the skeletal joints using a chain of Jacobians, so that at the top level, only the reduced degrees of freedom of the skeleton are used to completely drive both bones and musculotendons. Our approach can handle all commonly used musculotendon path types, including those with multiple path points and wrapping surfaces. For muscle paths involving wrapping surfaces, we use neural networks to model the Jacobians, trained using existing wrapping surface libraries, which allows us to effectively handle the Jacobian discontinuities that occur when musculotendon paths collide with wrapping surfaces. We demonstrate support for higher-order time integrators, complex joints, inverse dynamics, Hill-type muscle models, and differentiability. In the limit, as the muscle mass is reduced to zero, our approach gracefully degrades to traditional simulators without support for muscle inertia. Finally, it is possible to mix and match inertial and non-inertial musculotendons, depending on the application. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email