Part I: Leg Muscle Force Estimation in Swing

Model

A. Explore the model and the Forward Dynamics Tool

Objective: Visualize an inverse kinematic solution of the swing phase using the OpenSim GUI. Run a forward dynamic simulation of the swing phase with no muscle control.

Open the leg model, load the motion, and view the animation.

Use the File menu in the OpenSim GUI to open the 6-degree-of-freedom, 9-muscle leg model file (leg6dof9musc.osim) and load the swing motion file (leg69_IK_swing.mot from the Swing folder). The loaded motion will appear under the Motions tree in the Navigator panel. Animate the swing phase motion using the movie controls located above the view pane.

In the Navigator panel, rename the motion from Results to IK so it is distinct from the motions we will be generating below.

Stop the motion and reset the player to the first time step. In the Coordinates tab, set the current pose of the model (i.e., the initial pose of the inverse kinematics (IK) solution) as the default pose by clicking “Poses >” and selecting “Set Default”. Save the model after changing the default parameters.

Use the Forward Dynamics Tool to create a simulation of the system with no muscle controls.

Read the following set of pages on How Forward Dynamics Works.

Open the Forward Dynamics Tool by selecting “Forward Dynamics...” from the Tools menu. Configure the tool as follows:

• In the Input pane, check “Solve for equilibrium for actuator states”, which will equilibrate the fiber and tendon forces in each muscle at the initial time step.
• In the Time pane, specify the time range corresponding to the IK motion (0.117 to 0.617 seconds).
• In the Output pane, create a new directory in the Swing folder (e.g., <YourWorkingDir>\Swing\FWD_No_Controls).
• Save your settings to a file (e.g., leg69_Setup_Forward_No_Controls.xml).
• Run the simulation. Replay the resulting motion using the movie controls.
• Close the tool.
• In the Navigator, rename the motion to FWD_unlocked.

Lock the pelvis coordinates and plot simulation results.

In this forward simulation, the model simply falls because there are no forces to oppose gravity. In the Coordinates tab, lock the pelvis coordinates (pelvis_tilt, pelvis_tx, and pelvis_ty) to prevent the model from falling. Re-run the Forward Dynamics Tool, making sure to use "FWD_locked" in the names of the output folder, the settings file, and the resulting motion.

Open the plotter by selecting "Plot..." from the Tools menu. The X and Y quantities for the plot are selected in the Curve Add panel. Data can be loaded by clicking "Y-Quantity..." and selecting the "Load file..." option (note that the generalized coordinates will appear at the top of the list and the generalized speeds are suffixed with '/speed' ("_u" in versions before 4.0)).

B. Simulate swing phase with manually selected excitations

Objective: Using the Excitation Editor in the OpenSim GUI, find a set of muscle excitations that will match the given swing phase motion.

Make sure that the pelvis coordinates (pelvis_tilt, pelvis_tx, and pelvis_ty) are locked before proceeding with the remaining sections.

In the previous forward dynamic simulation, no controls were specified for the muscle actuators. As a result, no active muscle forces were generated.

Use the Excitation Editor to generate a set of muscle excitations that track the IK solution.

Read the page on how to use the Excitation Editor to learn how to select and change points on a muscle excitation curve.

Two helper files for setting up the excitation editing can be found in the Swing folder: leg69_Forward_Default_Controls.xml and leg69_Forward_Default_States.sto. The first file, leg69_Forward_Default_Controls.xml, adds control information for each muscle in the model. The controller is defined by a set of control nodes (time and value pairs) that are linearly interpolated over time to form the muscle excitation signal. These nodes can be edited using the Excitation Editor. The second file, leg69_Forward_Default_States.sto, contains the initial coordinates, generalized velocities, muscle activations, and muscle fiber lengths that completely specify the starting state used in forward simulations.

Open the Forward Tool and configure the simulation as follows:

• In the Input pane, for Controls, select the XML file that contains the controls of the forward simulation (e.g., leg69_Forward_Default_Controls.xml).
• In the Input pane, for Initial State, select the storage file that contains the initial states of the forward simulation (e.g., leg69_Forward_Default_States.sto).
• Open the Excitation Editor by clicking on the pencil next to the Controls textbox. From the list of muscles, select a muscle capable of generating hip flexion and click OK.
• Increase the excitation of the muscle as desired.
• Save the excitations to a controls file (e.g., leg69_Forward_Controls_<muscle_name>.xml).
• In the Input pane, for Controls, select the newly created XML file.
• Check the “Solve for equilibrium for actuator states” box.
• Check that the time range matches the original time range for IK.
• Specify a new output directory (e.g., Swing\FWD_<muscle_name>).
• Save your settings to a file (e.g., leg69_Setup_Forward_<muscle_name>.xml).
• Run the simulation and compare the values of the hip flexion, knee angle, and ankle angle coordinates to the inverse kinematics solution.
• Close the tool.

