Quantitative assessment of gait determinants during single stance via a three-dimensional model--Part 1. Normal gait

In this two-part paper, a variety of three-dimensional, dynamical models are constructed for simulating the single support phases of normal and pathological human gait. A major objective of this work is to quantify the influence of individual gait determinants on the ground reaction forces generated during normal, level walking. To this end, Part 1 presents a three-dimensional, seven degree-of-freedom model incorporating five of the six fundamental determinants of gait. On the basis of crude muscle-force and/or joint-moment trajectories, body-segmental motions and ground reaction forces are synthesized open loop. Through a quantitative comparison with experimental gait data, the model's predictions are evaluated. Our simulation results suggest that pelvic list is not as dominant a dynamical determinant as either stance knee flexion-extension or foot and knee interaction. Transverse pelvic rotation, however, makes an important contribution by limiting the magnitude of the horizontal ground reaction prior to opposite heel-strike.     

A Phase-Dependent Hypothesis for Locomotor Functions of Human Foot Complex

The human foot is a very complex structure comprising numerous bones, muscles, ligaments and synovial joints. As the only component in contact with the ground, the foot complex delivers a variety of biomechanical functions during human locomotion, e.g. body support and propulsion, stability maintenance and impact absorption. These need the human foot to be rigid and damped to transmit ground reaction forces to the upper body and maintain body stability, and also to be compliant and resilient to moderate risky impacts and save energy. How does the human foot achieve these apparent conflicting functions？ In this study, we propose a phase-dependent hypothesis for the overall locomotor functions of the human foot complex based on in-vivo measurements of human natural gait and simulation results of a mathematical foot model. We propse that foot functions are highly dependent on gait phase, which is a major characteristics of human locomotion. In early stance just after heel strike, the foot mainly works as a shock absorber by moderating high impacts using the viscouselastic heel pad in both vertical and horizontal directions. In mid-stance phase （-80% of stance phase）, the foot complex can be considered as a springy rocker, reserving external mechanical work using the foot arch whilst moving ground contact point forward along a curved path to maintain body stability. In late stance after heel off, the foot complex mainly serves as a force modulator like a gear box, modulating effective mechanical advantages of ankle plantiflexor muscles using metatarsal-phalangeal joints. A sound under- standing of how diverse functions are implemented in a simple foot segment during human locomotion might be useful to gain insight into the overall foot locomotor functions and hence to facilitate clinical diagnosis, rehabilitation product design and humanoid robot development.     

Mechanical energy and effective foot mass during impact loading of walking and running

The human heel pad is considered an important structure for attenuation of the transient force caused by heel-strike. Although the mechanical properties of heel pads are relatively well understood, the mechanical energy (Etot) absorbed by the heel pad during the impact phase has never been documented directly because data on the effective foot mass (Meff) was previously unavailable during normal forward locomotion. In this study, we use the impulse-momentum method (IMM) for calculating Meff from moving subjects. Mass-spring-damper models were developed to evaluate errors and to examine the effects of pad property, upper body mass, and effective leg spring on Meff. We simultaneously collected ground reaction forces, pad deformation, and lower limb kinematics during impact phase of barefoot walking, running, and crouched walking. The latter was included to examine the effect of knee angle on Meff. The magnitude of Meff as a percentage of body mass (M(B)) varies with knee angle at impact and significantly differs among gaits: 6.3%M(B) in walking, 5.3%M(B) in running, and 3.7%M(B) in crouched walking. Our modeling results suggested that Meff is insensitive to heel pad resilience and effective leg stiffness. At the instant prior to heel strike, Etot ranges from 0.24 to 3.99 J. The combination of video and forceplate data used in this study allows analyses of Etot and Etot as a function of heel-strike kinematics during normal locomotion. Relationship between Meff and knee angle provides insights into how changes in posture moderate impact transients at different gaits.     

A Dynamic Finite Element Analysis of Human Foot Complex in the Sagittal Plane during Level Walking

