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. ; ; (, 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
  2. ; ; ; ; (, 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
  3. ; ; ; ; (, 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
  4. ; (, Integrative And Comparative Biology)
    Synopsis Movement is integral to animal life, and most animal movement is actuated by the same engine: striated muscle. Muscle input is typically mediated by skeletal elements, resulting in musculoskeletal systems that are geared: at any instant, the muscle force and velocity are related to the output force and velocity only via a proportionality constant G, the “mechanical advantage”. The functional analysis of such “simple machines” has traditionally centered around this instantaneous interpretation, such that a small vs large G is thought to reflect a fast vs forceful system, respectively. But evidence is mounting that a comprehensive analysis ought to also consider the mechanical energy output of a complete contraction. Here, we approach this task systematically, and deploy the theory of physiological similarity to study how gearing affects the flow of mechanical energy in a minimalist model of a musculoskeletal system. Gearing influences the flow of mechanical energy in two key ways: it can curtail muscle work output, because it determines the ratio between the characteristic muscle kinetic energy and work capacity; and it defines how each unit of muscle work is partitioned into different system energies, that is, into kinetic vs “parasitic” energy such as heat. As a consequence of both effects, delivering maximum work in minimum time and with maximum output speed generally requires a mechanical advantage of intermediate magnitude. This optimality condition can be expressed in terms of two dimensionless numbers that reflect the key geometric, physiological, and physical properties of the interrogated musculoskeletal system, and the environment in which the contraction takes place. Illustrative application to exemplar musculoskeletal systems predicts plausible mechanical advantages in disparate biomechanical scenarios, yields a speculative explanation for why gearing is typically used to attenuate the instantaneous force output ($$G_{\text{opt}} \lt 1)$$, and predicts how G needs to vary systematically with animal size to optimize the delivery of mechanical energy, in superficial agreement with empirical observations. A many-to-one mapping from musculoskeletal geometry to mechanical performance is identified, such that differences in G alone do not provide a reliable indicator for specialization for force vs speed—neither instantaneously, nor in terms of mechanical energy output. The energy framework presented here can be used to estimate an optimal mechanical advantage across variable muscle physiology, anatomy, mechanical environment, and animal size, and so facilitates investigation of the extent to which selection has made efficient use of gearing as a degree of freedom in musculoskeletal “design.” 
    more » « less
  5. ; ; ; ; ; ; ; ; (, Biological Cybernetics)
    Abstract The coordination of complex behavior requires knowledge of both neural dynamics and the mechanics of the periphery. The feeding system ofAplysia californicais an excellent model for investigating questions in soft body systems’ neuromechanics because of its experimental tractability. Prior work has attempted to elucidate the mechanical properties of the periphery by using a Hill-type muscle model to characterize the force generation capabilities of the key protractor muscle responsible for movingAplysia’s grasper anteriorly, the I2 muscle. However, the I1/I3 muscle, which is the main driver of retractions ofAplysia’s grasper, has not been characterized. Because of the importance of the musculature’s properties in generating functional behavior, understanding the properties of muscles like the I1/I3 complex may help to create more realistic simulations of the feeding behavior ofAplysia, which can aid in greater understanding of the neuromechanics of soft-bodied systems. To bridge this gap, in this work, the I1/I3 muscle complex was characterized using force-frequency, length-tension, and force-velocity experiments and showed that a Hill-type model can accurately predict its force-generation properties. Furthermore, the muscle’s peak isometric force and stiffness were found to exceed those of the I2 muscle, and these results were analyzed in the context of prior studies on the I1/I3 complex’s kinematics in vivo. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email