Multibody Muscle Driven Model of an Instrumented Prosthetic Knee During Squat and Toe Rise Motions

[+] Author and Article Information
Antonis P. Stylianou

Research Associate
e-mail: stylianoua@umkc.edu

Trent M. Guess

Associate Professor
e-mail: guesstr@umkc.edu

Mohammad Kia

Research Associate
e-mail: kiam@umkc.edu
Department of Civil & Mechanical Engineering,
Musculoskeletal Biomechanics Research Lab,
University of Missouri—Kansas City,
350K Robert H Flarsheim Hall,
5110 Rockhill Rd.,
Kansas City, MO 64110

1Corresponding author.

Contributed by the Bioengineering Division of ASME for publication in the JOURNAL OF BIOMECHANICAL ENGINEERING. Manuscript received August 22, 2012; final manuscript received March 1, 2013; accepted manuscript posted March 8, 2013; published online April 5, 2013. Assoc. Editor: Brian D. Stemper.

J Biomech Eng 135(4), 041008 (Apr 05, 2013) (14 pages) Paper No: BIO-12-1369; doi: 10.1115/1.4023982 History: Received August 22, 2012; Revised March 01, 2013

Detailed knowledge of knee joint kinematics and dynamic loading is essential for improving the design and outcomes of surgical procedures, tissue engineering applications, prosthetics design, and rehabilitation. The need for dynamic computational models that link kinematics, muscle and ligament forces, and joint contacts has long been recognized but such body-level forward dynamic models do not exist in recent literature. A main barrier in using computational models in the clinic is the validation of the in vivo contact, muscle, and ligament loads. The purpose of this study was to develop a full body, muscle driven dynamic model with subject specific leg geometries and validate it during squat and toe-rise motions. The model predicted loads were compared to in vivo measurements acquired with an instrumented knee implant. Data for this study were provided by the “Grand Challenge Competition to Predict In-Vivo Knee Loads” for the 2012 American Society of Mechanical Engineers Summer Bioengineering Conference. Data included implant and bone geometries, ground reaction forces, EMG, and the instrumented knee implant measurements. The subject specific model was developed in the multibody framework. The knee model included three ligament bundles for the lateral collateral ligament (LCL) and the medial collateral ligament (MCL), and one bundle for the posterior cruciate ligament (PCL). The implanted tibia tray was segmented into 326 hexahedral elements and deformable contacts were defined between the elements and the femoral component. The model also included 45 muscles on each leg. Muscle forces were computed for the muscle driven simulation by a feedback controller that used the error between the current muscle length in the forward simulation and the muscle length recorded during a kinematics driven inverse simulation. The predicted tibia forces and torques, ground reaction forces, electromyography (EMG) patterns, and kinematics were compared to the experimentally measured values to validate the model. Comparisons were done graphically and by calculating the mean average deviation (MAD) and root mean squared deviation (RMSD) for all outcomes. The MAD value for the tibia vertical force was 279 N for the squat motion and 325 N for the toe-rise motion, 45 N and 53 N for left and right foot ground reaction forces during the squat and 94 N and 82 N for toe-rise motion. The maximum MAD value for any of the kinematic outcomes was 7.5 deg for knee flexion-extension during the toe-rise motion.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Fregly, B. J., Besier, T. F., Lloyd, D. G., Delp, S. L., Banks, S. A., Pandy, M. G., and D'Lima, D. D., 2012, “Grand challenge competition to predict in vivo knee loads,” J. Orthop. Res., 30(4), pp. 503–513. [CrossRef] [PubMed]
Mundermann, A., Dyrby, C. O., D'Lima, D. D., Colwell, C. W., Jr., and Andriacchi, T. P., 2008, “In vivo knee loading characteristics during activities of daily living as measured by an instrumented total knee replacement,” J. Orthop. Res., 26(9), pp. 1167–1172. [CrossRef] [PubMed]
Heinlein, B., Kutzner, I., Graichen, F., Bender, A., Rohlmann, A., Halder, A. M., Beier, A., and Bergmann, G., 2009, “ESB Clinical Biomechanics Award 2008: Complete data of total knee replacement loading for level walking and stair climbing measured in vivo with a follow-up of 6-10 months,” Clin. Biomech. (Bristol, Avon), 24(4), pp. 315–326. [CrossRef] [PubMed]
Andriacchi, T. P., Mundermann, A., Smith, R. L., Alexander, E. J., Dyrby, C. O., and Koo, S., 2004, “A framework for the in vivo pathomechanics of osteoarthritis at the knee,” Ann. Biomed. Eng., 32(3), pp. 447–457. [CrossRef] [PubMed]
Wimmer, M. A., and Andriacchi, T. P., 1997, “Tractive forces during rolling motion of the knee: implications for wear in total knee replacement,” J. Biomech., 30(2), pp. 131–137. [CrossRef] [PubMed]
Sathasivam, S., and Walker, P. S., 1998, “Computer model to predict subsurface damage in tibial inserts of total knees,” J. Orthop. Res., 16(5), pp. 564–571. [CrossRef] [PubMed]
Blankevoort, L., Kuiper, J. H., Huiskes, R., and Grootenboer, H. J., 1991, “Articular contact in a three-dimensional model of the knee,” J. Biomech., 24(11), pp. 1019–1031. [CrossRef] [PubMed]
Pandy, M. G., Sasaki, K., and Kim, S., 1998, “A Three-Dimensional Musculoskeletal Model of the Human Knee Joint. Part 1: Theoretical Construct,” Comput. Methods Biomech. Biomed. Engin., 1(2), pp. 87–108. [CrossRef] [PubMed]
Abdel-Rahman, E. M., and Hefzy, M. S., 1998, “Three-dimensional dynamic behaviour of the human knee joint under impact loading,” Med. Eng. Phys., 20(4), pp. 276–290. [CrossRef] [PubMed]
Cohen, Z. A., Roglic, H., Grelsamer, R. P., Henry, J. H., Levine, W. N., Mow, V. C., and Ateshian, G. A., 2001, “Patellofemoral stresses during open and closed kinetic chain exercises. An analysis using computer simulation,” Am. J. Sports Med., 29(4), pp. 480–487. [PubMed]
Dhaher, Y. Y., and Kahn, L. E., 2002, “The effect of vastus medialis forces on patello-femoral contact: a model-based study,” J. Biomech. Eng., 124(6), pp. 758–767. [CrossRef] [PubMed]
Chao, E. Y., 2003, “Graphic-based musculoskeletal model for biomechanical analyses and animation,” Med. Eng. Phys., 25(3), pp. 201–212. [CrossRef] [PubMed]
Elias, J. J., Wilson, D. R., Adamson, R., and Cosgarea, A. J., 2004, “Evaluation of a computational model used to predict the patellofemoral contact pressure distribution,” J. Biomech., 37(3), pp. 295–302. [CrossRef] [PubMed]
Piazza, S. J., and Delp, S. L., 2001, “Three-dimensional dynamic simulation of total knee replacement motion during a step-up task,” J. Biomech. Eng., 123(6), pp. 599–606. [CrossRef] [PubMed]
Bei, Y., and Fregly, B. J., 2004, “Multibody dynamic simulation of knee contact mechanics,” Med. Eng. Phys., 26(9), pp. 777–789. [CrossRef] [PubMed]
Halloran, J. P., Petrella, A. J., and Rullkoetter, P. J., 2005, “Explicit finite element modeling of total knee replacement mechanics,” J. Biomech., 38(2), pp. 323–331. [CrossRef] [PubMed]
Kim, H. J., Fernandez, J. W., Akbarshahi, M., Walter, J. P., Fregly, B. J., and Pandy, M. G., 2009, “Evaluation of predicted knee-joint muscle forces during gait using an instrumented knee implant,” J. Orthop. Res., 27(10), pp. 1326–1331. [CrossRef] [PubMed]
Lundberg, H. J., Foucher, K. C., Andriacchi, T. P., and Wimmer, M. A., 2012, “Direct comparison of measured and calculated total knee replacement force envelopes during walking in the presence of normal and abnormal gait patterns,” J. Biomech., 45(6), pp. 990–996. [CrossRef] [PubMed]
Lundberg, H. J., Foucher, K. C., and Wimmer, M. A., 2009, “A parametric approach to numerical modeling of TKR contact forces,” J. Biomech., 42(4), pp. 541–545. [CrossRef] [PubMed]
Godest, A. C., Beaugonin, M., Haug, E., Taylor, M., and Gregson, P. J., 2002, “Simulation of a knee joint replacement during a gait cycle using explicit finite element analysis,” J. Biomech., 35(2), pp. 267–275. [CrossRef] [PubMed]
Baldwin, M. A., Clary, C. W., Fitzpatrick, C. K., Deacy, J. S., Maletsky, L. P., and Rullkoetter, P. J., 2012, “Dynamic finite element knee simulation for evaluation of knee replacement mechanics,” J. Biomech., 45(3), pp. 474–483. [CrossRef] [PubMed]
Lanovaz, J. L., and Ellis, R. E., 2009, “A cadaverically evaluated dynamic FEM model of closed-chain TKR mechanics,” J. Biomech. Eng., 131(5), p. 051002. [CrossRef]
Zelle, J., Heesterbeek, P. J., DeWaal Malefijt, M., and Verdonschot, N., 2010, “Numerical analysis of variations in posterior cruciate ligament properties and balancing techniques on total knee arthroplasty loading,” Med. Eng. Phys., 32(7), pp. 700–707. [CrossRef] [PubMed]
Fregly, B. J., Bei, Y., and Sylvester, M. E., 2003, “Experimental evaluation of an elastic foundation model to predict contact pressures in knee replacements,” J. Biomech., 36(11), pp. 1659–1668. [CrossRef] [PubMed]
Lifemodeler, I., 2010, “Lifemod Manual.”, Lifemodeler Inc., San Clemente, CA.
Wismans, J., Veldpaus, F., Janssen, J., Huson, A., and Struben, P., 1980, “A three-dimensional mathematical model of the knee-joint,” J. Biomech., 13(8), pp. 677–685. [CrossRef] [PubMed]
DeFrate, L. E., Gill, T. J., and Li, G., 2004, “In vivo function of the posterior cruciate ligament during weightbearing knee flexion,” Am. J. Sports Med., 32(8), pp. 1923–1928. [CrossRef] [PubMed]
Guess, T. M., Liu, H., Bhashyam, S., and Thiagarajan, G., 2013, “A multibody knee model with discrete cartilage prediction of tibio-femoral contact mechanics,” Comput. Meth. Biomech. Biomed. Eng., 16(3), pp. 256–270. [CrossRef]
Guess, T. M., Thiagarajan, G., Kia, M., and Mishra, M., 2010, “A subject specific multibody model of the knee with menisci,” Med. Eng. Phys., 32(5), pp. 505–515. [CrossRef] [PubMed]
Guess, T. M., and Stylianou, A., 2012, “Simulation of anterior cruciate ligament deficiency in a musculoskeletal model with anatomical knees,” Open Biomed. Eng. J., 6, pp. 23–32. [CrossRef] [PubMed]
Machado, M., Moreir, P., Flores, P., and Lankarani, H. M., 2012, “Compliant contact force models in multibody dynamics:Evolution of the Hertz contact theory.,” Mechanism and Machine Theory, 53, pp. 99–121. [CrossRef]
Guess, T. M., and Maletsky, L. P., 2005, “Computational modelling of a total knee prosthetic loaded in a dynamic knee simulator,” Med. Eng. Phys., 27(5), pp. 357–367. [CrossRef] [PubMed]
Li, L., Tong, K., Song, R., and Koo, T. K., 2007, “Is maximum isometric muscle stress the same among prime elbow flexors?,” Clin. Biomech. (Bristol, Avon), 22(8), pp. 874–883. [CrossRef] [PubMed]
Fukunaga, T., Roy, R. R., Shellock, F. G., Hodgson, J. A., and Edgerton, V. R., 1996, “Specific tension of human plantar flexors and dorsiflexors,” J. Appl. Physiol., 80(1), pp. 158–165. [PubMed]
McMahon, T. A., 1984, Muscles, reflexes, and locomotion, Princeton University Press, Princeton, NJ.
Smith, G., 1989, “Padding point extrapolation techniques for the Butterworth digital filter,” J. Biomech., 22(8–9), pp. 967–971. [CrossRef] [PubMed]
Winter, D. A., 1990, Biomechanics and motor control of human movement, Wiley, New York, NY.
Shelburne, K. B., Torry, M. R., and Pandy, M. G., 2006, “Contributions of muscles, ligaments, and the ground-reaction force to tibiofemoral joint loading during normal gait,” J. Orthop. Res., 24(10), pp. 1983–1990. [CrossRef] [PubMed]
Lin, Y. C., Walter, J. P., Banks, S. A., Pandy, M. G., and Fregly, B. J., 2010, “Simultaneous prediction of muscle and contact forces in the knee during gait,” J. Biomech., 43(5), pp. 945–952. [CrossRef] [PubMed]
Bloemker, K. H., Guess, T. M., Maletsky, L., and Dodd, K., 2012, “Computational knee ligament modeling using experimentally determined zero-load lengths,” Open Biomed. Eng. J., 6, pp. 33–41. [CrossRef] [PubMed]
Guess, T. M., 2012, “Forward Dynamics Simulation Using a Natural Knee with Menisci in the Multibody Framework,” Multibody Syst. Dyn., 28, pp. 37–53. [CrossRef]
Lloyd, D. G., and Buchanan, T. S., 2001, “Strategies of muscular support of varus and valgus isometric loads at the human knee,” J. Biomech., 34(10), pp. 1257–1267. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Full body multibody model with artificial left knee and discretized tibia insert

Grahic Jump Location
Fig. 2

Model predicted and measured tibial component forces and torques during squat motion

Grahic Jump Location
Fig. 3

Model predicted and measured ground reaction forces for left and right feet during squat motion

Grahic Jump Location
Fig. 4

Computed joint kinematics during muscle driven simulation (FD) versus inverse dynamics simulation (ID) during squat motion

Grahic Jump Location
Fig. 5

Model predicted muscle activations versus measured normalized EMG for the primary muscles involved in squat motion

Grahic Jump Location
Fig. 6

Contact pressure (N/mm2) distribution on tibia component during one cycle of squat motion

Grahic Jump Location
Fig. 7

Model predicted and measured tibial component forces and torques during toe-rise motion

Grahic Jump Location
Fig. 8

Model predicted and measured ground reaction forces for left and right feet during toe-rise motion

Grahic Jump Location
Fig. 9

Computed joint kinematics during muscle driven simulation (FD) versus inverse dynamics simulation (ID) during toe-rise motion

Grahic Jump Location
Fig. 10

Model predicted muscle activations versus measured normalized EMG for the primary muscles involved in toe-rise motion

Grahic Jump Location
Fig. 11

Contact pressure (N/mm2) distribution on tibia component during one cycle of toe-rise motion



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In