Paired nonlinear behavior of active and passive joint torques associated with preparation for walk-to-run gait transition

This study investigated the lower extremity torque's active and passive features during the walk-to-run gait transition with continuously increased walking speed. Fourteen volunteers participated in the experiment. Kinematic and kinetic data were collected synchronously. Five strides leading up the gait transition were examined. Peaks of the passive (e.g., contact) and active (e.g., generalized muscle torques), along with net joint torque, and time to peak torques exhibited significant differences at the last stride before gait transition, compared to the first four strides, at the ankle, knee, and hip joints, respectively. Selected peak joint active and passive torques showed significant and opposite trends at critical events within a stride cycle: such ankle joint right after heel-contact, knee joint during weight acceptance, and both hip and knee joints right before toe-off. The magnitude and the corresponding time to active and passive peak torque changed in a nonlinear pattern before the transition from walk to run. The lower extremity segment-interaction during gait transition appeared to be an active reorganization exemplified by the interaction between the lower extremity's active and passive torque components.     

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.     

The scaling or ontogeny of human gait kinetics and walk-run transition: The implications of work vs. peak power minimization

A simple model is developed to find vertical force profiles and stance durations that minimize either limb mechanical work or peak power demands during bipedal locomotion. The model predicts that work minimization is achieved with a symmetrical vertical force profile, consistent with previous models and observations of adult humans, and data for 487 participants (predominantly 11–18 years old) required to walk at a range of speeds at a Science Fair. Work minimization also predicts the discrete walk-run transition, familiar for adult humans. In contrast, modeled peak limb mechanical power demands are minimized with an early skew in vertical ground reaction force that increases with speed, and stance durations that decrease steadily with speed across the work minimizing walk-run transition speed. The peak power minimization model therefore predicts a continuous walk-run gait transition that is quantitatively consistent with measurements of younger children (1.1–4.7 years) required to locomote at a range of speeds but free to select their own gaits.     

Spatio-temporal gait characteristics during transitions from trot to canter in horses

Gaits can be defined based upon specific interlimb coordination patterns characteristic of a limited range of speeds, with one or more defining variables changing discontinuously at a transition. With changing speed, horses perform a repertoire of gaits (walk, trot, canter and gallop), with transitions between them. Knowledge of the series of kinematic events necessary to realize a gait is essential for understanding the proximate mechanisms as well as the control underlying gait transitions. We studied the kinematics of the actual transition from trot to canter in miniature horses. The kinematics were characterized at three different levels: the whole-body level, the spatio-temporal level of the foot falls and the level of basic limb kinematics. This concept represents a hierarchy: the horse's center of mass (COM) moves forward by means of the coordinated action of the limbs and changes in the latter are the result of alterations in the basic limb kinematics. Early and short placement of the fore limb was observed before the dissociation of the footfalls of one of the diagonal limb pairs when entering the canter. Dissociation coincided with increased amplitude and wavelength of the oscillations of the trunk in the sagittal plane. The increased amplitude cannot be explained solely by the passive effects of acceleration or by neck and head movements which are inconsistent with the timing of the transition. We propose that the transition is initiated by the fore limb followed by subsequent changes in the hind limbs in a series of kinematic events that take about 2.5 strides to complete.     

High day-to-day repeatability of lower extremity muscle activation patterns and joint biomechanics of dual-belt treadmill gait: A reliability study in healthy young adults

PurposeThe reliability of lower extremity muscle activation patterns has not been clearly studied in a dual-belt instrumented treadmill environment. The primary study objective was to quantify the day-to-day reliability of quadriceps, hamstrings, gastrocnemius and gluteus medius activation patterns in healthy young adult gait. Secondarily, the reliability of spatiotemporal, and knee/hip motion and moment-based gait outcomes was assessed.Scope: 20 young adults were recruited and tested on two separate days. Using standardized procedures, participants were prepared for surface electromyography and lower extremity motion capture. All individuals walked on a dual-belt instrumented treadmill while muscle activation, segment motions and ground reaction forces were recorded. Sagittal plane motion and net external sagittal and frontal plane moments were calculated. Discrete biomechanical and muscle activation measures were calculated, and non-negative matrix factorization extracted amplitude and temporal muscle activation features. Intraclass Correlation Coefficients, Standard Error of Measurement and Minimum Detectable Change were calculated.ConclusionsHigh to excellent Intraclass correlation coefficients were found between visits for most primary and secondary outcomes. The absolute and relative reliability, including Minimum Detectable Change values, provided in this study support the use of dual-belt instrumented treadmill walking as an acceptable medium to collect biomechanical and lower extremity EMG outcomes for future studies.     

