Evaluation of a computational model used to predict the patellofemoral contact pressure distribution

One possible cause of patellofemoral pain syndrome is excessive lateral force acting on the patella. Although several treatment methods focus on decreasing the lateral force acting on the patella, the relationship between the lateral force and the patellofemoral contact pressure distribution is unclear. A computational model has been developed to determine how loading variations alter the patellofemoral force and pressure distributions for individual knees. The model allows variation in the quadriceps and patella tendon forces, and calculates the predicted contact pressure distribution using the discrete element analysis technique. To characterize the accuracy of the model, four cadaver knees were flexed on a knee simulator with three initial Q-angles, while recording the force and pressure distributions with a pressure sensor. A model of each knee was created from CT data. Using the external force applied to the knee, the geometry of the knee, and the quadriceps origin as input, the pressure distribution was calculated during flexion. Similar trends were noted for the computational and experimental results. The percentage of the total force applied to the lateral cartilage increased with the Q-angle. The maximum contact pressure increased during flexion. The maximum lateral contact pressure increased with the Q-angle for three knees. For the other knee, increasing the Q-angle decreased the maximum lateral pressure. The maximum medial contact pressure decreased as the Q-angle increased. By characterizing the influence of patellofemoral loading on the force and pressure distributions, the computational model could be used to evaluate treatment methods prescribed for patellofemoral pain.     

Finite element analysis to characterize how varying patellar loading influences pressure applied to cartilage: model evaluation

A finite element analysis (FEA) modeling technique has been developed to characterize how varying the orientation of the patellar tendon influences the patellofemoral pressure distribution. To evaluate the accuracy of the technique, models were created from MRI images to represent five knees that were previously tested in vitro to determine the influence of hamstrings loading on patellofemoral contact pressures. Hamstrings loading increased the lateral and posterior orientation of the patellar tendon. Each model was loaded at 40°, 60°, and 80° of flexion with quadriceps force vectors representing the experimental loading conditions. The orientation of the patellar tendon was represented for the loaded and unloaded hamstrings conditions based on experimental measures of tibiofemoral alignment. Similar to the experimental data, simulated loading of the hamstrings within the FEA models shifted the center of pressure laterally and increased the maximum lateral pressure. Significant (p < 0.05) differences were identified for the center of pressure and maximum lateral pressure from paired t-tests carried out at the individual flexion angles. The ability to replicate experimental trends indicates that the FEA models can be used for future studies focused on determining how variations in the orientation of the patellar tendon related to anatomical or loading variations or surgical procedures influence the patellofemoral pressure distribution.     

In vivo patellar tracking induced by individual quadriceps components in individuals with patellofemoral pain

Patellofemoral pain is a common knee disorder with a multi-factorial etiology related to abnormal patellar tracking. Our hypothesis was that the pattern of three-dimensional rotation and translation of the patella induced by selective activation of individual quadriceps components would differ between subjects with patellofemoral pain and healthy subjects. Nine female subjects with patellofemoral pain and seven healthy female subjects underwent electrical stimulation to selectively activate individual quadriceps components (vastus medialis obliquus, VMO; vastus medialis lateralis, VML; vastus lateralis, VL) with the knee at 0° and 20° flexion, while three-dimensional patellar tracking was recorded. Normalized direction of rotation and direction of translation characterized the relative amplitudes of each component of patellar movement. VMO activation in patellofemoral pain caused greater medial patellar rotation (distal patellar pole rotates medially in frontal plane) at both knee positions (p<0.01), and both VMO and VML activation caused increased anterior patellar translation (p<0.001) in patellofemoral pain compared to healthy subjects at 20° knee flexion. VL activation caused more lateral patellar translation (p<0.001) in patellofemoral pain compared to healthy subjects. In healthy subjects the 3-D mechanical action of the VMO is actively modulated with knee flexion angle while such modulation was not observed in PFP subjects. This could be due to anatomical differences in the VMO insertion on the patella and medial quadriceps weakness. Quantitative evaluation of the influence of individual quadriceps components on patellar tracking will aid understanding of the knee extensor mechanism and provide insight into the etiology of patellofemoral pain.     

Comparing two estimations of the quadriceps force distribution for use during patellofemoral simulation

