title: Generating optimal control simulations of musculoskeletal movement using OpenSim and MATLAB.
abstract: Computer modeling, simulation and optimization are powerful tools that have seen increased use in biomechanics research. Dynamic optimizations can be categorized as either data-tracking or predictive problems. The data-tracking approach has been used extensively to address human movement problems of clinical relevance. The predictive approach also holds great promise, but has seen limited use in clinical applications. Enhanced software tools would facilitate the application of predictive musculoskeletal simulations to clinically-relevant research. The open-source software OpenSim provides tools for generating tracking simulations but not predictive simulations. However, OpenSim includes an extensive application programming interface that permits extending its capabilities with scripting languages such as MATLAB. In the work presented here, we combine the computational tools provided by MATLAB with the musculoskeletal modeling capabilities of OpenSim to create a framework for generating predictive simulations of musculoskeletal movement based on direct collocation optimal control techniques. In many cases, the direct collocation approach can be used to solve optimal control problems considerably faster than traditional shooting methods. Cyclical and discrete movement problems were solved using a simple 1 degree of freedom musculoskeletal model and a model of the human lower limb, respectively. The problems could be solved in reasonable amounts of time (several seconds to 1-2 hours) using the open-source IPOPT solver. The problems could also be solved using the fmincon solver that is included with MATLAB, but the computation times were excessively long for all but the smallest of problems. The performance advantage for IPOPT was derived primarily by exploiting sparsity in the constraints Jacobian. The framework presented here provides a powerful and flexible approach for generating optimal control simulations of musculoskeletal movement using OpenSim and MATLAB. This should allow researchers to more readily use predictive simulation as a tool to address clinical conditions that limit human mobility.
authors.affiliation: Department of Kinesiology, University of Massachusetts Amherst , Amherst, MA , United States.
authors.forename: Leng-Feng
authors.initials: LF
authors.lastname: Lee
authors.affiliation: Department of Kinesiology, University of Massachusetts Amherst , Amherst, MA , United States.
authors.forename: Brian R
authors.initials: BR
authors.lastname: Umberger
keywords: Dynamics; Musculoskeletal model; Optimization; Predictive simulation
references.citation: Procedia IUTAM. 2011;2:212-232
references.pmid: 25893160
references.citation: Gait Posture. 2016 Jan;43:24-30
references.pmid: 26669947
references.citation: J Biomech. 2010 Apr 19;43(6):1055-60
references.pmid: 20074736
references.citation: J Biomech. 2001 Aug;34(8):1077-83
references.pmid: 11448699
references.citation: J Biomech. 2012 May 11;45(8):1517-21
references.pmid: 22464351
references.citation: J Biomech Eng. 2012 Nov;134(11):111005
references.pmid: 23387787
references.citation: J Neurophysiol. 2014 Oct 15;112(8):1815-24
references.pmid: 24944215
references.citation: IEEE Trans Neural Syst Rehabil Eng. 2010 Jun;18(3):263-73
references.pmid: 20378480
references.citation: Comput Methods Biomech Biomed Engin. 2015 ;18(13):1445-58
references.pmid: 24995410
references.citation: J Biomech Eng. 1992 Nov;114(4):450-60
references.pmid: 1487896
references.citation: PLoS One. 2015 Apr 01;10(4):e0121407
references.pmid: 25830913
references.citation: Stud Health Technol Inform. 2011;163:156-62
references.pmid: 21335781
references.citation: J Biomech. 2006;39(10):1769-77
references.pmid: 16046223
references.citation: J Biomech. 2003 Aug;36(8):1111-6
references.pmid: 12831736
references.citation: J Biomech. 1987;20(2):187-201
references.pmid: 3571299
references.citation: Procedia IUTAM. 2011 Jan 1;2(2011):297-316
references.pmid: 22102983
references.citation: Proc Biol Sci. 2012 Apr 22;279(1733):1498-505
references.pmid: 22072601
references.citation: Procedia IUTAM. 2011;2:241-261
references.pmid: 25866705
references.citation: J Biomech. 2001 Nov;34(11):1387-98
references.pmid: 11672713
references.citation: IEEE Trans Biomed Eng. 2007 Nov;54(11):1940-50
references.pmid: 18018689
references.citation: IEEE Trans Biomed Eng. 1990 Aug;37(8):757-67
references.pmid: 2210784
references.citation: J Biomech. 2006;39(6):1107-15
references.pmid: 16023125
references.citation: Annu Rev Biomed Eng. 2001;3:245-73
references.pmid: 11447064
references.citation: J Biomech Eng. 2001 Oct;123(5):381-90
references.pmid: 11601721
references.citation: J R Soc Interface. 2010 Sep 6;7(50):1329-40
references.pmid: 20356877