Walk-, run- and gallop-like gait patterns in human sideways locomotion

The purpose of this study is to examine the characteristics of gait patterns in human preferred sideways locomotion at increasing speeds. Fifteen healthy young males were asked to step sideways on a treadmill at various speeds of 1.3–6.1 km/h. The times of foot contact and take-off were analyzed. Three gait patterns were observed. At slow speeds, all of the subjects performed a walk-like pattern. When the treadmill speed exceeded approximately 3.5 km/h, the subjects preferred gait patterns with a flight phase. Most of the subjects performed an asymmetric gait pattern that was similar to a forward gallop, whereas only two out of fifteen subjects performed a run-like gait pattern. Because the left and right legs are positioned along the movement direction, it might be more efficient to divide roles between the leading and trailing limbs at high speeds: the leading limb functions to produces breaking and vertical force, and the trailing limb mainly absorbs the impact of foot contact and generates propulsive forces.     

The influence of step-down technique on lower extremity mechanics during curb descent

When stepping down from a curb, individuals typically make initial ground contact with either their rearfoot or forefoot. The purpose of this study was to compare vertical ground reaction forces, lower extremity mechanics, and intra-limb work distribution when individuals adopt a rearfoot technique vs. a forefoot technique, during simulated curb descent. Sixteen subjects stepped down from a platform with both a rearfoot and a forefoot technique. Vertical ground reaction forces and sagittal plane joint kinematics and kinetics were examined for the lead limb during the step-down task. Paired t-tests were used for comparison. Subjects demonstrated greater ankle joint power and negative work, and less hip joint power and negative work, with the forefoot technique vs. the rearfoot technique. Total lower extremity negative work was greater for the forefoot technique vs. the rearfoot technique. The percent contribution to the total negative work was greater for the ankle joint, and less for the hip and knee joints, with the forefoot technique vs. the rearfoot technique. The results of this study may provide insight into how curb descent technique can be modified to alter lower extremity loading.     

The effect of a novel gait retraining device on lower limb kinematics and muscle activation in healthy adults

The Re-Link Trainer (RLT) is a modified walking frame with a linkage system designed to apply a non-individualized kinematic constraint to normalize gait trajectory of the left limb. The premise behind the RLT is that a user’s lower limb is constrained into a physiologically normal gait pattern, ideally generating symmetry across gait cycle parameters and kinematics. This pilot study investigated adaptations in the natural gait pattern of healthy adults when using the RLT compared to normal overground walking. Bilateral lower limb kinematic and electromyography data were collected while participants walked overground at a self-selected speed, followed by walking in the RLT. A series of 2-way analyses of variance examined between-limb and between-condition differences. Peak hip extension and knee flexion were reduced bilaterally when walking in the RLT. Left peak hip extension occurred earlier in the gait cycle when using the RLT, but later for the right limb. Peak hip flexion was significantly increased and occurred earlier for the constrained limb, while peak plantarflexion was significantly reduced. Peak knee flexion and plantarflexion in the right limb occurred later when using the RLT. Significant bilateral reductions in peak electromyography amplitude were evident when walking in the RLT, along with a significant shift in when peak muscle activity was occurring. These findings suggest that the RLT does impose a significant constraint, but generates asymmetries in lower limb kinematics and muscle activity patterns. The large interindividual variation suggests users may utilize differing motor strategies to adapt their gait pattern to the imposed constraint.     

A general model for estimating lower extremity inertial properties of individuals with transtibial amputation

