The foot as a shock absorber

A mathematical analysis of the deformation of the foot is developed to determine the role that stretch of ligaments and tendons plays in absorbing shock following impact. Our analysis is based on an anatomical biomechanical model that includes each of the bones of the foot. We calculate the time course of the deflection of the joints and the elongation of the ligaments and tendons and determine the ground reaction force acting on the heel. Quasi-linear viscoelastic theory is used for soft tissue constitutive relationships. With biomechanical data selected from the literature, we obtain a vertical force impact peak of 8000 N, occurring at 16 ms following heel strike. This is of higher magnitude and shorter duration than is found experimentally, as is to be expected, since we did not include the heel pad in our model and we assumed that the impact surface was ideally rigid.     

The potential influence of the heel counter on internal stress during static standing: a combined finite element and positional MRI investigation

Confinement of the heel due to the counter of the shoe is believed to influence heel pad biomechanics. Using a two-dimensional finite element model of the heel pad and shoe during a simulation of static standing, the aim of this study was to quantify the potential effect of confinement on internal heel pad stress. Non-weightbearing MRI and weightbearing MRI with plantar pressure and ground reaction force data were recorded for a single subject. The non-weightbearing MRI was used to create two FE models of the heel pad, using either homogeneous or composite material properties. The composite model included a distinction in material properties between fat pad and skin. Vertical and medial-lateral forces, as measured on the subject's heel, were applied to the models and vertical compressive strains for both models were comparable with those observed by weightbearing MRI. However, only for the composite model was the predicted plantar pressure distribution comparable with measured data. The composite model was therefore used in further analyses. In this composite model, the internal stresses were located mainly in the skin and were predominantly tensile in nature, whereas the stress state in the fat pad approached hydrostatic conditions. A representation of a running shoe, including an insole, midsole and heel counter was then added to the composite heel pad to form the shod model. In order to investigate the counter effect, the load was applied to the shod model with and without the heel counter. The effect of the counter on peak stress was to elevate compression (0-50%), reduce tension (22-34%) and reduce shear (22-28%) in the skin. In addition, the counter reduced both compressive (20-40%) and shear (58-80%) stress in the fat pad and tension in the fat pad remained negligible. Taken together the results indicate that a well-fitted counter works in sympathy with the internal structure of the heel pad and could be an effective reducer of heel pad stress. However, further research needs to be undertaken to assess the long-term effects on the soft-tissues, practicalities of achieving good fit and behavior under dynamic events.     

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.     

Greater toe grip and gentler heel strike are the strategies to adapt to slippery surface

This study investigated the plantar pressure distribution during gait on wooden surface with different slipperiness in the presence of contaminants. Fifteen Chinese males performed 10 walking trials on a 5-m wooden walkway wearing cloth shoe in four contaminated conditions (dry, sand, water, oil). A pressure insole system was employed to record the plantar pressure data at 50Hz. Peak pressure and time-normalized pressure-time integral were evaluated in nine regions. In comparing walking on slippery to non-slippery surfaces, results showed a 30% increase of peak pressure beneath the hallux (from 195.6 to 254.1kPa), with a dramatic 79% increase in the pressure time integral beneath the hallux (from 63.8 to 114.3kPa) and a 34% increase beneath the lateral toes (from 35.1 to 47.2kPa). In addition, the peak pressure beneath the medial and lateral heel showed significant 20-24% reductions, respectively (from 233.6-253.5 to 204.0-219.0kPa). These findings suggested that greater toe grip and gentler heel strike are the strategies to adapt to slippery surface. Such strategies plantarflexed the ankle and the metatarsals to achieve a flat foot contact with the ground, especially at heel strike, in order to shift the ground reaction force to a more vertical direction. As the vertical ground reaction force component increased, the available ground friction increased and the floor became less slippery. Therefore, human could walk without slip on slippery surfaces with greater toe grip and gentler heel strike as adaptation strategies.     

