Ajay Seth

Muscle contributions to propulsion and support during running

Muscles actuate running by developing forces that propel the body forward while supporting the bodys weight. To understand how muscles contribute to propulsion (i.e., forward acceleration of the mass center) and support (i.e., upward acceleration of the mass center) during running we developed a three-dimensional muscle-actuated simulation of the running gait cycle. The simulation is driven by 92 musculotendon actuators of the lower extremities and torso and includes the dynamics of arm motion. We analyzed the simulation to determine how each muscle contributed to the acceleration of the body mass center. During the early part of the stance phase, the quadriceps muscle group was the largest contributor to braking (i.e., backward acceleration of the mass center) and support. During the second half of the stance phase, the soleus and gastrocnemius muscles were the greatest contributors to propulsion and support. The arms did not contribute substantially to either propulsion or support, generating less than 1% of the peak mass center acceleration. However, the arms effectively counterbalanced the vertical angular momentum of the lower extremities. Our analysis reveals that the quadriceps and plantarflexors are the major contributors to acceleration of the body mass center during running.

Flexing computational muscle: modeling and simulation of musculotendon dynamics.

Muscle-driven simulations of human and animal motion are widely used to complement physical experiments for studying movement dynamics. Musculotendon models are an essential component of muscle-driven simulations, yet neither the computational speed nor the biological accuracy of the simulated forces has been adequately evaluated. Here we compare the speed and accuracy of three musculotendon models: two with an elastic tendon (an equilibrium model and a damped equilibrium model) and one with a rigid tendon. Our simulation benchmarks demonstrate that the equilibrium and damped equilibrium models produce similar force profiles but have different computational speeds. At low activation, the damped equilibrium model is 29 times faster than the equilibrium model when using an explicit integrator and 3 times faster when using an implicit integrator; at high activation, the two models have similar simulation speeds. In the special case of simulating a muscle with a short tendon, the rigid-tendon model produces forces that match those generated by the elastic-tendon models, but simulates 2-54 times faster when an explicit integrator is used and 6-31 times faster when an implicit integrator is used. The equilibrium, damped equilibrium, and rigid-tendon models reproduce forces generated by maximally-activated biological muscle with mean absolute errors less than 8.9%, 8.9%, and 20.9% of the maximum isometric muscle force, respectively. When compared to forces generated by submaximally-activated biological muscle, the forces produced by the equilibrium, damped equilibrium, and rigid-tendon models have mean absolute errors less than 16.2%, 16.4%, and 18.5%, respectively. To encourage further development of musculotendon models, we provide implementations of each of these models in OpenSim version 3.1 and benchmark data online, enabling others to reproduce our results and test their models of musculotendon dynamics.

Muscle contributions to support and progression during single-limb stance in crouch gait

Pathological movement patterns like crouch gait are characterized by abnormal kinematics and muscle activations that alter how muscles support the body weight during walking. Individual muscles are often the target of interventions to improve crouch gait, yet the roles of individual muscles during crouch gait remain unknown. The goal of this study was to examine how muscles contribute to mass center accelerations and joint angular accelerations during single-limb stance in crouch gait, and compare these contributions to unimpaired gait. Subject-specific dynamic simulations were created for ten children who walked in a mild crouch gait and had no previous surgeries. The simulations were analyzed to determine the acceleration of the mass center and angular accelerations of the hip, knee, and ankle generated by individual muscles. The results of this analysis indicate that children walking in crouch gait have less passive skeletal support of body weight and utilize substantially higher muscle forces to walk than unimpaired individuals. Crouch gait relies on the same muscles as unimpaired gait to accelerate the mass center upward, including the soleus, vasti, gastrocnemius, gluteus medius, rectus femoris, and gluteus maximus. However, during crouch gait, these muscles are active throughout single-limb stance, in contrast to the modulation of muscle forces seen during single-limb stance in an unimpaired gait. Subjects walking in crouch gait rely more on proximal muscles, including the gluteus medius and hamstrings, to accelerate the mass center forward during single-limb stance than subjects with an unimpaired gait.

Is My Model Good Enough? Best Practices for Verification and Validation of Musculoskeletal Models and Simulations of Movement

