Two efficient static optimization algorithms that account for muscle-tendon equilibrium: approaching the constraint Jacobian via a constant or a cubic spline function. Issue 11 (17th August 2020)

• Record Type: Journal Article
• Authors: Michaud, Benjamin; Begon, Mickaël
• Is Part Of: Computer methods in biomechanics and biomedical engineering. Volume 23:Issue 11(2020)
• Page Start: 703
• Page End: 709
• Publication Date: 2020-08-17
• Subjects: Static optimization -- linearizing constraint -- constraint Jacobian; Biomechanics -- Data processing -- Periodicals; Biomedical engineering -- Periodicals; Biomechanics -- Periodicals; Biomedical Engineering -- methods -- Periodicals; Computing Methodologies -- Periodicals; 612.7
• DOI: 10.1080/10255842.2020.1759042
• Languages: English
• ISSNs: 1025-5842

Abstract

Abstract: Despite recent advances in algorithms for estimating muscle activities, static optimization remains the most used. Static optimization estimates muscle activations required to obtain a particular set of estimated kinematics. Although fast, static optimization may require considerable time for long trials. Improvements have been proposed in the past, but the current implementations are either accurate and slow (such as the most traditional implementation) or fast but less accurate (such as the linearized at maximal activations method used by OpenSim ). Two innovative algorithms are proposed to improve both optimization time and accuracy of the static optimization. The first, designed to be fast, linearizes the constraint—i.e., constructing a constant constraint Jacobian—at the activations of the previous frame. The second, designed to be as accurate as possible, approaches the constraint Jacobian by cubic splines. Their performance and accuracy are compared to the traditional and OpenSim implementations. The linearized method performed as fast as the OpenSim implementation and was more accurate (0.3% of RMSE versus 5.9%). The spline method had excellent accuracy (0.1% of RMSE), but was 2 X slower than the linearized approaches. Nevertheless, it was 100 X faster than the traditional implementation. Our linearized method is therefore recommended when fast computation is needed, such as real-time applications, while the spline method is recommended otherwise.