Friedl De Groote

The flexion synergy, mother of all synergies and father of new models of gait

Recently there has been a growing interest in the modular organization of leg movements, in particular those related to locomotion. One of the basic modules involves the flexion of the leg during swing and it was shown that this module is already present in neonates (Dominici et al., 2011). In this paper, we question how these finding build upon the original work by Sherrington, who proposed that the flexor reflex is the basic building block of flexion during swing phase. Similarly, the relation between the flexor reflex and the withdrawal reflex modules of Schouenborg and Weng (1994) will be discussed. It will be argued that there is large overlap between these notions on modules and the older concepts of reflexes. In addition, it will be shown that there is a great flexibility in the expression of some of these modules during gait, thereby allowing for a phase-dependent modulation of the appropriate responses. In particular, the end of the stance phase is a period when the flexor synergy is facilitated. It is proposed that this is linked to the activation of circuitry that is responsible for the generation of locomotor patterns (CPG, “central pattern generator”). More specifically, it is suggested that the responses in that period relate to the activation of a flexor burst generator. The latter structure forms the core of a new asymmetric model of the CPG. This activation is controlled by afferent input (facilitation by a broad range of afferents, suppression by load afferent input). Meanwhile, many of these physiologic features have found their way in the control of very flexible walking bipedal robots.

Evaluation of Direct Collocation Optimal Control Problem Formulations for Solving the Muscle Redundancy Problem.

Estimation of muscle forces during motion involves solving an indeterminate problem (more unknown muscle forces than joint moment constraints), frequently via optimization methods. When the dynamics of muscle activation and contraction are modeled for consistency with muscle physiology, the resulting optimization problem is dynamic and challenging to solve. This study sought to identify a robust and computationally efficient formulation for solving these dynamic optimization problems using direct collocation optimal control methods. Four problem formulations were investigated for walking based on both a two and three dimensional model. Formulations differed in the use of either an explicit or implicit representation of contraction dynamics with either muscle length or tendon force as a state variable. The implicit representations introduced additional controls defined as the time derivatives of the states, allowing the nonlinear equations describing contraction dynamics to be imposed as algebraic path constraints, simplifying their evaluation. Problem formulation affected computational speed and robustness to the initial guess. The formulation that used explicit contraction dynamics with muscle length as a state failed to converge in most cases. In contrast, the two formulations that used implicit contraction dynamics converged to an optimal solution in all cases for all initial guesses, with tendon force as a state generally being the fastest. Future work should focus on comparing the present approach to other approaches for computing muscle forces. The present approach lacks some of the major limitations of established methods such as static optimization and computed muscle control while remaining computationally efficient.

Task constraints and minimization of muscle effort result in a small number of muscle synergies during gait

Finding muscle activity generating a given motion is a redundant problem, since there are many more muscles than degrees of freedom. The control strategies determining muscle recruitment from a redundant set are still poorly understood. One theory of motor control suggests that motion is produced through activating a small number of muscle synergies, i.e., muscle groups that are activated in a fixed ratio by a single input signal. Because of the reduced number of input signals, synergy-based control is low dimensional. But a major criticism on the theory of synergy-based control of muscles is that muscle synergies might reflect task constraints rather than a neural control strategy. Another theory of motor control suggests that muscles are recruited by optimizing performance. Optimization of performance has been widely used to calculate muscle recruitment underlying a given motion while assuming independent recruitment of muscles. If synergies indeed determine muscle recruitment underlying a given motion, optimization approaches that do not model synergy-based control could result in muscle activations that do not show the synergistic muscle action observed through electromyography (EMG). If, however, synergistic muscle action results from performance optimization and task constraints (joint kinematics and external forces), such optimization approaches are expected to result in low-dimensional synergistic muscle activations that are similar to EMG-based synergies. We calculated muscle recruitment underlying experimentally measured gait patterns by optimizing performance assuming independent recruitment of muscles. We found that the muscle activations calculated without any reference to synergies can be accurately explained by on average four synergies. These synergies are similar to EMG-based synergies. We therefore conclude that task constraints and performance optimization explain synergistic muscle recruitment from a redundant set of muscles.

Gait alterations to effectively reduce hip contact forces

Patients with hip pathology present alterations in gait which have an effect on joint moments and loading. In knee osteoarthritic patients, the relation between medial knee contact forces and the knee adduction moment are currently being exploited to define gait retraining strategies to effectively reduce pain and disease progression. However, the relation between hip contact forces and joint moments has not been clearly established. Therefore, this study aims to investigate the effect of changes in hip and pelvis kinematics during gait on internal hip moments and contact forces which is calculated using muscle driven simulations. The results showed that frontal plane kinetics have the largest effect on hip contact forces. Given the high correlation between the change in hip adduction moment and contact force at initial stance (R2 = 0.87), this parameter can be used to alter kinematics and predict changes in contact force. At terminal stance the hip adduction and flexion moment can be used to predict changes in contact force (R2 = 0.76). Therefore, gait training that focuses on decreasing hip adduction moments, a wide base gait pattern, has the largest potential to reduce hip contact forces.

