Development of a Subject-Specific Foot-Ground Contact Model for Walking

Computational evaluation of psoas muscle influence on walking function following internal hemipelvectomy with reconstruction

An emerging option for internal hemipelvectomy surgery is custom prosthesis reconstruction. This option typically recapitulates the resected pelvic bony anatomy with the goal of maximizing post-surgery walking function while minimizing recovery time. However, the current custom prosthesis design process does not account for the patient's post-surgery prosthesis and bone loading patterns, nor can it predict how different surgical or rehabilitation decisions (e.g., retention or removal of the psoas muscle, strengthening the psoas) will affect prosthesis durability and post-surgery walking function. These factors may contribute to the high observed failure rate for custom pelvic prostheses, discouraging orthopedic oncologists from pursuing this valuable treatment option. One possibility for addressing this problem is to simulate the complex interaction between surgical and rehabilitation decisions, post-surgery walking function, and custom pelvic prosthesis design using patient-specific neuromusculoskeletal models. As a first step toward developing this capability, this study used a personalized neuromusculoskeletal model and direct collocation optimal control to predict the impact of ipsilateral psoas muscle strength on walking function following internal hemipelvectomy with custom prosthesis reconstruction. The influence of the psoas muscle was targeted since retention of this important muscle can be surgically demanding for certain tumors, requiring additional time in the operating room. The post-surgery walking predictions emulated the most common surgical scenario encountered at MD Anderson Cancer Center in Houston. Simulated post-surgery psoas strengths included 0% (removed), 50% (weakened), 100% (maintained), and 150% (strengthened) of the pre-surgery value. However, only the 100% and 150% cases successfully converged to a complete gait cycle. When post-surgery psoas strength was maintained, clinical gait features were predicted, including increased stance width, decreased stride length, and increased lumbar bending towards the operated side. Furthermore, when post-surgery psoas strength was increased, stance width and stride length returned to pre-surgery values. These results suggest that retention and strengthening of the psoas muscle on the operated side may be important for maximizing post-surgery walking function. If future studies can validate this computational approach using post-surgery experimental walking data, the approach may eventually influence surgical, rehabilitation, and custom prosthesis design decisions to meet the unique clinical needs of pelvic sarcoma patients.

Inertial Sensor-to-Segment Calibration for Accurate 3D Joint Angle Calculation for Use in OpenSim

Inertial capture (InCap) systems combined with musculoskeletal (MSK) models are an attractive option for monitoring 3D joint kinematics in an ecological context. However, the primary limiting factor is the sensor-to-segment calibration, which is crucial to estimate the body segment orientations. Walking, running, and stair ascent and descent trials were measured in eleven healthy subjects with the Xsens InCap system and the Vicon 3D motion capture (MoCap) system at a self-selected speed. A novel integrated method that combines previous sensor-to-segment calibration approaches was developed for use in a MSK model with three degree of freedom (DOF) hip and knee joints. The following were compared: RMSE, range of motion (ROM), peaks, and R² between InCap kinematics estimated with different calibration methods and gold standard MoCap kinematics. The integrated method reduced the RSME for both the hip and the knee joints below 5°, and no statistically significant differences were found between MoCap and InCap kinematics. This was consistent across all the different analyzed movements. The developed method was integrated on an MSK model workflow, and it increased the sensor-to-segment calibration accuracy for an accurate estimate of 3D joint kinematics compared to MoCap, guaranteeing a clinical easy-to-use approach.

Evaluation of Optimal Control Approaches for Predicting Active Knee-Ankle-Foot-Orthosis Motion for Individuals With Spinal Cord Injury