Lower extremity joint moment magnitudes during swing are dependent on the inertial properties of the prosthesis and residual limb of individuals with transtibial amputation (TTA). Often, intact limb inertial properties (INTACT) are used for prosthetic limb values in an inverse dynamics model even though these values overestimate the amputated limb’s inertial properties. The purpose of this study was to use subject-specific (SPECIFIC) measures of prosthesis inertial properties to generate a general model (GENERAL) for estimating TTA prosthesis inertial properties. Subject-specific mass, center of mass, and moment of inertia were determined for the shank and foot segments of the prosthesis (n = 11) using an oscillation technique and reaction board. The GENERAL model was derived from the means of the SPECIFIC model. Mass and segment lengths are required GENERAL model inputs. Comparisons of segment inertial properties and joint moments during walking were made using three inertial models (unique sample; n = 9): (1) SPECIFIC, (2) GENERAL, and (3) INTACT. Prosthetic shank inertial properties were significantly smaller with the SPECIFIC and GENERAL model than the INTACT model, but the SPECIFIC and GENERAL model did not statistically differ. Peak knee and hip joint moments during swing were significantly smaller for the SPECIFIC and GENERAL model compared with the INTACT model and were not significantly different between SPECIFIC and GENERAL models. When subject-specific measures are unavailable, using the GENERAL model produces a better estimate of prosthetic side inertial properties resulting in more accurate joint moment measurements for individuals with TTA than the INTACT model.     

The effect on conventional gait model kinematics and kinetics of hip joint centre equations in adult healthy gait

The Conventional Gait Model (CGM) needs to benefit from large investigations on localization of the hip joint centre (HJC). Incorrect positions from the native equations were demonstrated (Sangeux et al., 2014; Harrington et al., 2007). More accurate equations were proposed but their impact on kinematics and kinetic CGM outputs was never evaluated. This short communication aims at examining if adoption of new HJC equations would alter standard CGM outputs. Sixteen able bodied participants underwent a full 3-D optoelectronic gait analysis followed by a 3-D ultrasound localization of their hips. Data were processed through the open source python package pyCGM2 replicating kinematic and kinetic processing of the native CGM. Compared with 3D ultrasound location, Hara equations improved the accuracy of sagittal plane kinematics (0.6°) and kinetics (0.02 N m kg⁻¹) for the hip. The worst case participant exhibited Harrington’s equations reached a deviation of 3° for the sagittal kinematics. In the coronal plane, Hara and Harrington equations presented similar differences (1°) for the hip whilst Davis equations had the largest deviation for hip abduction (2.7°) and hip abductor moment (0.10 N m kg⁻¹).Both Harrington and Hara equations improved the CGM location of the HJC. Hara equations improved results in the sagittal plane, plus utilise a single anthropometrics measurement, leg length, that may be more robust. However, neither set of equations had significant effect on kinematics. We reported some effects on kinetics, particularly in the coronal plane, which warrant caution in interpreting outputs using different sets of equations.     

Ground reaction forces and lower-limb joint kinetics of turning gait in typically developing children

Turning is a common locomotor task essential to daily activity; however, very little is known about the forces and moments responsible for the kinematic adaptations occurring relative to straight-line gait in typically developing children. Thus, the aims of this study were to analyse ground reaction forces (GRFs), ground reaction free vertical torque (T_Z), and the lower-limb joint kinetics of 90° outside (step) and inside (spin) limb turns. Step, spin, and straight walking trials from fifty-four typically developing children were analysed. All children were fit with the Plug-in Gait and Oxford Foot Model marker sets while walking over force plates embedded in the walkway. Net internal joint moments and power were computed via a standard inverse dynamics approach. All dependent variables were statistically analysed over the entire curves using the mean difference 95% bootstrap confidence band approach. GRFs were directed medially for step turns and laterally for spin turns during the turning phase. Directions were reversed and magnitudes decreased during the approach phase. Step turns showed reduced ankle power generation, while spin turns showed large T_Z. Both strategies required large knee and hip coronal and transverse plane moments during swing. These kinetic differences highlight adaptations required to maintain stability and reorient the body towards the new walking direction during turning. From a clinical perspective, turning gait may better reveal weaknesses and motor control deficits than straight walking in pathological populations, such as children with cerebral palsy, and could potentially be implemented in standard gait analysis sessions.     

