It is estimated that 27 million adults in the United States have clinical osteoarthritis with knee osteoarthritis being the most prevalent form affecting 28% of adults over age 45 and 37% of adults over age 60 [1RC Lawrence, DT Felson, CG Helmick, LM Arnold, H Choi, RA Deyo, S Gabriel, R Hirsch, MC Hochberg, GG Hunder, J M Jordan, JN Katz, HM Kremers, and F Wolfe, ""Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II"", Arthritis. Rheum, vol. 58, pp. 26-35, Jan 2008.]. Although the prevalence and debilitating nature of osteoarthritis is well documented, the etiology of this chronic disease is not completely understood. The most significant cause of osteoarthritis in young adults is joint trauma including tears of the menisci or ligaments. For example, 12 years after rupture of the anterior cruciate ligament (ACL), 51% of a female population with a mean age of 31 years had radiographic evidence of osteoarthritis [2LS Lohmander, A Ostenberg, M Englund, and H Roos, ""High prevalence of knee osteoarthritis, pain, and functional limitations in female soccer players twelve years after anterior cruciate ligament injury"", Arthritis Rheum, vol. 50, pp. 3145-52, Oct 2004.]. Similarly, 43% of patients with intact ACLs and isolated limited meniscectomy had radiographic evidence of osteoarthritis after 16 years [3M Englund, EM Roos, and LS Lohmander, ""Impact of type of meniscal tear on radiographic and symptomatic knee osteoarthritis: a sixteen-year followup of meniscectomy with matched controls"", Arthritis Rheum , vol. 48, pp. 2178-87, Aug 2003.]. In addition, the severity of osteoarthritis correlates with the severity of meniscal injury and the amount of tissue removed [4M Englund, and LS Lohmander, ""Risk factors for symptomatic knee osteoarthritis fifteen to twenty-two years after meniscectomy"", Arthritis. Rheum, vol. 50, pp. 2811-9, Sep 2004.].
Abnormal knee kinematics and meniscus injury resulting from ACL deficiency are often implicated in joint degeneration [5LS Lohmander, PM Englund, LL Dahl, and EM Roos, ""The long-term consequence of anterior cruciate ligament and meniscus injuries: osteoarthritis"", Am. J. Sports. Med, vol. 35, pp. 1756-69, Oct 2007.] while ACL reconstruction surgeries attempt to restore normal tibio-femoral motion and prevent long-term consequences [6C van Eck, Z Working, and F Fu, ""Current concepts in anatomic single- and double-bundle anterior cruciate ligament reconstruction"", Phys. Sportsmed, vol. 39, pp. 140-8, May 2011.]. The observed differences in tibio-femoral kinematics between normal and ACL deficient knees are small, typically only a few millimeters [7LE Defrate, R Papannagari, TJ Gill, JM Moses, NP Pathare, and G and Li, ""The 6 degrees of freedom kinematics of the knee after anterior cruciate ligament deficiency: an in vivo imaging analysis"", Am. J. Sports. Med, vol. 34, pp. 1240-6, Aug 2006., 8AD Georgoulis, A Papadonikolakis, CD Papageorgiou, A Mitsou, and N Stergiou, ""Three-dimensional tibiofemoral kinematics of the anterior cruciate ligament-deficient and reconstructed knee during walking"", Am. J. Sports. Med, vol. 31, pp. 75-9, Jan-Feb 2003.]. Studies have also shown that ACL reconstruction does not significantly reduce the incidence of early onset osteoarthritis [2LS Lohmander, A Ostenberg, M Englund, and H Roos, ""High prevalence of knee osteoarthritis, pain, and functional limitations in female soccer players twelve years after anterior cruciate ligament injury"", Arthritis Rheum, vol. 50, pp. 3145-52, Oct 2004.]. Increased knowledge of tibiofemoral joint contact mechanics would increase our understanding of the effects of ACL injury and help guide ACL reconstruction methods. Medical imaging has been used to quantify motion and cartilage contact in normal and ACL deficient knees including knee loading in a magnetic resonance imaging (MRI) tunnel [9JM Scarvell, PN Smith, KM Refshauge, H Galloway, and K Woods, ""Comparison of kinematics in the healthy and ACL injured knee using MRI"", J. Biomech, vol. 38, pp. 255-62, Feb 2005.] and MRI combined with fluoroscopic images captured during a weight-bearing lunge [10SK Van de Velde, JT Bingham, A Hosseini, M Kozanek, LE DeFrate, TJ Gill, and G Li, ""Increased tibiofemoral cartilage contact deformation in patients with anterior cruciate ligament deficiency"", Arthritis. Rheum, vol. 60, pp. 3693-702, Dec 2009.]. Computational models capable of predicting tibio-menisco-femoral contact mechanics during movement would also provide a valuable tool for increasing our understanding of normal joint mechanics as well as knee injury and joint degeneration.
The main objective of this study was to develop a muscle driven forward dynamics model of the lower extremities in a dual limb squat, validate the computational model against experimental data and then simulate the same motion in two different ACL deficient conditions. The model combined a cadaver anatomical knee with anthropometric and motion data from a female subject of similar height and weight as the cadaver donor. The model was developed in the multibody framework and it was validated against measurements of the ground reaction forces and surface electromyography (EMG) taken during the squat motion. The multibody model included representation of the menisci with an intact ACL, transected posterolateral ACL bundle, and full ACL transection. The contact pressures on the tibia plateau for the intact and ACL deficient cases were compared during a muscle driven forward dynamics simulation.
A previously developed anatomical computational knee model was used for this study. The multibody model was created in MD Adams (MSC Software Corporation, Santa Ana, CA) and has been previously described [11TM Guess, H Liu, S Bhashyam, and G Thiagarajan, ""A multibody knee model with discrete cartilage prediction of tibiofemoral contact mechanics"", Computer Methods in Biomechanics and Biomedical. Engineering, pp. 1-15, 2011.]. Only a brief summary is provided here. Knee geometries (tibia, femur, articular cartilage and menisci) were derived from magnetic resonance images (MRI) of a fresh frozen cadaver knee (55 years old, female, left knee, 170 cm height, 72 kg mass). After MRI, the cadaver knee was mounted in a dynamic knee simulator (Kansas Knee Simulator, University of Kansas, Lawrence, KS) and manipulated through a walk cycle. During testing, the forces applied by the actuators of the machine were recorded and bone motion was measured using rigid body markers attached to the femur, tibia, and patella with an Optotrak 3020 system (Northern Digital Inc., Waterloo, Ontario). For calculation of ligament zero-load lengths, the motion of the tibia and femur rigid body markers was recorded while the joint was manually moved through its full range of motion with a minimum force applied (as judged by the experimenter). After testing, the knee was disarticulated and the insertion sites for the anterior cruciate ligament (ACL), posterior cruciate ligament (PCL), medial collateral ligament (MCL), lateral collateral ligament (LCL), patellar tendon, menisci horn attachments, and menisci transverse ligament were measured with an Optotrak digitizing probe.
where k is a stiffness parameter, εl is a spring parameter assumed to be 0.03, ε is the engineering strain of each ligament bundle, l is ligament bundle length, and l0 is the zero-load length. Values of k for each ligament bundle came from Wismans et al. and Blankevoort et al. [15L Blankevoort, JH Kuiper, R Huiskes, and HJ Grootenboer, "" Articular contact in a three-dimensional model of the knee"", J. Biomech, vol. 24, pp. 1019-31, 1991., 16J Wismans, F Veldpaus, J Janssen, A Huson, and P Struben, ""A three-dimensional mathematical model of the knee-joint"", J. Biomech, vol. 13, pp. 677-85, 1980.]. Each 1-D spring also included a parallel damper and a damping coefficient of 0.5 Ns/mm was used for each ligament bundle.
where Fc is the contact force, δ is the interpenetration of geometries, δ. is the velocity of interpenetration, k is a spring constant, n is the compliance exponent, and B(δ) is a damping coefficient. The contact parameters defined between each cartilage rigid body and the femur cartilage were derived by optimizing the contact parameters to match contact pressure predictions of a finite element model [11TM Guess, H Liu, S Bhashyam, and G Thiagarajan, ""A multibody knee model with discrete cartilage prediction of tibiofemoral contact mechanics"", Computer Methods in Biomechanics and Biomedical. Engineering, pp. 1-15, 2011.]. Specifically, the contact parameters were optimized to minimize the difference between pressure predictions of the multibody model and that of an identically loaded and constrained finite element model. The resulting parameters were k = 140 N/mm1/1.3, n = 1.3, and B = 5 Ns/mm.
where Fθ,r,z and Tθ,r,z are the translational and torsional forces between elements acting in the circumferential, radial, and axial directions. Kθ, Kr, Kz, Kθr, Kθz, Krz, Tθ, Tr, and Tz are the stiffness matrix parameters, θ, r, and z are relative translational displacements and a, b, and c are relative rotational displacements. The stiffness matrix parameters are shown in Table 1. Values for the stiffness parameters were derived through an optimization process that minimized the displacement error between identically loaded finite element and multibody menisci models [18TM Guess, G Thiagarajan, M Kia, and M Mishra, ""A subject specific multibody model of the knee with menisci"", Med. Eng. Phys, vol. 32, pp. 505-15, Jun 2010.].
Each meniscus was attached to the tibia via 4 horn attachments (2 posterior and 2 medial per meniscus) modeled as simple tension only springs with a spring constant of 1000 N/mm [19TL Donahue, ML Hull, MM Rashid, and CR Jacobs, ""A finite element model of the human knee joint for the study of tibiofemoral contact"", J. Biomech. Eng, vol. 124, pp. 273-80, Jun 2002.]. The transverse ligament was modeled as a tension only spring with a stiffness of 200 N/mm. A parallel damper with damping coefficient of 0.5 Ns/mm was included for each horn attachment and the transverse ligament spring. Deformable contacts using Eq. 3 were defined between each meniscus element and the femur cartilage. In addition, a deformable contact (k = 80 N/mm1/2.5, n = 2.5) was defined between each meniscus element and every tibia cartilage element that might come in contact with the meniscus element during movement (122 contacts medial side and 134 contacts lateral side). A deformable contact was also defined between the patella cartilage and femur cartilage geometries (Fig. 1 ). To replicate the anatomical knee force-displacement relationships for the right leg, a right knee was created by mirroring the geometries and ligament insertion sites of the left knee. As the focus of the study was on the cadaver based left knee, the menisci were not included in the mirrored right knee.
Components of the multi-body knee model. The arrows indicate geometries connected through deformable contacts.
Left knee during the muscle-driven forward dynamics simulation. The small red arrows represent the forces from the contacts and ligaments defined in the model.
Measured ground reaction forces for the right and left feet (a.) and predicted ground reaction forces for the intact knee simulation (b.). The x-axis is oriented in the medial-lateral direction, the y-axis is in the anterior-posterior direction (positive is posterior and negative is anterior), and z is vertical.
Predicted muscle tension and measured EMG for the dual limb squat.
Flexion angle of the left hip during the dual limb squat.
Tibia plateau for the three model versions at extension (a.) and flexion (b.) during the forward dynamics simulation. The yellow circles show the location of the peak tibia cartilage contact pressure on the medial and lateral plateaus.
Contact pressure on the medial tibia plateau at extension and flexion.
Contact pressure on the lateral tibia plateau at extension and flexion.
Motion of the lateral meniscus relative to the tibia plateau for the three model versions. Shown is the displacement of the center of mass of the center meniscus rigid body (indicated with a white dot in the figure) relative to the center of the lateral tibia cartilage plateau. The posterior and lateral directions are positive. As a reference, hip flexion angle for the intact model is also shown.
Tibia plateau for the three model versions at 0.25 second intervals during the forward dynamics simulation. The shaded areas result from the interpenetration of geometries and are the contact patch.
Internal forces in the menisci rigid bodies for the three model versions at 0.25 second intervals during the forward dynamics simulation.
With informed consent, motion, ground reaction forces, and surface electromyography (EMG) were collected on a young female subject (age 27) of similar height (175 cm) and weight (64 kg) as that of the cadaver knee donor. Markers were placed on the subject according to the plug-in-gait protocol and a dual-limb squat activity was performed where the subject’s right and left feet were isolated on different force plates. During the squat the heels were always in contact with the ground and the subject was instructed to start with the knee slightly flexed. A squat activity with an initial flexion angle matching the flexion of the computational model was chosen (Fig. 2 ). The squat began with a knee flexion angle of 16° (0.25 s, extension) and the deepest flexion angle was 87° (1.625 s, flexion).
LifeMOD (LifeModeler, Inc., San Clemente, CA) was used for development and simulation of the musculoskeletal model (Fig. 2 ). LifeMOD is a virtual human modeling and simulation software add-on to MD ADAMS. The program is used for development of musculoskeletal structures and dynamic movement simulations based on subject anthropometrics and captured motion data. For forward dynamics simulations, LifeMOD predicts the muscle forces that will reproduce measured body motion. The model developed for this project included the pelvis, right lower limb, and left lower limb. The right and left hip were modeled as spherical joints as well as the right and left ankle. Anthropometric measurements from the test subject were used to scale the model. The previously developed left knee and mirrored right knee were placed in the musculoskeletal model with the femur and tibia of the anatomical knees rigidly attached to the musculoskeletal model femur and tibia. Measured motion data from the gait testing was imported into the model and knee flexion-extension motion was used to position the anatomical knees relative to the musculoskeletal model. This was accomplished by systematically repositioning the anatomical knees until the force required to move the leg though the prescribed motion was minimized. The forty-five muscle LifeMOD leg muscle set was then added to both the right and left legs based on relative positions of the hip, knee and ankle. The default attachments of the left quadriceps muscles were modified to insert on the patella based on experimental cadaver measurements (Fig. 3 ).
The measured motion from the dual-limb squat trial was then used in an inverse dynamics simulation to move the lower body as constrained by the hips, ankles, and anatomical knees. During this step the length of each muscle, through their respective via points, was recorded. Next, the motion constraints were removed and each foot was connected to ground with a 6-axis spring-damper. During forward dynamics simulations, muscle forces were generated (within physiological constraints) via feedback control loops that reproduced the muscle lengths recorded during the inverse dynamics training.
Muscle driven forward dynamics squat simulations were run on three versions of the model (Table 2). Joint angles, contact forces, muscles forces, meniscus displacement, and ground-reaction forces were recorded during simulation of the three models. Contact pressures on the left knee tibia plateaus were estimated by dividing the contact force on each element by its contact surface area (typically ~ 17 mm2).
Fig. (4 ) provides the measured and predicted ground reaction forces during the forward simulation for both the right and left foot for the intact model. The predicted ground reaction forces did not significantly change for the partial and full model versions. Fig. (5 ) shows the predicted muscle tension and measured EMG (raw signal rectified and low pass filtered) for four muscles on the left leg where EMG was collected (Peroneus Longus, Gluteus Medius, Vastus Lateralis, Medial Gastrocnemius). Between model versions the magnitude of muscle tension changes, but the overall pattern is generally consistent.
The tibio-femoral compressive force at flexion was 830 N on the lateral side and 2020 N on the medial side at flexion. The total compressive force at the joint was 2850 N or 4.5 times body weight of the subject used in the gait measurements.
During the squat cycle, all three model versions maintained the same hip flexion angle (Fig. 6 ). Hip flexion angle varied from 8° at the start of the squat (0.25 s) to 50° at the deepest portion of the squat (1.625 s) and the knee flexion angle varied from approximately 16° to 87° during the same time frame. Table 3 provides the peak predicted contact pressures on the medial and lateral plateaus at extension and flexion for all three model versions. At flexion (87°), peak contact pressure increased on the medial plateau with partial and full ACL transection. On the lateral side, the peak cartilage contact pressure occurred underneath the meniscus at flexion for both the partially and fully transected ACL. At extension, the peak contact force occurred under the lateral meniscus for the full ACL transection.
On the medial plateau, the location of the peak contact pressure and the contact location moved posterior and lateral with partial and full ACL transection at extension (Fig. 7 and Fig. 8 ). At flexion, the location of the peak contact force did not move, but the magnitude of the peak contact force increased. On the lateral plateau, the peak contact location moved to underneath the meniscus at extension for the full transection (Fig. 7 and Fig. 9 ). At flexion, the location of peak contact pressure moved from uncovered cartilage to underneath the meniscus for both the partial and full transection cases (Fig. 7 and Fig. 9 ).
Fig. (10 ) provides the motion of the center lateral meniscus rigid body relative to the tibia plateau. The meniscus moves in the posterior direction during extension prior to beginning the squat for all three models. The posterior motion was 3.5 mm greater for both the partial and full transection cases than for the intact knee. The lateral meniscus was pulled in the posterior direction by movement of the femur. The posterior motion occurred faster for the full transection compared to the partial transection case. For all three models, the menisci start in the same position at time zero. There was also a 1 mm difference in medial motion for the transected cases compared to the intact knee. The medial motion was created by tension in the posterior horn attachments and deformation of the entire meniscus as it was stretched in the anterior-posterior direction. Fig. (11 ) shows the tibial plateau for all three models at 0.25 second intervals during the squat cycle. At 0.50 seconds the greater posterior motion of the lateral tibia can be seen for the transected cases. In Fig. (12 ) the internal force experienced by the menisci rigid bodies is shown at 0.25 second intervals during the motion. After 0.75 seconds the internal forces experienced by the lateral meniscus is altered for the partial and full transection cases.
During flexion, the lateral meniscus becomes “wedged” between the tibia and femur for both partial and full ACL transection as indicated by the cartilage contact pressure on elements under the menisci and the changes in the forces experienced by the menisci rigid bodies. This “wedging” effect does not occur for the intact knee. The location of the highest “wedging” force corresponds to the periphery of the central-anterior region of the lateral meniscus. The wedging also occurred at extension for the full transection case at the interior of the central-posterior region of the lateral meniscus (Fig. 7 and Fig. 9 ). The internal forces in the lateral meniscus are also altered for the partial and full transection cases. In particular, the circumferential forces decrease in the posterior lateral meniscus as the knee flexes (Fig. 12 ) and the internal shearing forces increase in the anterior lateral meniscus, where the wedging occurs.
A study evaluating the location of meniscus tears on 476 patients with ACL injury showed that tears were equally distributed between the lateral and medial side and more tears occurred on the posterior of the lateral meniscus than the anterior [20JP3rd Smith, and GR Barrett, ""Medial and lateral meniscal tear patterns in anterior cruciate ligament- deficient knees. A prospective analysis of 575 tears"", Am. J. Sports. Med, vol. 29, pp. 415-9, 2001.]. The location of meniscus injury evaluated on 159 patients within 3 days of ACL injury during skiing revealed that 83 % of the tears occurred on the lateral side [21JB Duncan, R Hunter, M Purnell, and J Freeman, ""Meniscal injuries associated with acute anterior cruciate ligament tears in alpine skiers"", Am. J. Sports. Med, vol. 23, pp. 170-2, Mar-Apr 1995.]. In the current study, only a simple dual limb squat was simulated and contact pressures were greatest during flexion. More complicated activities such as gait and stair climbing would result in different tibio-femoral loading and motion possibly causing the meniscus “wedging” to occur at different locations.
The tibio-femoral compressive forces in this model were higher than the ones reported by Mundermann et al., [22A Mundermann, CO Dyrby, DD D'Lima, CW Colwell Jr, and TP Andriacchi, ""In vivo knee loading characteristics during activities of daily living as measured by an instrumented total knee replacement", J. Orthop. Res, vol. 26, pp. 1167-72, Sep 2008.] for a similar squat motion. In Mundermann’s study the compressive force was measured directly using an instrumented tibial insert after knee reconstruction surgery and it was 2.5 times body weight at flexion. The discrepancy between the natural knee model compressive force and the measured prosthetic compressive force may be explained by the fact that the musculoskeletal model was a combination of a generic skeletal model based on gait subject height and weight and an MRI derived anatomical knee. The stitching of the geometries could have resulted in slight misalignments of the muscle lines of actions that in turn caused higher loads on the cartilage contacts. This issue will be investigated further in future models where the geometries of all bones will be derived from one subject’s MRI data.
Peak contact pressure and contact locations are similar for the partial tear (posterolateral bundle transection) and complete ACL transection cases at flexion, particularly on the lateral side. For the squat simulations, the tibio-femoral contact location shifts to the posterior and lateral direction with ACL transection. This displacement is only a few millimeters and is larger at extension. The change in peak contact and contact patch location on the tibia plateau is consistent with experimental measurements and observations [9JM Scarvell, PN Smith, KM Refshauge, H Galloway, and K Woods, ""Comparison of kinematics in the healthy and ACL injured knee using MRI"", J. Biomech, vol. 38, pp. 255-62, Feb 2005., 10SK Van de Velde, JT Bingham, A Hosseini, M Kozanek, LE DeFrate, TJ Gill, and G Li, ""Increased tibiofemoral cartilage contact deformation in patients with anterior cruciate ligament deficiency"", Arthritis. Rheum, vol. 60, pp. 3693-702, Dec 2009., 23JM Scarvell, PN Smith, KM Refshauge, HR Galloway, and K.R Woods, ""Association between abnormal kinematics and degenerative change in knees of people with chronic anterior cruciate ligament deficiency: a magnetic resonance imaging study"", Aust. J. Physiother, vol. 51, pp. 233-40, 2005.]. For both the partial and full transection simulations, small changes in femur motion relative to the tibia push the lateral meniscus in the posterior direction, exposing it to “wedging”. This “wedging” effect increases the loads on the meniscus and it is possible that the increased loading damages the tissue and increases the risk of cartilage degeneration during ACL deficiency [24A Meunier, M Odensten, and L Good, ""Long-term results after primary repair or non-surgical treatment of anterior cruciate ligament rupture: a randomized study with a 15-year follow-up"", Scand. J. Med. Sci. Sports, vol. 17, pp. 230-7, Jun 2007., 25JK Seon, EK Song, and SJ Park, ""Osteoarthritis after anterior cruciate ligament reconstruction using a patellar tendon autograft"", Int. Orthop , vol. 30, pp. 94-8, Apr 2006.].
For all simulations only the lower body was modeled. As seen in Fig. (2 ), the upper body of the test subject leans forward during the squat, allowing the body’s center of pressure to stay between its base of support. During the forward dynamics simulations, a LifeMOD tracking agent was used to maintain balance during the squat. The tracking agent is a 6-axis spring located between a “dummy” rigid body and the pelvis. The “dummy” rigid body is driven by a motion constraint that follows pelvis motion measured during the inverse dynamics simulation. If the forward dynamics pelvis motion closely follows the inverse dynamics pelvis motion, the tracking agent will have minimal influence on the model. Also, the tracking agent imparts no force in the vertical direction, regardless of pelvis motion. Measured ground reaction forces are symmetrical in the medial-lateral and anterior-posterior directions, but the simulation results include asymmetry in both axes (Fig. 4 ). For example, both the right and left knee have posterior shear forces while the knee was flexing and an anterior ground reaction shear force when the knee was extending. The model pushes against the tracking agent when these asymmetries occur and these forces will have an influence on muscle force prediction. Misalignment of the knee relative to the hip and ankle and slight misalignment between the left and right sides may be responsible for the predicted ground reaction force asymmetries. Previous work that used a hinge joint for the right knee showed significant asymmetries in the vertical ground reaction forces during the squat [26T Guess, ""Forward dynamics simulation using a natural knee with menisci in the multibody framework"", Multibody. System. Dynamics, pp. 1-17, .]. Changing the right knee from a hinge joint to an anatomical joint eliminated the vertical asymmetry. In addition, the vertical ground reaction force predictions may be further improved by including representation of the upper body. In the current model simulations, vertical ground reaction forces were greater for the left foot than the right. Experimentally, the right foot had a greater vertical reaction force (Fig. 4 ). Ongoing work in our lab has shown that including representation of the upper body improves prediction of ground reaction forces, particularly the vertical force distribution between the right and left leg. Future musculoskeletal models will include representation of the upper body.
Combining a cadaver knee with in vivo gait measurements from a subject of similar height and weight allows in vitro validation of the knee model, but it is also a limitation of the study. Additionally, the same measured motion from a subject with an intact ACL was used for the inverse dynamics solution of all three model versions. ACL deficient subjects demonstrate altered gait patterns [27R Ferber, LR Osternig, MH Woollacott, NJ Wasielewski, and JH Lee, ""Gait mechanics in chronic ACL deficiency and subsequent repair"", Clin. Biomech. (Bristol., Avon), vol. 17, pp. 274-85, May 2002.] and it is possible that the kinematics of the dual limb squat would also be different. This altered motion could affect simulation predictions at the knee. Finally, only a limited number of muscle activations were measured with EMG and maximum voluntary muscle contraction was not measured.
The primary objective of this study was to create a musculoskeletal model of the lower extremities that included anatomical representation of the knees, validate the model with experimental data, and then simulate two ACL deficient conditions for a simple motion (squat). The current study is a first step in the development of anatomically correct subject specific musculoskeletal models of the lower extremities to study joint function. Future refinements of this modeling method will include finer segmentations of the menisci and cartilage structures and validation of the models in other dynamic motions.
The authors have no conflict of interested with the presented work.
Development of the cadaver knee model was funded by the National Science Foundation, Grant Number 506297, under the IMAG program for Multiscale Modeling. The authors gratefully acknowledge the work of Mohammad Kia and Gavin Paiva in the development of the macros used to generate the multibody tissue models for the menisci and cartilage. The authors also acknowledge the contributions of researchers in the Experimental Joint Biomechanics Research Lab at the University of Kansas, Lawrence KS.
 RC Lawrence, DT Felson, CG Helmick, LM Arnold, H Choi, RA Deyo, S Gabriel, R Hirsch, MC Hochberg, GG Hunder, J M Jordan, JN Katz, HM Kremers, and F Wolfe, ""Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II"", Arthritis. Rheum, vol. 58, pp. 26-35, Jan 2008.
 LS Lohmander, A Ostenberg, M Englund, and H Roos, ""High prevalence of knee osteoarthritis, pain, and functional limitations in female soccer players twelve years after anterior cruciate ligament injury"", Arthritis Rheum, vol. 50, pp. 3145-52, Oct 2004.
 M Englund, EM Roos, and LS Lohmander, ""Impact of type of meniscal tear on radiographic and symptomatic knee osteoarthritis: a sixteen-year followup of meniscectomy with matched controls"", Arthritis Rheum , vol. 48, pp. 2178-87, Aug 2003.
 M Englund, and LS Lohmander, ""Risk factors for symptomatic knee osteoarthritis fifteen to twenty-two years after meniscectomy"", Arthritis. Rheum, vol. 50, pp. 2811-9, Sep 2004.
 LS Lohmander, PM Englund, LL Dahl, and EM Roos, ""The long-term consequence of anterior cruciate ligament and meniscus injuries: osteoarthritis"", Am. J. Sports. Med, vol. 35, pp. 1756-69, Oct 2007.
 C van Eck, Z Working, and F Fu, ""Current concepts in anatomic single- and double-bundle anterior cruciate ligament reconstruction"", Phys. Sportsmed, vol. 39, pp. 140-8, May 2011.
 LE Defrate, R Papannagari, TJ Gill, JM Moses, NP Pathare, and G and Li, ""The 6 degrees of freedom kinematics of the knee after anterior cruciate ligament deficiency: an in vivo imaging analysis"", Am. J. Sports. Med, vol. 34, pp. 1240-6, Aug 2006.
 AD Georgoulis, A Papadonikolakis, CD Papageorgiou, A Mitsou, and N Stergiou, ""Three-dimensional tibiofemoral kinematics of the anterior cruciate ligament-deficient and reconstructed knee during walking"", Am. J. Sports. Med, vol. 31, pp. 75-9, Jan-Feb 2003.
 JM Scarvell, PN Smith, KM Refshauge, H Galloway, and K Woods, ""Comparison of kinematics in the healthy and ACL injured knee using MRI"", J. Biomech, vol. 38, pp. 255-62, Feb 2005.
 SK Van de Velde, JT Bingham, A Hosseini, M Kozanek, LE DeFrate, TJ Gill, and G Li, ""Increased tibiofemoral cartilage contact deformation in patients with anterior cruciate ligament deficiency"", Arthritis. Rheum, vol. 60, pp. 3693-702, Dec 2009.
 TM Guess, H Liu, S Bhashyam, and G Thiagarajan, ""A multibody knee model with discrete cartilage prediction of tibiofemoral contact mechanics"", Computer Methods in Biomechanics and Biomedical. Engineering, pp. 1-15, 2011.
 VB Duthon, C Barea, S Abrassart, JH Fasel, D Fritschy, and J Menetrey, ""Anatomy of the anterior cruciate ligament"", Knee. Surg Sports. Traumatol. Arthrosc, vol. 14, pp. 204-13, Mar 2006.
 AA Amis, AM Bull, CM Gupte, I Hijazi, A Race, and JR Robinson, ""Biomechanics of the PCL and related structures: posterolateral, posteromedial and meniscofemoral ligaments"", Knee Surg. Sports. Traumatol. Arthrosc, vol. 11, pp. 271-81, Sep 2003.
 SE Park, LE DeFrate, JF Suggs, TJ Gill, HE Rubash, and G Li, ""The change in length of the medial and lateral collateral ligaments during in vivo knee flexion"", Knee, vol. 12, pp. 377-82, Oct 2005.
 L Blankevoort, JH Kuiper, R Huiskes, and HJ Grootenboer, "" Articular contact in a three-dimensional model of the knee"", J. Biomech, vol. 24, pp. 1019-31, 1991.
 J Wismans, F Veldpaus, J Janssen, A Huson, and P Struben, ""A three-dimensional mathematical model of the knee-joint"", J. Biomech, vol. 13, pp. 677-85, 1980.
 DC Fithian, MA Kelly, and VC Mow, ""Material properties and structure-function relationships in the menisci"", Clin. Orthop, pp. 9-31, 1990.
 TM Guess, G Thiagarajan, M Kia, and M Mishra, ""A subject specific multibody model of the knee with menisci"", Med. Eng. Phys, vol. 32, pp. 505-15, Jun 2010.
 TL Donahue, ML Hull, MM Rashid, and CR Jacobs, ""A finite element model of the human knee joint for the study of tibiofemoral contact"", J. Biomech. Eng, vol. 124, pp. 273-80, Jun 2002.
 JP3rd Smith, and GR Barrett, ""Medial and lateral meniscal tear patterns in anterior cruciate ligament- deficient knees. A prospective analysis of 575 tears"", Am. J. Sports. Med, vol. 29, pp. 415-9, 2001.
 JB Duncan, R Hunter, M Purnell, and J Freeman, ""Meniscal injuries associated with acute anterior cruciate ligament tears in alpine skiers"", Am. J. Sports. Med, vol. 23, pp. 170-2, Mar-Apr 1995.
 A Mundermann, CO Dyrby, DD D'Lima, CW Colwell Jr, and TP Andriacchi, ""In vivo knee loading characteristics during activities of daily living as measured by an instrumented total knee replacement", J. Orthop. Res, vol. 26, pp. 1167-72, Sep 2008.
 JM Scarvell, PN Smith, KM Refshauge, HR Galloway, and K.R Woods, ""Association between abnormal kinematics and degenerative change in knees of people with chronic anterior cruciate ligament deficiency: a magnetic resonance imaging study"", Aust. J. Physiother, vol. 51, pp. 233-40, 2005.
 A Meunier, M Odensten, and L Good, ""Long-term results after primary repair or non-surgical treatment of anterior cruciate ligament rupture: a randomized study with a 15-year follow-up"", Scand. J. Med. Sci. Sports, vol. 17, pp. 230-7, Jun 2007.
 JK Seon, EK Song, and SJ Park, ""Osteoarthritis after anterior cruciate ligament reconstruction using a patellar tendon autograft"", Int. Orthop , vol. 30, pp. 94-8, Apr 2006.
 T Guess, ""Forward dynamics simulation using a natural knee with menisci in the multibody framework"", Multibody. System. Dynamics, pp. 1-17, .
 R Ferber, LR Osternig, MH Woollacott, NJ Wasielewski, and JH Lee, ""Gait mechanics in chronic ACL deficiency and subsequent repair"", Clin. Biomech. (Bristol., Avon), vol. 17, pp. 274-85, May 2002.