The objective of this study is to develop a computational framework for investigating the dynamic behavior and the internal loading conditions of the human foot complex during locomotion. A subject-specific dynamic finite element model in the sagittal plane was constructed based on anatomical structures segmented from medical CT scan images. Three-dimensional gait measurements were conducted to support and validate the model. Ankle joint forces and moment derived from gait measurements were used to drive the model. Explicit finite element simulations were conducted, covering the entire stance phase from heel-strike impact to toe-off. The predicted ground reaction forces, center of pressure, foot bone motions and plantar surface pressure showed reasonably good agreement with the gait measurement data over most of the stance phase. The prediction discrepancies can be explained by the assumptions and limitations of the model. Our analysis showed that a dynamic FE simulation can improve the prediction accuracy in the peak plantar pressures at some parts of the foot complex by 10%–33% compared to a quasi-static FE simulation. However, to simplify the costly explicit FE simulation, the proposed model is confined only to the sagittal plane and has a simplified representation of foot structure. The dynamic finite element foot model proposed in this study would provide a useful tool for future extension to a fully muscle-driven dynamic three-dimensional model with detailed representation of all major anatomical structures, in order to investigate the structural dynamics of the human foot musculoskeletal system during normal or even pathological functioning.     

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.     

Optimal foot shape for a passive dynamic biped

Passive walking dynamics describe the motion of a biped that is able to "walk" down a shallow slope without any actuation or control. Instead, the walker relies on gravitational and inertial effects to propel itself forward, exhibiting a gait quite similar to that of humans. These purely passive models depend on potential energy to overcome the energy lost when the foot impacts the ground. Previous research has demonstrated that energy loss at heel-strike can vary widely for a given speed, depending on the nature of the collision. The point of foot contact with the ground (relative to the hip) can have a significant effect: semi-circular (round) feet soften the impact, resulting in much smaller losses than point-foot walkers. Collisional losses are also lower if a single impulse is broken up into a series of smaller impulses that gradually redirect the velocity of the center of mass rather than a single abrupt impulse. Using this principle, a model was created where foot-strike occurs over two impulses, "heel-strike" and "toe-strike," representative of the initial impact of the heel and the following impact as the ball of the foot strikes the ground. Having two collisions with the flat-foot model did improve efficiency over the point-foot model. Representation of the flat-foot walker as a rimless wheel helped to explain the optimal flat-foot shape, driven by symmetry of the virtual spoke angles. The optimal long period foot shape of the simple passive walking model was not very representative of the human foot shape, although a reasonably anthropometric foot shape was predicted by the short period solution.     

An efficient probabilistic methodology for incorporating uncertainty in body segment parameters and anatomical landmarks in joint loadings estimated from inverse dynamics

Inverse dynamics is a standard approach for estimating joint loadings in the lower extremity from kinematic and ground reaction data for use in clinical and research gait studies. Variability in estimating body segment parameters and uncertainty in defining anatomical landmarks have the potential to impact predicted joint loading. This study demonstrates the application of efficient probabilistic methods to quantify the effect of uncertainty in these parameters and landmarks on joint loading in an inverse-dynamics model, and identifies the relative importance of the parameters and landmarks to the predicted joint loading. The inverse-dynamics analysis used a benchmark data set of lower-extremity kinematics and ground reaction data during the stance phase of gait to predict the three-dimensional intersegmental forces and moments. The probabilistic analysis predicted the 1-99 percentile ranges of intersegmental forces and moments at the hip, knee, and ankle. Variabilities, in forces and moments of up to 56% and 156% of the mean values were predicted based on coefficients of variation less than 0.20 for the body segment parameters and standard deviations of 2 mm for the anatomical landmarks. Sensitivity factors identified the important parameters for the specific joint and component directions. Anatomical landmarks affected moments to a larger extent than body segment parameters. Additionally, for forces, anatomical landmarks had a larger effect than body segment parameters, with the exception of segment masses, which were important to the proximal-distal joint forces. The probabilistic modeling approach predicted the range of possible joint loading, which has implications in gait studies, clinical assessments, and implant design evaluations.     

