Research Papers

A Nonlinear Leg Damping Model for the Prediction of Running Forces and Stability

[+] Author and Article Information
Ian Abraham

Department of Mechanical Engineering,
Rutgers University,
Piscataway Township, NJ 08854
e-mail: ianabraham21@gmail.com

ZhuoHua Shen

School of Mechanical Engineering,
Purdue University,
West Lafayette, IN 47907
e-mail: Shen38@purdue.edu

Justin Seipel

School of Mechanical Engineering,
Purdue University,
West Lafayette, IN 47907
e-mail: jseipel@purdue.edu

Manuscript received March 8, 2014; final manuscript received October 2, 2014; published online April 2, 2015. Assoc. Editor: José L. Escalona.

J. Comput. Nonlinear Dynam 10(5), 051008 (Apr 02, 2015) (8 pages) Paper No: CND-14-1068; doi: 10.1115/1.4028751 History: Received March 08, 2014

Despite the neuromechanical complexity underlying animal locomotion, the steady-state center-of-mass motions and ground reaction forces of animal running can be predicted by simple spring-mass models such as the canonical spring-loaded inverted pendulum (SLIP) model. Such SLIP models have been useful for the fields of biomechanics and robotics in part because ground reaction forces are commonly measured and readily available for comparing with model predictions. To better predict the stability of running, beyond the canonical conservative SLIP model, more recent extensions have been proposed and investigated with hip actuation and linear leg damping (e.g., hip-actuated SLIP). So far, these attempts have gained improved prediction of the stability of locomotion but have led to a loss of the ability to accurately predict ground reaction forces. Unfortunately, the linear damping utilized in current models leads to an unrealistic prediction of damping force and ground reaction force with a large nonzero magnitude at touchdown (TD). Here, we develop a leg damping model that is bilinear in leg length and velocity in order to yield improved damping force and ground reaction force prediction. We compare the running ground reaction forces, small and large perturbation stability, parameter sensitivity, and energetic cost resulting from both the linear and bilinear damping models. We found that bilinear damping helps to produce more realistic, smooth vertical ground reaction forces, thus fixing the current problem with the linear damping model. Despite large changes in the damping force and power loss profile during the stance phase, the overall dynamics and energetics on a stride-to-stride basis of the two models are largely the same, implying that the integrated effect of damping over a stride is what matters most to the stability and energetics of running. Overall, this new model, an actuated SLIP model with bilinear damping, can provide significantly improved prediction of ground reaction forces as well as stability and energetics of locomotion.

Copyright © 2015 by ASME
Topics: Stability , Damping
Your Session has timed out. Please sign back in to continue.


