In vivo patellofemoral forces in high flexion total knee arthroplasty

This study compares the in vivo patellofemoral contact forces generated in high flexion fixed bearing posterior cruciate retaining Nexgen CR-Flex (PCR) and high flexion posterior stabilized Nexgen LPS-Flex (LPS) TKAs with that of normal knees from full knee extension to maximum weight bearing flexion. Ten patients with the PCR total knee arthroplasty (TKA), ten with the LPS TKA and seven patients having normal knees were fluoroscoped while performing a deep knee bend activity. In vivo femorotibial kinematics, obtained from 3D-to-2D registration technique, and patellar kinematics obtained by direct measurements from the fluoroscopic images were entered into a 3D inverse dynamics mathematical model to determine the in vivo contact forces at the knee. The variation in the patellofemoral and quadriceps forces with flexion were found to be similar across the three groups-increasing from full extension to 90 degrees of flexion, reaching a maximum between 90 degrees and 120 degrees of flexion and then decreasing until maximum flexion. At maximum knee flexion, these forces were found to be significantly lower in the normal knees than in the TKAs. The patellar ligament to quadriceps force ratio decreased with the increase in knee flexion while the patellofemoral to quadriceps force ratio increased. A strong correlation was found to exist between the patellofemoral forces, the femorotibial contact forces and the forces in the extensor mechanism. The PCR TKA in this study exhibited greater resemblance to the normal patients with respect to the patellofemoral forces than the LPS TKA though significant differences in the two implant types were not observed.     

The importance of quadriceps and hamstring muscle loading on knee kinematics and in-situ forces in the ACL

This study investigated the effect of hamstring co-contraction with quadriceps on the kinematics of the human knee joint and the in-situ forces in the anterior cruciate ligament (ACL) during a simulated isometric extension motion of the knee. Cadaveric human knee specimens (n = 10) were tested using the robotic universal force moment sensor (UFS) system and measurements of knee kinematics and in-situ forces in the ACL were based on reference positions on the path of passive flexion/extension motion of the knee. With an isolated 200 N quadriceps load, the knee underwent anterior and lateral tibial translation as well as internal tibial rotation with respect to the femur. Both translation and rotation increased when the knee was flexed from full extension to 30 of flexion; with further flexion, these motion decreased. The addition of 80 N antagonistic hamstrings load significantly reduced both anterior and lateral tibial translation as well as internal tibial rotation at knee flexion angles tested except at full extension. At 30 of flexion, the anterior tibial translation, lateral tibial translation, and internal tibial rotation were significantly reduced by 18, 46, and 30%, respectively (p<0.05). The in-situ forces in the ACL under the quadriceps load were found to increase from 27.8+/-9.3 N at full extension to a maximum of 44.9+/-13.8 N at 15 of flexion and then decrease to 10 N beyond 60 of flexion. The in-situ force at 15 was significantly higher than that at other flexion angles (p<0.05). The addition of the hamstring load of 80 N significantly reduced the in-situ forces in the ACL at 15, 30 and 60 of flexion by 30, 43, and 44%, respectively (p<0.05). These data demonstrate that maximum knee motion may not necessarily correspond to the highest in-situ forces in the ACL. The data also suggest that hamstring co-contraction with quadriceps is effective in reducing excessive forces in the ACL particularly between 15 and 60 of knee flexion.     

Forces acting on the patella during maximal voluntary contraction of the quadriceps femoris muscle at different knee flexion/extension angles

From knee extension moments measured with a dynamometer, the quadriceps muscle force, the patellar ligament force and the reaction force in the patellofemoral joint at various knee angles (0-90 degrees) were estimated. The information needed to calculate the combined effect of both patellofemoral and tibiofemoral joint on the mechanical advantage of the muscle was obtained from lateral-view radiographs of autopsy knees. The results show that the smallest quadriceps force (2,000 N) is exerted at maximal extension, and the largest force (8,000 N) at about 75 degrees of flexion. The patellar ligament force reaches a maximum (5,000 N) at 60 degrees. The reaction force in the patellofemoral joint is the smallest (1,000 N) at extension and is of the same values as the muscle force in a range from 75 to 90 degrees. Especially at large flexion angles, the value of the estimated forces is considerably larger (by 100%) than reported in the literature. This difference is attributed to the influence of the patellofemoral joint on the mechanical advantage of the muscle, which has not been taken into account in other studies.     