Determining the centre of pressure during walking and running using an instrumented treadmill

In this paper, a new method of determining spatial and temporal gait parameters by using centre of pressure (CoP) data is presented. A treadmill is used which was developed to overcome limitations of regular methods for the analysis of spatio-temporal gait parameters and ground reaction forces during walking and running. The design of the treadmill is based on the use of force transducers underneath a separate left and right plate, which together form the treadmill walking surface. The results of test procedures and measurements show that accurate recordings of vertical ground reaction force can be obtained. These recordings enable a separate analysis of vertical ground reaction forces during double support phases in walking, and the analysis of changes in the centre of pressure (CoP) position during subsequent foot placements. From the CoP data, temporal gait parameters (e.g. duration of left/right support and swing phases) and spatial gait parameters (i.e. left/right step lengths and widths) can be derived.     

A closed-loop cadaveric foot and ankle loading model

Investigations of human foot and ankle biomechanics rely chiefly on cadaver experiments. The application of proper force magnitudes to the cadaver foot and ankle is essential to obtain valid biomechanical data. Data for external ground reaction forces are readily available from human motion analysis. However, determining appropriate forces for extrinsic foot and ankle muscles is more problematic. A common approach is the estimation of forces from muscle physiological cross-sectional areas and electromyographic data. We have developed a novel approach for loading the Achilles and posterior tibialis tendons that does not prescribe predetermined muscle forces. For our loading model, these muscle forces are determined experimentally using independent plantarflexion and inversion angle feedback control. The independent (input) parameters -- calcaneus plantarflexion, calcaneus inversion, ground reaction forces, and peroneus forces -- are specified. The dependent (output) parameters -- Achilles force, posterior tibialis force, joint motion, and spring ligament strain -- are functions of the independent parameters and the kinematics of the foot and ankle. We have investigated the performance of our model for a single, clinically relevant event during the gait cycle. The instantaneous external forces and foot orientation determined from human subjects in a motion analysis laboratory were simulated in vitro using closed-loop feedback control. Compared to muscle force estimates based on physiological cross-sectional area data and EMG activity at 40% of the gait cycle, the posterior tibialis force and Achilles force required when using position feedback control were greater.     

The influence of foot position on body dynamics

During locomotion, the human body exhibits inherent dynamic properties such as mass (M), stiffness (K) and damping (B). During the gait cycle, foot contact with the ground progresses from the heel to the toe. Contact forces between the foot and ground are defined as ground reaction forces (GRF). It is unclear how body dynamics are affected by foot landing position. If the shape of GRF is indicative of body dynamics, our understanding of gait patterns in normal and pathologic conditions may improve. The aims of this study were to determine:(1) whether foot landing position affects the inherent dynamics of the human body and (2) the extent to which the GRF curve reflects the response of inherent body dynamics to sudden loading.Eight non-disabled control volunteers performed a series of small jumps and landed on one leg with a fully extended knee in three foot landing positions: heel, mid-foot, and toe. They then walked at self-paced velocity over force plates. For each foot landing position, values of K, B and the dimensionless damping coefficient, ξ, were calculated from the period of vertical body oscillations, T, and compared with an ANOVA test. In addition, the time between the two peaks of the vertical GRF, T_GRF, was compared with T. We found that that K and B decreased and ξ did not change (p<0.01) between heel to toe-landing positions. T_GRF was not different than T for the toe-landing position, which suggests that the dynamic body response has major impact on the shape of GRF.     

Pedestrian locomotion energetics and gait characteristics of a diving bird, the great cormorant, Phalacrocorax carbo