Excite three muscles to minimize the error of the hip flexion, knee angle, and ankle angle coordinates.

• By exciting a single muscle in early swing, reduce the error of the hip flexion coordinate in the forward simulation.
• By exciting a second muscle along with your previous excitation, reduce the error of the hip flexion and knee angle coordinates.
• By exciting a third muscle, reduce the error of the hip flexion, knee angle, and ankle angle coordinates.

C. Simulate swing phase with activations from Static Optimization (SO)

Objective: Using the Static Optimization Tool, find a set of muscle controls for the swing phase. Use the resulting controls in a forward dynamic simulation and compare the results.

Use the Static Optimization Tool to find a set of muscle controls that approximately matches the IK motion.

Read the following set of pages on How Static Optimization Works.

Open the Static Optimization Tool and configure it as follows:

• In the Input pane, load the leg69_IK_swing.mot file as the Motion.
• In the Input pane, check the box to filter coordinates and enter a cutoff frequency of 6 Hz.
• In the Objective Function pane, set the muscle activation exponent to 2.0 and check “Use muscle force-length-velocity relation”.
• Set the time range to between 0.117 and 0.617 seconds (if it is not already set).
• Specify the output directory (e.g., Swing\SO).
• Save your settings to a file (e.g., leg69_Setup_SO.xml).
• Make sure the pelvis coordinates (pelvis_tilt, pelvis_tx, and pelvis_ty) are locked.
• Run the tool; close the tool when you are finished.

Use the Forward Dynamic Tool to examine the motion using the estimated muscle controls from Static Optimization.

Run a forward dynamic simulation of the swing phase using the controls generated by Static Optimization. Configure the Forward Dynamics Tool as follows:

• In the Input pane, for Controls, navigate to the folder containing the results from Static Optimization and select the controls file.
• In the Input pane, for Initial State, select the storage file that contains the desired initial states (e.g., leg69_Forward_Default_States.sto).
• Check the “Solve for equilibrium for actuator states” box.
• Specify the same time range as before.
• Specify the output directory (e.g., Swing\FWD_SO).
• Save your settings to a file (e.g., leg69_Setup_Forward_SO.xml).
• Run the tool; close the tool when you are finished.

D. Simulate swing phase with excitations from Computed Muscle Control (CMC)

Objective: Using the CMC Tool, find a set of muscle controls for the swing phase. Use the resulting controls in a forward dynamic simulation and compare the results.

Use the Computed Muscle Control (CMC) Tool to find a set of muscle controls that approximately matches the IK motion.

Read the following set of pages on How Computed Muscle Control Works.

Open the CMC Tool and configure it as follows:

• In the Input pane, specify leg69_IK_swing.mot as the desired kinematics.
• In the Input pane, specify leg69_CMC_Swing_Tracking_Tasks.xml as the tracking task.
• Check the box to filter the coordinates and set the cutoff frequency to 6 Hz.
• Include actuator control constraints that define the muscle control limits (e.g., leg69_muscles_control_limits.xml).
• Set the time range to between 0.117 and 0.617 seconds (if it is not already set).
• Choose a “CMC look-ahead window” that is the approximate time in which muscles can change their output forces in response to a change in input controls (use 0.01s).
• Specify the output directory (e.g., Swing\CMC).
• Save your settings to a file (e.g., leg69_Setup_CMC.xml).
• Make sure the pelvis coordinates (pelvis_tilt, pelvis_tx, and pelvis_ty) are locked.
• Run the tool; close the tool when you are finished.

Use the Forward Dynamics Tool to examine the motion using the estimated muscle controls from CMC.

Run a forward dynamic simulation of the swing phase using the controls generated by CMC. Configure the tool as follows:

• In the Input pane, for Controls, select the controls file generated by CMC.
• In the Input pane, for Initial State, select the states file generated by CMC.
• Check the “Solve for equilibrium for actuator states” box.
• Specify the output directory (e.g., Swing\FWD_CMC).
• Save your settings to a file (e.g., leg69_Setup_Forward_CMC.xml).
• Run the tool; close the tool when you are finished.