Computational modeling and simulation of neuromusculoskeletal (NMS) systems enables researchers and clinicians to study the complex dynamics underlying human and animal movement. NMS models use equations derived from physical laws and biology to help solve challenging real-world problems, from designing prosthetics that maximize running speed to developing exoskeletal devices that enable walking after a stroke. NMS modeling and simulation has proliferated in the biomechanics research community over the past 25 years, but the lack of verification and validation standards remains a major barrier to wider adoption and impact. The goal of this paper is to establish practical guidelines for verification and validation of NMS models and simulations that researchers, clinicians, reviewers, and others can adopt to evaluate the accuracy and credibility of modeling studies. In particular, we review a general process for verification and validation applied to NMS models and simulations, including careful formulation of a research question and methods, traditional verification and validation steps, and documentation and sharing of results for use and testing by other researchers. Modeling the NMS system and simulating its motion involves methods to represent neural control, musculoskeletal geometry, muscle-tendon dynamics, contact forces, and multibody dynamics. For each of these components, we review modeling choices and software verification guidelines; discuss variability, errors, uncertainty, and sensitivity relationships; and provide recommendations for verification and validation by comparing experimental data and testing robustness. We present a series of case studies to illustrate key principles. In closing, we discuss challenges the community must overcome to ensure that modeling and simulation are successfully used to solve the broad spectrum of problems that limit human mobility.

Are Subject-Specific Musculoskeletal Models Robust to the Uncertainties in Parameter Identification?

Subject-specific musculoskeletal modeling can be applied to study musculoskeletal disorders, allowing inclusion of personalized anatomy and properties. Independent of the tools used for model creation, there are unavoidable uncertainties associated with parameter identification, whose effect on model predictions is still not fully understood. The aim of the present study was to analyze the sensitivity of subject-specific model predictions (i.e., joint angles, joint moments, muscle and joint contact forces) during walking to the uncertainties in the identification of body landmark positions, maximum muscle tension and musculotendon geometry. To this aim, we created an MRI-based musculoskeletal model of the lower limbs, defined as a 7-segment, 10-degree-of-freedom articulated linkage, actuated by 84 musculotendon units. We then performed a Monte-Carlo probabilistic analysis perturbing model parameters according to their uncertainty, and solving a typical inverse dynamics and static optimization problem using 500 models that included the different sets of perturbed variable values. Model creation and gait simulations were performed by using freely available software that we developed to standardize the process of model creation, integrate with OpenSim and create probabilistic simulations of movement. The uncertainties in input variables had a moderate effect on model predictions, as muscle and joint contact forces showed maximum standard deviation of 0.3 times body-weight and maximum range of 2.1 times body-weight. In addition, the output variables significantly correlated with few input variables (up to 7 out of 312) across the gait cycle, including the geometry definition of larger muscles and the maximum muscle tension in limited gait portions. Although we found subject-specific models not markedly sensitive to parameter identification, researchers should be aware of the model precision in relation to the intended application. In fact, force predictions could be affected by an uncertainty in the same order of magnitude of its value, although this condition has low probability to occur.

Contributions of muscles to mediolateral ground reaction force over a range of walking speeds

Impaired control of mediolateral body motion during walking is an important health concern. Developing treatments to improve mediolateral control is challenging, partly because the mechanisms by which muscles modulate mediolateral ground reaction force (and thereby modulate mediolateral acceleration of the body mass center) during unimpaired walking are poorly understood. To investigate this, we examined mediolateral ground reaction forces in eight unimpaired subjects walking at four speeds and determined the contributions of muscles, gravity, and velocity-related forces to the mediolateral ground reaction force by analyzing muscle-driven simulations of these subjects. During early stance (0-6% gait cycle), peak ground reaction force on the leading foot was directed laterally and increased significantly (p<0.05) with walking speed. During early single support (14-30% gait cycle), peak ground reaction force on the stance foot was directed medially and increased significantly (p<0.01) with speed. Muscles accounted for more than 92% of the mediolateral ground reaction force over all walking speeds, whereas gravity and velocity-related forces made relatively small contributions. Muscles coordinate mediolateral acceleration via an interplay between the medial ground reaction force contributed by the abductors and the lateral ground reaction forces contributed by the knee extensors, plantarflexors, and adductors. Our findings show how muscles that contribute to forward progression and body-weight support also modulate mediolateral acceleration of the body mass center while weight is transferred from one leg to another during double support.

Muscle contributions to vertical and fore-aft accelerations are altered in subjects with crouch gait