Gait restoration of individuals with spinal cord injury can be partially achieved using active orthoses or exoskeletons. To improve the walking ability of each patient as much as possible, it is important to personalize the parameters that define the device actuation. This study investigates whether using an optimal control-based predictive simulation approach to personalize pre-defined knee trajectory parameters for an active knee-ankle-foot orthosis (KAFO) used by spinal cord injured (SCI) subjects could potentially be an alternative to the current trial-and-error approach. We aimed to find the knee angle trajectory that produced an improved orthosis-assisted gait pattern compared to the one with passive support (locked knee). We collected experimental data from a healthy subject assisted by crutches and KAFOs (with locked knee and with knee flexion assistance) and from an SCI subject assisted by crutches and KAFOs (with locked knee). First, we compared different cost functions and chose the one that produced results closest to experimental locked knee walking for the healthy subject (angular coordinates mean RMSE was 5.74°). For this subject, we predicted crutch-orthosis-assisted walking imposing a pre-defined knee angle trajectory for different maximum knee flexion parameter values, and results were evaluated against experimental data using that same pre-defined knee flexion trajectories in the real device. Finally, using the selected cost function, gait cycles for different knee flexion assistance were predicted for an SCI subject. We evaluated changes in four clinically relevant parameters: foot clearance, stride length, cadence, and hip flexion ROM. Simulations for different values of maximum knee flexion showed variations of these parameters that were consistent with experimental data for the healthy subject (e.g., foot clearance increased/decreased similarly in experimental and predicted motions) and were reasonable for the SCI subject (e.g., maximum parameter values were found for moderate knee flexion). Although more research is needed before this method can be applied to choose optimal active orthosis controller parameters for specific subjects, these findings suggest that optimal control prediction of crutch-orthosis-assisted walking using biomechanical models might be used in place of the trial-and-error method to select the best maximum knee flexion angle during gait for a specific SCI subject.

A Computational Gait Model With a Below-Knee Amputation and a Semi-Active Variable-Stiffness Foot Prosthesis

INTRODUCTION

Simulations based on computational musculoskeletal models are powerful tools for evaluating the effects of potential biomechanical interventions, such as implementing a novel prosthesis. However, the utility of simulations to evaluate the effects of varied prosthesis design parameters on gait mechanics has not been fully realized due to the lack of a readily-available limb loss-specific gait model and methods for efficiently modeling the energy storage and return dynamics of passive foot prostheses. The purpose of this study was to develop and validate a forward simulation-capable gait model with lower-limb loss and a semi-active variable-stiffness foot (VSF) prosthesis.

METHODS

A seven-segment 28-DoF gait model was developed and forward kinematics simulations, in which experimentally observed joint kinematics were applied and the resulting contact forces under the prosthesis evolved accordingly, were computed for four subjects with unilateral below-knee amputation walking with a VSF.

RESULTS

Model-predicted resultant ground reaction force (GRFR) matched well under trial-specific optimized parameter conditions (mean R2: 0.97, RMSE: 7.7% body weight (BW)) and unoptimized (subject-specific, but not trial-specific) parameter conditions (mean R2: 0.93, RMSE: 12% BW). Simulated anterior-posterior center of pressure demonstrated a mean R2 = 0.64 and RMSE = 14% foot length. Simulated kinematics remained consistent with input data (0.23 deg RMSE, R2 > 0.99) for all conditions.

CONCLUSIONS

These methods may be useful for simulating gait among individuals with lower-limb loss and predicting GRFR arising from gait with novel VSF prostheses. Such data are useful to optimize prosthesis design parameters on a user-specific basis.

A Reduced-Order Computational Model of a Semi-Active Variable-Stiffness Foot Prosthesis

Passive energy storage and return (ESR) feet are current performance standard in lower limb prostheses. A recently developed semi-active variable-stiffness foot (VSF) prosthesis balances the simplicity of a passive ESR device with the adaptability of a powered design. The purpose of this study was to model and simulate the ESR properties of the VSF prosthesis. The ESR properties of the VSF were modeled as a lumped parameter overhung beam. The overhung length is variable, allowing the model to exhibit variable ESR stiffness. Foot-ground contact was modeled using sphere-to-plane contact models. Contact parameters were optimized to represent the geometry and dynamics of the VSF and its foam base. Static compression tests and gait were simulated. Simulation outcomes were compared to corresponding experimental data. Stiffness of the model matched that of the physical VSF (R2: 0.98, root-mean-squared error (RMSE): 1.37 N/mm). Model-predicted resultant ground reaction force (GRFR) matched well under optimized parameter conditions (R2: 0.98, RMSE: 5.3% body weight,) and unoptimized parameter conditions (R2: 0.90, mean RMSE: 13% body weight). Anterior-posterior center of pressure matched well with R2 > 0.94 and RMSE < 9.5% foot length in all conditions. The ESR properties of the VSF were accurately simulated under benchtop testing and dynamic gait conditions. These methods may be useful for predicting GRFR arising from gait with novel prostheses. Such data are useful to optimize prosthesis design parameters on a user-specific basis.