Estimating the complete ground reaction forces with pressure insoles in walking

This study presented a method to estimate the complete ground reaction forces from pressure insoles in walking. Five male subjects performed 10 walking trials in a laboratory. The complete ground reaction forces were collected during a right foot stride by a force plate at 1000Hz. Simultaneous plantar pressure data were collected at 100Hz by a pressure insole system with 99 sensors covering the whole plantar area. Stepwise linear regressions were performed to individually reconstruct the complete ground reaction forces in three directions from the 99 individual pressure data until redundancy among the predictors occurred. An additional linear regression was performed to reconstruct the vertical ground reaction force by the sum of the value of the 99 pressure sensors. Five other subjects performed the same walking test for validation. Estimated ground reaction forces in three directions were calculated with the developed regression models, and were compared with the real data recorded from force plate. Accuracy was represented by the correlation coefficient and the root mean square error. Results showed very good correlation in anterior-posterior (0.928) and vertical (0.989) directions, and reasonable correlation in medial-lateral direction (0.719). The root mean square error was about 12%, 5% and 28% of the peak recorded value. Future studies should aim to generalize the methods or to establish specific methods to other subjects, patients, motions, footwear and floor conditions. The method gives an extra option to study an estimation of the complete ground reaction forces in any environment without the constraints from the number and location of force plates.     

Force measurements during running on different instrumented treadmills

One method to determine the forces produced during running is to conduct extensive kinematic and kinetic analysis. These analyses can be performed by having an individual perform repeated over-ground running trials or simply run continuously on an instrumented treadmill. The forces produced during over-ground running may not be the same as the forces during treadmill running and these differences could be attributed to a number of factors, including the design of the instrumented treadmill. The purpose of this paper was to determine whether there are differences in force measurements on different instrumented treadmill setups in comparison to over-ground running and to correct for any of these differences using a theoretical model. 11 participants ran on three different treadmills and performed over-ground running at 2.7, 3.6, and 4.5 m/s. Ground reaction forces were measured via force plates and an instrumented pressure insole. We found that the magnitude of the vertical ground reaction force differed between the three treadmills and over-ground running. The difference in ground reaction forces estimated by the pressure insole and the treadmill-force-plate system or instrumented treadmill can be explained by a three degree of freedom mechanical model of a person running on a treadmill and this model could potentially be used to correct for errors in force measurement from instrumented treadmills. The model included a force plate, a treadmill, and a wobbling mass with varying natural frequencies and damping characteristics, and constant masses. These findings provide researchers a method to correct forces from an instrumented treadmill set-up to determine a close approximation of the actual forces experienced by a participant during treadmill running.     

Correlation between ground reaction force and tibial acceleration in vertical jumping

Modern electronics allow for the unobtrusive measurement of accelerations outside the laboratory using wireless sensor nodes. The ability to accurately measure joint accelerations under unrestricted conditions, and to correlate them with jump height and landing force, could provide important data to better understand joint mechanics subject to real-life conditions. This study investigates the correlation between peak vertical ground reaction forces, as measured by a force plate, and tibial axial accelerations during free vertical jumping. The jump heights calculated from force-plate data and accelerometer measurements are also compared. For six male subjects participating in this study, the average coefficient of determination between peak ground reaction force and peak tibial axial acceleration is found to be 0.81. The coefficient of determination between jump height calculated using force plate and accelerometer data is 0.88. Data show that the landing forces could be as high as 8 body weights of the jumper. The measured peak tibial accelerations ranged up to 42 g. Jump heights calculated from force plate and accelerometer sensors data differed by less than 2.5 cm. It is found that both impact accelerations and landing forces are only weakly correlated with jump height (the average coefficient of determination is 0.12). This study shows that unobtrusive accelerometers can be used to determine the ground reaction forces experienced in a jump landing. Whereas the device also permitted an accurate determination of jump height, there was no correlation between peak ground reaction force and jump height.     