Diurnal variation in gait characteristics and transition speed

The aim of this study was to investigate the effect of time-of-day on Preferred Transition Speed (PTS) and spatiotemporal organization of walking and running movements. Twelve active male subjects participated in the study (age: 27.2?±?4.9 years; height: 177.9?±?5.4?cm; body mass: 75.9?±?5.86?kg). First, PTS was determined at 08:00?h and 18:00?h. The mean of the two PTS recorded at the two times-of-day tested was used as a reference (PTS_m). Then, subjects were asked to walk and run on a treadmill at three imposed speeds (PTS_m, PTS_m?+?0.3?m.s^?1, and PTS_m???0.3?m.s^?1) at 08:00?h and 18:00?h. Mean stride length, temporal stride, spatial stride variability, and temporal stride variability were used for gait analysis. The PTS observed at 08:00?h (2.10?±?0.17?m.s^?1) tends to be lower (p?=?0.077) than that recorded at 18:00?h (2.14?±?0.19?m.s^?1). Stride lengths recorded while walking (p?=?0.038) and running (p?=?0.041) were shorter at 08:00?h than 18:00?h. No time-of-day effect was observed for stride frequency during walking and running trials. When walking, spatial stride variability (p?=?0.020) and temporal stride variability (p?=?0.028) were lower at 08:00?h than at 18:00?h. When running, no diurnal variation of spatial stride variability or temporal stride variability was detected.     

The influence of estimated body segment parameters on predicted joint kinetics during diplegic cerebral palsy gait

Inverse Dynamic calculations are routinely used in joint moment and power estimates during gait with anthropometric data often taken from published sources. Many biomechanical analyses have highlighted the need to obtain subject-specific anthropometric data (e.g. Mass, Centre of Mass, Moments of Inertia) yet the types of imaging techniques required to achieve this are not always available in the clinical setting. Differences in anthropometric sets have been shown to affect the reactive force and moment calculations in normal subjects but the effect on a paediatric diplegic cerebral palsy group has not been investigated. The aim of this study was to investigate the effect of using different anthropometric sets on predicted sagittal plane moments during normal and diplegic cerebral palsy gait. Three published anthropometric sets were applied to the reactive force and moment calculations of 14 Cerebral Palsy and 14 Control subjects. Statistically significant differences were found when comparing the different anthropometric sets but variability in the resulting sagittal plane moment calculations between sets was low (0.01–0.07 Nm/kg). In addition, the GDI-Kinetic, used as an outcome variable to assess whether differences were clinically meaningful, indicated no clinically meaningful difference between sets. The results suggest that the effects of using different anthropometric sets on the kinetic profiles of normal and diplegic cerebral palsy subjects are clinically insignificant.     

Kalman smoothing improves the estimation of joint kinematics and kinetics in marker-based human gait analysis

We developed a Kalman smoothing algorithm to improve estimates of joint kinematics from measured marker trajectories during motion analysis. Kalman smoothing estimates are based on complete marker trajectories. This is an improvement over other techniques, such as the global optimisation method (GOM), Kalman filtering, and local marker estimation (LME), where the estimate at each time instant is only based on part of the marker trajectories. We applied GOM, Kalman filtering, LME, and Kalman smoothing to marker trajectories from both simulated and experimental gait motion, to estimate the joint kinematics of a ten segment biomechanical model, with 21 degrees of freedom. Three simulated marker trajectories were studied: without errors, with instrumental errors, and with soft tissue artefacts (STA). Two modelling errors were studied: increased thigh length and hip centre dislocation. We calculated estimation errors from the known joint kinematics in the simulation study. Compared with other techniques, Kalman smoothing reduced the estimation errors for the joint positions, by more than 50% for the simulated marker trajectories without errors and with instrumental errors. Compared with GOM, Kalman smoothing reduced the estimation errors for the joint moments by more than 35%. Compared with Kalman filtering and LME, Kalman smoothing reduced the estimation errors for the joint accelerations by at least 50%. Our simulation results show that the use of Kalman smoothing substantially improves the estimates of joint kinematics and kinetics compared with previously proposed techniques (GOM, Kalman filtering, and LME) for both simulated, with and without modelling errors, and experimentally measured gait motion.     