Role of gastrocnemius activation in knee joint biomechanics: gastrocnemius acts as an ACL antagonist

Gastrocnemius is a premier muscle crossing the knee, but its role in knee biomechanics and on the anterior cruciate ligament (ACL) remains less clear when compared to hamstrings and quadriceps. The effect of changes in gastrocnemius force at late stance when it peaks on the knee joint response and ACL force was initially investigated using a lower extremity musculoskeletal model driven by gait kinematics—kinetics. The tibiofemoral joint under isolated isometric contraction of gastrocnemius was subsequently analyzed at different force levels and flexion angles (0°–90°). Changes in gastrocnemius force at late stance markedly influenced hamstrings forces. Gastrocnemius acted as ACL antagonist by substantially increasing its force. Simulations under isolated contraction of gastrocnemius confirmed this role at all flexion angles, in particular, at extreme knee flexion angles (0° and 90°). Constraint on varus/valgus rotations substantially decreased this effect. Although hamstrings and gastrocnemius are both knee joint flexors, they play opposite roles in respectively protecting or loading ACL. Although the quadriceps is also recognized as antagonist of ACL, at larger joint flexion and in contrast to quadriceps, activity in gastrocnemius substantially increased ACL forces (anteromedial bundle). The fact that gastrocnemius is an antagonist of ACL should help in effective prevention and management of ACL injuries.     

A detailed and validated three dimensional dynamic model of the patellofemoral joint

A detailed 3D anatomical model of the patellofemoral joint was developed to study the tracking, force, contact and stability characteristics of the joint. The quadriceps was considered to include six components represented by 15 force vectors. The patellar tendon was modeled using four bundles of viscoelastic tensile elements. Each of the lateral and medial retinaculum was modeled by a three-bundle nonlinear spring. The femur and patella were considered as rigid bodies with their articular cartilage layers represented by an isotropic viscoelastic material. The geometrical and tracking data needed for model simulation, as well as validation of its results, were obtained from an in vivo experiment, involving MR imaging of a normal knee while performing isometric leg press against a constant 140 N force. The model was formulated within the framework of a rigid body spring model and solved using forth-order Runge-Kutta, for knee flexion angles between zero and 50 degrees. Results indicated a good agreement between the model predictions for patellar tracking and the experimental results with RMS deviations of about 2 mm for translations (less than 0.7 mm for patellar mediolateral shift), and 4 degrees for rotations (less than 3 degrees for patellar tilt). The contact pattern predicted by the model was also consistent with the results of the experiment and the literature. The joint contact force increased linearly with progressive knee flexion from 80 N to 210 N. The medial retinaculum experienced a peak force of 18 N at full extension that decreased with knee flexion and disappeared entirely at 20 degrees flexion. Analysis of the patellar time response to the quadriceps contraction suggested that the muscle activation most affected the patellar shift and tilt. These results are consistent with the recent observations in the literature concerning the significance of retinaculum and quadriceps in the patellar stability.     

Knee muscle forces during walking and running in patellofemoral pain patients and pain-free controls

One proposed mechanism of patellofemoral pain, increased stress in the joint, is dependent on forces generated by the quadriceps muscles. Describing causal relationships between muscle forces, tissue stresses, and pain is difficult due to the inability to directly measure these variables in vivo. The purpose of this study was to estimate quadriceps forces during walking and running in a group of male and female patients with patellofemoral pain (n=27, 16 female; 11 male) and compare these to pain-free controls (n=16, 8 female; 8 male). Subjects walked and ran at self-selected speeds in a gait laboratory. Lower limb kinematics and electromyography (EMG) data were input to an EMG-driven musculoskeletal model of the knee, which was scaled and calibrated to each individual to estimate forces in 10 muscles surrounding the joint. Compared to controls, the patellofemoral pain group had greater co-contraction of quadriceps and hamstrings (p=0.025) and greater normalized muscle forces during walking, even though the net knee moment was similar between groups. Muscle forces during running were similar between groups, but the net knee extension moment was less in the patellofemoral pain group compared to controls. Females displayed 30–50% greater normalized hamstring and gastrocnemius muscle forces during both walking and running compared to males (p<0.05). These results suggest that some patellofemoral pain patients might experience greater joint contact forces and joint stresses than pain-free subjects. The muscle force data are available as supplementary material.     

Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces

