The study objective was to investigate the influence of coronal plane alignment and ligament properties on total knee replacement (TKR) contact loads during walking. We created a subject-specific knee model of an 83-year-old male who had an instrumented TKR. The knee model was incorporated into a lower extremity musculoskeletal model and included deformable contact, ligamentous structures, and six degrees-of-freedom (DOF) tibiofemoral and patellofemoral joints. A novel numerical optimization technique was used to simultaneously predict muscle forces, secondary knee kinematics, ligament forces, and joint contact pressures from standard gait analysis data collected on the subject. The nominal knee model predictions of medial, lateral, and total contact forces during gait agreed well with TKR measures, with root-mean-square (rms) errors of 0.23, 0.22, and 0.33 body weight (BW), respectively. Coronal plane component alignment did not affect total knee contact loads, but did alter the medial-lateral load distribution, with 4 deg varus and 4 deg valgus rotations in component alignment inducing þ17% and À23% changes in the first peak medial tibiofemoral contact forces, respectively. A Monte Carlo analysis showed that uncertainties in ligament stiffness and reference strains induce 60.2 BW uncertainty in tibiofemoral force estimates over the gait cycle. Ligament properties had substantial influence on the TKR load distributions, with the medial collateral ligament and iliotibial band (ITB) properties having the largest effects on medial and lateral compartment loading, respectively. The computational framework provides a viable approach for virtually designing TKR components, considering parametric uncertainty and predicting the effects of joint alignment and soft tissue balancing procedures on TKR function during movement.

Articular geometry in the knee varies widely among people which has implications for risk of injury and pathology. The goals of this work were to develop a framework to systematically vary geometry in a multibody knee model and to use this framework to investigate the effect of morphological features on dynamic knee kinematics and contact mechanics. A statistical shape model of the tibiofemoral and patellofemoral joints was created from magnetic resonance images of 14 asymptomatic knees. The shape model was then used to generate 37 unique multibody knee models based on −3 to +3 standard deviations of the scores for the first six principal components identified. Each multibody model was then incorporated into a lower extremity musculoskeletal model and the Concurrent Optimization of Muscle Activations and Kinematics (COMAK) routine was used to simulate knee mechanics for overground walking. Changes in articular geometry affected knee function, resulting in differences up to 17° in orientation, 8 mm in translation, 0.7 BW in contact force, and 2.0 MPa in mean cartilage contact pressure. Understanding the relationship between shape and function in a joint could provide insight into the mechanisms behind injury and pathology and the variability in response to treatment.

Interventions used to treat patellofemoral pain in runners are often designed to alter patellofemoral mechanics. This study used a computational model to investigate the influence of two interventions, step rate manipulation and quadriceps strengthening, on patellofemoral contact pressures during running. Running mechanics were analyzed using a lower extremity musculoskeletal model that included a knee with six degree-of-freedom tibiofemoral and patellofemoral joints. An elastic foundation model was used to compute articular contact pressures. The lower extremity model was scaled to anthropometric dimensions of 22 healthy adults, who ran on an instrumented treadmill at 90%, 100% and 110% of their preferred step rate. Numerical optimization was then used to predict the muscle forces, secondary tibiofemoral kinematics and all patellofemoral kinematics that would generate the measured hip, knee and ankle joint accelerations. Mean and peak patella contact pressures reached 5.0 and 9.7 MPa during the midstance phase of running. Increasing step rate by 10% significantly reduced mean contact pressures by 10.4% and contact area by 7.4%, but had small effects on lateral patella translation and tilt. Enhancing vastus medialis strength did not substantially affect pressure magnitudes or lateral patella translation, but did shift contact pressure medially toward the patellar median ridge. Thus, the model suggests that step rate tends to primarily modulate the magnitude of contact pressure and contact area, while vastus medialis strengthening has the potential to alter mediolateral pressure locations. These results are relevant to consider in the design of interventions used to prevent or treat patellofemoral pain in runners.

Computational knee models provide a powerful platform to investigate the effects of injury and surgery on functional knee behavior. The objective of this study was to use a multibody knee model to investigate the influence of ligament properties on tibiofemoral kinematics and cartilage contact pressures in the stance phase of walking. The knee model included 14 ligament bundles and articular cartilage contact acting across the tibiofemoral and patellofemoral joints. The knee was incorporated into a lower extremity musculoskeletal model and used to simulate knee mechanics during the stance phase of normal walking. A Monte Carlo approach was employed to assess the influence ligament stiffness and reference strains on knee mechanics. The ACL, MCL and posterior capsule properties exhibited significant influence on anterior tibial translation at heel strike, with the ACL acting as the primary restraint to anterior translation in mid-stance. The MCL and LCL exhibited the greatest influence on tibial rotation from heel strike through mid-stance. Simulated tibial plateau contact location was dependent on the ACL, MCL and LCL properties, while pressure magnitudes were most dependent on the ACL. A decrease in ACL stiffness or reference strain significantly increased average contact pressure in mid-stance, with the pressure migrating posteriorly on the medial tibial plateau. These ligament-dependent shifts in tibiofemoral cartilage contact during walking are potentially relevant to consider when investigating the causes of early onset osteoarthritis following knee ligament injury and surgical treatment.

Background: Graft placement is a modifiable and often discussed surgical factor in anterior cruciate ligament (ACL) reconstruction (ACLR). However, the sensitivity of functional knee mechanics to variability in graft placement is not well understood. Purpose: To (1) investigate the relationship of ACL graft tunnel location and graft angle with tibiofemoral kinematics in patients with ACLR, (2) compare experimentally measured relationships with those observed with a computational model to assess the predictive capabilities of the model, and (3) use the computational model to determine the effect of varying ACL graft tunnel placement on tibiofemoral joint mechanics during walking. Study Design: Controlled laboratory study. Methods: Eighteen participants who had undergone ACLR were tested. Bilateral ACL footprint location and graft angle were assessed using magnetic resonance imaging (MRI). Bilateral knee laxity was assessed at the completion of rehabilitation. Dynamic MRI was used to measure tibiofemoral kinematics and cartilage contact during active knee flexion-extension. Additionally, a total of 500 virtual ACLR models were created from a nominal computational knee model by varying ACL footprint locations, graft stiffness, and initial tension. Laxity tests, active knee extension, and walking were simulated with each virtual ACLR model. Linear regressions were performed between internal knee mechanics and ACL graft tunnel locations and angles for the patients with ACLR and the virtual ACLR models. Results: Static and dynamic MRI revealed that a more vertical graft in the sagittal plane was significantly related ( P < .05) to a greater laxity compliance index ( R2 = 0.40) and greater anterior tibial translation and internal tibial rotation during active knee extension ( R2 = 0.22 and 0.23, respectively). Similarly, knee extension simulations with the virtual ACLR models revealed that a more vertical graft led to greater laxity compliance index, anterior translation, and internal rotation ( R2 = 0.56, 0.26, and 0.13). These effects extended to simulations of walking, with a more vertical ACL graft inducing greater anterior tibial translation, ACL loading, and posterior migration of contact on the tibial plateaus. Conclusion: This study provides clinical evidence from patients who underwent ACLR and from complementary modeling that functional postoperative knee mechanics are sensitive to graft tunnel locations and graft angle. Of the factors studied, the sagittal angle of the ACL was particularly influential on knee mechanics. Clinical Relevance: Early-onset osteoarthritis from altered cartilage loading after ACLR is common. This study shows that postoperative cartilage loading is sensitive to graft angle. Therefore, variability in graft tunnel placement resulting in small deviations from the anatomic ACL angle might contribute to the elevated risk of osteoarthritis after ACLR.