Forward dynamic simulation of Japanese macaque bipedal locomotion demonstrates better energetic economy in a virtualised plantigrade posture

A plantigrade foot with a large robust calcaneus is regarded as a distinctive morphological feature of the human foot; it is presumably the result of adaptation for habitual bipedal locomotion. The foot of the Japanese macaque, on the other hand, does not have such a feature, which hampers it from making foot-ground contact at the heel during bipedal locomotion. Understanding how this morphological difference functionally affects the generation of bipedal locomotion is crucial for elucidating the evolution of human bipedalism. In this study, we constructed a forward dynamic simulation of bipedal locomotion in the Japanese macaque based on a neuromusculoskeletal model to evaluate how virtual manipulation of the foot structure from digitigrade to plantigrade affects the kinematics, dynamics, and energetics of bipedal locomotion in a nonhuman primate whose musculoskeletal anatomy is not adapted to bipedalism. The normal bipedal locomotion generated was in good agreement with that of actual Japanese macaques. If, as in human walking, the foot morphology was altered to allow heel contact, the vertical ground reaction force profile became double-peaked and the cost of transport decreased. These results suggest that evolutionary changes in the foot structure were important for the acquisition of human-like efficient bipedal locomotion.

A Conceptual Blueprint for Making Neuromusculoskeletal Models Clinically Useful

The ultimate goal of most neuromusculoskeletal modeling research is to improve the treatment of movement impairments. However, even though neuromusculoskeletal models have become more realistic anatomically, physiologically, and neurologically over the past 25 years, they have yet to make a positive impact on the design of clinical treatments for movement impairments. Such impairments are caused by common conditions such as stroke, osteoarthritis, Parkinson’s disease, spinal cord injury, cerebral palsy, limb amputation, and even cancer. The lack of clinical impact is somewhat surprising given that comparable computational technology has transformed the design of airplanes, automobiles, and other commercial products over the same time period. This paper provides the author’s personal perspective for how neuromusculoskeletal models can become clinically useful. First, the paper motivates the potential value of neuromusculoskeletal models for clinical treatment design. Next, it highlights five challenges to achieving clinical utility and provides suggestions for how to overcome them. After that, it describes clinical, technical, collaboration, and practical needs that must be addressed for neuromusculoskeletal models to fulfill their clinical potential, along with recommendations for meeting them. Finally, it discusses how more complex modeling and experimental methods could enhance neuromusculoskeletal model fidelity, personalization, and utilization. The author hopes that these ideas will provide a conceptual blueprint that will help the neuromusculoskeletal modeling research community work toward clinical utility.

A probabilistic method to estimate gait kinetics in the absence of ground reaction force measurements

Human joint torques during gait are usually computed using inverse dynamics. This method requires a skeletal model, kinematics and measured ground reaction forces and moments (GRFM). Measuring GRFM is however only possible in a controlled environment. This paper introduces a probabilistic method based on probabilistic principal component analysis to estimate the joint torques for healthy gait without measured GRFM. A gait dataset of 23 subjects was obtained containing kinematics, measured GRFM and joint torques from inverse dynamics in order to obtain a probabilistic model. This model was then used to estimate the joint torques of other subjects without measured GRFM. Only kinematics, a skeletal model and timing of gait events are needed. Estimation only takes 0.28 ms per time instant. Using cross-validation, the resulting root mean square estimation errors for the lower-limb joint torques are found to be approximately 0.1 Nm/kg, which is 6-18% of the range of the ground truth joint torques. Estimated joint torque and GRFM errors are up to two times smaller than model-based state-of-the-art methods. Model-free artificial neural networks can achieve lower errors than our method, but are less repeatable, do not contain uncertainty information on the estimates and are difficult to use in situations which are not in the learning set. In contrast, our method performs well in a new situation where the walking speed is higher than in the learning dataset. The method can for example be used to estimate the kinetics during overground walking without force plates, during treadmill walking without (separate) force plates and during ambulatory measurements.

