Implicit methods for efficient musculoskeletal simulation and optimal control.

The ordinary differential equations for musculoskeletal dynamics are often numerically stiff and highly nonlinear. Consequently, simulations require small time steps, and optimal control problems are slow to solve and have poor convergence. In this paper, we present an implicit formulation of musculoskeletal dynamics, which leads to new numerical methods for simulation and optimal control, with the expectation that we can mitigate some of these problems. A first order Rosenbrock method was developed for solving forward dynamic problems using the implicit formulation. It was used to perform real-time dynamic simulation of a complex shoulder arm system with extreme dynamic stiffness. Simulations had an RMS error of only 0.11 degrees in joint angles when running at real-time speed. For optimal control of musculoskeletal systems, a direct collocation method was developed for implicitly formulated models. The method was applied to predict gait with a prosthetic foot and ankle. Solutions were obtained in well under one hour of computation time and demonstrated how patients may adapt their gait to compensate for limitations of a specific prosthetic limb design. The optimal control method was also applied to a state estimation problem in sports biomechanics, where forces during skiing were estimated from noisy and incomplete kinematic data. Using a full musculoskeletal dynamics model for state estimation had the additional advantage that forward dynamic simulations, could be done with the same implicitly formulated model to simulate injuries and perturbation responses. While these methods are powerful and allow solution of previously intractable problems, there are still considerable numerical challenges, especially related to the convergence of gradient-based solvers.

[1]  Benno M. Nigg,et al.  A Video Technique for Obtaining 3-D Coordinates in Alpine Skiing , 1996 .

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

[3]  K N An,et al.  Determination of muscle orientations and moment arms. , 1984, Journal of biomechanical engineering.

[4]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[5]  Dan Simon,et al.  Optimal State Estimation: Kalman, H∞, and Nonlinear Approaches , 2006 .

[6]  F.C.T. van der Helm,et al.  A finite element musculoskeletal model of the shoulder mechanism. , 1994 .

[7]  Antonie J. van den Bogert,et al.  Analysis and simulation of mechanical loads on the human musculoskeletal system: a methodological overview. , 1994 .

[8]  Antonie J. van den Bogert,et al.  A Real-Time, 3-D Musculoskeletal Model for Dynamic Simulation of Arm Movements , 2009, IEEE Transactions on Biomedical Engineering.

[9]  S. McLean,et al.  Development and validation of a 3-D model to predict knee joint loading during dynamic movement. , 2003, Journal of biomechanical engineering.

[10]  F. V. D. van der Helm A finite element musculoskeletal model of the shoulder mechanism. , 1994, Journal of biomechanics.

[11]  S. McLean,et al.  Sagittal plane biomechanics cannot injure the ACL during sidestep cutting. , 2004, Clinical biomechanics.

[12]  Andrew P. Witkin,et al.  Spacetime constraints , 1988, SIGGRAPH.

[13]  R. Crowninshield,et al.  THE PREDICTION OF FORCES IN JOINT STRUCTURES: DISTRIBUTION OF INTERSEGMENTAL RESULTANTS , 1981, Exercise and sport sciences reviews.

[14]  Marko Ackermann,et al.  Optimality principles for model-based prediction of human gait. , 2010, Journal of biomechanics.

[15]  Richard R Neptune,et al.  Biomechanics and muscle coordination of human walking: part II: lessons from dynamical simulations and clinical implications. , 2003, Gait & posture.

[16]  C. Spoor,et al.  Measuring muscle and joint geometry parameters of a shoulder for modeling purposes. , 1999, Journal of biomechanics.

[17]  John T. Betts,et al.  Practical Methods for Optimal Control and Estimation Using Nonlinear Programming , 2009 .

[18]  E. Otten Concepts and Models of Functional Architecture in Skeletal Muscle , 1988, Exercise and sport sciences reviews.

[19]  Manindra R. Kaphle,et al.  Optimality in forward dynamics simulations. , 2008, Journal of biomechanics.

[20]  R. Wyn Morton,et al.  Basic Biomechanics of the Musculoskeletal System (2nd ed.) , 1990 .

[21]  R. Brand,et al.  The biomechanics and motor control of human gait: Normal, elderly, and pathological , 1992 .

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

[23]  Walter Herzog,et al.  Model-based estimation of muscle forces exerted during movements. , 2007, Clinical biomechanics.

[24]  O. von Stryk,et al.  Efficient forward dynamics simulation and optimization of human body dynamics , 2006 .

[25]  M. Roche,et al.  Rosenbrock methods for Differential Algebraic Equations , 1987 .

[26]  C. Spoor,et al.  Knee muscle moment arms from MRI and from tendon travel. , 1992, Journal of biomechanics.

[27]  Anil V. Rao,et al.  Practical Methods for Optimal Control Using Nonlinear Programming , 1987 .

[28]  A J van den Bogert,et al.  Multi-body modelling and simulation of animal locomotion. , 1993, Acta anatomica.

[29]  Christopher L Vaughan,et al.  Theories of bipedal walking: an odyssey. , 2003, Journal of biomechanics.

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

[31]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[32]  H. Hatze The complete optimization of a human motion , 1976 .

[33]  Richard A. Brand,et al.  The biomechanics and motor control of human gait: Normal, elderly, and pathological , 1992 .

[34]  Jeffrey A. Reinbolt,et al.  From the Selectedworks of Jeffrey A. Reinbolt Design of Patient-specific Gait Modifications for Knee Osteoarthritis Rehabilitation Design of Patient-specific Gait Modifications for Knee Osteoarthritis Rehabilitation , 2022 .

[35]  R R Neptune,et al.  Optimization algorithm performance in determining optimal controls in human movement analyses. , 1999, Journal of biomechanical engineering.

[36]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[37]  W. Sellers,et al.  Stride lengths, speed and energy costs in walking of Australopithecus afarensis: using evolutionary robotics to predict locomotion of early human ancestors , 2005, Journal of The Royal Society Interface.

[38]  Ayman Habib,et al.  OpenSim: Open-Source Software to Create and Analyze Dynamic Simulations of Movement , 2007, IEEE Transactions on Biomedical Engineering.

[39]  J. Heegaard,et al.  Predictive algorithms for neuromuscular control of human locomotion. , 2001, Journal of biomechanics.