Human Leg Model Predicts Muscle Forces, States, and Energetics during Walking

Humans employ a high degree of redundancy in joint actuation, with different combinations of muscle and tendon action providing the same net joint torque. Both the resolution of these redundancies and the energetics of such systems depend on the dynamic properties of muscles and tendons, particularly their force-length relations. Current walking models that use stock parameters when simulating muscle-tendon dynamics tend to significantly overestimate metabolic consumption, perhaps because they do not adequately consider the role of elasticity. As an alternative, we posit that the muscle-tendon morphology of the human leg has evolved to maximize the metabolic efficiency of walking at self-selected speed. We use a data-driven approach to evaluate this hypothesis, utilizing kinematic, kinetic, electromyographic (EMG), and metabolic data taken from five participants walking at self-selected speed. The kinematic and kinetic data are used to estimate muscle-tendon lengths, muscle moment arms, and joint moments while the EMG data are used to estimate muscle activations. For each subject we perform an optimization using prescribed skeletal kinematics, varying the parameters that govern the force-length curve of each tendon as well as the strength and optimal fiber length of each muscle while seeking to simultaneously minimize metabolic cost and maximize agreement with the estimated joint moments. We find that the metabolic cost of transport (MCOT) values of our participants may be correctly matched (on average 0.36±0.02 predicted, 0.35±0.02 measured) with acceptable joint torque fidelity through application of a single constraint to the muscle metabolic budget. The associated optimal muscle-tendon parameter sets allow us to estimate the forces and states of individual muscles, resolving redundancies in joint actuation and lending insight into the potential roles and control objectives of the muscles of the leg throughout the gait cycle.

[1]  M L Audu,et al.  A dynamic optimization technique for predicting muscle forces in the swing phase of gait. , 1987, Journal of biomechanics.

[2]  Daniel P. Ferris,et al.  Mechanics and energetics of level walking with powered ankle exoskeletons , 2008, Journal of Experimental Biology.

[3]  R. F. Goldman,et al.  Predicting energy expenditure with loads while standing or walking very slowly. , 1977, Journal of applied physiology: respiratory, environmental and exercise physiology.

[4]  Hartmut Geyer,et al.  A Muscle-Reflex Model That Encodes Principles of Legged Mechanics Produces Human Walking Dynamics and Muscle Activities , 2010, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[5]  M. Pandy,et al.  Dynamic optimization of human walking. , 2001, Journal of biomechanical engineering.

[6]  H. Hatze,et al.  A myocybernetic control model of skeletal muscle , 1977, Biological Cybernetics.

[7]  Kalyanmoy Deb,et al.  Multi-objective optimization using evolutionary algorithms , 2001, Wiley-Interscience series in systems and optimization.

[8]  D W Moran,et al.  A computationally efficient method for solving the redundant problem in biomechanics. , 1995, Journal of biomechanics.

[9]  E L Morin,et al.  Sampling, noise-reduction and amplitude estimation issues in surface electromyography. , 2002, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[10]  O. Schmitt The heat of shortening and the dynamic constants of muscle , 2017 .

[11]  T. Roberts The integrated function of muscles and tendons during locomotion. , 2002, Comparative biochemistry and physiology. Part A, Molecular & integrative physiology.

[12]  S L Delp,et al.  A graphics-based software system to develop and analyze models of musculoskeletal structures. , 1995, Computers in biology and medicine.

[13]  Hartmut Geyer,et al.  Generalization of a muscle-reflex control model to 3D walking , 2013, 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[14]  D. Lloyd,et al.  An EMG-driven musculoskeletal model to estimate muscle forces and knee joint moments in vivo. , 2003, Journal of biomechanics.

[15]  J. Rubenson,et al.  On the ascent: the soleus operating length is conserved to the ascending limb of the force–length curve across gait mechanics in humans , 2012, Journal of Experimental Biology.

[16]  W. Herzog,et al.  The EMG-force relationship of the cat soleus muscle and its association with contractile conditions during locomotion. , 1995, The Journal of experimental biology.

[17]  F. Zajac Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. , 1989, Critical reviews in biomedical engineering.

[18]  Thomas A. McMahon,et al.  Muscles, Reflexes, and Locomotion , 1984 .

[19]  Y. Istefanopulos,et al.  IEEE Engineering in Medicine and Biology Society , 2019, IEEE Transactions on Biomedical Engineering.

[20]  Robert W. Mann,et al.  Myoelectric Signal Processing: Optimal Estimation Applied to Electromyography - Part I: Derivation of the Optimal Myoprocessor , 1980, IEEE Transactions on Biomedical Engineering.

[21]  David G Lloyd,et al.  Estimation of muscle forces and joint moments using a forward-inverse dynamics model. , 2005, Medicine and science in sports and exercise.