Computational Design of FastFES Treatment to Improve Propulsive Force Symmetry During Post-stroke Gait: A Feasibility Study

Stroke is a leading cause of long-term disability worldwide and often impairs walking ability. To improve recovery of walking function post-stroke, researchers have investigated the use of treatments such as fast functional electrical stimulation (FastFES). During FastFES treatments, individuals post-stroke walk on a treadmill at their fastest comfortable speed while electrical stimulation is delivered to two muscles of the paretic ankle, ideally to improve paretic leg propulsion and toe clearance. However, muscle selection and stimulation timing are currently standardized based on clinical intuition and a one-size-fits-all approach, which may explain in part why some patients respond to FastFES training while others do not. This study explores how personalized neuromusculoskeletal models could potentially be used to enable individual-specific selection of target muscles and stimulation timing to address unique functional limitations of individual patients post-stroke. Treadmill gait data, including EMG, surface marker positions, and ground reactions, were collected from an individual post-stroke who was a non-responder to FastFES treatment. The patient's gait data were used to personalize key aspects of a full-body neuromusculoskeletal walking model, including lower-body joint functional axes, lower-body muscle force generating properties, deformable foot-ground contact properties, and paretic and non-paretic leg neural control properties. The personalized model was utilized within a direct collocation optimal control framework to reproduce the patient's unstimulated treadmill gait data (verification problem) and to generate three stimulated walking predictions that sought to minimize inter-limb propulsive force asymmetry (prediction problems). The three predictions used: (1) Standard muscle selection (gastrocnemius and tibialis anterior) with standard stimulation timing, (2) Standard muscle selection with optimized stimulation timing, and (3) Optimized muscle selection (soleus and semimembranosus) with optimized stimulation timing. Relative to unstimulated walking, the optimal control problems predicted a 41% reduction in propulsive force asymmetry for scenario (1), a 45% reduction for scenario (2), and a 64% reduction for scenario (3), suggesting that non-standard muscle selection may be superior for this patient. Despite these predicted improvements, kinematic symmetry was not noticeably improved for any of the walking predictions. These results suggest that personalized neuromusculoskeletal models may be able to predict personalized FastFES training prescriptions that could improve propulsive force symmetry, though inclusion of kinematic requirements would be necessary to improve kinematic symmetry as well.

A Quick Turn of Foot: Rigid Foot-Ground Contact Models for Human Motion Prediction

Computer simulation can be used to predict human walking motions as a tool of basic science, device design, and for surgical planning. One the challenges of predicting human walking is accurately synthesizing both the movements and ground forces of the stance foot. Though the foot is commonly modeled as a viscoelastic element, rigid foot-ground contact models offer some advantages: fitting is reduced to a geometric problem, and the numerical stiffness of the equations of motion is similar in both swing and stance. In this work, we evaluate two rigid-foot ground contact models: the ellipse-foot (a single-segment foot), and the double-circle foot (a two-segment foot). To evaluate the foot models we use three different comparisons to experimental data: first we compare how accurately the kinematics of the ankle frame fit those of the model when it is forced to track the measured center-of-pressure (CoP) kinematics; second, we compare how each foot affects how accuracy of a sagittal plane gait model that tracks a subjects walking motion; and third, we assess how each model affects a walking motion prediction. For the prediction problem we consider a unique cost function that includes terms related to both muscular effort and foot-ground impacts. Although the ellipse-foot is superior to the double-circle foot in terms of fit and the accuracy of the tracking OCP solution, the predictive simulation reveals that the ellipse-foot is capable of producing large force transients due to its geometry: when the ankle quickly traverses its u-shaped trajectory, the body is accelerated the body upwards, and large ground forces result. In contrast, the two-segment double-circle foot produces ground forces that are of a similar magnitude to the experimental subject because the additional forefoot segment plastically contacts the ground, arresting its motion, similar to a human foot.