Great cormorants Phalacrocorax carbo are foot propelled diving birds that seem poorly suited to locomotion on land. They have relatively short legs, which are presumably adapted for the generation of high forces during the power stroke of aquatic locomotion, and walk with a pronounced "clumsy waddle". We hypothesise (1) that the speed, independent minimum cost of locomotion (C min, ml O2 m(-1)) will be high for cormorants during treadmill exercise, and (2) that cormorants will have a relatively limited speed range in comparison to more cursorial birds. We measured the rate of oxygen consumption (V02) of cormorants during pedestrian locomotion on a treadmill, and filmed them to determine duty factor (the fraction of stride period that the foot is in contact with the ground), foot contact time (tc), stride frequency (f), swing phase duration and stride length. C min was 2.1-fold higher than that predicted by their body mass and phylogenetic position, but was not significantly different from the C min of runners (Galliformes and Struthioniformes). The extrapolated gamma-intercept of the relationship between V02 and speed was 1.9-fold higher than that predicted by allometry. Again, cormorants were not significantly different from runners. Contrary to our hypothesis, we therefore conclude that cormorants do not have high pedestrian transport costs. Cormorants were observed to use a grounded gait with two double support phases at all speeds measured, and showed an apparent gait transition between 0.17 and 0.25 m s(-1). This transition occurs at a Froude number between 0.016 and 0.037, which is lower than the value of approximately 0.5 observed for many other species. However, despite the use of a limited speed range, and a gait transition at relatively low speed, we conclude that the pedestrian locomotion of these foot propelled diving birds is otherwise generally similar to that of cursorial birds at comparable relative velocities.     

Dynamic plantar pressure distribution during terrestrial locomotion of bonobos (Pan paniscus)

We collected high-resolution plantar pressure distributions of seven bonobos during terrestrial bipedal and quadrupedal locomotion (N = 146). Functional foot length, degree of hallux abduction, and total contact time were determined, and plots, showing pressure as a function of time for six different foot regions, were generated. We also studied five adult humans for comparison (N = 13). Both locomotion types of the bonobo show a large variation in plantar pressure distributions, which could be due to the interference of instantaneous behavior with locomotion and differences in walking speed and body dimensions. The heel and the lateral midfoot typically touch down simultaneously at initial ground contact in bipedal and quadrupedal walking of bonobos, in contrast with the typical heel-strike of human bipedalism. The center of pressure follows a curved course during quadrupedalism, as a consequence of the medial weight transfer during mid-stance. Bipedal locomotion of bonobos is characterized by a more plantar positioning of the feet and by a shorter contact time than during quadrupedal walking, according to a smaller stride and step length at a higher frequency. We observed a varus position of the foot with an abducted hallux, which likely possesses an important sustaining and stabilizing function during terrestrial locomotion.     

Pattern of anterior cruciate ligament force in normal walking

The goal of this study was to calculate and explain the pattern of anterior cruciate ligament (ACL) loading during normal level walking. Knee-ligament forces were obtained by a two-step procedure. First, a three-dimensional (3D) model of the whole body was used together with dynamic optimization theory to calculate body-segmental motions, ground reaction forces, and leg-muscle forces for one cycle of gait. Joint angles, ground reaction forces, and muscle forces obtained from the gait simulation were then input into a musculoskeletal model of the lower limb that incorporated a 3D model of the knee. The relative positions of the femur, tibia, and patella and the forces induced in the knee ligaments were found by solving a static equilibrium problem at each instant during the simulated gait cycle. The model simulation predicted that the ACL bears load throughout stance. Peak force in the ACL (303 N) occurred at the beginning of single-leg stance (i.e., contralateral toe off). The pattern of ACL force was explained by the shear forces acting at the knee. The balance of muscle forces, ground reaction forces, and joint contact forces applied to the leg determined the magnitude and direction of the total shear force acting at the knee. The ACL was loaded whenever the total shear force pointed anteriorly. In early stance, the anterior shear force from the patellar tendon dominated the total shear force applied to the leg, and so maximum force was transmitted to the ACL at this time. ACL force was small in late stance because the anterior shear forces supplied by the patellar tendon, gastrocnemius, and tibiofemoral contact were nearly balanced by the posterior component of the ground reaction.     

