Massimo Sartori

EMG-Driven Forward-Dynamic Estimation of Muscle Force and Joint Moment about Multiple Degrees of Freedom in the Human Lower Extremity

This work examined if currently available electromyography (EMG) driven models, that are calibrated to satisfy joint moments about one single degree of freedom (DOF), could provide the same musculotendon unit (MTU) force solution, when driven by the same input data, but calibrated about a different DOF. We then developed a novel and comprehensive EMG-driven model of the human lower extremity that used EMG signals from 16 muscle groups to drive 34 MTUs and satisfy the resulting joint moments simultaneously produced about four DOFs during different motor tasks. This also led to the development of a calibration procedure that allowed identifying a set of subject-specific parameters that ensured physiological behavior for the 34 MTUs. Results showed that currently available single-DOF models did not provide the same unique MTU force solution for the same input data. On the other hand, the MTU force solution predicted by our proposed multi-DOF model satisfied joint moments about multiple DOFs without loss of accuracy compared to single-DOF models corresponding to each of the four DOFs. The predicted MTU force solution was (1) a function of experimentally measured EMGs, (2) the result of physiological MTU excitation, (3) reflected different MTU contraction strategies associated to different motor tasks, (4) coordinated a greater number of MTUs with respect to currently available single-DOF models, and (5) was not specific to an individual DOF dynamics. Therefore, our proposed methodology has the potential of producing a more dynamically consistent and generalizable MTU force solution than was possible using single-DOF EMG-driven models. This will help better address the important scientific questions previously approached using single-DOF EMG-driven modeling. Furthermore, it might have applications in the development of human-machine interfaces for assistive devices.

Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces.

Estimating tibiofemoral joint contact forces is important for understanding the initiation and progression of knee osteoarthritis. However, tibiofemoral contact force predictions are influenced by many factors including muscle forces and anatomical representations of the knee joint. This study aimed to investigate the influence of subject-specific geometry and knee joint kinematics on the prediction of tibiofemoral contact forces using a calibrated EMG-driven neuromusculoskeletal model of the knee. One participant fitted with an instrumented total knee replacement walked at a self-selected speed while medial and lateral tibiofemoral contact forces, ground reaction forces, whole-body kinematics, and lower-limb muscle activity were simultaneously measured. The combination of generic and subject-specific knee joint geometry and kinematics resulted in four different OpenSim models used to estimate muscle-tendon lengths and moment arms. The subject-specific geometric model was created from CT scans and the subject-specific knee joint kinematics representing the translation of the tibia relative to the femur was obtained from fluoroscopy. The EMG-driven model was calibrated using one walking trial, but with three different cost functions that tracked the knee flexion/extension moments with and without constraint over the estimated joint contact forces. The calibrated models then predicted the medial and lateral tibiofemoral contact forces for five other different walking trials. The use of subject-specific models with minimization of the peak tibiofemoral contact forces improved the accuracy of medial contact forces by 47% and lateral contact forces by 7%, respectively compared with the use of generic musculoskeletal model.

A musculoskeletal model of human locomotion driven by a low dimensional set of impulsive excitation primitives.

Human locomotion has been described as being generated by an impulsive (burst-like) excitation of groups of musculotendon units, with timing dependent on the biomechanical goal of the task. Despite this view being supported by many experimental observations on specific locomotion tasks, it is still unknown if the same impulsive controller (i.e., a low-dimensional set of time-delayed excitastion primitives) can be used as input drive for large musculoskeletal models across different human locomotion tasks. For this purpose, we extracted, with non-negative matrix factorization, five non-negative factors from a large sample of muscle electromyograms in two healthy subjects during four motor tasks. These included walking, running, sidestepping, and crossover cutting maneuvers. The extracted non-negative factors were then averaged and parameterized to obtain task-generic Gaussian-shaped impulsive excitation curves or primitives. These were used to drive a subject-specific musculoskeletal model of the human lower extremity. Results showed that the same set of five impulsive excitation primitives could be used to predict the dynamics of 34 musculotendon units and the resulting hip, knee and ankle joint moments (i.e., NRMSE = 0.18 ± 0.08, and R2 = 0.73 ± 0.22 across all tasks and subjects) without substantial loss of accuracy with respect to using experimental electromyograms (i.e., NRMSE = 0.16 ± 0.07, and R2 = 0.78 ± 0.18 across all tasks and subjects). Results support the hypothesis that biomechanically different motor tasks might share similar neuromuscular control strategies. This might have implications in neurorehabilitation technologies such as human-machine interfaces for the torque-driven, proportional control of powered prostheses and orthoses. In this, device control commands (i.e., predicted joint torque) could be derived without direct experimental data but relying on simple parameterized Gaussian-shaped curves, thus decreasing the input drive complexity and the number of needed sensors.

