Ability of the planar spring-mass model to predict mechanical parameters in running humans

The planar spring-mass model is a simple mathematical model of bouncing gaits, such as running, trotting and hopping. Although this model has been widely used in the study of locomotion, its accuracy in predicting locomotor mechanics has not been systematically quantified. We determined the percent error of the model in predicting 10 locomotor parameters in running humans by comparing the model predictions to experimental data from humans running in normal gravity and simulated reduced gravity. We tested the hypotheses that the model would overestimate horizontal impulse and the change in mechanical energy of the centre of mass (COM) during stance. The model provided good predictions of stance time, vertical impulse, contact length, duty factor, relative stride length and relative peak force. All predictions of these parameters were within 20% of measured values and at least 90% of predictions of each parameter were within 10% of measured values (median absolute errors: <7%). This suggests that the model incorporates all features of running humans that have a significant influence upon these six parameters. As simulated gravity level decreased, the magnitude of the errors in predicting each of these parameters either decreased or stayed constant, indicating that this is a good model of running in simulated reduced gravity. As hypothesised, horizontal impulse and change in mechanical energy of the COM during stance were overestimated (median absolute errors: 43.6% and 26.2%, respectively). Aerial time and peak vertical COM displacement during stance were also systematically overestimated (median absolute errors: 17.7% and 22.9%, respectively). Care should be taken to ensure that the model is used only to investigate parameters which it can predict accurately. It would be useful to extend this analysis to other species and gaits.     

A simple method for measuring stiffness during running

The spring-mass model, representing a runner as a point mass supported by a single linear leg spring, has been a widely used concept in studies on running and bouncing mechanics. However, the measurement of leg and vertical stiffness has previously required force platforms and high-speed kinematic measurement systems that are costly and difficult to handle in field conditions. We propose a new "sine-wave" method for measuring stiffness during running. Based on the modeling of the force-time curve by a sine function,this method allows leg and vertical stiffness to be estimated from just a few simple mechanical parameters: body mass, forward velocity, leg length, flight time, and contact time. We compared this method to force-platform-derived stiffness measurements for treadmill dynamometer and overground running conditions, at velocities ranging from 3.33 m.s-1 to maximal running velocity in both recreational and highly trained runners. Stiffness values calculated with the proposed method ranged from 0.67 % to 6.93 % less than the force platform method, and thus were judged to be acceptable. Furthermore, significant linear regressions (p < 0.01) close to the identity line were obtained between force platform and sine-wave model values of stiffness. Given the limits inherent in the use of the spring-mass model, it was concluded that this sine-wave method allows leg and stiffness estimates in running on the basis of a few mechanical parameters, and could be useful in further field measurements.     

Gender differences in active musculoskeletal stiffness. Part II. Quantification of leg stiffness during functional hopping tasks.

Leg stiffness was compared between age-matched males and females during hopping at preferred and controlled frequencies. Stiffness was defined as the linear regression slope between the vertical center of mass (COM) displacement and ground-reaction forces recorded from a force plate during the stance phase of the hopping task. Results demonstrate that subjects modulated the vertical displacement of the COM during ground contact in relation to the square of hopping frequency. This supports the accuracy of the spring-mass oscillator as a representative model of hopping. It also maintained peak vertical ground-reaction load at approximately three times body weight. Leg stiffness values in males (33.9+/-8.7 kN/m) were significantly (p<0.01) greater than in females (26.3+/-6.5 kN/m) at each of three hopping frequencies, 3.0, 2.5 Hz, and a preferred hopping rate. In the spring-mass oscillator model leg stiffness and body mass are related to the frequency of motion. Thus male subjects necessarily recruited greater leg stiffness to drive their heavier body mass at the same frequency as the lighter female subjects during the controlled frequency trials. However, in the preferred hopping condition the stiffness was not constrained by the task because frequency was self-selected. Nonetheless, both male and female subjects hopped at statistically similar preferred frequencies (2.34+/-0.22 Hz), therefore, the females continued to demonstrate less leg stiffness. Recognizing the active muscle stiffness contributes to biomechanical stability as well as leg stiffness, these results may provide insight into the gender bias in risk of musculoskeletal knee injury.     

Stance leg control: variation of leg parameters supports stable hopping