Estimating tibiofemoral joint contact forces is important for understanding the initiation and progression of knee osteoarthritis. However, tibiofemoral contact force predictions are influenced by many factors including muscle forces and anatomical representations of the knee joint. This study aimed to investigate the influence of subject-specific geometry and knee joint kinematics on the prediction of tibiofemoral contact forces using a calibrated EMG-driven neuromusculoskeletal model of the knee. One participant fitted with an instrumented total knee replacement walked at a self-selected speed while medial and lateral tibiofemoral contact forces, ground reaction forces, whole-body kinematics, and lower-limb muscle activity were simultaneously measured. The combination of generic and subject-specific knee joint geometry and kinematics resulted in four different OpenSim models used to estimate muscle–tendon lengths and moment arms. The subject-specific geometric model was created from CT scans and the subject-specific knee joint kinematics representing the translation of the tibia relative to the femur was obtained from fluoroscopy. The EMG-driven model was calibrated using one walking trial, but with three different cost functions that tracked the knee flexion/extension moments with and without constraint over the estimated joint contact forces. The calibrated models then predicted the medial and lateral tibiofemoral contact forces for five other different walking trials. The use of subject-specific models with minimization of the peak tibiofemoral contact forces improved the accuracy of medial contact forces by 47% and lateral contact forces by 7%, respectively compared with the use of generic musculoskeletal model.     

Patellofemoral joint forces

In this review of patellofemoral joint forces as they might apply to implant design, methodologies for estimating forces on the patella and estimates of the forces, as reported in the literature, are summarized. Two methodologies exist for studying joint loads; one that measures kinematics in-vivo and uses analysis to estimate the joint loads and another that measures ground reaction forces and uses analysis to estimate the joint loads. In both these analyses many assumptions are required with varying degrees of uncertainty; here, those assumptions are examined with data from the published literature. The topics covered include: relationships between quadriceps forces and patellofemoral forces or patella ligament forces, relationships between knee joint moments and quadriceps forces, knee joint moments in various gaits, relationships between patellofemoral forces and lateral subluxation forces, and relationships between patella forces and inferior-superior forces. In many cases, there is little data on patella forces during normal activities, in other cases, there are some discrepancies in reported patella forces, i.e. during squat.     

The coupled motion of the femur and patella during in vivo weightbearing knee flexion

The movement of the knee joint consists of a coupled motion between the tibiofemoral and patellofemoral articulations. This study measured the six degrees-of-freedom kinematics of the tibia, femur, and patella using dual-orthogonal fluoroscopy and magnetic resonance imaging. Ten normal knees from ten living subjects were investigated during weightbearing flexion from full extension to maximum flexion. The femoral and the patellar motions were measured relative to the tibia. The femur externally rotated by 12.9 deg and the patella tilted laterally by 16.3 deg during the full range of knee flexion. Knee flexion was strongly correlated with patellar flexion (R(2)=0.91), posterior femoral translation was strongly correlated to the posterior patellar translation (R(2)=0.87), and internal-external rotation of the femur was correlated to patellar tilt (R(2)=0.73) and medial-lateral patellar translation (R(2)=0.63). These data quantitatively indicate a kinematic coupling between the tibia, femur, and patella, and provide base line information on normal knee joint kinematics throughout the full range of weightbearing flexion. The data also suggest that the kinematic coupling of tibia, femur, and patella should be considered when investigating patellar pathologies and when developing surgical techniques to treat knee joint diseases.     

Regulation of the patellofemoral contact area: An essential mechanism in patellofemoral joint mechanics?

Although the relationship between contact area and pressure under physiological loading has been described in the feline patellofemoral joint, this interaction has only been examined under simplified loading conditions and/or considerably lower forces than those occurring during demanding activities in humans. We hypothesized that patellofemoral contact area increases non-linearly under an increasing joint reaction force to regulate patellofemoral pressure. Eight human cadaveric knees were ramp loaded with muscle forces representative of the stance phase of stair climbing at 30° knee flexion. Continuous pressure data were acquired with a pressure sensitive film that was positioned within the patellofemoral joint. While pressure was linearly dependent upon the resulting joint reaction force, contact area asymptotically approached a maximum value and reached 95% of this maximum at patellofemoral forces of 349–723 N (95% CI). Our findings indicate that the regulatory influence of increasing contact area to protect against high patellofemoral pressure is exhausted at relatively low loads.     