CEINMS: A toolbox to investigate the influence of different neural control solutions on the prediction of muscle excitation and joint moments during dynamic motor tasks

Personalized neuromusculoskeletal (NMS) models can represent the neurological, physiological, and anatomical characteristics of an individual and can be used to estimate the forces generated inside the human body. Currently, publicly available software to calculate muscle forces are restricted to static and dynamic optimisation methods, or limited to isometric tasks only. We have created and made freely available for the research community the Calibrated EMG-Informed NMS Modelling Toolbox (CEINMS), an OpenSim plug-in that enables investigators to predict different neural control solutions for the same musculoskeletal geometry and measured movements. CEINMS comprises EMG-driven and EMG-informed algorithms that have been previously published and tested. It operates on dynamic skeletal models possessing any number of degrees of freedom and musculotendon units and can be calibrated to the individual to predict measured joint moments and EMG patterns. In this paper we describe the components of CEINMS and its integration with OpenSim. We then analyse how EMG-driven, EMG-assisted, and static optimisation neural control solutions affect the estimated joint moments, muscle forces, and muscle excitations, including muscle co-contraction.

Modeling and simulating the neuromuscular mechanisms regulating ankle and knee joint stiffness during human locomotion

This work presents an electrophysiologically and dynamically consistent musculoskeletal model to predict stiffness in the human ankle and knee joints as derived from the joints constituent biological tissues (i.e., the spanning musculotendon units). The modeling method we propose uses electromyography (EMG) recordings from 13 muscle groups to drive forward dynamic simulations of the human leg in five healthy subjects during overground walking and running. The EMG-driven musculoskeletal model estimates musculotendon and resulting joint stiffness that is consistent with experimental EMG data as well as with the experimental joint moments. This provides a framework that allows for the first time observing 1) the elastic interplay between the knee and ankle joints, 2) the individual muscle contribution to joint stiffness, and 3) the underlying co-contraction strategies. It provides a theoretical description of how stiffness modulates as a function of muscle activation, fiber contraction, and interacting tendon dynamics. Furthermore, it describes how this differs from currently available stiffness definitions, including quasi-stiffness and short-range stiffness. This work offers a theoretical and computational basis for describing and investigating the neuromuscular mechanisms underlying human locomotion.

Neural Data-Driven Musculoskeletal Modeling for Personalized Neurorehabilitation Technologies

Objectives: The development of neurorehabilitation technologies requires the profound understanding of the mechanisms underlying an individuals motor ability and impairment. A major factor limiting this understanding is the difficulty of bridging between events taking place at the neurophysiologic level (i.e., motor neuron firings) with those emerging at the musculoskeletal level (i.e. joint actuation), in vivo in the intact moving human. This review presents emerging model-based methodologies for filling this gap that are promising for developing clinically viable technologies. Methods: We provide a design overview of musculoskeletal modeling formulations driven by recordings of neuromuscular activity with a critical comparison to alternative model-free approaches in the context of neurorehabilitation technologies. We present advanced electromyography-based techniques for interfacing with the human nervous system and model-based techniques for translating the extracted neural information into estimates of motor function. Results: We introduce representative application areas where modeling is relevant for accessing neuromuscular variables that could not be measured experimentally. We then show how these variables are used for designing personalized rehabilitation interventions, biologically inspired limbs, and human-machine interfaces. Conclusion: The ability of using electrophysiological recordings to inform biomechanical models enables accessing a broader range of neuromechanical variables than analyzing electrophysiological data or movement data individually. This enables understanding the neuromechanical interplay underlying in vivo movement function, pathology, and robot-assisted motor recovery. Significance: Filling the gap between our understandings of movement neural and mechanical functions is central for addressing one of the major challenges in neurorehabilitation: personalizing current technologies and interventions to an individuals anatomy and impairment.

Human-like compliant locomotion: state of the art of robotic implementations

This review paper provides a synthetic yet critical overview of the key biomechanical principles of human bipedal walking and their current implementation in robotic platforms. We describe the functional role of human joints, addressing in particular the relevance of the compliant properties of the different degrees of freedom throughout the gait cycle. We focused on three basic functional units involved in locomotion, i.e. the ankle-foot complex, the knee, and the hip-pelvis complex, and their relevance to whole-body performance. We present an extensive review of the current implementations of these mechanisms into robotic platforms, discussing their potentialities and limitations from the functional and energetic perspectives. We specifically targeted humanoid robots, but also revised evidence from the field of lower-limb prosthetics, which presents innovative solutions still unexploited in the current humanoids. Finally, we identified the main critical aspects of the process of translating human principles into actual machines, providing a number of relevant challenges that should be addressed in future research.