Geyer, H. , Seyfarth, A. , and Blickhan, R. , 2005, “Spring-Mass Running: Simple Approximate Solution and Application to Gait Stability,” J. Theor. Biol., 232(3), pp. 315–328.
Blickhan, R. , and Full, R. , 1992, “Similarity in Multilegged Locomotion: Bouncing Like a Monopode,” J. Comp. Physiol., A, 173(5), pp. 509–517.
Full, R. , and Koditschek, D. , 1999, “Templates and Anchors: Neuromechanical Hypotheses of Legged Locomotion on Land,” J. Exp. Biol., 202(Pt 23), pp. 3325–3332. [PubMed]
Schwind, W. , 1998, “Spring Loaded Inverted Pendulum Running: A Plant Model,” Ph.D. thesis, University of Michigan, Ann Arbor, MI.
Ghigliazza, R. M. , Altendorfer, R. , Holmes, P. , and Koditschek, D. , 2003, “A Simply Stabilized Running Model,” SIAM J. Appl. Dyn. Syst., 2(2), pp. 187–218. [CrossRef]
Seyfarth, A. , Geyer, H. , and Herr, H. , 2003, “Swing-Leg Retraction: A Simple Control Model for Stable Running,” J. Exp. Biol., 206, pp. 2547–2555. [CrossRef] [PubMed]
Seyfarth, A. , Geyer, H. , Gunther, M. , and Blickhan, R. , 2002, “A Movement Criterion for Running,” J. Biomech., 35(5), pp. 649–655. [CrossRef] [PubMed]
Shen, Z. H. , and Seipel, J. , 2012, “A Fundamental Mechanism of Legged Locomotion With Hip Torque and Leg Damping,” Bioinspiration Biomimetics, 7(4), p. 046010. [CrossRef] [PubMed]
Seyfarth, A. , Geyer, H. , and Herr, H. , 2003, “Swing-Leg Retraction: A Simple Control Model for Stable Running,” J. Exp. Biol., 206(15), pp. 2547–2555. [CrossRef] [PubMed]
Blum, Y. , Lipfert, S. , Rummel, J. , and Seyfarth, A. , 2010, “Swing Leg Control in Human Running,” Bioinspiration Biomimetics, 5(2), p. 026006. [CrossRef] [PubMed]
Ankarali, M. M. , and Saranli, U. , 2010, “Stride-to-Stride Energy Regulation for Robust Self-Stability of a Torque-Actuated Dissipative Spring-Mass Hopper,” Chaos, 20(3), p. 033121. [CrossRef] [PubMed]
Seipel, J. , and Holmes, P. , 2007, “A Simple Model for Clock-Actuated Legged Locomotion,” Regul. Chaotic Dyn., 12(5), pp. 502–520. [CrossRef]
Saranli, U. , Buehler, M. , and Koditschek, D. , 2001, “RHex: A Simple and Highly Mobile Hexapod Robot,” Int. J. Rob. Res., 20(1), pp. 616–631. [CrossRef]
Holmes, P. , Full, R. , Koditschek, D. , and Guckenheimer, J. , 2006, “The Dynamics of Legged Locomotion: Models, Analyses, and Challenges,” SIAM Rev., 48, pp. 207–304. [CrossRef]
Schmitt, J. , 2006, “A Simple Stabilizing Control for Sagittal Plane Locomotion,” ASME J. Comput. Nonlinear Dyn., 1(4), pp. 348–357. [CrossRef]
Maus, H. M. , Rummel, J. , and Seyfarth, A. , 2008, “Stable Upright Walking and Running Using a Simple Pendulum Based Control Scheme,” CLAWAR, pp. 623–629.
Sreenath, K. , Park, H. W. , Poulakakis, I. , and Grizzle, J. , 2013, “Embedding Active Force Control Within the Compliant Hybrid Zero Dynamics to Achieve Stable, Fast Running on MABEL,” Int. J. Rob. Res.(IJRR), 32(3), pp. 324–345. [CrossRef]
Altendorfer, R. , Koditschek, D. , Komsuoglu, H. , Buehler, M. , Moore, N. , and McMordie, D. , 2001, “Evidence for Spring Loaded Inverted Pendulum Running in a Hexapod Robot,” Experimental Robotics VII, Springer Verlag, Berlin, Heidelberg Germany, pp. 291–302.
Komsuoglu, H. , Majumdar, A. , Aydin, Y. , and Koditschek, D. , 2010, “Characterization of Dynamic Behaviors in a Hexapod Robot,” Experimental Robotics, Springer, Berling, Heidelberg, Germany, pp. 667–684.
Kim, K. Y. , Kwon, O. , Yeon, J. S. , and Park, J. H. , 2006, “Elliptic Trajectory Generation for Galloping Quadruped Robots,” IEEE International Conference on Robotics and Biomimetics, Kunming, Dec. 17–20, pp. 103–108.
Mordatch, I. , de Lasa, M. , and Hertzmann, A. , 2010, “Robust Physics-Based Locomotion Using Low-Dimensional Planning,” ACM Trans. Graphics, 29(3), p. 71.
Keilson, S. E. , Teich, M. C. , and Kahnna, S. M. , 1989, “Models of Nonlinear Vibrations. I. Oscillator With Bilinear Resistance,” Acta Otolaryngol., 108(s467), pp. 241–248. [CrossRef]
Chen, J. T. , and You, D. W. , 1997, “Hysteretic Damping Revisited,” Adv. Eng. Software, 28(3), pp. 165–171. [CrossRef]
de Kraker, B. , 2004, A Numerical-Experimental Approach in Structural Dynamics, Shaker Publishing.
Den Hartog, J. P. , 1985, Mechanical Vibrations, Courier Dover Publications.
Hill, A. V. , 1938, “The Heat of Shortening and the Dynamic Constants of Muscle,” Proc. R. Soc. London, 126(843), pp. 136–195. [CrossRef]
Hill, A. V. , 1949, “The Abrupt Transition From Rest to Activity in Muscle,” Proc. R. Soc. London, 136(884), pp. 399–420. [CrossRef]
Schmidt-Nielsen, K. , 1972, “Locomotion: Energy Cost of Swimming, Flying, and Running,” Science, 177(4045), pp. 222–228. [CrossRef] [PubMed]
McGowan, C. P. , Grabowski, A. M. , McDermott, W. J. , Herr, H. M. , and Kram, R. , 2012, “Leg Stiffness of Sprinters Using Running-Specific Prostheses,” J. R. Soc. Interface, 9(73), pp. 1975–1982. [CrossRef] [PubMed]
Alexander, R. M. , 2005, “Models and the Scaling of Energy Costs for Locomotion,” J. Exp. Biol., 208(9), pp. 1645–1652. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