A methodology to accurately quantify patellofemoral cartilage contact kinematics by combining 3D image shape registration and cine-PC MRI velocity data

Patellofemoral osteoarthritis and its potential precursor patellofemoral pain syndrome (PFPS) are common, costly, and debilitating diseases. PFPS has been shown to be associated with altered patellofemoral joint mechanics; however, an actual variation in joint contact stresses has not been established due to challenges in accurately quantifying in vivo contact kinematics (area and location). This study developed and validated a method for tracking dynamic, in vivo cartilage contact kinematics by combining three magnetic resonance imaging (MRI) techniques, cine-phase contrast (CPC), multi-plane cine (MPC), and 3D high-resolution static imaging. CPC and MPC data were acquired from 12 healthy volunteers while they actively extended/flexed their knee within the MRI scanner. Since no gold standard exists for the quantification of in vivo dynamic cartilage contact kinematics, the accuracy of tracking a single point (patellar origin relative to the femur) represented the accuracy of tracking the kinematics of an entire surface. The accuracy was determined by the average absolute error between the PF kinematics derived through registration of MPC images to a static model and those derived through integration of the CPC velocity data. The accuracy ranged from 0.47 mm to 0.77 mm for the patella and femur and from 0.68 mm to 0.86 mm for the patellofemoral joint. For purely quantifying joint kinematics, CPC remains an analytically simpler and more accurate (accuracy <0.33 mm) technique. However, for application requiring the tracking of an entire surface, such as quantifying cartilage contact kinematics, this combined imaging approach produces accurate results with minimal operator intervention.     

A method of calculating physiologically relevant joint reaction forces during forward dynamic simulations of movement from an existing knee model

In the commonly used SIMM software, which includes a complete musculoskeletal model of the lower limbs, the reaction forces at the knee are computed. These reaction forces represent the bone-on-bone contact forces and the soft tissue forces (e.g. ligaments) other than muscles acting at the joint. In the knee model integrated into this software, a patellotibial joint rather than a patellofemoral joint is defined, and a force acting along the direction of the patellar ligament is not included. Although this knee model results in valid kinematics and muscle moment arms, the reaction forces at the knee calculated do not represent physiologic knee joint reaction forces. Hence our objectives were to develop a method of calculating physiologic knee joint reaction forces using the knee model incorporated into the SIMM software and to demonstrate the differences in the forces returned by SIMM and the physiologic forces in an example. Our method converts the anatomically fictional patellotibial joint into a patellofemoral joint and computes the force in an inextensible patellar ligament. In our example, the rectus femoris was fully excited isometrically, with the knee and hip flexed to 90 degrees . The resulting SIMM tibiofemoral joint reaction force was primarily shear, because the quadriceps force was applied to the tibia via the fictional patellotibial joint. In contrast the physiologic tibiofemoral joint reaction force was primarily compression, because the quadriceps force was applied through the patellar ligament. This result illustrates that the physiologic knee joint reaction forces are profoundly different than the forces returned by SIMM. However physiologic knee joint reaction forces can be computed with postprocessing of SIMM results.     

A mathematical model of the patellofemoral joint

A mathematical model of the patellofemoral joint taking into account movements and forces in the sagittal plane is described. The system parameters of the model are the locations of the attachments of the quadriceps muscle and the patellar ligament, the length of the patellar ligament, the dimensions of the patella and the geometry of the articulating surfaces. They were obtained from ten autopsy knees. The model enables calculation of the relative position of the patella, patellar ligament and quadriceps tendon, the location of the patellofemoral contact point and the magnitude of the patellofemoral compression force and the force in the patellar ligament as a function of the location of the tibial tuberosity at different flexion-extension angles of the knee. The model is validated by comparing model data with experimentally determined data.     

Computational biodynamics of human knee joint in gait: from muscle forces to cartilage stresses

