Table of Contents | ||
---|---|---|
|
Info |
---|
The tutorial below is designed for use with OpenSim version 4.0 and later. A version of the tutorial compatible with OpenSim version 3.3 is available here. |
I. Objectives
Purpose
The purpose of this tutorial is to demonstrate how OpenSim solves an inverse kinematics and dynamics problem using experimental data. To diagnose movement disorders and study human movement, biomechanists frequently ask human subjects to perform movements in a motion capture laboratory and use computational tools to analyze these movements. A common step in analyzing a movement is to compute the joint angles and joint moments of the subject during movement. OpenSim has tools for computing these quantities:
...
Inverse kinematics computes the joint angles for a musculoskeletal model that best reproduce the motion of a subject. Inverse dynamics then uses joint use joint angles, angular velocities, and angular accelerations of the model, together with the experimental ground reaction forces and moments, to solve for the net reaction forces and net moments at each of the joints. The schematic below shows an overview of the inverse kinematics and inverse dynamics problems.
In this tutorial, you will:
...
Each section of the tutorial guides you in using certain tools within and asks you to answer a few questions. The menu titles and option names you must select and any commands you must type to run OpenSim will appear in bold face boldface. The questions can be answered based on information from OpenSim and basic knowledge of the human musculoskeletal system. After you complete the tutorial, feel free to explore OpenSim and the other analysis tools further on your own. Depending on the amount of exploration you do, this tutorial should take 1-2 hours to complete.
...
III. Scaling A Musculoskeletal Model
Subject-specific modeling involves (i) scaling a generic musculoskeletal model to modify the anthropometry, or physical dimensions, of the generic model so that it matches the anthropometry of a particular subject and (ii) registering the markers placed on the model to match the locations on the subject. Scaling and Registration are the most important steps in solving inverse kinematics and inverse dynamics problems because IK and ID solutions are sensitive to the accuracy of the scaling and registration. To scale the generic model and register the markers:
This Scale Setup file is an xml file that contains pre-configured settings to scale and register the generic gait2354 musculoskeletal model to the dimensions of a subject that we have experimental data for. A detailed explanation of the Scale Tool can be found on the Scaling page of the documentation. | |
Model ScalingIn OpenSim, the scaling step adjusts both the dimensions of the body segments, as well as the mass properties (mass and inertia tensor). Scaling can be performed using a combination of two methods: (1) Manual Scaling: Scaling that allows the user to scale a segment based on some predetermined scale factor. Manual scaling is sometimes necessary when suitable data are not available, or if the scale factors were determined using an alternative algorithm. (2) Measurement-based Scaling: Scaling that determines scale factors for a body segment by comparing distance measurements between specified landmarks on the model, known as model markers, and the corresponding experimental marker positions.
Marker RegistrationIn , OpenSim, the registration step adjusts the location of model markers to match the location of markers on the subject. To do this, you must first estimate a pose for the model that closely resembles the pose of the subject during the experimental static trial. | |
To complete the scale step:
|
...
where q is the vector of generalized coordinates (e.g., joint angles), xiexp is the position of experimental marker i, xi(q) is the position of the corresponding model marker i (which depends on q), and wi is the weight associated with marker i.
To setup set up an inverse kinematics analysis:
subject01_Setup_IK.xml contains pre-configured settings for the inverse kinematics tool. Notice the text boxes in the dialog window are now are filled with values. A detailed explanation of the Inverse Kinematics Tool can be found on the Inverse Kinematics page of the documentation.
| |
Navigate to the Weights tab.
| |
To perform inverse kinematics:
|
...
4. In the inverse kinematics Tool dialog window, click the Weights tab and scroll through the list of markers in the top half of the weights tab. Which markers have weighting values less than one? Why?
Hint: Think about joints that have not been modelledmodeled.
5. Based on information in the Messages window, what is the root-mean-squared (RMS) error of all the markers in the last frame of the motion? Include units. Does this seem reasonable? Explain.
...
To setup an inverse dynamics analysis:
A detailed explanation of the Inverse Dynamics Tool can be found on the Inverse Dynamics page of the documentation. | |
To perform inverse dynamics:
| |
When completed, examine the results of the inverse dynamics solution by plotting the net moments at the left and right ankles:
| |
In solving the inverse dynamics problem, both kinematic data and force plate data were used, making this an over-determined problem. In other words, the problem has more equations than unknowns (i.e., degrees of freedom). Due to errors in the experimental motion data and inaccuracies in the musculoskeletal model, it turns out that Newton's second law is violated, or [3]. One method to handle this inconsistency is to compute and apply residual forces and moments to a particular body segment in the model, such that Newton's second law becomes: † An analogous equation relates the ground reaction moment, to the residual moment, . In this musculoskeletal model, the residuals are applied to the pelvis segment. To see the residuals from the inverse dynamics solution, in a new plot window, Plot pelvis_tx_force, pelvis_ty_force, and pelvis_tz_force versus time. Using this plot, answer question 9. While applying residual forces and moments makes the model's motion dynamically consistent with the external forces , this strategy is undesirable because the residuals can be large. More advanced strategies have been developed to deal with the problem of residuals and dynamic inconsistencies, such as least-squares optimization [3], the Residual Elimination Algorithm (REA) [5], and the Residual Reduction Algorithm (RRA) [6]. OpenSim implements a Residual Reduction Algorithm as part of its workflow for generating muscle-actuated simulations [6]. A detailed explanation of the Residual Reduction Algorithm (RRA) can be found on the Residual Reduction Algorithm page of the documentation. For additional information on these strategies, please also see [3], [5], [6], and [7]. |
...
9. What are the maximum magnitudes of the residual forces? Using the mass of the subject from Question 1, what fraction of body weight are the maximum residual forces?
...
Acknowledgments
The experimental gait data were collected by Jill Higginson and Chand John in the Neuromuscular Biomechanics Lab at the University of Delaware [8]. The data include marker trajectories and ground reaction forces for an adult male walking at a self-selected speed on an instrumented split-belt treadmill. Please note that the data distributed with OpenSim is from a different subject than the one described in the paper. Data collection protocols were the same for both subjects.
...