Ultrasonic motion analysis system--measurement of temporal and spatial gait parameters

The duration of stance and swing phase and step and stride length are important parameters in human gait. In this technical note a low-cost ultrasonic motion analysis system is described that is capable of measuring these temporal and spatial parameters while subjects walk on the floor. By using the propagation delay of sound when transmitted in air, this system is able to record the position of the subjects' feet. A small ultrasonic receiver is attached to both shoes of the subject while a transmitter is placed stationary on the floor. Four healthy subjects were used to test the device. Subtracting positions of the foot with zero velocity yielded step and stride length. The duration of stance and swing phase was calculated from heel-strike and toe-off. Comparison with data obtained from foot contact switches showed that applying two relative thresholds to the speed graph of the foot could reliably generate heel-strike and toe-off. Although the device is tested on healthy subjects in this study, it promises to be extremely valuable in examining pathological gait. When gait is asymmetrical, walking speed is not constant or when patients do not completely lift their feet, most existing devices will fail to correctly assess the proper gait parameters. Our device does not have this shortcoming and it will accurately demonstrate asymmetries and variations in the patient's gait. As an example, the recording of a left hemiplegic patient is presented in the discussion.     

Change in knee contact force with simulated change in body weight

The relationship between obesity, weight gain and progression of knee osteoarthritis is well supported, suggesting that excessive joint loading may be a mechanism responsible for cartilage deterioration. Examining the influence of weight gain on joint compressive forces is difficult, as both muscles and ground reaction forces can have a significant impact on the forces experienced during gait. While previous studies have examined the relationship between body weight and knee forces, these studies have used models that were not validated using experimental data. Therefore, the objective of this study was to evaluate the relationship between changes in body weight and changes in knee joint contact forces for an individual's gait pattern using musculoskeletal modeling that is validated against known internal compressive forces. Optimal weighting constants were determined for three subjects to generate valid predictions of knee contact forces (KCFs) using in vivo data collection with instrumented total knee arthroplasty. A total of five simulations per walking trial were generated for each subject, from 80% to 120% body weight in 10% increments, resulting in 50 total simulations. The change in peak KCF with respect to body weight was found to be constant and subject-specific, predominantly determined by the peak force during the baseline condition at 100% body weight. This relationship may be further altered by any change in kinematics or body mass distribution that may occur as a result of a change in body weight or exercise program.     

Measured and estimated ground reaction forces for multi-segment foot models

Accurate measurement of ground reaction forces under discrete areas of the foot is important in the development of more advanced foot models, which can improve our understanding of foot and ankle function. To overcome current equipment limitations, a few investigators have proposed combining a pressure mat with a single force platform and using a proportionality assumption to estimate subarea shear forces and free moments. In this study, two adjacent force platforms were used to evaluate the accuracy of the proportionality assumption on a three segment foot model during normal gait. Seventeen right feet were tested using a targeted walking approach, isolating two separate joints: transverse tarsal and metatarsophalangeal. Root mean square (RMS) errors in shear forces up to 6% body weight (BW) were found using the proportionality assumption, with the highest errors (peak absolute errors up to 12% BW) occurring between the forefoot and toes in terminal stance. The hallux exerted a small braking force in opposition to the propulsive force of the forefoot, which was unaccounted for by the proportionality assumption. While the assumption may be suitable for specific applications (e.g. gait analysis models), it is important to understand that some information on foot function can be lost. The results help highlight possible limitations of the assumption. Measured ensemble average subarea shear forces during normal gait are also presented for the first time.     

An MRI-compatible foot-loading device for assessment of internal tissue deformation