Using a validated finite element model of the intact knee joint we aim to compute muscle forces and joint response in the stance phase of gait. The model is driven by reported in vivo kinematics-kinetics data and ground reaction forces in asymptomatic subjects. Cartilage layers and menisci are simulated as depth-dependent tissues with collagen fibril networks. A simplified model with less refined mesh and isotropic depth-independent cartilage is also considered to investigate the effect of model accuracy on results. Muscle forces and joint detailed response are computed following an iterative procedure yielding results that satisfy kinematics/kinetics constraints while accounting at deformed configurations for muscle forces and passive properties. Predictions confirm that muscle forces and joint response alter substantially during the stance phase and that a simplified joint model may accurately be used to estimate muscle forces but not necessarily contact forces/areas, tissue stresses/strains, and ligament forces. Predictions are in general agreement with results of earlier studies. Performing the analyses at 6 periods from beginning to the end (0%, 5%, 25%, 50%, 75% and 100%), hamstrings forces peaked at 5%, quadriceps forces at 25% whereas gastrocnemius forces at 75%. ACL Force reached its maximum of 343 N at 25% and decreased thereafter. Contact forces reached maximum at 5%, 25% and 75% periods with the medial compartment carrying a major portion of load and experiencing larger relative movements and cartilage strains. Much smaller contact stresses were computed at the patellofemoral joint. This novel iterative kinematics-driven model is promising for the joint analysis in altered conditions.     

Biomechanics of changes in ACL and PCL material properties or prestrains in flexion under muscle force-implications in ligament reconstruction

The effects of changes in cruciate ligament material and prestrain on knee joint biomechanics following ligament reconstruction surgery by a tendon are not adequately known. A 3D nonlinear finite element model of the entire knee joint was used to investigate the joint response at different flexion angles under a quadriceps force while varying ACL and PCL initial strains or material properties. The ACL and PCL forces as well as tibiofemoral contact forces/areas substantially increased with greater ACL or PCL initial strains or stiffness. The patellofemoral contact force slightly increased whereas the tibial extensor moment slightly decreased with tenser or stiffer ACL. Reverse trends were predicted with slacker ACL. Results confirm the hypotheses that changes in the prestrain of one cruciate ligament substantially influence the force in the other cruciate ligament and the entire joint and that the use of the patellar tendon (PT) as a replacement for cruciate ligaments markedly alters the joint biomechanics with trends similar to those predicted when increasing prestrains. Forces in both ACL and PCL ligaments increased as one of them became tenser or stiffer and diminished as it became slacker. These results have important consequences in joint biomechanics following ligament injuries or replacement and tend to recommend the use of grafts with smaller prestrains (i.e. slacker than intact) when using the PT as the replacement material with stiffness greater than that of replaced ligament itself.     

A planar model of the knee joint to characterize the knee extensor mechanism

A simple planar static model of the knee joint was developed to calculate effective moment arms for the quadriceps muscle. A pathway for the instantaneous center of rotation was chosen that gives realistic orientations of the femur relative to the tibia. Using the model, nonlinear force and moment equilibrium equations were solved at one degree increments for knee flexion angles from 0 (full extension) to 90 degrees, yielding patellar orientation, patellofemoral contact force and patellar ligament force and direction with respect to both the tibial insertion point and the tibiofemoral contact point. The computer-derived results from this two-dimensional model agree with results from more complex models developed previously from experimentally obtained data. Due to our model's simplicity, however, the operation of the patellar mechanism as a lever as well as a spacer is clearly illustrated. Specifically, the thickness of the patella was found to increase the effective moment arm significantly only at flexions below 35 degrees even though the actual moment arm exhibited an increase throughout the flexion range. Lengthening either the patella or the patellar ligament altered the force transmitted from the quadriceps to the patellar ligament, significantly increasing the effective moment arm at flexions greater than 25 degrees. We conclude that the levering action of the patella is an essential mechanism of knee joint operation at moderate to high flexion angles.     

Biomechanics of changes in ACL and PCL material properties or prestrains in flexion under muscle force-implications in ligament reconstruction