The spring-loaded inverted pendulum describes the planar center-of-mass dynamics of legged locomotion. This model features linear springs with constant parameters as legs. In biological systems, however, spring-like properties of limbs can change over time. Therefore, in this study, it is asked how variation of spring parameters during ground contact would affect the dynamics of the spring-mass model. Neglecting damping initially, it is found that decreasing stiffness and increasing rest length of the leg during a stance phase are required for orbitally stable hopping. With damping, stable hopping is found for a larger region of rest-length rates and stiffness rates. Here, also increasing stiffness and decreasing rest length can result in stable hopping. Within the predicted range of leg parameter variations for stable hopping, there is no need for precise parameter tuning. Since hopping gaits form a subset of the running gaits (with vanishing horizontal velocity), these results may help to improve leg design in robots and prostheses.     

Compliant bipedal model with the center of pressure excursion associated with oscillatory behavior of the center of mass reproduces the human gait dynamics

Although the compliant bipedal model could reproduce qualitative ground reaction force (GRF) of human walking, the model with a fixed pivot showed overestimations in stance leg rotation and the ratio of horizontal to vertical GRF. The human walking data showed a continuous forward progression of the center of pressure (CoP) during the stance phase and the suspension of the CoP near the forefoot before the onset of step transition. To better describe human gait dynamics with a minimal expense of model complexity, we proposed a compliant bipedal model with the accelerated pivot which associated the CoP excursion with the oscillatory behavior of the center of mass (CoM) with the existing simulation parameter and leg stiffness. Owing to the pivot acceleration defined to emulate human CoP profile, the arrival of the CoP at the limit of the stance foot over the single stance duration initiated the step-to-step transition. The proposed model showed an improved match of walking data. As the forward motion of CoM during single stance was partly accounted by forward pivot translation, the previously overestimated rotation of the stance leg was reduced and the corresponding horizontal GRF became closer to human data. The walking solutions of the model ranged over higher speed ranges (~1.7 m/s) than those of the fixed pivoted compliant bipedal model (~1.5 m/s) and exhibited other gait parameters, such as touchdown angle, step length and step frequency, comparable to the experimental observations. The good matches between the model and experimental GRF data imply that the continuous pivot acceleration associated with CoM oscillatory behavior could serve as a useful framework of bipedal model.     

The effect of speed on leg stiffness and joint kinetics in human running

The goals of this study were to examine the following hypotheses: (a) there is a difference between the theoretically calculated (McMahon and Cheng, 1990. Journal of Biomechanics 23, 65-78) and the kinematically measured length changes of the spring-mass model and (b) the leg spring stiffness, the ankle spring stiffness and the knee spring stiffness are influenced by running speed. Thirteen athletes took part in this study. Force was measured using a "Kistler" force plate (1000 Hz). Kinematic data were recorded using two high-speed (120 Hz) video cameras. Each athlete completed trials running at five different velocities (approx. 2.5, 3.5, 4.5, 5.5 and 6.5 m/s). Running velocity influences the leg spring stiffness, the effective vertical spring stiffness and the spring stiffness at the knee joint. The spring stiffness at the ankle joint showed no statistical difference (p < 0.05) for the five velocities. The theoretically calculated length change of the spring-mass model significantly (p < 0.05) overestimated the actual length change. For running velocities up to 6.5 m/s the leg spring stiffness is influenced mostly by changes in stiffness at the knee joint.     

Running in the real world: adjusting leg stiffness for different surfaces.

A running animal coordinates the actions of many muscles, tendons, and ligaments in its leg so that the overall leg behaves like a single mechanical spring during ground contact. Experimental observations have revealed that an animal''s leg stiffness is independent of both speed and gravity level, suggesting that it is dictated by inherent musculoskeletal properties. However, if leg stiffness was invariant, the biomechanics of running (e.g. peak ground reaction force and ground contact time) would change when an animal encountered different surfaces in the natural world. We found that human runners adjust their leg stiffness to accommodate changes in surface stiffness, allowing them to maintain similar running mechanics on different surfaces. These results provide important insight into mechanics and control of animal locomotion and suggest that incorporating an adjustable leg stiffness in the design of hopping and running robots is important if they are to match the agility and speed of animals on varied terrain.     

The spring-mass model for running and hopping