A predictive model of muscle excitations based on muscle modularity for a large repertoire of human locomotion conditions

Humans can efficiently walk across a large variety of terrains and locomotion conditions with little or no mental effort. It has been hypothesized that the nervous system simplifies neuromuscular control by using muscle synergies, thus organizing multi-muscle activity into a small number of coordinative co-activation modules. In the present study we investigated how muscle modularity is structured across a large repertoire of locomotion conditions including five different speeds and five different ground elevations. For this we have used the non-negative matrix factorization technique in order to explain EMG experimental data with a low-dimensional set of four motor components. In this context each motor components is composed of a non-negative factor and the associated muscle weightings. Furthermore, we have investigated if the proposed descriptive analysis of muscle modularity could be translated into a predictive model that could: (1) Estimate how motor components modulate across locomotion speeds and ground elevations. This implies not only estimating the non-negative factors temporal characteristics, but also the associated muscle weighting variations. (2) Estimate how the resulting muscle excitations modulate across novel locomotion conditions and subjects. The results showed three major distinctive features of muscle modularity: (1) the number of motor components was preserved across all locomotion conditions, (2) the non-negative factors were consistent in shape and timing across all locomotion conditions, and (3) the muscle weightings were modulated as distinctive functions of locomotion speed and ground elevation. Results also showed that the developed predictive model was able to reproduce well the muscle modularity of un-modeled data, i.e., novel subjects and conditions. Muscle weightings were reconstructed with a cross-correlation factor greater than 70% and a root mean square error less than 0.10. Furthermore, the generated muscle excitations matched well the experimental excitation with a cross-correlation factor greater than 85% and a root mean square error less than 0.09. The ability of synthetizing the neuromuscular mechanisms underlying human locomotion across a variety of locomotion conditions will enable solutions in the field of neurorehabilitation technologies and control of bipedal artificial systems. Open-access of the model implementation is provided for further analysis at https://simtk.org/home/p-mep/.

MOtoNMS: A MATLAB toolbox to process motion data for neuromusculoskeletal modeling and simulation

BackgroundNeuromusculoskeletal modeling and simulation enable investigation of the neuromusculoskeletal system and its role in human movement dynamics. These methods are progressively introduced into daily clinical practice. However, a major factor limiting this translation is the lack of robust tools for the pre-processing of experimental movement data for their use in neuromusculoskeletal modeling software.ResultsThis paper presents MOtoNMS (matlab MOtion data elaboration TOolbox for NeuroMusculoSkeletal applications), a toolbox freely available to the community, that aims to fill this lack. MOtoNMS processes experimental data from different motion analysis devices and generates input data for neuromusculoskeletal modeling and simulation software, such as OpenSim and CEINMS (Calibrated EMG-Informed NMS Modelling Toolbox). MOtoNMS implements commonly required processing steps and its generic architecture simplifies the integration of new user-defined processing components. MOtoNMS allows users to setup their laboratory configurations and processing procedures through user-friendly graphical interfaces, without requiring advanced computer skills. Finally, configuration choices can be stored enabling the full reproduction of the processing steps. MOtoNMS is released under GNU General Public License and it is available at the SimTK website and from the GitHub repository. Motion data collected at four institutions demonstrate that, despite differences in laboratory instrumentation and procedures, MOtoNMS succeeds in processing data and producing consistent inputs for OpenSim and CEINMS.ConclusionsMOtoNMS fills the gap between motion analysis and neuromusculoskeletal modeling and simulation. Its support to several devices, a complete implementation of the pre-processing procedures, its simple extensibility, the available user interfaces, and its free availability can boost the translation of neuromusculoskeletal methods in daily and clinical practice.

Bone remodelling in the natural acetabulum is influenced by muscle force‐induced bone stress

A modelling framework using the international Physiome Project is presented for evaluating the role of muscles on acetabular stress patterns in the natural hip. The novel developments include the following: (i) an efficient method for model generation with validation; (ii) the inclusion of electromyography-estimated muscle forces from gait; and (iii) the role that muscles play in the hip stress pattern. The 3D finite element hip model includes anatomically based muscle area attachments, material properties derived from Hounsfield units and validation against an Instron compression test. The primary outcome from this study is that hip loading applied as anatomically accurate muscle forces redistributes the stress pattern and reduces peak stress throughout the pelvis and within the acetabulum compared with applying the same net hip force without muscles through the femur. Muscle forces also increased stress where large muscles have small insertion sites. This has implications for the hip where bone stress and strain are key excitation variables used to initiate bone remodelling based on the strain-based bone remodelling theory. Inclusion of muscle forces reduces the predicted sites and degree of remodelling. The secondary outcome is that the key muscles that influenced remodelling in the acetabulum were the rectus femoris, adductor magnus and iliacus.