EMG analysis has indicated that the vastus lateralis and vastus medialis contribute less to the quadriceps moment during knee extension than the physiological cross-sectional areas (PCSA's) of the muscles indicate. Both PCSA- and EMG-based quadriceps force distributions were utilized while computationally simulating knee extension. For both distributions, a 10 degrees increase in the Q-angle and a 50% decrease in the force applied by the vastus medialis were simulated, and the influence of these changes on the resultant force and moment applied by the quadriceps muscles and the patella tendon was quantified. For both quadriceps force distributions, increasing the Q-angle increased the lateral force and the moment acting to rotate the distal patella laterally. Due to the relatively large forces initially attributed to the vastus medialis and vastus lateralis for the PCSA-based quadriceps force distribution, decreasing the vastus medialis force created a large force imbalance between these two muscles. For the PCSA-based quadriceps force distribution, decreasing the vastus medialis force increased the lateral rotation moment and the moment acting to tilt the patella laterally. For the EMG-based quadriceps force distribution, decreasing the vastus medialis force produced relatively little change in the tilt and rotation moments. For both quadriceps force distributions, increasing the Q-angle increased the maximum and mean cartilage pressure during flexion, but decreasing the vastus medialis force only increased the cartilage pressures for the PCSA-based quadriceps distribution. The choice of the initial quadriceps distribution can influence the outcome of patellofemoral simulation when manipulating quadriceps muscle forces.     

Dynamic simulation of tibial tuberosity realignment: model evaluation

This study was performed to evaluate a dynamic multibody model developed to characterize the influence of tibial tuberosity realignment procedures on patellofemoral motion and loading. Computational models were created to represent four knees previously tested at 40°, 60°, and 80° of flexion with the tibial tuberosity in a lateral, medial and anteromedial positions. The experimentally loaded muscles, major ligaments of the knee, and patellar tendon were represented. A repeated measures ANOVA with post-hoc testing was performed at each flexion angle to compare data between the three positions of the tibial tuberosity. Significant experimental trends for decreased patella flexion due to tuberosity anteriorization and a decrease in the lateral contact force due to tuberosity medialization were reproduced computationally. The dynamic multibody modeling technique will allow simulation of function for symptomatic knees to identify optimal surgical treatment methods based on parameters related to knee pathology and pre-operative kinematics.     

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.     

Stair climbing results in more challenging patellofemoral contact mechanics and kinematics than walking at early knee flexion under physiological-like quadriceps loading

The mechanical environment during stair climbing has been associated with patellofemoral pain, but the contribution of loading to this condition is not clearly understood. It was hypothesized that the loading conditions during stair climbing induce higher patellofemoral pressures, a more lateral force distribution on the trochlea and a more lateral shift and tilt of the patella compared to walking at early knee flexion. Optical markers for kinematic measurements were attached to eight cadaveric knees, which were loaded with muscle forces at instances of walking and stair climbing cycles at 12° and 30° knee flexion. Contact mechanics were determined using a pressure sensitive film. At 12° knee flexion, stair climbing loads resulted in higher peak pressure (p=0.012) than walking, more lateral force distribution (p=0.012) and more lateral tilt (p=0.012), whilst mean pressure (p=0.069) and contact area (p=0.123) were not significantly different. At 30° knee flexion, although stair climbing compared to walking loads resulted in significantly higher patellofemoral mean (p=0.012) and peak pressures (p=0.012), contact area (p=0.025), as well as tilt (p=0.017), the medial–lateral force distribution (p=0.674) was not significantly different. No significant differences were observed in patellar shift between walking and stair climbing at either 12° (p=0.093) or 30° (p=0.575) knee flexion. Stair climbing thus leads to more challenging patellofemoral contact mechanics and kinematics than level walking at early knee flexion. The increase in patellofemoral pressure, lateral force distribution and lateral tilt during stair climbing provides a possible biomechanical explanation for the patellofemoral pain frequently experienced during this activity.     

The effect of vastus medialis forces on patello-femoral contact: a model-based study

A mathematical model of the patello-femoral joint was introduced to investigate the impact of the vastus medialis (longus, obliquus) forces on the lateral contact force levels. In the model, the quadriceps were represented as five separate forces: vastus lateralis, vastus intermedius, rectus femoris, vastus medialis longus (VML), and obliquus (VMO). By varying the relative force generation ratios of the quadriceps heads, the patello-femoral contact forces were estimated. We sought to analytically determine the range of forces in the VMO and VML that cause a reduction or an increase of lateral contact forces, often the cause of patello-femoral pain. Our results indicated that increased contact forces are more dependent on combinations of muscle forces than solely VMO weakness. Moreover, our simulation data showed that the contact force levels are also highly dependent on the knee flexion angle. These findings suggest that training targeted to reduce contact forces through certain joint angles could actually result in a significant increase of the contact forces through other joint angles.     

Patellofemoral cartilage stresses are most sensitive to variations in vastus medialis muscle forces

The purpose of this study was to evaluate the effects of variations in quadriceps muscle forces on patellofemoral stress. We created subject-specific finite element models for 21 individuals with chronic patellofemoral pain and 16 pain-free control subjects. We extracted three-dimensional geometries from high resolution magnetic resonance images and registered the geometries to magnetic resonance images from an upright weight bearing squat with the knees flexed at 60°. We estimated quadriceps muscle forces corresponding to 60° knee flexion during a stair climb task from motion analysis and electromyography-driven musculoskeletal modelling. We applied the quadriceps muscle forces to our finite element models and evaluated patellofemoral cartilage stress. We quantified cartilage stress using an energy-based effective stress, a scalar quantity representing the local stress intensity in the tissue. We used probabilistic methods to evaluate the effects of variations in quadriceps muscle forces from five trials of the stair climb task for each subject. Patellofemoral effective stress was most sensitive to variations in forces in the two branches of the vastus medialis muscle. Femur cartilage effective stress was most sensitive to variations in vastus medialis forces in 29/37 (78%) subjects, and patella cartilage effective stress was most sensitive to variations in vastus medialis forces in 21/37 (57%) subjects. Femur cartilage effective stress was more sensitive to variations in vastus medialis longus forces in subjects classified as maltrackers compared to normal tracking subjects (p?=?0.006). This study provides new evidence of the importance of the vastus medialis muscle in the treatment of patellofemoral pain.     

The effect on patellofemoral joint stability of selective cutting of lateral retinacular and capsular structures

Patient selection for lateral retinacular release (LRR) and its efficacy are controversial. Iatrogenic medial subluxation can occur with inappropriate LRR. The aim of this study was to determine the reduction in patellofemoral stability with progressively more extensive LRR. The force required to displace the patella 10 mm medially and laterally in nine cadaveric knees was measured with and without loading of the quadriceps and iliotibial band. The knee was tested intact, then after progressive release beginning proximal to the patella (PR), the mid-level between the proximal and distal limit of the patella (MR) where the fibres are more transverse, then distally till Gerdy's tubercle (DR) and finally the joint capsule (CR). Both medial and lateral stability decreased with progressive releases, larger for the medial. The MR caused a significant reduction of lateral stability between 30° and 90° of knee flexion. There was an 8% reduction in medial stability at 0° flexion with a complete LRR (DR). A comparable reduction in medial stability in the loaded knee at 20° and 30° flexion was obtained with MR alone, with no further reduction after DR. A capsular release caused a further reduction in medial stability at 0° and 20° and this was marked in the unloaded knee. In extension, the main lateral restraint was the joint capsule. At 30° flexion, the transverse fibres were the main contributor to the lateral restraint.     

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 geometric approach to study the contact mechanisms in the patellofemoral joint of normal versus patellofemoral pain syndrome subjects

The biomechanics of the patellofemoral (PF) joint is complex in nature, and the aetiology of such manifestations of PF instability as patellofemoral pain syndrome (PFPS) is still unclear. At this point, the particular factors affecting PFPS have not yet been determined. This study has two objectives: (1) The first is to develop an alternative geometric method using a three-dimensional (3D) registration technique and linear mapping to investigate the PF joint contact stress using an indirect measure: the depth of virtual penetration (PD) of the patellar cartilage surface into the femoral cartilage surface. (2) The second is to develop 3D PF joint models using the finite element analysis (FEA) to quantify in vivo cartilage contact stress and to compare the peak contact stress location obtained from the FE models with the location of the maximum PD. Magnetic resonance images of healthy and PFPS subjects at knee flexion angles of 15°, 30° and 45° during isometric loading have been used to develop the geometric models. The results obtained from both approaches demonstrated that the subjects with PFPS show higher PD and contact stresses than the normal subjects. Maximum stress and PD increase with flexion angle, and occur on the lateral side in healthy and on the medial side in PFPS subjects. It has been concluded that the alternative geometric method is reliable in addition to being computationally efficient compared with FEA, and has the potential to assess the mechanics of PFPS with an accuracy similar to the FEA.     

Activation of the VMO and VL during dynamic mini-squat exercises with and without isometric hip adduction

Objective: the purpose of this study was to compare vastus medialis obliquus (VMO) and vastus lateralis (VL) activity while performing a mini-squat with and without isometric hip adduction.

Design and setting: a repeated measures within subjects design was used. Subjects performed two sets of three repetitions of a traditional mini-squat and a mini-squat with concurrent hip adduction (squeeze).

Subjects: 20 recreationally active subjects (10 men, 10 women AGE=28.10±5.91 years, HEIGHT=170.94±11.03 cm, MASS=72.32±16.66 kg) with no history of patellofemoral pain (PFP), quadriceps injury, or other knee injury participated in the study.

Measurements: the EMG signal of the VMO and VL was recorded bilaterally during both exercises. EMG data were normalized to the maximal voluntary isometric contraction (MVIC) of the quadriceps produced during seated, isometric knee extension.

Results: results of repeated measures ANOVA's revealed that the squeeze squat produced significantly greater VMO and VL activity than the traditional squat (p=0.02). For both the traditional and squeeze squats, intrasession reliability from the first to the second set was calculated using intraclass correlation coefficient (ICC) formula (3:1) bilaterally for both the VMO and the VL. All ICC values were greater than 0.9.

Conclusion: combining isometric hip adduction with a mini-squat exercise significantly increases the activity of the quadriceps. Performing mini-squats with isometric hip adduction will be beneficial to patellofemoral patients as they increase quadriceps activity, however, based on our data we cannot conclude that this exercise preferentially recruits the VMO. Further research is needed to determine the exact mechanism by which quadriceps function is altered.     

Combined measurement and modeling of specimen-specific knee mechanics for healthy and ACL-deficient conditions

Quantifying the mechanical environment at the knee is crucial for developing successful rehabilitation and surgical protocols. Computational models have been developed to complement in vitro studies, but are typically created to represent healthy conditions, and may not be useful in modeling pathology and repair. Thus, the objective of this study was to create finite element (FE) models of the natural knee, including specimen-specific tibiofemoral (TF) and patellofemoral (PF) soft tissue structures, and to evaluate joint mechanics in intact and ACL-deficient conditions. Simulated gait in a whole joint knee simulator was performed on two cadaveric specimens in an intact state and subsequently repeated following ACL resection. Simulated gait was performed using motor-actuated quadriceps, and loads at the hip and ankle. Specimen-specific FE models of these experiments were developed in both intact and ACL-deficient states. Model simulations compared kinematics and loading of the experimental TF and PF joints, with average RMS differences [max] of 3.0° [8.2°] and 2.1° [8.4°] in rotations, and 1.7 [3.0] and 2.5 [5.1] mm in translations, for intact and ACL-deficient states, respectively. The timing of peak quadriceps force during stance and swing phase of gait was accurately replicated within 2° of knee flexion and with an average error of 16.7% across specimens and pathology. Ligament recruitment patterns were unique in each specimen; recruitment variability was likely influenced by variations in ligament attachment locations. ACL resections demonstrated contrasting joint mechanics in the two specimens with altered knee motion shown in one specimen (up to 5 mm anterior tibial translation) while increased TF joint loading was shown in the other (up to 400 N).     

Development and validation of a weight-bearing finite element model for total knee replacement

Total knee arthroplasty (TKA) is a successful procedure for osteoarthritis. However, some patients (19%) do have pain after surgery. A finite element model was developed based on boundary conditions of a knee rig. A 3D-model of an anatomical full leg was generated from magnetic resonance image data and a total knee prosthesis was implanted without patella resurfacing. In the finite element model, a restarting procedure was programmed in order to hold the ground reaction force constant with an adapted quadriceps muscle force during a squat from 20° to 105° of flexion. Knee rig experimental data were used to validate the numerical model in the patellofemoral and femorotibial joint. Furthermore, sensitivity analyses of Young’s modulus of the patella cartilage, posterior cruciate ligament (PCL) stiffness, and patella tendon origin were performed. Pearson’s correlations for retropatellar contact area, pressure, patella flexion, and femorotibial ap-movement were near to 1. Lowest root mean square error for retropatellar pressure, patella flexion, and femorotibial ap-movement were found for the baseline model setup with Young’s modulus of 5 MPa for patella cartilage, a downscaled PCL stiffness of 25% compared to the literature given value and an anatomical origin of the patella tendon. The results of the conducted finite element model are comparable with the experimental results. Therefore, the finite element model developed in this study can be used for further clinical investigations and will help to better understand the clinical aspects after TKA with an unresurfaced patella.     

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.     

Iliotibial band tension affects patellofemoral and tibiofemoral kinematics

The iliotibial band (ITB) has an important role in knee mechanics and tightness can cause patellofemoral maltracking. This study investigated the effects of increasing ITB tension on knee kinematics. Nine fresh-frozen cadaveric knees had the components of the quadriceps loaded with 175 N. A Polaris optical tracking system was used to acquire joint kinematics during extension from 100° to 0° flexion. This was repeated after the following ITB loads: 30, 60 and 90 N. There was no change with 30 N load for patellar translation. On average, at 60 and 90 N, the patella translated laterally by 0.8 and 1.4 mm in the mid flexion range compared to the ITB unloaded condition. The patella became more laterally tilted with increasing ITB loads by 0.7°, 1.2° and 1.5° for 30, 60 and 90 N, respectively. There were comparable increases in patellar lateral rotation (distal patella moves laterally) towards the end of the flexion cycle. Increased external rotation of the tibia occurred from early flexion onwards and was maximal between 60° and 75° flexion. The increase was 5.2°, 9.5° and 13° in this range for 30, 60 and 90 N, respectively. Increased tibial abduction with ITB loads was not observed. The combination of increased patellar lateral translation and tilt suggests increased lateral cartilage pressure. Additionally, the increased tibial external rotation would increase the Q angle. The clinical consequences and their relationship to lateral retinacular releases may be examined, now that the effects of a tight ITB are known.     

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.