The goals of this study were to determine if the muscle contributions to vertical and fore-aft acceleration of the mass center differ between crouch gait and unimpaired gait and if these muscle contributions change with crouch severity. Examining muscle contributions to mass center acceleration provides insight into the roles of individual muscles during gait and can provide guidance for treatment planning. We calculated vertical and fore-aft accelerations using musculoskeletal simulations of typically developing children and children with cerebral palsy and crouch gait. Analysis of these simulations revealed that during unimpaired gait the quadriceps produce large upward and backward accelerations during early stance, whereas the ankle plantarflexors produce large upward and forward accelerations later in stance. In contrast, during crouch gait, the quadriceps and ankle plantarflexors produce large, opposing fore-aft accelerations throughout stance. The quadriceps force required to accelerate the mass center upward was significantly larger in crouch gait than in unimpaired gait and increased with crouch severity. The gluteus medius accelerated the mass center upward during midstance in unimpaired gait; however, during crouch gait the upward acceleration produced by the gluteus medius was significantly reduced. During unimpaired gait the quadriceps and ankle plantarflexors accelerate the mass center at different times, efficiently modulating fore-aft accelerations. However, during crouch gait, the quadriceps and ankle plantarflexors produce fore-aft accelerations at the same time and the opposing fore-aft accelerations generated by these muscles contribute to the inefficiency of crouch gait.

A Biomechanical Model of the Scapulothoracic Joint to Accurately Capture Scapular Kinematics during Shoulder Movements

The complexity of shoulder mechanics combined with the movement of skin relative to the scapula makes it difficult to measure shoulder kinematics with sufficient accuracy to distinguish between symptomatic and asymptomatic individuals. Multibody skeletal models can improve motion capture accuracy by reducing the space of possible joint movements, and models are used widely to improve measurement of lower limb kinematics. In this study, we developed a rigid-body model of a scapulothoracic joint to describe the kinematics of the scapula relative to the thorax. This model describes scapular kinematics with four degrees of freedom: 1) elevation and 2) abduction of the scapula on an ellipsoidal thoracic surface, 3) upward rotation of the scapula normal to the thoracic surface, and 4) internal rotation of the scapula to lift the medial border of the scapula off the surface of the thorax. The surface dimensions and joint axes can be customized to match an individual’s anthropometry. We compared the model to “gold standard” bone-pin kinematics collected during three shoulder tasks and found modeled scapular kinematics to be accurate to within 2mm root-mean-squared error for individual bone-pin markers across all markers and movement tasks. As an additional test, we added random and systematic noise to the bone-pin marker data and found that the model reduced kinematic variability due to noise by 65% compared to Euler angles computed without the model. Our scapulothoracic joint model can be used for inverse and forward dynamics analyses and to compute joint reaction loads. The computational performance of the scapulothoracic joint model is well suited for real-time applications; it is freely available for use with OpenSim 3.2, and is customizable and usable with other OpenSim models.

A rolling constraint reproduces ground reaction forces and moments in dynamic simulations of walking, running, and crouch gait

Recent advances in computational technology have dramatically increased the use of muscle-driven simulation to study accelerations produced by muscles during gait. Accelerations computed from muscle-driven simulations are sensitive to the model used to represent contact between the foot and ground. A foot-ground contact model must be able to calculate ground reaction forces and moments that are consistent with experimentally measured ground reaction forces and moments. We show here that a rolling constraint can model foot-ground contact and reproduce measured ground reaction forces and moments in an induced acceleration analysis of muscle-driven simulations of walking, running, and crouch gait. We also illustrate that a point constraint and a weld constraint used to model foot-ground contact in previous studies produce inaccurate reaction moments and lead to contradictory interpretations of muscle function. To enable others to use and test these different constraint types (i.e., rolling, point, and weld constraints) we have included them as part of an induced acceleration analysis in OpenSim, a freely-available biomechanics simulation package.

Simulating Ideal Assistive Devices to Reduce the Metabolic Cost of Running.

Tools have been used for millions of years to augment the capabilities of the human body, allowing us to accomplish tasks that would otherwise be difficult or impossible. Powered exoskeletons and other assistive devices are sophisticated modern tools that have restored bipedal locomotion in individuals with paraplegia and have endowed unimpaired individuals with superhuman strength. Despite these successes, designing assistive devices that reduce energy consumption during running remains a substantial challenge, in part because these devices disrupt the dynamics of a complex, finely tuned biological system. Furthermore, designers have hitherto relied primarily on experiments, which cannot report muscle-level energy consumption and are fraught with practical challenges. In this study, we use OpenSim to generate muscle-driven simulations of 10 human subjects running at 2 and 5 m/s. We then add ideal, massless assistive devices to our simulations and examine the predicted changes in muscle recruitment patterns and metabolic power consumption. Our simulations suggest that an assistive device should not necessarily apply the net joint moment generated by muscles during unassisted running, and an assistive device can reduce the activity of muscles that do not cross the assisted joint. Our results corroborate and suggest biomechanical explanations for similar effects observed by experimentalists, and can be used to form hypotheses for future experimental studies. The models, simulations, and software used in this study are freely available at simtk.org and can provide insight into assistive device design that complements experimental approaches.