Effects of kinematic complexity and number of muscles on musculoskeletal model robustness to muscle dysfunction

The role of the deep cervical extensor muscles in multi-directional isometric neck strength

Clinical management of whiplash-associated disorders is challenging and often unsuccessful, with over a third of whiplash injuries progressing to chronic neck pain. Previous imaging studies have identified muscle fat infiltration, indicative of muscle weakness, in the deep cervical extensor muscles (multifidus and semispinalis cervicis). Yet, kinematic and muscle redundancy prevent the direct assessment of individual neck muscle strength, making it difficult to determine the role of these muscles in motor dysfunction. The purpose of this study was to determine the effects of deep cervical extensor muscle weakness on multi-directional neck strength and muscle activation patterns. Maximum isometric forces and associated muscle activation patterns were computed in 25 test directions using a 3-joint, 24-muscle musculoskeletal model of the head and neck. The computational approach accounts for differential torques about the upper and lower cervical spine. To facilitate clinical translation, the test directions were selected based on locations where resistance could realistically be applied to the head during clinical strength assessments. Simulation results reveal that the deep cervical extensor muscles are active and contribute to neck strength in directions with an extension component. Weakness of this muscle group leads to complex compensatory muscle activation patterns characterized primarily by increased activation of the superficial extensors and deep upper cervical flexors, and decreased activation of the deep upper cervical extensors. These results provide a biomechanistic explanation for movement dysfunction that can be used to develop targeted diagnostics and treatments for chronic neck pain in whiplash-associated disorders.

Comparison of strength profile representations using musculoskeletal models and their applications in robotics

Musculoskeletal models provide an approach towards simulating the ability of the human body in a variety of human-robot applications. A promising use for musculoskeletal models is to model the physical capabilities of the human body, for example, estimating the strength at the hand. Several methods of modelling and representing human strength with musculoskeletal models have been used in ergonomic analysis, human-robot interaction and robotic assistance. However, it is currently unclear which methods best suit modelling and representing limb strength. This paper compares existing methods for calculating and representing the strength of the upper limb using musculoskeletal models. It then details the differences and relative advantages of the existing methods, enabling the discussion on the appropriateness of each method for particular applications.

Muscle redundancy is greatly reduced by the spatiotemporal nature of neuromuscular control

Animals must control numerous muscles to produce forces and movements with their limbs. Current theories of motor optimization and synergistic control are predicated on the assumption that there are multiple highly diverse feasible activations for any motor task ("muscle redundancy"). Here, we demonstrate that the dimensionality of the neuromuscular control problem is greatly reduced when adding the temporal constraints inherent to any sequence of motor commands: the physiological time constants for muscle activation-contraction dynamics. We used a seven-muscle model of a human finger to fully characterize the seven-dimensional polytope of all possible motor commands that can produce fingertip force vector in any direction in 3D, in alignment with the core models of Feasibility Theory. For a given sequence of seven force vectors lasting 300 ms, a novel single-step extended linear program finds the 49-dimensional polytope of all possible motor commands that can produce the sequence of forces. We find that muscle redundancy is severely reduced when the temporal limits on muscle activation-contraction dynamics are added. For example, allowing a generous ± 12% change in muscle activation within 50 ms allows visiting only ∼ 7% of the feasible activation space in the next time step. By considering that every motor command conditions future commands, we find that the motor-control landscape is much more highly structured and spatially constrained than previously recognized. We discuss how this challenges traditional computational and conceptual theories of motor control and neurorehabilitation for which muscle redundancy is a foundational assumption.

Biomechanical modeling for the estimation of muscle forces: toward a common language in biomechanics, medical engineering, and neurosciences

Different research fields, such as biomechanics, medical engineering or neurosciences take part in the development of biomechanical models allowing for the estimation of individual muscle forces involved in motor action. The heterogeneity of the terminology used to describe these models according to the research field is a source of confusion and can hamper collaboration between the different fields. This paper proposes a common language based on lexical disambiguation and a synthesis of the terms used in the literature in order to facilitate the understanding of the different elements of biomechanical modeling for force estimation, without questioning the relevance of the terms used in each field or the different model components or their interest. We suggest that the description should start with an indication of whether the muscle force estimation problem is solved following the physiological movement control (from the nervous drive to the muscle force production) or in the opposite direction. Next, the suitability of the model for force production estimation at a given time or for monitoring over time should be specified. Authors should pay particular attention to the method description used to find solutions, specifying whether this is done during or after data collection, with possible method adaptations during processing. Finally, the presence of additional data must be specified by indicating whether they are used to drive, assist, or calibrate the model. Describing and classifying models in this way will facilitate the use and application in all fields where the estimation of muscle forces is of real, direct, and concrete interest.