The actuated-SLIP model. A stride includes stance from TD to LO, and then flight until the next TD. CoM velocity is measured at every TD event.

Grahic Jump Location
Fig. 2

A free body diagram of the (a) main body during stance, (b) leg during stance, (c) main body in flight

Grahic Jump Location
Fig. 3

Here, we show the trajectories of hip-SLIP with linear damping and bilinear damping. The parameters used for hip-SLIP with bilinear damping are τ˜ = 0.40, krel = 20, v˜n = 1.4344, δ = 0.1266, β = 67 deg, and ζ = 10.32 and for hip-SLIP with linear damping, τ˜ = 0.40, krel = 20, v˜n = 1.4344, δ = 0.1449, β = 67 deg, and ζ = 0.4657. The open circles superimposed on the trajectories indicate TD and LO events.

Grahic Jump Location
Fig. 4

Here, we show leg length, leg velocity, damping force, and the power lost through damping for hip-SLIP with both linear and bilinear damping. The parameters used for hip-SLIP with bilinear damping are τ˜ = 0.40, krel = 20, v˜n = 1.4344, δ = 0.1266, β = 67 deg, and ζ = 10.32 and for hip-SLIP with linear damping, τ˜ = 0.40, krel = 20, v˜n = 1.4344, δ = 0.1449, β = 67 deg, and ζ = 0.4657.

Grahic Jump Location
Fig. 5

The vertical ground reaction force of human runners [29] is compared with that of the linear and bilinear cases of damping in the hip-SLIP model. The following human-representative parameters are used for the comparison: τ˜ = 0.4, β = 67 deg, krel = 20, v˜n = 1.4344. For the linear damping case, ζ = 0.4657 and δ = 0.1446. For the bilinear damping case, ζ = 10.32 and δ = 0.1266. In this figure, all three forces are normalized by body weight.

Grahic Jump Location
Fig. 6

Here, we compare how bilinear damping changes the maximum eigenvalue of hip-SLIP for each fixed point solution as hip torque is varied from 0 (SLIP) all the way up to τ = 2, though values above τ = 0.5 are not likely to be commonly seen in animal locomotion. We compare the eigenvalues of hip-SLIP with bilinear damping to that of hip-SLIP with linear damping. We fix the parameters v˜n = 1.4344 and krel = 12.7004, while varying τ˜ and searching for fixed point solutions and solving for the eigenvalues via small perturbation analysis.

Grahic Jump Location
Fig. 7

Here, we look at the effects that bilinear damping has on the basin of attraction for hip-SLIP for four different levels of hip actuation: (a) τ˜ = 0 (equivalent to the canonical SLIP model), (b) τ˜ = 0.10, (c) τ˜ = 0.40, and (d) τ˜ = 1.0. The constant parameters for the (top row) hip-SLIP with linear damping are krel = 12.7004, ζ = 0.3678, and β = 60 deg. The parameters for the (bottom row) hip-SLIP with bilinear damping are krel = 12.7004, ζ = 4.817, and β = 60 deg.

Grahic Jump Location
Fig. 8

Here, we show the speed of all stable periodic hip-SLIP locomotion solutions found for varying levels of torque, for the cases (a) linear leg damping and (b) bilinear leg damping. Fixed point solutions are searched for within the range of δn from −30 deg to the fixed touch down angle β and v˜n from 0 to 5.

Grahic Jump Location
Fig. 9

Here, we show the maximum eigenvalue magnitude of hip-SLIP with (a) linear damping and (b) bilinear damping as relative stiffness krel changes. In the top row, the parameters that are held constant are β = 60 deg and v˜n = 1.4344. Each branch for plots (a) and (b) represents distinct incremental τ˜ values of 0.5, 0.8, and 1.2. In the row below, the parameters that are held constant are τ˜ = 1 and v˜n = 1.4344. Each branch for plots (a) and (b) represents distinct incremental β values of 60, 65, 70 deg.

Grahic Jump Location
Fig. 10

Here, we determine how bilinear damping effects the cost of transport of hip-SLIP. For hip-SLIP, we fixed krel = 20, β = 1.1694, and varied τ˜. The region in gray shows the biologically relevant region for human locomotion.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In