Forward dynamic simulation provides a powerful framework for characterizing internal loads and for predicting changes in movement due to injury, impairment or surgical intervention. However, the computational challenge of generating simulations has greatly limited the use and application of forward dynamic models for simulating human gait. In this study, we introduce an optimal estimation approach to efficiently solve for generalized accelerations that satisfy the overall equations of motion and best agree with measured kinematics and ground reaction forces. The estimated accelerations are numerically integrated to enforce dynamic consistency over time, resulting in a forward dynamic simulation. Numerical optimization is then used to determine a set of initial generalized coordinates and speeds that produce a simulation that is most consistent with the measured motion over a full cycle of gait. The proposed method was evaluated with synthetically created kinematics and force plate data in which both random noise and bias errors were introduced. We also applied the method to experimental gait data collected from five young healthy adults walking at a preferred speed. We show that the proposed residual elimination algorithm (REA) converges to an accurate solution, reduces the detrimental effects of kinematic measurement errors on joint moments, and eliminates the need for residual forces that arise in standard inverse dynamics. The greatest improvements in joint kinetics were observed proximally, with the algorithm reducing joint moment errors due to marker noise by over 20% at the hip and over 50% at the low back. Simulated joint angles were generally within $1deg$ of recorded values when REA was used to generate a simulation from experimental gait data. REA can thus be used as a basis for generating accurate simulations of subject-specific gait dynamics.

1.
Zajac
,
F. E.
,
Neptune
,
R. R.
, and
Kautz
,
S. A.
, 2003, “
Biomechanics and Muscle Coordination of Human Walking. Part II. Lessons From Dynamical Simulations and Clinical Implications
,”
Gait and Posture
0966-6362,
14
, pp.
1
17
.
2.
Delp
,
S. L.
,
Anderson
,
F. C.
,
Arnold
,
A. S.
,
,
P.
,
Habib
,
A.
,
Chand
,
J.
,
Guendelman
,
E.
, and
Thelen
,
D. G.
, 2007, “
Opensim: Open-Source Software to Create and Analyze Dynamic Simulations of Movement
,”
IEEE Trans. Biomed. Eng.
0018-9294,
54
, pp.
1940
1950
.
3.
Winter
,
D.
, 1990, “
Kinetics: Forces and Moments of Force
,”
Biomechanics and Motor Control of Human Movement
,
Wiley-Interscience
,
New York
, pp.
75
102
.
4.
Cahout
,
V.
,
Luc
,
M.
, and
David
,
A.
, 2002, “
Static Optimal Estimation of Joint Accelerations for Inverse Dynamics Problem Solution
,”
J. Biomech.
0021-9290,
35
, pp.
1507
1513
.
5.
Kuo
,
A.
, 1998, “
A Least-Squares Estimation Approach to Improving the Precision of Inverse Dynamics Computations
,”
ASME J. Biomech. Eng.
0148-0731,
120
, pp.
148
159
.
6.
van den Bogert
,
A. J.
, and
Su
,
A.
, 2008, “
A Weighted Least Squares Method for Inverse Dynamic Analysis
,”
Comput. Methods Biomech. Biomed. Eng.
1025-5842,
11
, pp.
3
9
.
7.
Chao
,
E.
, and
Rim
,
K.
, 1973, “
Application of Optimization Principles in Determining the Applied Moments in Human Leg Joints During Gait
,”
J. Biomech.
0021-9290,
6
, pp.
497
510
.
8.
Anderson
,
F.
, and
Pandy
,
M.
, 2001, “
Dynamic Optimization of Human Walking
,”
ASME J. Biomech. Eng.
0148-0731,
135
(
5
), pp.
381
390
.
9.
Neptune
,
R.
,
Kautz
,
S.
, and
Zajec
,
F.
, 2001, “
Contributions of the Individual Ankle Plantar Flexors to Support, Forward Progression and Swing Initiation During Walking
,”
J. Biomech.
0021-9290,
34
, pp.
1387
1398
.
10.
Greenwood
,
D.
, 1988, “
Dynamics of a Rigid Body
,”
Principles of Dynamics
,
Prentice Hall
,
Englewood Cliffs, NJ
, pp.
389
468
.
11.
Yamaguchi
,
G. T.
, 2002, “
Properties of the Pseudoinverse
,”
Dynamic Modeling of Musculoskeletal Motion
,
,
Dordrecht
, pp.
225
226
.
12.
Thelen
,
D.
, and
Anderson
,
F.
, 2006, “
Using Computed Muscle Control to Generate Forward Dynamic Simulations of Human Walking From Experimental Data
,”
J. Biomech.
0021-9290,
39
, pp.
1107
1115
.
13.
McConville
,
J. T.
,
Churchill
,
T. D.
,
Kaleps
,
I.
,
Clauser
,
C. E.
, and
Cuzzi
,
J.
, 1980, “
Anthropometric Relationships of Body and Body Segment Moments of Inertia
,”
Aerospace Medical Research Laboratory, Wright-Patterson Air Force Base
, Technical Report No. AMRL-TR-80-119.
14.
Lu
,
T. W.
, and
O’Connor
,
J. J.
, 1999, “
Bone Position Estimation From Skin Marker Coordinates Using Global Optimisation With Joint Constraints
,”
J. Biomech.
0021-9290,
32
, pp.
129
134
.
15.
Seth
,
A.
, and
Pandy
,
M.
, 2007, “
A Neuromusculoskeletal Tracking Method for Estimating Individual Muscle Forces in Human Movement
,”
J. Biomech.
0021-9290,
40
, pp.
356
366
.
16.
Siegel
,
K. L.
,
Kepple
,
T. M.
, and
Stanhope
,
S. J.
, 2006, “
Using Induced Accelerations to Understand Knee Stability During Gait of Individuals With Muscle Weakness
,”
Gait and Posture
0966-6362,
23
, pp.
435
440
.
17.
Cappozzo
,
A.
,
Leo
,
T.
, and
Pedotti
,
A.
, 1975, “
A General Computation Method for the Analysis of Human Locomotion
,”
J. Biomech.
0021-9290,
8
, pp.
307
320
.
18.
Challis
,
J. H.
, and
Kerwin
,
D. G.
, 1996, “
Quantification of the Uncertainties in Resultant Joint Moments Computed in a Dynamic Activity
,”
J. Sports Sci.
0264-0414,
14
, pp.
219
231
.
19.
Vaughan
,
C.
,
Andrews
,
J.
, and
Hay
,
J.
, 1982, “
Selection of Body Segment Parameters by Optimization Methods
,”
ASME J. Biomech. Eng.
0148-0731,
104
, pp.
38
44
.
20.
Reinbolt
,
J. A.
,
Schutte
,
B. J.
,
Fregly
,
B. J.
,
Koh
,
B. I.
,
George
,
A. D.
, and
Mitchell
,
K. H.
, 2005, “
Determination of Patient-Specific Multi-Joint Kinematic Models Through Two-Level Optimization
,”
J. Biomech.
0021-9290,
38
, pp.
621
626
.
21.
Reinbolt
,
J. A.
,
Haftka
,
R. T.
,
Chmielewski
,
T. L.
, and
Fregly
,
B. J.
, 2007, “
Are Patient-Specific Joint and Inertial Parameters Necessary for Accurate Inverse Dynamics Analyses of Gait?
,”
IEEE Trans. Biomed. Eng.
0018-9294,
54
, pp.
782
793
.