A simple spring-mass model consisting of a massless spring attached to a point mass describes the interdependency of mechanical parameters characterizing running and hopping of humans as a function of speed. The bouncing mechanism itself results in a confinement of the free parameter space where solutions can be found. In particular, bouncing frequency and vertical displacement are closely related. Only a few parameters, such as the vector of the specific landing velocity and the specific leg length, are sufficient to determine the point of operation of the system. There are more physiological constraints than independent parameters. As constraints limit the parameter space where hopping is possible, they must be tuned to each other in order to allow for hopping at all. Within the range of physiologically possible hopping frequencies, a human hopper selects a frequency where the largest amount of energy can be delivered and still be stored elastically. During running and hopping animals use flat angles of the landing velocity resulting in maximum contact length. In this situation ground reaction force is proportional to specific contact time and total displacement is proportional to the square of the step duration. Contact time and hopping frequency are not simply determined by the natural frequency of the spring-mass system, but are influenced largely by the vector of the landing velocity. Differences in the aerial phase or in the angle of the landing velocity result in the different kinematic and dynamic patterns observed during running and hopping. Despite these differences, the model predicts the mass specific energy fluctuations of the center of mass per distance to be similar for runners and hoppers and similar to empirical data obtained for animals of various size.     

An actuated dissipative spring-mass walking model: Predicting human-like ground reaction forces and the effects of model parameters

Simple models are widely used to understand the mechanics of human walking. The optimization-based minimal biped model and spring-loaded-inverted-pendulum (SLIP) model are two popular models that can achieve human-like walking patterns. However, ground reaction forces (GRF) from these two models still deviate from experimental data. In this paper, we proposed an actuated dissipative spring-mass model by integrating these two models to realize more human-like GRF patterns. We first explored the function of stiffness, damping, and weights of both energy cost and force cost in the objective function and found that these parameters have distinctly different influences on the optimized gait and GRF profiles. The stiffness and objective weight affect the number and size of peaks in the vertical GRF and stance time. The damping changes the relative size of the peaks but has little influence on stance time. Based on these observations, these parameters were manually tuned at three different speeds to approach experimentally measured vertical GRF and the highest correlation coefficient can reach 0.983. These results indicate that the stiffness, damping, and proper objective functions are all important factors in achieving human-like motion for this simple walking model. These findings can facilitate the understanding of human walking dynamics and may be applied in future biped models.     

Effective leg stiffness in running