Simple and complex models for studying muscle function in walking

While simple models can be helpful in identifying basic features of muscle function, more complex models are needed to discern the functional roles of specific muscles in movement. In this paper, two very different models of walking, one simple and one complex, are used to study how muscle forces, gravitational forces and centrifugal forces (i.e. forces arising from motion of the joints) combine to produce the pattern of force exerted on the ground. Both the simple model and the complex one predict that muscles contribute significantly to the ground force pattern generated in walking; indeed, both models show that muscle action is responsible for the appearance of the two peaks in the vertical force. The simple model, an inverted double pendulum, suggests further that the first and second peaks are due to net extensor muscle moments exerted about the knee and ankle, respectively. Analyses based on a much more complex, muscle-actuated simulation of walking are in general agreement with these results; however, the more detailed model also reveals that both the hip extensor and hip abductor muscles contribute significantly to vertical motion of the centre of mass, and therefore to the appearance of the first peak in the vertical ground force, in early single-leg stance. This discrepancy in the model predictions is most probably explained by the difference in model complexity. First, movements of the upper body in the sagittal plane are not represented properly in the double-pendulum model, which may explain the anomalous result obtained for the contribution of a hip-extensor torque to the vertical ground force. Second, the double-pendulum model incorporates only three of the six major elements of walking, whereas the complex model is fully 3D and incorporates all six gait determinants. In particular, pelvic list occurs primarily in the frontal plane, so there is the potential for this mechanism to contribute significantly to the vertical ground force, especially during early single-leg stance when the hip abductors are activated with considerable force.     

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.     

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.     

Towards a footwear design tool: Influence of shoe midsole properties and ground stiffness on the impact force during running

Several spring–damper–mass models of the human body have been developed in order to reproduce the measured ground vertical reaction forces during human running (McMahon and Cheng, 1990; Ferris et al., 1999; Liu and Nigg, 2000). In particular, Liu and Nigg introduced at the lower level of their model, i.e. at the interface between the human body and the ground, a nonlinear element representing simultaneously the shoe midsoles and the ground flexibility. The ground reaction force is modelled as the force supported by this nonlinear element, whose parameters are identified from several sets of experimental data. This approach proved to be robust and quite accurate. However, it does not explicitly take into account the shoe and the ground properties. It turns out to be impossible to study the influence of shoe materials on the impact force, for instance for footwear design purposes. In this paper, a modification of the Liu and Nigg's model is suggested, where the original nonlinear element is replaced with a bi-layered spring–damper–mass model: the first layer represents the shoe midsole and the second layer is associated with the ground.Ground is modelled as an infinite elastic half-space. We have assumed a viscoelastic behaviour of the shoe material, so the damping of shoe material is taken into account. A methodology for the shoe-soles characterization is proposed and used together with the proposed model. A parametric study is then conducted and the influence of the shoe properties on the impact force is quantified. Moreover, it is shown that impact forces are strongly affected by the ground stiffness, which should therefore be considered as an essential parameter in the footwear design.     

Modelling, simulation and optimisation of a human vertical jump.

This paper describes an efficient biomechanical model of the human lower limb with the aim of simulating a real human jump movement consisting of an upword propulsion, a flying and a landing phase. A multiphase optimal control technique is used to solve the muscle force sharing problem. To understand how intermuscular control coordinates limb muscle excitations, the human body is reduced to a single lower limb consisting of three rigid bodies. The biomechanical system is activated by nine muscle-tendon actuators representing the basic properties of muscles during force generation. For the calculation of the minimal muscle excitations of the jump movement, the trajectory of the hip joint is given as a rheonomic constraint and the contact forces (ground reaction forces) are determined by force plates. Based on the designed musculoskeletal model and on the differential equations of the multibody system, muscle excitations and muscle forces necessary for a vertical jump movement are calculated. The validity of the system is assessed comparing the calculated muscle excitations with the registered surface electromyogramm (EMG) of the muscles. The achieved results indicate a close relationship between the predicted and the measured parameters.     

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.     