The effects of changes in cruciate ligament material and prestrain on knee joint biomechanics following ligament reconstruction surgery by a tendon are not adequately known. A 3D nonlinear finite element model of the entire knee joint was used to investigate the joint response at different flexion angles under a quadriceps force while varying ACL and PCL initial strains or material properties. The ACL and PCL forces as well as tibiofemoral contact forces/areas substantially increased with greater ACL or PCL initial strains or stiffness. The patellofemoral contact force slightly increased whereas the tibial extensor moment slightly decreased with tenser or stiffer ACL. Reverse trends were predicted with slacker ACL. Results confirm the hypotheses that changes in the prestrain of one cruciate ligament substantially influence the force in the other cruciate ligament and the entire joint and that the use of the patellar tendon (PT) as a replacement for cruciate ligaments markedly alters the joint biomechanics with trends similar to those predicted when increasing prestrains. Forces in both ACL and PCL ligaments increased as one of them became tenser or stiffer and diminished as it became slacker. These results have important consequences in joint biomechanics following ligament injuries or replacement and tend to recommend the use of grafts with smaller prestrains (i.e. slacker than intact) when using the PT as the replacement material with stiffness greater than that of replaced ligament itself.     

Influence of patella alta on knee extensor mechanics

The purpose of this study was to compare the knee extensor mechanics in persons with and without patella alta. Thirteen subjects with patella alta and 14 subjects with normal patellar position participated in the study. Sagittal and axial MR images of the knee were acquired at 0°, 20°, 40°, and 60° of knee flexion. Measurements of actual moment arm, patellar ligament/quadriceps tendon force ratio, quadriceps effective moment arm, and joint reaction force/quadriceps force ratio were obtained. There were no differences between groups in terms of actual moment arm. However, subjects with patella alta had significantly larger patellar ligament/quadriceps tendon force ratios (1.04±0.02 vs. 0.92±0.02) and quadriceps effective moment arms (4.40±0.09 vs. 4.00±0.09 cm) when averaged across the range of knee flexion angles tested. There was no difference in the joint reaction force/quadriceps force ratio between groups. The observed differences in knee extensor mechanics suggest that individuals with patella alta have a more efficient knee extensor mechanism and would be expected to generate similar joint reaction forces per unit quadriceps force compared to subjects with normal patellar position. Therefore, persons with patella alta may experience less patellofemoral joint reaction force to overcome the same knee flexion moment in the range of 0°–60° of knee flexion.     

Normalized patellofemoral joint reaction force is greater in individuals with patellofemoral pain

Patellofemoral pain is a disabling, highly prevalent pathology. Altered patellofemoral contact forces are theorized to contribute to this pain. Musculoskeletal modeling has been employed to better understand the etiology of patellofemoral pain. Currently, there are no data on the effective quadriceps moment arm for individuals with patellofemoral pain, forcing researchers to apply normative values when modeling such individuals. In addition, the ratio of patellofemoral reaction force to quadriceps force is often used as a surrogate for patellofemoral joint contact force, ignoring the fact that the quadriceps efficiency can vary with pathology and intervention. Thus, the purposes of this study were to: (1) quantify the effective quadriceps moment arm in individuals with patellofemoral pain and compare this value to a control cohort and (2) develop a novel methodology for quantifying the normalized patellofemoral joint reaction force in vivo during dynamic activities. Dynamic MR data were captured as subjects with patellofemoral pain (30F/3M) cyclically flexed their knee from 10° to 40°. Data for control subjects (29F/9M) were taken from a previous study. The moment arm data acquired across a large cohort of individuals with patellofemoral pain should help advance musculoskeletal modeling. The primary finding of this study was an increased mean normalized patellofemoral reaction force of 14.9% (maximum values at a knee angle of 10°) in individuals with patellofemoral pain. Understanding changes in the normalized patellofemoral reaction force with pathology may lead to improvements in clinical decision making, and consequently treatments, by providing a more direct measure of altered patellofemoral joint forces.     

A mathematical model for the evaluation of the behaviour during flexion of condylar-type knee prostheses

A 3D knee model was developed in order to evaluate the mechanical behaviour during flexion of condylar-type knee prosthesis. Based on the total energy minimization principle, it takes into account the articular surfaces (the tibial surface being deformable), the body weight, and the patello femoral joint. It generates the kinematics of the joint, the motion of the centre of contact, the quadriceps forces, the pressure distribution on the tibial plateau, and ligament lengths and forces between 0 and 120 degrees of flexion. The results for ten digitized knees and the commercially available prostheses are presented. They are in general agreement with experimental results published in the literature. It is concluded that this computer program may be, within its limitations, a useful tool in the preliminary evaluation of new condylar-type knee prosthesis designs.