Influence of gait speed on the control of mediolateral dynamic stability during gait initiation

This study investigated the influence of gait speed on the control of mediolateral dynamic stability during gait initiation. Thirteen healthy young adults initiated gait at three self-selected speeds: Slow, Normal and Fast. The results indicated that the duration of anticipatory postural adjustments (APA) decreased from Slow to Fast, i.e. the time allocated to propel the centre of mass (COM) towards the stance-leg side was shortened. Likely as an attempt at compensation, the peak of the anticipatory centre of pressure (COP) shift increased. However, COP compensation was not fully efficient since the results indicated that the mediolateral COM shift towards the stance-leg side at swing foot-off decreased with gait speed. Consequently, the COM shift towards the swing-leg side at swing heel-contact increased from Slow to Fast, indicating that the mediolateral COM fall during step execution increased as gait speed rose. However, this increased COM fall was compensated by greater step width so that the margin of stability (the distance between the base-of-support boundary and the mediolateral component of the “extrapolated centre of mass”) at heel-contact remained unchanged across the speed conditions. Furthermore, a positive correlation between the mediolateral extrapolated COM position at heel-contact and step width was found, indicating that the greater the mediolateral COM fall, the greater the step width. Globally, these results suggest that mediolateral APA and step width are modulated with gait speed so as to maintain equivalent mediolateral dynamical stability at the time of swing heel-contact.     

The influence of knee alignment on lower extremity kinetics during squats

The squat is an assessment of lower extremity alignment during movement, however there is little information regarding altered joint kinetics during poorly performed squats. The purpose of this study was to examine changes in joint kinetics and power from altered knee alignment during a squat. Thirty participants completed squats while displacing the knee medially, anteriorly, and with neutral alignment (control). Sagittal and frontal plane torques at the ankle, knee, and hip were altered in the descending and ascending phase of the squat in both the medial and anterior malaligned squat compared to the control squat. Ankle and trunk power increased and hip power decreased in the medial malaligned squat compared to the control squat. Ankle, knee, and trunk power increased and hip power decreased in the anterior malaligned squat compared to the control squat. Changes in joint torques and power during malaligned squats suggest that altered knee alignment increases ankle and trunk involvement to execute the movement. Increased anterior knee excursion during squatting may also lead to persistent altered loading of the ankle and knee. Sports medicine professionals using the squat for quadriceps strengthening must consider knee alignment to reduce ankle and trunk involvement during the movement.     

A lower extremity model for muscle-driven simulation of activity using explicit finite element modeling

A key strength of computational modeling is that it can provide estimates of muscle, ligament, and joint loads, stresses, and strains through non-invasive means. However, simulations that can predict the forces in the muscles during activity while maintaining sufficient complexity to realistically represent the muscles and joint structures can be computationally challenging. For this reason, the current state of the art is to apply separate rigid-body dynamic and finite-element (FE) analyses in series. However, the use of two or more disconnected models often fails to capture key interactions between the joint-level and whole-body scales. Single framework MSFE models have the potential to overcome the limitations associated with disconnected models in series. The objectives of the current study were to create a multi-scale FE model of the human lower extremity that combines optimization, dynamic muscle modeling, and structural FE analysis in a single framework and to apply this framework to evaluate the mechanics of healthy knee specimens during two activities. Two subject-specific FE models (Model 1, Model 2) of the lower extremity were developed in ABAQUS/Explicit including detailed representations of the muscles. Muscle forces, knee joint loading, and articular contact were calculated for two activities using an inverse dynamics approach and static optimization. Quadriceps muscle forces peaked at the onset of chair rise (2174 N, 1962 N) and in early stance phase (510 N, 525 N), while gait saw peak forces in the hamstrings (851 N, 868 N) in midstance. Joint forces were similar in magnitude to available telemetric patient data. This study demonstrates the feasibility of detailed quasi-static, muscle-driven simulations in an FE framework.