skip to main content

Title: Evaluation of Direct Collocation Optimal Control Problem Formulations for Solving the Muscle Redundancy Problem
Estimation of muscle forces during motion involves solving an indeterminate problem (more unknown muscle forces than joint moment constraints), frequently via optimization methods. When the dynamics of muscle activation and contraction are modeled for consistency with muscle physiology, the resulting optimization problem is dynamic and challenging to solve. This study sought to identify a robust and computationally efficient formulation for solving these dynamic optimization problems using direct collocation optimal control methods. Four problem formulations were investigated for walking based on both a two and three dimensional model. Formulations differed in the use of either an explicit or implicit representation of contraction dynamics with either muscle length or tendon force as a state variable. The implicit representations introduced additional controls defined as the time derivatives of the states, allowing the nonlinear equations describing contraction dynamics to be imposed as algebraic path constraints, simplifying their evaluation. Problem formulation affected computational speed and robustness to the initial guess. The formulation that used explicit contraction dynamics with muscle length as a state failed to converge in most cases. In contrast, the two formulations that used implicit contraction dynamics converged to an optimal solution in all cases for all initial guesses, with tendon force as a state generally being the fastest. Future work should focus on comparing the present approach to other approaches for computing muscle forces. The present approach lacks some of the major limitations of established methods such as static optimization and computed muscle control while remaining computationally efficient.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Annals of Biomedical Engineering
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Continuum robots suffer large deflections due to internal and external forces. Accurate modeling of their passive compliance is necessary for accurate environmental interaction, especially in scenarios where direct force sensing is not practical. This paper focuses on deriving analytic formulations for the compliance of continuum robots that can be modeled as Kirchhoff rods. Compared to prior works, the approach presented herein is not subject to the constant-curvature assumptions to derive the configuration space compliance, and we do not rely on computationally-expensive finite difference approximations to obtain the task space compliance. Using modal approximations over curvature space and Lie group integration, we obtain closed-form expressions for the task and configuration space compliance matrices of continuum robots, thereby bridging the gap between constant-curvature analytic formulations of configuration space compliance and variable curvature task space compliance. We first present an analytic expression for the compliance of aingle Kirchhoff rod.We then extend this formulation for computing both the task space and configuration space compliance of a tendon-actuated continuum robot. We then use our formulation to study the tradeoffs between computation cost and modeling accuracy as well as the loss in accuracy from neglecting the Jacobian derivative term in the compliance model. Finally, we experimentally validate the model on a tendon-actuated continuum segment, demonstrating the model’s ability to predict passive deflections with error below 11.5% percent of total arc length. 
    more » « less
  2. null (Ed.)
    Optimal Transport (OT) distances such as Wasserstein have been used in several areas such as GANs and domain adaptation. OT, however, is very sensitive to outliers (samples with large noise) in the data since in its objective function, every sample, including outliers, is weighed similarly due to the marginal constraints. To remedy this issue, robust formulations of OT with unbalanced marginal constraints have previously been proposed. However, employing these methods in deep learning problems such as GANs and domain adaptation is challenging due to the instability of their dual optimization solvers. In this paper, we resolve these issues by deriving a computationally-efficient dual form of the robust OT optimization that is amenable to modern deep learning applications. We demonstrate the effectiveness of our formulation in two applications of GANs and domain adaptation. Our approach can train state-of-the-art GAN models on noisy datasets corrupted with outlier distributions. In particular, our optimization computes weights for training samples reflecting how difficult it is for those samples to be generated in the model. In domain adaptation, our robust OT formulation leads to improved accuracy compared to the standard adversarial adaptation methods. 
    more » « less
  3. Abstract Introduction

    Optimal frequency modulation during functional electrical stimulation (FES) may minimize or delay the onset of FES‐induced muscle fatigue.


    An offline dynamic optimization method, constrained to a modified Hill‐Huxley model, was used to determine the minimum number of pulses that would maintain a constant desired isometric contraction force.


    Six able‐bodied participants were recruited for the experiments, and their quadriceps muscles were stimulated while they sat on a leg extension machine. The force–time (F–T) integrals and peak forces after the pulse train was delivered were found to be statistically significantly greater than the force–time integrals and peak forces obtained after a constant frequency train was delivered.


    Experimental results indicated that the optimized pulse trains induced lower levels of muscle fatigue compared with constant frequency pulse trains. This could have a potential advantage over current FES methods that often choose a constant frequency stimulation train.Muscle Nerve57: 634–641, 2018

    more » « less
  4. Abstract Cerebrovascular accidents like a stroke can affect the lower limb as well as upper extremity joints (i.e., shoulder, elbow, or wrist) and hinder the ability to produce necessary torque for activities of daily living. In such cases, muscles’ ability to generate forces reduces, thus affecting the joint’s torque production. Understanding how muscles generate forces is a key element to injury detection. Researchers have developed several computational methods to obtain muscle forces and joint torques. Electromyography (EMG) driven modeling is one of the approaches to estimate muscle forces and obtain joint torques from muscle activity measurements. Musculoskeletal models and EMG-driven models require necessary muscle-specific parameters for the calculation. The focus of this study is to investigate the EMG-driven approach along with an upper extremity musculoskeletal model to determine muscle forces of two major muscle groups, biceps brachii and triceps brachii, consisting of seven muscle-tendon units. Estimated muscle forces are used to determine the elbow joint torque. Experimental EMG signals and motion capture data are collected for a healthy subject. The musculoskeletal model is scaled to match the geometric parameters of the subject. Then, the approach calculates muscle forces and joint moment for two tasks: simple elbow flexion extension and triceps kickback. Individual muscle forces and net joint torques for both tasks are estimated. The study also has compared the effect of muscle-tendon parameters (optimal fiber length and tendon slack length) on the estimated results. 
    more » « less
  5. There has been an increasing interest in the use of autonomous underwater robots to monitor freshwater and marine environments. In particular, robots that propel and maneuver themselves like fish, often known as robotic fish, have emerged as mobile sensing platforms for aquatic environments. Highly nonlinear and often under-actuated dynamics of robotic fish present significant challenges in control of these robots. In this work, we propose a nonlinear model predictive control (NMPC) approach to path-following of a tail-actuated robotic fish that accommodates the nonlinear dynamics and actuation constraints while minimizing the control effort. Considering the cyclic nature of tail actuation, the control design is based on an averaged dynamic model, where the hydrodynamic force generated by tail beating is captured using Lighthill's large-amplitude elongated-body theory. A computationally efficient approach is developed to identify the model parameters based on the measured swimming and turning data for the robot. With the tail beat frequency fixed, the bias and amplitude of the tail oscillation are treated as physical variables to be manipulated, which are related to the control inputs via a nonlinear map. A control projection method is introduced to accommodate the sector-shaped constraints of the control inputs while minimizing the optimization complexity in solving the NMPC problem. Both simulation and experimental results support the efficacy of the proposed approach. In particular, the advantages of the control projection method are shown via comparison with alternative approaches. 
    more » « less