The effect of muscle weakness on the capability gap during gross motor function: a simulation study supporting design criteria for exoskeletons of the lower limb

BackgroundEnabling persons with functional weaknesses to perform activities of daily living (ADL) is one of the main challenges for the aging society. Powered orthoses, or exoskeletons, have the potential to support ADL while promoting active participation of the user. For this purpose, assistive devices should be designed and controlled to deliver assistance as needed (AAN). This means that the level of assistance should bridge the capability gap, i.e. the gap between the capabilities of the subjects and the task requirements. However, currently the actuators of exoskeletons are mainly designed using inverse dynamics (ID) based calculations of joint moments. The goal of the present study is to calculate the capability gap for the lower limb during ADL when muscle weakness is present, which is needed for appropriate selection of actuators to be integrated in exoskeletons.MethodsA musculoskeletal model (MM) is used to calculate the joint kinematics, joint kinetics and muscle forces of eight healthy subjects during ADL (gait, sit-to-stand, stand-to-sit, stair ascent, stair descent). Muscle weakness was imposed to the MM by a stepwise decrease in maximal isometric force imposed to all muscles. Muscle forces were calculated using static optimization. In order to compensate for muscle weakness, ideal moment actuators that represent the motors of an exoskeleton in the simulation were added to deliver AAN required to perform the task.ResultsThe ID approach overestimates the required assistance since it relies solely on the demands of the task, whereas the AAN approach incorporates the capabilities of the subject. Furthermore, the ID approach delivers continuous support whereas the AAN approach targets the period where a capability gap occurs. The level of muscle weakness for which the external demands imposed by ADL can no longer be met by active muscle force production, is respectively 40%, 70%, 80% and 30%.ConclusionsThe present workflow allows estimating the AAN during ADL for different levels of muscle weakness, which can be used in the mechatronic design and control of powered exoskeletons. The AAN approach is a more physiological approach than the ID approach, since the MM accounts for the subject-specific capabilities of the user.

Altering length and velocity feedback during a neuro-musculoskeletal simulation of normal gait contributes to hemiparetic gait characteristics.

BackgroundSpasticity is an important complication after stroke, especially in the anti-gravity muscles, i.e. lower limb extensors. However the contribution of hyperexcitable muscle spindle reflex loops to gait impairments after stroke is often disputed. In this study a neuro-musculoskeletal model was developed to investigate the contribution of an increased length and velocity feedback and altered reflex modulation patterns to hemiparetic gait deficits.MethodsA musculoskeletal model was extended with a muscle spindle model providing real-time length and velocity feedback of gastrocnemius, soleus, vasti and rectus femoris during a forward dynamic simulation (neural control model). By using a healthy subject’s base muscle excitations, in combination with increased feedback gains and altered reflex modulation patterns, the effect on kinematics was simulated. A foot-ground contact model was added to account for the interaction effect between the changed kinematics and the ground. The qualitative effect i.e. the directional effect and the specific gait phases where the effect is present, on the joint kinematics was then compared with hemiparetic gait deviations reported in the literature.ResultsOur results show that increased feedback in combination with altered reflex modulation patterns of soleus, vasti and rectus femoris muscle can contribute to excessive ankle plantarflexion/inadequate dorsiflexion, knee hyperextension/inadequate flexion and increased hip extension/inadequate flexion during dedicated gait cycle phases. Increased feedback of gastrocnemius can also contribute to excessive plantarflexion/inadequate dorsiflexion, however in combination with excessive knee and hip flexion. Increased length/velocity feedback can therefore contribute to two types of gait deviations, which are both in accordance with previously reported gait deviations in hemiparetic patients. Furthermore altered modulation patterns, in particular the reduced suppression of the muscle spindle feedback during swing, can contribute largely to an increased plantarflexion and knee extension during the swing phase and consequently to hampered toe clearance.ConclusionsOur results support the idea that hyperexcitability of length and velocity feedback pathways, especially in combination with altered reflex modulation patterns, can contribute to deviations in hemiparetic gait. Surprisingly, our results showed only subtle temporal differences between length and velocity feedback. Therefore, we cannot attribute the effects seen in kinematics to one specific type of feedback.

The effect of perturbing body segment parameters on calculated joint moments and muscle forces during gait

This study examined the effect of body segment parameter (BSP) perturbations on joint moments calculated using an inverse dynamics procedure and muscle forces calculated using computed muscle control (CMC) during gait. BSP (i.e. segment mass, center of mass location (com) and inertia tensor) of the left thigh, shank and foot of a scaled musculoskeletal model were perturbed. These perturbations started from their nominal value and were adjusted to ±40% in steps of 10%, for both individual as well as combined perturbations in BSP. For all perturbations, an inverse dynamics procedure calculated the ankle, knee and hip moments based on an identical inverse kinematics solution. Furthermore, the effect of applying a residual reduction algorithm (RRA) was investigated. Muscle excitations and resulting muscle forces were calculated using CMC. The results show only a limited effect of an individual parameter perturbation on the calculated moments, where the largest effect is found when perturbing the shank com (MS(com,shank), the ratio of absolute difference in torque and relative parameter perturbation, is maximally -7.81 N m for hip flexion moment). The additional influence of perturbing two parameters simultaneously is small (MS(mass+com,thigh) is maximally 15.2 N m for hip flexion moment). RRA made small changes to the model to increase the dynamic consistency of the simulation (after RRA MS(com,shank) is maximally 5.01 N m). CMC results show large differences in muscle forces when BSP are perturbed. These result from the underlying forward integration of the dynamic equations.