It is well known that mechanical forces acting within the soft tissues of the foot can contribute to the formation of neuropathic ulcers in people with diabetes. Presently, only surface measurements of plantar pressure are used clinically to estimate risk status due to mechanical loading. It is currently not known how surface measurements relate to the three-dimensional (3-D) internal stress/strain state of the foot. This article describes the development of a foot-loading device that allows for the direct observation of the internal deformation of foot tissues under known forces. Ground reaction forces and plantar pressure distributions during normal walking were measured in ten healthy young adults. One instant in the gait cycle, when pressure under the metatarsal heads reached a peak, was extracted for simulation in an MR imager. T1-weighted 3-D gradient echo MRI sets were collected as the simulated walking ground reaction force was incrementally applied to the foot by the novel foot-loading device. The sub-metatarsal head soft-tissue thickness decreased rapidly at first and then reached a plateau. Peak plantar pressure measurements collected within the loading device (161+/-75kPa) were lower in magnitude and less focal than pressures measured during walking (492+/-91kPa). This finding implies that although the device successfully applied full peak walking ground reaction forces to the foot, they were not distributed in the same manner as during walking. Although not representative of gait, the data collected from this in vivo mechanical test are suitable for determination of foot tissue material properties or, when combined with finite element modeling, to examine the relationship between surface loading and internal stress.     

Peak dynamic force in human gait

Studies were made of the forces generated at heel stroke in human gait using both force plates having a high resonant frequencies (capable of picking up high frequency components in the contact force) as well as a force transducer inserted into the heel of the shoe of the subjects. The output traces were analyzed for the existence of high frequency impulsive loads during a normal walking cycle. The effect of the complicance of the foot and floor was studied with the force transducers. The results showed that during normal human gait the lower limb is subjected to a high frequency impulsive load at heel strike. The severity of this impulse varied with the individual, the velocity and angle with which the limb aproached the ground and the compliance of the two materials coming in contact at heel strike. The magnitude of this peak force varied from 0.5 to 1.25 times body weight and its frequency components from 10 to 75 Hz.     

Biomechanical analysis of the three-dimensional foot structure during gait: a basic tool for clinical applications

A novel three-dimensional numerical model of the foot, incorporating, for the first time in the literature, realistic geometric and material properties of both skeletal and soft tissue components of the foot, was developed for biomechanical analysis of its structural behavior during gait. A system of experimental methods, integrating the optical Contact Pressure Display (CPD) method for plantar pressure measurements and a Digital Radiographic Fluoroscopy (DRF) instrument for acquisition of skeletal motion during gait, was also developed in this study and subsequently used to build the foot model and validate its predictions. Using a Finite Element solver, the stress distribution within the foot structure was obtained and regions of elevated stresses for six subphases of the stance (initial-contact, heel-strike, midstance, forefoot-contact, push-off, and toe-off) were located. For each of these subphases, the model was adapted according to the corresponding fluoroscopic data, skeletal dynamics, and active muscle force loading. Validation of the stress state was achieved by comparing model predictions of contact stress distribution with respective CPD measurements. The presently developed measurement and numerical analysis tools open new approaches for clinical applications, from simulation of the development mechanisms of common foot disorders to pre- and post-interventional evaluation of their treatment.     

Assessment of the influence of foot orthoses in the hip loading conditions during walking: a single case study

Despite their large clinical application, the understanding of the effects of foot orthoses on the lower limb kinematics and kinetics is limited. In this context, we propose an advanced musculoskeletal model to assess the influence of foot orthoses in the loading conditions within an osteoarthritic hip joint during gait. Experimental data are collected for a single pathological subject presenting a coxarthrosis (with and without orthoses), and a healthy subject during walking. An inverse dynamic approach coupled with an optimisation method evaluates the forces developed by 14 muscles and the hip contact reaction force. Contact reaction and muscular force magnitudes are closed whether the patient is walking with or without foot orthoses. Nevertheless, contact reaction amplitudes and orientations show differences in relation to those calculated for the healthy subject. The results obtained allow us to formulate some assumptions concerning the causes of coxarthrosis evolution and treatment.