[22]  Jack M. Winters,et al.  Multiple Muscle Systems: Biomechanics and Movement Organization , 2011 .

[23]  G. Lichtwark,et al.  Muscle fascicle and series elastic element length changes along the length of the human gastrocnemius during walking and running. , 2007, Journal of biomechanics.

[24]  Emery N. Brown,et al.  Human Leg Model Predicts Ankle Muscle-Tendon Morphology, State, Roles and Energetics in Walking , 2011, PLoS Comput. Biol..

[25]  Edmund Y. S. Chao,et al.  GRAPHICAL INTERPRETATION OF THE SOLUTION TO THE REDUNDANT PROBLEM IN BIOMECHANICS. , 1978 .

[26]  J. Winters Hill-Based Muscle Models: A Systems Engineering Perspective , 1990 .

[27]  M. Alexander,et al.  Principles of Neural Science , 1981 .

[28]  S. Schultz Principles of Neural Science, 4th ed. , 2001 .

[29]  Alan M. Wilson,et al.  Optimal muscle fascicle length and tendon stiffness for maximising gastrocnemius efficiency during human walking and running. , 2008, Journal of theoretical biology.

[30]  J. Brockway Derivation of formulae used to calculate energy expenditure in man. , 1987, Human nutrition. Clinical nutrition.

[31]  Scott L. Delp,et al.  Fibre operating lengths of human lower limb muscles during walking , 2011, Philosophical Transactions of the Royal Society B: Biological Sciences.

[32]  P. Komi,et al.  Muscle-tendon interaction and elastic energy usage in human walking. , 2005, Journal of applied physiology.

[33]  Richard R Neptune,et al.  The effect of walking speed on muscle function and mechanical energetics. , 2008, Gait & posture.

[34]  F. Zajac,et al.  Muscle coordination of maximum-speed pedaling. , 1997, Journal of biomechanics.

[35]  T. Fukunaga,et al.  In vivo behaviour of human muscle tendon during walking , 2001, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[36]  Ken Endo,et al.  A Model of Muscle–Tendon Function in Human Walking at Self-Selected Speed , 2014, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[37]  Neville Hogan,et al.  Myoelectric Signal Processing: Optimal Estimation Applied to Electromyography - Part II: Experimental Demonstration of Optimal Myoprocessor Performance , 1980, IEEE Transactions on Biomedical Engineering.

[38]  Dario Farina,et al.  EMG-Driven Forward-Dynamic Estimation of Muscle Force and Joint Moment about Multiple Degrees of Freedom in the Human Lower Extremity , 2012, PloS one.

[39]  Hartmut Geyer,et al.  Integration of an adaptive swing control into a neuromuscular human walking model , 2013, 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[40]  Gary S Chleboun,et al.  Fascicle length change of the human tibialis anterior and vastus lateralis during walking. , 2007, The Journal of orthopaedic and sports physical therapy.

[41]  Ken Endo,et al.  A quasi-passive model of human leg function in level-ground walking , 2006, 2006 IEEE/RSJ International Conference on Intelligent Robots and Systems.

[42]  A. Hill The efficiency of mechanical power development during muscular shortening and its relation to load , 1964, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[43]  Philip E. Martin,et al.  A Model of Human Muscle Energy Expenditure , 2003, Computer methods in biomechanics and biomedical engineering.

[44]  J. Cram,et al.  Introduction to Surface Electromyography , 1998 .

[45]  Hartmut Geyer,et al.  Regulating speed and generating large speed transitions in a neuromuscular human walking model , 2012, 2012 IEEE International Conference on Robotics and Automation.

[46]  Terence D Sanger,et al.  Bayesian filtering of myoelectric signals. , 2007, Journal of neurophysiology.

[47]  D. Farris,et al.  Human medial gastrocnemius force–velocity behavior shifts with locomotion speed and gait , 2012, Proceedings of the National Academy of Sciences.

[48]  M. Pandy,et al.  In vivo behavior of the human soleus muscle with increasing walking and running speeds. , 2015, Journal of applied physiology.

[49]  M. Latash,et al.  Electromechanical delay: An experimental artifact. , 1992, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[50]  J. He,et al.  Feedback gains for correcting small perturbations to standing posture , 1989, Proceedings of the 28th IEEE Conference on Decision and Control,.

[51]  B. R. Umberger,et al.  Stance and swing phase costs in human walking , 2010, Journal of The Royal Society Interface.

[52]  Vladlen Koltun,et al.  Optimizing locomotion controllers using biologically-based actuators and objectives , 2012, ACM Trans. Graph..

[53]  Samuel R. Hamner,et al.  How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds , 2013, Journal of Experimental Biology.