Use of pressure insoles to calculate the complete ground reaction forces

A method to calculate the complete ground reaction force (GRF) components from the vertical GRF measured with pressure insoles is presented and validated. With this approach it is possible to measure several consecutive steps without any constraint on foot placement and compute a standard inverse dynamics analysis with the estimated GRF.     

The ground reaction force pattern from the hindlimb of the horse simulated by a spring model

A model consisting of a spring loaded by a time-dependent mass is presented simulating the vertical and longitudinal horizontal ground reaction force patterns obtained from the hindlimb of a walking horse.     

The oscillatory behavior of the CoM facilitates mechanical energy balance between push-off and heel strike

Humans use equal push-off and heel strike work during the double support phase to minimize the mechanical work done on the center of mass (CoM) during the gait. Recently, a step-to-step transition was reported to occur over a period of time greater than that of the double support phase, which brings into question whether the energetic optimality is sensitive to the definition of the step-to-step transition. To answer this question, the ground reaction forces (GRFs) of seven normal human subjects walking at four different speeds (1.1-2.4 m/s) were measured, and the push-off and heel strike work for three differently defined step-to-step transitions were computed based on the force, work, and velocity. To examine the optimality of the work and the impulse data, a hybrid theoretical-empirical analysis is presented using a dynamic walking model that allows finite time for step-to-step transitions and incorporates the effects of gravity within this period. The changes in the work and impulse were examined parametrically across a range of speeds. The results showed that the push-off work on the CoM was well balanced by the heel strike work for all three definitions of the step-to-step transition. The impulse data were well matched by the optimal impulse predictions (R(2)>0.7) that minimized the mechanical work done on the CoM during the gait. The results suggest that the balance of push-off and heel strike energy is a consistent property arising from the overall gait dynamics, which implies an inherited oscillatory behavior of the CoM, possibly by spring-like leg mechanics.     

The equations of motion for a standing human reveal three mechanisms for balance

The equations of motion for a standing multi-segment human model are derived. Output quantity of these equations is the horizontal acceleration of the whole-body centre of mass (CoM). There are three input terms and they can be identified as the three mechanisms by which balance can be maintained: (1) by moving the centre of pressure with respect to the vertical projection of the CoM, (2) by counter-rotating segments around the CoM, and (3) by applying an external force, other than the ground reaction force. For the first two mechanisms the respective contributions to CoM acceleration can be obtained from force plate recordings. This is illustrated by some example data from experiments, which show that the contribution from mechanism 2 can be considerable, e.g. in one-legged standing.     

The influence of soft tissue movement on ground reaction forces, joint torques and joint reaction forces in drop landings

The aim of this study was to determine the effects that soft tissue motion has on ground reaction forces, joint torques and joint reaction forces in drop landings. To this end a four body-segment wobbling mass model was developed to reproduce the vertical ground reaction force curve for the first 100 ms of landing. Particular attention was paid to the passive impact phase, while selecting most model parameters a priori, thus permitting examination of the rigid body assumption on system kinetics. A two-dimensional wobbling mass model was developed in DADS (version 9.00, CADSI) to simulate landing from a drop of 43 cm. Subject-specific inertia parameters were calculated for both the rigid links and the wobbling masses. The magnitude and frequency response of the soft tissue of the subject to impulsive loading was measured and used as a criterion for assessing the wobbling mass motion. The model successfully reproduced the vertical ground reaction force for the first 100 ms of the landing with a peak vertical ground reaction force error of 1.2% and root mean square errors of 5% for the first 15 ms and 12% for the first 40 ms. The resultant joint forces and torques were lower for the wobbling mass model compared with a rigid body model, up to nearly 50% lower, indicating the important contribution of the wobbling masses on reducing system loading.     

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.