Neuromuscular control: from a biomechanist's perspective

Understanding neural control of movement necessitates a collaborative approach between many disciplines, including biomechanics, neuroscience, and motor control. Biomechanics grounds us to the laws of physics that our musculoskeletal system must obey. Neuroscience reveals the inner workings of our nervous system that functions to control our body. Motor control investigates the coordinated motor behaviours we display when interacting with our environment. The combined efforts across the many disciplines aimed at understanding human movement has resulted in a rich and rapidly growing body of literature overflowing with theories, models, and experimental paradigms. As a result, gathering knowledge and drawing connections between the overlapping but seemingly disparate fields can be an overwhelming endeavour. This review paper evolved as a need for us to learn of the diverse perspectives underlying current understanding of neuromuscular control. The purpose of our review paper is to integrate ideas from biomechanics, neuroscience, and motor control to better understand how we voluntarily control our muscles. As biomechanists, we approach this paper starting from a biomechanical modelling framework. We first define the theoretical solutions (i.e., muscle activity patterns) that an individual could feasibly use to complete a motor task. The theoretical solutions will be compared to experimental findings and reveal that individuals display structured muscle activity patterns that do not span the entire theoretical solution space. Prevalent neuromuscular control theories will be discussed in length, highlighting optimality, probabilistic principles, and neuromechanical constraints, that may guide individuals to families of muscle activity solutions within what is theoretically possible. Our intention is for this paper to serve as a primer for the neuromuscular control scientific community by introducing and integrating many of the ideas common across disciplines today, as well as inspire future work to improve the representation of neural control in biomechanical models.

Number of synergies impacts sensitivity of gait to weakness and contracture

Muscle activity during gait can be described by a small set of synergies, weighted groups of muscles, that are theorized to reflect underlying neural control. For people with neurologic injuries, like cerebral palsy or stroke, even fewer synergies are required to explain muscle activity during gait. This reduction in synergies is thought to reflect altered control and is associated with impairment severity and treatment outcomes. Individuals with neurologic injuries also develop secondary musculoskeletal impairments, like weakness or contracture, that can impact gait. Yet, the combined impacts of altered control and musculoskeletal impairments on gait remains unclear. In this study, we use a two-dimensional musculoskeletal model constrained to synergy control to simulate unimpaired gait. We vary the number of synergies, while simulating muscle weakness and contracture to examine how altered control impacts sensitivity to musculoskeletal impairment while tracking unimpaired gait. Results demonstrate that reducing the number of synergies increases sensitivity to weakness and contracture for specific muscle groups. For example, simulations using five-synergy control tolerated 40% and 51% more knee extensor weakness than those using four- or three-synergy control, respectively. Furthermore, when constrained to four- or three-synergy control, the model was increasingly sensitive to contracture and weakness of proximal muscles, such as the hamstring and hip flexors. Contrastingly, neither the amount of generalized nor plantarflexor weakness tolerated was affected by the number of synergies. These findings highlight the interactions between altered control and musculoskeletal impairments, emphasizing the importance of measuring and incorporating both in future simulation and experimental studies.

Effects of kinematic complexity and number of muscles on musculoskeletal model robustness to muscle dysfunction

The robustness of motor outputs to muscle dysfunction has been investigated using musculoskeletal modeling, but with conflicting results owing to differences in model complexity and motor tasks. Our objective was to systematically study how the number of kinematic degrees of freedom, and the number of independent muscle actuators alter the robustness of motor output to muscle dysfunction. We took a detailed musculoskeletal model of the human leg and systematically varied the model complexity to create six models with either 3 or 7 kinematic degrees of freedom and either 14, 26, or 43 muscle actuators. We tested the redundancy of each model by quantifying the reduction in sagittal plane feasible force set area when a single muscle was removed. The robustness of feasible force set area to the loss of any single muscle, i.e. general single muscle loss increased with the number of independent muscles and decreased with the number of kinematic degrees of freedom, with the robust area varying from 1% and 52% of the intact feasible force set area. The maximum sensitivity of the feasible force set to the loss of any single muscle varied from 75% to 26% of the intact feasible force set area as the number of muscles increased. Additionally, the ranges of feasible muscle activation for maximum force production were largely unconstrained in many cases, indicating ample musculoskeletal redundancy even for maximal forces. We propose that ratio of muscles to kinematic degrees of freedom can be used as a rule of thumb for estimating musculoskeletal redundancy in both simulated and real biomechanical systems.