How gravity and muscle action control mediolateral center of mass excursion during slow walking: a simulation study.

Maintaining mediolateral (ML) balance is very important to prevent falling during walking, especially at very slow speeds. The effect of walking speed on support and propulsion of the center of mass (COM) has been focus of previous studies. However, the influence of speed on ML COM control and the associated coupling with sagittal plane control remains unclear. Simulations of walking at very slow and normal speeds were generated for twelve healthy subjects. Our results show that gluteus medius (GMED) contributions to ML stability decrease, while its contributions to sagittal plane accelerations increase during very slow compared to normal walking. Simultaneously the destabilizing influence of gravity increases in ML direction at a very slow walking speed. This emphasizes the need for a tight balance between gravity and gluteus medius action to ensure ML stability. When walking speed increases, GMED has a unique role in controlling ML acceleration and therefore stabilizing ML COM excursion. Contributions of other muscles decrease in all directions during very slow speed. Increased contributions of these muscles are therefore required to provide for both stability and propulsion when walking speed increases.

Sensitivity of predicted muscle forces during gait to anatomical variability in musculotendon geometry.

Scaled generic musculoskeletal models are commonly used to drive dynamic simulations of motions. It is however, acknowledged that not accounting for variability in musculoskeletal geometry and musculotendon parameters may confound the simulation results, even when analysing control subjects. This study documents the three-dimensional anatomical variability of musculotendon origins and insertions of 33 lower limb muscles determined based on magnetic resonance imaging in six subjects. This anatomical variability was compared to the musculotendon point location in a generic musculoskeletal model. Furthermore, the sensitivity of muscle forces during gait, calculated using static optimization, to perturbations of the musculotendon point location was analyzed with a generic model. More specific, a probabilistic approach was used: for each analyzed musculotendon point, the three-dimensional location was re-sampled with a uniform Latin hypercube method within the anatomical variability and the static optimization problem was then re-solved for all perturbations. We found that musculotendon point locations in the generic model showed only variable correspondences with the anatomical variability. The anatomical variability of musculotendon point location did affect the calculated muscle forces: muscles most sensitive to perturbations within the anatomical variability are iliacus and psoas. Perturbation of the gluteus medius anterior, iliacus and psoas induces the largest concomitant changes in muscle forces of the unperturbed muscles. Therefore, when creating subject-specific musculoskeletal models, these attachment points should be defined accurately. In addition, the size of the anatomical variability of the musculotendon point location was not related to the sensitivity of the calculated muscle forces.

A new method for estimating subject-specific muscle-tendon parameters of the knee joint actuators: a simulation study.

A new method for the estimation of subject-specific muscle-tendon parameters of the knee actuators based on dynamometry experiments is presented. The algorithm aims at estimating the tendon slack length and the optimal muscle fiber length by minimizing the difference between experimentally reproduced and model-based joint moments. The key innovative features are as follows: (i) the inclusion of a priori physiological knowledge to define a physiologically feasible set, the hot start for the optimization, and constraints for the optimization and (ii) the introduction of a new (affine) transformation of the muscle-tendon parameters, which greatly improves the numerical condition of the optimization. The influence of the initial guess and of measurement noise was studied in a simulation environment, and the performance was compared with that of the method presented earlier by Garner and Pandy for the upper limb. The tendon slack length was estimated for 97.5/63% (extensors/flexors) of all initial guesses within 2% of the ground truth. The optimal fiber length was estimated for 89/90% (extensors/flexors) of all initial guesses within 2% of the ground truth. When 10 Nm measurement noise was added, the mean value of the estimated tendon slack length deviated at most 1.9/1.6% (extensors/flexors) from the ground truth whereas the standard deviations were at most 5.1/3.9%. The mean value of the estimated optimal fiber length deviated at most 4.3/3.0% (extensors/flexors) from the ground truth whereas the standard deviations were at most 10.2/15.5%. In comparison, mean values resulting from the method of Garner and Pandy deviated up to 181% ( ± 123%) and 119% ( ± 30%) from the ground truth for, respectively, optimal fiber length and tendon slack length of rectus femoris. We concluded that the presented method had a low dependency on the initial guess and outperformed the method of Garner and Pandy in terms of accuracy by at least one order of magnitude when parameters were estimated from noisy data. The improvements open new perspectives for subject-specific modelling of muscles and tendons, which is beneficial for the accuracy of human motion simulations.