Three-dimensional data-tracking dynamic optimization simulations of human locomotion generated by direct collocation

The aim of this study was to perform full-body three-dimensional (3D) dynamic optimization simulations of human locomotion by driving a neuromusculoskeletal model toward in vivo measurements of body-segmental kinematics and ground reaction forces. Gait data were recorded from 5 healthy participants who walked at their preferred speeds and ran at 2m/s. Participant-specific data-tracking dynamic optimization solutions were generated for one stride cycle using direct collocation in tandem with an OpenSim-MATLAB interface. The body was represented as a 12-segment, 21-degree-of-freedom skeleton actuated by 66 muscle-tendon units. Foot-ground interaction was simulated using six contact spheres under each foot. The dynamic optimization problem was to find the set of muscle excitations needed to reproduce 3D measurements of body-segmental motions and ground reaction forces while minimizing the time integral of muscle activations squared. Direct collocation took on average 2.7±1.0h and 2.2±1.6h of CPU time, respectively, to solve the optimization problems for walking and running. Model-computed kinematics and foot-ground forces were in good agreement with corresponding experimental data while the calculated muscle excitation patterns were consistent with measured EMG activity. The results demonstrate the feasibility of implementing direct collocation on a detailed neuromusculoskeletal model with foot-ground contact to accurately and efficiently generate 3D data-tracking dynamic optimization simulations of human locomotion. The proposed method offers a viable tool for creating feasible initial guesses needed to perform predictive simulations of movement using dynamic optimization theory. The source code for implementing the model and computational algorithm may be downloaded at http://simtk.org/home/datatracking.

Muscle Synergies Facilitate Computational Prediction of Subject-Specific Walking Motions

Researchers have explored a variety of neurorehabilitation approaches to restore normal walking function following a stroke. However, there is currently no objective means for prescribing and implementing treatments that are likely to maximize recovery of walking function for any particular patient. As a first step toward optimizing neurorehabilitation effectiveness, this study develops and evaluates a patient-specific synergy-controlled neuromusculoskeletal simulation framework that can predict walking motions for an individual post-stroke. The main question we addressed was whether driving a subject-specific neuromusculoskeletal model with muscle synergy controls (5 per leg) facilitates generation of accurate walking predictions compared to a model driven by muscle activation controls (35 per leg) or joint torque controls (5 per leg). To explore this question, we developed a subject-specific neuromusculoskeletal model of a single high-functioning hemiparetic subject using instrumented treadmill walking data collected at the subject's self-selected speed of 0.5 m/s. The model included subject-specific representations of lower-body kinematic structure, foot-ground contact behavior, electromyography-driven muscle force generation, and neural control limitations and remaining capabilities. Using direct collocation optimal control and the subject-specific model, we evaluated the ability of the three control approaches to predict the subject's walking kinematics and kinetics at two speeds (0.5 and 0.8 m/s) for which experimental data were available from the subject. We also evaluated whether synergy controls could predict a physically realistic gait period at one speed (1.1 m/s) for which no experimental data were available. All three control approaches predicted the subject's walking kinematics and kinetics (including ground reaction forces) well for the model calibration speed of 0.5 m/s. However, only activation and synergy controls could predict the subject's walking kinematics and kinetics well for the faster non-calibration speed of 0.8 m/s, with synergy controls predicting the new gait period the most accurately. When used to predict how the subject would walk at 1.1 m/s, synergy controls predicted a gait period close to that estimated from the linear relationship between gait speed and stride length. These findings suggest that our neuromusculoskeletal simulation framework may be able to bridge the gap between patient-specific muscle synergy information and resulting functional capabilities and limitations.