Leg stiffness is a common parameter used to characterize leg function during bouncing gaits, like running and hopping. In the literature, different methods to approximate leg stiffness based on kinetic and kinematic parameters are described. A challenging point in estimating leg stiffness is the definition of leg compression during contact. In this paper four methods (methods A–D) based on ground reaction forces (GRF) and one method (method E) relying on temporal parameters are described. Leg stiffness calculated by these five methods is compared with running patterns, predicted by the spring mass model.The best and simplest approximation of leg stiffness is method E. It requires only easily accessible parameters (contact time, flight time, resting leg length, body mass and the leg's touch down angle). Method D is of similar quality but additionally requires the time-dependent progression of the GRF. The other three methods show clear differences from the model predictions by over- or underestimating leg stiffness, especially at slow speeds.Leg stiffness is derived from a conceptual model of legged locomotion and does not exist without this model. Therefore, it is important to prove which experimental method is suited best for approximating the stiffness in a specific task. This will help to interpret the predictions of the conceptual model in comparison with experimental data.     

Changes in running mechanics and spring-mass behavior induced by a mountain ultra-marathon race

Changes in running mechanics and spring-mass behavior due to fatigue induced by a mountain ultra-marathon race (MUM, 166km, total positive and negative elevation of 9500m) were studied in 18 ultra-marathon runners. Mechanical measurements were undertaken pre- and 3h post-MUM at 12km h(-1) on a 7m long pressure walkway: contact (t(c)), aerial (t(a)) times, step frequency (f), and running velocity (v) were sampled and averaged over 5-8 steps. From these variables, spring-mass parameters of peak vertical ground reaction force (F(max)), vertical downward displacement of the center of mass (Δz), leg length change (ΔL), vertical (k(vert)) and leg (k(leg)) stiffness were computed. After the MUM, there was a significant increase in f (5.9±5.5%; P<0.001) associated with reduced t(a) (-18.5±17.4%; P<0.001) with no change in t(c), and a significant decrease in both Δz and F(max) (-11.6±10.5 and -6.3±7.3%, respectively; P<0.001). k(vert) increased by 5.6±11.7% (P=0.053), and k(leg) remained unchanged. These results show that 3h post-MUM, subjects ran with a reduced vertical oscillation of their spring-mass system. This is consistent with (i) previous studies concerning muscular structure/function impairment in running and (ii) the hypothesis that these changes in the running pattern could be associated with lower overall impact (especially during the braking phase) supported by the locomotor system at each step, potentially leading to reduced pain during running.     

How well can spring-mass-like telescoping leg models fit multi-pedal sagittal-plane locomotion data?

Idealized mathematical models of animals, with point-mass bodies and spring-like legs, have been used by researchers to study various aspects of terrestrial legged locomotion. Here, we fit a bipedal spring-mass model to the ground reaction forces of human running, a horse trotting, and a cockroach running. We find that, in all three cases, while the model captures center-of-mass motions and vertical force variations well, horizontal forces are less well reproduced, primarily due to variations in net force vector directions that the model cannot accommodate. The fits result in different apparent leg stiffnesses in the three animals. Assuming a simple fixed leg-angle touch-down strategy, we find that the gaits of these models are stable in different speed-step length regimes that overlap with those used by humans and horses, but not with that used by cockroaches.     

A bipedal compliant walking model generates periodic gait cycles with realistic swing dynamics

A simple spring mechanics model can capture the dynamics of the center of mass (CoM) during human walking, which is coordinated by multiple joints. This simple spring model, however, only describes the CoM during the stance phase, and the mechanics involved in the bipedality of the human gait are limited. In this study, a bipedal spring walking model was proposed to demonstrate the dynamics of bipedal walking, including swing dynamics followed by the step-to-step transition. The model consists of two springs with different stiffnesses and rest lengths representing the stance leg and swing leg. One end of each spring has a foot mass, and the other end is attached to the body mass. To induce a forward swing that matches the gait phase, a torsional hip joint spring was introduced at each leg. To reflect the active knee flexion for foot clearance, the rest length of the swing leg was set shorter than that of the stance leg, generating a discrete elastic restoring force. The number of model parameters was reduced by introducing dependencies among stiffness parameters. The proposed model generates periodic gaits with dynamics-driven step-to-step transitions and realistic swing dynamics. While preserving the mimicry of the CoM and ground reaction force (GRF) data at various gait speeds, the proposed model emulated the kinematics of the swing leg. This result implies that the dynamics of human walking generated by the actuations of multiple body segments is describable by a simple spring mechanics.     

The correlation of segment accelerations and impact forces with knee angle in jump landing

Impact forces and shock deceleration during jumping and running have been associated with various knee injury etiologies. This study investigates the influence of jump height and knee contact angle on peak ground reaction force and segment axial accelerations. Ground reaction force, segment axial acceleration, and knee angles were measured for 6 male subjects during vertical jumping. A simple spring-mass model is used to predict the landing stiffness at impact as a function of (1) jump height, (2) peak impact force, (3) peak tibial axial acceleration, (4) peak thigh axial acceleration, and (5) peak trunk axial acceleration. Using a nonlinear least square fit, a strong (r = 0.86) and significant (p < or = 0.05) correlation was found between knee contact angle and stiffness calculated using the peak impact force and jump height. The same model also showed that the correlation was strong (r = 0.81) and significant (p < or = 0.05) between knee contact angle and stiffness calculated from the peak trunk axial accelerations. The correlation was weaker for the peak thigh (r = 0.71) and tibial (r = 0.45) axial accelerations. Using the peak force but neglecting jump height in the model, produces significantly worse correlation (r = 0.58). It was concluded that knee contact angle significantly influences both peak ground reaction forces and segment accelerations. However, owing to the nonlinear relationship, peak forces and segment accelerations change more rapidly at smaller knee flexion angles (i.e., close to full extension) than at greater knee flexion angles.     

Vertical stiffness during one-legged hopping with and without using a running-specific prosthesis

Although athletes with unilateral below-the-knee amputations (BKAs) generally use their affected leg, including their prosthesis, as their take-off leg for the long jump, little is known about the spring-like leg behavior and stiffness regulation of the affected leg. The purpose of this study was to investigate vertical stiffness during one-legged hopping in an elite-level long jump athlete with a unilateral BKA. We used the spring-mass model to calculate vertical stiffness, which equals the ratio of maximum vertical ground reaction force to maximum center of mass displacement, while the athlete with a BKA hopped on one leg at a range of frequencies. Then, we compared the vertical stiffness of this athlete to seven non-amputee elite-level long-jumpers. We found that from 1.8 to 3.4 Hz, the vertical stiffness of the unaffected leg for an athlete with a BKA increases with faster hopping frequencies, but the vertical stiffness of the affected leg remains nearly constant across frequencies. The athlete with a BKA attained the desired hopping frequencies at 2.2 and 2.6 Hz, but was unable to match the lowest (1.8 Hz) and two highest frequencies (3.0 and 3.4 Hz) using his affected leg. We also found that at 2.5 Hz, unaffected leg vertical stiffness was 15% greater than affected leg vertical stiffness, and the vertical stiffness of non-amputee long-jumpers was 32% greater than the affected leg vertical stiffness of an athlete with a BKA. The results of the present study suggest that the vertical stiffness regulation strategy of an athlete with a unilateral BKA is not the same in the unaffected versus affected legs, and compared to non-amputees.     

Dynamics of the long jump

A mechanical model is proposed which quantitatively describes the dynamics of the centre of gravity (c.g.) during the take-off phase of the long jump. The model entails a minimal but necessary number of components: a linear leg spring with the ability of lengthening to describe the active peak of the force time curve and a distal mass coupled with nonlinear visco-elastic elements to describe the passive peak. The influence of the positions and velocities of the supported body and the jumper's leg as well as of systemic parameters such as leg stiffness and mass distribution on the jumping distance were investigated. Techniques for optimum operation are identified: (1) There is a minimum stiffness for optimum performance. Further increase of the stiffness does not lead to longer jumps. (2) For any given stiffness there is always an optimum angle of attack. (3) The same distance can be achieved by different techniques. (4) The losses due to deceleration of the supporting leg do not result in reduced jumping distance as this deceleration results in a higher vertical momentum. (5) Thus, increasing the touch-down velocity of the jumper's supporting leg increases jumping distance.     

Runners adjust leg stiffness for their first step on a new running surface.

Human runners adjust the stiffness of their stance leg to accommodate surface stiffness during steady state running. This adjustment allows runners to maintain similar center of mass movement (e.g., ground contact time and stride frequency) regardless of surface stiffness. When runners encounter abrupt transitions in the running surface, they must either make a rapid adjustment or allow the change in the surface stiffness to disrupt their running mechanics. Our goal was to determine how quickly runners adjust leg stiffness when they encounter an abrupt but expected change in surface stiffness that they have encountered previously. Six human subjects ran at 3 m s(-1) on a rubber track with two types of rubber surfaces: a compliant "soft" surface (ksurf = 21.3 kN m(-1) and a non-compliant "hard" surface (ksurf = 533 kN m(-1). We found that runners completely adjusted leg stiffness for their first step on the new surface after the transition. For example, runners decreased leg stiffness by 29% between the last step on the soft surface and the first step on the hard surface (from 10.7 kN m(-1) to 7.6 kN m(-1), respectively). As a result, the vertical displacement of the center of mass during stance ( approximately 7 cm) did not change at the transition despite a reduction in surface compression from 6 cm to less than 0.25 cm. By rapidly adjusting leg stiffness, each runner made a smooth transition between surfaces so that the path of the center of mass was unaffected by the change in surface stiffness.     

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.     

Mechanical analysis of the landing phase in heel-toe running.

Results of mechanical analyses of running may be helpful in the search for the etiology of running injuries. In this study a mechanical analysis was made of the landing phase of three trained heel-toe runners, running at their preferred speed and style. The body was modeled as a system of seven linked rigid segments, and the positions of markers defining these segments were monitored using 200 Hz video analysis. Information about the ground reaction force vector was collected using a force plate. Segment kinematics were combined with ground reaction force data for calculation of the net intersegmental forces and moments. The vertical component of the ground reaction force vector Fz was found to reach a first peak approximately 25 ms after touch-down. This peak occurs because, in the support leg, the vertical acceleration of the knee joint is not reduced relative to that of the ankle joint by rotation of the lower leg, so that the support leg segments collide with the floor. Rotation of the support upper leg, however, reduces the vertical acceleration of the hip joint relative to that of the knee joint, and thereby plays an important role in limiting the vertical forces during the first 40 ms. Between 40 and 100 ms after touch-down, the vertical forces are mainly limited by rotation of the support lower leg. At the instant that Fz reaches its first peak, net moments about ankle, knee and hip joints of the support leg are virtually zero. The net moment about the knee joint changed from -100 Nm (flexion) at touch-down to +200 Nm (extension) 50 ms after touch-down. These changes are too rapid to be explained by variations in the muscle activation levels and were ascribed to spring-like behavior of pre-activated knee flexor and knee extensor muscles. These results imply that the runners investigated had no opportunity to control the rotations of body segments during the first part of the contact phase, other than by selecting a certain geometry of the body and muscular (co-)activation levels prior to touch-down.